Gianluca De Santis

Patient-specific computational fluid dynamics: structured mesh generation from coronary angiography

Patient-specific simulations are widely used to investigate the local hemodynamics within realistic morphologies. However, pre-processing and mesh generation are time consuming, operator dependent, and the quality of the resulting mesh is often suboptimal. Therefore, a semi-automatic methodology for patient-specific reconstruction and structured meshing of a left coronary tree from biplane angiography is presented. Seven hexahedral grids have been generated with the new method (50,000–3,200,000 cells) and compared to nine unstructured tetrahedral grids with prismatic boundary layer (150,000–3,100,000 cells). Steady-state blood flow simulation using Computational Fluid Dynamics (CFD) has been used to calculate the Wall Shear Stress (WSS). Our results (99 percentile, area-weighted and local WSS values along a line) demonstrate that hexahedral meshes with respect to tetrahedral/prismatic meshes converge better, and for the same accuracy of the result, six times less cells and 14 times less computational time are required. Hexahedral meshes are superior to tetrahedral/prismatic meshes and should be preferred for the calculation of the WSS.

An Integrated Framework to Quantitatively Link Mouse-Specific Hemodynamics to Aneurysm Formation in Angiotensin II-infused ApoE −/− mice

Locally disturbed flow has been suggested to play a (modulating) role in abdominal aortic aneurysm (AAA) formation, but no longitudinal studies have been performed yet due to (a.o.) a lack of human data prior to AAA formation. In this study we made use of recent advances in small animal imaging technology in order to set up entirely mouse-specific computational fluid dynamics (CFD) simulations of the abdominal aorta in an established ApoE −/− mouse model of AAA formation, combining (i) in vivo contrast-enhanced micro-CT scans (geometrical model) and (ii) in vivo high-frequency ultrasound scans (boundary conditions). Resulting areas of disturbed flow at baseline were compared to areas of AAA at end-stage. Qualitative results showed that AAA dimension is maximal in areas that are situated proximal to those areas that experience most disturbed flow in three out of four S developing an AAA. Although further quantitative analysis did not reveal any obvious relationship between areas that experience most disturbed flow and the end-stage AAA dimensions, we cannot exclude that hemodynamics play a role in the initial phases of AAA formation. Due to its mouse-specific and in vivo nature, the presented methodology can be used in future research to link detailed and animal-specific (baseline) hemodynamics to (end-stage) arterial disease in longitudinal studies in mice.

Impact of carotid stent cell design on vessel scaffolding: A case study comparing experimental investigation and numerical simulations

Purpose To quantitatively evaluate the impact of carotid stent cell design on vessel scaffolding by using patient-specific finite element analysis of carotid artery stenting (CAS). Methods The study was organized in 2 parts: (1) validation of a patient-specific finite element analysis of CAS and (2) evaluation of vessel scaffolding. Micro-computed tomography (CT) images of an open-cell stent deployed in a patient-specific silicone mock artery were compared with the corresponding finite element analysis results. This simulation was repeated for the closed-cell counterpart. In the second part, the stent strut distribution, as reflected by the inter-strut angles, was evaluated for both cell types in different vessel cross sections as a measure of scaffolding. Results The results of the patient-specific finite element analysis of CAS matched well with experimental stent deployment both qualitatively and quantitatively, demonstrating the reliability of the numerical approach. The measured inter-strut angles suggested that the closed-cell design provided superior vessel scaffolding compared to the open-cell counterpart. However, the full strut interconnection of the closed-cell design reduced the stents ability to accommodate to the irregular eccentric profile of the vessel cross section, leading to a gap between the stent surface and the vessel wall. Conclusion Even though this study was limited to a single stent design and one vascular anatomy, the study confirmed the capability of dedicated computer simulations to predict differences in scaffolding by open- and closed-cell carotid artery stents. These simulations have the potential to be used in the design of novel carotid stents or for procedure planning.

The Impact of Simplified Boundary Conditions and Aortic Arch Inclusion on CFD Simulations in the Mouse Aorta: A Comparison With Mouse-specific Reference Data

Computational fluid dynamics (CFD) simulations allow for calculation of a detailed flow field in the mouse aorta and can thus be used to investigate a potential link between local hemodynamics and disease development. To perform these simulations in a murine setting, one often needs to make assumptions (e.g. when mouse-specific boundary conditions are not available), but many of these assumptions have not been validated due to a lack of reference data. In this study, we present such a reference data set by combining high-frequency ultrasound and contrast-enhanced micro-CT to measure (in vivo) the time-dependent volumetric flow waveforms in the complete aorta (including seven major side branches) of 10 male ApoE -/- deficient mice on a C57Bl/6 background. In order to assess the influence of some assumptions that are commonly applied in literature, four different CFD simulations were set up for each animal: (i) imposing the measured volumetric flow waveforms, (ii) imposing the average flow fractions over all 10 animals, presented as a reference data set, (iii) imposing flow fractions calculated by Murrays law, and (iv) restricting the geometrical model to the abdominal aorta (imposing measured flows). We found that - even if there is sometimes significant variation in the flow fractions going to a particular branch - the influence of using average flow fractions on the CFD simulations is limited and often restricted to the side branches. On the other hand, Murrays law underestimates the fraction going to the brachiocephalic trunk and strongly overestimates the fraction going to the distal aorta, influencing the outcome of the CFD results significantly. Changing the exponential factor in Murrays law equation from 3 to 2 (as suggested by several authors in literature) yields results that correspond much better to those obtained imposing the average flow fractions. Restricting the geometrical model to the abdominal aorta did not influence the outcome of the CFD simulations. In conclusion, the presented reference dataset can be used to impose boundary conditions in the mouse aorta in future studies, keeping in mind that they represent a subsample of the total population, i.e., relatively old, non-diseased, male C57Bl/6 ApoE -/- mice.

Patient-Specific Computer Modeling to Predict Aortic Regurgitation after Transcatheter Aortic Valve Replacement

Outcome of transcatheter aortic valve replacement (TAVR) depends on a combination of patient-, procedure-, and operator-related variables. Specific device–host-related interactions may also be involved and may result in, for instance, incomplete and/or nonuniform frame expansion that in turn may

Patient-specific computer modelling of coronary bifurcation stenting: the John Doe programme.

John Doe, an 81-year-old patient with a significant distal left main (LM) stenosis, was treated using a provisional stenting approach. As part of an European Bifurcation Club (EBC) project, the complete stenting procedure was repeated using computational modelling. First, a tailored three-dimensional (3D) reconstruction of the bifurcation anatomy was created by fusion of multislice computed tomography (CT) imaging and intravascular ultrasound. Second, finite element analysis was employed to deploy and post-dilate the stent virtually within the generated patient-specific anatomical bifurcation model. Finally, blood flow was modelled using computational fluid dynamics. This proof-of-concept study demonstrated the feasibility of such patient-specific simulations for bifurcation stenting and has provided unique insights into the bifurcation anatomy, the technical aspects of LM bifurcation stenting, and the positive impact of adequate post-dilatation on blood flow patterns. Potential clinical applications such as virtual trials and preoperative planning seem feasible but require a thorough clinical validation of the predictive power of these computer simulations.

A Computational Study of the Hemodynamic Impact of Open- Versus Closed-Cell Stent Design in Carotid Artery Stenting

The aim of this study is to analyze the shape and flow changes of a patient-specific carotid artery after carotid artery stenting (CAS) performed using an open-cell (stent-O) or a closed-cell (stent-C) stent design. First, a stent reconstructed from micro-computed tomography (microCT) is virtually implanted in a left carotid artery reconstructed from CT angiography. Second, an objective analysis of the stent-to-vessel apposition is used to quantify the lumen cross-sectional area and the incomplete stent apposition (ISA). Third, the carotid artery lumen is virtually perfused in order to quantify its resistance to flow and its exposure to atherogenic or thrombogenic hemodynamic conditions. After CAS, the minimum cross-sectional area of the internal carotid artery (ICA) (external carotid artery [ECA]) changes by +54% (-12%) with stent-O and +78% (-17%) with stent-C; the resistance to flow of the ICA (ECA) changes by -21% (+13%) with stent-O and -26% (+18%) with stent-C. Both stent designs suffer from ISA but the malapposed stent area is larger with stent-O than stent-C (29.5 vs. 14.8 mm(2) ). The untreated vessel is not exposed to atherogenic flow conditions whereas an area of 67.6 mm(2) (104.9) occurs with stent-O (stent-C). The area of the stent surface exposed to thrombogenic risk is 5.42 mm(2) (7.7) with stent-O (stent-C). The computer simulations of stenting in a patients carotid artery reveal a trade-off between cross-sectional size and flow resistance of the ICA (enlarged and circularized) and the ECA (narrowed and ovalized). Such a trade-off, together with malapposition, atherogenic risk, and thrombogenic risk is stent-design dependent.

How an Eukaryotic Cell Senses the Substrate Stiffness? An Exploration Using a Finite Element Model With Cytoskeleton Modelled as Tensegrity Structure

Mammalian cells in vivo are connected to the ECM (or other substrate, SS) or other cells having elastic moduli ranging from 10 to 10000 Pa. Several experimental evidences relate cell processes (e.g. changes of spreading area, cytoskeletal filament assembling, focal adhesion complex (FAC)) to SS stiffness, but how a passive substrate affects these cell processes is still unknown [1,2].Copyright

A Novel Tram Stent Method in the Treatment of Coronary Bifurcation Lesions – Finite Element Study

A novel stent was designed for the treatment of coronary bifurcation lesion, and it was investigated for its performance by finite element analysis. This study was performed in search of a novel method of treatment of bifurcation lesion with provisional stenting. A bifurcation model was created with the proximal vessel of 3.2 mm diameter, and the distal vessel after the side branch (2.3 mm) was 2.7 mm. A novel stent was designed with connection links that had a profile of a tram. Laser cutting and shape setting of the stent was performed, and thereafter it was crimped and deployed over a balloon. The contact pressure, stresses on the arterial wall, stresses on the stent, the maximal principal log strain of the main artery and the side-branch were studied. The study was performed in Abaqus, Simulia. The stresses on the main branch and the distal branch were minimally increased after deployment of this novel stent. The side branch was preserved, and the stresses on the side branch were lesser; and at the confluence of bifurcation on either side of the side branch origin the von-Mises stress was marginally increased. The stresses and strain at the bifurcation were significantly lesser than the stresses and strain of the currently existing techniques used in the treatment of bifurcation lesions though the study was primarily focused only on the utility of the new technology. There is a potential for a novel Tram-stent method in the treatment of coronary bifurcation lesions.

Patient-specific computer simulation to elucidate the role of contact pressure in the development of new conduction abnormalities after catheter-based implantation of a self-expanding aortic valve

Background— The extent to which pressure generated by the valve on the aortic root plays a role in the genesis of conduction abnormalities after transcatheter aortic valve replacement (TAVR) is unknown. This study elucidates the role of contact pressure and contact pressure area in the development of conduction abnormalities after TAVR using patient-specific computer simulations. Methods and Results— Finite-element computer simulations were performed to simulate TAVR of 112 patients who had undergone TAVR with the self-expanding CoreValve/Evolut R valve. On the basis of preoperative multi-slice computed tomography, a patient-specific region of the aortic root containing the atrioventricular conduction system was determined by identifying the membranous septum. Contact pressure and contact pressure index (percentage of area subjected to pressure) were quantified and compared in patients with and without new conduction abnormalities. Sixty-two patients (55%) developed a new left bundle branch block or a high-degree atrioventricular block after TAVR. Maximum contact pressure and contact pressure index (median [interquartile range]) were significantly higher in patients with compared with those without new conduction abnormalities (0.51 MPa [0.43–0.70 MPa] and 33% [22%–44%], respectively, versus 0.29 MPa [0.06–0.50 MPa] and 12% [1%–28%]). By multivariable regression analysis, only maximum contact pressure (odds ratio, 1.35; confidence interval, 1.1–1.7; P=0.01) and contact pressure index (odds ratio, 1.52; confidence interval, 1.1–2.1; P=0.01) were identified as independent predictors for conduction abnormalities, but not implantation depth. Conclusions— Patient-specific computer simulations revealed that maximum contact pressure and contact pressure index are both associated with new conduction abnormalities after CoreValve/Evolut R implantation and can predict which patient will have conduction abnormalities.