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

On the feasibility of ultrasound Doppler-based personalized hemodynamic modeling of the abdominal aorta

BACKGROUND

Personalized modeling is a promising tool to improve abdominal aortic aneurysm (AAA) rupture risk assessment. Computed tomography (CT) and quantitative flow (Q-flow) magnetic resonance imaging (MRI) are widely regarded as the gold standard for acquiring patient-specific geometry and velocity profiles, respectively. However, their frequent utilization is hindered by various drawbacks. Ultrasound is used extensively in current clinical practice and offers a safe, rapid and cost-effective method to acquire patient-specific geometries and velocity profiles. This study aims to extract and validate patient-specific velocity profiles from Doppler ultrasound and to examine the impact of the velocity profiles on computed hemodynamics.

METHODS

Pulsed-wave Doppler (PWD) and color Doppler (CD) data were successfully obtained for six volunteers and seven patients and employed to extract the flow pulse and velocity profile over the cross-section, respectively. The US flow pulses and velocity profiles as well as generic Womersley profiles were compared to the MRI velocities and flows. Additionally, CFD simulations were performed to examine the combined impact of the velocity profile and flow pulse.

RESULTS

Large discrepancies were found between the US and MRI velocity profiles over the cross-sections, with differences for US in the same range as for the Womersley profile. Differences in flow pulses revealed that US generally performs best in terms of maximum flow, forward flow and ratios between forward and backward flow, whereas it often overestimates the backward flow. Both spatial patterns and magnitude of the computed hemodynamics were considerably affected by the prescribed velocity boundary conditions. Larger errors and smaller differences between the US and generic CFD cases were observed for patients compared to volunteers.

CONCLUSION

These results show that it is feasible to acquire the patient-specific flow pulse from PWD data, provided that the PWD acquisition could be performed proximal to the aneurysm region, and resulted in a triphasic flow pattern. However, obtaining the patient-specific velocity profile over the cross-section using CD data is not reliable. For the volunteers, utilizing the US flow profile instead of the generic flow profile generally resulted in improved performance, whereas this was the case in more than half of the cases for the patients.

Effect of Intraluminal Thrombus Burden on the Risk of Abdominal Aortic Aneurysm Rupture

Abdominal aortic aneurysm (AAA) is a critical health disorder, where the abdominal aorta dilates more than 50% of its normal diameter. Enlargement in abdominal aorta alters the hemodynamics and flow-induced forces on the AAA wall. Depending on the flow conditions, the hemodynamic forces on the wall may result in excessive mechanical stresses that lead to AAA rupture. The risk of rupture can be predicted using advanced computational techniques such as computational fluid dynamics (CFD) and fluid-structure interaction (FSI). For a reliable rupture risk assessment, formation of intraluminal thrombus (ILT) and uncertainty in arterial material properties should be taken into account, mainly due to the patient-specific differences and unknowns in AAAs. In this study, AAA models are computationally investigated by performing CFD simulations combined with FSI analysis. Various levels of ILT burdens are artificially generated in a realistic AAA geometry, and the peak effective stresses are evaluated to elucidate the effect of material models and ILT formation. The results indicate that increasing the ILT burden leads to lowered effective stresses on the AAA wall. The material properties of the artery and ILT are also effective on the stresses; however, these effects are limited compared to the effect of ILT volume in the AAA sac.

The Impact of a Limited Field-of-View on Computed Hemodynamics in Abdominal Aortic Aneurysms: Evaluating the Feasibility of Completing Ultrasound Segmentations with Parametric Geometries

To improve abdominal aortic aneurysm (AAA) rupture risk assessment, a large, longitudinal study on AAA hemodynamics and biomechanics is necessary, using personalized fluid-structure interaction (FSI) modeling. 3-dimensional, time-resolved ultrasound (3D+t US) is the preferred image modality to obtain the patient-specific AAA geometry for such a study, since it is safe, affordable and provides temporal information. However, the 3D+t US field-of-view (FOV) is limited and therefore often fails to capture the inlet and aorto-iliac bifurcation geometry. In this study, a framework was developed to add parametric inlet and bifurcation geometries to the abdominal aortic aneurysm geometry by employing dataset statistics and parameters of the AAA geometry. The impact of replacing the patient-specific inlet and bifurcation geometries, acquired using computed tomography (CT) scans, by parametric geometries was evaluated by examining the differences in hemodynamics (systolic and time-averaged wall shear stress and oscillatory shear index) in the aneurysm region. The results show that the inlet geometry has a larger effect on the AAA hemodynamics (median differences of 7.5 to 18.8%) than the bifurcation geometry (median differences all below 1%). Therefore, it is not feasible to replace the patient-specific inlet geometry by a generic one. Future studies should investigate the possibilities of extending the proximal FOV of 3D+t US. However, this study did show the feasibility of adding a parametric bifurcation geometry to the aneurysm geometry. After extending the proximal FOV, the obtained framework can be used to extract AAA geometries from 3D+t US for FSI simulations, despite the absence of the bifurcation geometry.

Prediction of wall stress and oxygen flow in patient-specific abdominal aortic aneurysms: the role of intraluminal thrombus

In this study, the biomechanical role of intraluminal thrombus (ILT) in an abdominal aortic aneurysm (AAA) is investigated. The implications of ILT in AAA are controversial in literature. Previous studies have demonstrated that ILT provides a biomechanical advantage by decreasing wall stress, whereas other studies have associated ILT with inhibiting oxygen transport and inducing aortic wall weakening. Therefore, we sought to explore the connection between ILT, mechanical stresses, and oxygen flow in different geometries of patient-specific aneurysms with varying ILT morphologies. The objective is to investigate the extent to which ILT influences the prediction of aneurysmal wall stresses that are associated with rupture, as well as oxygen concentrations to measure tissue oxygen deprivation. Three patient-specific AAA geometries are considered, and two models, one with ILT and one without ILT, are created for each patient to assess the effect of ILT presence. A fluid-structure interaction approach is used to couple the blood flow, wall deformation, and oxygen mass transport. Results are presented for hemodynamics patterns, wall stress measures, and oxygen metrics within the arterial wall. While ILT is found to reduce wall stress, simulations confirm that ILT decreases oxygen transport within the tissue significantly, leading to wall hypoxia.

Numerical analysis of biothermal-fluids and cardiac thermal pulse of abdominal aortic aneurysm

Abdominal aortic aneurysms are serious and difficult to detect, conditions can be deadly if they rupture. In this study, the heat transfer and flow physics of Abdominal Aortic Aneurysm (AAA) were discussed and associated with cardiac cycle to illustrate the cardiac thermal pulse (CTP) of AAA. A CTP and infrared thermography (IRT) evaluation-based on AAA and abdomen skin surface detection method was proposed, respectively. Infrared thermography (IRT) is a promising imaging technique that may detect AAA quicker and cheaper than other imaging techniques (as biomarker). From CFD rigid-wall and FSI Analysis, the transient bioheat transfer effect resulted in a distinct thermal signature (circular thermal elevation) on the temperature profile of midriff skin surface, at both regular body temperature and supine position, under normal clinical temperature. However, it is important to note that thermography is not a perfect technology, and it does have some limitations, such as lack of clinical trials. There is still work to be done to improve this imaging technique and make it a more viable and accurate method for detecting abdominal aortic aneurysms. However, thermography is currently one of the most convenient technologies in this field, and it has the potential to detect abdominal aortic aneurysms earlier than other techniques. CTP, on the other hand, was used to examine the thermal physics of AAA. In CFD rigid-wall Analysis, AAA had a CTP that only responded to systolic phase at regular body temperature. In contrast, a healthy abdominal aorta displayed a CTP that responded to the full cardiac cycle, including diastolic phase at all simulated cases. Besides, the findings from FSI Analysis suggest the influence of numerical simulation techniques on the prediction of thermal physics behaviours of AAA and abdominal skin surface. Lastly, this study correlated the relationship between natural convective heat transfer coefficient with AAA and provided reference for potential clinical diagnostic using IRT in clinical implications.

Assessment of the systolic rise time by photoplethysmography in peripheral arterial diseases: a comparative study with ultrasound Doppler

Aims

Peripheral arterial disease (PAD) is a major public health burden requiring more intensive population screening. Ankle brachial index (ABI) using arm and ankle cuffs is considered as the reference method for the detection of PAD. Although it requires a rigorous methodology by trained operators, it remains time-consuming and more technically difficult in patients with diabetes due to mediacalcosis. Techniques based on the study of hemodynamic, such as the systolic rise time (SRT), appear promising but need to be validated. We retrospectively compared the reliability and accuracy of SRT using a photoplethysmography (PPG) technique to the SRT measured by ultrasound doppler (UD) in PAD patients diagnosed with the ABI (137 patients, 200 lower limbs).

Methods and results

There was a significant correlation between SRT measured with UD (SRTud) compared with that with PPG (SRTppg, r = 0.25; P = 0.001). Best correlation was found in patients without diabetes (r = 0.40; P = 0.001). Bland and Altman analysis showed a good agreement between the SRTud and SRTppg. In contrast, there was no significant correlation between UD and PPG in diabetes patients. Furthermore, patients with diabetes exhibited a significant increase of SRTppg (P = 0.02) compared with patients without diabates but not with the SRTud (P = 0.18). The SRTppg was significantly linked to the arterial velocity waveforms, the type of arterial lesion but not vascular surgery revascularization technique.

Conclusion

This monocentric pilot study shows that SRT measured with the PPG signal reliably correlates with SRT recorded with UD. The PPG is an easy to use technique in the hand of non-expert with a potential interest for general screening of PAD, especially in diabetes patients, due to its ease to use.

Micromechanical and Ultrastructural Properties of Abdominal Aortic Aneurysms

Abdominal aortic aneurysms are a common condition of uncertain pathogenesis that can rupture if left untreated. Current recommended thresholds for planned repair are empirical and based entirely on diameter. It has been observed that some aneurysms rupture before reaching the threshold for repair whilst other larger aneurysms do not rupture. It is likely that geometry is not the only factor influencing rupture risk. Biomechanical indices aiming to improve and personalise rupture risk prediction require, amongst other things, knowledge of the material properties of the tissue and realistic constitutive models. These depend on the composition and organisation of the vessel wall which has been shown to undergo drastic changes with aneurysmal degeneration, with loss of elastin, smooth muscle cells, and an accumulation of isotropically arranged collagen. Most aneurysms are lined with intraluminal thrombus, which has an uncertain effect on the underlying vessel wall, with some authors demonstrating a reduction in wall stress and others a reduction in wall strength. The majority of studies investigating biomechanical properties of ex vivo abdominal aortic aneurysm tissues have used low-resolution techniques, such as tensile testing, able to measure the global material properties at the macroscale. High-resolution engineering techniques such as nanoindentation and atomic force microscopy have been modified for use in soft biological tissues and applied to vascular tissues with promising results. These techniques have the potential to advance the understanding and improve the management of abdominal aortic aneurysmal disease.

3D finite-element modeling of vascular adaptation after endovascular aneurysm repair

Aneurysm shrinkage is clinically observed after successful endovascular aortic aneurysm repair (EVAR). However, global understanding of post-operative aneurysm evolutions remains weak. In this work, we propose to study these effects using numerical simulation. We set up a 3D finite-element model of post-EVAR vascular adaptation within an open-source finite-element code, which was initially developed for growth and remodeling (G&R). We modeled the endograft with a set of uniaxial prestrained springs that apply radial forces on the inner surface of the artery. Constitutive equations, momentum balance equations, and equations related to the mechanobiology of the artery were formulated based on the homogenized constrained mixture theory. We performed a sensitivity analysis by varying different selected parameters, namely oversizing and compliance of the stent-graft, gain parameters related to collagen G&R, and the residual pressure in the aneurysm sac. This permitted us to evaluate how each factor influences post-EVAR vascular adaptation. It was found that oversizing, compliance or gain parameters have a limited influence compared to that of the residual pressure in the aneurysm sac, which was found to play a critical role in the stability of aneurysm after stent-graft implantation. An excessive residual pressure larger than 50 mmHg can induce a continuous expansion of the aneurysm while a moderate residual pressure below this critical threshold yields continuous shrinkage of the aneurysm. Moreover, it was found that elderly patients, with relatively lower amounts of remnant elastin in the arterial wall, are more sensitive to the effect of residual pressure. Therefore, these results show that elderly patients may present a higher potential risk of aortic sac expansion due to intra-aneurysm sac pressure after EVAR than younger patients.

Hemodynamic simulation of abdominal aortic aneurysm on idealised models: Investigation of stress parameters during disease progression

BACKGROUND AND OBJECTIVE

Analysis and prediction of rupture risk of abdominal aortic aneurysms (AAA) facilitates planning for surgical interventions and assessment of plausible treatment modalities. Present approach of using maximum diameter criterion, is giving way to hemodynamic and bio-mechanical based predictors in conjunction with Computational fluid dynamic (CFD) simulations. Detailed studies on hemodynamic and bio-mechanical parameters at the stage of maximum growth/rupture is of practical importance to the clinical community. However, understanding the changes in these parameters at different stages of growth, will be useful for clinicians, in planning routine monitoring to reduce the risk of sudden rupture. This is particularly useful in medical resource starved nations. Present study investigates the hemodynamic and bio-mechanical changes occurring during the growth stages of aortic aneurysms using fluid structure interaction (FSI) studies.

METHOD

Six idealized fusiform aneurysm models spanning high (shorter) and low (longer) values of the shape index (DHr), have been analysed at three different stages of growth viz, a D_max of 3.5cm, 4.25cm, 5cm. Pulsatile Newtonian blood flow, passing through an elastic arterial vessel wall with uniform thickness is assumed. Two-way coupled fluid structure interaction have been employed for the numerical simulation of blood flow dynamics and arterial wall mechanics.

RESULTS

Wall shear stress (WSS) parameters and vonmises stress indicators, co-relating rupture and thrombus formation, have been extracted and reported, at each growth stage. When the aneurysm progresses in diameter, the areas recording abnormally low TAWSS, as well as areas of high/low OSI were found to increase at different rates for shorter and longer aneurysms. Moreover, drastic increase in the maximum wall stresses (MWS) and wall displacement were observed as the aneurysm approached the critical diameter.

CONCLUSION

Hemodynamic predictors were found to be highly dependent on the shape index (DHr), when the aneurysm was small, whereas significant influence of DHr on the wall stresses happens, as the aneurysm approaches the critical diameter. Inconsistent variation of these indicators exhibited by shorter aneurysms (high DHr) at different growth stages, demands routine monitoring (using scans), of such aneurysms, to prevent unexpected rupture.

Numerical Investigation of the Fetal Left Heart Hemodynamics During Gestational Stages

Flow-driven hemodynamic forces on the cardiac tissues have critical importance, and have a significant role in the proper development of the heart. These mechanobiological mechanisms govern the cellular responses for the growth and remodeling of the heart, where the altered hemodynamic environment is believed to be a major factor that is leading to congenital heart defects (CHDs). In order to investigate the mechanobiological development of the normal and diseased hearts, identification of the blood flow patterns and wall shear stresses (WSS) on these tissues are required for an accurate hemodynamic assessment. In this study, we focus on the left heart hemodynamics of the human fetuses throughout the gestational stages. Computational fetal left heart models are created for the healthy fetuses using the ultrasound images at various gestational weeks. Realistic inflow boundary conditions are implemented in the models using the Doppler ultrasound measurements for resolving the specific blood flow waveforms in the mitral valve. Obtained results indicate that WSS and vorticity levels in the fetal left heart decrease with the development of the fetus. The maximum WSS around the mitral valve is determined around 36 Pa at the gestational week of 16. This maximum WSS decreases to 11 Pa at the gestational week of 27, indicating nearly three-times reduction in the peak shear stress. These findings reveal the highly dynamic nature of the left heart hemodynamics throughout the development of the human fetus and shed light into the relevance of hemodynamic environment and development of CHDs.

Ultrasound-Based Fluid-Structure Interaction Modeling of Abdominal Aortic Aneurysms Incorporating Pre-stress

Currently, the prediction of rupture risk in abdominal aortic aneurysms (AAAs) solely relies on maximum diameter. However, wall mechanics and hemodynamics have shown to provide better risk indicators. Patient-specific fluid-structure interaction (FSI) simulations based on a non-invasive image modality are required to establish a patient-specific risk indicator. In this study, a robust framework to execute FSI simulations based on time-resolved three-dimensional ultrasound (3D+t US) data was obtained and employed on a data set of 30 AAA patients. Furthermore, the effect of including a pre-stress estimation (PSE) to obtain the stresses present in the measured geometry was evaluated. The established workflow uses the patient-specific 3D+t US-based segmentation and brachial blood pressure as input to generate meshes and boundary conditions for the FSI simulations. The 3D+t US-based FSI framework was successfully employed on an extensive set of AAA patient data. Omitting the pre-stress results in increased displacements, decreased wall stresses, and deviating time-averaged wall shear stress and oscillatory shear index patterns. These results underline the importance of incorporating pre-stress in FSI simulations. After validation, the presented framework provides an important tool for personalized modeling and longitudinal studies on AAA growth and rupture risk.

Boundary conditions investigation to improve computer simulation of cerebrospinal fluid dynamics in hydrocephalus patients

Three-D head geometrical models of eight healthy subjects and 11 hydrocephalus patients were built using their CINE phase-contrast MRI data and used for computer simulations under three different inlet/outlet boundary conditions (BCs). The maximum cerebrospinal fluid (CSF) pressure and the ventricular system volume were more effective and accurate than the other parameters in evaluating the patients' conditions. In constant CSF pressure, the computational patient models were 18.5% more sensitive to CSF volume changes in the ventricular system under BC "C". Pulsatile CSF flow rate diagrams were used for inlet and outlet BCs of BC "C". BC "C" was suggested to evaluate the intracranial compliance of the hydrocephalus patients. The results suggested using the computational fluid dynamic (CFD) method and the fully coupled fluid-structure interaction (FSI) method for the CSF dynamic analysis in patients with external and internal hydrocephalus, respectively.

Computational Modeling of Blood Flow Hemodynamics for Biomechanical Investigation of Cardiac Development and Disease

The heart is the first functional organ in a developing embryo. Cardiac development continues throughout developmental stages while the heart goes through a serious of drastic morphological changes. Previous animal experiments as well as clinical observations showed that disturbed hemodynamics interfere with the development of the heart and leads to the formation of a variety of defects in heart valves, heart chambers, and blood vessels, suggesting that hemodynamics is a governing factor for cardiogenesis, and disturbed hemodynamics is an important source of congenital heart defects. Therefore, there is an interest to image and quantify the flowing blood through a developing heart. Flow measurement in embryonic fetal heart can be performed using advanced techniques such as magnetic resonance imaging (MRI) or echocardiography. Computational fluid dynamics (CFD) modeling is another approach especially useful when the other imaging modalities are not available and in-depth flow assessment is needed. The approach is based on numerically solving relevant physical equations to approximate the flow hemodynamics and tissue behavior. This approach is becoming widely adapted to simulate cardiac flows during the embryonic development. While there are few studies for human fetal cardiac flows, many groups used zebrafish and chicken embryos as useful models for elucidating normal and diseased cardiogenesis. In this paper, we explain the major steps to generate CFD models for simulating cardiac hemodynamics in vivo and summarize the latest findings on chicken and zebrafish embryos as well as human fetal hearts.

On the importance of tunica intima in the aging aorta: a three-layered in silico model for computing wall stresses in abdominal aortic aneurysms

Layer-specific experimental data for human aortic tissue suggest that, in aged arteries and arteries with non-atherosclerotic intimal thickening, the innermost layer of the aorta increases significantly its stiffness and thickness, becoming load-bearing. However, there are very few computational studies of abdominal aortic aneurysms (AAAs) that take into account the mechanical contribution of the three layers that comprise the aneurysmal tissue. In this paper, a three-layered finite element model is proposed from the simplest uniaxial stress state to geometrically parametrized models of AAAs with different asymmetry values. Comparisons are made between a three-layered artery wall and a mono-layered intact artery, which represents the complex behavior of the aggregate adventitia-media-intima in a single layer with averaged mechanical properties. Likewise, the response of our idealized geometries is compared with similar experimental and numerical models. Finally, the mechanical contributions of adventitia, media and intima are analyzed for the three-layered aneurysms through the evaluation of the mean stress absorption percentage. Results show the relevance and necessity of considering the inclusion of tunica intima in multi-layered models of AAAs for getting accurate results in terms of peak wall stresses and displacements.

Targeted Drug Delivery of Microbubble to Arrest Abdominal Aortic Aneurysm Development: A Simulation Study Towards Optimized Microbubble Design

Abdominal aortic aneurysm (AAA) is an irreversible bulge in the artery with higher prevalence among the elderlies. Increase of the aneurysm diameter by time is a fatal phenomenon which will lead to its sidewall rupture. Invasive surgical treatments are vital in preventing from AAA development. These approaches however have considerable side effects. Targeted drug delivery using microbubbles (MBs) has been recently employed to suppress the AAA growth. The present study is aimed to investigate the surface adhesion of different types of drug-containing MBs to the inner wall of AAA through ligand-receptor binding, using fluid-structure interaction (FSI) simulation by using a patient CT-scan images of the vascular system. The effect of blood flow through AAA on MBs delivery to the intended surface was also addressed. For this purpose, the adherence of four types of MBs with three different diameters to the inner surface wall of AAA was studied in a patient with 40-mm diameter aneurysm. The effects of the blood mechanical properties on the hematocrit (Hct) percentage of patients suffering from anemia and diabetes were studied. Moreover, the impact of variations in the artery inlet velocity on blood flow was addressed. Simulation results demonstrated the dependency of the surface density of MBs (SDM) adhered on the AAA lumen to the size and the type of MBs. It was observed that the amount of SDM due to adhesion on the AAA lumen for one of the commercially-approved MBs (Optison) with a diameter of 4.5 μm was higher than the other MBs. Furthermore, we have shown that the targeted drug delivery to the AAA lumen is more favorable in healthy individuals (45% Hct) compared to the patients with diabetes and anemia. Also, it was found that the targeted drug delivery method using MBs on the patients having AAA with complicated aneurysm shape and negative inlet blood flow velocity can be severely affected.

Numerical investigation of the effects of blood rheology and wall elasticity in abdominal aortic aneurysm under pulsatile flow conditions

BACKGROUND

Previous studies on aneurysm modeling have focused on the blood rheology and vessel elasticity separately. The combined effects of blood shear thinning properties and wall elasticity need to be revealed.

OBJECTIVE

To provide insights on how pulsatile hemodynamics vary with blood rheology and vessel elasticity for a developed abdominal aortic aneurysm (AAA).

METHOD

An Arbitrary Lagrangian-Eulerian fluid-solid interaction method is adopted with the Newtonian and the shear thinning Carreau constitutive models for the fluid with the linearly elastic and the hyperelastic Yeoh models for the vessel. Finite element based numerical solver is used to simulate the blood flow in the AAA.

RESULTS

Newtonian model overestimates the velocity values compared to the Carreau model and the difference in the velocity field increases as the shear rate decreases at the instances of the cardiac cycle. The rigid walled simulations display higher deviations in the velocity and wall shear stress with the fluid rheology. The risk indicators show that Newtonian assumption combined with the linearly elastic model may overlook degeneration risk of arterial tissue.

CONCLUSIONS

Newtonian assumption for the blood as well as modelling the arterial wall as linearly elastic lead to significant differences in oscillatory hemodynamic properties with respect to the use of Carreau fluid together with hyperelastic vessel model, even in large vessel aneurysms.

Fluid–Solid Interaction Analysis on Iliac Bifurcation Artery: A Numerical Study

Atherosclerosis, which is commonly seen at regions with low wall shear stress (WSS) level in bifurcations, is a kind of fibro-fatty plaque accumulated on arterial walls. Aortic and iliac bifurcations have the highest proportion of patients among all atherosclerosis cases, thus it is necessary to numerically analyze the flow distribution and predict plaque positions in these bifurcations. Furthermore, using fluid–solid interaction (FSI) method could obtain a more exact flow pattern in arteries. In this study, a patient-specific model of aortic and iliac bifurcations was simulated with both FSI and rigid-wall cases. We analyzed the vessel deformation, WSS and flow distribution of this model. Computed tomography (CT) angiography was used in our study to create patient-specific model of aorto-iliac arteries. Real material properties and pulsatile fluid boundary conditions were applied in solid and fluid zones, respectively. We performed FSI and ordinary computational fluid dynamics (CFD) simulations with AYSYS 15.0 software (ANSYS Inc., Canonsburg, PA), and compared the diameter change, WSS and flow field between these two results. The diameter change between systolic phase and diastolic phase is 8–9% on abdominal aorta, and 3% on external and internal iliac arteries. The compliance of vessels corresponds to in-vivo observations. At peak systolic phase, the average WSS obtained in FSI simulation is 10% lower than in rigid-wall result, area of low-WSS region ([Formula: see text]) also increases by 78%. Wall deformation has a greater impact on WSS of those vessels with larger diameter, but hardly changes the shear level in smaller branches. Our result also shows that iliac bifurcations reveal more complicated secondary flow in systolic phase, comparing to other vessels, and stenosed iliac artery has more severe secondary flow than healthy artery. We obtained a feasible method for hemodynamic FSI research. The material parameters, boundary conditions and mesh could be used for further simulations, while the WSS and flow distribution may support clinical diagnosis and treatment. We concluded that compliance is a must-consider factor for simulating an accurate wall shear stress, because the vessel deformation in FSI simulation will significantly change the distribution of low-WSS zones. Moreover, more complicated secondary flow is detected in iliac arteries because it may interact between bifurcations. Stenosis in artery may also have a blocking effect on downstream blood flow.

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.

Fluid-Structure Interaction in Abdominal Aortic Aneurysms: Effect of Haematocrit

The Abdominal Aortic Aneurysm (AAA) is a local dilation of the abdominal aorta and it is a cause for serious concern because of the high mortality associated with its rupture. Consequently, the understanding of the phenomena related to the creation and the progression of an AAA is of crucial importance. In this work, the complicated interaction between the blood flow and the AAA wall is numerically examined using a fully coupled Fluid-Structure Interaction (FSI) method. The study investigates the possible link between the dynamic behavior of an AAA and the blood viscosity variations attributed to the haematocrit value, while it also incorporates the pulsatile blood flow, the non-Newtonian behavior of blood and the hyperelasticity of the arterial wall. It was found that blood viscosity has no significant effect on von Mises stress magnitude and distribution, whereas there is a close relation between the haematocrit value and the Wall Shear Stress (WSS) magnitude in AAAs. This WSS variation can possibly alter the mechanical properties of the arterial wall and increase its growth rate or even its rupture possibility. The relationship between haematocrit and dynamic behavior of an AAA can be helpful in designing a patient specific treatment.

Regional pulse wave velocity and stress in aneurysmal arch-shaped aorta

The pulse wave velocity (PWV) has been shown to be associated with the properties of blood vessel and a cardiovascular risk factor such as aneurysm. The global PWV estimation is applied in conventional clinical diagnosis. However, the geometry of blood vessel changes along the wave traveling path and the global PWV estimation may not always detect regional wall changes resulting from cardiovascular diseases. In this study, a fluid structure interaction (FSI) analysis was applied on arch-shaped aortas with and without aneurysm aimed at determining the effects of the number of aneurysm, aneurysm size and the modulus ratio (aneurysm to wall modulus) on the pulse wave propagation and velocity. The characterization for each stage of aneurysmal aorta was simulated by progressively increasing aortic stiffness and aneurysm size. The pulse wave propagations and velocities were estimated from the two-dimensional spatial-temporal plot of the normalized wall displacement based on elastic deformation. The descending forward and arch reflected PWVs of aneurysmal aortic arch models were found up to 9.7% and 122.8%, respectively, deviate from the PWV of non-aneurysmal aortic arch model. The PWV patterns and magnitudes can be used to distinguish the characterization of the normal and aneurysmal aortic walls and shown to be relevant regional markers utilized in clinical diagnosis.

Fluid-structure interaction modeling of aneurysmal arteries under steady-state and pulsatile blood flow: a stability analysis

Tortuous aneurysmal arteries are often associated with a higher risk of rupture but the mechanism remains unclear. The goal of this study was to analyze the buckling and post-buckling behaviors of aneurysmal arteries under pulsatile flow. To accomplish this goal, we analyzed the buckling behavior of model carotid and abdominal aorta with aneurysms by utilizing fluid-structure interaction (FSI) method with realistic waveforms boundary conditions. FSI simulations were done under steady-state and pulsatile flow for normal (1.5) and reduced (1.3) axial stretch ratios to investigate the influence of aneurysm, pulsatile lumen pressure and axial tension on stability. Our results indicated that aneurysmal artery buckled at the critical buckling pressure and its deflection nonlinearly increased with increasing lumen pressure. Buckling elevates the peak stress (up to 118%). The maximum aneurysm wall stress at pulsatile FSI flow was (29%) higher than under static pressure at the peak lumen pressure of 130 mmHg. Buckling results show an increase in lumen shear stress at the inner side of the maximum deflection. Vortex flow was dramatically enlarged with increasing lumen pressure and artery diameter. Aneurysmal arteries are more susceptible than normal arteries to mechanical instability which causes high stresses in the aneurysm wall that could lead to aneurysm rupture.

The influence of fiber dispersion on the mechanical response of aortic tissues in health and disease: a computational study

Changes in the structural components of aortic tissues have been shown to play a significant role in the pathogenesis of aortic degeneration. Therefore, reliable stress analyses require a suitable and meaningful constitutive model that captures micro-structural changes. As recent data show, in-plane and out-of-plane collagen fiber dispersions vary significantly between healthy and aneurysmatic aortic walls. The aim of this study is to computationally investigate the influence of fiber dispersion on the mechanical response of aortic tissues in health and disease. In particular, the influence of three different fiber dispersions is studied: (i) non-rotationally symmetric dispersion, the most realistic assumption for aortic tissues; (ii) transversely isotropic dispersion, a special case; (iii) perfectly aligned fibers (no dispersion in either plane), another special case. Explicit expressions for the stress and elasticity tensors as needed for the implementation in a finite element code are provided. Three representative numerical examples are studied: planar biaxial extension, inflation of residually stressed and pre-stretched aortic segments and inflation of an idealized abdominal aortic aneurysm (AAA) geometry. For the AAA geometry the case of isotropic dispersion is additionally analyzed. Documented structural and mechanical parameters are taken from human aortas (healthy media/adventitia and AAA). The influence of fiber dispersions upon magnitudes and distributions of stresses and deformations are presented and analyzed. Stresses vary significantly, especially in the AAA case, where material stiffening is significantly influenced by fiber dispersion. The results highlight the need to incorporate the structural differences into finite element simulations to obtain more accurate stress predictions. Additionally, results show the capability of one constitutive model to represent different scenarios of aortic micro-structures allowing future studies of collagen reorientation during disease progression.

INFLUENCE OF THE LOCAL MEAN CURVATURE ON THE ABDOMINAL AORTIC ANEURYSM STRESS DISTRIBUTION

An abdominal aortic aneurysm (AAA) is a permanent focal dilatation of the abdominal aorta of at least 1.5 times its normal diameter. Although the criterion of maximum diameter is still used in clinical practice to decide when to proceed with surgical intervention, numerical studies have demonstrated the importance of other geometric factors. In this work, the influence of the local mean curvature on AAA stress distribution has been analyzed in synthetic AAA geometries via finite element analysis. The results show a significant correlation between this geometric parameter and stress, suggesting that local mean curvature should also be considered along with the diameter criterion when making decisions about surgery.

The Obsolete Maximum Diameter Criterion, the Evident Role of Biomechanical (Pressure) Indices, the New Role of Hemodynamic (Flow) Indices, and the Multi-Modal Approach to the Rupture Risk Assessment of Abdominal Aortic Aneurysms

Although the therapeutic management of abdominal aortic aneurysms (AAAs) is currently based on the maximum diameter criterion, this has often proved inaccurate and misleading. Conversely, the biomechanical approach, which takes into account the pressure-induced wall stress exerted at every point throughout the aneurysmal surface, has been proven superior in predicting the rupture risk of AAAs, and its value is being increasingly recognized among physicians. More recently, hemodynamic indices, such as flow-induced wall shear stresses, have been indicated as potentially significant determinants of AAA natural history. Ultimately, a statistical model that takes into account all these factors may be relevant for making a sound prediction of the rupture risk of aneurysms and optimizing the management of these patients.

Deformation and dynamic response of abdominal aortic aneurysm sealing

Endovascular sealing is a new technique for the repair of abdominal aortic aneurysms. Commercially available in Europe since 2013, it takes a revolutionary approach to aneurysm repair through minimally invasive techniques. Although aneurysm sealing may be thought as more stable than conventional endovascular stent graft repairs, post-implantation movement of the endoprosthesis has been described, potentially leading to late complications. The paper presents for the first time a model, which explains the nature of forces, in static and dynamic regimes, acting on sealed abdominal aortic aneurysms, with references to real case studies. It is shown that elastic deformation of the aorta and of the endoprosthesis induced by static forces and vibrations during daily activities can potentially promote undesired movements of the endovascular sealing structure.

Deformation and dynamic response of abdominal aortic aneurysm sealing

Endovascular sealing is a new technique for the repair of abdominal aortic aneurysms. Commercially available in Europe since 2013, it takes a revolutionary approach to aneurysm repair through minimally invasive techniques. Although aneurysm sealing may be thought as more stable than conventional endovascular stent graft repairs, post-implantation movement of the endoprosthesis has been described, potentially leading to late complications. The paper presents for the first time a model, which explains the nature of forces, in static and dynamic regimes, acting on sealed abdominal aortic aneurysms, with references to real case studies. It is shown that elastic deformation of the aorta and of the endoprosthesis induced by static forces and vibrations during daily activities can potentially promote undesired movements of the endovascular sealing structure.

A Methodology for Verifying Abdominal Aortic Aneurysm Wall Stress

An abdominal aortic aneurysm (AAA) is a permanent focal dilatation of the abdominal aorta of at least 1.5 times its normal diameter. Although the criterion of maximum diameter is still used in clinical practice to decide on a timely intervention, numerical studies have demonstrated the importance of other geometric factors. However, the major drawback of numerical studies is that they must be validated experimentally before clinical implementation. This work presents a new methodology to verify wall stress predicted from the numerical studies against the experimental testing. To this end, four AAA phantoms were manufactured using vacuum casting. The geometry of each phantom was subject to microcomputed tomography (μCT) scanning at zero and three other intraluminal pressures: 80, 100, and 120 mm Hg. A zero-pressure geometry algorithm was used to calculate the wall stress in the phantom, while the numerical wall stress was calculated with a finite-element analysis (FEA) solver based on the actual zero-pressure geometry subjected to 80, 100, and 120 mm Hg intraluminal pressure loading. Results demonstrate the moderate accuracy of this methodology with small relative differences in the average wall stress (1.14%). Additionally, the contribution of geometric factors to the wall stress distribution was statistically analyzed for the four phantoms. The results showed a significant correlation between wall thickness and mean curvature (MC) with wall stress.

A Methodology for the Derivation of Unloaded Abdominal Aortic Aneurysm Geometry With Experimental Validation

In this work, we present a novel method for the derivation of the unloaded geometry of an abdominal aortic aneurysm (AAA) from a pressurized geometry in turn obtained by 3D reconstruction of computed tomography (CT) images. The approach was experimentally validated with an aneurysm phantom loaded with gauge pressures of 80, 120, and 140 mm Hg. The unloaded phantom geometries estimated from these pressurized states were compared to the actual unloaded phantom geometry, resulting in mean nodal surface distances of up to 3.9% of the maximum aneurysm diameter. An in-silico verification was also performed using a patient-specific AAA mesh, resulting in maximum nodal surface distances of 8 μm after running the algorithm for eight iterations. The methodology was then applied to 12 patient-specific AAA for which their corresponding unloaded geometries were generated in 5-8 iterations. The wall mechanics resulting from finite element analysis of the pressurized (CT image-based) and unloaded geometries were compared to quantify the relative importance of using an unloaded geometry for AAA biomechanics. The pressurized AAA models underestimate peak wall stress (quantified by the first principal stress component) on average by 15% compared to the unloaded AAA models. The validation and application of the method, readily compatible with any finite element solver, underscores the importance of generating the unloaded AAA volume mesh prior to using wall stress as a biomechanical marker for rupture risk assessment.

Endoleak Assessment Using Computational Fluid Dynamics and Image Processing Methods in Stented Abdominal Aortic Aneurysm Models

Endovascular aortic aneurysm repair (EVAR) is a predominant surgical procedure to reduce the risk of aneurysm rupture in abdominal aortic aneurysm (AAA) patients. Endoleak formation, which eventually requires additional surgical reoperation, is a major EVAR complication. Understanding the etiology and evolution of endoleak from the hemodynamic perspective is crucial to advancing the current posttreatments for AAA patients who underwent EVAR. Therefore, a comprehensive flow assessment was performed to investigate the relationship between endoleak and its surrounding pathological flow fields through computational fluid dynamics and image processing. Six patient-specific models were reconstructed, and the associated hemodynamics in these models was quantified three-dimensionally to calculate wall stress. To provide a high degree of clinical relevance, the mechanical stress distribution calculated from the models was compared with the endoleak positions identified from the computed tomography images of patients through a series of imaging processing methods. An endoleak possibly forms in a location with high local wall stress. An improved stent graft (SG) structure is conceived accordingly by increasing the mechanical strength of the SG at peak wall stress locations. The presented analytical paradigm, as well as numerical analysis using patient-specific models, may be extended to other common human cardiovascular surgeries.

FSI Simulations of Pulse Wave Propagation in Human Abdominal Aortic Aneurysm: The Effects of Sac Geometry and Stiffness

This study aims to quantify the effects of geometry and stiffness of aneurysms on the pulse wave velocity (PWV) and propagation in fluid–solid interaction (FSI) simulations of arterial pulsatile flow. Spatiotemporal maps of both the wall displacement and fluid velocity were generated in order to obtain the pulse wave propagation through fluid and solid media, and to examine the interactions between the two waves. The results indicate that the presence of abdominal aortic aneurysm (AAA) sac and variations in the sac modulus affect the propagation of the pulse waves both qualitatively (eg, patterns of change of forward and reflective waves) and quantitatively (eg, decreasing of PWV within the sac and its increase beyond the sac as the sac stiffness increases). The sac region is particularly identified on the spatiotemporal maps with a region of disruption in the wave propagation with multiple short-traveling forward/reflected waves, which is caused by the change in boundary conditions within the saccular region. The change in sac stiffness, however, is more pronounced on the wall displacement spatiotemporal maps compared to those of fluid velocity. We conclude that the existence of the sac can be identified based on the solid and fluid pulse waves, while the sac properties can also be estimated. This study demonstrates the initial findings in numerical simulations of FSI dynamics during arterial pulsations that can be used as reference for experimental and in vivo studies. Future studies are needed to demonstrate the feasibility of the method in identifying very mild sacs, which cannot be detected from medical imaging, where the material property degradation exists under early disease initiation.

Compliance mismatch and compressive wall stresses drive anomalous remodelling of pulmonary trunks reinforced with Dacron grafts

Synthetic grafts are often satisfactory employed in cardiac and vascular surgery, including expanded poly(ethylene terephthalate) or expanded poly(tetrafluoroethylene). However, accumulating evidences suggest the emergence of worrisome issues concerning the long-term fate of prosthetic grafts as large vessel replacement. Disadvantages related to the use of synthetic grafts can be traced in their inability of mimicking the elasto-mechanical characteristics of the native vascular tissue, local suture overstress leading to several prosthesis-related complications and retrograde deleterious effects on valve competence, cardiac function and perfusion. Motivated by this, in the present work it is analyzed - by means of both elemental biomechanical paradigms and more accurate in silico Finite Element simulations - the physical interaction among aorta, autograft and widely adopted synthetic (Dacron) prostheses utilized in transposition of pulmonary artery, highlighting the crucial role played by somehow unexpected stress fields kindled in the vessel walls and around suture regions, which could be traced as prodromal to the triggering of anomalous remodelling processes and alterations of needed surgical outcomes. Theoretical results are finally compared with histological and surgical data related to a significant experimental animal campaign conducted by performing pulmonary artery transpositions in 30 two-month old growing lambs, followed up during growth for six months. The in vivo observations demonstrate the effectiveness of the proposed biomechanical hypothesis and open the way for possible engineering-guided strategies to support and optimize surgical procedures.

Patient-specific modelling of abdominal aortic aneurysms: The influence of wall thickness on predicted clinical outcomes

Rupture of abdominal aortic aneurysms (AAAs) is linked to aneurysm morphology. This study investigates the influence of patient-specific (PS) AAA wall thickness on predicted clinical outcomes. Eight patients under surveillance for AAAs were selected from the MA(3)RS clinical trial based on the complete absence of intraluminal thrombus. Two finite element (FE) models per patient were constructed; the first incorporated variable wall thickness from CT (PS_wall), and the second employed a 1.9mm uniform wall (Uni_wall). Mean PS wall thickness across all patients was 1.77±0.42mm. Peak wall stress (PWS) for PS_wall and Uni_wall models was 0.6761±0.3406N/mm(2) and 0.4905±0.0850N/mm(2), respectively. In 4 out of 8 patients the Uni_wall underestimated stress by as much as 55%; in the remaining cases it overestimated stress by up to 40%. Rupture risk more than doubled in 3 out of 8 patients when PS_wall was considered. Wall thickness influenced the location and magnitude of PWS as well as its correlation with curvature. Furthermore, the volume of the AAA under elevated stress increased significantly in AAAs with higher rupture risk indices. This highlights the sensitivity of standard rupture risk markers to the specific wall thickness strategy employed.

A Review of Computational Methods to Predict the Risk of Rupture of Abdominal Aortic Aneurysms

Computational methods have played an important role in health care in recent years, as determining parameters that affect a certain medical condition is not possible in experimental conditions in many cases. Computational fluid dynamics (CFD) methods have been used to accurately determine the nature of blood flow in the cardiovascular and nervous systems and air flow in the respiratory system, thereby giving the surgeon a diagnostic tool to plan treatment accordingly. Machine learning or data mining (MLD) methods are currently used to develop models that learn from retrospective data to make a prediction regarding factors affecting the progression of a disease. These models have also been successful in incorporating factors such as patient history and occupation. MLD models can be used as a predictive tool to determine rupture potential in patients with abdominal aortic aneurysms (AAA) along with CFD-based prediction of parameters like wall shear stress and pressure distributions. A combination of these computer methods can be pivotal in bridging the gap between translational and outcomes research in medicine. This paper reviews the use of computational methods in the diagnosis and treatment of AAA.

Computational Growth and Remodeling of Abdominal Aortic Aneurysms Constrained by the Spine

Abdominal aortic aneurysms (AAAs) evolve over time, and the vertebral column, which acts as an external barrier, affects their biomechanical properties. Mechanical interaction between AAAs and the spine is believed to alter the geometry, wall stress distribution, and blood flow, although the degree of this interaction may depend on AAAs specific configurations. In this study, we use a growth and remodeling (G&R) model, which is able to trace alterations of the geometry, thus allowing us to computationally investigate the effect of the spine for progression of the AAA. Medical image-based geometry of an aorta is constructed along with the spine surface, which is incorporated into the computational model as a cloud of points. The G&R simulation is initiated by local elastin degradation with different spatial distributions. The AAA-spine interaction is accounted for using a penalty method when the AAA surface meets the spine surface. The simulation results show that, while the radial growth of the AAA wall is prevented on the posterior side due to the spine acting as a constraint, the AAA expands faster on the anterior side, leading to higher curvature and asymmetry in the AAA configuration compared to the simulation excluding the spine. Accordingly, the AAA wall stress increases on the lateral, posterolateral, and the shoulder regions of the anterior side due to the AAA-spine contact. In addition, more collagen is deposited on the regions with a maximum diameter. We show that an image-based computational G&R model not only enhances the prediction of the geometry, wall stress, and strength distributions of AAAs but also provides a framework to account for the interactions between an enlarging AAA and the spine for a better rupture potential assessment and management of AAA patients.

Patient-specific computational modeling of blood flow in the pulmonary arterial circulation

Computational fluid dynamics (CFD) modeling of the pulmonary vasculature has the potential to reveal continuum metrics associated with the hemodynamic stress acting on the vascular endothelium. It is widely accepted that the endothelium responds to flow-induced stress by releasing vasoactive substances that can dilate and constrict blood vessels locally. The objectives of this study are to examine the extent of patient specificity required to obtain a significant association of CFD output metrics and clinical measures in models of the pulmonary arterial circulation, and to evaluate the potential correlation of wall shear stress (WSS) with established metrics indicative of right ventricular (RV) afterload in pulmonary hypertension (PH). Right Heart Catheterization (RHC) hemodynamic data and contrast-enhanced computed tomography (CT) imaging were retrospectively acquired for 10 PH patients and processed to simulate blood flow in the pulmonary arteries. While conducting CFD modeling of the reconstructed patient-specific vasculatures, we experimented with three different outflow boundary conditions to investigate the potential for using computationally derived spatially averaged wall shear stress (SAWSS) as a metric of RV afterload. SAWSS was correlated with both pulmonary vascular resistance (PVR) (R(2)=0.77, P<0.05) and arterial compliance (C) (R(2)=0.63, P<0.05), but the extent of the correlation was affected by the degree of patient specificity incorporated in the fluid flow boundary conditions. We found that decreasing the distal PVR alters the flow distribution and changes the local velocity profile in the distal vessels, thereby increasing the local WSS. Nevertheless, implementing generic outflow boundary conditions still resulted in statistically significant SAWSS correlations with respect to both metrics of RV afterload, suggesting that the CFD model could be executed without the need for complex outflow boundary conditions that require invasively obtained patient-specific data. A preliminary study investigating the relationship between outlet diameter and flow distribution in the pulmonary tree offers a potential computationally inexpensive alternative to pressure based outflow boundary conditions.

The role of wall shear stress in the assessment of right ventricle hydraulic workload

Pulmonary hypertension (PH) is a devastating disease affecting approximately 15-50 people per million, with a higher incidence in women. PH mortality is mostly attributed to right ventricle (RV) failure, which results from RV hypotrophy due to an overburdened hydraulic workload. The objective of this study is to correlate wall shear stress (WSS) with hemodynamic metrics that are generally accepted as clinical indicators of RV workload and are well correlated with disease outcome. Retrospective right heart catheterization data for 20 PH patients were analyzed to derive pulmonary vascular resistance (PVR), arterial compliance (C), and an index of wave reflections (Γ). Patient-specific contrast-enhanced computed tomography chest images were used to reconstruct the individual pulmonary arterial trees up to the seventh generation. Computational fluid dynamics analyses simulating blood flow at peak systole were conducted for each vascular model to calculate WSS distributions on the endothelial surface of the pulmonary arteries. WSS was found to be decreased proportionally with elevated PVR and reduced C. Spatially averaged WSS (SAWSS) was positively correlated with PVR (R (2) = 0.66), C (R (2) = 0.73), and Γ (R (2) = 0.5) and also showed promising preliminary correlations with RV geometric characteristics. Evaluating WSS at random cross sections in the proximal vasculature (main, right, and left pulmonary arteries), the type of data that can be acquired from phase-contrast magnetic resonance imaging, did not reveal the same correlations. In conclusion, we found that WSS has the potential to be a viable and clinically useful noninvasive metric of PH disease progression and RV health. Future work should be focused on evaluating whether SAWSS has prognostic value in the management of PH and whether it can be used as a rapid reactivity assessment tool, which would aid in selection of appropriate therapies.

An animal-specific FSI model of the abdominal aorta in anesthetized mice

Recent research has revealed that angiotensin II-induced abdominal aortic aneurysm in mice can be related to medial ruptures occurring in the vicinity of abdominal side branches. Nevertheless a thorough understanding of the biomechanics near abdominal side branches in mice is lacking. In the current work we present a mouse-specific fluid-structure interaction (FSI) model of the abdominal aorta in ApoE(-/-) mice that incorporates in vivo stresses. The aortic geometry was based on contrast-enhanced in vivo micro-CT images, while aortic flow boundary conditions and material model parameters were based on in vivo high-frequency ultrasound. Flow waveforms predicted by FSI simulations corresponded better to in vivo measurements than those from CFD simulations. Peak-systolic principal stresses at the inner and outer aortic wall were locally increased caudal to the celiac and left lateral to the celiac and mesenteric arteries. Interestingly, these were also the locations at which a tear in the tunica media had been observed in previous work on angiotensin II-infused mice. Our preliminary results therefore suggest that local biomechanics play an important role in the pathophysiology of branch-related ruptures in angiotensin-II infused mice. More elaborate follow-up research is needed to demonstrate the role of biomechanics and mechanobiology in a longitudinal setting.

Investigation of material modeling in fluid-structure interaction analysis of an idealized three-layered abdominal aorta: aneurysm initiation and fully developed aneurysms

Different material models for an idealized three-layered abdominal aorta are compared using computational techniques to study aneurysm initiation and fully developed aneurysms. The computational model includes fluid-structure interaction (FSI) between the blood vessel and the blood. In order to model aneurysm initiation, the medial region was degenerated to mimic the medial loss occurring in the inception of an aneurysm. Various cases are considered in order to understand their effects on the initiation of an abdominal aortic aneurysm. The layers of the blood vessel were modeled using either linear elastic materials or Mooney-Rivlin (otherwise known as hyperelastic) type materials. The degenerated medial region was also modeled in either linear elastic or hyperelastic-type materials and assumed to be in the shape of an arc with a thin width or a circular ring with different widths. The blood viscosity effect was also considered in the initiation mechanism. In addition, dynamic analysis of the blood vessel was performed without interaction with the blood flow by applying time-dependent pressure inside the lumen in a three-layered abdominal aorta. The stresses, strains, and displacements were compared for a healthy aorta, an initiated aneurysm and a fully developed aneurysm. The study shows that the material modeling of the vessel has a sizable effect on aneurysm initiation and fully developed aneurysms. Different material modeling of degeneration regions also affects the stress-strain response of aneurysm initiation. Additionally, the structural analysis without considering FSI (called noFSI) overestimates the peak von Mises stress by 52% at the interfaces of the layers.

Effects of intraluminal thrombus on patient-specific abdominal aortic aneurysm hemodynamics via stereoscopic particle image velocity and computational fluid dynamics modeling

The pathology of the human abdominal aortic aneurysm (AAA) and its relationship to the later complication of intraluminal thrombus (ILT) formation remains unclear. The hemodynamics in the diseased abdominal aorta are hypothesized to be a key contributor to the formation and growth of ILT. The objective of this investigation is to establish a reliable 3D flow visualization method with corresponding validation tests with high confidence in order to provide insight into the basic hemodynamic features for a better understanding of hemodynamics in AAA pathology and seek potential treatment for AAA diseases. A stereoscopic particle image velocity (PIV) experiment was conducted using transparent patient-specific experimental AAA models (with and without ILT) at three axial planes. Results show that before ILT formation, a 3D vortex was generated in the AAA phantom. This geometry-related vortex was not observed after the formation of ILT, indicating its possible role in the subsequent appearance of ILT in this patient. It may indicate that a longer residence time of recirculated blood flow in the aortic lumen due to this vortex caused sufficient shear-induced platelet activation to develop ILT and maintain uniform flow conditions. Additionally, two computational fluid dynamics (CFD) modeling codes (Fluent and an in-house cardiovascular CFD code) were compared with the two-dimensional, three-component velocity stereoscopic PIV data. Results showed that correlation coefficients of the out-of-plane velocity data between PIV and both CFD methods are greater than 0.85, demonstrating good quantitative agreement. The stereoscopic PIV study can be utilized as test case templates for ongoing efforts in cardiovascular CFD solver development. Likewise, it is envisaged that the patient-specific data may provide a benchmark for further studying hemodynamics of actual AAA, ILT, and their convolution effects under physiological conditions for clinical applications.

The influence of computational assumptions on analysing abdominal aortic aneurysm haemodynamics

The variation in computational assumptions for analysing abdominal aortic aneurysm haemodynamics can influence the desired output results and computational cost. Such assumptions for abdominal aortic aneurysm modelling include static/transient pressures, steady/transient flows and rigid/compliant walls. Six computational methods and these various assumptions were simulated and compared within a realistic abdominal aortic aneurysm model with and without intraluminal thrombus. A full transient fluid-structure interaction was required to analyse the flow patterns within the compliant abdominal aortic aneurysms models. Rigid wall computational fluid dynamics overestimates the velocity magnitude by as much as 40%-65% and the wall shear stress by 30%-50%. These differences were attributed to the deforming walls which reduced the outlet volumetric flow rate for the transient fluid-structure interaction during the majority of the systolic phase. Static finite element analysis accurately approximates the deformations and von Mises stresses when compared with transient fluid-structure interaction. Simplifying the modelling complexity reduces the computational cost significantly. In conclusion, the deformation and von Mises stress can be approximately found by static finite element analysis, while for compliant models a full transient fluid-structure interaction analysis is required for acquiring the fluid flow phenomenon.