Ender A. Finol

Fluid-structure interaction in abdominal aortic aneurysms: effects of asymmetry and wall thickness

BackgroundAbdominal aortic aneurysm (AAA) is a prevalent disease which is of significant concern because of the morbidity associated with the continuing expansion of the abdominal aorta and its ultimate rupture. The transient interaction between blood flow and the wall contributes to wall stress which, if it exceeds the failure strength of the dilated arterial wall, will lead to aneurysm rupture. Utilizing a computational approach, the biomechanical environment of virtual AAAs can be evaluated to study the affects of asymmetry and wall thickness on this stress, two parameters that contribute to increased risk of aneurysm rupture.MethodsTen virtual aneurysm models were created with five different asymmetry parameters ranging from β = 0.2 to 1.0 and either a uniform or variable wall thickness to study the flow and wall dynamics by means of fully coupled fluid-structure interaction (FSI) analyses. The AAA wall was designed to have a (i) uniform 1.5 mm thickness or (ii) variable thickness ranging from 0.5 – 1.5 mm extruded normally from the boundary surface of the lumen. These models were meshed with linear hexahedral elements, imported into a commercial finite element code and analyzed under transient flow conditions. The method proposed was then compared with traditional computational solid stress techniques on the basis of peak wall stress predictions and cost of computational effort.ResultsThe results provide quantitative predictions of flow patterns and wall mechanics as well as the effects of aneurysm asymmetry and wall thickness heterogeneity on the estimation of peak wall stress. These parameters affect the magnitude and distribution of Von Mises stresses; varying wall thickness increases the maximum Von Mises stress by 4 times its uniform thickness counterpart. A pre-peak systole retrograde flow was observed in the AAA sac for all models, which is due to the elastic energy stored in the compliant arterial wall and the expansion force of the artery during systole.ConclusionBoth wall thickness and geometry asymmetry affect the stress exhibited by a virtual AAA. Our results suggest that an asymmetric AAA with regional variations in wall thickness would be exposed to higher mechanical stresses and an increased risk of rupture than a more fusiform AAA with uniform wall thickness. Therefore, it is important to accurately reproduce vessel geometry and wall thickness in computational predictions of AAA biomechanics.

Wall stress and flow dynamics in abdominal aortic aneurysms: finite element analysis vs. fluid–structure interaction

Abdominal aortic aneurysm (AAA) rupture is the clinical manifestation of an induced force exceeding the resistance provided by the strength of the arterial wall. This force is most frequently assumed to be the product of a uniform luminal pressure acting along the diseased wall. However fluid dynamics is a known contributor to the pathogenesis of AAAs, and the dynamic interaction of blood flow and the arterial wall represents the in vivo environment at the macro-scale. The primary objective of this investigation is to assess the significance of assuming an arbitrary estimated peak fluid pressure inside the aneurysm sac for the evaluation of AAA wall mechanics, as compared with the non-uniform pressure resulting from a coupled fluid–structure interaction (FSI) analysis. In addition, a finite element approach is utilised to estimate the effects of asymmetry and wall thickness on the wall stress and fluid dynamics of ten idealised AAA models and one non-aneurysmal control. Five degrees of asymmetry with uniform and variable wall thickness are used. Each was modelled under a static pressure-deformation analysis, as well as a transient FSI. The results show that the inclusion of fluid flow yields a maximum AAA wall stress up to 20% higher compared to that obtained with a static wall stress analysis with an assumed peak luminal pressure of 117 mmHg. The variable wall models have a maximum wall stress nearly four times that of a uniform wall thickness, and also increasing with asymmetry in both instances. The inclusion of an axial stretch and external pressure to the computational domain decreases the wall stress by 17%.

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patients average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.

Blood Flow Dynamics in Saccular Aneurysm Models of the Basilar Artery

Blood flow dynamics under physiologically realistic pulsatile conditions plays an important role in the growth, rupture, and surgical treatment of intracranial aneurysms. The temporal and spatial variations of wall pressure and wall shear stress in the aneurysm are hypothesized to be correlated with its continuous expansion and eventual rupture. In addition, the assessment of the velocity field in the aneurysm dome and neck is important for the correct placement of endovascular coils. This paper describes the flow dynamics in two representative models of a terminal aneurysm of the basilar artery under Newtonian and non-Newtonian fluid assumptions, and compares their hemodynamics with that of a healthy basilar artery. Virtual aneurysm models are investigated numerically, with geometric features defined by beta = 0 deg and beta = 23.2 deg, where beta is the tilt angle of the aneurysm dome with respect to the basilar artery. The intra-aneurysmal pulsatile flow shows complex ring vortex structures for beta = 0 deg and single recirculation regions for beta = 23.2 deg during both systole and diastole. The pressure and shear stress on the aneurysm wall exhibit large temporal and spatial variations for both models. When compared to a non-Newtonian fluid, the symmetric aneurysm model (beta = 0 deg) exhibits a more unstable Newtonian flow dynamics, although with a lower peak wall shear stress than the asymmetric model (beta = 23.2 deg). The non-Newtonian fluid assumption yields more stable flows than a Newtonian fluid, for the same inlet flow rate. Both fluid modeling assumptions, however, lead to asymmetric oscillatory flows inside the aneurysm dome.

Three-Dimensional Geometrical Characterization of Abdominal Aortic Aneurysms: Image-Based Wall Thickness Distribution

The clinical assessment of abdominal aortic aneurysm (AAA) rupture risk is based on the quantification of AAA size by measuring its maximum diameter from computed tomography (CT) images and estimating the expansion rate of the aneurysm sac over time. Recent findings have shown that geometrical shape and size, as well as local wall thickness may be related to this risk; thus, reliable noninvasive image-based methods to evaluate AAA geometry have a potential to become valuable clinical tools. Utilizing existing CT data, the three-dimensional geometry of nine unruptured human AAAs was reconstructed and characterized quantitatively. We propose and evaluate a series of 1D size, 2D shape, 3D size, 3D shape, and second-order curvature-based indices to quantify AAA geometry, as well as the geometry of a size-matched idealized fusiform aneurysm and a patient-specific normal abdominal aorta used as controls. The wall thickness estimation algorithm, validated in our previous work, is tested against discrete point measurements taken from a cadaver tissue model, yielding an average relative difference in AAA wall thickness of 7.8%. It is unlikely that any one of the proposed geometrical indices alone would be a reliable index of rupture risk or a threshold for elective repair. Rather, the complete geometry and a positive correlation of a set of indices should be considered to assess the potential for rupture. With this quantitative parameter assessment, future research can be directed toward statistical analyses correlating the numerical values of these parameters with the risk of aneurysm rupture or intervention (surgical or endovascular). While this work does not provide direct insight into the possible clinical use of the geometric parameters, we believe it provides the foundation necessary for future efforts in that direction.

Flow-induced Wall Shear Stress in Abdominal Aortic Aneurysms: Part I - Steady Flow Hemodynamics

Numerical predictions of blood flow patterns and hemodynamic stresses in Abdominal Aortic Aneurysms (AAAs) are performed in a two-aneurysm, axisymmetric, rigid wall model using the spectral element method. Homogeneous, Newtonian blood flow is simulated under steady conditions for the range of Reynolds numbers 10 h Re h 2265. Flow hemodynamics are quantified by calculating the distributions of wall pressure (p w ), wall shear stress ( w ), Wall Shear Stress Gradient (WSSG). A correlation between maximum values of hemodynamic stresses and Reynolds number is established, and the spatial distribution of WSSG is considered as a hemodynamic force that may cause damage to the arterial wall at an intermediate stage of AAA growth. The temporal distribution of hemodynamic stresses in pulsatile flow and their physical implications in AAA rupture are discussed in Part II of this paper.

The Effect of Material Model Formulation in the Stress Analysis of Abdominal Aortic Aneurysms

A reliable estimation of wall stress in Abdominal Aortic Aneurysms (AAAs), requires performing an accurate three-dimensional reconstruction of the medical image-based native geometry and modeling an appropriate constitutive law for the aneurysmal tissue material characterization. A recent study on the biaxial mechanical behavior of human AAA tissue specimens demonstrates that aneurysmal tissue behaves mechanically anisotropic. Results shown in this communication show that the peak wall stress is highly sensitive to the anisotropic model used for the stress analysis. In addition, the present investigation indicates that structural parameters (e.g., collagen fiber orientation) should be determined independently and not by means of non-linear fitting to stress–strain test data. Fiber orientation identified in this manner could lead to overestimated peak wall stresses.

A Relation Between Near-Wall Particle-Hemodynamics and Onset of Thrombus Formation in Abdominal Aortic Aneurysms

A novel computational particle-hemodynamics analysis of key criteria for the onset of an intraluminal thrombus (ILT) in a patient-specific abdominal aortic aneurysm (AAA) is presented. The focus is on enhanced platelet and white blood cell residence times as well as their elevated surface-shear loads in near-wall regions of the AAA sac. The generalized results support the hypothesis that a patient’s AAA geometry and associated particle-hemodynamics have the potential to entrap activated blood particles, which will play a role in the onset of ILT. Although the ILT history of only a single patient was considered, the modeling and simulation methodology provided allow for the development of an efficient computational tool to predict the onset of ILT formation in complex patient-specific cases.

Wall apposition assessment and performance comparison of distal protection filters

Purpose: To assess the wall apposition of 3 distal protection filters used in carotid artery stenting (CAS) for cerebral protection and quantify the effect on the in vitro capture efficiency of the filters under simulated physiological flow conditions. Methods: The 3 distal protection filters (Angioguard XP, FilterWire EZ, and RX Accunet) were deployed in silicone flow models of 5.0-, 5.5-, and 6.0-mm inner diameter and were tested at a mean flow rate characteristic of the human internal carotid artery while injecting polydispersed microspheres simulating plaque emboli. The injected microspheres had a diameter larger than the pore size of the devices tested, so it was conjectured that any microspheres missed by the device traveled between the device basket and the vessel wall. Varying the diameter of the vessel phantom within the recommended vessel diameter treatment range for each device simulated the variability of vessel diameter in vivo, allowing the quantification of device wall apposition. Results: None of the devices tested completely prevented distal embolization in the flow model. The RX Accunet device has the best overall wall apposition, yielding gaps of 0.075% of the vessel cross-sectional area. The FilterWire EZ device had the best overall average filtration rate, failing to capture only 0.8% of plaque particles. There were no statistically significant differences in the wall apposition assessment or the capture efficiency of the RX Accunet and FilterWire EZ devices. Conclusion: Several complications related to apposition of the filter basket on the vessel wall and device retrieval were detected in all the devices. It is inferred that the adaptability of the filter basket to conform to the vessel cross section at the site of deployment is the primary design variable responsible for distal embolization during CAS with cerebral protection.

FSI Analysis of a Healthy and a Stenotic Human Trachea Under Impedance-Based Boundary Conditions

In this work, a fluid-solid interaction (FSI) analysis of a healthy and a stenotic human trachea was studied to evaluate flow patterns, wall stresses, and deformations under physiological and pathological conditions. The two analyzed tracheal geometries, which include the first bifurcation after the carina, were obtained from computed tomography images of healthy and diseased patients, respectively. A finite element-based commercial software code was used to perform the simulations. The tracheal wall was modeled as a fiber reinforced hyperelastic solid material in which the anisotropy due to the orientation of the fibers was introduced. Impedance-based pressure waveforms were computed using a method developed for the cardiovascular system, where the resistance of the respiratory system was calculated taking into account the entire bronchial tree, modeled as binary fractal network. Intratracheal flow patterns and tracheal wall deformation were analyzed under different scenarios. The simulations show the possibility of predicting, with FSI computations, flow and wall behavior for healthy and pathological tracheas. The computational modeling procedure presented herein can be a useful tool capable of evaluating quantities that cannot be assessed in vivo, such as wall stresses, pressure drop, and flow patterns, and to derive parameters that could help clinical decisions and improve surgical outcomes.