A literature review of the numerical analysis of abdominal aortic aneurysms treated with endovascular stent grafts

Patient-Specific Numerical Simulations of Endovascular Procedures in Complex Aortic Pathologies: Review and Clinical Perspectives

The endovascular technique is used in the first line treatment in many complex aortic pathologies. Its clinical outcome is mostly determined by the appropriate selection of a stent-graft for a specific patient and the operator's experience. New tools are still needed to assist practitioners with decision making before and during procedures. For this purpose, numerical simulation enables the digital reproduction of an endovascular intervention with various degrees of accuracy. In this review, we introduce the basic principles and discuss the current literature regarding the use of numerical simulation for endovascular management of complex aortic diseases. Further, we give the future direction of everyday clinical applications, showing that numerical simulation is about to revolutionize how we plan and carry out endovascular interventions.

Volumetric lattice Boltzmann method for wall stresses of image-based pulsatile flows

Image-based computational fluid dynamics (CFD) has become a new capability for determining wall stresses of pulsatile flows. However, a computational platform that directly connects image information to pulsatile wall stresses is lacking. Prevailing methods rely on manual crafting of a hodgepodge of multidisciplinary software packages, which is usually laborious and error-prone. We present a new computational platform, to compute wall stresses in image-based pulsatile flows using the volumetric lattice Boltzmann method (VLBM). The novelty includes: (1) a unique image processing to extract flow domain and local wall normality, (2) a seamless connection between image extraction and VLBM, (3) an en-route calculation of strain-rate tensor, and (4) GPU acceleration (not included here). We first generalize the streaming operation in the VLBM and then conduct application studies to demonstrate its reliability and applicability. A benchmark study is for laminar and turbulent pulsatile flows in an image-based pipe (Reynolds number: 10 to 5000). The computed pulsatile velocity and shear stress are in good agreements with Womersley's analytical solutions for laminar pulsatile flows and concurrent laboratory measurements for turbulent pulsatile flows. An application study is to quantify the pulsatile hemodynamics in image-based human vertebral and carotid arteries including velocity vector, pressure, and wall-shear stress. The computed velocity vector fields are in reasonably well agreement with MRA (magnetic resonance angiography) measured ones. This computational platform is good for image-based CFD with medical applications and pore-scale porous media flows in various natural and engineering systems.

Analytical methods for braided stents design and comparison with FEA

Braiding technology is nowadays commonly adopted to build stent-like devices. Indeed, these endoprostheses, thanks to their typical great flexibility and kinking resistance, find several applications in mini-invasive treatments, involving but not limiting to the cardiovascular field. The design process usually involves many efforts and long trial and error processes before identifying the best combination of manufacturing parameters. This paper aims to provide analytical tools to support the design and optimization phases: the developed equations, based on few geometrical parameters commonly used for describing braided stents and material stiffness, are easily implementable in a worksheet and allow predicting the radial rigidity of braided stents, also involving complex features such as multiple twists and looped ends, and the diameter variation range. Finite element simulations, previously validated with respect to experimental tests, were used as a comparator to prove the reliability of the analytical results. The illustrated tools can assess the impact of each selected parameter modification and are intended to guide the optimal selection of geometrical and mechanical stent proprieties to obtain the desired radial rigidity, deliverability (minimum diameter), and, if forming processes are planned to modify the shape of the stent, the required diameter variations (maximum and minimum diameters).

Impact of Thoracic Endografting on the Hemodynamics of the Native Aorta: Pre- and Postoperative Assessments of Wall Shear Stress and Vorticity Using Computational Fluid Dynamics

PURPOSE

To quantify the hemodynamic consequences of thoracic endovascular aortic repair (TEVAR) by comparing the preoperative and postoperative wall shear stress (WSS) and vorticity profiles on computational fluid dynamics (CFD) simulations.

MATERIALS AND METHODS

The pre- and postoperative computed tomography (CT) scans from 20 consecutive patients (median age 69 years, range 20-87) treated for different thoracic aortic pathologies (11 aneurysms, 5 false aneurysms, 3 penetrating ulcers, and 1 traumatic aortic rupture) were segmented to construct patient-specific CFD models using a meshless code. The simulations were run over the cardiac cycle, and the WSS and vorticity values measured at the proximal and distal landing zones were compared.

RESULTS

The CFD runs provided 4-dimensional simulations of blood flow in all patients. WSS and vorticity profiles at the proximal landing zone (located in zones 0-3 in 15 patients) varied in 18 and 20 of the cases, respectively; WSS was increased in 11 cases and the vorticity in 9. Pre- and postoperative WSS median values were 4.19 and 4.90 Pa, respectively. Vorticity median values were 40.38 and 39.17 Hz, respectively.

CONCLUSION

TEVAR induces functional alterations in the native thoracic aorta, though the prognostic significance of these changes is still unknown. CFD appears to be a valuable tool to explore aortic hemodynamics, and its application in a larger series would help define a predictive role for these hemodynamic assessments.

Anthropomorphic and biomechanical mockup for abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) is an asymptomatic condition due to the dilation of abdominal aorta along with progressive wall degeneration, where rupture of AAA is life-threatening. Failures of AAA endovascular repair (EVAR) reflect our inadequate knowledge about the complex interaction between the aortic wall and medical devices. In this regard, we are presenting a hydrogel-based anthropomorphic mockup (AMM) to better understand the biomechanical constraints during EVAR. By adjusting the cryogenic treatments, we tailored the hydrogel to mimic the mechanical behavior of human AAA wall, thrombus and abdominal fat. A specific molding sequence and a pressurizing system were designed to reproduce the geometrical and diseased characteristics of AAA. A mechanically, anatomically and pathologically realistic AMM for AAA was developed for the first time, EVAR experiments were then performed with and without the surrounding fat. Substantial displacements of the aortic centerlines and vessel expansion were observed in the case without surrounding fat, revealing an essential framework created by the surrounding fat to account for the interactions with medical devices. In conclusion, the importance to consider surrounding tissue for the global deformation of AAA during EVAR was highlighted. Furthermore, potential use of this AMM for medical training was also suggested.

Advanced Radial Basis Functions Mesh Morphing for High Fidelity Fluid-Structure Interaction with Known Movement of the Walls: Simulation of an Aortic Valve

High fidelity Fluid-Structure Interaction (FSI) can be tackled by means of non-linear Finite Element Models (FEM) suitable to capture large deflections of structural parts interacting with fluids and by means of detailed Computational Fluid Dynamics (CFD). High fidelity is gained thanks to the spatial resolution of the computational grids and a key enabler to have a proper exchange of information between the structural solver and the fluid one is the management of the interfaces. A class of applications consists in problems where the complex movement of the walls is known in advance or can be computed by FEM and has to be transferred to the CFD solver. The aforementioned approach, known also as one-way FSI, requires effective methods for the time marching adaption of the computation grid of the CFD model. A versatile and well established approach consists in a continuum update of the mesh that is regenerated so to fit the evolution of the moving walls. In this study, an innovative method based on Radial Basis Functions (RBF) mesh morphing is proposed, allowing to keep the same mesh topology suitable for a continuum update of the shape. A set of key configurations are exactly guaranteed whilst time interpolation is adopted between frames. The new framework is detailed and then demonstrated, adopting as a reference the established approach based on remeshing, for the study of a Polymeric-Prosthetic Heart Valve (P-PHV).

Fatigue lifetime prediction of a reduced-diameter dental implant system: Numerical and experimental study

OBJECTIVE

To validate the fatigue lifetime of a reduced-diameter dental implant system predicted by three-dimensional finite element analysis (FEA) by testing physical implant specimens using an accelerated lifetime testing (ALT) strategy with the apparatus specified by ISO 14801.

METHODS

A commercially-available reduced-diameter titanium dental implant system (Straumann Standard Plus NN) was digitized using a micro-CT scanner. Axial slices were processed using an interactive medical image processing software (Mimics) to create 3D models. FEA analysis was performed in ABAQUS, and fatigue lifetime was predicted using fe-safe^® software. The same implant specimens (n=15) were tested at a frequency of 2Hz on load frames using apparatus specified by ISO 14801 and ALT. Multiple step-stress load profiles with various aggressiveness were used to improve testing efficiency. Fatigue lifetime statistics of physical specimens were estimated in a reliability analysis software (ALTA PRO). Fractured specimens were examined using SEM with fractographic technique to determine the failure mode.

RESULTS

FEA predicted lifetime was within the 95% confidence interval of lifetime estimated by experimental results, which suggested that FEA prediction was accurate for this implant system. The highest probability of failure was located at the root of the implant body screw thread adjacent to the simulated bone level, which also agreed with the failure origin in physical specimens.

SIGNIFICANCE

Fatigue lifetime predictions based on finite element modeling could yield similar results in lieu of physical testing, allowing the use of virtual testing in the early stages of future research projects on implant fatigue.

A new approach for the pre-clinical optimization of a spatial configuration of bifurcated endovascular prosthesis placed in abdominal aortic aneurysms

Complexity of the spatial configuration of an aortic implant with bifurcation in the distal part is related to changes in blood hemodynamic in the area of bifurcation which may disturb blood flow and lead to thrombus formation. This study was designed to characterize parameters which define spatial configuration of an aortic implant for which the risk of thrombus formation is the smallest. We used AngioCT data from 74 patients, aged 55 ±10 years, after endovascular procedure to prepare 3D geometries of stent-grafts. Computational Fluid Dynamics (CFD) simulations were used to reconstruct blood hemodynamic and simulate thrombus formation. Next, geometric parameters of stent-grafts included the ratio of volume of upper part to the bifurcations, the relation of inlet and outlet diameters of a stent-graft and deformations in the iliac part of the stent-graft were analyzed. We also analyzed tortuosities (spiral twisting of the flow around the flow direction) and bends (the largest angulation in distal part of a stent-graft). The CFD results were confronted with AngioCT data to verify if computer generated thrombus appeared in particular patient. Additionally, geometric parameters of analyzed stent-grafts were used to propose a mathematical tool for prediction of thrombus appearance. The results showed that tortuosities and bends of a stent-graft had the highest impact on thrombus formation. Formation of thrombi was observed in 22% to 31% of cases (at blood hematocrit Hct = 40%) even for small values of tortuosities and bends indicating that these parameters are dominant in determining blood clotting. Our calculated results overlapped with clinical data in 80% to 91%. Therefore, we conclude that tortuosities and bends have high impact on thrombus formation and should be under special attention during stent-graft recommendation and patients’ follow-ups.

Hemodynamic Profile of Two Aortic Endografts Accounting for Their Postimplantation Position

Endovascular aneurysm repair (EVAR) is a clinically effective technique for treating anatomically eligible abdominal aortic aneurysms (AAAs), involving the deployment of an endograft (EG) that is designed to prevent blood leakage in the aneurysmal sac. While most EGs have equivalent operating principles, the hemodynamic environment established by different EGs is not necessarily the same. So, to unveil the post-EVAR hemodynamic properties, we need an EG-specific computational approach that currently lacks from the literature. Endurant and Excluder are two EGs with similar pre-installation designs. We assumed that the flow conditions in the particular EGs do not vary significantly. The hypothesis was tested combining image reconstructions, computational fluid dynamics (CFD), and statistics, taking into account the postimplantation position of the EGs. Ten patients with Endurant EGs and ten patients with Excluder EGs were included in this study. The two groups were matched with respect to the preoperative morphological characteristics of the AAAs. The EG models are derived from image reconstructions of postoperative computed tomography scans. Wall shear stress (WSS), displacement force, velocity, and helicity were calculated in regions of interest within the EG structures, i.e., the main body, the upper and lower part of the limbs. Excluder generated higher WSS compared to Endurant, especially on the lower part of the limbs (p = 0.001). Spatial fluctuations of WSS were observed on the upper part of the Excluder limbs. Higher blood velocity was induced by Excluder in all the regions of interest (p = 0.04, p = 0.01, and p = 0.004). Focal points of secondary flow were detected in the main body of Endurant and the limbs of Excluder. The displacement force acting on the lower part of the Excluder limbs was stronger compared to the Endurant one (p = 0.03). The results showed that two similar EGs implanted in similar AAAs can induce significantly different flow properties. The delineation of the hemodynamic features associated with the various commercially available EGs could further promote the personalization of treatment offered to aneurysmal patients and inspire ideas for the improvement of EG designs in the future.

A 3D model of thrombus formation in a stent-graft after implantation in the abdominal aorta

Here we present a 3D kinetic model of thrombus formation in an endovascular prosthesis after implantation in an abdominal aorta with an aneurysm. The computational fluid dynamic technique (CFD) was used to determine the process of thrombus formation and growth in the stent-graft on the basis of the medical data from computed tomography angiography and Doppler ultrasound examination of 10 patients. The Quemada model was used to describe rheological properties of blood. Results of the CFD simulations were validated based on actual data from patients with diagnosed thrombi in aortic implants. The results show that the elaborated CFD model correctly predicted thrombus formation, shape and deposition site in an endovascular prosthesis. The developed CFD model of thrombus growth can be applied to predict the risk of thrombus formation in stent-grafts and assist in selection of geometry of the endovascular prosthesis to reduce possible complications after stent-graft implantation using only basic medical data.

Improving the efficiency of abdominal aortic aneurysm wall stress computations

An abdominal aortic aneurysm is a pathological dilation of the abdominal aorta, which carries a high mortality rate if ruptured. The most commonly used surrogate marker of rupture risk is the maximal transverse diameter of the aneurysm. More recent studies suggest that wall stress from models of patient-specific aneurysm geometries extracted, for instance, from computed tomography images may be a more accurate predictor of rupture risk and an important factor in AAA size progression. However, quantification of wall stress is typically computationally intensive and time-consuming, mainly due to the nonlinear mechanical behavior of the abdominal aortic aneurysm walls. These difficulties have limited the potential of computational models in clinical practice. To facilitate computation of wall stresses, we propose to use a linear approach that ensures equilibrium of wall stresses in the aneurysms. This proposed linear model approach is easy to implement and eliminates the burden of nonlinear computations. To assess the accuracy of our proposed approach to compute wall stresses, results from idealized and patient-specific model simulations were compared to those obtained using conventional approaches and to those of a hypothetical, reference abdominal aortic aneurysm model. For the reference model, wall mechanical properties and the initial unloaded and unstressed configuration were assumed to be known, and the resulting wall stresses were used as reference for comparison. Our proposed linear approach accurately approximates wall stresses for varying model geometries and wall material properties. Our findings suggest that the proposed linear approach could be used as an effective, efficient, easy-to-use clinical tool to estimate patient-specific wall stresses.

Morphologic evaluation of ruptured and symptomatic abdominal aortic aneurysm by three-dimensional modeling

OBJECTIVE

To identify geometric indices of abdominal aortic aneurysms (AAAs) on computed tomography that are associated with higher risk of rupture.

METHODS

This retrospective case-control, institutional review board-approved study involved 63 cases with ruptured or symptomatic AAA and 94 controls with asymptomatic AAA. Three-dimensional models were generated from computed tomography segmentation and used for the calculation of 27 geometric indices. On the basis of the results of univariate analysis and multivariable sequential logistic regression analyses with a forward stepwise model selection based on likelihood ratios, a traditional model based on gender and maximal diameter (Dmax) was compared with a model that also incorporated geometric indices while adjusting for gender and Dmax. Receiver operating characteristic (ROC) curves were calculated for these two models to evaluate their classification accuracy.

RESULTS

Univariate analysis revealed that gender (P = .024), Dmax (P = .001), and 14 other geometric indices were associated with AAA rupture at P < .05. In the multivariable analysis, adjusting for gender and Dmax, the AAA with a higher bulge location (P = .020) and lower mean averaged area (P = .005) were associated with AAA rupture. With these two geometric indices, the area under the ROC curve showed an improvement from 0.67 (95% confidence interval, 0.58-0.77) to 0.75 (95% confidence interval, 0.67-0.83; P < .001). Our predictive model showed comparable sensitivity (64% vs 60%) and specificity (79% vs 77%) with current treatment criteria based on gender and diameter at the point optimizing the Youden index (sensitivity + specificity - 1) on the ROC curve.

CONCLUSIONS

Two geometric indices derived from AAA three-dimensional modeling were independently associated with AAA rupture. The addition of these indices in a predictive model based on current treatment criteria modestly improved the accuracy to detect aneurysm rupture.