Francesco Iannaccone

Our capricious vessels: The influence of stent design and vessel geometry on the mechanics of intracranial aneurysm stent deployment.

There is a growing interest in virtual tools to assist clinicians in evaluating different procedures and devices for endovascular treatment. In the present study we use finite element analysis to investigate the influence of stent design and vessel geometry for stent assisted coiling of intracranial aneurysms. Nine virtual stenting procedures were performed: three nitinol stent designs ((i) an open cell stent resembling the Neuroform, (ii) a generic stiff and (iii) a more flexible closed cell design), were deployed in three patient-specific cerebral aneurysmatic vessels. We investigated the percentage of strut area covering the aneurysm neck, the straightening induced on the cerebrovasculature by the stent placement (quantified by the reduction in tortuosity), and stent apposition to the wall (quantified as the percentage of struts within 0.2mm of the vessel). The results suggest that the open cell design better covers the aneurysm neck (11.0±1.1%) compared to both the stiff (7.8±1.6%) and flexible (8.7±1.6%) closed cell stents, and induces less straightening of the vessel (-5.1±1.6% vs. -42.9±9.8% and -26.9±11.9% ). The open cell design has, however, less struts apposing well to the vessel wall (56.0±6.4%) compared to the flexible (73.4±4.6%) and stiff (70.4±5.1%) closed cell design. With the presented study, we hope to contribute to and improve aneurysm treatment, using a novel patient specific environment as a possible pre-operative tool to evaluate mechanical stent behavior in different vascular geometries.

Virtual evaluation of stent graft deployment: a validated modeling and simulation study.

The presented study details the virtual deployment of a bifurcated stent graft (Medtronic Talent) in an Abdominal Aortic Aneurysm model, using the finite element method. The entire deployment procedure is modeled, with the stent graft being crimped and bent according to the vessel geometry, and subsequently released. The finite element results are validated in vitro with placement of the device in a silicone mock aneurysm, using high resolution CT scans to evaluate the result. The presented work confirms the capability of finite element computer simulations to predict the deformed configuration after endovascular aneurysm repair (EVAR). These simulations can be used to quantify mechanical parameters, such as neck dilations, radial forces and stresses in the device, that are difficult or impossible to obtain from medical imaging.

Fluid-Structure Interaction Simulation of Prosthetic Aortic Valves: Comparison between Immersed Boundary and Arbitrary Lagrangian-Eulerian Techniques for the Mesh Representation.

In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations’ outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results.

Synthetic dataset generation for the analysis and the evaluation of image-based hemodynamics of the human aorta

Here, we consider the issue of generating a suitable controlled environment for the evaluation of phase contrast (PC) MRI measurements. The computational framework, tailored to build synthetic datasets, is based on a two-step approach, i.e., define and implement (1) an accurate CFD model and (2) an image generator able to mime the overall outcomes of a PC MRI acquisition starting from datasets retrieved by the computational model. About 20 different datasets were built by changing relevant image parameters (pixel size, slice thickness, time frames per cardiac cycle). Focusing our attention on the thoracic aorta, synthetic images were processed in order to: (1) verify to which extent the fluid dynamics into the aortic arch is influenced by the image parameters; (2) establish the effect of spatial and temporal interpolation. Our study demonstrates that the integral scale of the aortic bulk flow could be described satisfactorily even when using images which are nowadays acquirable with MRI scanners. However, attention must be paid to near-wall velocities that can be affected by large inaccuracy. In detail, in bulk flow regions error values are well bounded (below 5% for most of the analyzed resolutions), while errors greater than 100% are systematically present at the vessel’s wall. Moreover, also the data interpolation process can be responsible for large inaccuracies in new data generation, due to the inherent complexity of the flow field in some connected regions.

The influence of vascular anatomy on carotid artery stenting: A parametric study for damage assessment

Carotid artery stenting is emerging as an alternative technique to surgery for the treatment of symptomatic severe carotid stenosis. Clinical and experimental evidence demonstrates that both plaque morphology and biomechanical changes due to the device implantation can be possible causes of an unsuccessful treatment. In order to gain further insights of the endovascular intervention, a virtual environment based on structural finite element simulations was built to emulate the stenting procedure on generalized atherosclerotic carotid geometries which included a damage model to quantify the injury of the vessel. Five possible lesion scenarios were simulated by changing both material properties and vascular geometrical features to cover both presumed vulnerable and stable plaques. The results were analyzed with respect to lumen gain and wall stresses which are potentially related to the failure of the procedure according to previous studies. Our findings show that an elliptic lumen shape and a thinner fibrous cap with an underlying lipid pool result in higher stenosis reduction, while large calcifications and fibrotic tissue are more prone to recoil. The shielding effect of a thicker fibrous cap helps to reduce local compressive stresses in the soft plaque. The presence of a soft plaque reduces the damage in the healthy vascular structures. Contrarily, the presence of hard plaque promotes less damage volume in the fibrous cap and reduces stress peaks in this region, but they seem to increase stresses in the media-intima layer. Finally the reliability of the achieved results was put into clinical perspective.

Filling the void: A coalescent numerical and experimental technique to determine aortic stent graft mechanics

The presented study details a combined experimental and computational method to assess and compare the mechanical behavior of the main body of 4 different stent graft designs. The mechanical response to a flat plate compression and radial crimping of the devices is derived and related to geometrical and material features of different stent designs. The finite element modeling procedure is used to complement the experimental results and conduct a solution sensitivity study. Finite element evaluations of the mechanical behavior match well with experimental findings and are used as a quantitative basis to discuss design characteristics of the different devices.

What if you stretch the IFU? A mechanical insight into stent graft Instructions For Use in angulated proximal aneurysm necks

Endovascular treatment for patients with a proximal neck anatomy outside instructions for use is an ongoing topic of debate in endovascular aneurysm repair. This paper employs the finite element method to offer insight into possible adverse effects of deploying a stent graft into an angulated geometry. The effect of angulation, straight neck length and device oversize was investigated in a full factorial parametric analysis. Stent apposition, area reduction of the graft, asymmetry of contact forces and the ability to find a good seal were investigated. Most adverse effects are expected for combinations of high angulation and short straight landing zones. Higher oversize has a beneficiary effect, but not enough to compensate the adverse effects of (very) short and angulated angles. Our analysis shows that for an angle between the suprarenal aorta and proximal neck above 60°, proximal kinking of the device can occur. The method used offers a engineering view on the morphological limits of EVAR for a clinically used device.

Assessment of shear stress related parameters in the carotid bifurcation using mouse-specific FSI simulations

The ApoE(-)(/)(-) mouse is a common small animal model to study atherosclerosis, an inflammatory disease of the large and medium sized arteries such as the carotid artery. It is generally accepted that the wall shear stress, induced by the blood flow, plays a key role in the onset of this disease. Wall shear stress, however, is difficult to derive from direct in vivo measurements, particularly in mice. In this study, we integrated in vivo imaging (micro-Computed Tomography-µCT and ultrasound) and fluid-structure interaction (FSI) modeling for the mouse-specific assessment of carotid hemodynamics and wall shear stress. Results were provided for 8 carotid bifurcations of 4 ApoE(-)(/)(-) mice. We demonstrated that accounting for the carotid elasticity leads to more realistic flow waveforms over the complete domain of the model due to volume buffering capacity in systole. The 8 simulated cases showed fairly consistent spatial distribution maps of time-averaged wall shear stress (TAWSS) and relative residence time (RRT). Zones with reduced TAWSS and elevated RRT, potential indicators of atherosclerosis-prone regions, were located mainly at the outer sinus of the external carotid artery. In contrast to human carotid hemodynamics, no flow recirculation could be observed in the carotid bifurcation region.

2-D Versus 3-D Cross-Correlation-Based Radial and Circumferential Strain Estimation Using Multiplane 2-D Ultrafast Ultrasound in a 3-D Atherosclerotic Carotid Artery Model

Three-dimensional (3-D) strain estimation might improve the detection and localization of high strain regions in the carotid artery (CA) for identification of vulnerable plaques. This paper compares 2-D versus 3-D displacement estimation in terms of radial and circumferential strain using simulated ultrasound (US) images of a patient-specific 3-D atherosclerotic CA model at the bifurcation embedded in surrounding tissue generated with ABAQUS software. Global longitudinal motion was superimposed to the model based on the literature data. A Philips L11-3 linear array transducer was simulated, which transmitted plane waves at three alternating angles at a pulse repetition rate of 10 kHz. Interframe (IF) radio-frequency US data were simulated in Field II for 191 equally spaced longitudinal positions of the internal CA. Accumulated radial and circumferential displacements were estimated using tracking of the IF displacements estimated by a two-step normalized cross-correlation method and displacement compounding. Least-squares strain estimation was performed to determine accumulated radial and circumferential strain. The performance of the 2-D and 3-D methods was compared by calculating the root-mean-squared error of the estimated strains with respect to the reference strains obtained from the model. More accurate strain images were obtained using the 3-D displacement estimation for the entire cardiac cycle. The 3-D technique clearly outperformed the 2-D technique in phases with high IF longitudinal motion. In fact, the large IF longitudinal motion rendered it impossible to accurately track the tissue and cumulate strains over the entire cardiac cycle with the 2-D technique.

Healthy and diseased coronary bifurcation geometries influence near-wall and intravascular flow: A computational exploration of the hemodynamic risk

Local hemodynamics has been identified as one main determinant in the onset and progression of atherosclerotic lesions at coronary bifurcations. Starting from the observation that atherosensitive hemodynamic conditions in arterial bifurcation are majorly determined by the underlying anatomy, the aim of the present study is to investigate how peculiar coronary bifurcation anatomical features influence near-wall and intravascular flow patterns. Different bifurcation angles and cardiac curvatures were varied in population-based, idealized models of both stenosed and unstenosed bifurcations, representing the left anterior descending (LAD) coronary artery with its diagonal branch. Local hemodynamics was analyzed in terms of helical flow and exposure to low/oscillatory shear stress by performing computational fluid dynamics simulations. Results show that bifurcation angle impacts lowly hemodynamics in both stenosed and unstenosed cases. Instead, curvature radius influences the generation and transport of helical flow structures, with smaller cardiac curvature radius associated to higher helicity intensity. Stenosed bifurcation models exhibit helicity intensity values one order of magnitude higher than the corresponding unstenosed cases. Cardiac curvature radius moderately affects near-wall hemodynamics of the stenosed cases, with smaller curvature radius leading to higher exposure to low shear stress and lower exposure to oscillatory shear stress. In conclusion, the proposed controlled benchmark allows investigating the effect of various geometrical features on local hemodynamics at the LAD/diagonal bifurcation, highlighting that cardiac curvature influences near wall and intravascular hemodynamics, while bifurcation angle has a minor effect.