Derivation of flow related risk indices for stenosed left anterior descending coronary arteries with the use of computer simulations

Ascending Aortic Volume: A Feasible Indicator for Ascending Aortic Aneurysm Elective Surgery?

Diameter-based criterion have been widely adopted for preventive surgery of ascending thoracic aortic aneurysm (ATAA). However, recent and growing evidence has shown that diameter-based methods may not be sufficient for identifying patients who are at risk of an ATAA. In this study, fluid-structure interaction (FSI) analysis was performed on one-hundred ATAA geometries reconstructed from clinical data to examine the relationship between hemodynamic conditions, ascending aortic volume (AAV), ascending aortic curvature, and aortic ratios measured from the reconstructed 3D models. The simulated hemodynamic and biomechanical parameters were compared among different groups of ATAA geometries classified based on AAV. The ATAAs with enlarged AAV showed significantly compromised hemodynamic conditions and higher mechanical wall stress. The maximum oscillatory shear index (OSI), particle residence time (PRT) and wall stress (WS) were significantly higher in enlarged ATAAs compared with controls (0.498 [0.497, 0.499] vs 0.499 [0.498, 0.499], p = 0.002, 312.847 [207.445, 519.391] vs 996.047 [640.644, 1573.140], p < 0.001, 769.680 [668.745, 879.795] vs 1072.000 [873.060, 1280.000] kPa, p < 0.001, respectively). Values were reported as median with interquartile range (IQR). AAV was also found to be more strongly correlated with these parameters compared to maximum diameter. The correlation coefficient between AAV and average WS was as high as 0.92 (p < 0.004), suggesting that AAV might be a feasible risk identifier for ATAAs. STATEMENT OF SIGNIFICANCE: Ascending thoracic aortic aneurysm is associated with the risk of dissection or rupture, creating life-threatening conditions. Current surgical intervention guidelines are purely diameter based. Recently, many studies proposed to incorporate other morphological parameters into the current clinical guidelines to better prevent severe adverse aortic events like rupture or dissection. The purpose of this study is to gain a better understanding of the relationship between morphological parameters and hemodynamic parameters in ascending aortic aneurysms using fluid-solid-interaction analysis on patient-specific geometries. Our results suggest that ascending aortic volume may be a better indicator for surgical intervention as it shows a stronger association with pathogenic hemodynamic conditions.

Numerical simulation in the abdominal aorta and the visceral arteries with or without stenosis based on 2D PCMRI

It is important to obtain accurate boundary conditions (BCs) in hemodynamic simulations. This article aimed to improve the accuracy of BCs in computational fluid dynamics (CFD) simulation and analyze the differences in hemodynamics between healthy volunteers and patients with visceral arterial stenosis (VAS). The geometric models of seven cases were reconstructed using the magnetic resonance angiogram (MRA) or computed tomography angiogram (CTA) imaging data. The physiological flow waveforms obtained from 2D Phase Contrast Magnetic Resonance Imaging (PCMRI) were imposed on the aortic inlet and the visceral arteries' outlets. The individualized RCR values of the three-element Windkessel model were imposed on the aortic outlet. CFD simulations were run in the open-source software: svSolver. Two specific time points were selected to compare the hemodynamics of healthy volunteers and patients with VAS. The results suggested that blood in the stenotic visceral arteries flowed at high speed throughout the cardiac cycle. The low pressure is distributed at stenotic lesions. The wall shear stress (WSS) reached 4 Pa near stenotic locations. The low time average wall shear stress (TAWSS), high oscillatory shear index (OSI), and high relative residence time (RRT) concentrated in the abdominal aorta. Besides, the ratios of the areas with low TAWSS, high OSI, and high RRT to the computational domain were higher in patients with VAS than which in the healthy volunteers. The individualized BCs were used for hemodynamic simulations and results suggest that patients with stenosis have a higher risk of blood retention and atherosclerosis formation in the abdominal aorta.

Patient-Specific Computational Analysis of Hemodynamics in Adult Pulmonary Hypertension

Pulmonary hypertension (PH) is a progressive disease characterized by elevated pressure and vascular resistance in the pulmonary arteries. Nearly 250,000 hospitalizations occur annually in the US with PH as the primary or secondary condition. A definitive diagnosis of PH requires right heart catheterization (RHC) in addition to a chest computed tomography, a walking test, and others. While RHC is the gold standard for diagnosing PH, it is invasive and posseses inherent risks and contraindications. In this work, we characterized the patient-specific pulmonary hemodynamics in silico for diverse PH WHO groups. We grouped patients on the basis of mean pulmonary arterial pressure (mPAP) into three disease severity groups: at-risk ([Formula: see text], denoted with A), mild ([Formula: see text], denoted with M), and severe ([Formula: see text], denoted with S). The pulsatile flow hemodynamics was simulated by evaluating the three-dimensional Navier-Stokes system of equations using a flow solver developed by customizing OpenFOAM libraries (v5.0, The OpenFOAM Foundation). Quasi patient-specific boundary conditions were implemented using a Womersley inlet velocity profile and transient resistance outflow conditions. Hemodynamic indices such as spatially averaged wall shear stress ([Formula: see text]), wall shear stress gradient ([Formula: see text]), time-averaged wall shear stress ([Formula: see text]), oscillatory shear index ([Formula: see text]), and relative residence time ([Formula: see text]), were evaluated along with the clinical metrics pulmonary vascular resistance ([Formula: see text]), stroke volume ([Formula: see text]) and compliance ([Formula: see text]), to assess possible spatiotemporal correlations. We observed statistically significant decreases in [Formula: see text], [Formula: see text], and [Formula: see text], and increases in [Formula: see text] and [Formula: see text] with disease severity. [Formula: see text] was moderately correlated with [Formula: see text] and [Formula: see text] at the mid-notch stage of the cardiac cycle when these indices were computed using the global pulmonary arterial geometry. These results are promising in the context of a long-term goal of identifying computational biomarkers that can serve as surrogates for invasive diagnostic protocols of PH.

Coupling hemodynamics with mechanobiology in patient-specific computational models of ascending thoracic aortic aneurysms

BACKGROUND AND OBJECTIVE

The prevention of ascending thoracic aortic aneurysms (ATAAs), which affect thousands of persons every year worldwide, remains a major issue. ATAAs may be caused by anything that weakens the aortic wall. Altered hemodynamics, which concerns a majority of patients with bicuspid aortic valves, has been shown to be related to such weakening and to contribute to ATAA development and progression. However the underlying mechanisms remain unclear and computational modeling in this field could help significantly to elucidate how hemodynamics and mechanobiology interact in ATAAs.

METHODS

Accordingly, we propose a numerical framework combining computational fluid dynamics and 4D flow magnetic resonance imaging (MRI) coupled with finite element (FE) analyses to simulate growth and remodeling (G&R) occurring in patient-specific aortas in relation with altered hemodynamics. The geometries and the blood velocities obtained from 4D flow MRI are used as boundary conditions for CFD simulations. CFD simulations provide an estimation of the wall shear stress (WSS) and relative residence time (RRT) distribution across the luminal surface of the wall. An initial insult is then applied to the FE model of the aortic wall, assuming that the magnitude of the insult correlates spatially with the normalized RRT distribution obtained from CFD simulations. G&R simulations are then performed. The material behavior of each Gauss point in these FE models is evolved continuously to compensate for the deviation of the actual wall stress distribution from the homeostatic state after the initial insult. The whole approach is illustrated on two healthy and two diseased subjects. The G&R parameters are calibrated against previously established statistical models of ATAA growth rates.

RESULTS

Among the variety of results provided by G&R simulations, the analysis focused especially on the evolution of the wall stiffness, which was shown to be a major risk factor for ATAAs. It was shown that the G&R parameters, such as for instance the rate of collagen production or cell mechanosensitivity, play a critical role in ATAA progression and remodeling.

CONCLUSIONS

These preliminary findings show that patient-specific computational modeling coupling hemodynamics with mechanobiology is a promising approach to explore aneurysm progression.

A computational study of the hemodynamics of bioprosthetic aortic valves with reduced leaflet motion

We employ a reduced degree-of-freedom aortic valve model to investigate the flow physics associated with early-stage reduced leaflet motion in bioprosthetic aortic valves. The model is coupled with a sharp-interface immersed boundary based incompressible flow solver to efficiently simulate the fluid-structure interaction. A total of 19 cases of flow through aortic valves with varying degrees of reduced leaflet motion (RLM) are considered. The characteristics of the aortic jet and the consequent aorta wall loading patterns are analyzed. Our results show that asymmetric RLM tilts the aortic jet and leads to large reverse and recirculating flow regions downstream from leaflets with restricted mobility. The changes in flow patterns increase wall pressure and shear stress fluctuations, and result in asymmetric oscillating shear on the aorta wall. These findings have implications for auscultation based diagnosis of this condition as well as the health of the aorta.

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.

HEMODYNAMIC EFFECT OF OBSTRUCTION TO RENAL ARTERIES CAUSED BY STENT GRAFTS IN PATIENTS WITH ABDOMINAL AORTIC ANEURYSMS

To investigate the hemodynamic effects of partial obstruction to the renal orifice caused by inappropriate stent-graft location. Pre-operative and deployment models of a stent graft with various degrees of obstruction to the renal orifice are constructed based on medical images of abdominal aortic aneurysm. Hemodynamics, including flow pattern, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), as well as relative residence time (RRT) are analyzed using numerical simulations. Flow rate distributions are assessed and verified by in vitro experiments. Results show that partial blockage to the renal branch orifice leads to flow recirculation and vortices with low wall shear stress around the renal ostia, whereas OSI and RRT on the renal arteries increase with the degree of obstruction. While the in vitro experiment indicates a decreasing flow rate to the bilateral renal arteries as renal artery ostia are obstructed. In conclusion, obstruction to the renal arteries induced by an inappropriate stent graft location causes stenosis in the renal artery in the long term. This study reveals a possible pathological mechanism of renal complications due to the implantation of a stent graft.

Modelling of thrombin generation under flow in realistic left anterior descending geometries

Currently there are no available methods for prediction of thrombotic complications in Coronary Artery disease. Additionally, blood coagulation tests are mainly performed in a steady system while coagulation in vivo occurs under flow conditions. In this work, a phenomenological model for coagulation up-to thrombin generation is proposed; the model is mainly based on the results of thrombin generation assays and therefore it can account for the variation of the coagulability that is observed in different individuals. The model is applied on 3 cases of left anterior descending arteries (LAD) with 50% maximum stenosis placed at a different location and have been statistically assessed as of different complication risk. The simulations showed that parameters of thrombin generation assays obtain different values when they refer to thrombin generation under realistic coronary flow conditions. The flow conditions prevailing locally because of the geometric differences among the arterial trees can lead to different initiation times and thrombin production rates and it also alters the spatial distribution of the coagulation products. Similarly, small changes of the coagulation characteristics of blood under identical flow conditions can allow or prevent the initiation of coagulation. The results indicate that combined consideration of geometry and coagulation characteristics of blood can lead to entirely different conclusions compared to independent assessment of each factor.

Vascular Remodelling Relates to an Elevated Oscillatory Shear Index and Relative Residence Time in Spontaneously Hypertensive Rats

Haemodynamic disorders are common clinical findings in hypertension and lead to adverse cardiovascular events. However, the haemodynamic conditions in hypertension models are poorly understood. This study aimed to observe the characteristics of haemodynamics in spontaneously hypertensive rats (SHRs) and antihypertensive-treated SHRs. Twenty-four adult male SHRs and Wistar-Kyoto rats (WKYs) were randomly divided into four groups and treated for 7 days as follows: WKY-CON (WKYs + saline), WKY-NIF (WKYs + nifedipine, 50 mg/kg/day), SHR-CON (SHRs + saline), and SHR-NIF (SHRs + nifedipine). Aortic computational fluid dynamics (CFD) models were simulated to obtain the haemodynamic parameters. We found that in the hypertensive (SHR-CON) and blood pressure-controlled (SHR-NIF) groups, the oscillatory shear index (OSI) and relative residence time (RRT), which are key haemodynamics indices, were markedly elevated. Furthermore, there was a correlation between both the elevated OSI and RRT with the vascular wall thickening in regions near the inner wall of the aortic arch. Our research demonstrates that haemodynamics remains disturbed even if the blood pressure is normalized. In addition, vascular remodelling may play an important role in maintaining elevated OSI and RRT values.