Biomechanical Investigation of Disturbed Hemodynamics-Induced Tissue Degeneration in Abdominal Aortic Aneurysms Using Computational and Experimental Techniques.

Abdominal aortic aneurysm (AAA) is the dilatation of the aorta beyond 50% of the normal vessel diameter. It is reported that 4-8% of men and 0.5-1% of women above 50 years of age bear an AAA and it accounts for ~15,000 deaths per year in the United States alone. If left untreated, AAA might gradually expand until rupture; the most catastrophic complication of the aneurysmal disease that is accompanied by a striking overall mortality of 80%. The precise mechanisms leading to AAA rupture remains unclear. Therefore, characterization of disturbed hemodynamics within AAAs will help to understand the mechanobiological development of the condition which will contribute to novel therapies for the condition. Due to geometrical complexities, it is challenging to directly quantify disturbed flows for AAAs clinically. Two other approaches for this investigation are computational modeling and experimental flow measurement. In computational modeling, the problem is first defined mathematically, and the solution is approximated with numerical techniques to get characteristics of flow. In experimental flow measurement, once the setup providing physiological flow pattern in a phantom geometry is constructed, velocity measurement system such as particle image velocimetry (PIV) enables characterization of the flow. We witness increasing number of applications of these complimentary approaches for AAA investigations in recent years. In this paper, we outline the details of computational modeling procedures and experimental settings and summarize important findings from recent studies, which will help researchers for AAA investigations and rupture mechanics.

Introduction

Abdominal aortic aneurysm (AAA) is the dilatation of the abdominal aorta beyond 50% of the normal vessel diameter, due to degeneration of the arterial wall (McGloughlin Timothy and Doyle Barry, 2010). It is reported that 4–8% of men and 0.5–1% of women above 50 years of age bear an AAA and it accounts for ~15,000 deaths per year in the United States alone (Sakalihasan et al., 2005; Kontopodis et al., 2014). AAA may result in rupture of the vessel wall, which is a fatal surgical emergency because of reduced blood flow to vital organs and hematocele (i.e., swelling caused by blood collecting in a body cavity). AAA rupture is the most catastrophic complication of the aneurysmal disease that is accompanied by a striking overall mortality of 80% (Bengtsson and Bergqvist, 1993). Unfortunately, clinical symptoms for impending AAA rupture are not observed for nearly 75% of the patients (Anjum et al., 2012).

Current clinical practice is surgical repair if the diameter of the AAA is >5.5 cm or the growth rate is over 1 cm/year and to follow up AAAs with diameters < 5.5 cm at 6-month intervals (Cosford and Leng, 2007). It is reported that only 25% of AAAs ruptured in the patient's lifetime (Darling et al., 1977), and 13% of AAAs with a diameter of < 5 cm are ruptured. For this reason, AAA diameter is not the sole indicator for the rupture and other factors should also be considered for the rupture risk assessment, together with the size of AAA.

The precise mechanisms leading to AAA rupture remains unclear. From a purely mechanical point of view, rupture of an AAA occurs when the mechanical stresses (i.e., internal forces per unit area) acting on the aneurysm exceeds the ability of the wall tissue to withstand these stresses (i.e., the wall's failure strength). Wall shear stress (i.e., WSS, the frictional force exerted by blood flow on the luminal surface) is also thought to play an important role (Peattie et al., 1996). While the blood flow in normal aorta (diameter within 2–2.5 cm) is mainly anterograde with high WSS, during AAA, circulatory flows emerge within the vessel. Hence, flow becomes disturbed with oscillatory characteristics leading to low WSS (Tanweer et al., 2014). The disturbance of flow in AAA is believed to contribute to the progression of the disease by activating the inflammatory markers of the endothelial cells lining the vessel wall which might lead to degeneration and weakening of the vessel wall (Franck et al., 2013). Therefore, prediction of the growth course of an AAA and precise determination of the rupture risk is a very challenging problem and attracts substantial attention.

Currently, there is no established conclusive approach for AAA rupture risk assessment, rather than the limited clinical guidelines for the AAA size. Ideally, the problem needs to be investigated from a biomechanical point of view, for a comprehensive wall stress, wall strength, and hemodynamic analyses. Computational modeling based on finite element analysis (FEA) is a powerful technique for estimating mechanical behavior of materials (i.e., mechanical behavior is the amount of stresses and deformations in a substance under external forces). This is especially beneficial in medical research where in most cases direct measurement of in-vivo mechanical response of diseased tissue is not possible. In this technique, first an accurate representation of the problem geometry is generated. The second step is defining the conditions in the model boundaries. The third step is defining different zones in the geometry, and assigning appropriate material properties for each zone and the fourth step is dividing zones into smaller elements (meshing). The last step is numerically solving appropriate governing equations in each small element to obtain a solution field of interest in the entire geometry. Each of these steps are critically important for solution accuracy.

Computational modeling has been used frequently in recent years to investigate rupture mechanics for AAA tissue. Comparison of ruptured and electively repaired AAA cases via FEA, revealed elevated wall stress levels for ruptured tissue (Fillinger et al., 2002). Several observations on excised tissue suggest that AAA formation is accompanied by an increase in wall stress as well as a corresponding decrease in wall strength (Vorp et al., 1996; Raghavan et al., 2000) and rupture point on AAA is usually coincides with peak wall stress locations (Fillinger et al., 2002). In very few studies where previous medical images for ruptured AAA tissue were available, patient-specific FEA could successfully identify known exact future rupture locations as high wall stress regions (Doyle et al., 2014). The FEA results prove the potential significant contribution of computational analysis for rupture risk assessment for AAA.

As computational methods have widely been used in investigation of hemodynamics and mechanical behavior of arterial tissue, the experimental techniques are also utilized in characterization of flow dynamics through AAAs. Both approaches are crucial and complement each other with offering in depth analysis where the level of its intensity and accuracy depend on the assumptions in computational models and simplifications in experimental methods. In line with the aforementioned computational studies, many experimental investigations of the hemodynamics through AAAs have been conducted in literature, and various qualitative and quantitative flow measurement techniques have been utilized for that purpose. A typical experimental set up for the analysis of hemodynamics contains physiological flow circulatory system including pump, piping, and pressure compliance, test section that contains artery model (phantoms), blood mimicking fluid, and flow measurement systems (like particle image velocimetry) to track the movement of fluid particles to calculate the velocities first, and then WSS levels on the phantoms. Previously, this approach has been adapted by several researchers to investigate AAA rupture mechanism (Tanweer et al., 2014; Wang et al., 2016).

In this paper, we will summarize the techniques utilized in experimental and computational studies of AAAs. We will briefly explain important aspects of the modeling procedure and how to set up and run relevant simulations and also, how to generate a flow circulation set up with particle image velocimetry (PIV). In addition, we will summarize important findings for experimental and computational biomechanical assessments of AAAs.

Computational Investigation of Biomechanics of AAAs

Computational models enable researchers to approximate biomechanical behavior of tissue stress and blood flow hemodynamics under realistic conditions. The models are generated by defining the conditions on the boundaries and numerically solving governing equations in fluid and solid domains. Generation of the problem geometry, setting up solid and fluid models, coupling the solutions through FSI procedure are important steps which contribute to enhance the reliability and accuracy of the biomechanical assessment as will be mentioned in the forthcoming sections.

Mathematical Details of the Problem

Geometry of the Problem

AAA has unique, patient-specific and complex geometry with a wall thickness around 1.5 mm (Raghavan et al., 2004) and mostly does not represent an axisymmetric form. Healthy abdominal aortic diameter is around 2 cm and in case of dilatation, AAA diameter can expand up to 9 cm. Proximal AAA neck angle, distal iliac bifurcation angle, heterogeneous wall thickness are patient-specific geometric parameters. Left and right renal arteries, left and right iliac arteries and superior mesenteric artery are the main branching arteries in abdominal aorta. In Figure 1, different patient-specific AAA geometries, branching arteries and disturbed flow in AAA sac are presented.

Figure 1

Six different patient-specific medical images and corresponding reconstructed AAA geometries. The branching arteries and disturbed AAA flow in the 3rd patient-specific medical image are described at the right side. Compared to undilated vessel, the flow hemodynamics in the AAA sac changed significantly. The artery labels are as follows: H, Hepatic artery; S, Splenic artery; RR, Right renal artery; ARR, Accessory right renal artery; SM, Superior mesenteric artery; LR, Left renal artery; AAA, Abdominal aortic aneurysm; RI, Right iliac artery; LI, Left iliac artery; REI, Right external iliac artery; RII, Right internal iliac artery; LII, Left internal iliac artery; LEI, Left external iliac artery [The figure is adapted from Les et al. (2010) and used with permission].

Main geometric parameters for a patient-specific AAA are presented in Figure 2. The AAA asymmetry is defined using the relation defined by Vorp et al. (1998). β is the parameter of asymmetry defined according to the central axis of undilated diameter, as shown in Figure 2. In Equation (1), r and R are the radii measured from center of the undilated portion to the posterior and anterior walls, respectively.

Figure 2

Main geometric parameters for a patient-specific AAA geometry. LAAA, Aneurysm length; DAMAX, Maximum aneurysm diameter; dproximal_neck, Inlet diameter of AAA sac; ddistal_neck, Outlet diameter of AAA sac; dabdominal_aorta, Normal abdominal aorta diameter; θ, Proximal neck angle; ϕ, Iliac bifurcation angle; L, Absolute length of tortuous vessel; τ, Imaginary straight line starting at the center of normal abdominal aorta and ending at the iliac bifurcation. The midsection at the location of the maximum AAA diameter (DAMAX) is represented at the right side. At the maximum diameter midsection, r and R are the radii measured from center of the undilated portion (i.e., normal abdominal aorta) to the posterior and anterior walls, respectively [The figure is adapted from Soudah et al. (2013) and used with permission].

In addition to the maximum aneurysm diameter, tortuosity, curvature, proximal neck angle, iliac bifurcation angle, and complex shape of the aneurysm itself are influencing factors on the wall stress (Stringfellow et al., 1987; Vorp et al., 1998; Hua and Mower, 2001). Up to 80% of AAA rupture are observed on the posterior wall (Darling et al., 1977), demonstrating the importance of asymmetry and curved AAA geometry in rupture. Loss of curvature on the posterior wall and AAA asymmetry are suggested as dominant wall stress increasing factors (Scotti et al., 2008).

Generation of Model Geometry

Magnetic resonance imaging (MRI), three-dimensional (3D) ultrasound, computerized tomography (CT) are the most common medical imaging techniques used to extract the realistic geometry of AAAs. At present, primary imaging method is CT. However, due to its high cost, significant radiation dose, injection of ionated contrast medium associated with nephrotoxity and relatively less availability, alternative methods such as 3D ultrasound can be used to widespread the biomechanical assessments.

After obtaining patient-specific medical images in DICOM (Digital Imaging and Communications in Medicine) format, 3D AAA models can be reconstructed using segmentation software such as MIMICS (Materialize, Leuven, Belgium), VESSEG (Carniege Mellon University, Pittsburgh, PA), ImFusion Suite (ImFusion GmbH, Munich, Germany), and open-source software such as SimVascular, VMTK, and ITK-SNAP. Lumen and AAA wall are segmented separately to differentiate these zones. Segmented 3D model might include protrusions and tight internal corners, which requires application of smoothing algorithms, to prepare a suitable geometry for FEA simulations. Autodesk MeshMixer, AngioLab, and MeshLab softwares can be used for further smoothing and mesh optimization processes. For smoothing process, different algorithms can be applied such as Laplacian (Field, 1988), HC Laplacian (Vollmer et al., 2001), and Taubin's low pass filter (Taubin, 1995). In Figure 3, the reconstruction procedure of a patient-specific AAA model is given.

Figure 3

Reconstruction procedure for a patient-specific AAA model [The figure is adapted from Les et al. (2010) and used with permission].

Governing Equations in Solid Domain

The governing equation in solid domain is the momentum conservation given in Equation (2) (Scotti et al., 2008).

In Equation (2), τ is solid stress tensor; is body forces term per unit volume; ρ is AAA wall mass density; and is local acceleration of solid. Lagrangian description is typically used to track the deformation of the solid domain (Donea et al., 1982; Bathe and Zhang, 2004). Neglecting the gravitational forces () is a commonly used approach due to their insignificant effect on wall stresses (Wolters et al., 2005).

Important Stress Parameters for Solid Domain

For the solid domain, principal wall stresses (σ1, σ2, σ3), Von Mises stress and wall displacements are critical parameters for AAA rupture risk assessment. The peak wall stresses are better indicators when compared to the maximum AAA diameter, since rupture is the mechanical failure of the wall where the strength of wall is not sufficient to withstand the peak wall stress (Fillinger et al., 2003; Wolters et al., 2005). Von Mises stress is a measure used for failure prediction based on three principal stresses as given in Equation (3).

In Equation (3), the term in the left side is square of Von Mises stress, σ is the local principal stress and σ is the uniaxial failure strength of the wall (Scotti et al., 2008).

Incorporation of Intraluminal Thrombus (ILT) to the Models

In solid domain, the presence of intraluminal thrombus (ILT) is an influencing factor for the wall stress depending on its shape, size, and material properties (Mower et al., 1997; Di Martino et al., 1998; Di Martino and Vorp, 2003). Flow hemodynamics is also closely related with the ILT formation on AAA wall. ILT contains immune and inflammatory response agents affecting the evolution of the disease (Adolph et al., 1997). In a related clinical study, ILT is observed near the site of rupture for 80% of the autopsies (Simão da Silva et al., 2000), as being a source of proteolytic activity, local wall thinning, wall weakening, and hypoxia (Swedenborg and Eriksson, 2006; Houard et al., 2007). ILT morphology can exhibit a layered structural material behavior or fluid-like homogeneous medium. It becomes a poroelastic material at a certain state of maturity which is affecting blood transport (Tong and Holzapfel, 2015). Flow stagnation and associated low WSS, high residence time of platelets and monocytes in AAA sac contribute to high potential of wall-cell adhesion forming ILT (Kelsey et al., 2016). In a recent study, Di Achille et al. (2016) predicted ILT formation and progression sites in patient-specific AAA models, using a phenomenological metric of thrombus deposition potential which is indicating a balance between the thrombogenicity and hemodynamic shear forces on the wall. The relative overlap between predicted and actual thrombus covered areas in six patient-specific AAA models was reported around 80 ± 16% (Di Achille et al., 2016). Therefore, ILT formations on vessel walls should be assigned with appropriate material constants with accurate geometric representations for model accuracy since ILT significantly affects both the wall mechanics and flow hemodynamics, indicating the importance of interactions between solid and fluid domains.

Wall Material Properties

Aorta consists of three layers: the intima, media, and adventitia (Lasheras, 2006). Intima and adventitia are the inner and outer layers, respectively. The distribution of layer thicknesses has a ratio about 20:47:33 for intima: media: adventitia (Humphrey and Holzapfel, 2012). Elastic modulus ratio for intima: media: adventitia is ~1:3:2 (Khanafer and Berguer, 2009; Gao et al., 2013). The highest stresses are observed in the media layer due to the normal pressure of cyclic flow (Simsek and Kwon, 2015). The maximum wall deformation can reach up to 2.2 mm at the systolic phase (Canchi et al., 2018). One of the possible reasons of AAA formation is reported as the medial loss in the arterial wall with degeneration of smooth muscle cells (Humphrey and Taylor, 2008). In case of AAA wall weakening depending on the media loss, elastic modulus of the media layer can be reduced by 20 times (Feng et al., 2008). Modeling AAA wall as a single-layered structure is a commonly used approach for simplification of the problem.

Governing Equations in Fluid Domain

Flow velocity, pressure and wall shear stress (WSS) exerted by fluid viscous forces are determined by solving the Navier-Stokes and continuity equations. For an incompressible and homogeneous fluid, the Navier-Stokes equations can be defined using the Arbitrary Lagrangian Eulerian (ALE) description given in Equation (4) and continuity equation given in Equation (5) (Donea et al., 1982; Zhang et al., 2003; Scotti and Finol, 2007).

The fluid velocity vector is denoted by v; time is denoted by t; the velocity of the fluid domain (i.e., moving coordinate velocity primarily due to FSI) is denoted by w; fluid stress tensor is denoted by τ; and body forces term is denoted by . The fluid stress tensor (τ) is defined in Equation (6), in terms of fluid pressure (p), Kronecker delta (δ), dynamic viscosity (μ), and strain rate (ε). The strain rate can be written in terms of velocity vector (v), as given in Equation (7). The effect of gravitational acceleration is not critical and the body forces on the fluid () can be neglected (Scotti et al., 2008).

Important Hemodynamic Parameters for Fluid Domain

In fluid domain, wall shear stress (WSS), oscillatory shear index (OSI), intraluminal pressure, flow path, and flow velocity are of main interest. Due to complex AAA geometries, WSS has spatially and temporally complicated distributions (Arzani and Shadden, 2015). WSS is a measure of flow-driven tangential forces per unit area on the AAA wall. The magnitude of WSS at the fluid-structure interface can be determined by multiplying the viscosity (μ) with the local shear rate () (Wolters et al., 2005). In large arteries, WSS amplitudes typically range from 1 to 5 Pa, therefore WSS values smaller than 1 Pa are evaluated as low WSS (Ene-Iordache and Remuzzi, 2012; Qiu et al., 2018). Time averaged WSS (TAWSS) can be determined as given in Equation (8) where T is the integration period (Arzani et al., 2014; Arzani, 2018).

OSI is a measure defining the unidirectionality of shear stress which is sensitive to turbulence and it can be determined as given in Equation (9) (Arzani et al., 2017).

When OSI is zero, it is indicating that shear stress is unidirectional. If OSI is 0.5, it means that time average of shear stress is zero. Endothelial cell activation potential (ECAP) is the ratio of OSI and TAWSS as given in Equation (10). ECAP is used to characterize ILT susceptibility on AAA wall (Di Achille et al., 2014). Critical threshold value of ECAP is stated as 1.4 Pa−1, where above this value, there is a high potential of ILT formation (Kelsey et al., 2016). Therefore, as ECAP increases due to emergence of circulatory flows in diseased states, there is a higher chance of wall inflammation.

Fluid Properties and Flow Regime

Blood has non-Newtonian characteristics where the viscosity decreases with increased shear rate. AAA flow is inherently pulsatile and has a turbulent nature particularly at mid-diastolic phase. At the systolic phase, peak blood flow velocity in AAA can reach up to 20 cm/s (Casciaro et al., 2018). Larger aneurysmal diameters lead to higher turbulence intensity and large recirculating vortices increasing the blood residence time in the aneurysm sac, particularly at the diastolic phase of the physiological pulsatile flow (Khanafer et al., 2007).

In many studies, AAA flow is considered as unsteady laminar flow (Scotti et al., 2008; Chandra et al., 2013; Morris et al., 2013; Soudah et al., 2013; Owen et al., 2016) due to not exceeding the threshold Reynolds number (2000–2300) for transition to turbulence in pipe flow. However, complex flow geometry, sudden lumen expansion in AAA sac and pulsatile nature of flow are the factors that might trigger the turbulence in low Reynolds numbers (Poelma et al., 2015). In order to resolve the turbulent effects accurately, Les et al. (2010) and Arzani et al. (2014) applied more demanding direct numerical simulation (DNS) approach, and Khanafer et al. (2007) used k-ω turbulence model.

Biochemical Transport

Oxygen (Sun et al., 2009), low density lipoproteins (Choudhury et al., 2019), and chemical species transported by blood flow have an influence on the AAA progression and ILT formation on the wall. Shear rate is an Eulerian measure that cannot quantify particle transport (Shadden and Arzani, 2015). Locally high surface concentration of chemical species are not always corresponds to the regions with low WSS that are determined using Eulerian mass transport (Choudhury et al., 2019). The residence time and flow path of the micro-particles can be determined using Lagrangian mass transport (Arzani et al., 2014, 2017). By this way, the particles in the blood flow can be tracked and their interaction with AAA wall can be modeled using convection, diffusion, and reaction equations coupled with the computed flow field (Biasetti et al., 2012). These fluid-chemical models provide better insight to understand AAA and ILT pathophysiology.

Setting Up and Solution of the Solid Domain via FEA

Defining Material Constants

AAA wall has hyperelastic (i.e., non-linear stress-strain behavior), viscoelastic (i.e., time varying response due to relaxation) and anisotropic (i.e., direction dependent) material properties due to its layered and fiber-oriented structure (Vande Geest et al., 2006a,b). In order to determine the wall deformation and stress in a reliable manner, constitutive equations are needed to model wall stress-strain behavior. As a first approximation, AAA wall can be modeled as a linearly elastic, homogeneous medium using elastic modulus of 2.7 MPa, Poisson's ratio of 0.45, and mass density of 2,000 kg/m3 (Di Martino et al., 2001). For improved accuracy, non-linear stress-strain behavior should be taken into account using hyperelastic models (Scotti et al., 2008). Raghavan and Vorp (2000) obtained an experimental fit for non-linear stress-strain curve of AAA wall. Experimental data can be represented using Mooney-Rivlin hyperelastic material model considering a strain energy density function (W) given in Equation (11) (Rivlin and Saunders, 1951).

In Equation (11), I is the ith invariant of the left Cauchy-Green tensor, c is the material parameter for fitting the experimental data. Chandra et al. (2013) used second order Mooney-Rivlin approach for modeling the AAA wall and ILT as given in Equations (12,13), respectively.

Considering the population averages (Raghavan and Vorp, 2000; Vande Geest et al., 2006a,b), material parameters in Equations (12) and (13) are determined as c10 = 17.4 N/cm2, c20 = 188.1 N/cm2, c01 = 7.98 N/cm2, and c02 = 8.71 N/cm2. Mass densities of AAA wall and ILT are assumed to be 1.2 and 1.1 g/cm3, respectively.

Boundary Conditions in Solid Domain

In general, AAA tissue is isolated by applying zero translation and zero rotation fixed on the inlet and outlet of the solid domain (Scotti et al., 2005). The surrounding branch arteries of aorta produce tethering effect on the AAA model. This tethering effect is reflected to the wall either by using previously described zero displacement and rotation fixed on the inlet and outlet boundaries; or alternatively, an axial stretch can be defined on these solid inlet and outlet boundaries as suggested by Tang et al. (2005), since arterial wall is physiologically under tension in reality (Holzapfel et al., 2000). Scotti et al. (2008) applied 5% axial stretch on the inlet and outlet boundaries of the solid for modeling the tethering effect.

In addition to inlet and outlet boundaries, intra-abdominal pressure and contact of AAA with the surrounding tissues and organs can also be considered for defining boundary conditions (BC) in solid domain. Scotti et al. (2008) applied intra-abdominal pressure of 12 mmHg (Hinnen et al., 2005) on the outer AAA wall. In reality, AAA is located along the spinal column and there is a certain physical contact between AAA and spinal column which provides additional support on the posterior wall. This can be modeled with an additional BC on the posterior wall by fixing the contact surface with zero displacement (Scotti et al., 2008).

If solid FEA model is not coupled with a flow domain, then intraluminal pressure BC is applied on the inner AAA wall as a steady or transient load. This analysis is also known as computational solid stress (CSS) approach. In case of a FSI simulation, in which the solid FEA model is coupled with a flow domain, inner AAA wall is set as a FSI boundary surface.

Setting Up and Solution of the Fluid Domain Via Computational Fluid Dynamics

Modeling Non-Newtonian Fluid Viscosity

Computational fluid dynamics (CFD) simulations are performed to numerically determine the flow variables such as flow velocity, pressure and WSS by modeling the fluid viscosity. Newtonian fluids have constant viscosity, but for non-Newtonian fluids such as blood, the viscosity changes as function of shear rate. For large arteries, the shear strain rate in the flow exceeds 50 s−1 and the viscosity demonstrates nearly constant behavior depending on the high shear rate (Young, 1979). Therefore, modeling the blood as a Newtonian fluid is a common approach for the AAA flow. Recently, Arzani (2018) reported that blood residence time affects the viscosity and traditional non-Newtonian models may exaggerate the shear-dependent viscosity behavior of blood, thus it is recommended to use Newtonian models if minor portion of AAA is accompanied by high blood residence time.

Blood mass density is usually taken as 1.05 g/cm3 with a constant dynamic viscosity (μ) of 0.035 Poise (Chandra et al., 2013). Khanafer et al. (2006) compared Newtonian and non-Newtonian fluid models for the same AAA geometry and stated that the maximum pressure and maximum WSS differences were 2.53 and 26.7%, respectively. The higher difference in maximum WSS is due to the effect of near-wall turbulence which is underestimated by Newtonian models. Non-Newtonian behavior of blood can be modeled using Carreau-Yasuda viscosity model parameters given in Equation (14).

In Equation (14), is shear rate, η is shear rate dependent viscosity, η∞ is viscosity at high shear rate, η0 is viscosity at low shear rate, and λ, a, n are empirical constants. By approximating the blood viscosity measurements of Thurston (1979), empirical parameters yield as η∞ = 0.00476 Pa s, η0 = 0.0519 Pa s, a = 0.409, n = 0.191, and λ = 0.438s.

Fluid Domain Boundary Conditions

Inlet and outlet BC are defined considering the physiological pulsatile flow. Commonly, time-varying flow velocity given in Figure 4A is defined at the inlet of fluid domain and time-varying intraluminal pressure given in Figure 4B is prescribed at the outlet (Scotti et al., 2008). No slip BC is applied on the wall considering the viscous blood flow.

Figure 4

Waveforms of applied boundary conditions at the inlet and the outlet of the fluid domain. Dash lines show the moment of peak values in the waveforms. (A) Sample inlet flow velocity profile. (B) Sample outlet pressure profile [The figure is adapted from Scotti et al. (2008) and used with permission].

The time-dependent inlet flow rate given in Figure 4A can be applied considering different velocity profiles. Applying plug velocity profile is the simplest approach where the axial flow velocity is uniform through the inlet boundary surface (Bluestein et al., 2009). Alternatively, fully developed velocity profile can be used considering Womersley profile which has zero velocity on the wall and has a parabolic distribution with maximum flow velocity at the midpoint of the inlet surface, which is more realistic due to including the effect of viscous boundary layer (Womersley, 1955; Papaharilaou et al., 2007). To enhance the accuracy of the model, patient-specific inlet velocity profiles can also be applied by obtaining realistic profiles non-invasively using phase-contrast magnetic resonance images (PC-MRI) (Kose et al., 2006).

After determining patient-specific inlet BC, distribution of mass flow rate to the branches of mesenteric, renal, and iliac arteries should also be considered, since they are significantly affecting the mass flow rate at the outlet of AAA. One-dimensional arterial tree model can be used to determine appropriate inflow and outflow conditions to better reflect the reality for biological relevance (Formaggia et al., 2003; Wolters et al., 2005). Les et al. (2010) modeled the effect of branch arteries at the downstream vasculature using three-element Windkessel model considering the capacitance, proximal resistance and distal resistance at the downstream. Using PC-MRI scanning, mean flow rates at supraceliac and infrarenal locations are measured as 3.51 and 1.31 L/min, respectively. The difference of 2.2 L/min is distributed to the arterial branches between supraceliac (1 cm above celiac artery) and infrarenal (1 cm below the most distal renal artery) levels (Les et al., 2010). It is reported that using zero-pressure outlet BC is not sufficient, and Windkessel BC provides more realistic flow and pressure features.

For determining patient-specific outlet pressure BC, a catheter can be placed inside AAA sac, however this invasive procedure is not preferred by clinicians during AAA repair (Chandra et al., 2013). Alternatively, a non-invasive method can be applied to obtain patient-specific outlet pressure. van ‘t Veer et al. (2008) compared non-invasive brachial cuff blood pressure measurements with invasive catheter pressure measurements inside AAA sac. Brachial cuff pressure measurements resulted in 5% underestimation for systolic phase, and 12% overestimation for diastolic phase compared to the intraluminal pressure in AAA. Using these correlations, a patient-specific estimation of fluid outlet pressure can be predicted non-invasively.

Solution and Mesh Independence for Solid and Fluid Domains

Solid and fluid domains are spatially discretized using a generated mesh composed of tetrahedral, hexahedral or polyhedral elements as given in Figure 5. In fluid domain, mesh density should be relatively higher in the boundary layers close to the wall and in regions that are expected to have high velocity gradients as well. Obtaining a mesh-independent solution is a requirement prior to the in depth analysis of AAA hemodynamics. Typically, if a relative difference of 2% is achieved for variables obtained by different mesh sizes, the results can be considered as mesh-independent (Kelsey et al., 2016). The quality of CFD mesh resolution can also be checked using the minimum Kolmogorov length scale. For large eddy simulation approach, the ratio of mesh element size to the Kolmogorov length scale should be at most 20 and 40 for fine and coarse fluid meshes, respectively (Celik et al., 2009).

Figure 5

Discretization of the problem domain using finite element meshes. (A) Structured hexahedral fluid mesh using idealized AAA model. (B) Structured hexahedral solid mesh using idealized AAA model. (C) Unstructured tetrahedral fluid mesh using patient-specific AAA model. [A,B are adapted from Scotti et al. (2008) and used with permission, (C) is adapted from Wolters et al. (2005) and used with permission].

In order to check the mesh-independency of results, solution of first three cardiac cycles should be ignored. Because, periodically converged results are generally obtained after the third cycle. A time step of 0.001 s (~1/1,000th of cardiac cycle) is suitable for obtaining unsteady solutions (Les et al., 2010). If stability problems arise in the analysis, the time step can be decreased to achieve convergence. Courant number (Co ≡ u(t)/x) should be equal or lower than one for the increased stability. The Courant number shows how much information traverses (i.e., Flow velocity, u) the length of a mesh element (x) within a time step (t).

Fluid-Structure Interaction (FSI)

To simulate deformation of AAA tissue under the effect of blood hemodynamics accurately, FSI approach needs to be adapted. This is because of the strong interactions between flowing blood and vessel walls. Blood flow generates unsteady forces on vessel walls that causes deformation of the walls. These deformations in turn influence blood flow patterns. Therefore, vessel wall behavior cannot be predicted accurately if these counter interacting forces are ignored. The commercial software packages such as ANSYS, ADINA, ABAQUS, COMSOL are commonly used for FSI modeling.

FSI modeling can be performed using three different numerical approaches, which are 1-way uncoupled (also known as FEA or CSS), explicit 2-way coupled and implicit 2-way coupled methods. 1-way uncoupled method is a static investigation where spatially uniform or non-uniform intraluminal pressure load is applied to the inner AAA wall without considering the effects of dynamic flow and only provides solid domain parameters such as wall stresses and wall displacements. In a related study (Scotti et al., 2008), it is stated that 1-way uncoupled approach underestimates the maximum wall stress due to neglecting the hemodynamic effects. To have a better understanding, FSI analysis should be performed using 2-way coupled methods. In a numerical investigation, 1-way uncoupled and 2-way coupled FSI simulations are compared using the same patient-specific AAA model and it is reported that 1-way uncoupled FSI resulted in 14, 4, and 18% difference in peak principle stress (σ1), principle strain (ε1), and WSS, respectively (Chandra et al., 2013).

In 2-way coupled methods, a FSI boundary is defined between the blood and AAA wall. On this FSI boundary surface, solutions of solid and fluid domains are coupled considering displacement compatibility and traction equilibrium given in Equations (15,16) (Bathe et al., 1999). On FSI boundary, solid and fluid displacement vectors are denoted by d and d; solid and fluid unit normal vectors are denoted by and , respectively.

Two-way coupled FSI can be performed using explicit or implicit approaches depending on the level of physical interaction between the fluid and solid. When a strong interaction exists as present in AAA, explicit approach will be insufficient, and implicit approach should be preferred for enhanced accuracy (Amindari et al., 2017). For implicit approach, there are two options as iterative implicit or fully coupled. In fully coupled implicit method, all governing fluid and solid equations are solved simultaneously, requiring high computational memory and leading to an excessive solution time. However, iterative implicit approach requires less computational power, and consists of a number of coupling iterations until reaching interaction convergence at each time step. Using smaller time steps lead to more stable results. Increasing the number of coupling iterations and using relaxation factors will help to overcome the stability problems.

Recent Findings on Growth and Rupture Mechanics of AAAs From Computational Studies

It is now widely accepted that biomechanical factors play role in degeneration and eventual rupture of AAA tissue. There is controversy regarding which biomechanical parameters are important in rupture. Below we summarize the important findings by categorizing the effects of various parameters.

Effect of Wall Thickness and AAA Asymmetry

Scotti et al. (2005) investigated the effect of varying wall thicknesses by performing fully coupled FSI simulations using idealized AAA geometries based on medical images. The thickness of AAA wall has an average of 1.45 mm, but it can decrease to 0.23 mm near the rupture site (Raghavan et al., 2004, 2006). A variable thickness between 0.5 and 1.5 mm was distributed along the AAA wall as inversely proportional to the cross-sectional diameter. When variable wall thickness model was compared with a uniform wall thickness (1.5 mm) model, it was observed that peak wall stress was increased by four-times for the variable wall thickness case.

The same research group also investigated the effect of AAA asymmetry (Scotti et al., 2008). Different levels of asymmetry were modeled using idealized geometries with the same patient-specific inlet velocity BC. For the most asymmetric model with β = 0.2 (see Equation (1) for definition of β), AAA diameter expanded by 15.2% at the peak systolic pressure. When an axisymmetric (β = 1) AAA model is considered, the diameter expansion is observed as 12.8%, implying that increasing asymmetry resulted in higher AAA deformation and therefore higher peak wall stress on the wall. As AAA becomes more asymmetric, location of peak wall stress shifted from anterior to the posterior wall.

Effect of Proximal Neck and Iliac Bifurcation Angles

Drewe et al. (2017) performed FSI simulations to investigate the effects of proximal neck and lateral iliac bifurcation angles considering idealized AAA models. Recent morphological comparisons showed that AAAs with large iliac bifurcation angle have a lower rupture risk (Drewe et al., 2017). Proximal neck angle has less impact on AAA hemodynamics compared to the iliac bifurcation angle. When the iliac bifurcation angle increased from 30 to 150°, peak WSS increased more than 2-fold (from 2.91 to 6.19 Pa), peak von Mises wall stress increased about 30% (from 0.186 to 0.243 MPa) and ECAP decreased by 57% (from 11.41 to 7.25). Larger iliac bifurcation angle was more protective of ILT formation and AAA expansion due to provoking high WSS and low ECAP conditions (Drewe et al., 2017). However, excessive load on the iliac arteries with increased bifurcation angle may increase the risk of an iliac artery aneurysm initiation (Xenos et al., 2010).

Effect of AAA Diameter and Wall Stress

The maximum AAA diameter is the first indicator for the treatment. For the current practice, when maximum AAA diameter exceeds 5–6 cm or diameter growth rate is higher than 1 cm per year, open surgery or endovascular treatment methods are performed considering the life expectancy of the patient (Scott et al., 2002; Longo and Upchurch, 2005; Chandra et al., 2013). Canchi et al. (2018) performed a comparative FSI study, considering two patient-specific AAAs with maximum aneurysmal diameters of 3.5 and 7 cm. Maximum principle stresses were determined as 0.30 and 0.22 MPa for 3.5 and 7 cm AAA diameters, respectively, implying that size of AAA is not the sole determinant for the rupture risk. It is reported that maximum wall stress was 12% more accurate for predicting rupture compared to using maximum AAA diameter alone as an indicator (Fillinger et al., 2002, 2003).

The maximum wall stress is generally observed at the transition of sac to neck of AAA wall as shown in Figure 6 (Chandra et al., 2013; Doyle et al., 2014). FEA investigations showed that peak wall stress on the posterior AAA wall was within the range of 290 to 450 kPa (Raghavan et al., 2000), while non-aneurysmal aorta had a peak stress around 120 kPa. Scotti et al. (2008) compared FSI and FEA approaches on the same AAA model, and showed that FSI resulted in 25% increased wall stresses (peak values between 275 and 398 kPa) compared to FEA. Chandra et al. (2013) also performed FSI and obtained relatively higher peak wall stresses between 750 and 870 kPa. The reason for the high peak wall stress was considered to be the non-uniform AAA wall thickness which could increase the stresses at low thickness areas.

Figure 6

Von Mises wall stress contour plot from different views, on a patient-specific AAA geometry at peak systolic pressure. The middle figure is a sectional view. High stress region exists in the posterior wall, however the peak wall stress is observed at the proximal neck of the aneurysm [The figure is adapted from Doyle et al. (2014) and used with permission].

Effect of WSS and OSI

WSS and OSI are important hemodynamic parameters in turbulent AAA flow. Les et al. (2010) investigated patient-specific AAA hemodynamics using high-resolution CFD simulations (about 8-million mesh elements), hypothesizing that physical lower limb exercise might decrease the growth rate of AAA, since exercise resulted in a high amplitude WSS pattern and lowered OSI on AAA wall. Moderate turbulence was observed in AAA during exercise, while resting conditions led to mild turbulence. Rest-to-exercise TAWSS changes were found to be statistically significant. For example, at supraceliac level, TAWSS at rest was 3.6 dyn/cm2, and it increased to 9.2 dyn/cm2 during the exercise. At mid-aneurysm level, TAWSS at rest was 7.3 dyn/cm2, while at exercise it was 21.7 dyn/cm2. OSI values at rest were 0.28 and 0.27 for supraceliac and mid-aneurysm locations, respectively. These OSI values decreased to 0.18 and 0.21 during the exercise.

Qiu et al. (2018) performed CFD simulations using three ruptured and one non-ruptured patient-specific AAA models. The rupture sites were found to be near the fluid stagnation regions which have nearly zero WSS with high WSS gradients (WSSG). In Figure 7, the distribution of WSS and OSI are provided on a patient-specific model. Most researchers agree that the locations with low WSS, high OSI and high ECAP are prone to thrombus formation and have a higher risk of rupture (Les et al., 2010; Kelsey et al., 2016). On the other hand, in some studies reporting controversial results (O'Rourke et al., 2012; Arzani et al., 2014; Mohamied et al., 2015; Singh et al., 2018), it is stated that low WSS and high OSI regions do not coincide with thrombus deposition and atherosclerosis sites. Therefore, the exact effect of these hemodynamic parameters on the rupture mechanism is not yet fully understood.

Figure 7

TAWSS (Time averaged WSS), OSI, and ECAP contour plots on a patient specific AAA model. ILT formation is more likely to be observed in high ECAP regions. Rupture risk and aneurysm growth rate increase at regions with high OSI and low TAWSS [The figure is adapted from Kelsey et al. (2016) and used with permission].

Effect of Vascular Growth and Remodeling

Arterial growth and remodeling are investigated and modeled in several studies (Humphrey and Rajagopal, 2003; Watton et al., 2004; Gleason and Humphrey, 2005; Valentín et al., 2011, 2013; Humphrey and Holzapfel, 2012; Karšaj and Humphrey, 2012). In a recent study, Wu and Shadden (2015) presented a computational framework by coupling blood flow hemodynamics with vascular growth and remodeling (G&R). The wall is treated as a constrained mixture consisting of anisotropic collagen, elastin and smooth muscle fibers. Depending on the mechanical stimuli on the wall, the vessel has the ability of adaption for returning back to its homeostatic state by means of removal of old constituents and production of new ones. The lifespan of collagen is within 70–80 days (Wilson et al., 2013), where the flow simulation cycle has a time scale around 1 s. For computational efficiency, flow simulations only performed in cases when G&R caused severe change in geometry and boundary conditions, meaning that the last flow simulation results were used for all time in between until a new one is performed. The deviations of wall tension and WSS were used as the influencing factors for mass production rate on the vessel to return its homeostatic state. WSS changed about one order of magnitude in the aneurysm site and contributed to G&R. Stresses on the wall converged to homeostatic values after 110 G&R time steps (equivalent to about 800 days) and significant wall expansion was observed at the entire mass loss introduced regions.

Effect of Surgical and Endovascular Treatment Methods

Sughimoto et al. (2014) evaluated AAA surgery by comparing pre-operative and post-operative hemodynamics on CT-based geometries. AAA was replaced with a 30 mm straight graft under cardiopulmonary bypass. They proposed a novel parameter of pulsatile energy loss index (PELI) which evaluates the energy loss between inlet and outlet of interested artery. After AAA treatment, PELI decreased from 0.986 (pre-operative value for the whole aorta) to 0.820 (see Sughimoto et al. (2014) for calculation of PELI), indicating that the grafting procedure improved the energy efficiency of blood flow delivery and reduced the left ventricle afterload. For a young healthy adult, PELI was measured as 0.0215, which was significantly low when compared to the AAA patient.

Casciaro et al. (2018) performed CFD simulations before and after two different endovascular surgical treatments of AAA. The Nellix endograft (Endologix, Irvine, California) proposed a method of AAA treatment based on endovascular aneurysm sealing (EVAS). For conventional endovascular aneurysm repair (EVAR), endografts have proximal fixation mechanism. The lateral neck angle (Kandail et al., 2015) and implantation position (Raptis et al., 2017) of the endografts alter the hemodynamics of EVAR. On the other hand, EVAS consists of two stents which are separately engaged to left and right common iliac arteries. Expandable endobags surround the stents and provide continuous sealing along the internal surface of the aneurysm (Böckler et al., 2015). Due to significant anatomic modifications after EVAS, 2-fold pressure increase was seen at the level of renal arteries. In addition, the peak flow velocity inside EVAS endograft stents was 60 cm/s, which was three-times higher than the peak blood flow velocity after EVAR treatment. This 3-fold increase in peak flow velocity after EVAS, also resulted in 60% higher WSS compared to EVAR (Casciaro et al., 2018).

Experimental Investigation of Hemodynamics of AAAs

AAA hemodynamics is quite complex primarily due to irregular shape, flexible arteries, turbulent flow, and non-Newtonian behavior of blood. As computational methods have widely been used in investigation of hemodynamics and mechanical behavior of arterial tissue, experimental techniques are also utilized in characterization of flow dynamics through AAAs. Both approaches are crucial and complement each other with offering in depth analysis where the corresponding results depend on the assumptions in computational models and simplifications in experimental methods.

Various qualitative and quantitative flow measurement techniques have been utilized for experimental investigation of AAAs. In early studies, qualitative flow visualization techniques have been widely implemented, in which using localized injections, the patterns of dye as streaklines were generated (Ku et al., 1989). This simple method provides insights on overall behavior of flow structure in the region of interest. In terms of detailed velocity information in the AAAs, non-intrusive, and quantitative techniques such as MRI and Doppler Ultrasound Imaging (DUS) are appropriate for detailed anatomical analysis by means of tomographic slices. In MRI based system encoding the flow velocity is obtained by means of the changes in MR signal phase along a magnetic field gradient (Wang et al., 2016). Moreover, laser-based techniques including Laser Doppler Anemometry (LDA) and Particle Image Velocimetry (PIV) are also implemented for detailed velocity information. In LDA, pointwise velocity measurements are performed by means of Doppler effect using two laser beams for each velocity component. This method has very high spatial and temporal resolutions (~kHz) and allows to measure reverse flow and turbulent fluctuations, which is quite critical for understanding the growth and the rupture mechanics of AAAs (Yip and Yu, 2002). In PIV, two component velocity information on 2-D plane is obtained at relatively lower temporal resolution (~Hz) using two laser sheets and an advanced camera, where in stereoscopic PIV three-velocity components are obtained on 2-D plane using two cameras. This method is quite effective in understanding the key parameters of disturbed hemodynamics of AAAs since it provides detailed velocity field information and therefore have widely been used recently (Deplano et al., 2013; Chen et al., 2014). For the quantification of three-dimensional velocity field, PIV offers two recent state of the art versions: Tomographic PIV and Holographic PIV, which are expected to appear in the studies of disturbed hemodynamics through AAAs in near future. Tomographic PIV has recently been applied to the hemodynamics of intracranial aneurysms (Roloff et al., 2017, 2018).

Experimental studies on blood flow inside AAAs is primarily categorized and populated in two major groups: the properties of the aorta phantom and conditions of the flow. Major consideration for the aorta phantom are elasticity of arterial wall and geometry of the aorta (change in wall diameter, aneurysm shape, artery bifurcation, etc.) whereas major considerations in flow conditions, include flow waveform (steady or physiological), type of working fluid (Newtonian, non-Newtonian), and pulse rate (exercise or resting). Considering the overall flow structure in aorta and AAAs, fully developed flow at the entrance of the AAA expands and creates counter rotating recirculation along with a jet flow at the center, which are due to the adverse pressure gradient imposed at the aneurysm bulge. Location and strength of the recirculation vortex primarily depend on the inlet flow waveform, bulge shape, and elasticity of arterial wall. For the steady inlet flow condition, which is not realistic, the recirculating vortex is larger in extent in average and located closer to downstream of the bulge, whereas in physiological flow condition recirculation region is alternating and moving back and forth due to periodic flow condition. It is also indicated that the spatial extent, the location, and the strength of the recirculation vortex are significantly affected by the bulge shape and the wall rigidity (Egelhoff et al., 1999; Yu, 2000; Deplano et al., 2007; Meyer et al., 2011). Flow structure inside the AAAs varies significantly according to these considerations, as reported by many experimental studies (Egelhoff et al., 1999; Yu, 2000; Deplano et al., 2007; Meyer et al., 2011).

Therefore, because of the complexities in geometry and flow conditions, the experimental set up for AAA investigations need to be designed to mimic natural realistic in-vivo conditions within AAA as close as possible to obtain accurate results. A typical experimental set up for the analysis of hemodynamics contains flow circulatory system including pump, piping, pressure compliance, test section that contains artery model, blood mimicking fluid, and flow measurement systems as briefly mentioned previously such as MRI, LDA, or PIV. In this part of the paper, the components of experimental set up and the techniques that are used widely in literature are explained in detail, and the results of recent experimental studies are discussed.

Flow Circulatory System

In cardiovascular biomechanics, it is difficult to measure hemodynamic quantities clinically due to the constraints of the system. Conducting in vivo animal studies are also very challenging because of the geometrical complexities and ethical issues. Therefore, in recent years, simulating in vivo conditions by using computational models or experimental setups become very popular. In vitro studies can be performed by generating the physiological flows through anatomical geometries of interest to simulate exact in vivo conditions in laboratory test setup; i.e., flow rates and pressure at the inlet and the exit, geometry of the artery and working fluid of a specific artery should be replicated physiologically. In Figure 8, an example of a test setup is presented. The main component of such a setup is the flow source, which is generally a pump, supplying a fluid flow into downstream. At the downstream, there are resistances to simulate physiological pressures. This is generally accomplished by using compliance chambers; in other words, vascular simulators to mimic realistic arterial pressures for specific arteries. Besides mimicking the physiological flow conditions, arterial phantoms which are desired to be studied should also be modeled and manufactured (Doyle et al., 2008, 2009b; Corbett et al., 2009). Finally, a blood mimicking fluid should be selected that matches biomechanical behaviors of blood.

Figure 8

Circulatory flow loop model for flow measurement system of abdominal aortic aneurysm containing tank, physiological flow generator, flowmeters, pressure sensors, compliances, and AAA box [The figure is adapted from Deplano et al. (2013) and used with permission].

Pump

Blood flows through arteries in an unsteady manner, which is characterized by a rhythmic repeatability and called pulsatile, i.e., such flows repeat themselves with a certain frequency. The pulsatile physiological flow profile, is presented in Figure 4, and has a complex waveform. In earlier studies, although it is far from the exact physiological conditions, researchers use constant flow rate in hemodynamic analyses. This is primarily due to the simplicities both in generating the steady flow in a circulatory loop and in performing flow measurements in such systems. In that sense, many experimental studies have been conducted at steady flow conditions using different type of flow generators, such as variable speed centrifugal or gear pumps (Ku et al., 1989; Asbury et al., 1995; Bluestein et al., 1996; Boutsianis et al., 2008), steady submersible pumps (Chen et al., 2014, 2015) or sometimes, by only using the effect of gravity, which requires placing a head tank (O'Rourke and McCullough, 2010).

However, a better description of hemodynamics of AAA closer to the physiological conditions is needed for the complete understanding of growth and rupture behavior of AAA. Indeed, comparison of steady and pulsatile flow conditions reveals important differences in hemodynamics which are also influenced by other system parameters including wall rigidity and AAA bulge shape (Yu, 2000). For example, with steady flow conditions, a vortex ring is localized at the distal end of the bulge, creating oscillations in wall shear stresses proximal to the corresponding localized region. On the other hand, with pulsatile flow case, the vortex ring is not localized and appears at the proximal site at the early systolic phase proceeding through the downstream throughout the cardiac cycle (Yu, 2000). This means that extrapolating steady flow results in extreme conditions might underestimate certain key parameters in physiological flow conditions. Considering the steady flow generators, centrifugal pumps, gravity-driven systems, and gear pumps have been commonly used. Later on, researchers have adapted various components to generate pulsatility at the downstream of these steady flow generators, including computer controlled valves (Egelhoff et al., 1999), piston-cylinder arrangements (Stamatopoulos et al., 2010, 2011), and gear pumps (Mechoor et al., 2016). In addition, commercial pumps which can generate either steady or pulsatile flow waveforms for several arteries are also available. All these alternatives for pumps are discussed along with addressing the pros and cons of each of them.

Steady flow generators

Centrifugal pumps. The most widely used flow generators are steady flow pumps, where the centrifugal pumps (Asbury et al., 1995; Bluestein et al., 1996) and submersible pumps (Chen et al., 2014, 2015) are generally used for that purpose. Centrifugal pumps are cost effective and are easily available in market with a very wide product range, however the flow rate is highly affected by the pressure drop, and thus, is not very convenient for hemodynamic analyses.

Gear pumps. Gear pumps are positive displacement pumps and are convenient for providing quite stable and constant flow rates at broad range of pressure drop conditions. In earlier studies, variable speed gear pumps have been employed in hemodynamics studies to generate steady flow conditions (Asbury et al., 1995; Boutsianis et al., 2008).

Gravity-driven systems. In gravity driven systems, the effect of gravity is used as the driving force with locating a head tank of fluid at certain elevation above the test section to maintain pressure and velocity distribution in AAA. With the effect of gravity, steady mean flow can be generated, but in order to maintain the physiological flow pattern, an additional equipment such as computer-controlled valves (Egelhoff et al., 1999) or piston-cylinder arrangements (Yu, 2000; Stamatopoulos et al., 2010, 2011) are needed to be used at the downstream.

Pulsatility generators

In order to generate pulsatility, additional computer-controlled equipment should be adapted to circulatory systems. For this purpose, the most common integrated systems include computer-controlled valves, piston-cylinder arrangements, or gear pumps.

Valves. In general, centrifugal pumps are adapted only for steady circulations, where in a very few studies, a flapper nozzle valve is adapted to the flow loop to produce pulsatility (Moore et al., 1992; Moore and Ku, 1994). In general, the desired waveform including physiological pattern with the reverse flow cannot be maintained properly in centrifugal pump driven systems. An application of valves in physiological flow generation is utilized by Egelhoff et al. (1999), where a gravity driven pump with a head tank kept 234 cm above the AAA test section to provide desired entrance flow conditions and a computer controlled diverter valve at the downstream of AAA model have been used to generate the pulsatility. The flow waveform for the corresponding setup is quite similar to realistic case, but deviates from the exact physiological waveform. A similar gravity driven set up with computer controlled rotating spherical valve downstream of the test section can also be used (Nikolaidis and Mathioulakis, 2002). Relatively enhanced version of the gravity driven flow generator with valves has been developed by Peattie et al. (2004), where two computer-controlled valves, one in the forward direction while the other is in the retrograde direction, were embedded to the system and more realistic waveform with backflow region was obtained with the improved control mechanism compared to Egelhoff et al. (1999) and Nikolaidis and Mathioulakis (2002).

Piston-cylinder arrangements. In piston-driven pumps, a piston is used as linear actuator to provide periodicity to the mean flow that is generated by a gear pump or a head tank located at a certain elevation above the test section, as explained previously. Piston-driven pumps serve better performance in terms of controllability to generate desired flow rate compared to previous orientations, but still physiological waveform, especially backflow regions, cannot be generated properly. In Figure 9, a test rig with a piston-driven pump is presented (Yu, 2000). The mean flow in the test section is generated by a head tank, which is fed by a submersible pump. Pulsatility is generated by means of a pulse generating module, which contains a DC gear-motor fitted with the circular cam of a small piston-cylinder assembly. DC gear-motor drives the piston by means of the circular cam to create necessary pulsating flow conditions (Yu, 2000; Yu and Zhao, 2000). The flow waveform generated in such a flow loop is periodic in sinusoidal form but again not physiological. To obtain more realistic flow waveforms, in later studies, a variable speed electrical motor is used to drive circular cam of linear reciprocating piston (Yip and Yu, 2001, 2002; Stamatopoulos et al., 2010, 2011). In some orientations, motors drive the reciprocating pistons with screw mechanisms such as rack and pinion and lead screw, rather than circular cam (Duclaux et al., 2010; Ene et al., 2011; Morris et al., 2013; Wang et al., 2016). In Figure 10, an example for piston with lead screw mechanism is shown (Tsai and Savaş, 2010). In that case, a gear pump provides the steady mean component of the desired flow waveform and piston arrangement provides the oscillatory component of the flow waveform with the help of a motor driven lead screw. The gear pump is preloaded with a back-pressure valve at its discharge, providing that the gears remain in constant contact to prevent reverse flow. Although the system is capable of generating different pulsatile waveforms, great deviations from physiological cases are reported, especially for the regions of steep changes in the waveform (Tsai and Savaş, 2010). Sometimes two-sided cylinders are also employed with such screw mechanisms. In such a case, a computer controlled motor driven rack mounted piston, divides a cylinder into two parts, and they are connected with a valve which ensures that pump is completely empty on one side of the cylinder while the other side is refilled to keep the flow direction same. The problem for two-sided piston cylinder assembly is its limited operating range in terms of peak flow because of the decreased stroke volume (Frayne et al., 1992; Salsac et al., 2006).

Figure 9

Experimental test set up for AAA test model. The mean flow is generated by a head tank and pulsatility is added by piston cylinder assembly fitted with a circular cam, which is driven by a DC gear-motor [The figure is adapted from Yu (2000) and used with permission].

Figure 10

Sketch and picture of flow circulatory system using gear pump and piston-cylinder arrangement with lead screw as flow generator. Steady flow is supplied by gear pump, while piston and servo motor combination generates the pulsation [The figure is adapted from Tsai and Savaş (2010) and used with permission].

Gear pumps. In recent years, researchers have found that gear pumps are quite appropriate in physiological flow generation, where a controller needs to be integrated with a servo motor to run the gear pump. In that case, reverse flows can also be easily generated. An example of such system is presented in Figure 11 (Mechoor et al., 2016). Different from the aforementioned systems that generate pulsatility, computer-controlled gear pumps operate continuously without requiring any additional equipment and can replicate the physiological waveforms quite accurately (Gaillard and Deplano, 2005; Deplano et al., 2007, 2013, 2014; Mechoor et al., 2016). In addition, Mechoor et al. (2016) developed a closed loop system with integrating a feedback mechanism from output flow to the computer controlled gear pump system to better match the desired physiological flows and pressure waveforms.

Figure 11

Schematic representation of flow circulatory system using computer-controlled gear pump with a feedback mechanism to achieve desired physiological flow and pressure waveform [The figure is adapted from Mechoor et al. (2016) and used with permission].

Commercially available pumps

Besides these component-based in-house built systems, there are commercially available physiological flow generators. Harvard Pump (Harvard Apparatus, Holliston, MA) and SuperPump (Vivitro Labs, Victoria, BC, Canada) are commercially available pumps that are partially programmable and offer several number of arterial waveforms. A fully programmable commercial pump is CardioFlow (Shelley Medical Imaging Technologies, London, ON, Canada), which is again a two-sided piston pump that generates pulsatile flow, by dividing the cylinder into two chambers, each with fluid ports at the end. Some of the studies that utilize these pumps in hemodynamic analysis of different arteries are Pahlevan and Gharib (2013), Groves et al. (2014), and Najjari and Plesniak (2016). Resulting waveforms for these commercial pumps are very similar to physiological cases of predefined arteries, but the range of arteries that they can replicate is limited.

Piping and Pressure Compliance

As a component of flow circulatory systems, piping is critical in terms of obtaining correct inlet waveforms. The main concerns in piping are the general arrangement and the distance between the pump exit and test section inlet, which determine whether the fully developed flow is generated at the inlet of the test section. Some researchers specify this required distance as 100 times the diameter of the tubing before flow enters the test section (Moore et al., 1992). In a latter study, Durst et al. (2005) have suggested a formula, as given in Equation (17).

where D is the tube diameter (m), L is the entrance length (m) and Re is the Reynolds number. Tsai and Savaş (2010) states in the study of cerebral saccular aneurysms that the entrance length required for fully developed flow as given in Equation (18).

He and Ku (1994) have claimed that the maximum length of the inlet pipe to satisfy the requirement of fully developed flow conditions for pulsatile flow is smaller than the tube length for steady flow conditions.

In order to match the pressure waveform at the inlet of the test section, various methods have been implemented. Researchers frequently employing downstream vascular simulators, replicate the impedances along the arteries to generate pressure waveform at specific arteries, which are generally composed of flow resistances and compliance chambers. In some recent studies, resistor–capacitor–resistor (RCR) module is adapted, which is composed of proximal and distal resistors (pinch valves) with a compliance chamber (air chamber) (Deplano et al., 2013; Mechoor et al., 2016). Besides all of these, reproducing physiological pressures throughout the aneurysm sac can be accomplished very accurately by the lumped-parameter boundary condition module, which is suggested by Kung et al. (2011). Such module is adapted downstream of aneurysm model to replicate outlet boundary condition, consisting of an inductance (L), proximal resistance (Rp), capacitance (C), and distal resistance (Rd), which is, in total, called as four-element Windkessel model.

Test Section

Another component of a circulatory flow setup is the model of the artery of interest, or so called, anatomical phantom. In early studies, an aneurysm model, which has transparent rigid walls with simplified geometry, has been integrated to flow loop (Egelhoff et al., 1999). These simplified and rigid models cannot replicate the exact compliant nature of the vessels. However, they provide insight about the governing physical processes occurring inside the AAA, including but not limited to separation of the flow at the inlet of AAA, vortex generation and jet regions, transition to turbulence, oscillatory flow structure, and high or low wall shear stress regions. With such controllable geometries, the parametric studies regarding aspect ratio, bulge diameter, and bulge asymmetry can be characterized. In recent investigations, more realistic geometries have been implemented. Rigid aneurysm models with simplified geometries are in the form of straight tube with a concentric bulge, which can be described by an ellipse formula (Yu, 2000), and are generally manufactured from glass (Salsac et al., 2006), or resins (Biglino et al., 2013). On the other hand, in order to understand exact hemodynamics and biomechanical forces that are generated through the aneurysms, it is necessary to work with compliant arterial models (Deplano et al. (2007) to replicate in vivo condition by manufacturing patient-specific phantom.

Manufacturing the realistic and patient specific arterial model is quite challenging. Different techniques are applied in literature by using different materials like silicone, polyurethane, and latex, where in the most well-known technique polyjet printing is utilized to manufacture transparent, 3D, and compliant phantoms (Sulaiman et al., 2008; Biglino et al., 2013; Ionita et al., 2014). 3D patient-specific vascular imaging data generated by means of CT scanner or MR techniques are collected and exported into Standard Tessellation Language (STL) file format, and this STL file is imported into 3D polyjet printer. The printer head generates layers of liquid photopolymer and builds them up to create 3D models. This technique enables to utilize different materials with a high resolution. As another method 3D printed vessel lumens are used as molds of the real phantoms that are going to be casted (Ho et al., 2017). For non-transparent phantoms, other techniques with non-transparent materials are also available to manufacture patient specific artery models where Doppler Ultrasound measurements are conducted for flow quantification. As an example, embedding silicone-elastomer vessel model into the agar based tissue-mimicking material (TMM) having same acoustic properties with artery (Poepping et al., 2004), direct machining of phantom using numerically controlled milling rather than producing a mold (Wong et al., 2008) or some other techniques (Watts et al., 2007; Allard et al., 2013) can be given. Even though Doppler Ultrasound technique does not require transparent phantoms, which may be an advantage, it provides only the maximum velocity in a flow cross section. Therefore, detailed hemodynamics evaluation is not possible with Doppler Ultrasound technique.

Blood-Mimicking Fluid

Blood is a non-Newtonian fluid, where the viscosity changes with shear rate. In literature, most researchers have utilized Newtonian fluids as blood-mimicking fluid because blood behaves like a Newtonian fluid for certain shear rate values (Berger and Jou, 2000). However, its non-Newtonian behavior becomes critical when shear rate is low and Deplano et al. (2013) reports that wall shear rate in some regions of AAA drops to values of 1 s−1. In this case, the blood acts as a non-Newtonian shear thinning fluid because viscosity decreases with increasing shear rate (Mandal, 2005). Non-Newtonian behavior of blood can be modeled using Carreau-Yasuda viscosity model, which is given in Equation (14). Most of the experimental studies utilize the Newtonian assumption, and generally pure water or the mixture of glycerin and water with different volumetric ratios have been widely used as blood mimicking fluid (Boutsianis et al., 2008; Stamatopoulos et al., 2010, 2011). There is only few studies that use non-Newtonian fluids of which behavior are very similar to the blood, like Xanthane Gum solution with glycerin (Deplano et al., 2014; Najjari and Plesniak, 2016).

Flow Measurement System

In early studies, flow field is generally observed with the help of flow visualization technique (Budwig, 1994; Moore and Ku, 1994). The growing trend, which is the Particle Image Velocimetry (PIV) has been developed since the early 1980's to acquire velocity vector information of a whole flow field instantaneously with a high resolution (Adrian, 1991). Since 1993, PIV has been utilized in biologically important flows because of being a non-intrusive approach. In terms of the flow field information, different versions of PIV are available including 2D (planar) PIV, Stereoscopic PIV, Holographic PIV, and Tomographic PIV. In most of the AAA studies, planar PIV has been implemented, while Stereoscopic PIV has been utilized in a very few studies. As of authors knowledge, although Tomographic PIV has been applied to disturbed hemodynamics through intracranial aneurysms (Roloff et al., 2017, 2018), Holographic PIV and Tomographic PIV have not been applied to AAA studies yet.

Planar PIV is the mapping of average displacements of seeding particles within interrogation areas over a small time interval by means of two successive images of an illuminated plane in a fluid flow. Interrogation areas are small sub areas of each image, containing sufficient number of illuminated seeding particles. To find the displacements of particles for a short time interval, interrogation areas of two successive images are compared by means of cross correlation technique, resulting in determination of the most probable velocity vector for the particles inside of that interrogation areas, and this procedure continues until all interrogation areas are correlated and whole flow field velocity vectors are generated.

The fluid is seeded with neutrally buoyant particles, which should be sufficient in amount and size, typically varying from 5 to 50 μm, depending on the flow rate, magnification of the camera, and the field of view as well (Stamhuis, 2006). There are several types of seeding particles with respect to illumination characteristics, such as reflective, scattering, and fluorescent. Among all of these options, fluorescent particles have an advantage that fluorescent light can be distinguished from the illumination, which increases the visibility of the particles (Tsai and Savaş, 2010). Illumination of seeding particles with laser sheets and capturing the images with camera bring the refraction problem where the refractive indices of mediums and inclination angles of surfaces when lights passes through them are quite critical. Since the AAA models are curved in shape, in addition to the efforts given to match the refractive indices of working fluid and phantom, a box with flat surfaces perpendicular to the surface normal of camera lens plane covers the test section and is filled with working fluid to minimize refraction problems (Budwig, 1994).

In planar PIV, two component velocity information on 2-D plane is obtained using two laser sheets and an advanced camera, where in Stereoscopic PIV three-velocity components are obtained on 2-D plane using two cameras, as can be seen in Figure 12. Note that, lens planes and image planes of cameras are tilted to satisfy the Scheimpflug arrangement and plexiglass prisms in addition to the box for the test section are used to minimize refraction problems. In literature, there are many studies that utilize PIV technique to study the hemodynamics through AAA. In Table 1, PIV studies that have been conducted for the last 20 years are listed, considering the flow conditions at the inlet, type of blood mimicking fluids, wall types, geometry, and Reynolds number ranges.

Figure 12

Typical camera arrangement of the SPIV. δ is the angle between the lens plane and the image plane, and faces of plexiglass prism is perpendicular to corresponding camera, which is used to minimize refraction related problems [The figure is adapted from Deplano et al. (2016) and used with permission].

Table 1

PIV studies conducted in the field of hemodynamics in AAA for the last 20 years.

References	Flow generator	Blood mimicking fluid	Wall type	Geometry	Reynolds numbers
Yu, 2000	Gravity driven pump with piston cylinder arrangement (Unsteady-Sinusoidal)	Newtonian fluid: A solution mixture of glycerin and water	Rigid wall-Pyrex glass tubes	Straight Tube and Axisymmetric, elliptical shaped bulge	Steady Flow: 400–1400 Unsteady Flow, Peak Value: 1,274
Salsac et al., 2006	Two sided piston cylinder arrangement (Unsteady-Physiological)	Newtonian fluid: Pure water	Rigid wall-Glass	Axisymmetric bulges with different diameters	Unsteady flow, Peak Value: 2,700
Deplano et al., 2007	Computer controlled pump (Unsteady-Physiological)	Newtonian fluid: Aqueous glycerin solution, 60% water	Rigid wall-Glass Compliant wall-Molded polyurethane	Asymmetric bulge	–
Boutsianis et al., 2008	Variable speed gear pump (Steady)	Newtonian fluid: A mixture of 40% water, 60% glycerol	Compliant wall-Silicone phantom	Patient-specific aneurysm, gathered by CT scanning	Steady flow: 560
Stamatopoulos et al., 2010	Linear reciprocating piston cylinder arrangement (Unsteady-Sinusoidal)	Newtonian fluid: A water and glycerin solution (40:60 by volume)	Rigid wall-Elastomer material (Sylgard-184)	Axisymmetric, elliptical shaped bulge	Steady Flow: 105–690 Unsteady Flow, 105–690
Stamatopoulos et al., 2011	Linear reciprocating piston cylinder arrangement (Unsteady-Physiological)	Newtonian fluid: A water and glycerin solution (40:60 by volume)	Compliant wall-Liquid silicon elastomer	Patient-specific aneurysm, gathered by CT scanning	Unsteady flow, Peak Value: 541
Deplano et al., 2013	Computer controlled pump (Unsteady-Physiological)	Newtonian fluid: Aqueous glycerin solution, 60% water	Compliant wall-Molded polyurethane	Asymmetric bulge with symmetric and asymmetric iliac bifurcation	Unsteady Flow, Peak Value: 1,876
Chen et al., 2014	Steady Submersible Pump	Newtonian fluid: A mixture of NaI and water	Rigid wall	Patient-specific aneurysm, produced by rapid prototyping	Steady Flow: 2,234
Deplano et al., 2014	Computer controlled pump (Unsteady-Physiological)	Shear-thinning fluid: Aqueous solution of Xanthane gum (XG)	Compliant wall-Molded polyurethane	Asymmetric bulge with asymmetric iliac bifurcation	Unsteady Flow, Peak Value: 1,941
Deplano et al., 2016	Computer controlled pump (Unsteady-Physiological)	Shear-thinning fluid: Aqueous solution of Xanthane gum (XG)	Compliant wall-Molded polyurethane	Asymmetric bulge with asymmetric iliac bifurcation	Unsteady Flow, Peak Value: 2,298
Wang et al., 2016	Piston cylinder arrangement (Unsteady-Physiological)	Newtonian fluid: Pure water	Rigid wall-Glass	Straight Tube and Axisymmetric, elliptical shaped bulge

Scaling Parameters

Pulsatile flow through the arteries is generally characterized by two non-dimensional parameters, Reynolds and Womersley numbers, which are defined as the Equations (19,20).

where U is the velocity, D is the tube diameter, ν is the kinematic viscosity and ω is the frequency of periodicity (Womersley, 1955). Reynolds number is the ratio of inertia forces to viscous forces, which determines the flow regime and varies from 1 for very small arteries to ~4,000 for aorta. Womersley number is the ratio of unsteady inertia forces to viscous forces, which can also be represented in terms of Reynolds number and Strouhal number. For low Womersley number flows, viscous forces dominate and velocity profile inside the arteries become parabolic, and velocity at the centerline oscillates with a phase of driving pressure gradient. For high Womersley number flows (α > 10), unsteady inertia forces dominate and velocity profile has a nearly flat velocity profile, and in that case, flow represents situations with rapid acceleration and deceleration (Womersley, 1955; Ku, 1997). Womersley number value for femoral artery is nearly four, while for aorta, it reaches up to 22.

These non-dimensional parameters are required to be kept constant when scaling is conducted, enabling to perform similarity analysis. The artery models that are studied during the experiments might be scaled to a larger size because generally the vessels in human circulatory system are relatively small for macro scale investigations. In that case, the normal vessel is the prototype with exact dimensions and velocity values, while the artery model utilized in experiments might be the scaled version of the prototype. In that case, performing a similarity analysis by matching Reynolds and Womersley numbers of prototype and model will mimic physiological conditions and enables the simulation of performance of the prototype. By equating the Reynolds and Womersley numbers for the model and prototype, as given in Equations (21,22).

where subscript m is representing model and p is representing prototype and with the following ratio of diameters of arteries and kinematic viscosities, which is given in Equation (23).

where L is diameter ratio and ν is the kinematic viscosity ratio of prototype to model, the velocity and frequency of the model become

To increase the resolution of a small vessel by increasing its diameter, the velocity should be decreased by the ratio of 1/L to match the Reynolds number, while frequency should be decreased by the ratio of to match the Womersley number. In certain test setups, the frequency adjustment might be needed possibly due to the following two reasons: (1) Flow circulatory system could not reach the desired frequencies for the flow waveform, (2) Flow measurement system could not provide enough sample in one cycle. In that case, the desired velocity and diameter for the setup are determined accordingly. Similarly, if the flow rate in the setup is the limitation, then once the velocity is determined, diameter of the setup is adjusted according to the Reynolds number and finally the frequency of the waveform is set based on the Womersley number to ultimately reach the similarity.

Key Parameters to Study in Experimental Investigation of Hemodynamics of AAA

In experimental studies, researchers have generally reported the results of the hemodynamics in AAA in terms of velocity and vorticity parameters (Salsac et al., 2006; Stamatopoulos et al., 2010) obtained using several vortex identification techniques, such as swirling strength (Deplano et al., 2007). Swirling strength of a vortex is defined by λ, while is called enstrophy, which is energy of vorticity and utilized for decreasing background noise (Zhou et al., 1999).

In recent years, with the utilization of Stereoscopic PIV technique, some researchers have also presented the 3D vortex evaluation through aneurysm sac (Deplano et al., 2016). In this technique, Deplano et al. (2016) have focused on the effect of transverse velocity component (w) measurement by means of comparing its magnitude to other velocity components, such as axial (u), vertical (v), or total velocity vector (V), of which magnitude is

In addition, the magnitude of 2D velocity vector, (V), can also be written as

For each spatial point of different planes inside aneurysm, relative difference, RD between these two velocity magnitudes is calculated as

Once the RD values are calculated at each point, the mean value of the relative difference, RD, can be obtained on the corresponding planes. If RD, is lower than 10%, it can be interpreted as u and v components are sufficient to represent flow behavior, whereas for higher mean relative difference values, w component needs to be included (Deplano et al., 2016). The quantities such as contours and , , RD values, which underline the specific importance of w with respect to u and v components, can be utilized in order to investigate the importance of transverse velocity component (w) on hemodynamics and AAA progression. In addition, vortex ring and its strain dynamics which can also be obtained from Stereoscopic PIV data, are used to understand different flow instabilities and growth mechanism (Deplano et al., 2016). The rate of deformation tensor can be decomposed into rate of strain tensor, which is pure rate of deformation tensor S, and pure rate of rotation tensor Ω. The strain rate tensor S enables to visualize and quantify the direction of compression and stretching. If the eigenvalues of tensor S are positive, the fluid element is stretched while they are negative, the element is compressed (Bouremel et al., 2009).

Important Findings on Hemodynamics of AAAs From Experimental Studies

Hemodynamics in Undilated Aorta

In early studies, researchers have focused on hemodynamics inside undilated aorta to understand the effect of flow on dilation of arterial wall. Ku et al. (1989) investigated undilated abdominal aorta by using a realistic glass abdominal aorta model and observed steady flow inside that artery by injecting a dye into the model. Flow pattern at the infrarenal aorta, where AAA generation is observed, has more complex structure than in the suprarenal aorta because of the following reasons: (1) The aorta has a curvature at that section, which localizes transient separation to its posterior wall; (2) There are branches of arteries transferring blood to kidneys after the suprarenal and before the infrarenal segments. The blood flow transferred from these branches creates secondary flows at the upstream of infrarenal aorta; (3) The bifurcation at the downstream of infrarenal aorta is creating a horseshoe vortex. All these unsteady flow structures make the infrarenal segment appropriate to generate an aneurysmal bulge.

Hemodynamics in AAA

Effect of pulsatility

In early studies on AAA (Nerem, 1984; Egelhoff et al., 1999; Yu, 2000), steady flow analysis has been conducted and as reported under laminar conditions, a jet of fluid passing though the core of the aneurysm is surrounded by a circulating vortex. The flow velocity through the center of the aneurysm is 40 times higher than the velocity at the recirculation region, and WSS value has the largest peak at the reattachment point of this recirculating flow. The WSS magnitude in the recirculation region is about ten times less than the value in the upstream tube, because the core flow separation at the bulge causes the negative WSS values. Scherer (1973) reports that the transition to turbulence begins at Re = 2900. Asbury et al. (1995) studied AAA without any bifurcation and utilized Color Doppler Flow Imaging (CDFI) method to visualize the flow and Laser Doppler Velocimetry (LDV) method to quantify the flow in order to address the effect of bulge diameter and turbulence on hemodynamics. They reported a correlation between the aneurysm size and rupture because especially in larger aneurysms turbulent flow is observed, which causes higher WSS values throughout the bulge (Asbury et al., 1995).

The cycle of physiological flow with strong acceleration and deceleration of flow in AAA has strong effects on vortex dynamics, which moves from proximal to the distal end creating significant fluctuations in WSS and pressure though the aneurysm bulge. Fukushima et al. (1989) conducted the first experimental and numerical study on pulsatile flow conditions for glass models of axisymmetric and asymmetric AAAs. They observed that primary and secondary vortices were present inside the bulge, and reported that the peak negative shear stress occurs at the distal end. Yu (2000) performed the first PIV experiments for both steady and pulsatile flow waveforms in pyrex glass aneurysm models without any iliac bifurcation to see the effects on sinusoidal waveform on hemodynamics. Yu and Zhao (2000) compared the steady and pulsatile flow dynamics in a straight tube with a side bulge using PIV, and they reported significant differences between flow patterns. Stamatopoulos et al. (2010) reported the difference between steady and pulsatile, sinusoidal flow in AAA, where the flow separation at the proximal end and reattachment at the distal end of aneurysm bulge is witnessed for steady case, but for unsteady case, these locations vary for each instant of cardiac cycle. However, the peak WSS values are observed at the model exit for both flow conditions.

Effects of wall material properties, working fluid properties, and geometric parameters on AAA hemodynamics

In most of the experimental studies, rigid axisymmetric models were used to investigate the hemodynamics in aneurysm bulge (Yu, 2000; Salsac et al., 2006; Stamatopoulos et al., 2010). Deplano et al. (2007) and Meyer et al. (2011) have tested two aneurysm models one having rigid and the other with compliant walls. The rigid aneurysm was made of glass while the compliant one was made of molded polyurethane, and both had the identical geometry as a simple bulge for comparison purposes. As can be seen in Figure 13, there is significant difference between the flow fields inside both of the aneurysms. In rigid model, there is only viscous dissipation, while in compliant model, viscoelastic dissipation is also observed, which is indicated as more critical. During the accelerating phase of the cycle, compliant walls absorb kinetic energy in the potential energy form, leading the walls become expanded. While the flow is decelerating, this stored kinetic energy retracts the walls. This expansion and retraction contribute the progression of vortices at the distal end of aneurysm during the deceleration phase. They conclude that increasing wall compliance causes collision of vortices with the walls, increasing both the local pressures and wall stresses especially at the distal end (Deplano et al., 2007; Meyer et al., 2011). Producing models with exact in vivo wall material properties and thicknesses are also important to decide rupture location. In literature, there are several studies focusing on identifying rupture location in compliant AAA models by means of tensile tests (Doyle et al., 2009a, 2010) and biaxial tests (O'Leary et al., 2014). Mechanical behaviors of ILT and wall calcification and their effects on rupture have also been characterized by means of these tests (Ene et al., 2011; O'Leary et al., 2013, 2015).

Figure 13

Temporal evolution of the swirling strength λ at different time instants in (A) the rigid model and the compliant model in horizontal plane, and (B) the compliant model at horizontal and vertical planes, under exercise conditions. During the deceleration phase for the compliant model, vortices impact on the walls and swirling strength increases, including vortex shedding occurrence [The figure is adapted from Deplano et al. (2007) and used with permission].

The real aneurysms also contain iliac bifurcation. Deplano et al. (2013) have studied asymmetric aorto-iliac bifurcation for a simple bulge. Two artery models, one with a straight outlet tube and the other with iliac bifurcation have been studied to observe the effect of bifurcation on flow structure. The results indicate that for the model with iliac bifurcations, intensity of the vortex ring impact on the anterior wall is about 90% higher than without bifurcation. When the vortex ring impinges on the wall, the forces generated by this impingement might possibly cause the rupture (Chu et al., 1995). For a compliant aneurysm bulge with aorta iliac bifurcation, stereoscopic PIV technique has been implemented to observe vortex ring impingement inside the aneurysm as 3D quantification (Deplano et al., 2016). In Figure 14, contours of in six planes and four time instants of cardiac cycle is represented in order to understand the effect of transverse velocity component, w. During the deceleration phase, where T = 0.45, 0.55 and 0.85, the magnitude of the mean transverse velocity component, w, is 0.7 times the axial component and 1.7 times the vertical component, especially for C and D planes which are located within the AAA bulge.

Figure 14

Contours of from Stereoscopic PIV measurements in six planes and four time instants of cardiac cycle (A) T = 0.3; maximum flow rate, (B) T = 0.45; decelerating flow rate, (C) T = 0.55; minimum flow rate, (D) T = 0.85; flow rate is nearly zero. Velocity vector projection and isosurface λ = 8 (in white) is imposed on each plane. Flow stagnation area (in black) is represented by weak components of velocity vectors [The figure is adapted from Deplano et al. (2016) and used with permission].

The other important concern about AAA is intraluminal thrombus (ILT) developed inside the aneurysms, consisting of blood proteins, platelets, and cells. Platelets stick the vortex ring generated inside the aneurysm throughout the cardiac cycle, and are transported from proximal to distal end, released during vortex break up. They adhere to low wall shear stress sites, which will lead to thrombus formation (Biasetti et al., 2011; O'Rourke et al., 2012). The question that ILT formation decreases the rupture risk or not still has controversy. In an experimental and numerical study on AAA with and without ILT, it is reported that a recirculation region is generated inside AAA without thrombus while it is disappeared in AAA with ILT (Chen et al., 2014), and there are several studies that are seeking the governing hemodynamic parameters affecting platelet deposition (Deplano et al., 2013, 2014).

As the patient-specific studies, there are a few examples such as Particle Tracking Velocimetry (PTV) study of Boutsianis et al. (2008) and PIV study of Stamatopoulos et al. (2011), where they report that WSS have the peak values at the proximal and distal sites, like in simplified bulges, and the velocity waveform at the anterior site is highly disturbed because the vortex is generated and transmitted to downstream at this site, rather than posterior.

For mimicking the shear stress and shear strain rate relation of blood, Deplano et al. (2014) have used Xanthane Gum dissolved in aqueous solutions of glycerol, a shear thinning fluid, where they report that shear thinning fluid model imposes higher shear stress values throughout the aneurysm bulge than Newtonian model. They conclude that the rheology of working fluid affects the hemodynamics inside aneurysm model, and should be taken into account for more realistic experimentation. On the other hand, in another study which is performed on a curved artery model with again a shear thinning fluid, mixture of fluids like Xanthane Gum, glycerin, water, and sodium iodide (Najjari and Plesniak, 2016). The results are very similar to the studies for curved arteries with Newtonian fluids, concluding that for large arteries with large flow rates, rheology of the blood does not affect the flow structure. However, at that point it is important to note that, they have utilized comparatively high sodium iodide concentration in their mixture than required to resolve the refractive index problems, leading that their working fluid is less viscoelastic than the blood.

Conclusion

Experimental and computational studies indicate that a certain number of parameters can be interpreted as the indicators of the rupture. Typically, the peak values of WSS and wall stress (structural stress) are interpreted as the footprints of an upcoming rupture. Some key parameters that are derived from WSS, such as TAWSS, OSI, and ECAP are also utilized to detect candidate locations of the rupture (Arzani and Shadden, 2015; Arzani, 2018; Singh et al., 2018). Unfortunately, determining the exact location of rupture by deciding the most effective hemodynamic parameter is still a controversial issue. For improved rupture risk assessment, ILT and plaque formations, patient-specific boundary conditions, arterial thickness variation, and heterogeneous material properties of arterial wall should be considered in a comprehensive computational analysis. For the experimental studies, flow loop should exactly replicate the in-vivo characteristics using a programmable, real time physiological pump that can supply re-producible and realistic inflow conditions to test model (Mechoor et al., 2016). In addition, lumped-parameter outlet boundary condition modules should be used to produce physiological pressures together with a compliant artery model having a patient specific geometry. The main purpose of these studies, in the end, should be to obtain a clinically useful tool in decision making step to understand whether the repair is necessary or not. In addition to mechanical analyses, the genetic aspect of the problem should be taken into account, since the disturbed biomechanical environment inside the aneurysm may affect gene expression patterns which may result in altered growth at later stages (Krishna et al., 2010).

Author Contributions

HS and BR wrote the first draft and revised the manuscript. MY and HY structured, reviewed, and revised the manuscript. All authors read and approved the submitted version of the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.