Fluid-structure interaction within realistic three-dimensional models of the aneurysmatic aorta as a guidance to assess the risk of rupture of the aneurysm

Towards the mechanical characterisation of unruptured intracranial aneurysms: Numerical modelling of interactions between a deformation device and the aneurysm wall

Intracranial aneurysm is a critical pathology related to the arterial wall deterioration. This work is an essential aspect of a large scale project aimed at providing clinicians with a non-invasive patient-specific decision support tool regarding the rupture risk assessment. A machine learning algorithm links the aneurysm shape observed and a database of UIA clinical images associated with in vivo wall mechanical properties and rupture characterisation. The database constitution is derived from a device prototype coupled with medical imaging. It provides the mechanical characterisation of the aneurysm from the wall deformation obtained by inverse analysis based on the variation of luminal volume. Before performing in vivo tests of the device on small animals, a numerical model was built to quantify the device's impact on the aneurysm wall under natural blood flow conditions. As the clinician will never be able to precisely situate the device, several locations were considered. In preparation for the inverse analysis procedure, artery material laws of increasing complexity were studied (linear elastic, hyper elastic Fung-like). Considering all the device locations and material laws, the device induced relative displacements to the Systole peak (worst case scenario with the highest mechanical stimulus linked to the blood flow) ranging from 375 μm to 1.28 mm. The variation of luminal volume associated with the displacements was between 0.95 % and 4.3 % compared to the initial Systole volume of the aneurysm. Significant increase of the relative displacements and volume variations were found with the study of different cardiac cycle moments between the blood flow alone and the device application. For forthcoming animal model studies, Spectral Photon CT Counting, with a minimum spatial resolution of 250 μm, was selected as the clinical imaging technique. Based on this preliminary study, the displacements and associated volume variations (baseline for inverse analyse), should be observable and exploitable.

Design, manufacturing and testing of a green non-isocyanate polyurethane prosthetic heart valve

The sole effective treatment for most patients with heart valve disease is valve replacement by implantation of mechanical or biological prostheses. However, mechanical valves represent high risk of thromboembolism, and biological prostheses are prone to early degeneration. In this work, we aim to determine the potential of novel environmentally-friendly non-isocyanate polyurethanes (NIPUs) for manufacturing synthetic prosthetic heart valves. Polyhydroxyurethane (PHU) NIPUs are synthesized via an isocyanate-free route, tested in vitro, and used to produce aortic valves. PHU elastomers reinforced with a polyester mesh show mechanical properties similar to native valve leaflets. These NIPUs do not cause hemolysis. Interestingly, both platelet adhesion and contact activation-induced coagulation are strongly reduced on NIPU surfaces, indicating low thrombogenicity. Fibroblasts and endothelial cells maintain normal growth and shape after indirect contact with NIPUs. Fluid-structure interaction (FSI) allows modeling of the ideal valve design, with minimal shear stress on the leaflets. Injection-molded valves are tested in a pulse duplicator and show ISO-compliant hydrodynamic performance, comparable to clinically-used bioprostheses. Poly(tetrahydrofuran) (PTHF)-NIPU patches do not show any evidence of calcification over a period of 8 weeks. NIPUs are promising sustainable biomaterials for the manufacturing of improved prosthetic valves with low thrombogenicity.

An assessment of how the anterior cerebral artery anatomy impacts ACoA aneurysm formation based on CFD analysis

OBJECTIVE

The aim of this study was to identify independent anatomic, morphologic and hemodynamic features of the ACoA (anterior communicating artery) complex that serve as risk factors for the occurrence of ACoA aneurysms.

METHODS

Fifteen consecutive patients with 15 ACoA aneurysms were included. Computational fluid dynamics (CFD) simulations based on patient-specific models were carried out using 3D time-of-flight magnetic resonance angiography (3D-TOF-MRA) images. A reverse reconstruction technique was used to generate a pre-aneurysm vessel anatomy. Geometric parameters and hemodynamic changes were compared and evaluated.

RESULTS

The overall prevalence of symmetric, dysplastic, and absent A1 segments were 53.3%, 26.7%, and 20%. The mean wall shear stress (WSS) of the absent group (AG) was significantly higher than that of the symmetric group (SG) and dysplastic group (DG). The absolute mean A1 artery flow rate (410.2 ± 88 versus 439.4 ± 101 mL/min; p = .45) of the aneurysm side was similar between the SG and DG but significantly higher in the AG (528.1 ± 77 mL/min; p < .05). The A1-A2 angles of the aneurysm side showed no significant differences among the 3 groups (p = .32). However, the mean A1-A2 angle on the aneurysm side was smaller than the contralateral A1-A2 angle (101.9 ± 9.1˚ versus 120.3 ± 7.7˚; p <.05). A regression analysis demonstrated that high WSS was significantly associated with a large A1-A2 ratio (R2=0.52; p <.05).

CONCLUSIONS

ACoA aneurysms are a high-WSS pathology. Severe flow impingement and the anatomic vasculature structures play a role in triggering the occurrence of ACoA aneurysms.

In vivo mechanical characterization of arterial wall using an inverse analysis procedure: application on an animal model of intracranial aneurysm

Intracranial aneurysm is a pathology related to the deterioration of the arterial wall. This work is an essential part of a large-scale project aimed at providing clinicians with a non-invasive patient-specific decision support tool to facilitate the rupture risk assessment. It will lean on the link between the aneurysm shape clinically observed and a database derived from the in vivo mechanical characterization of aneurysms. To supply this database, a deformation device prototype of the arterial wall was developed. Its use coupled with medical imaging (spectral photon-counting computed tomography providing a spatial resolution down to 250 μm) is used to determine the in vivo mechanical properties of the wall based on the inverse analysis of the quantification of the wall deformation observed experimentally. This study presents the in vivo application of this original procedure to an animal model of aneurysm. The mechanical properties of the aneurysm wall identified were consistent with the literature, and the errors between the numerical and experimental results were less than 10%. Based on these parameters, this study allows the assessment of the aneurysm stress state for a known solicitation and points towards the definition of a rupture criterion.

Morphology of Abdominal Aortic Aneurysms and Correlation with Biomechanical Tests of Aneurysmal Wall Fragments

BACKGROUND

Evaluate how specific morphologic aspects of abdominal aortic aneurysms (AAAs), including asymmetries, curvatures, tortuosities, and angulations, among others can influence the intrinsic biomechanical properties of the AAA's wall. This study analyzed the correlation of geometric measurements (1-dimensional, 2-dimensional, and 3-dimensional) of preoperative tomographic images of AAA with uniaxial biomechanical tests of the arterial wall fragments of these AAA obtained in open surgical repair of aneurysms.

METHODS

It was a multicenter, experimental, and observational study, and initially 54 individuals were selected who underwent open surgical of AAA, with valid biomechanical tests of the anterior wall of the AAA. Seven individuals were excluded because they had poor preoperative quality computed tomography scans and/or artifacts that impeded image segmentation and extraction of AAA geometric indices. The aortic fragments were subjected to uniaxial biomechanical destructive tests to obtain the following data: maximum load, failure stress, failure tension, failure strain energy, strain, and fragment thickness. In the same patients, preoperative computed tomography scans were performed with the extraction of 26 geometric indices, subdivided into 9 1-dimensional indices, 6 2-dimensional indices, and 11 3-dimensional indices. Data were subjected to statistical analysis using SPSS version 28.

RESULTS

Comparing ruptured and unruptured AAA, no statistical difference was observed between the biomechanical and geometric parameters. The fragment thickness of the ruptured AAA was lower than that of the unruptured AAA (P < 0.05). By comparing tomographic geometric indices and biomechanical parameters of the aortic fragments using Pearson's coefficient, positive and linear correlations (P < 0.05) were observed between the geometric variable maximum diameter (Dmax) of the AAA with maximum load (r = 0.408), failure tension (r = 0.372), and failure stress (r = 0.360). Positive and linear correlations were also observed between the variable diameter/height ratio (DHr) and the maximum load (r = 0.360), failure tension (r = 0.354), and failure stress (r = 0.289). The geometric variable DHr was dependent and correlated with Dmax. Simple regression analysis showed that R2 varied between 8.3% and 16.7%, and all models were significant (P < 0.05).

CONCLUSIONS

Dmax and DHr were linearly and positively correlated with the resistance parameters (maximum load, failure tension, and failure stress) of the AAA fragments. The DHr variable is dependent and correlated with Dmax. There was no correlation between the other geometric indices and the biomechanical parameters of the AAA wall. The asymmetries did not globally influence the biomechanics of AAA wall; however, they may influence regionally. Larger AAAs were stronger than smaller ones. Therefore, such findings may point toward Dmax is still the main geometric parameter, which influences the anterior wall, and possibly globally in the AAA.

Computational Study of Abdominal Aortic Aneurysm Walls Accounting for Patient-Specific Non-Uniform Intraluminal Thrombus Thickness and Distinct Material Models: A Pre- and Post-Rupture Case

An intraluminal thrombus (ILT) is present in the majority of abdominal aortic aneurysms, playing a crucial role in their growth and rupture. Although most computational studies do not include the ILT, in the present study, this is taken into account, laying out the whole simulation procedure, namely, from computed tomography scans to medical image segmentation, geometry reconstruction, mesh generation, biomaterial modeling, finite element analysis, and post-processing, all carried out in open software. By processing the tomography scans of a patient's aneurysm before and after rupture, digital twins are reconstructed assuming a uniform aortic wall thickness. The ILT and the aortic wall are assigned different biomaterial models; namely, the first is modeled as an isotropic linear elastic material, and the second is modeled as the Mooney-Rivlin hyperelastic material as well as the transversely isotropic hyperelastic Holzapfel-Gasser-Ogden nonlinear material. The implementation of the latter requires the designation of local Cartesian coordinate systems in the aortic wall, suitably oriented in space, for the proper orientation of the collagen fibers. The composite aneurysm geometries (ILT and aortic wall structures) are loaded with normal and hypertensive static intraluminal pressure. Based on the calculated stress and strain distributions, ILT seems to be protecting the aneurysm from a structural point of view, as the highest stresses appear in the thrombus-free areas of the aneurysmal wall.

Heat transfer mechanism in idealized healthy and diseased aortas using fluid-structure interaction method

The heat transfer mechanism inside the human aorta may be related to the physiological function and lesion formation of the aortic wall. The objective of this study was to acquire the temperature distribution in the three-dimensional idealized aorta. An idealized healthy aortic geometry and three representative diseased aortas: aortic aneurysm, coarctation of the aorta, and aortic dissection were constructed. Advanced fluid-structure interaction (FSI) computational framework was applied to predict the aortic temperature distribution. The movement of the aortic root due to the heartbeat was also considered. The displacement distribution of the aortic vessel wall was consistent with clinical observation. The lesser curvature of the aortic arch, aneurysm body, coarctation region, and false lumen were all exposed to relatively high temperatures (over 310.006 K). We found that the rigid wall assumption slightly underestimated the magnitude of the whole aortic wall-averaged temperature while the changing trend and local temperature were like the results of the FSI method. Besides, the wall-averaged temperature would increase and the temperature inflection point would advance when the aortic vessel wall was loaded with a high heat flux. This pilot study revealed the aortic heat transfer mechanism and temperature distribution, and the findings may help to understand the physiological characteristics of the aortic vessel wall.

A two-scale numerical study on the mechanobiology of abdominal aortic aneurysms

Abdominal aortic aneurysms (AAAs) are a serious condition whose pathophysiology is related to phenomena occurring at different length scales. To gain a better understanding of the disease, this work presents a multi-scale computational study that correlates AAA progression with microstructural and mechanical alterations in the tissue. Macro-scale geometries of a healthy aorta and idealized aneurysms with increasing diameter are developed on the basis of existing experimental data and subjected to physiological boundary conditions. Subsequently, microscopic representative volume elements of the abluminal side of each macro-model are employed to analyse the local kinematics at the cellular scale. The results suggest that the formation of the aneurysm disrupts the micromechanics of healthy tissue, which could trigger collagen growth and remodelling by mechanosensing cells. The resulting changes to the macro-mechanics and microstructure of the tissue seem to establish a new homeostatic state at the cellular scale, at least for the diameter range investigated.

The role of aorta distal to stent in the occurrence of distal stent graft-induced new entry tear: A computational fluid dynamics and morphological study

Distal stent graft-induced new entry tear (dSINE) is an important complication of thoracic endovascular aortic repair (TEVAR) for the treatment of type B aortic dissection (TBAD). This study aims to explore whether the aorta distal to the stent plays an important role in the occurrence of dSINE. Sixty-nine patient-specific geometrical models of twenty-three enrolled patients were reconstructed from preoperative, postoperative, and predSINE computed tomography scans. Computational fluid dynamics (CFD) simulations were performed to calculate the von Mises stress in the CFD group. Meanwhile, morphological measurements were performed in all patients, including measurements of the inverted pyramid index at different follow-up time points and the postoperative true lumen volume change rate. In the CFD study, the time-averaged von Mises stress of the true lumen distal to the stent in dSINE patients was significantly higher than that in the CFD controls (20.42 kPa vs. 15.47 kPa). In the morphological study, a special aortic plane (plane A) with an extremely small area distal to the stent was observed in dSINE patients, which resulted in an inverted pyramid structure in the true lumen distal to the stent. This structure in dSINE patients became increasingly obvious during the follow-up period and finally reached the maximum value before dSINE occurred (mean, 3.91 vs. 1.23). At the same time, enlargement of the true lumen distal to the stent occurs before dSINE, manifesting as a continuous increase in the true lumen volume (mean, 0.70 vs. 013). A new theory of what causes dSINE to occur has been proposed: the inverted pyramid structure of the true lumen distal to the stent caused an increase in the von Mises stress in this region and aortic enlargement, which ultimately led to the occurrence of dSINE.

Bio-inspired microfluidics: A review

Biomicrofluidics, a subdomain of microfluidics, has been inspired by several ideas from nature. However, while the basic inspiration for the same may be drawn from the living world, the translation of all relevant essential functionalities to an artificially engineered framework does not remain trivial. Here, we review the recent progress in bio-inspired microfluidic systems via harnessing the integration of experimental and simulation tools delving into the interface of engineering and biology. Development of "on-chip" technologies as well as their multifarious applications is subsequently discussed, accompanying the relevant advancements in materials and fabrication technology. Pointers toward new directions in research, including an amalgamated fusion of data-driven modeling (such as artificial intelligence and machine learning) and physics-based paradigm, to come up with a human physiological replica on a synthetic bio-chip with due accounting of personalized features, are suggested. These are likely to facilitate physiologically replicating disease modeling on an artificially engineered biochip as well as advance drug development and screening in an expedited route with the minimization of animal and human trials.

A review on the biomechanical behaviour of the aorta

Large aortic aneurysm and acute and chronic aortic dissection are pathologies of the aorta requiring surgery. Recent advances in medical intervention have improved patient outcomes; however, a clear understanding of the mechanisms leading to aortic failure and, hence, a better understanding of failure risk, is still missing. Biomechanical analysis of the aorta could provide insights into the development and progression of aortic abnormalities, giving clinicians a powerful tool in risk stratification. The complexity of the aortic system presents significant challenges for a biomechanical study and requires various approaches to analyse the aorta. To address this, here we present a holistic review of the biomechanical studies of the aorta by categorising articles into four broad approaches, namely theoretical, in vivo, experimental and combined investigations. Experimental studies that focus on identifying mechanical properties of the aortic tissue are also included. By reviewing the literature and discussing drawbacks, limitations and future challenges in each area, we hope to present a more complete picture of the state-of-the-art of aortic biomechanics to stimulate research on critical topics. Combining experimental modalities and computational approaches could lead to more comprehensive results in risk prediction for the aortic system.

AI-powered assessment of biomarkers for growth prediction of abdominal aortic aneurysms

Objective

The purpose of this study was to employ biomechanics-based biomarkers to locally characterize abdominal aortic aneurysm (AAA) tissue and investigate their relation to local aortic growth by means of an artificial intelligence model.

Methods

The study focused on a population of 36 patients with AAAs undergoing serial monitoring with electrocardiogram-gated multiphase computed tomography angiography acquisitions. The geometries of the aortic lumen and wall were reconstructed from the baseline scans and used for the baseline assessment of regional aortic weakness with three functional biomarkers, time-averaged wall-shear stress, in vivo principal strain, and intra-luminal thrombus thickness. The biomarkers were encoded as regional averages on axial and circumferential sections perpendicularly to the aortic centerline. Local diametric growth was obtained as difference in diameter between baseline and follow-up at the level of each axial section. An artificial intelligence model was developed to predict accelerated aneurysmal growth with the Extra Trees algorithm used as a binary classifier where the positive class represented regions that grew more than 2.5 mm/year. Additional clinical biomarkers, such as maximum aortic diameter at baseline, were also investigated as predictors of growth.

Results

The area under the curve for the constructed receiver operating characteristic curve for the Extra Trees classifier showed a very good performance in predicting relevant aortic growth (area under the curve = 0.92), with the three biomechanics-based functional biomarkers being objectively selected as the main predictors of growth.

Conclusions

The use of features based on the functional and local characterization of the aortic tissue resulted in a superior performance in terms of growth prediction when compared with models based on geometrical assessments. With rapid growth linked to increasing risk for patients with AAAs, the ability to access functional information related to tissue weakening and disease progression at baseline has the potential to support early clinical decisions and improve disease management.

Three-Dimensional Characterization of Aortic Root Motion by Vascular Deformation Mapping

The aorta is in constant motion due to the combination of cyclic loading and unloading with its mechanical coupling to the contractile left ventricle (LV) myocardium. This aortic root motion has been proposed as a marker for aortic disease progression. Aortic root motion extraction techniques have been mostly based on 2D image analysis and have thus lacked a rigorous description of the different components of aortic root motion (e.g., axial versus in-plane). In this study, we utilized a novel technique termed vascular deformation mapping (VDM(D)) to extract 3D aortic root motion from dynamic computed tomography angiography images. Aortic root displacement (axial and in-plane), area ratio and distensibility, axial tilt, aortic rotation, and LV/Ao angles were extracted and compared for four different subject groups: non-aneurysmal, TAA, Marfan, and repair. The repair group showed smaller aortic root displacement, aortic rotation, and distensibility than the other groups. The repair group was also the only group that showed a larger relative in-plane displacement than relative axial displacement. The Marfan group showed the largest heterogeneity in aortic root displacement, distensibility, and age. The non-aneurysmal group showed a negative correlation between age and distensibility, consistent with previous studies. Our results revealed a strong positive correlation between LV/Ao angle and relative axial displacement and a strong negative correlation between LV/Ao angle and relative in-plane displacement. VDM(D)-derived 3D aortic root motion can be used in future studies to define improved boundary conditions for aortic wall stress analysis.

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 study on the impact of AAA wall motion on the hemodynamics based on 4D CT image data

Purpose: To analyze the effect of the physiological deformation of the vessel wall on the hemodynamics in the abdominal aortic aneurysm (AAA), this paper compared the hemodynamics in AAA based on the moving boundary (MB) simulation and the rigid wall (RW) simulation.

Method: Patient-specific models were reconstructed to generate mesh based on four-dimensional computed tomography angiography (4D CT) data. The dynamic mesh technique was used to achieve deformation of the vessel wall, surface mesh and volume mesh of the fluid domain were successively remeshed at each time step. Besides, another rigid wall simulation was performed. Hemodynamics obtained from these two simulations were compared.

Results: Flow field and wall shear stress (WSS) distribution are similar. When using the moving boundary method (MBM), mean time-averaged wall shear stress (TAWSS) is lower, mean oscillatory shear index (OSI) and mean relative residence time (RRT) are higher. When using the 10th and 20th percentile values for TAWSS and 80th and 90th percentile values for RRT, the ratios of areas with low TAWSS, high OSI and high RRT to the entire vessel wall are higher than those assuming the vessel as rigid. In addition, one overlapping region of low TAWSS, high OSI and high RRT by using the MBM is consistent with the location of thrombus obtained from the follow-up imaging data.

Conclusion: The hemodynamics results by using the MBM reflect a higher blood retention effect. This paper presents a potential tool to assess the risk of intraluminal thrombus (ILT) formation based on the MBM.

How does hemodynamics affect rupture tissue mechanics in abdominal aortic aneurysm: Focus on wall shear stress derived parameters, time-averaged wall shear stress, oscillatory shear index, endothelial cell activation potential, and relative residence time

An abdominal aortic aneurysm (AAA) is a critical health condition with a risk of rupture, where the diameter of the aorta enlarges more than 50% of its normal diameter. The incidence rate of AAA has increased worldwide. Currently, about three out of every 100,000 people have aortic diseases. The diameter and geometry of AAAs influence the hemodynamic forces exerted on the arterial wall. Therefore, a reliable assessment of hemodynamics is crucial for predicting the rupture risk. Wall shear stress (WSS) is an important metric to define the level of the frictional force on the AAA wall. Excessive levels of WSS deteriorate the remodeling mechanism of the arteries and lead to abnormal conditions. At this point, WSS-related hemodynamic parameters, such as time-averaged WSS (TAWSS), oscillatory shear index (OSI), endothelial cell activation potential (ECAP), and relative residence time (RRT) provide important information to evaluate the shear environment on the AAA wall in detail. Calculation of these parameters is not straightforward and requires a physical understanding of what they represent. In addition, computational fluid dynamics (CFD) solvers do not readily calculate these parameters when hemodynamics is simulated. This review aims to explain the WSS-derived parameters focusing on how these represent different characteristics of disturbed hemodynamics. A representative case is presented for spatial and temporal formulation that would be useful for interested researchers for practical calculations. Finally, recent hemodynamics investigations relating WSS-related parameters with AAA rupture risk assessment are presented. This review will be useful to understand the physical representation of WSS-related parameters in cardiovascular flows and how they can be calculated practically for AAA investigations.

Study of flow, Bioheat transfer and cardiac thermal pulse of aneurysm in the abdominal aortic

Abdominal aortic aneurysms (AAA) are serious and difficult to detect conditions that can be deadly if they rupture. Infrared thermography (IRT) is a promising imaging technique that can detect abdominal aortic aneurysms more quickly and less costly than other imaging techniques. A clinical biomarker of circular thermal elevation on the midriff skin surface of AAA patient at various scenarios was expected during diagnosis using IRT scanner. 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 in detecting abdominal aortic aneurysms. Nevertheless, thermography is currently one of the most convenient technologies in imaging, and it has the potential to detect abdominal aortic aneurysms earlier than other techniques. Cardiac thermal pulse (CTP), on the other hand, was used to examine the thermal physics of AAA. AAA had a CTP that only responded to systolic phase at regular body temperature. Whereas the AAA wall would establish thermal homeostasis with blood temperature following a quasi-linear relationship as the body experienced fever or stage-2 hypothermia. In contrast, a healthy abdominal aorta displayed a CTP that responded to the full cardiac cycle, including diastolic phase at all simulated scenarios.

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.

Experimental and theoretical investigation of aortic wall tissue in tensile tests

BACKGROUND

Understanding the mechanical properties of aortic tissue is essential for developing numerical computation tools and assessing the risk of aortic aneurysm fractures. Tensile tests using aortic wall specimens allow for the determination of stress and strain depending on the location and direction of the sample.

OBJECTIVE

The aim of this study was to perform a mechanical tensile test using canine aorta samples and create a numerical model of aortic tissue tension from the processed data.

METHODS

Dogbone-shaped samples were dissected from canine aortic segments. The initial measurements were made at zero tension and the tensile tests were conducted at 10 mm/min until rupture. Force and stretch data were used to obtain engineering and true stress-strain curves. The true stress-strain curves were taken until the maximum strength was obtained, after which they were smoothed and fitted using a logistic function with three coefficients. These curves were then used as material mechanical properties for a numerical model of the aortic tissue tension. A simplified rectangle form was used to mimic the middle of the dogbone-shaped portion of the tissue specimen. Experimental displacement data were collected for the boundary conditions of the finite element 3D model.

RESULTS

The experimental data processing revealed that the logistic function described the nonlinear behaviour of the aorta soft tissue with an accuracy of 95% from the start of the tension to the media layer rupture. By applying numerical simulations, we obtained a correspondence of the load curve with an RMSE = 0.069 for the theoretical and experimental external tension data.

CONCLUSION

The numerical investigation confirmed that the non-linear soft tissue was validated by applying a logistic function approach to the mechanical properties of the aortic wall.

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.

Computational Analysis of Flow Structures in Turbulent Ventricular Blood Flow Associated With Mitral Valve Intervention

Cardiac disease and clinical intervention may both lead to an increased risk for thrombosis events due to a modified blood flow in the heart, and thereby a change in the mechanical stimuli of blood cells passing through the chambers of the heart. Specifically, the degree of platelet activation is influenced by the level and type of mechanical stresses in the blood flow. In this article we analyze the blood flow in the left ventricle of the heart through a computational model constructed from patient-specific data. The blood flow in the ventricle is modelled by the Navier-Stokes equations, and the flow through the mitral valve by a parameterized model which represents the projected opening of the valve. A finite element method is used to solve the equations, from which a simulation of the velocity and pressure of the blood flow is constructed. The intraventricular blood flow is complex, in particular in diastole when the inflow jet from the atrium breaks down into turbulent flow on a range of scales. A triple decomposition of the velocity gradient tensor is then used to distinguish between rigid body rotational flow, irrotational straining flow, and shear flow. The triple decomposition enables the separation of three fundamentally different flow structures, that each generates a distinct type of mechanical stimulus on the blood cells in the flow. We compare the results in a simulation where a mitral valve clip intervention is modelled, which leads to a significant modification of the intraventricular flow. Further, we perform a sensitivity study of the results with respect to the positioning of the clip. It was found that the shear in the simulation cases treated with clips increased more compared to the untreated case than the rotation and strain did. A decrease in valve opening area of 64% in one of the cases led to a 90% increase in rotation and strain, but a 150% increase in shear. The computational analysis opens up for improvements in models of shear-induced platelet activation, by offering an algorithm to distinguish shear from other modalities in intraventricular blood flow.

Component quantification of aortic blood flow energy loss using computational fluid-structure interaction hemodynamics

BACKGROUND AND OBJECTIVES

The aorta serves as the main tube of the human blood circulation system. Energy loss (EL) occurs when blood flows through the aorta and there may be a potential correlation between EL and aortic diseases. However, the components of blood flow EL are still not fully understood. This study aims to quantitatively reveal the EL components in healthy and diseased aortas.

METHODS

We construct an idealized healthy aorta and three idealized representative diseased aortas: aortic aneurysm, coarctation of the aorta, and aortic dissection. Computational hemodynamic studies are carried out by using the fluid-structure interaction simulation framework.

RESULTS

Four kinds of EL components: viscous friction, turbulence dissipation, wall deformation, and local lesion are firstly acquired in healthy and diseased aortas based on the high-resolution blood flow information. Viscous friction contributes most to the EL (45.69%-57.22%). EL caused by the deformation of the aortic wall ranks second (15.18%-33.12%). The proportions of turbulence dissipation and local lesion depend on individual geometric characteristics. Besides, the buffering efficiency of the healthy and diseased aorta is about 80%.

CONCLUSIONS

This study quantitatively reports the components of blood flow EL in healthy and diseased aortas, the finding may provide novel insights into the pathogenesis of aortic diseases.

Risk evaluation of type B aortic dissection based on WSS-based indicators distribution in different types of aortic arch

BACKGROUND AND OBJECTIVE

The underlying mechanism of aortic dissection (AD) remains unclear and the onset of AD is still unpredictable. Although clinical study with statistical analysis has reported that type III aortic arch may have strong correlation with type B AD (TBD), the effects of different arch types on the wall shear stress (WSS) have not been clarified.

METHODS

As a complementary work, this study numerically investigated the distribution of five WSS-based indicators in thirty aortic arches without AD, which were classified into three groups based on the arch types.

RESULTS

The distribution of most WSS indicators, such as time averaged WSS (TAWSS), oscillatory shear index (OSI) and relative residence time (RRT) had no significant difference among different types of aortic arches (P>0.05). However, a multidirectional WSS index, namely CFI, was found its maximum value was positively correlated with type III aortic arch in proximal descending aorta (p<0.001, r = 0.65).

CONCLUSIONS

It can be concluded that the enhancement or oscillation of WSS may not be the main reason of TBD is prevalence in type III arches, while the multidirectional WSS distribution may be an important factor. It can be further referred that the CFI may have a potential to predict the onset of TBD.

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.

Parallel generation of extensive vascular networks with application to an archetypal human kidney model

Given the relevance of the inextricable coupling between microcirculation and physiology, and the relation to organ function and disease progression, the construction of synthetic vascular networks for mathematical modelling and computer simulation is becoming an increasingly broad field of research. Building vascular networks that mimic in vivo morphometry is feasible through algorithms such as constrained constructive optimization (CCO) and variations. Nevertheless, these methods are limited by the maximum number of vessels to be generated due to the whole network update required at each vessel addition. In this work, we propose a CCO-based approach endowed with a domain decomposition strategy to concurrently create vascular networks. The performance of this approach is evaluated by analysing the agreement with the sequentially generated networks and studying the scalability when building vascular networks up to 200 000 vascular segments. Finally, we apply our method to vascularize a highly complex geometry corresponding to the cortex of a prototypical human kidney. The technique presented in this work enables the automatic generation of extensive vascular networks, removing the limitation from previous works. Thus, we can extend vascular networks (e.g. obtained from medical images) to pre-arteriolar level, yielding patient-specific whole-organ vascular models with an unprecedented level of detail.

Fluid-structure interaction: Insights into biomechanical implications of endograft after thoracic endovascular aortic repair

Thoracic endovascular aortic repair (TEVAR) has developed to be the most effective treatment for aortic diseases. This study aims to evaluate the biomechanical implications of the implanted endograft after TEVAR. We present a novel image-based, patient-specific, fluid-structure computational framework. The geometries of blood, endograft, and aortic wall were reconstructed based on clinical images. Patient-specific measurement data was collected to determine the parameters of the three-element Windkessel. We designed three postoperative scenarios with rigid wall assumption, blood-wall interaction, blood-endograft-wall interplay, respectively, where a two-way fluid-structure interaction (FSI) method was applied to predict the deformation of the composite stent-wall. Computational results were validated with Doppler ultrasound data. Results show that the rigid wall assumption fails to predict the waveforms of blood outflow and energy loss (EL). The complete storage and release process of blood flow energy, which consists of four phases is captured by the FSI method. The endograft implantation would weaken the buffer function of the aorta and reduce mean EL by 19.1%. The closed curve area of wall pressure and aortic volume could indicate the EL caused by the interaction between blood flow and wall deformation, which accounts for 68.8% of the total EL. Both the FSI and endograft have a slight effect on wall shear stress-related-indices. The deformability of the composite stent-wall region is remarkably limited by the endograft. Our results highlight the importance of considering the interaction between blood flow, the implanted endograft, and the aortic wall to acquire physiologically accurate hemodynamics in post-TEVAR computational studies and the deformation of the aortic wall is responsible for the major EL of the blood flow.

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.

Mechano-chemo-biological Computational Models for Arteries in Health, Disease and Healing: From Tissue Remodelling to Drug-eluting Devices

This review aims to highlight urgent priorities for the computational biomechanics community in the framework of mechano-chemo-biological models. Recent approaches, promising directions and open challenges on the computational modelling of arterial tissues in health and disease are introduced and investigated, together with in silico approaches for the analysis of drug-eluting stents that promote pharmacological-induced healing. The paper addresses a number of chemo-biological phenomena that are generally neglected in biomechanical engineering models but are most likely instrumental for the onset and the progression of arterial diseases. An interdisciplinary effort is thus encouraged for providing the tools for an effective in silico insight into medical problems. An integrated mechano-chemo-biological perspective is believed to be a fundamental missing piece for crossing the bridge between computational engineering and life sciences, and for bringing computational biomechanics into medical research and clinical practice.

Comparative evaluation of ballet-type and conventional stent graft configurations for endovascular aneurysm repair: A CFD analysis

PURPOSE:

Cross limb stent graft (SG) configuration technique for endovascular aneurysm repair (EVAR) is employed for splayed aortic bifurcations to avoid device kinking and smoothen cannulation. The present study investigates three types of stent graft (SG) configurations for endovascular aneurysm repair (EVAR) in abdominal aortic aneurysm. A computational fluid dynamic analysis was performed on the pulsatile non-Newtonian flow characteristics in three ideally modeled geometries of abdominal aortic (AA) SG configurations.

METHODS:

The three planar and crosslimb SG configurations were ideally modeled, namely, top-down nonballet-type, top-down ballet-type, and bottom-up nonballet-type configurations. In top-down SG configuration, most of the device is deployed in the main body in the vicinity of renal artery and the limbs are extended to the iliac artery. While in the bottom-up configuration, some of the SG device is deployed in the main body, the limbs are deployed in aortic bifurcation, and the extra stent graft of the main body is extended to the proximal aorta until the below of the renal artery. The effects of non-Newtonian pulsatile flow on the wall stresses and flow patterns of the three models were investigated and compared. Moreover, the average wall shear stress (AWSS), oscillatory shear stress index (OSI), absolute helicity, pressure distribution, graft displacement and flow visualization plots were analyzed.

RESULTS:

The top-down ballet-type showed less branch blockage effect than the top-down nonballet-type models. Furthermore, the top-down ballet-type configuration showed an increased tendency to sustain high WSS and higher helicity characteristics than that of the bottom-up and top-down non-ballet type configurations. However, displacement forces of the top-down ballet-type configuration were 40% and 9.6% higher than those of the bottom-up and top-down nonballet-type configurations, respectively.

CONCLUSIONS:

Some complications such as graft tearing, thrombus formation, limb disconnection during long term follow up periods might be relevant to hemodynamic characteristics according to the configurations of EVAR. Hence, the reported data required to be validated with the clinical results.

Aortic hemodynamics assessment prior and after valve sparing reconstruction: A patient-specific 4D flow-based FSI model

INTRODUCTION

Valve-sparing root replacement (VSRR) of the ascending aorta is a life-saving procedure for the treatment of aortic aneurysms, but patients remain at risk for post-operative events involving the downstream native aorta, the mechanism for which is uncertain. It is possible that proximal graft replacement of the ascending aorta induces hemodynamics alterations in the descending aorta, which could trigger adverse events. Herein, we present a fluid-structure interaction (FSI) protocol, based on patient-specific geometry and boundary conditions, to assess impact of proximal aortic grafts on downstream aortic hemodynamics and distensibility.

METHODS

Cardiac magnetic resonance (CMR), including MRA, cine-CMR and 4D flow sequences, was performed prior and after VSRR on one subject. Central blood pressure was non-invasively acquired at the time of the CMR: data were used to reconstruct the pre- and post-VSRR model and derive patient-specific boundary conditions for the FSI and a computational fluid dynamic (CFD) analysis with the same settings. Results were validated comparing the predicted velocity field against 4D flow dataset, over four landmarks along the aorta, and the predicted distensibility against the cine-CMR derived value.

RESULTS

Instantaneous velocity magnitudes extracted from 4D flow and FSI were similar (p > 0.05), while CFD-predicted velocity was significantly higher (p < 0.001), especially in the descending aorta of the pre-VSRR model (v_max was 73 cm/s, 76 cm/s and 99 cm/s, respectively). As measured in cine-CMR, FSI predicted an increase in descending aorta distensibility after grafting (i.e., 4.02 to 5.79 10^-3 mmHg^-1). In the descending aorta, the post-VSRR model showed increased velocity, aortic distensibility, stress and strain and wall shear stress.

CONCLUSIONS

Our Results indicate that i) the distensibility of the wall cannot be neglected, and hence the FSI method is necessary to obtain reliable results; ii) graft implantation induces alterations in the hemodynamics and biomechanics along the thoracic aorta, that may trigger adverse vessel remodeling.

Application of Proper Orthogonal Decomposition to Study Coherent Flow Structures in a Saccular Aneurysm

Aneurysms are localized expansions of weakened blood vessels that can be debilitating or fatal upon rupture. Previous studies have shown that flow in an aneurysm exhibits complex flow structures that are correlated with its inflow conditions. Therefore, the objective of this study was to demonstrate the application of proper orthogonal decomposition (POD) to study the impact of different inflow conditions on energetic flow structures and their temporal behavior in an aneurysm. To achieve this objective, experiments were performed on an idealized rigid sidewall aneurysm model. A piston pump system was used for precise inflow control, i.e., peak Reynolds number (Rep) and Womersley number (α) were varied from 50 to 270 and 2 to 5, respectively. The velocity flow field measurements at the midplane location of the idealized aneurysm model were performed using particle image velocimetry (PIV). The results demonstrate the efficacy of POD in decomposing complex data, and POD was able to capture the energetic flow structures unique to each studied inflow condition. Furthermore, the time-varying coefficient results highlighted the interplay between the coefficients and their corresponding POD modes, which in turn helped explain how POD modes impact certain flow features. The low-order reconstruction results were able to capture the flow evolution and provide information on complex flow in an aneurysm. The POD and low-order reconstruction results also indicated that vortex formation, evolution, and convection varied with an increase in α, while vortex strength and formation of secondary structures were correlated with an increase in Rep.

NUMERICAL COUPLING ANALYSIS OF THE INFLUENCE OF BLOOD FLOW ON THE MECHANICAL RESPONSE FOR LIVER

As an important basis for determining the state of the liver, the mechanical responses are associated with many factors, and belong to a complex coupling system. Liver tissue has significantly complicated vascular channels. The vascular diameter, vascular deflection angle and vascular depth are defined as the key characteristic parameters. The influences of these parameters on the mechanical responses were analyzed. On the basis of the real mechanical parameters, the coupled numerical model of blood vessel, blood flow and liver tissue was established. The corresponding mechanical responses are obtained by utilizing the different vascular parameters. The effects of vascular parameters on the differences among the mechanical response difference and high strain modulus were analyzed. It was found that the blood vessels in the central area could reduce the liver mechanical response. The inner diameter parameter had main influences on the regions where the stain was more than 0.1. The mechanical difference is greater with larger inner diameter. The influences of vascular depth are greatest when the vascular depth was in the intermediate value, which would increase the liver mechanical responses. With the increment of vascular deflection angle, the liver mechanical response would also increase, and exceed the mechanical response without blood vessels. The findings after analyzing the influence of vascular parameters will provide a basis for the quantitative studies on the influence of blood vessels.

Linking Aortic Mechanical Properties, Gene Expression and Microstructure: A New Perspective on Regional Weakening in Abdominal Aortic Aneurysms

Background: Current clinical practice for the assessment of abdominal aortic aneurysms (AAA) is based on vessel diameter and does not account for the multifactorial, heterogeneous remodeling that results in the regional weakening of the aortic wall leading to aortic growth and rupture. The present study was conducted to determine correlations between a novel non-invasive surrogate measure of regional aortic weakening and the results from invasive analyses performed on corresponding ex vivo aortic samples. Tissue samples were evaluated to classify local wall weakening and the likelihood of further degeneration based on non-invasive indices.

Methods: A combined, image-based fluid dynamic and in-vivo strain analysis approach was used to estimate the Regional Aortic Weakness (RAW) index and assess individual aortas of AAA patients prior to elective surgery. Nine patients were treated with complete aortic resection allowing the systematic collection of tissue samples that were used to determine regional aortic mechanics, microstructure and gene expression by means of mechanical testing, microscopy and transcriptomic analyses.

Results: The RAW index was significantly higher for samples exhibiting lower mechanical strength (p = 0.035) and samples classified as low elastin content (p = 0.020). Samples with higher RAW index had the greatest number of genes differentially expressed compared to any constitutive metric. High RAW samples showed a decrease in gene expression for elastin and a down-regulation of pathways responsible for cell movement, reorganization of cytoskeleton, and angiogenesis.

Conclusions: This work describes the first AAA index free of assumptions for material properties and accounting for patient-specific mechanical behavior in relation to aneurysm strength. Use of the RAW index captured biomechanical changes linked to the weakening of the aorta and revealed changes in microstructure and gene expression. This approach has the potential to provide an improved tool to aid clinical decision-making in the management of aortic pathology.

Numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation

Background

The effects of arterial wall compliance on blood flow have been revealed using fluid-structure interaction in last decades. However, microcirculation is not considered in previous researches. In fact, microcirculation plays a key role in regulating blood flow. Therefore, it is very necessary to involve microcirculation in arterial hemodynamics.

Objective

The main purpose of the present study is to investigate how wall compliance affects the flow characteristics and to establish the comparisons of these flow variables with rigid wall when microcirculation is considered.

Methods

We present numerical modeling in arterial hemodynamics incorporating fluid-structure interaction and microcirculation. A novel outlet boundary condition is employed to prescribe microcirculation in an idealised model.

Results

The novel finding in this work is that wall compliance under the consideration of microcirculation leads to the increase of wall shear stress in contrast to rigid wall, contrary to the traditional result that wall compliance makes wall shear stress decrease when a constant or time dependent pressure is specified at an outlet.

Conclusions

This work provides the valuable study of hemodynamics under physiological and realistic boundary conditions and proves that wall compliance may have a positive impact on wall shear stress based on this model. This methodology in this paper could be used in real model simulations.

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.

Numerical modeling and simulations of type B aortic dissection treated by stent-grafts with different oversizing ratios

Retrograde type A dissection after thoracic endovascular aortic repair has been a major drawback of endovascular treatment. This study investigated the biomechanical mechanism of stent-graft-induced new lesions after implantation and analyzed the relationship between radial force and spring-back force of the stent-graft when it was implanted virtually under different oversizing ratios. Based on the computed tomography angiography images, a three-dimensional geometric model of a patient-specific aortic dissection was established. The stent was designed in CAD software and the stent-graft implantation procedure under different oversizing ratios was simulated in the finite element analysis software. Implantation simulations were performed six times for each stent-graft model under 0%, 3%, 6%, 9%, 12%, and 15% oversizing ratios and the peak stress of the aorta was compared among groups. It was observed that the peak stress of the aorta was located where the proximal bare stent interacted with aortic wall and its value was increased by 62.2% from 0% to 15% oversizing ratio. The conclusions are reached that the long-term higher stress in the aortic wall may lead to the emergence of new lesions in these areas, and the radial force plays a key role in the formation of a new entry in the real aorta model.

Computational prediction of hemodynamical and biomechanical alterations induced by aneurysm dilatation in patient-specific ascending thoracic aortas

The aim of the present work is to propose a robust computational framework combining computational fluid dynamics (CFD) and 4D flow MRI to predict the progressive changes in hemodynamics and wall rupture index (RPI) induced by aortic morphological evolutions in patients harboring ascending thoracic aortic aneurysms (ATAAs). An analytical equation has been proposed to predict the aneurysm progression based on age, sex, and body surface area. Parameters such as helicity, wall shear stress (WSS), time-averaged WSS, oscillatory shear index, relative residence time, and viscosity were evaluated for two patients at different stages of aneurysm growth, and compared with age-sex-matched healthy subjects. The study shows that evolution of hemodynamics and RPI, despite being very slow in ATAAs, is strongly affected by morphological alterations and, in turn could impact biomechanical factors and aortic mechanobiology. An aspect of the current work is that the patient-specific 4D MRI data sets were obtained with a follow-up of 1 year and the measured time-averaged velocity maps and flow eccentricity were compared with the CFD simulation for validation. The computational framework presented here is capable of capturing the blood flow patterns and the hemodynamic descriptors during the various stages of aneurysm growth. Further investigations will be conducted in order to verify these results on a larger cohort of patients and with long follow-up times to finally elucidate the link between deranged hemodynamics, AA geometry, and wall mechanical properties in ATAAs.

A novel combined fluid dynamic and strain analysis approach identified abdominal aortic aneurysm rupture

Clinical decision-making for surgical repair of abdominal aortic aneurysms based on maximum aortic diameter presents limitations as rupture can occur below threshold for some aneurysms, whereas others are stable at large sizes. The proposed approach combines hemodynamics and geometric indices with in vivo deformation analysis to assess local weakening of the aortic wall for a case of impending rupture that was confirmed during open surgical repair. A new combined index, the Regional Rupture Potential, is introduced to help the assessment of individual aneurysms and their rupture risk, providing a rationale for clinical decisions.

Computational analysis of Nitinol stent-graft for endovascular aortic repair (EVAR) of abdominal aortic aneurysm (AAA): Crimping, sealing and fluid-structure interaction (FSI)

OBJECTIVES

We evaluate the crimping strain, sealing stress and contact forces on a Nitinol stent deployed in the aorta during endovascular aortic (or aneurysm) repair (EVAR) procedures. Nitinol shape memory effect (SME) is used. We also study the fluid-structure interaction (FSI) of the blood flow on the stented aorta.

METHODS

We employ Solidworks to generate a closed-cell honeycomb stent structure used to treat abdominal aortic aneurysm (AAA). We use the commercial Abaqus/Simulia finite element (FEM) simulation package to study the displacements and stresses experienced by the stent during the crimping phase and deployment into the aortic segment. The Nitinol stent is covered with Dacron, a popular graft material. We implement our own user-material (UMAT) subroutines to model the shape memory effect (SME) of Nitinol. The effect of the stent geometry is analyzed. We use the FSI analysis in Abaqus/Simulia to understand the effect of hemodynamic loading on the stent.

RESULTS

Results indicate that the crimping strain increases as the stent strut spacing increases. This is also the case for the radius of curvature. Maximum strains developed on the stent during crimping are in the order of 10%. Stresses exerted by the stent needed to completely seal the aorta are found to be below the yield stress values of Nitinol (700 MPa). Wall shear stresses (WSS) on the stented aorta are close to WSS obtained on the aorta alone.

CONCLUSION

Using Nitinol's thermo-reactivity property as opposed to its superelasticity causes the stent-graft to deploy more gently.

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.

A novel method for describing biomechanical properties of the aortic wall based on the three-dimensional fluid-structure interaction model

OBJECTIVES

Our goal was to present a novel non-invasive approach for assessment of aortic wall displacement to describe its biomechanical properties during the cardiac cycle.

METHODS

The fluid-structure interaction (FSI) technique was used to reconstruct aortic wall displacement based on computed tomography angiography and 2-dimensional speckle-tracking technique (2DSTT) data collected from 20 patients [10 with healthy aortas (AA) and 10 with abdominal aortic aneurysms (AAAs)]. The mechanical properties of the wall of the aorta were described by the Yeoh hyperelastic materials model with α and β parameters, and wall displacement was determined with 2DSTT. The mechanical parameters of the wall of the aorta in the FSI model were automatically updated in the calculation loop until the calculated and clinically measured wall movements were the same.

RESULTS

Results showed 98% accuracy of FSI compared to 2DSTT for AA and AAA (P > 0.05). The mean wall deformation for AA was 2.45 ± 0.12 mm and 2.49 ± 0.10 mm for FSI and 2DSTT, respectively (P = 0.40), whereas that for AAA was 2.84 ± 0.44 mm and 2.88 ± 0.45 mm, respectively (P = 0.83). The FSI analysis indicated that the α and β parameters for AA were equal to 14.35 ± 1.30 N⋅cm-2 and 9.33 ± 1.08 N⋅cm-2, respectively; and for AAA, α was 11.00 ± 0.49 N⋅cm-2 and β was 79.46 ± 4.32 N⋅cm-2.

CONCLUSIONS

The FSI technique may be successfully applied to assess the mechanical parameters of patient-specific aortic walls using computed tomography angiographic and 2DSTT measurements.

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.

Computational study on hemodynamic changes in patient-specific proximal neck angulation of abdominal aortic aneurysm with time-varying velocity

Aneurysms are considered as a critical cardiovascular disease worldwide when they rupture. The clinical understanding of geometrical impact on the flow behaviour and biomechanics of abdominal aortic aneurysm (AAA) is progressively developing. Proximal neck angulations of AAAs are believed to influence the hemodynamic changes and wall shear stress (WSS) within AAAs. Our aim was to perform pulsatile simulations using computational fluid dynamics (CFD) for patient-specific geometry to investigate the influence of severe angular (≥ 60°) neck on AAA's hemodynamic and wall shear stress. The patient's geometrical characteristics were obtained from a computed tomography images database of AAA patients. The AAA geometry was reconstructed using Mimics software. In computational method, blood was assumed Newtonian fluid and an inlet varying velocity waveform in a cardiac cycle was assigned. The CFD study was performed with ANSYS software. The results of flow behaviours indicated that the blood flow through severe bending of angular neck leads to high turbulence and asymmetry of flows within the aneurysm sac resulting in blood recirculation. The high wall shear stress (WSS) occurred near the AAA neck and on surface of aneurysm sac. This study explained and showed flow behaviours and WSS progression within high angular neck AAA and risk prediction of abdominal aorta rupture. We expect that the visualization of blood flow and hemodynamic changes resulted from CFD simulation could be as an extra tool to assist clinicians during a decision making when estimation the risks of interventional procedures.