Computational simulations of hemodynamic changes within thoracic, coronary, and cerebral arteries following early wall remodeling in response to distal aortic coarctation

Influence of boundary conditions and blood rheology on indices of wall shear stress from IVUS-based patient-specific stented coronary artery simulations

The long-term clinical efficacy of coronary stents is limited by restenosis. Coronary stenting results in altered arterial geometry, local blood flow patterns, and wall shear stress (WSS), all of which can influence restenosis. Computational fluid dynamics (CFD) simulations employ assumptions about blood properties and boundary conditions, which also influence WSS alterations from stenting. Our objective was to evaluate three common assumptions applied with stented coronary artery CFD simulations (inlet velocity profile, outlet boundary conditions, and viscosity) to provide insight for future studies. A patient-specific right coronary artery was reconstructed from intravascular ultrasound and computed tomography imaging. Time-averaged WSS (TAWSS) and oscillatory shear index (OSI) were compared for CFD simulations using parabolic and Womersley (α = 2.5) velocity profiles and Newtonian versus non-Newtonian (Carreau-Yasuda) viscosity. TAWSS and OSI differences were quantified for 5-element lumped parameter network (LPN) outlet boundary conditions as compared to a 3-element Windkessel approach previously applied by neglecting ventricular contraction. Differences in velocity profiles were negligible beyond two diameters from the inlet. Differences between Womersley inlet 3-element and 5-element LPNs were ~ 5% for TAWSS and ~ 200% for OSI, and most pronounced near stent struts. OSI differences were due to pressure differences between outlet boundary conditions inducing different near-wall velocity gradients. TAWSS differences between viscosity models were greatest near struts (~ 100%). Collectively these results suggest certain assumptions commonly applied for simulations of stented coronary arteries have a greater impact on TAWSS and OSI than others, with outlet boundary conditions being paramount, followed by viscosity model and inlet velocity profile.

Reliability of characterising coronary artery flow with the flow-split outflow strategy: Comparison against the multiscale approach

BACKGROUND

In computational modelling of coronary haemodynamics, imposing patient-specific flow conditions is paramount, yet often impractical due to resource and time constraints, limiting the ability to perform a large number of simulations particularly for diseased cases.

OBJECTIVE

To compare coronary haemodynamics quantified using a simplified flow-split strategy with varying exponents against the clinically verified but computationally intensive multiscale simulations under both resting and hyperaemic conditions in arteries with varying degrees of stenosis.

METHODS

Six patient-specific left coronary artery trees were segmented and reconstructed, including three with severe (>70 %) and three with mild (<50 %) focal stenoses. Simulations were performed for the entire coronary tree to account for the flow-limiting effects from epicardial artery stenoses. Both a 0D-3D coupled multiscale model and a flow-split approach with four different exponents (2.0, 2.27, 2.33, and 3.0) were used. The resulting prominent haemodynamic metrics were statistically compared between the two methods.

RESULTS

Flow-split and multiscale simulations did not significantly differ under resting conditions regardless of the stenosis severity. However, under hyperaemic conditions, the flow-split method significantly overestimated the time-averaged wall shear stress by up to 16.8 Pa (p = 0.031) and underestimate the fractional flow reserve by 0.327 (p = 0.043), with larger discrepancies observed in severe stenoses than in mild ones. Varying the exponent from 2.0 to 3.0 within the flow-split methods did not significantly affect the haemodynamic results (p > 0.141).

CONCLUSIONS

Flow-split strategies with exponents between 2.0 and 3.0 are appropriate for modelling stenosed coronaries under resting conditions. Multiscale simulations are recommended for accurate modelling of hyperaemic conditions, especially in severely stenosed arteries.(247/250 words).

Hemodynamics of Proximal Coronary Lesions in Patients Undergoing Transcatheter Aortic Valve Implantation: Patient-Specific In Silico Study

Aortic stenosis (AS) frequently coexists with coronary artery disease (CAD), complicating revascularization decisions. The use of coronary physiology indices, such as the fractional flow reserve (FFR), instantaneous wave-free ratio (iFR), and coronary flow reserve (CFR), in AS patients remains debated, particularly after transcatheter aortic valve implantation (TAVI). In this study, we employ computational fluid dynamics (CFD) to evaluate coronary hemodynamics and assess changes in the wall shear stress (WSS) before and after TAVI. Our analysis demonstrates strong agreement between CFD-derived and invasive FFR measurements, confirming CFD's reliability as a non-invasive tool for coronary physiology assessment. Furthermore, our results show no significant changes in FFR (p=0.92), iFR (p=0.67), or CFR (p=0.34) post-TAVI, suggesting that these indices remain stable following aortic valve intervention. However, a significant reduction in high WSS exposure (59% to 40.8%, p<0.001) and the oscillatory shear index (OSI: 0.32 to 0.21, p<0.001) was observed, indicating improved hemodynamic stability. These findings suggest that coronary physiology indices remain reliable for revascularization guidance post-TAVI and highlight a potential beneficial effect of aortic stenosis treatment on plaque shear stress dynamics. Our study underscores the clinical utility of CFD modeling in CAD management, paving the way for further research into its prognostic implications.

Mechanical, structural, and morphological differences in the iliac arteries

Iliac arteries play a crucial role in peripheral blood circulation. They are susceptible to various diseases, including aneurysms and atherosclerosis. Structure, material properties, and biomechanical forces acting on different regions of the iliac vasculature may contribute to the localization and progression of these pathologies. We examined 33 arterial specimens from common iliac (CI), external iliac (EI), and internal iliac (II) arteries obtained from 11 human donors (62 ± 12 years). We conducted morphometric, mechanical, and structural analyses using planar biaxial tests, constitutive modeling, and bi-directional histology on transverse and axial sections. The iliac arteries exhibited increased tortuosity and varying disease distribution with age. CI and II arteries displayed non-uniform age-related disease progression around their circumference, while EI remained healthy even in older individuals. Trends in load-free and stress-free thickness varied along the iliac vasculature. Longitudinally, EI exhibited the highest compliance compared to other iliac vessels. In contrast, CI was stiffest longitudinally, and EI was the stiffest circumferentially. Material parameters for all iliac vessels are reported for four common constitutive relations. Elastin near the internal elastic lamina displayed greater waviness in EI and II compared to CI. Also, EI had the least glycosaminoglycans (GAGs) and the highest elastin content. Our findings highlight variations in the morphological, mechanical, and structural properties of iliac arteries along their length. This data can inform vascular disease development and computational studies, and guide the development of biomimetic repair materials and devices tailored to specific iliac locations, improving vascular repair strategies.

Establishing the distribution of cerebrovascular resistance using computational fluid dynamics and 4D flow MRI

Cerebrovascular resistance (CVR) regulates blood flow in the brain, but little is known about the vascular resistances of the individual cerebral territories. We present a method to calculate these resistances and investigate how CVR varies in the hemodynamically disturbed brain. We included 48 patients with stroke/TIA (29 with symptomatic carotid stenosis). By combining flow rate (4D flow MRI) and structural computed tomography angiography (CTA) data with computational fluid dynamics (CFD) we computed the perfusion pressures out from the circle of Willis, with which CVR of the MCA, ACA, and PCA territories was estimated. 56 controls were included for comparison of total CVR (tCVR). CVR were 33.8 ± 10.5, 59.0 ± 30.6, and 77.8 ± 21.3 mmHg s/ml for the MCA, ACA, and PCA territories. We found no differences in tCVR between patients, 9.3 ± 1.9 mmHg s/ml, and controls, 9.3 ± 2.0 mmHg s/ml (p = 0.88), nor in territorial CVR in the carotid stenosis patients between ipsilateral and contralateral hemispheres. Territorial resistance associated inversely to territorial brain volume (p < 0.001). These resistances may work as reference values when modelling blood flow in the circle of Willis, and the method can be used when there is need for subject-specific analysis.

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary anatomies combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. We propose an optimization-based method to personalize multiscale coronary flow simulations by assimilating clinical CT myocardial perfusion imaging and cardiac function measurements to yield patient-specific flow distributions and model parameters. Using this proof-of-concept study on a cohort of six patients, we reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based purely on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

Using Gaussian process for velocity reconstruction after coronary stenosis applicable in positron emission particle tracking: An in-silico study

Accurate velocity reconstruction is essential for assessing coronary artery disease. We propose a Gaussian process method to reconstruct the velocity profile using the sparse data of the positron emission particle tracking (PEPT) in a biological environment, which allows the measurement of tracer particle velocity to infer fluid velocity fields. We investigated the influence of tracer particle quantity and detection time interval on flow reconstruction accuracy. Three models were used to represent different levels of stenosis and anatomical complexity: a narrowed straight tube, an idealized coronary bifurcation with stenosis, and patient-specific coronary arteries with a stenotic left circumflex artery. Computational fluid dynamics (CFD), particle tracking, and the Gaussian process of kriging were employed to simulate and reconstruct the pulsatile flow field. The study examined the error and uncertainty in velocity profile reconstruction after stenosis by comparing particle-derived flow velocity with the CFD solution. Using 600 particles (15 batches of 40 particles) released in the main coronary artery, the time-averaged error in velocity reconstruction ranged from 13.4% (no occlusion) to 161% (70% occlusion) in patient-specific anatomy. The error in maximum cross-sectional velocity at peak flow was consistently below 10% in all cases. PEPT and kriging tended to overestimate area-averaged velocity in higher occlusion cases but accurately predicted maximum cross-sectional velocity, particularly at peak flow. Kriging was shown to be useful to estimate the maximum velocity after the stenosis in the absence of negative near-wall velocity.

Predictors of Myocardial Ischemia in Patients with Kawasaki Disease: Insights from Patient-Specific Simulations of Coronary Hemodynamics

Current treatments for patients with coronary aneurysms caused by Kawasaki disease (KD) are based primarily on aneurysm size. This ignores hemodynamic factors influencing myocardial ischemic risk. We performed patient-specific computational hemodynamics simulations for 15 KD patients, with parameters tuned to patients' arterial pressure and cardiac function. Ischemic risk was evaluated in 153 coronary arteries from simulated fractional flow reserve (FFR), wall shear stress, and residence time. FFR correlated weakly with aneurysm [Formula: see text]-scores (correlation coefficient, [Formula: see text]) but correlated better with the ratio of maximum-to-minimum aneurysmal lumen diameter ([Formula: see text]). FFR dropped more rapidly distal to aneurysms, and this correlated more with the lumen diameter ratio ([Formula: see text]) than [Formula: see text]-score ([Formula: see text]). Wall shear stress correlated better with the diameter ratio ([Formula: see text]), while residence time correlated more with [Formula: see text]-score ([Formula: see text]). Overall, the maximum-to-minimum diameter ratio predicted ischemic risk better than [Formula: see text]-score. Although FFR immediately distal to aneurysms was nonsignificant, its rapid rate of decrease suggests elevated risk.

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary models combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. Personalized flow distributions and model parameters are informed by clinical CT myocardial perfusion imaging and cardiac function using surrogate-based optimization. We reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

Impact of TAVR on coronary artery hemodynamics using clinical measurements and image-based patient-specific in silico modeling

In recent years, transcatheter aortic valve replacement (TAVR) has become the leading method for treating aortic stenosis. While the procedure has improved dramatically in the past decade, there are still uncertainties about the impact of TAVR on coronary blood flow. Recent research has indicated that negative coronary events after TAVR may be partially driven by impaired coronary blood flow dynamics. Furthermore, the current technologies to rapidly obtain non-invasive coronary blood flow data are relatively limited. Herein, we present a lumped parameter computational model to simulate coronary blood flow in the main arteries as well as a series of cardiovascular hemodynamic metrics. The model was designed to only use a few inputs parameters from echocardiography, computed tomography and a sphygmomanometer. The novel computational model was then validated and applied to 19 patients undergoing TAVR to examine the impact of the procedure on coronary blood flow in the left anterior descending (LAD) artery, left circumflex (LCX) artery and right coronary artery (RCA) and various global hemodynamics metrics. Based on our findings, the changes in coronary blood flow after TAVR varied and were subject specific (37% had increased flow in all three coronary arteries, 32% had decreased flow in all coronary arteries, and 31% had both increased and decreased flow in different coronary arteries). Additionally, valvular pressure gradient, left ventricle (LV) workload and maximum LV pressure decreased by 61.5%, 4.5% and 13.0% respectively, while mean arterial pressure and cardiac output increased by 6.9% and 9.9% after TAVR. By applying this proof-of-concept computational model, a series of hemodynamic metrics were generated non-invasively which can help to better understand the individual relationships between TAVR and mean and peak coronary flow rates. In the future, tools such as these may play a vital role by providing clinicians with rapid insight into various cardiac and coronary metrics, rendering the planning for TAVR and other cardiovascular procedures more personalized.

Recasting Current Knowledge of Human Fetal Circulation: The Importance of Computational Models

Computational hemodynamic simulations are becoming increasingly important for cardiovascular research and clinical practice, yet incorporating numerical simulations of human fetal circulation is relatively underutilized and underdeveloped. The fetus possesses unique vascular shunts to appropriately distribute oxygen and nutrients acquired from the placenta, adding complexity and adaptability to blood flow patterns within the fetal vascular network. Perturbations to fetal circulation compromise fetal growth and trigger the abnormal cardiovascular remodeling that underlies congenital heart defects. Computational modeling can be used to elucidate complex blood flow patterns in the fetal circulatory system for normal versus abnormal development. We present an overview of fetal cardiovascular physiology and its evolution from being investigated with invasive experiments and primitive imaging techniques to advanced imaging (4D MRI and ultrasound) and computational modeling. We introduce the theoretical backgrounds of both lumped-parameter networks and three-dimensional computational fluid dynamic simulations of the cardiovascular system. We subsequently summarize existing modeling studies of human fetal circulation along with their limitations and challenges. Finally, we highlight opportunities for improved fetal circulation models.

Evaluation of Different Cannulation Strategies for Aortic Arch Surgery Using a Cardiovascular Numerical Simulator

Aortic disease has a significant impact on quality of life. The involvement of the aortic arch requires the preservation of blood supply to the brain during surgery. Deep hypothermic circulatory arrest is an established technique for this purpose, although neurological injury remains high. Additional techniques have been used to reduce risk, although controversy still remains. A three-way cannulation approach, including both carotid arteries and the femoral artery or the ascending aorta, has been used successfully for aortic arch replacement and redo procedures. We developed circuits of the circulation to simulate blood flow during this type of cannulation set up. The CARDIOSIM© cardiovascular simulation platform was used to analyse the effect on haemodynamic and energetic parameters and the benefit derived in terms of organ perfusion pressure and flow. Our simulation approach based on lumped-parameter modelling, pressure-volume analysis and modified time-varying elastance provides a theoretical background to a three-way cannulation strategy for aortic arch surgery with correlation to the observed clinical practice.

Long-term prognostic impact of paravalvular leakage on coronary artery disease requires patient-specific quantification of hemodynamics

Transcatheter aortic valve replacement (TAVR) is a frequently used minimally invasive intervention for patient with aortic stenosis across a broad risk spectrum. While coronary artery disease (CAD) is present in approximately half of TAVR candidates, correlation of post-TAVR complications such as paravalvular leakage (PVL) or misalignment with CAD are not fully understood. For this purpose, we developed a multiscale computational framework based on a patient-specific lumped-parameter algorithm and a 3-D strongly-coupled fluid-structure interaction model to quantify metrics of global circulatory function, metrics of global cardiac function and local cardiac fluid dynamics in 6 patients. Based on our findings, PVL limits the benefits of TAVR and restricts coronary perfusion due to the lack of sufficient coronary blood flow during diastole phase (e.g., maximum coronary flow rate reduced by 21.73%, 21.43% and 21.43% in the left anterior descending (LAD), left circumflex (LCX) and right coronary artery (RCA) respectively (N = 6)). Moreover, PVL may increase the LV load (e.g., LV load increased by 17.57% (N = 6)) and decrease the coronary wall shear stress (e.g., maximum wall shear stress reduced by 20.62%, 21.92%, 22.28% and 25.66% in the left main coronary artery (LMCA), left anterior descending (LAD), left circumflex (LCX) and right coronary artery (RCA) respectively (N = 6)), which could promote atherosclerosis development through loss of the physiological flow-oriented alignment of endothelial cells. This study demonstrated that a rigorously developed personalized image-based computational framework can provide vital insights into underlying mechanics of TAVR and CAD interactions and assist in treatment planning and patient risk stratification in patients.

STUDY ON THE NONLINEAR DEFORMATION MECHANICAL PROPERTIES OF A NEW TYPE OF VENA CAVA FILTER

In order to analyze the mechanical properties of a new type of retrievable vena cava filter, the Denali and Aeginsy retrievable filters are selected, and the in vivo and in vitro mechanical properties of the three types of retrievable vena cava filters are compared. The nonlinear deformation process of the filter which is pressed into the transport sheath, flicked and released into the vascular vessel is simulated by using the finite element method, and the mechanical behaviors of the filter and the blood vessel wall are analyzed. During the crimping process, the elastic strain peak value of the filter wall increases with increasing amount of grip. The maximum elastic strain peak value of the three filters is less than 8%, and the safety factor is greater than 2, which meets the strength requirement. Equivalent stress and strain on the vascular wall decreases along the direction of blood flow after the filter is opened. The peak value of equivalent stress and strain on the vascular wall corresponding to Aegisy filter is the largest, and the impact on the vascular wall is the largest. The total displacement of the three filters under the action of thrombus and blood fluids is not more than 80[Formula: see text][Formula: see text]m. The impact of the new filter on the vascular wall is smaller than that of Aegisy filter, and the total displacement after implantation is smaller than that of Denali filter. The results are valuable for the evaluation and clinical application of the new vena cava filter.

Computational Assessment of Hemodynamic Significance in Patients With Intramural Anomalous Aortic Origin of the Coronary Artery Using Virtually Derived Fractional Flow Reserve and Downstream Microvascular Resistance

Anomalous aortic origin of a coronary artery (AAOCA) is the second most common cause of sudden cardiac death in young athletes. One of the hypothesized mechanisms of ischemia in these patients is the lateral compression of the anomalous artery with an intramural or interarterial course. The presence of a narrowing in the anomalous artery will cause physiologic changes in downstream resistance that should be included for computational assessment of possible clinical ramifications. In this study, we created different compression levels, i.e., proximal narrowing, in the intramural course of a representative patient model and calculated hyperemic stenosis resistance (HSR) as well as virtual fractional flow reserve (vFFR). Models also included the effect of the distal hyperemic microvascular resistance (HMR) on vFFR. Our results agreed with similar FFR studies indicating that FFR is increased with increasing HMR and that different compression levels could have similar FFR depending on the HMR. For example, vFFR at HSR: 1.0-1.3 and HMR: 2.30 mmHg/cm/s is 0.68 and close to vFFR at HSR: 0.6-0.7 and HMR: 1.6 mmHg/cm/s, which is 0.7. The current findings suggest that functional assessment of anomalous coronary arteries through FFR should consider the vascular resistance distal to the narrowing in addition to the impact of a proximal narrowing and provides computational approaches for implementation of these important considerations.

Personalized Pre- and Post-Operative Hemodynamic Assessment of Aortic Coarctation from 3D Rotational Angiography

PURPOSE

Coarctation of Aorta (CoA) is a congenital disease consisting of a narrowing that obstructs the systemic blood flow. This proof-of-concept study aimed to develop a framework for automatically and robustly personalizing aortic hemodynamic computations for the assessment of pre- and post-intervention CoA patients from 3D rotational angiography (3DRA) data.

METHODS

We propose a framework that combines hemodynamic modelling and machine learning (ML) based techniques, and rely on 3DRA data for non-invasive pressure computation in CoA patients. The key features of our framework are a parameter estimation method for calibrating inlet and outlet boundary conditions, and regional mechanical wall properties, to ensure that the computational results match the patient-specific measurements, and an improved ML based pressure drop model capable of predicting the instantaneous pressure drop for a wide range of flow conditions and anatomical CoA variations.

RESULTS

We evaluated the framework by investigating 6 patient datasets, under pre- and post-operative setting, and, since all calibration procedures converged successfully, the proposed approach is deemed robust. We compared the peak-to-peak and the cycle-averaged pressure drop computed using the reduced-order hemodynamic model with the catheter based measurements, before and after virtual and actual stenting. The mean absolute error for the peak-to-peak pressure drop, which is the most relevant measure for clinical decision making, was 2.98 mmHg for the pre- and 2.11 mmHg for the post-operative setting. Moreover, the proposed method is computationally efficient: the average execution time was of only [Formula: see text] minutes on a standard hardware configuration.

CONCLUSION

The use of 3DRA for hemodynamic modelling could allow for a complete hemodynamic assessment, as well as virtual interventions or surgeries and predictive modeling. However, before such an approach can be used routinely, significant advancements are required for automating the workflow.

The effect of coarctation degrees on wall shear stress indices

Coarctation of the aorta (CoA) is a congenital tightening of the proximal descending aorta. Flow quantification can be immensely valuable for an early and accurate diagnosis. However, there is a lack of appropriate diagnostic approaches for a variety of cardiovascular diseases, such as CoA. An accurate understanding of the disease depends on measurements of the global haemodynamics (criteria for heart function) and also the local haemodynamics (detailed data on the dynamics of blood flow). Playing a significant role in clinical processes, wall shear stress (WSS) cannot be measured clinically; thus, computation tools are needed to give an insight into this crucial haemodynamic parameter. In the present study, in order to enable the progress of non-invasive approaches that quantify global and local haemodynamics for different CoA severities, innovative computational blueprint simulations that include fluid-solid interaction models are developed. Since there is no clear approach for managing the CoA regarding its severity, this study proposes the use of WSS indices and pressure gradient to better establish a framework for treatment procedures in CoA patients with different severities. This provides a platform for improving CoA therapy on a patient-specific level, in which physicians can perform treatment methods based on WSS indices on top of using a mere experience. Results show how severe CoA affects the aorta in comparison to the milder cases, which can give the medical community valuable information before and after any intervention.

CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation

In this work, we describe the CRIMSON (CardiovasculaR Integrated Modelling and SimulatiON) software environment. CRIMSON provides a powerful, customizable and user-friendly system for performing three-dimensional and reduced-order computational haemodynamics studies via a pipeline which involves: 1) segmenting vascular structures from medical images; 2) constructing analytic arterial and venous geometric models; 3) performing finite element mesh generation; 4) designing, and 5) applying boundary conditions; 6) running incompressible Navier-Stokes simulations of blood flow with fluid-structure interaction capabilities; and 7) post-processing and visualizing the results, including velocity, pressure and wall shear stress fields. A key aim of CRIMSON is to create a software environment that makes powerful computational haemodynamics tools accessible to a wide audience, including clinicians and students, both within our research laboratories and throughout the community. The overall philosophy is to leverage best-in-class open source standards for medical image processing, parallel flow computation, geometric solid modelling, data assimilation, and mesh generation. It is actively used by researchers in Europe, North and South America, Asia, and Australia. It has been applied to numerous clinical problems; we illustrate applications of CRIMSON to real-world problems using examples ranging from pre-operative surgical planning to medical device design optimization.

[4D FLOW Magnetic Resonance Imaging in the Study of Blood Flow in Patients With Aortic Coarctation in the Long-Term After Surgery]

Aim      Comprehensive evaluation of blood flow in the thoracic aorta using a software for 4D processing of magnetic resonance (MR) images of the heart and blood vessels (4D Flow) in patients with aortic coarctation in the late postoperative period.Materials and methods The MR study of the heart was performed for 10 patients (7 boys and 3 girls) aged 8 to 13 years (median, 9.5 [8.3; 10.8] years) who underwent resection with end-to-end anastomosis for aortic coarctation at age of 2 weeks to 10 months. MR tomography was performed on a 1.5 T MR scanner using a multichannel surface coil for scanning, electrocardiographic synchronization, and a specialized package of pulse sequences for scanning of the heart. Blood flow was evaluated with a 4D data handling software for processing of MR images of heart and blood vessels (4D Flow). The following blood flow parameters were analyzed: blood flow volume per second, peak blood flow velocity, peak and minimum blood flow area at the levels of ascending aorta, arch, isthmus, and descending aorta, and pressure gradient at the level of maximum narrowing of the aorta. 3D-MR images were used for evaluation of aortic geometry. Blood flow formation, distribution, and trajectories were analyzed by maps of vectors, particle trace, and stream lines. Statistical analysis was performed with a Statistica (v. 6.0 StatSoft Inc.) package.Results Accelerated flow in the region of residual aortic stenosis in systole was observed in all patients; 4 patients had an additional vortex flow below the aortic stenosis and a spiral flow in the descending aorta. The pressure gradient on the aortic isthmus was directly correlated with the left ventricular myocardial mass index (r=0.65; р=0.04) and indexes of blood flow in the ascending and descending aorta (р=0.03; р=0.026). No significant correlation was found for blood flow indexes and geometry of the aortic arch (H / L). Delayed contrast enhancement MR imaging did not detect any fibrotic changes in the myocardium in only one patient. The fibrosis severity inversely correlated with the right ventricular ejection fraction (r=0.65; р=0.04) and directly correlated with the pressure gradient at the aortic isthmus (r=0.63; p=0.05).Conclusion      The 4D MR image processing software for the heart and blood vessels allows studying the blood flow in detail under natural conditions, provides potential advantages in comprehensive evaluation of patients with aortic coarctation during a dynamic follow-up. For a definitive conclusion about the relationship between the altered blood flow in the thoracic aorta and markers of residual, post-correction pathology, larger studies are required as well as long-term follow-up of patients with documented pathological patterns of blood flow (changes in blood flow velocity and volume throughout the entire thoracic aorta in combination with disorders in the normal flow geometry during the cardiac cycle).

Patient-Specific Numerical Analysis of Coronary Flow in Children With Intramural Anomalous Aortic Origin of Coronary Arteries

Unroofing surgery for anomalous aortic origin of a coronary artery (AAOCA) alters coronary anatomy by opening the intramural segment so that the anomalous coronary orifice arises perpendicularly from appropriate aortic sinus. Computational fluid dynamics modeling (CFD) allows for quantification of hemodynamics linked to morbidity such as wall shear stress (WSS), relative to patient-specific features like the angle of origin (AO). We hypothesize that CFD will reveal abnormal WSS indices in unroofed arteries that are related to AO. Six AAOCA patients (3 left, 3 right) status post unroofing (median = 13.5 years, range 9-17) underwent cardiac magnetic resonance imaging. CFD models were created from pre (n = 2) and postunroofing (n = 6) cardiac magnetic resonance imaging data, for the anomalous and contralateral normally-arising arteries. Downstream vasculature was represented by lumped parameter networks. Time-averaged WSS (TAWSS) and oscillatory shear index (OSI) were quantified relative to AO and measured hemodynamics. TAWSS was elevated along the outer wall of the normally-arising left vs right coronary arteries, as well as along unroofed left vs right coronary arteries (n = 6/group). No significant differences were noted when comparing unroofed and same-sided normally-arising coronaries. TAWSS was reduced after unroofing (eg, 276 ± 28 dyne/cm² vs 91 ± 15 dyne/cm²; n = 2/group). Models with more acute preoperative AO indicated lower TAWSS at the proximity of ostium. Differences in OSI were not significant. Different flow patterns exist natively between right and left coronary arteries. Unroofing may normalize TAWSS but with variance related to the AO. This study suggests CFD may help stratify risk in AAOCA.

The effects of clinically-derived parametric data uncertainty in patient-specific coronary simulations with deformable walls

Cardiovascular simulations are increasingly used for noninvasive diagnosis of cardiovascular disease, to guide treatment decisions, and in the design of medical devices. Quantitative assessment of the variability of simulation outputs due to input uncertainty is a key step toward further integration of cardiovascular simulations in the clinical workflow. In this study, we present uncertainty quantification in computational models of the coronary circulation to investigate the effect of uncertain parameters, including coronary pressure waveform, intramyocardial pressure, morphometry exponent, and the vascular wall Young's modulus. We employ a left coronary artery model with deformable vessel walls, simulated via an Arbitrary-Lagrangian-Eulerian framework for fluid-structure interaction, with a prescribed inlet pressure and open-loop lumped parameter network outlet boundary conditions. Stochastic modeling of the uncertain inputs is determined from intra-coronary catheterization data or gathered from the literature. Uncertainty propagation is performed using several approaches including Monte Carlo, Quasi Monte Carlo sampling, stochastic collocation, and multi-wavelet stochastic expansion. Variabilities in the quantities of interest, including branch pressure, flow, wall shear stress, and wall deformation are assessed. We find that uncertainty in inlet pressures and intramyocardial pressures significantly affect all resulting QoIs, while uncertainty in elastic modulus only affects the mechanical response of the vascular wall. Variability in the morphometry exponent used to distribute the total downstream vascular resistance to the single outlets, has little effect on coronary hemodynamics or wall mechanics. Finally, we compare convergence behaviors of statistics of QoIs using several uncertainty propagation methods on three model benchmark problems and the left coronary simulations. From the simulation results, we conclude that the multi-wavelet stochastic expansion shows superior accuracy and performance against Quasi Monte Carlo and stochastic collocation methods.

Hemodynamics and tissue biomechanics of the thoracic aorta with a trileaflet aortic valve at different phases of valve opening

In a normal cardiac cycle, the trileaflet aortic valve opening is progressive, which correlates with the phasic blood flow. Therefore, we aimed to determine the impact of including an anatomically accurate reconstructed trileaflet aortic valve within a fluid-structure interaction (FSI) simulation model and determine the cyclical hemodynamic forces imposed on the thoracic aortic walls from aortic valve opening to closure. A pediatric patient with a normal trileaflet valve was recruited. Using the Cardiac Magnetic Resonance Data (CMR), a 3D model of the aortic valve and thoracic aorta was reconstructed. FSI simulations were employed to assess the tissue stress during a cardiac cycle as the result of changes in the valve opening. The blood flow was simulated as a mixture of blood plasma and red blood cells to account for non-Newtonian effects. The computation was validated with phase-contrast CMR. Windkessel boundary conditions were employed to ensure physiological pressures during the cardiac cycle. The leaflets' dynamic motion during the cardiac cycle was defined with an analytic grid velocity function. At the beginning of the valve opening a thin jet is developing. From mid-open towards full opening the stress level increases where the jet impinges the convex wall. At peak systole two counter-rotating Dean-like vortex cores manifest in the ascending aorta, which correlates with increased integrated mean stress levels. An accurate trileaflet aortic valve is needed for capturing of both primary and secondary flow features that impact the forces on the thoracic aorta wall. Omitting the aortic valve underestimates the biomechanical response.

Multilevel and multifidelity uncertainty quantification for cardiovascular hemodynamics

Standard approaches for uncertainty quantification in cardiovascular modeling pose challenges due to the large number of uncertain inputs and the significant computational cost of realistic three-dimensional simulations. We propose an efficient uncertainty quantification framework utilizing a multilevel multifidelity Monte Carlo (MLMF) estimator to improve the accuracy of hemodynamic quantities of interest while maintaining reasonable computational cost. This is achieved by leveraging three cardiovascular model fidelities, each with varying spatial resolution to rigorously quantify the variability in hemodynamic outputs. We employ two low-fidelity models (zero- and one-dimensional) to construct several different estimators. Our goal is to investigate and compare the efficiency of estimators built from combinations of these two low-fidelity model alternatives and our high-fidelity three-dimensional models. We demonstrate this framework on healthy and diseased models of aortic and coronary anatomy, including uncertainties in material property and boundary condition parameters. Our goal is to demonstrate that for this application it is possible to accelerate the convergence of the estimators by utilizing a MLMF paradigm. Therefore, we compare our approach to single fidelity Monte Carlo estimators and to a multilevel Monte Carlo approach based only on three-dimensional simulations, but leveraging multiple spatial resolutions. We demonstrate significant, on the order of 10 to 100 times, reduction in total computational cost with the MLMF estimators. We also examine the differing properties of the MLMF estimators in healthy versus diseased models, as well as global versus local quantities of interest. As expected, global quantities such as outlet pressure and flow show larger reductions than local quantities, such as those relating to wall shear stress, as the latter rely more heavily on the highest fidelity model evaluations. Similarly, healthy models show larger reductions than diseased models. In all cases, our workflow coupling Dakota's MLMF estimators with the SimVascular cardiovascular modeling framework makes uncertainty quantification feasible for constrained computational budgets.

Towards non-invasive computational-mechanics and imaging-based diagnostic framework for personalized cardiology for coarctation

Coarctation of the aorta (COA) is a congenital narrowing of the proximal descending aorta. Although accurate and early diagnosis of COA hinges on blood flow quantification, proper diagnostic methods for COA are still lacking because fluid-dynamics methods that can be used for accurate flow quantification are not well developed yet. Most importantly, COA and the heart interact with each other and because the heart resides in a complex vascular network that imposes boundary conditions on its function, accurate diagnosis relies on quantifications of the global hemodynamics (heart-function metrics) as well as the local hemodynamics (detailed information of the blood flow dynamics in COA). In this study, to enable the development of new non-invasive methods that can quantify local and global hemodynamics for COA diagnosis, we developed an innovative fast computational-mechanics and imaging-based framework that uses Lattice Boltzmann method and lumped-parameter modeling that only need routine non-invasive clinical patient data. We used clinical data of patients with COA to validate the proposed framework and to demonstrate its abilities to provide new diagnostic analyses not possible with conventional diagnostic methods. We validated this framework against clinical cardiac catheterization data, calculations using the conventional finite-volume method and clinical Doppler echocardiographic measurements. The diagnostic information, that the framework can provide, is vitally needed to improve clinical outcomes, to assess patient risk and to plan treatment.

A protocol for automated a posteriori adaptive meshing with SimVascular: a test case

OBJECTIVE

Operational details regarding the use of the adaptive meshing (AM) algorithm available in the SimVascular package are scarce despite its application in several studies. Lacking these details, novice users of the AM algorithm may experience undesirable outcomes post-adaptation such as increases in mesh error metrics, unpredictable increases in mesh size, and losses in geometric fidelity. Here we present a test case using our proposed iterative protocol that will help prevent these undesirable outcomes and enhance the utility of the AM algorithm. We present three trials (conservative, moderate, and aggressive settings) applied to a scenario modelling a Fontan junction with a patient-specific geometry and physiologically realistic boundary conditions.

RESULTS

In all three trials, an overall reduction in mesh error metrics is observed (range 47%-86%). The increase in the number of elements through each adaptation never exceeded the mesh size of the pre-adaptation mesh by one order of magnitude. In all three trials, the protocol resulted in consistent, repeatable improvements in mesh error metrics, no losses of geometric fidelity and steady increments in the number of elements in the mesh. Our proposed protocol prevented the aforementioned undesirable outcomes and can potentially save new users considerable effort and computing resources.

Effects of Normal Variation in the Rotational Position of the Aortic Root on Hemodynamics and Tissue Biomechanics of the Thoracic Aorta

PURPOSE

Variation in the rotational position of the aortic root relative to the left ventricle is present in normal trileaflet aortic valves. Its impact on the resulting fluid mechanics of blood flow in the thoracic aorta and structural mechanics in the aortic wall are unknown. We aimed to determine the regional hemodynamic and biomechanical differences in different rotational positions of the normal aortic root (clockwise, central, and counterclockwise positions).

METHOD

Cardiac magnetic resonance imaging (CMR) data was acquired from a normal pediatric patient. These were used for reconstruction of the aortic valve and thoracic aorta 3D model. Fluid-structure interaction (FSI) simulations were employed to study the influence of the root rotation with a central position as compared to observed extreme variations. Patient-specific phase-encoding CMR data were used to assess the validity of computed blood flow. The 3D FSI model was coupled with Windkessel boundary conditions that were tuned for physiological pressures. A grid velocity function was adopted for the valve motion during the systolic period.

RESULTS

The largest wall shear stress level is detected in the clockwise positioned aortic root at the sinutubular junction. Two counter-rotating vortex cores are formed within the aortic root of both the central and extreme root configurations, however, in the clockwise root the vortex system becomes more symmetric. This also coincides with more entrainment of the valve jet and more turbulence production along the shear layer.

CONCLUSION

A clockwise rotational position of the aortic root imparts an increased wall shear stress at the sinutubular junction and proximal ascending aorta in comparison to other root rotation positions. This may pose increased risk for dilation of the sinutubular junction and ascending aorta in the patient with a clockwise positioned aortic root compared to other normal positional configurations.

Uncertainty quantification of simulated biomechanical stimuli in coronary artery bypass grafts

Coronary artery bypass graft surgery (CABG) is performed on more than 400,000 patients annually in the U.S. However, saphenous vein grafts (SVGs) implanted during CABG exhibit poor patency compared to arterial grafts, with failure rates up to 40% within 10 years after surgery. Differences in mechanical stimuli are known to play a role in driving maladaptation and have been correlated with endothelial damage and thrombus formation. As these quantities are difficult to measure in vivo, multi-scale coronary models offer a way to quantify them, while accounting for complex coronary physiology. However, prior studies have primarily focused on deterministic evaluations, without reporting variability in the model parameters due to uncertainty. This study aims to assess confidence in multi-scale predictions of wall shear stress and wall strain while accounting for uncertainty in peripheral hemodynamics and material properties. Boundary condition distributions are computed by assimilating uncertain clinical data, while spatial variations of vessel wall stiffness are obtained through approximation by a random field. We developed a stochastic submodeling approach to mitigate the computational burden of repeated multi-scale model evaluations to focus exclusively on the bypass grafts. This produces a two-level decomposition of quantities of interest into submodel contributions and full model/submodel discrepancies. We leverage these two levels in the context of forward uncertainty propagation using a previously proposed multi-resolution approach. The time- and space-averaged wall shear stress is well estimated with a coefficient of variation of <35%, but ignorance about the spatial distribution on the wall elastic modulus and thickness lead to large variations in an objective measure of wall strain, with coefficients of variation up to 100%. Sensitivity analysis reveals how the interactions between the flow and material parameters contribute to output variability.

Hemodynamic variables in aneurysms are associated with thrombotic risk in children with Kawasaki disease

BACKGROUND

Thrombosis is a major adverse outcome associated with coronary artery aneurysms (CAAs) resulting from Kawasaki disease (KD). Clinical guidelines recommend initiation of anticoagulation therapy with maximum CAA diameter (D_max) ≥8 mm or Z-score ≥ 10. Here, we investigate the role of aneurysm hemodynamics as a superior method for thrombotic risk stratification in KD patients.

METHODS AND RESULTS

We retrospectively studied ten KD patients with CAAs, including five patients who developed thrombosis. We constructed patient-specific anatomic models from cardiac magnetic resonance images and performed computational hemodynamic simulations using SimVascular. Our simulations incorporated pulsatile flow, deformable arterial walls and boundary conditions automatically tuned to match patient-specific arterial pressure and cardiac output. From simulation results, we derived local hemodynamic variables including time-averaged wall shear stress (TAWSS), low wall shear stress exposure, and oscillatory shear index (OSI). Local TAWSS was significantly lower in CAAs that developed thrombosis (1.2 ± 0.94 vs. 7.28 ± 9.77 dynes/cm², p = 0.006) and the fraction of CAA surface area exposed to low wall shear stress was larger (0.69 ± 0.17 vs. 0.25 ± 0.26%, p = 0.005). Similarly, longer residence times were obtained in branches where thrombosis was confirmed (9.07 ± 6.26 vs. 2.05 ± 2.91 cycles, p = 0.004). No significant differences were found for OSI or anatomical measurements such us D_max and Z-score. Assessment of thrombotic risk according to hemodynamic variables had higher sensitivity and specificity compared to standard clinical metrics (D_max, Z-score).

CONCLUSIONS

Hemodynamic variables can be obtained non-invasively via simulation and may provide improved thrombotic risk stratification compared to current diameter-based metrics, facilitating long-term clinical management of KD patients with persistent CAAs.

Modeling Patient-Specific Magnetic Drug Targeting Within the Intracranial Vasculature

Drug targeting promises to substantially enhance future therapies, for example through the focussing of chemotherapeutic drugs at the site of a tumor, thus reducing the exposure of healthy tissue to unwanted damage. Promising work on the steering of medication in the human body employs magnetic fields acting on nanoparticles made of paramagnetic materials. We develop a computational tool to aid in the optimization of the physical parameters of these particles and the magnetic configuration, estimating the fraction of particles reaching a given target site in a large patient-specific vascular system for different physiological states (heart rate, cardiac output, etc.). We demonstrate the excellent computational performance of our model by its application to the simulation of paramagnetic-nanoparticle-laden flows in a circle of Willis geometry obtained from an MRI scan. The results suggest a strong dependence of the particle density at the target site on the strength of the magnetic forcing and the velocity of the background fluid flow.

Image-based computational assessment of vascular wall mechanics and hemodynamics in pulmonary arterial hypertension patients

Pulmonary arterial hypertension (PAH) is a disease characterized by an elevated pulmonary arterial (PA) pressure. While several computational hemodynamic models of the pulmonary vasculature have been developed to understand PAH, they are lacking in some aspects, such as the vessel wall deformation and its lack of calibration against measurements in humans. Here, we describe a computational modeling framework that addresses these limitations. Specifically, computational models describing the coupling of hemodynamics and vessel wall mechanics in the pulmonary vasculature of a PAH patient and a normal subject were developed. Model parameters, consisting of linearized stiffness E of the large vessels and Windkessel parameters for each outflow branch, were calibrated against in vivo measurements of pressure, flow and vessel wall deformation obtained, respectively, from right-heart catheterization, phase-contrast and cine magnetic resonance images. Calibrated stiffness E of the proximal PA was 2.0 and 0.5 MPa for the PAH and normal models, respectively. Calibrated total compliance C_T and resistance R_T of the distal vessels were, respectively, 0.32 ml/mmHg and 11.3 mmHg∗min/l for the PAH model, and 2.93 ml/mmHg and 2.6 mmHg∗min/l for the normal model. These results were consistent with previous findings that the pulmonary vasculature is stiffer with more constricted distal vessels in PAH patients. Individual effects on PA pressure due to remodeling of the distal and proximal compartments of the pulmonary vasculature were also investigated in a sensitivity analysis. The analysis suggests that the remodeling of distal vasculature contributes more to the increase in PA pressure than the remodeling of proximal vasculature.

RESEARCH ON BIOMECHANICS PROPERTIES AND HEMODYNAMICS PERFORMANCE OF THE CONVERTIBLE VENA CAVA FILTER

This paper reveals the biomechanical properties and hemodynamics performance of the convertible vena cava filter while it was in the process of implantation and conversion in the blood vessels. This paper uses the finite element method and computational fluid dynamics to analyze the interaction mechanism and influence of the convertible vena cava filter while the filter was in the process of implantation and conversion in the blood vessels. Additionally, six pigs were used as experimental samples to verify the effects of the filter in the blood vessels. The computer-aided simulation results showed that it was easier to cause damage to the vessel wall prior to the filter being converted into support. On the contrary, the stress, the peak value of the blood vessel’s stress, the outlet velocity, and the supporting stiffness were reduced after the conversion. Due to the intimal hyperplasia, the supporting element was easy to fix on the inner surface of the blood vessels, which was helpful for the correct positioning of the filter after the conversion. Meanwhile, the animal experiments proved that the surface of inferior vena cava wall was relatively smooth, and the filter did not cause vascular wall rupture. The computer-aided simulation and animal experiments proved the reasonability of the structure of the filter design, and that the filter has good biomechanical properties. The results will provide more scientific reference for the clinical treatment and design of the filter.

Fiber-reinforced computational model of the aortic root incorporating thoracic aorta and coronary structures

Cardiovascular diseases are still the leading causes of death in the developed world. The decline in the mortality associated with circulatory system diseases is accredited to development of new diagnostic and prognostic tools. It is well known that there is an inter relationship between the aortic valve impairment and pathologies of the aorta and coronary vessels. However, due to the limitations of the current tools, the possible link is not fully elucidated. Following our previous model of the aortic root including the coronaries, in this study, we have further developed the global aspect of the model by incorporating the anatomical structure of the thoracic aorta. This model is different from all the previous studies in the sense that inclusion of the coronary structures and thoracic aorta into the natural aortic valve introduces the notion of globality into the model enabling us to explore the possible link between the regional pathologies. The developed model was first validated using the available data in the literature under physiological conditions. Then, to provide a support for the possible association between the localized cardiovascular pathologies and global variations in hemodynamic conditions, we simulated the model for two pathological conditions including moderate and severe aortic valve stenoses. The findings revealed that malformations of the aortic valve are associated with development of low wall shear stress regions and helical blood flow in thoracic aorta that are considered major contributors to aortic pathologies.

Computational Fluid Dynamics and Aortic Thrombus Formation Following Thoracic Endovascular Aortic Repair

BACKGROUND

We present the possible utility of computational fluid dynamics in the assessment of thrombus formation and virtual surgical planning illustrated in a patient with aortic thrombus in a kinked ascending aortic graft following thoracic endovascular aortic repair.

METHODS

A patient-specific three-dimensional model was built from computed tomography. Additionally, we modeled 3 virtual aortic interventions to assess their effect on thrombosis potential: (1) open surgical repair, (2) conformable endografting, and (3) single-branched endografting. Flow waveforms were extracted from echocardiography and used for the simulations. We used the computational index termed platelet activation potential (PLAP) representing accumulated shear rates of fluid particles within a fluid domain to assess thrombosis potential.

RESULTS

The baseline model revealed high PLAP in the entire arch (119.8 ± 42.5), with significantly larger PLAP at the thrombus location (125.4 ± 41.2, p < 0.001). Surgical repair showed a 37% PLAP reduction at the thrombus location (78.6 ± 25.3, p < 0.001) and a 24% reduction in the arch (91.6 ± 28.9, p < 0.001). Single-branched endografting reduced PLAP in the thrombus region by 20% (99.7 ± 24.6, p < 0.001) and by 14% in the arch (103.8 ± 26.1, p < 0.001), whereas a more conformable endograft did not have a profound effect, resulting in a modest 4% PLAP increase (130.6 ± 43.7, p < 0.001) in the thrombus region relative to the baseline case.

CONCLUSIONS

Regions of high PLAP were associated with aortic thrombus. Aortic repair resolved pathologic flow patterns, reducing PLAP. Branched endografting also relieved complex flow patterns reducing PLAP. Computational fluid dynamics may assist in the prediction of aortic thrombus formation in hemodynamically complex cases and help guide repair strategies.

SimVascular: An Open Source Pipeline for Cardiovascular Simulation

Patient-specific cardiovascular simulation has become a paradigm in cardiovascular research and is emerging as a powerful tool in basic, translational and clinical research. In this paper we discuss the recent development of a fully open-source SimVascular software package, which provides a complete pipeline from medical image data segmentation to patient-specific blood flow simulation and analysis. This package serves as a research tool for cardiovascular modeling and simulation, and has contributed to numerous advances in personalized medicine, surgical planning and medical device design. The SimVascular software has recently been refactored and expanded to enhance functionality, usability, efficiency and accuracy of image-based patient-specific modeling tools. Moreover, SimVascular previously required several licensed components that hindered new user adoption and code management and our recent developments have replaced these commercial components to create a fully open source pipeline. These developments foster advances in cardiovascular modeling research, increased collaboration, standardization of methods, and a growing developer community.

Assessment of CardiOvascular Remodelling following Endovascular aortic repair through imaging and computation: the CORE prospective observational cohort study protocol

INTRODUCTION

Thoracic aortic stent grafts are orders of magnitude stiffer than the native aorta. These devices have been associated with acute hypertension, elevated pulse pressure, cardiac remodelling and reduced coronary perfusion. However, a systematic assessment of such cardiovascular effects of thoracic endovascular aortic repair (TEVAR) is missing. The CardiOvascular Remodelling following Endovascular aortic repair (CORE) study aims to (1) quantify cardiovascular remodelling following TEVAR and compare echocardiography against MRI, the reference method; (2) validate computational modelling of cardiovascular haemodynamics following TEVAR using clinical measurements, and virtually assess the impact of more compliant stent grafts on cardiovascular haemodynamics; and (3) investigate diagnostic accuracy of ECG and serum biomarkers for cardiac remodelling compared to MRI.

METHODS AND ANALYSIS

This is a prospective, nonrandomised, observational cohort study. We will use MRI, CT, echocardiography, intraluminal pressures, ECG, computational modelling and serum biomarkers to assess cardiovascular remodelling in two groups of patients with degenerative thoracic aneurysms or penetrating aortic ulcers: (1) patients managed with TEVAR and (2) control patients managed with medical therapy alone. Power analysis revealed a minimum total sample size of 20 patients (α=0.05, power=0.97) to observe significant left ventricular mass increase following TEVAR after 1 year. Consequently, we will include 12 patients in both groups. Advanced MRI sequences will be used to assess myocardial and aortic strain and distensibility, myocardial perfusion and aortic flow. ECG, echocardiography and serum biomarkers will be collected and compared against the imaging data. Computational models will be constructed from each patient imaging data, analysed and validated. All measurements will be collected at baseline (prior to TEVAR) and 1-year follow-up. The expected study period is 3 years.

ETHICS AND DISSEMINATION

This study has been approved by the University of Michigan IRB. The results will be disseminated through scientific journals and conference presentations.

TRIAL REGISTRATION NUMBER

NCT02735720.

Patient-Specific Simulations Reveal Significant Differences in Mechanical Stimuli in Venous and Arterial Coronary Grafts

Mechanical stimuli are key to understanding disease progression and clinically observed differences in failure rates between arterial and venous grafts following coronary artery bypass graft surgery. We quantify biologically relevant mechanical stimuli, not available from standard imaging, in patient-specific simulations incorporating non-invasive clinical data. We couple CFD with closed-loop circulatory physiology models to quantify biologically relevant indices, including wall shear, oscillatory shear, and wall strain. We account for vessel-specific material properties in simulating vessel wall deformation. Wall shear was significantly lower (p = 0.014*) and atheroprone area significantly higher (p = 0.040*) in venous compared to arterial grafts. Wall strain in venous grafts was significantly lower (p = 0.003*) than in arterial grafts while no significant difference was observed in oscillatory shear index. Simulations demonstrate significant differences in mechanical stimuli acting on venous vs. arterial grafts, in line with clinically observed graft failure rates, offering a promising avenue for stratifying patients at risk for graft failure.

Old Myths, New Concerns: the Long-Term Effects of Ascending Aorta Replacement with Dacron Grafts. Not All That Glitters Is Gold

Synthetic grafts are widely used in cardiac and vascular surgery since the mid-1970s. Despite their general good performance, inability of mimicking the elastomechanical characteristics of the native arterial tissue, and the consequent lack of adequate compliance, leads to a cascade of hemodynamic and biological alterations deeply affecting cardiovascular homeostasis. Those concerns have been reconsidered in more contemporaneous surgical and experimental reports which also triggered some research efforts in the tissue engineering field towards the realization of biomimetic arterial surrogates. The present review focuses on the significance of the “compliance mismatch” phenomenon occurring after aortic root or ascending aorta replacement with prosthetic grafts and discusses the clinical reflexes of this state of tissue incompatibility, as the loss of the native elastomechanical properties of the aorta can translate into detrimental effects on the normal efficiency of the aortic root complex with impact in the long-term results of patients undergoing aortic replacement.

The aortic reservoir-wave as a paradigm for arterial haemodynamics: insights from three-dimensional fluid-structure interaction simulations in a model of aortic coarctation

BACKGROUND

The reservoir-wave paradigm considers aortic pressure as the superposition of a 'reservoir pressure', directly related to changes in reservoir volume, and an 'excess' component ascribed to wave dynamics. The change in reservoir pressure is assumed to be proportional to the difference between aortic inflow and outflow (i.e. aortic volume changes), an assumption that is virtually impossible to validate in vivo. The aim of this study is therefore to apply the reservoir-wave paradigm to aortic pressure and flow waves obtained from three-dimensional fluid-structure interaction simulations in a model of a normal aorta, aortic coarctation (narrowed descending aorta) and stented coarctation (stiff segment in descending aorta).

METHOD AND RESULTS

We found no unequivocal relation between the intraaortic volume and the reservoir pressure for any of the simulated cases. When plotted in a pressure-volume diagram, hysteresis loops are found that are looped in a clockwise way indicating that the reservoir pressure is lower than the pressure associated with the change in volume. The reservoir-wave analysis leads to very high excess pressures, especially for the coarctation models, but to surprisingly little changes of the reservoir component despite the impediment of the buffer capacity of the aorta.

CONCLUSION

With the observation that reservoir pressure is not related to the volume in the aortic reservoir in systole, an intrinsic assumption in the wave-reservoir concept is invalidated and, consequently, also the assumption that the excess pressure is the component of pressure that can be attributed to wave travel and reflection.

VALIDATION OF LOSS-COEFFICIENT-BASED OUTLET BOUNDARY CONDITIONS FOR SIMULATING AORTIC FLOW

Flow in a polyurethane model of a human aorta, driven by a heart-lung machine, is analyzed experimentally and computationally for antegrade and retrograde perfusion. The purpose of the analysis is the validation of the previously proposed loss-coefficient-based outlet boundary condition for aortic branches. This model is claimed to be commonly applicable to different perfusion modes of the aorta, unlike the alternative straightforward constant-pressure outlet boundary condition. First, the antegrade perfusion is analyzed computationally and experimentally. This step delivers the loss-coefficients that are to be used in any other perfusion mode of the aorta. Subsequently, a retrograde perfusion is applied to the same aorta, where the flow rates at the outlets of the aortic branches are measured and predicted by applying the loss-coefficient-based outlet boundary conditions. A very good agreement of the predictions with the measurements is observed. The predictions delivered by the standard constant-pressure outlet boundary condition are observed, on the contrary, to be highly in error. Thus, the advocated loss-coefficient-based outlet boundary condition is experimentally validated. It is shown that it is applicable to different perfusion modes with a quite good accuracy, which is much higher compared to the straightforward constant-pressure outlet boundary condition.

Three-dimensional hemodynamics analysis of the circle of Willis in the patient-specific nonintegral arterial structures

The hemodynamic alteration in the cerebral circulation caused by the geometric variations in the cerebral circulation arterial network of the circle of Wills (CoW) can lead to fatal ischemic attacks in the brain. The geometric variations due to impairment in the arterial network result in incomplete cerebral arterial structure of CoW and inadequate blood supply to the brain. Therefore, it is of great importance to understand the hemodynamics of the CoW, for efficiently and precisely evaluating the status of blood supply to the brain. In this paper, three-dimensional computational fluid dynamics of the main CoW vasculature coupled with zero-dimensional lumped parameter model boundary condition for the CoW outflow boundaries is developed for analysis of the blood flow distribution in the incomplete CoW cerebral arterial structures. The geometric models in our study cover the arterial segments from the aorta to the cerebral arteries, which can allow us to take into account the innate patient-specific resistance of the arterial trees. Numerical simulations of the governing fluid mechanics are performed to determine the CoW arterial structural hemodynamics, for illustrating the redistribution of the blood flow in CoW due to the structural variations. We have evaluated our coupling methodology in five patient-specific cases that were diagnosed with the absence of efferent vessels or impairment in the connective arteries in their CoWs. The velocity profiles calculated by our approach in the segments of the patient-specific arterial structures are found to be very close to the Doppler ultrasound measurements. The accuracy and consistency of our hemodynamic results have been improved (to [Formula: see text] %) compared to that of the pure-resistance boundary conditions (of 43.5 [Formula: see text] 28 %). Based on our grouping of the five cases according to the occurrence of unilateral occlusion in vertebral arteries, the inter-comparison has shown that (i) the flow reduction in posterior cerebral arteries is the consequence of the unilateral vertebral arterial occlusion, and (ii) the flow rate in the anterior cerebral arteries is correlated with the posterior structural variations. This study shows that our coupling approach is capable of providing comprehensive information of the hemodynamic alterations in the pathological CoW arterial structures. The information generated by our methodology can enable evaluation of both the functional and structural status of the clinically significant symptoms, for assisting the treatment decision-making.

Non-invasive pressure difference estimation from PC-MRI using the work-energy equation

Pressure difference is an accepted clinical biomarker for cardiovascular disease conditions such as aortic coarctation. Currently, measurements of pressure differences in the clinic rely on invasive techniques (catheterization), prompting development of non-invasive estimates based on blood flow. In this work, we propose a non-invasive estimation procedure deriving pressure difference from the work-energy equation for a Newtonian fluid. Spatial and temporal convergence is demonstrated on in silico Phase Contrast Magnetic Resonance Image (PC-MRI) phantoms with steady and transient flow fields. The method is also tested on an image dataset generated in silico from a 3D patient-specific Computational Fluid Dynamics (CFD) simulation and finally evaluated on a cohort of 9 subjects. The performance is compared to existing approaches based on steady and unsteady Bernoulli formulations as well as the pressure Poisson equation. The new technique shows good accuracy, robustness to noise, and robustness to the image segmentation process, illustrating the potential of this approach for non-invasive pressure difference estimation.

Differential impact of local stiffening and narrowing on hemodynamics in repaired aortic coarctation: an FSI study

Even after successful treatment of aortic coarctation, a high risk of cardiovascular morbidity and mortality remains. Uncertainty exists on the factors contributing to this increased risk among which are the presence of (1) a residual narrowing leading to an additional resistance and (2) a less distensible zone disturbing the buffer function of the aorta. As the many interfering factors and adaptive physiological mechanisms present in vivo prohibit the study of the isolated impact of these individual factors, a numerical fluid-structure interaction model is developed to predict central hemodynamics in coarctation treatment. The overall impact of a stiffening on the hemodynamics is limited, with a small increase in systolic pressure (up to 8 mmHg) proximal to the stiffening which is amplified with increasing stiffening and length. A residual narrowing, on the other hand, affects the hemodynamics significantly. For a short segment (10 mm), the combination of a stiffening and narrowing (coarctation index 0.5) causes an increase in systolic pressure of 58 mmHg, with 31 mmHg due to narrowing and an additional 27 mmHg due to stiffening. For a longer segment (25 mm), an increase in systolic pressure of 50 mmHg is found, of which only 9 mmHg is due to stiffening.

Quantitative Assessment of Turbulence and Flow Eccentricity in an Aortic Coarctation: Impact of Virtual Interventions

Turbulence and flow eccentricity can be measured by magnetic resonance imaging (MRI) and may play an important role in the pathogenesis of numerous cardiovascular diseases. In the present study, we propose quantitative techniques to assess turbulent kinetic energy (TKE) and flow eccentricity that could assist in the evaluation and treatment of stenotic severities. These hemodynamic parameters were studied in a pre-treated aortic coarctation (CoA) and after several virtual interventions using computational fluid dynamics (CFD), to demonstrate the effect of different dilatation options on the flow field. Patient-specific geometry and flow conditions were derived from MRI data. The unsteady pulsatile flow was resolved by large eddy simulation including non-Newtonian blood rheology. Results showed an inverse asymptotic relationship between the total amount of TKE and degree of dilatation of the stenosis, where turbulent flow proximal the constriction limits the possible improvement by treating the CoA alone. Spatiotemporal maps of TKE and flow eccentricity could be linked to the characteristics of the jet, where improved flow conditions were favored by an eccentric dilatation of the CoA. By including these flow markers into a combined MRI-CFD intervention framework, CoA therapy has not only the possibility to produce predictions via simulation, but can also be validated pre- and immediate post treatment, as well as during follow-up studies.

Simulations reveal adverse hemodynamics in patients with multiple systemic to pulmonary shunts

For newborns diagnosed with pulmonary atresia or severe pulmonary stenosis leading to insufficient pulmonary blood flow, cyanosis can be mitigated with placement of a modified Blalock-Taussig shunt (MBTS) between the innominate and pulmonary arteries. In some clinical scenarios, patients receive two systemic-to-pulmonary connections, either by leaving the patent ductus arteriosus (PDA) open or by adding an additional central shunt (CS) in conjunction with the MBTS. This practice has been motivated by the thinking that an additional source of pulmonary blood flow could beneficially increase pulmonary flow and provide the security of an alternate pathway in case of thrombosis. However, there have been clinical reports of premature shunt occlusion when more than one shunt is employed, leading to speculation that multiple shunts may in fact lead to unfavorable hemodynamics and increased mortality. In this study, we hypothesize that multiple shunts may lead to undesirable flow competition, resulting in increased residence time (RT) and elevated risk of thrombosis, as well as pulmonary overcirculation. Computational fluid dynamics-based multiscale simulations were performed to compare a range of shunt configurations and systematically quantify flow competition, pulmonary circulation, and other clinically relevant parameters. In total, 23 cases were evaluated by systematically changing the PDA/CS diameter, pulmonary vascular resistance (PVR), and MBTS position and compared by quantifying oxygen delivery (OD) to the systemic and coronary beds, wall shear stress (WSS), oscillatory shear index (OSI), WSS gradient (WSSG), and RT in the pulmonary artery (PA), and MBTS. Results showed that smaller PDA/CS diameters can lead to flow conditions consistent with increased thrombus formation due to flow competition in the PA, and larger PDA/CS diameters can lead to insufficient OD due to pulmonary hyperfusion. In the worst case scenario, it was found that multiple shunts can lead to a 160% increase in RT and a 10% decrease in OD. Based on the simulation results presented in this study, clinical outcomes for patients receiving multiple shunts should be critically investigated, as this practice appears to provide no benefit in terms of OD and may actually increase thrombotic risk.

Simulation of short-term pressure regulation during the tilt test in a coupled 3D-0D closed-loop model of the circulation

Short-term fluctuations in arterial pressures arising from normal physiological function are buffered by a negative feedback system known as the arterial baroreflex. Initiated by altered biomechanical stretch in the vessel wall, the baroreflex coordinates a systemic response that alters heart rate, cardiac contractility and peripheral vessel vasoconstriction. In this work, a coupled 3D–0D formulation for the short-term pressure regulation of the systemic circulation is presented. Including the baroreflex feedback mechanisms, a patient-specific model of the large arteries is subjected to a simulated head up tilt test. Comparative simulations with and without baroreflex control highlight the critical role that the baroreflex has in regulating variations in pressures within the systemic circulation.

Perspective on CFD studies of coronary artery disease lesions and hemodynamics: a review

Coronary artery disease (CAD) is the most common cardiovascular disease. Early diagnosis of CAD's physiological significance is of utmost importance for guiding individualized risk-tailored treatment strategies. In this paper, we first review the state-of-the-art clinical diagnostic indices to quantify the severity of CAD and the associated invasive and noninvasive imaging technologies in order to quantify the anatomical parameters of diameter stenosis, area stenosis, and hemodynamic indices of coronary flow reserve and fractional flow reserve. With the development of computational technologies and CFD methods, tremendous progress has been made in applying image-based CFD simulation techniques to elucidate the effects of hemodynamics in vascular pathophysiology toward the initialization and progression of CAD. So then, we review the advancements of CFD technologies in patient-specific modeling, involving the development of geometry reconstruction, boundary conditions, and fluid-structure interaction. Next, we review the applications of CFD to stenotic sites, in order to compute their hemodynamic parameters and study the relationship between the hemodynamic conditions and the clinical indices, to thereby assess the amount of viable myocardium and candidacy for percutaneous coronary intervention. Finally, we review the strengths and limitations of current researches of applying CFD to CAD studies.

Quantification of regional differences in aortic stiffness in the aging human

There has been a growing awareness over the past decade that stiffening of the aorta, and its attendant effects on hemodynamics, is both an indicator and initiator of diverse cardiovascular, neurovascular, and renovascular diseases. Although different clinical metrics of arterial stiffness have been proposed and found useful in particular situations, there remains a need to understand better the complex interactions between evolving aortic stiffness and the hemodynamics. Computational fluid-solid-interaction (FSI) models are amongst the most promising means to understand such interactions for one can parametrically examine effects of regional variations in material properties and arterial geometry on local and systemic blood pressure and flow. Such models will not only increase our understanding, they will also serve as important steps towards the development of fluid-solid-growth (FSG) models that can further examine interactions between the evolving wall mechanics and hemodynamics that lead to arterial adaptations or disease progression over long periods. In this paper, we present a consistent quantification and comparison of regional nonlinear biaxial mechanical properties of the human aorta based on 19 data sets available in the literature and we calculate associated values of linearized stiffness over the cardiac cycle that are useful for initial large-scale FSI and FSG simulations. It is shown, however, that there is considerable variability amongst the available data and consequently that there is a pressing need for more standardized biaxial testing of the human aorta to collect data as a function of both location and age, particularly for young healthy individuals who serve as essential controls.