A patient-specific computational model of fluid–structure interaction in abdominal aortic aneurysms

Calibration of the Mechanical Boundary Conditions for a Patient-Specific Thoracic Aorta Model Including the Heart Motion Effect

OBJECTIVE

We propose a procedure for calibrating 4 parameters governing the mechanical boundary conditions (BCs) of a thoracic aorta (TA) model derived from one patient with ascending aortic aneurysm. The BCs reproduce the visco-elastic structural support provided by the soft tissue and the spine and allow for the inclusion of the heart motion effect.

METHODS

We first segment the TA from magnetic resonance imaging (MRI) angiography and derive the heart motion by tracking the aortic annulus from cine-MRI. A rigid-wall fluid-dynamic simulation is performed to derive the time-varying wall pressure field. We build the finite element model considering patient-specific material properties and imposing the derived pressure field and the motion at the annulus boundary. The calibration, which involves the zero-pressure state computation, is based on purely structural simulations. After obtaining the vessel boundaries from the cine-MRI sequences, an iterative procedure is performed to minimize the distance between them and the corresponding boundaries derived from the deformed structural model. A strongly-coupled fluid-structure interaction (FSI) analysis is finally performed with the tuned parameters and compared to the purely structural simulation.

RESULTS AND CONCLUSION

The calibration with structural simulations allows to reduce maximum and mean distances between image-derived and simulation-derived boundaries from 8.64 mm to 6.37 mm and from 2.24 mm to 1.83 mm, respectively. The maximum root mean square error between the deformed structural and FSI surface meshes is 0.19 mm. This procedure may prove crucial for increasing the model fidelity in replicating the real aortic root kinematics.

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.

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.

Inverse material parameter estimation of patient-specific finite element models at the carotid bifurcation: The impact of excluding the zero-pressure configuration and residual stress

The carotid bifurcation experiences a complex loading environment due to its anatomical structure. Previous in-vivo material parameter estimation methods often use simplified model geometries, isotropic hyperelastic constitutive equations or neglect key aspects of the vessel, such as the zero-pressure configuration or residual stress, all of which have independently been shown to alter the stress environment of the vessel wall. Characterizing the location of high stress in the vessel wall has often been proposed as a potential indicator of structural weakness. However, excluding the afore-mentioned zero-pressure configuration, residual stress and patient-specific material parameters can lead to an incorrect estimation of the true stress values observed, meaning that stress alone as a risk indicator of rupture is insufficient. In this study, we investigate how the estimated material parameters and overall stress distributions in geometries of carotid bifurcations, extracted from in-vivo MR images, alter with the inclusion of the zero-pressure configuration and residual stress. This approach consists of the following steps: (1) geometry segmentation and hexahedral meshing from in-vivo magnetic resonance images (MRI) at two known phases; (2) computation of the zero-pressure configuration and the associated residual stresses; (3) minimization of an objective function built on the difference between the stress states of an "almost true" stress field at two known phases and a "deformed" stress field by altering the input material parameters to determine patient-specific material properties; and (4) comparison of the stress distributions throughout these carotid bifurcations for all cases with estimated material parameters. This numerical approach provides insights into the need for estimation of both the zero-pressure configuration and residual stress for accurate material property estimation and stress analysis for the carotid bifurcation, establishing the reliability of stress as a rupture risk metric.

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.

Studying the effect of stent thickness and porosity on post-stent implantation hemodynamics

This study investigates the effect of stent thickness and stent porosity which are important factors determining the post-treatment intra-aneurysmal hemodynamics. The study uses computational fluid dynamics (CFD) to estimate the hemodynamic behaviour: flow velocity, pressure distributions, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), besides relative residence time (RRT) blood flow distribution in a proposed stent and three other commercially available stents. The hemodynamic behaviour is compared between four different cases. In each case, each stent has the specific thickness and porosity values. The results show that the velocity magnitude inside the sac declined in thinner stents and lower porosity stents, TAWSS on the aneurysmal wall declined linearly in lower porosity stents, OSI and RRT increased obviously in thicker stents and higher porosity stents. Finally, the results conclude that the stent with the lowest thickness and porosity has the best performance that leads to post-stent thrombus formation and healing. However, the proposed stent design, a more porous construct, has a higher RRT compared to the used commercially available stents, which helps promote the thrombus growth inside the aneurysm sac.

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.

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.

Fluid Structure Interaction on Paravalvular Leakage of Transcatheter Aortic Valve Implantation Related to Aortic Stenosis: A Patient-Specific Case

This study investigated the impact of paravalvular leakage (PVL) in relation to the different valve openings of the transcatheter aortic valve implantation (TAVI) valve using the fluid structure interaction (FSI) approach. Limited studies were found on the subject of FSI with regards to TAVI-PVL condition, which involves both fluid and structural responses in coupling interaction. Hence, further FSI simulation with the two-way coupling method is implemented to investigate the effects of hemodynamics blood flow along the patient-specific aorta model subjected to the interrelationship between PVL and the different valve openings using the established FSI software ANSYS 16.1. A 3D patient-specific aorta model is constructed using MIMICS software. The TAVI valve identical to Edward SAPIEN XT 26 (Edwards Lifesciences, Irvine, California), at different Geometrical Orifice Areas (GOAs), is implanted into the patient's aortic annulus. The leaflet opening of the TAVI valve is drawn according to severity of GOA opening represented in terms of 100%, 80%, 60%, and 40% opening, respectively. The result proved that the smallest percentage of GOA opening produced the highest possibility of PVL, increased the recirculatory flow proximally to the inner wall of the ascending aorta, and produced lower backflow velocity streamlines through the side area of PVL region. Overall, 40% GOA produced 89.17% increment of maximum velocity magnitude, 19.97% of pressure drop, 65.70% of maximum WSS magnitude, and a decrement of 33.62% total displacement magnitude with respect to the 100% GOA.

Finite Element Framework for Computational Fluid Dynamics in FEBio

The mechanics of biological fluids is an important topic in biomechanics, often requiring the use of computational tools to analyze problems with realistic geometries and material properties. This study describes the formulation and implementation of a finite element framework for computational fluid dynamics (CFD) in FEBio, a free software designed to meet the computational needs of the biomechanics and biophysics communities. This formulation models nearly incompressible flow with a compressible isothermal formulation that uses a physically realistic value for the fluid bulk modulus. It employs fluid velocity and dilatation as essential variables: The virtual work integral enforces the balance of linear momentum and the kinematic constraint between fluid velocity and dilatation, while fluid density varies with dilatation as prescribed by the axiom of mass balance. Using this approach, equal-order interpolations may be used for both essential variables over each element, contrary to traditional mixed formulations that must explicitly satisfy the inf-sup condition. The formulation accommodates Newtonian and non-Newtonian viscous responses as well as inviscid fluids. The efficiency of numerical solutions is enhanced using Broyden's quasi-Newton method. The results of finite element simulations were verified using well-documented benchmark problems as well as comparisons with other free and commercial codes. These analyses demonstrated that the novel formulation introduced in FEBio could successfully reproduce the results of other codes. The analogy between this CFD formulation and standard finite element formulations for solid mechanics makes it suitable for future extension to fluid-structure interactions (FSIs).

Morphology and Hemodynamics in Isolated Common Iliac Artery Aneurysms Impacts Proximal Aortic Remodeling

Objective—

Isolated common iliac artery aneurysms (CIAA) are rare. Their prognosis and influence on aortoiliac blood flow and remodeling are unclear. We evaluated the hypotheses that morphology at and distal to the aortic bifurcation, together with the associated hemodynamic changes, influence both the natural history of CIAA and proximal aortic remodeling.

Approach and Results—

Twenty-five isolated CIAAs (15 intact, 10 ruptured), in 23 patients were reconstructed and analyzed with computational fluid dynamics: all showed abnormal flow. Then we studied a series of 24 hypothetical aortoiliac geometries in silico with varying abdominal aortic deflection and aortic bifurcation angles: key findings were assessed in an independent validation cohort of 162 patients. Wall shear stress in isolated unilateral CIAAs was lower than the contralateral common iliac artery, 0.38±0.33 Pa versus 0.61±0.24 Pa, inversely associated with CIAA diameter ( P <0.001) and morphology (high shear stress in variants distal to a sharp kink). Rupture usually occurred in regions of elevated low and oscillatory shear with a wide aortic bifurcation angle. Abdominal aortas deflected towards the CIAA for most unilateral isolated CIAAs (14/21). In silico, wider bifurcation angles created high focal regions of low and oscillatory shear in the common iliac artery. The associations of unilateral CIAA with aortic deflection and common iliac artery diameter with bifurcation angle were confirmed in the validation cohort.

Conclusions—

Decreasing wall shear stress is strongly associated with CIAA progression (larger aneurysms and rupture), whereas abnormal blood flow in the CIAA seems to promote proximal aortic remodeling, with adaptive lateral deflection of the abdominal aorta towards the aneurysmal side.

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.

Uncertainty in cardiac myofiber orientation and stiffnesses dominate the variability of left ventricle deformation response

Computational cardiac modelling is a mature area of biomedical computing and is currently evolving from a pure research tool to aiding in clinical decision making. Assessing the reliability of computational model predictions is a key factor for clinical use, and uncertainty quantification (UQ) and sensitivity analysis are important parts of such an assessment. In this study, we apply UQ in computational heart mechanics to study uncertainty both in material parameters characterizing global myocardial stiffness and in the local muscle fiber orientation that governs tissue anisotropy. The uncertainty analysis is performed using the polynomial chaos expansion (PCE) method, which is a nonintrusive meta-modeling technique that surrogates the original computational model with a series of orthonormal polynomials over the random input parameter space. In addition, in order to study variability in the muscle fiber architecture, we model the uncertainty in orientation of the fiber field as an approximated random field using a truncated Karhunen-Loéve expansion. The results from the UQ and sensitivity analysis identify clear differences in the impact of various material parameters on global output quantities. Furthermore, our analysis of random field variations in the fiber architecture demonstrate a substantial impact of fiber angle variations on the selected outputs, highlighting the need for accurate assignment of fiber orientation in computational heart mechanics models.

A Formulation for Fluid Structure-Interactions in FEBio Using Mixture Theory

Many physiological systems involve strong interactions between fluids and solids, posing a signicant challenge when modeling biomechanics. The objective of this study was to implement a fluid-structure interaction (FSI) solver in the free, open-source finite element code FEBio (febio.org), that combined the existing solid mechanics and rigid body dynamics solver with a recently-developed computational fluid dynamics (CFD) solver. A novel Galerkin-based finite element FSI formulation was introduced based on mixture theory, where the FSI domain was described as a mixture of fluid and solid constituents that have distinct motions. The mesh was defined on the solid domain, specialized to have zero mass, negligible stiffness and zero frictional interactions with the fluid, whereas the fluid was modeled as isothermal and compressible. The mixture framework provided the foundation for evaluating material time derivatives in a material frame for the solid and in a spatial frame for the fluid. Similar to our recently reported CFD solver, our FSI formulation did not require stabilization methods to achieve good convergence, producing a compact set of equations and code implementation. The code was successfully verified against benchmark problems and an analytical solution for squeeze-film lubrication. It was validated against experimental measurements of the flow rate in a peristaltic pump, and illustrated using non-Newtonian blood flow through a bifurcated carotid artery with a thick arterial wall. The successful formulation and implementation of this FSI solver enhances the multiphysics modeling capabilities in FEBio relevant to the biomechanics and biophysics communities.

Patient-specific CFD modelling in the thoracic aorta with PC-MRI-based boundary conditions: A least-square three-element Windkessel approach

The increasing use of computational fluid dynamics for simulating blood flow in clinics demands the identification of appropriate patient-specific boundary conditions for the customization of the mathematical models. These conditions should ideally be retrieved from measurements. However, finite resolution of devices as well as other practical/ethical reasons prevent the construction of complete data sets necessary to make the mathematical problems well posed. Available data need to be completed by modelling assumptions, whose impact on the final solution has to be carefully addressed. Focusing on aortic vascular districts and related pathologies, we present here a method for efficiently and robustly prescribing phase contrast MRI-based patient-specific data as boundary conditions at the domain of interest. In particular, for the outlets, the basic idea is to obtain pressure conditions from an appropriate elaboration of available flow rates on the basis of a 3D/0D dimensionally heterogeneous modelling. The key point is that the parameters are obtained by a constrained optimization procedure. The rationale is that pressure conditions have a reduced impact on the numerical solution compared with velocity conditions, yielding a simulation framework less exposed to noise and inconsistency of the data, as well as to the arbitrariness of the underlying modelling assumptions. Numerical results confirm the reliability of the approach in comparison with other patient-specific approaches adopted in the literature.

Structural modelling of the cardiovascular system

Computational modelling of the cardiovascular system offers much promise, but represents a truly interdisciplinary challenge, requiring knowledge of physiology, mechanics of materials, fluid dynamics and biochemistry. This paper aims to provide a summary of the recent advances in cardiovascular structural modelling, including the numerical methods, main constitutive models and modelling procedures developed to represent cardiovascular structures and pathologies across a broad range of length and timescales; serving as an accessible point of reference to newcomers to the field. The class of so-called hyperelastic materials provides the theoretical foundation for the modelling of how these materials deform under load, and so an overview of these models is provided; comparing classical to application-specific phenomenological models. The physiology is split into components and pathologies of the cardiovascular system and linked back to constitutive modelling developments, identifying current state of the art in modelling procedures from both clinical and engineering sources. Models which have originally been derived for one application and scale are shown to be used for an increasing range and for similar applications. The trend for such approaches is discussed in the context of increasing availability of high performance computing resources, where in some cases computer hardware can impact the choice of modelling approach used.

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.

An atlas- and data-driven approach to initializing reaction-diffusion systems in computer cardiac electrophysiology

The cardiac electrophysiology (EP) problem is governed by a nonlinear anisotropic reaction-diffusion system with a very rapidly varying reaction term associated with the transmembrane cell current. The nonlinearity associated with the cell models requires a stabilization process before any simulation is performed. More importantly, when used in a 3-dimensional (3D) anatomy, it is not sufficient to perform this stabilization on the basis of isolated cells only, since the coupling of the different cells through the tissue greatly modulates the dynamics of the system. Therefore, stabilization of the system must be performed on the entire 3D model. This work develops a novel procedure for the initialization of reaction-diffusion systems for numerical simulations of cardiac EP from steady-state conditions. We exploit surface point correspondence to establish volumetric point correspondence. Upon introduction of a new 3D anatomy with surface point correspondence, a prediction of the cell model steady states is derived from the set of earlier biophysical simulations. We show that the prediction error is typically less than 10% for all model variables, with most variables showing even greater accuracy. When initializing simulations with the predicted model states, it is demonstrated that simulation times can be cut by at least two-thirds and potentially more, which saves hours or days of high-performance computing. Overall, these results increase the clinical applicability of detailed computational EP studies on personalized anatomies.

A FSI computational framework for vascular physiopathology: A novel flow-tissue multiscale strategy

A novel fluid-structure computational framework for vascular applications is herein presented. It is developed by combining the double multi-scale nature of vascular physiopathology in terms of both tissue properties and blood flow. Addressing arterial tissues, they are modelled via a nonlinear multiscale constitutive rationale, based only on parameters having a clear histological and biochemical meaning. Moreover, blood flow is described by coupling a three-dimensional fluid domain (undergoing physiological inflow conditions) with a zero-dimensional model, which allows to reproduce the influence of the downstream vasculature, furnishing a realistic description of the outflow proximal pressure. The fluid-structure interaction is managed through an explicit time-marching approach, able to accurately describe tissue nonlinearities within each computational step for the fluid problem. A case study associated to a patient-specific aortic abdominal aneurysmatic geometry is numerically investigated, highlighting advantages gained from the proposed multiscale strategy, as well as showing soundness and effectiveness of the established framework for assessing useful clinical quantities and risk indexes.

A comparison of hemodynamic metrics and intraluminal thrombus burden in a common iliac artery aneurysm

Aneurysms of the common iliac artery (CIAA) are typically found in association with an abdominal aortic aneurysm (AAA). Isolated CIAAs, in the absence of an AAA, are uncommon. Similar to AAAs, CIAA may develop intraluminal thrombus (ILT). As isolated CIAAs have a contralateral common iliac artery for comparison, they provide an opportunity to study the hemodynamic mechanisms behind ILT formation. In this study, we compared a large isolated CIAA and the contralateral iliac artery using computational fluid dynamics to determine if hemodynamic metrics correlate with the location of ILT. We performed a comprehensive computational fluid dynamics study and investigated the residence time of platelets and monocytes, velocity fields, time-averaged wall shear stress, oscillatory shear index, and endothelial cell activation potential. We then correlated these data to ILT burden determined with computed tomography. We found that high cell residence times, low time-averaged wall shear stress, high oscillatory shear index, and high endothelial cell activation potential all correlate with regions of ILT development. Our results show agreement with previous hypotheses of thrombus formation in AAA and provide insights into the computational hemodynamics of iliac artery aneurysms.

Hemodynamics in Transplant Renal Artery Stenosis and its Alteration after Stent Implantation Based on a Patient-specific Computational Fluid Dynamics Model

Background:

Accumulating studies on computational fluid dynamics (CFD) support the involvement of hemodynamic factors in artery stenosis. Based on a patient-specific CFD model, the present study aimed to investigate the hemodynamic characteristics of transplant renal artery stenosis (TRAS) and its alteration after stent treatment.

Methods:

Computed tomography angiography (CTA) data of kidney transplant recipients in a single transplant center from April 2013 to November 2014 were reviewed. The three-dimensional geometry of transplant renal artery (TRA) was reconstructed from the qualified CTA images and categorized into three groups: the normal, stenotic, and stented groups. Hemodynamic parameters including pressure distribution, velocity, wall shear stress (WSS), and mass flow rate (MFR) were extracted. The data of hemodynamic parameters were expressed as median (interquartile range), and Mann–Whitney U-test was used for analysis.

Results:

Totally, 6 normal, 12 stenotic, and 6 stented TRAs were included in the analysis. TRAS presented nonuniform pressure distribution, adverse pressure gradient across stenosis throat, flow vortex, and a separation zone at downstream stenosis. Stenotic arteries had higher maximal velocity and maximal WSS (2.94 [2.14, 3.30] vs. 1.06 [0.89, 1.15] m/s, 256.5 [149.8, 349.4] vs. 41.7 [37.8, 45.3] Pa at end diastole, P = 0.001; 3.25 [2.67, 3.56] vs. 1.65 [1.18, 1.72] m/s, 281.3 [184.3, 364.7] vs. 65.8 [61.2, 71.9] Pa at peak systole, P = 0.001) and lower minimal WSS and MFRs (0.07 [0.03, 0.13] vs. 0.52 [0.45, 0.67] Pa, 1.5 [1.0, 3.0] vs. 11.0 [8.0, 11.3] g/s at end diastole, P = 0.001; 0.08 [0.03, 0.19] vs. 0.70 [0.60, 0.81] Pa, 2.0 [1.3, 3.3] vs. 16.5 [13.0, 20.3] g/s at peak systole, P = 0.001) as compared to normal arteries. Stent implantation ameliorated all the alterations of the above hemodynamic factors except low WSS.

Conclusions:

Hemodynamic factors were significantly changed in severe TRAS. Stent implantation can restore or ameliorate deleterious change of hemodynamic factors except low WSS at stent regions.

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.

A guide to uncertainty quantification and sensitivity analysis for cardiovascular applications

As we shift from population-based medicine towards a more precise patient-specific regime guided by predictions of verified and well-established cardiovascular models, an urgent question arises: how sensitive are the model predictions to errors and uncertainties in the model inputs? To make our models suitable for clinical decision-making, precise knowledge of prediction reliability is of paramount importance. Efficient and practical methods for uncertainty quantification (UQ) and sensitivity analysis (SA) are therefore essential. In this work, we explain the concepts of global UQ and global, variance-based SA along with two often-used methods that are applicable to any model without requiring model implementation changes: Monte Carlo (MC) and polynomial chaos (PC). Furthermore, we propose a guide for UQ and SA according to a six-step procedure and demonstrate it for two clinically relevant cardiovascular models: model-based estimation of the fractional flow reserve (FFR) and model-based estimation of the total arterial compliance (CT ). Both MC and PC produce identical results and may be used interchangeably to identify most significant model inputs with respect to uncertainty in model predictions of FFR and CT . However, PC is more cost-efficient as it requires an order of magnitude fewer model evaluations than MC. Additionally, we demonstrate that targeted reduction of uncertainty in the most significant model inputs reduces the uncertainty in the model predictions efficiently. In conclusion, this article offers a practical guide to UQ and SA to help move the clinical application of mathematical models forward. Copyright © 2015 John Wiley & Sons, Ltd.

Numerical investigation of the hydrodynamic parameters of blood flow through stenotic descending aorta

In this study, the blood flow passing through a three-dimensional geometrically realistic stenosis is investigated both experimentally and numerically. Although the blood flow in stenotic arteries has been extensively studied in the past few decades, not much work has been focused on irregularity of stenosis. Thus, a model of an irregular stenotic descending aorta is used in this work. Due to the irregularity of stenosis model, the governing differential equations for continuity and momentum are solved numerically using finite-volume/finite-difference techniques in the generalized body-fitted coordinates. In order to verify the numerical results, the experimental measured pressure drops are compared with the numerical result. In addition, an improved method for nearly orthogonal grid generation is presented in numerical study. The grid generating system is based on the solution of a set of partial differential equations with finite difference discretization. Numerical calculations are performed to examine the effect of 55% (ratio between cross-sectional area at upstream and stenosis) irregular stenosis on the hemodynamic characteristics such as flow separation zone, wall shear stress and pressure drop. The maximum calculated wall shear stress is related to the maximum velocity gradient due to minimum cross-sectional area at the neck of stenosis. In addition, the pressure is shown as an important characteristic that is effecting on the resistance against the flow in the artery. Based on our results, the 55% irregular constriction is considered critical unlike the studies that have believed the reduction, which is greater than 75% become significant.

Hemodynamic analysis in an idealized artery tree: differences in wall shear stress between Newtonian and non-Newtonian blood models

Development of many conditions and disorders, such as atherosclerosis and stroke, are dependent upon hemodynamic forces. To accurately predict and prevent these conditions and disorders hemodynamic forces must be properly mapped. Here we compare a shear-rate dependent fluid (SDF) constitutive model, based on the works by Yasuda et al in 1981, against a Newtonian model of blood. We verify our stabilized finite element numerical method with the benchmark lid-driven cavity flow problem. Numerical simulations show that the Newtonian model gives similar velocity profiles in the 2-dimensional cavity given different height and width dimensions, given the same Reynolds number. Conversely, the SDF model gave dissimilar velocity profiles, differing from the Newtonian velocity profiles by up to 25% in velocity magnitudes. This difference can affect estimation in platelet distribution within blood vessels or magnetic nanoparticle delivery. Wall shear stress (WSS) is an important quantity involved in vascular remodeling through integrin and adhesion molecule mechanotransduction. The SDF model gave a 7.3-fold greater WSS than the Newtonian model at the top of the 3-dimensional cavity. The SDF model gave a 37.7-fold greater WSS than the Newtonian model at artery walls located immediately after bifurcations in the idealized femoral artery tree. The pressure drop across arteries reveals arterial sections highly resistive to flow which correlates with stenosis formation. Numerical simulations give the pressure drop across the idealized femoral artery tree with the SDF model which is approximately 2.3-fold higher than with the Newtonian model. In atherosclerotic lesion models, the SDF model gives over 1 Pa higher WSS than the Newtonian model, a difference correlated with over twice as many adherent monocytes to endothelial cells from the Newtonian model compared to the SDF model.

Progression of abdominal aortic aneurysm towards rupture: refining clinical risk assessment using a fully coupled fluid-structure interaction method

Rupture of abdominal aortic aneurysm (AAA) is associated with high mortality rates. Risk of rupture is multi-factorial involving AAA geometric configuration, vessel tortuosity, and the presence of intraluminal pathology. Fluid structure interaction (FSI) simulations were conducted in patient based computed tomography scans reconstructed geometries in order to monitor aneurysmal disease progression from normal aortas to non-ruptured and contained ruptured AAA (rAAA), and the AAA risk of rupture was assessed. Three groups of 8 subjects each were studied: 8 normal and 16 pathological (8 non-ruptured and 8 rAAA). The AAA anatomical structures segmented included the blood lumen, intraluminal thrombus (ILT), vessel wall, and embedded calcifications. The vessel wall was described with anisotropic material model that was matched to experimental measurements of AAA tissue specimens. A statistical model for estimating the local wall strength distribution was employed to generate a map of a rupture potential index (RPI), representing the ratio between the local stress and local strength distribution. The FSI simulations followed a clear trend of increasing wall stresses from normal to pathological cases. The maximal stresses were observed in the areas where the ILT was not present, indicating a potential protective effect of the ILT. Statistically significant differences were observed between the peak systolic stress and the peak stress at the mean arterial pressure between the three groups. For the ruptured aneurysms, where the geometry of intact aneurysm was reconstructed, results of the FSI simulations clearly depicted maximum wall stress at the a priori known location of rupture. The RPI mapping indicated several distinct regions of high RPI coinciding with the actual location of rupture. The FSI methodology demonstrates that the aneurysmal disease can be described by numerical simulations, as indicated by a clear trend of increasing aortic wall stresses in the studied groups, (normal aortas, AAAs and rAAAs). Ultimately, the results demonstrate that FSI wall stress mapping and RPI can be used as a tool for predicting the potential rupture of an AAA by predicting the actual rupture location, complementing current clinical practice by offering a predictive diagnostic tool for deciding whether to intervene surgically or spare the patient from an unnecessary risky operation.

Microstructural modelling of cerebral aneurysm evolution through effective stress mediated destructive remodelling

Recently, researchers have shown an increased interest in the biomechanical modelling of cerebral aneurysm development. In the present study a fluid-solid-growth model for the formation of a fusiform aneurysm has been presented in an axi-symmetric geometry of the internal carotid artery. This model is the result of two parallel mechanisms: first, defining arterial wall as a living tissue with the ability of degradation, growth and remodelling and second, full coupling of the wall and the blood flow. Here for the first time the degradation of elastin has been defined as a function of vascular wall effective stress to take into account the shear dependent nature of degradation and the mural-cell-mediated destructive activities. The model has been stabilized in size and mechanical properties and is consistent with other computational or clinical studies. Furthermore, the evolving microstructural properties of the wall during the evolution process have been predicted.

Derivation of aortic distensibility and pulse wave velocity by image registration with a physics-based regularisation term

Analysis of the cardiovascular system represents a classical problem in which the solid and fluid phases interact intimately, and so is a rich field of application for state-of-the-art fluid-solid interaction (FSI) analyses. In this paper, we focus on the human aorta. Solution of the full FSI problem requires knowledge of the material properties of the wall and information on vessel support. We show that variation of distensibility along the aorta can be obtained from four-dimensional image data using image registration. If pressure data at one point in the vessel are available, these can be converted to absolute values. Alternatively, values of pulse wave velocity along the vessel can be obtained. The quality of the extracted data is improved by the incorporation into the registration of a regularisation term based on the one-dimensional wave equation. The method has been validated using simulated data. For idealised vessels, the accuracy with which the distensibility and wave velocity can be extracted is high (1%-2%). The method is applied to six clinical datasets from patients with mild coarctation, for which it is shown that wave velocity along the aorta is relatively constant.

Fluid-structure interaction modeling of abdominal aortic aneurysms: the impact of patient-specific inflow conditions and fluid/solid coupling

Rupture risk assessment of abdominal aortic aneurysms (AAA) by means of biomechanical analysis is a viable alternative to the traditional clinical practice of using a critical diameter for recommending elective repair. However, an accurate prediction of biomechanical parameters, such as mechanical stress, strain, and shear stress, is possible if the AAA models and boundary conditions are truly patient specific. In this work, we present a complete fluid-structure interaction (FSI) framework for patient-specific AAA passive mechanics assessment that utilizes individualized inflow and outflow boundary conditions. The purpose of the study is two-fold: (1) to develop a novel semiautomated methodology that derives velocity components from phase-contrast magnetic resonance images (PC-MRI) in the infrarenal aorta and successfully apply it as an inflow boundary condition for a patient-specific fully coupled FSI analysis and (2) to apply a one-way-coupled FSI analysis and test its efficiency compared to transient computational solid stress and fully coupled FSI analyses for the estimation of AAA biomechanical parameters. For a fully coupled FSI simulation, our results indicate that an inlet velocity profile modeled with three patient-specific velocity components and a velocity profile modeled with only the axial velocity component yield nearly identical maximum principal stress (σ1), maximum principal strain (ε1), and wall shear stress (WSS) distributions. An inlet Womersley velocity profile leads to a 5% difference in peak σ1, 3% in peak ε1, and 14% in peak WSS compared to the three-component inlet velocity profile in the fully coupled FSI analysis. The peak wall stress and strain were found to be in phase with the systolic inlet flow rate, therefore indicating the necessity to capture the patient-specific hemodynamics by means of FSI modeling. The proposed one-way-coupled FSI approach showed potential for reasonably accurate biomechanical assessment with less computational effort, leading to differences in peak σ1, ε1, and WSS of 14%, 4%, and 18%, respectively, compared to the axial component inlet velocity profile in the fully coupled FSI analysis. The transient computational solid stress approach yielded significantly higher differences in these parameters and is not recommended for accurate assessment of AAA wall passive mechanics. This work demonstrates the influence of the flow dynamics resulting from patient-specific inflow boundary conditions on AAA biomechanical assessment and describes methods to evaluate it through fully coupled and one-way-coupled fluid-structure interaction analysis.

An integrated fluid-chemical model toward modeling the formation of intra-luminal thrombus in abdominal aortic aneurysms

Abdominal Aortic Aneurysms (AAAs) are frequently characterized by the presence of an Intra-Luminal Thrombus (ILT) known to influence their evolution biochemically and biomechanically. The ILT progression mechanism is still unclear and little is known regarding the impact of the chemical species transported by blood flow on this mechanism. Chemical agonists and antagonists of platelets activation, aggregation, and adhesion and the proteins involved in the coagulation cascade (CC) may play an important role in ILT development. Starting from this assumption, the evolution of chemical species involved in the CC, their relation to coherent vortical structures (VSs) and their possible effect on ILT evolution have been studied. To this end a fluid-chemical model that simulates the CC through a series of convection-diffusion-reaction (CDR) equations has been developed. The model involves plasma-phase and surface-bound enzymes and zymogens, and includes both plasma-phase and membrane-phase reactions. Blood is modeled as a non-Newtonian incompressible fluid. VSs convect thrombin in the domain and lead to the high concentration observed in the distal portion of the AAA. This finding is in line with the clinical observations showing that the thickest ILT is usually seen in the distal AAA region. The proposed model, due to its ability to couple the fluid and chemical domains, provides an integrated mechanochemical picture that potentially could help unveil mechanisms of ILT formation and development.

Simulation model for assessing quality of ultrasound strain estimation in abdominal aortic aneurysm

The purpose of this study was to develop a simulation model for evaluating methods for ultrasound strain estimation in abdominal aortic aneurysms. Wall geometry was obtained from a real ultrasound image and wall motion was simulated applying realistic blood pressures to a nonlinear viscoelastic wall model. The ultrasound simulation included speckle, absorption and angle dependent reflection. Gaussian white noise was added to simulate various noise levels. Despite not fully replicating real ultrasound images, the model simulated realistic circumferential variations in intensity and realistic speckle patterns and has potential for initial evaluation of strain estimation methods.

Stereoscopically observed deformations of a compliant abdominal aortic aneurysm model

A new experimental setup has been implemented to precisely measure the deformations of an entire model abdominal aortic aneurysm (AAA). This setup addresses a gap between the computational and experimental models of AAA that have aimed at improving the limited understanding of aneurysm development and rupture. The experimental validation of the deformations from computational approaches has been limited by a lack of consideration of the large and varied deformations that AAAs undergo in response to physiologic flow and pressure. To address the issue of experimentally validating these calculated deformations, a stereoscopic imaging system utilizing two cameras was constructed to measure model aneurysm displacement in response to pressurization. The three model shapes, consisting of a healthy aorta, an AAA with bifurcation, and an AAA without bifurcation, were also evaluated with computational solid mechanical modeling using finite elements to assess the impact of differences between material properties and for comparison against the experimental inflations. The device demonstrated adequate accuracy (surface points were located to within 0.07 mm) for capturing local variation while allowing the full length of the aneurysm sac to be observed at once. The experimental model AAA demonstrated realistic aneurysm behavior by having cyclic strains consistent with reported clinical observations between pressures 80 and 120 mm Hg. These strains are 1-2%, and the local spatial variations in experimental strain were less than predicted by the computational models. The three different models demonstrated that the asymmetric bifurcation creates displacement differences but not cyclic strain differences within the aneurysm sac. The technique and device captured regional variations of strain that are unobservable with diameter measures alone. It also allowed the calculation of local strain and removed rigid body motion effects on the strain calculation. The results of the computations show that an asymmetric aortic bifurcation created displacement differences but not cyclic strain differences within the aneurysm sac.

Importance of initial aortic properties on the evolving regional anisotropy, stiffness and wall thickness of human abdominal aortic aneurysms

Complementary advances in medical imaging, vascular biology and biomechanics promise to enable computational modelling of abdominal aortic aneurysms to play increasingly important roles in clinical decision processes. Using a finite-element-based growth and remodelling model of evolving aneurysm geometry and material properties, we show that regional variations in material anisotropy, stiffness and wall thickness should be expected to arise naturally and thus should be included in analyses of aneurysmal enlargement or wall stress. In addition, by initiating the model from best-fit material parameters estimated for non-aneurysmal aortas from different subjects, we show that the initial state of the aorta may influence strongly the subsequent rate of enlargement, wall thickness, mechanical behaviour and thus stress in the lesion. We submit, therefore, that clinically reliable modelling of the enlargement and overall rupture-potential of aneurysms may require both a better understanding of the mechanobiological processes that govern the evolution of these lesions and new methods of determining the patient-specific state of the pre-aneurysmal aorta (or correlation to currently unaffected portions thereof) through knowledge of demographics, comorbidities, lifestyle, genetics and future non-invasive or minimally invasive tests.

Mechanics, mechanobiology, and modeling of human abdominal aorta and aneurysms

Biomechanical factors play fundamental roles in the natural history of abdominal aortic aneurysms (AAAs) and their responses to treatment. Advances during the past two decades have increased our understanding of the mechanics and biology of the human abdominal aorta and AAAs, yet there remains a pressing need for considerable new data and resulting patient-specific computational models that can better describe the current status of a lesion and better predict the evolution of lesion geometry, composition, and material properties and thereby improve interventional planning. In this paper, we briefly review data on the structure and function of the human abdominal aorta and aneurysmal wall, past models of the mechanics, and recent growth and remodeling models. We conclude by identifying open problems that we hope will motivate studies to improve our computational modeling and thus general understanding of AAAs.

A fluid-structure interaction-based numerical investigation on the evolution of stress, strength and rupture potential of an abdominal aortic aneurysm

An abdominal aortic aneurysm (AAA) is an irreversible dilation of the abdominal artery. Once an aneurysm is detected by doctors, clinical intervention is usually recommended. The interventions involve traditional open surgery repair and endovascular aneurysm repair with a stent graft. Both types of prophylactic procedures are expensive and not without any risk to the patient. It is very difficult to balance the risk of aneurysm repair and the chance of rupture. The reason lies in that the changing trend of characteristic physical quantities with the evolution of AAA and the mechanisms that give rise to it are still not completely clear. In this study, computational 3D patient-specific model for investigating AAA development was established based on computed tomography (CT) images. Results showed that as the aneurysm evolved, peak wall stress and time-averaged wall shear stress distribution patterns changed. The expansion of AAA wall resulted in the increment of peak stress. The AAA wall compliance not only showed different magnitudes at different cross-sections of the aneurismal body, but also changed with the development of the aneurysm. Furthermore, minimum wall strength and rupture potential index during the three stages of AAA evolution were also investigated in detail. This study might provide valuable information on how to further explore the mechanical basis and the rupture potential during AAA evolution, and that it may assist clinical diagnostic procedures and avoid the potential risk of unnecessary surgical intervention.