Effect of Wall Elasticity on Hemodynamics and Wall Shear Stress in Patient-Specific Simulations in the Coronary Arteries

Hemodynamics and wall shear metrics in a pulmonary autograft: Comparing a fluid-structure interaction and computational fluid dynamics approach

OBJECTIVE

In young patients, aortic valve disease is often treated by placement of a pulmonary autograft (PA) which adapts to its new environment through growth and remodeling. To better understand the hemodynamic forces acting on the highly distensible PA in the acute phase after surgery, we developed a fluid-structure interaction (FSI) framework and comprehensively compared hemodynamics and wall shear-stress (WSS) metrics with a computational fluid dynamic (CFD) simulation.

METHODS

The FSI framework couples a prestressed non-linear hyperelastic arterial tissue model with a fluid model using the in-house coupling code CoCoNuT. Geometry, material parameters and boundary conditions are based on in-vivo measurements. Hemodynamics, time-averaged WSS (TAWSS), oscillatory shear index (OSI) and topological shear variation index (TSVI) are evaluated qualitatively and quantitatively for 3 different sheeps.

RESULTS

Despite systolic-to-diastolic volumetric changes of the PA in the order of 20 %, the point-by-point correlation of TAWSS and OSI obtained through CFD and FSI remains high (r > 0.9, p < 0.01) for TAWSS and (r > 0.8, p < 0.01) for OSI). Instantaneous WSS divergence patterns qualitatively preserve similarities, but large deformations of the PA leads to a decrease of the correlation between FSI and CFD resolved TSVI (r < 0.7, p < 0.01). Moderate co-localization between FSI and CFD is observed for low thresholds of TAWSS and high thresholds of OSI and TSVI.

CONCLUSION

FSI might be warranted if we were to use the TSVI as a mechano-biological driver for growth and remodeling of PA due to varying intra-vascular flow structures and near wall hemodynamics because of the large expansion of the PA.

Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use

PURPOSE

Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential.

METHODS

This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed.

RESULTS

Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data.

CONCLUSIONS

Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.

Integrating Computational and Biological Hemodynamic Approaches to Improve Modeling of Atherosclerotic Arteries

Atherosclerosis is the primary cause of cardiovascular disease, resulting in mortality, elevated healthcare costs, diminished productivity, and reduced quality of life for individuals and their communities. This is exacerbated by the limited understanding of its underlying causes and limitations in current therapeutic interventions, highlighting the need for sophisticated models of atherosclerosis. This review critically evaluates the computational and biological models of atherosclerosis, focusing on the study of hemodynamics in atherosclerotic coronary arteries. Computational models account for the geometrical complexities and hemodynamics of the blood vessels and stenoses, but they fail to capture the complex biological processes involved in atherosclerosis. Different in vitro and in vivo biological models can capture aspects of the biological complexity of healthy and stenosed vessels, but rarely mimic the human anatomy and physiological hemodynamics, and require significantly more time, cost, and resources. Therefore, emerging strategies are examined that integrate computational and biological models, and the potential of advances in imaging, biofabrication, and machine learning is explored in developing more effective models of atherosclerosis.

Predicting Lipid-Rich Plaque Progression in Coronary Arteries Using Multimodal Imaging and Wall Shear Stress Signatures

BACKGROUND

Plaque composition and wall shear stress (WSS) magnitude act as well-established players in coronary plaque progression. However, WSS magnitude per se does not completely capture the mechanical stimulus to which the endothelium is subjected, since endothelial cells experience changes in the WSS spatiotemporal configuration on the luminal surface. This study explores WSS profile and lipid content signatures of plaque progression to identify novel biomarkers of coronary atherosclerosis.

METHODS

Thirty-seven patients with acute coronary syndrome underwent coronary computed tomography angiography, near-infrared spectroscopy intravascular ultrasound, and optical coherence tomography of at least 1 nonculprit vessel at baseline and 1-year follow-up. Baseline coronary artery geometries were reconstructed from intravascular ultrasound and coronary computed tomography angiography and combined with flow information to perform computational fluid dynamics simulations to assess the timeaveraged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, which can be assessed by the topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage atheroma volume change at 1-year follow-up (Δplaque atheroma volume). Plaque composition information was extracted from near-infrared spectroscopy and optical coherence tomography.

RESULTS

Exposure to high TSVI and low TAWSS was associated with higher plaque progression (4.00±0.69% and 3.60±0.62%, respectively). Plaque composition acted synergistically with TSVI or TAWSS, resulting in the highest plaque progression (≥5.90%) at locations where lipid-rich content is exposed to high TSVI or low TAWSS.

CONCLUSIONS

Luminal exposure to high TSVI, solely or combined with a lipid-rich plaque phenotype, is associated with enhanced plaque progression at 1-year follow-up. Where plaque progression occurred, low TAWSS was also observed. These findings suggest TSVI, in addition to low TAWSS, as a potential biomechanical predictor for plaque progression, showing promise for clinical translation to improve patient prognosis.

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.

The Influence of Aortic Valve Disease on Coronary Hemodynamics: A Computational Model-Based Study

Aortic valve disease (AVD) often coexists with coronary artery disease (CAD), but whether and how the two diseases are correlated remains poorly understood. In this study, a zero-three dimensional (0-3D) multi-scale modeling method was developed to integrate coronary artery hemodynamics, aortic valve dynamics, coronary flow autoregulation mechanism, and systemic hemodynamics into a unique model system, thereby yielding a mathematical tool for quantifying the influences of aortic valve stenosis (AS) and aortic valve regurgitation (AR) on hemodynamics in large coronary arteries. The model was applied to simulate blood flows in six patient-specific left anterior descending coronary arteries (LADs) under various aortic valve conditions (i.e., control (free of AVD), AS, and AR). Obtained results showed that the space-averaged oscillatory shear index (SA-OSI) was significantly higher under the AS condition but lower under the AR condition in comparison with the control condition. Relatively, the overall magnitude of wall shear stress was less affected by AVD. Further data analysis revealed that AS induced the increase in OSI in LADs mainly through its role in augmenting the low-frequency components of coronary flow waveform. These findings imply that AS might increase the risk or progression of CAD by deteriorating the hemodynamic environment in coronary arteries.

Abnormal Wall Shear Stress Area is Correlated to Coronary Artery Bypass Graft Remodeling 1 Year After Surgery

Coronary artery bypass graft surgery is a common intervention for coronary artery disease; however, it suffers from graft failure, and the underlying mechanisms are not fully understood. To better understand the relation between graft hemodynamics and surgical outcomes, we performed computational fluid dynamics simulations with deformable vessel walls in 10 study participants (24 bypass grafts) based on CT and 4D flow MRI one month after surgery to quantify lumen diameter, wall shear stress (WSS), and related hemodynamic measures. A second CT acquisition was performed one year after surgery to quantify lumen remodeling. Compared to venous grafts, left internal mammary artery grafts experienced lower abnormal WSS (< 1 Pa) area one month after surgery (13.8 vs. 70.1%, p = 0.001) and less inward lumen remodeling one year after surgery (- 2.4% vs. - 16.1%, p = 0.027). Abnormal WSS area one month post surgery correlated with percent change in graft lumen diameter one year post surgery (p = 0.030). This study shows for the first time prospectively a correlation between abnormal WSS area one month post surgery and graft lumen remodeling 1 year post surgery, suggesting that shear-related mechanisms may play a role in post-operative graft remodeling and might help explain differences in failure rates between arterial and venous grafts.

Impact of cyclic bending on coronary hemodynamics

It remains unknown that the degree of bias in computational fluid dynamics results without considering coronary cyclic bending. This study aims to investigate the influence of different rates of coronary cyclic bending on coronary hemodynamics. To model coronary bending, a multi-ring-controlled fluid-structural interaction model was designed. A coronary artery was simulated with various cyclic bending rates (0.5, 0.75 and 1 s, corresponding to heart rates of 120, 80 and 60 bpm) and compared against a stable model. The simulated results show that the hemodynamic parameters of vortex Q-criterion, temporal wall shear stress (WSS), time-averaged WSS (TaWSS) and oscillatory shear index (OSI) were sensitive to the changes in cyclic rate. A higher heart rate resulted in higher magnitude and larger variance in the hemodynamic parameters. Whereas, the values and distributions of flow velocity and relative residence time (RRT) did not show significant differences between different bending periods. This study suggests that a stable coronary model is not sufficient to represent the hemodynamics in a bending coronary artery. Different heart rate conditions were found to have significant impact on the hemodynamic parameters. Thus, cyclic bending should be considered to mimic the realistic hemodynamics in future patient-specific coronary hemodynamics studies.

Colocalization of Coronary Plaque with Wall Shear Stress in Myocardial Bridge Patients

PURPOSE

Patients with myocardial bridges (MBs) have a higher prevalence of atherosclerosis. Wall shear stress (WSS) has previously been correlated with plaque in coronary artery disease patients, but such correlations have not been investigated in symptomatic MB patients. The aim of this paper was to use a multi-scale computational fluid dynamics (CFD) framework to simulate hemodynamics in MB patient, and investigate the co-localization of WSS and plaque.

METHODS

We identified N = 10 patients from a previously reported cohort of 50 symptomatic MB patients, all of whom had plaque in the proximal vessel. Dynamic 3D models were reconstructed from coronary computed tomography angiography (CCTA), intravascular ultrasound (IVUS) and catheter angiograms. CFD simulations were performed to compute WSS proximal to, within and distal to the MB. Plaque was quantified from IVUS images in 2 mm segments and registered to CFD model. Plaque area was compared to absolute and patient-normalized WSS.

RESULTS

WSS was lower in the proximal segment compared to the bridge segment (6.1 ± 2.9 vs. 16.0 ± 7.1 dynes/cm², p value < 0.01). Plaque area and plaque burden measured from IVUS peaked at 1-3 cm proximal to the MB entrance, coinciding with the first diagonal branch. Normalized WSS showed a statistically significant moderate correlation with plaque area (r = 0.41, p < 0.01).

CONCLUSION

WSS may be obtained non-invasively in MB patients and provides a surrogate marker of plaque area. Using CFD, it may be possible to non-invasively assess the extent of plaque area, and identify patients who could benefit from frequent monitoring or medical management.

Numerical study of the effect of hemodynamic variables on LDL concentration through the single layer of the Left Anterior Descending coronary artery (LAD) under the heart pulse

Heart attack is one of the most common causes of death in the world. Coronary artery disease is the most recognized cause of heart attack whose onset and progression have been attributed to low-density lipoprotein (LDL) passing through the wall of the artery. In this paper, hemodynamic variables as well as the concentration of LDL through the coronary porous artery at the Left Anterior Descending coronary artery (LAD), and its first diagonal branch (D1) under the heart motion investigated using computational simulation. The geometry that has been studied in this paper is the first bifurcation of Left Anterior Descending (LAD) that has been placed on a perimeter of hypothetical sphere representative of the heart geometry. Sinusoidal variations of sphere radii, simulated pulsating movement of the heart. Blood has been considered as a Newtonian and incompressible flow with pulsatile flow rate and real physiological profile. The plasma filtration boundary condition used over the walls in order to simulate the concentration of LDL to a one-layer artery wall. Variations in the concentration of LDL on the artery wall and its relation to oscillation on shear stress on the artery wall under different conditions are presented. Moreover, the effects of the pulsating inlet flow and dynamic movement of the artery are explored. The results declared that minimum shear stress and maximum LDL concentration take place at the bifurcation and on the myocardial wall which is in complete agreement with clinical studies. Furthermore, it has been shown that the heart pulse has a slight effect on the average time of concentration (0.1% increase); however, by analyzing all time steps, one could observe that the maximum concentration rises in some time steps; where this increases the possibility of LDL presence and helps them diffuse inside the artery wall.

Framework for lumen-based nonrigid tomographic coregistration of intravascular images

Purpose: Modern medical imaging enables clinicians to effectively diagnose, monitor, and treat diseases. However, clinical decision-making often relies on combined evaluation of either longitudinal or disparate image sets, necessitating coregistration of multiple acquisitions. Promising coregistration techniques have been proposed; however, available methods predominantly rely on time-consuming manual alignments or nontrivial feature extraction with limited clinical applicability. Addressing these issues, we present a fully automated, robust, nonrigid registration method, allowing for coregistering of multimodal tomographic vascular image datasets using luminal annotation as the sole alignment feature. Approach: Registration is carried out by the use of the registration metrics defined exclusively for lumens shapes. The framework is primarily broken down into two sequential parts: longitudinal and rotational registration. Both techniques are inherently nonrigid in nature to compensate for motion and acquisition artifacts in tomographic images. Results: Performance was evaluated across multimodal intravascular datasets, as well as in longitudinal cases assessing pre-/postinterventional coronary images. Low registration error in both datasets highlights method utility, with longitudinal registration errors-evaluated throughout the paired tomographic sequences-of 0.29 ± 0.14    mm ( < 2 longitudinal image frames) and 0.18 ± 0.16    mm ( < 1 frame) for multimodal and interventional datasets, respectively. Angular registration for the interventional dataset rendered errors of 7.7 ° ± 6.7 ° , and 29.1 ° ± 23.2 ° for the multimodal set. Conclusions: Satisfactory results across datasets, along with additional attributes such as the ability to avoid longitudinal over-fitting and correct nonlinear catheter rotation during nonrigid rotational registration, highlight the potential wide-ranging applicability of our presented coregistration method.

INVESTIGATION OF FLOW IN AN IDEALIZED CURVED ARTERY: COMPARATIVE STUDY USING CFD AND FSI WITH NEWTONIAN AND NON-NEWTONIAN FLUIDS

Artery curvatures, where disturbed flow patterns are expected, are preferred sites of formation of atherosclerosis. Experimental studies have shown that low and oscillating wall shear stress (WSS) plays an important role in the development and progression of atherosclerosis. Accurate estimation of these biomechanical parameters is important to assess the risk of atherosclerosis formation. The coupled effects of non-Newtonian behavior of blood and artery wall flexibility for the transient blood flow through an idealized curved coronary artery are investigated using computational fluid dynamics (CFD) as well as fluid–structure interaction (FSI) simulations. The choice of fluid model, Carreau and Newtonian, was found to impact the time averaged and minimum WSS values. The effects of wall deformation on time averaged wall shear tress were negligible. However, a comparison of temporal minima of WSS along the curvature showed significant variations between CFD and FSI simulations. Since low WSS values are crucial in the prediction of atherosclerosis development, it is concluded that both the non-Newtonian behavior of blood and the wall flexibility should be considered for computational studies.

Numerical predictions of shear stress and cyclic stretch in pulmonary hypertension due to left heart failure

Isolated post-capillary pulmonary hypertension (Ipc-PH) occurs due to left heart failure, which contributes to 1 out of every 9 deaths in the United States. In some patients, through unknown mechanisms, Ipc-PH transitions to combined pre-/post-capillary PH (Cpc-PH) and is associated with a dramatic increase in mortality. Altered mechanical forces and subsequent biological signaling in the pulmonary vascular bed likely contribute to the transition from Ipc-PH to Cpc-PH. However, even in a healthy pulmonary circulation, the mechanical forces in the smallest vessels (the arterioles, capillary bed, and venules) have not been quantitatively defined. This study is the first to examine this question via a computational fluid dynamics model of the human pulmonary arteries, arterioles, venules, and veins. Using this model, we predict temporal and spatial dynamics of cyclic stretch and wall shear stress with healthy and diseased hemodynamics. In the normotensive case for large vessels, numerical simulations show that large arteries have higher pressure and flow than large veins, as well as more pronounced changes in area throughout the cardiac cycle. In the microvasculature, shear stress increases and cyclic stretch decreases as vessel radius decreases. When we impose an increase in left atrial pressure to simulate Ipc-PH, shear stress decreases and cyclic stretch increases as compared to the healthy case. Overall, this model predicts pressure, flow, shear stress, and cyclic stretch that providing a way to analyze and investigate hypotheses related to disease progression in the pulmonary circulation.

A Pilot Study Characterizing Flow Patterns in the Thoracic Aorta of Patients With Connective Tissue Disease: Comparison to Age- and Gender-Matched Controls via Fluid Structure Interaction

Prior computational and imaging studies described changes in flow patterns for patients with Marfan syndrome, but studies are lacking for related populations. This pilot study addresses this void by characterizing wall shear stress (WSS) indices for patients with Loeys-Dietz and undifferentiated connective tissue diseases. Using aortic valve-based velocity profiles from magnetic resonance imaging as input to patient-specific fluid structure interaction (FSI) models, we determined local flow patterns throughout the aorta for four patients with various connective tissue diseases (Loeys-Dietz with the native aorta, connective tissue disease of unclear etiology with native aorta in female and male patients, and an untreated patient with Marfan syndrome, as well as twin patients with Marfan syndrome who underwent valve-sparing root replacement). FSI simulations used physiological boundary conditions and material properties to replicate available measurements. Time-averaged WSS (TAWSS) and oscillatory shear index (OSI) results are presented with localized comparison to age- and gender-matched control participants. Ascending aortic dimensions were greater in almost all patients with connective tissue diseases relative to their respective control. Differences in TAWSS and OSI were driven by local morphological differences and cardiac output. For example, the model for one twin had a more pronounced proximal descending aorta in the vicinity of the ductus ligamentum that impacted WSS indices relative to the other. We are optimistic that the results of this study can serve as a foundation for larger future studies on the connective tissue disorders presented in this article.

Sex-related differences in plaque characteristics and endothelial shear stress related plaque-progression in human coronary arteries

BACKGROUND AND AIMS

Clinical atherosclerosis manifestations are different in women compared to men. Since endothelial shear stress (ESS) is known to play a critical role in coronary atherosclerosis development, we investigated differences in anatomical characteristics and endothelial shear stress (ESS)-related plaque growth in human coronary arteries in men compared to women.

METHODS

1183 coronary arteries (male/female: 944/239) from the PREDICTION study were studied for differences in artery/plaque and ESS characteristics, and ESS-related plaque progression (6-10 months follow-up) among men and women and after stratification for age. All characteristics were derived from IVUS-based vascular profiling and reported per 3 mm-segments (13,030 3-mm-segments (male/female: 10,465/2,565)).

RESULTS

Coronary arteries and plaques were significantly smaller in females compared to males; but no important differences were observed in plaque burden, ESS and rate of plaque progression. Change in plaque burden was inversely related to ESS (p<0.001) with no difference between women versus men (β: -0.62 ± 0.13 vs -0.68 ± 0.05, p=0.62). However, stratification for age demonstrated that ESS-related plaque growth was more marked in young women compared to men (<55 years, β: -2.02 ± 0.61 vs -0.33 ± 0.10, p=0.007), reducing in magnitude over the age-categories up till 75 years.

CONCLUSIONS

Coronary artery and plaque size are smaller in women compared to men, but ESS and ESS- related plaque progression were similar. Sex-related differences in ESS-related plaque growth were evident after stratification for age. These observations suggest that although the fundamental processes of atherosclerosis progression are similar in men versus women, plaque progression may be influenced by age within gender.

An optimal control approach to determine resistance-type boundary conditions from in-vivo data for cardiovascular simulations

The choice of appropriate boundary conditions is a fundamental step in computational fluid dynamics (CFD) simulations of the cardiovascular system. Boundary conditions, in fact, highly affect the computed pressure and flow rates, and consequently haemodynamic indicators such as wall shear stress (WSS), which are of clinical interest. Devising automated procedures for the selection of boundary conditions is vital to achieve repeatable simulations. However, the most common techniques do not automatically assimilate patient-specific data, relying instead on expensive and time-consuming manual tuning procedures. In this work, we propose a technique for the automated estimation of outlet boundary conditions based on optimal control. The values of resistive boundary conditions are set as control variables and optimized to match available patient-specific data. Experimental results on four aortic arches demonstrate that the proposed framework can assimilate 4D-Flow MRI data more accurately than two other common techniques based on Murray's law and Ohm's law.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Morphological and hemodynamic analysis of the patient-specific renal cell carcinoma models

Although the morbidity of renal cell carcinoma (RCC) has been increasing as the seventh most common tumours, to our knowledge, there is few studies foucsing on the hemodynamics of the renal artery (RA) with RCC. The objective of this study is to perform morphological and hemodynamic analysis of the RA and abdominal aorta artery (AAA) in the control healthy and RCC patient groups. Three-dimensional (3D) geometries are reconstructed from 18 control healthy subjects and 15 RCC patients based on Computed Tomography Angiography (CTA) images. There is higer in the lumen diameter of the RA (6.21 ± 0.89 mm) and curvature of the RA (1.2 ± 0.07) in the RCC patient group compared with the control healthy group (4.29 ± 1.08 mm, 1.1 ± 0.1), respectively. In the hemodynamic analysis, the surface area ratio (%) of low time-averaged wall shear stress (SAR-TAWSS) at the RA (10.65 ± 11.65) and AAA (48.49 ± 12.79) in the RCC patient group is significantly higher than that in the control healthy group (0.23 ± 0.22, 21.57 ± 20.5), respectively. It is found that RCC altered the morphology of the RA in the RCC patient group, which could deteriorate the hemodynamic environment of the RA and AAA. The finding in this study could enhance us to understand the progression of vascular disease caused by RCC.

Effect of stenosis and dilatation on the hemodynamic parameters associated with left coronary artery

BACKGROUND AND OBJECTIVE

The main objective of the work is to examine the curvature effects of stenosis/dilatation region pertaining to left coronary artery. The hemodymamic features during the cardiac cycle is thoroughly examined.

METHODS

A numerical fluid structure interaction model incorporating multi- layered elastic artery wall, non-Newtonian blood viscosity and pulsating boundary conditions is developed. The composite arterial wall consists of a thin layer tunica intima, atheroma and a thick wall. Higher stiffness of atheroma is captured by using higher Young's modulus. The CFD and FSI models are validated with available experimental and analytical data. Computations are done with five different non-Newtonian models and arterial wall with various elasticity levels. The local and time averaged WSS, velocity contours downstream of stenosis, wall pressure and pressure drop during various phases of cardiac cycle are provided in detail.

RESULTS

The influence of non-Newtonian effects of blood viscosity is found to be significant especially at stenosis regions. The flexible wall caused wall deformation and the associated flow and pressure wave propagation affecting WSS and pressure drop compared to the rigid wall. Flow recirculation is noticed at stenosis downstream locations and its strength increases with increased severity of the stenosis. A stenosis is characterised by a sudden drop in wall pressure and a slower two stage recovery during peak velocity periods of the cardiac cycle.

CONCLUSIONS

The pressure drop, local WSS at stenosis centre, and radial velocity increase are significantly higher for stenosis cases and the effect is severe during peak diastole. The variation in hemodynamic parameters is found to be less significant for dilatation. Significantly lower WSS is noticed for the recirculation regions downstream of stenosis which can enhance the tendency for monocytes to attach to the endothelium. The radius of curvature of the stenosis is found to be the most sensitive parameter affecting the hemodynamic characteristics rather than the detailed geometry of the stenosis. The main effect of variation of artery wall stiffness is noted at recirculation regions present downstream of stenosis. The results from the study may be useful for predicting wall shear stress signatures associated with stenosis/dilatation changes and the management of specific cases.

Non-invasive characterization of complex coronary lesions

Conventional invasive diagnostic imaging techniques do not adequately resolve complex Type B and C coronary lesions, which present unique challenges, require personalized treatment and result in worsened patient outcomes. These lesions are often excluded from large-scale non-invasive clinical trials and there does not exist a validated approach to characterize hemodynamic quantities and guide percutaneous intervention for such lesions. This work identifies key biomarkers that differentiate complex Type B and C lesions from simple Type A lesions by introducing and validating a coronary angiography-based computational fluid dynamic (CFD-CA) framework for intracoronary assessment in complex lesions at ultrahigh resolution. Among 14 patients selected in this study, 7 patients with Type B and C lesions were included in the complex lesion group including ostial, bifurcation, serial lesions and lesion where flow was supplied by collateral bed. Simple lesion group included 7 patients with lesions that were discrete, [Formula: see text] long and readily accessible. Intracoronary assessment was performed using CFD-CA framework and validated by comparing to clinically measured pressure-based index, such as FFR. Local pressure, endothelial shear stress (ESS) and velocity profiles were derived for all patients. We validates the accuracy of our CFD-CA framework and report excellent agreement with invasive measurements ([Formula: see text]). Ultra-high resolution achieved by the model enable physiological assessment in complex lesions and quantify hemodynamic metrics in all vessels up to 1mm in diameter. Importantly, we demonstrate that in contrast to traditional pressure-based metrics, there is a significant difference in the intracoronary hemodynamic forces, such as ESS, in complex lesions compared to simple lesions at both resting and hyperemic physiological states [n = 14, [Formula: see text]]. Higher ESS was observed in the complex lesion group ([Formula: see text] Pa) than in simple lesion group ([Formula: see text] Pa). Complex coronary lesions have higher ESS compared to simple lesions, such differential hemodynamic evaluation can provide much the needed insight into the increase in adverse outcomes for such patients and has incremental prognostic value over traditional pressure-based indices, such as FFR.

Application of physics-based flow models in cardiovascular medicine: Current practices and challenges

Personalized physics-based flow models are becoming increasingly important in cardiovascular medicine. They are a powerful complement to traditional methods of clinical decision-making and offer a wealth of physiological information beyond conventional anatomic viewing using medical imaging data. These models have been used to identify key hemodynamic biomarkers, such as pressure gradient and wall shear stress, which are associated with determining the functional severity of cardiovascular diseases. Importantly, simulation-driven diagnostics can help researchers understand the complex interplay between geometric and fluid dynamic parameters, which can ultimately improve patient outcomes and treatment planning. The possibility to compute and predict diagnostic variables and hemodynamics biomarkers can therefore play a pivotal role in reducing adverse treatment outcomes and accelerate development of novel strategies for cardiovascular disease management.

A Hemodynamic Comparison of Myocardial Bridging and Coronary Atherosclerotic Stenosis: A Computational Model With Experimental Evaluation

Myocardial bridging (MB) and coronary atherosclerotic stenosis can impair coronary blood flow and may cause myocardial ischemia or even heart attack. It remains unclear how MB and stenosis are similar or different regarding their impacts on coronary hemodynamics. The purpose of this study was to compare the hemodynamic effects of coronary stenosis and MB using experimental and computational fluid dynamics (CFD) approaches. For CFD modeling, three MB patients with different levels of lumen obstruction, mild, moderate, and severe were selected. Patient-specific left anterior descending (LAD) coronary artery models were reconstructed from biplane angiograms. For each MB patient, the virtually healthy and stenotic models were also simulated for comparison. In addition, an in vitro flow-loop was developed, and the pressure drop was measured for comparison. The CFD simulations results demonstrated that the difference between MB and stenosis increased with increasing MB/stenosis severity and flowrate. Experimental results showed that increasing the MB length (by 140%) only had significant impact on the pressure drop in the severe MB (39% increase at the exercise), but increasing the stenosis length dramatically increased the pressure drop in both moderate and severe stenoses at all flow rates (31% and 93% increase at the exercise, respectively). Both CFD and experimental results confirmed that the MB had a higher maximum and a lower mean pressure drop in comparison with the stenosis, regardless of the degree of lumen obstruction. A better understanding of MB and atherosclerotic stenosis may improve the therapeutic strategies in coronary disease patients and prevent acute coronary syndromes.

Computational investigation of left ventricular hemodynamics following bioprosthetic aortic and mitral valve replacement

The left ventricle of the heart is a fundamental structure in the human cardiac system that pumps oxygenated blood into the systemic circulation. Several valvular conditions can cause the aortic and mitral valves associated with the left ventricle to become severely diseased and require replacement. However, the clinical outcomes of such operations, specifically the postoperative ventricular hemodynamics of replacing both valves, are not well understood. This work uses computational fluid-structure interaction (FSI) to develop an improved understanding of this effect by modeling a left ventricle with the aortic and mitral valves replaced with bioprostheses. We use a hybrid Arbitrary Lagrangian-Eulerian/immersogeometric framework to accommodate the analysis of cardiac hemodynamics and heart valve structural mechanics in a moving fluid domain. The motion of the endocardium is obtained from a cardiac biomechanics simulation and provided as an input to the proposed numerical framework. The results from the simulations in this work indicate that the replacement of the native mitral valve with a tri-radially symmetric bioprosthesis dramatically changes the ventricular hemodynamics. Most significantly, the vortical motion in the left ventricle is found to reverse direction after mitral valve replacement. This study demonstrates that the proposed computational FSI framework is capable of simulating complex multiphysics problems and can provide an in-depth understanding of the cardiac mechanics.

The Story of Wall Shear Stress in Coronary Artery Atherosclerosis: Biochemical Transport and Mechanotransduction

Coronary artery atherosclerosis is a local, multifactorial, complex disease, and the leading cause of death in the US. Complex interactions between biochemical transport and biomechanical forces influence disease growth. Wall shear stress (WSS) affects coronary artery atherosclerosis by inducing endothelial cell mechanotransduction and by controlling the near-wall transport processes involved in atherosclerosis. Each of these processes is controlled by WSS differently and therefore has complicated the interpretation of WSS in atherosclerosis. In this paper, we present a comprehensive theory for WSS in atherosclerosis. First, a short review of shear stress-mediated mechanotransduction in atherosclerosis was presented. Next, subject-specific computational fluid dynamics (CFD) simulations were performed in ten coronary artery models of diseased and healthy subjects. Biochemical-specific mass transport models were developed to study low-density lipoprotein, nitric oxide, adenosine triphosphate, oxygen, monocyte chemoattractant protein-1, and monocyte transport. The transport results were compared with WSS vectors and WSS Lagrangian coherent structures (WSS LCS). High WSS magnitude protected against atherosclerosis by increasing the production or flux of atheroprotective biochemicals and decreasing the near-wall localization of atherogenic biochemicals. Low WSS magnitude promoted atherosclerosis by increasing atherogenic biochemical localization. Finally, the attracting WSS LCS's role was more complex where it promoted or prevented atherosclerosis based on different biochemicals. We present a summary of the different pathways by which WSS influences coronary artery atherosclerosis and compare different mechanotransduction and biotransport mechanisms.

Multiscale modeling of blood flow to assess neurological complications in patients supported by venoarterial extracorporeal membrane oxygenation

Computational blood flow models in large arteries elucidate valuable relationships between cardiovascular diseases and hemodynamics, leading to improvements in treatment planning and clinical decision making. One such application with potential to benefit from simulation is venoarterial extracorporeal membrane oxygenation (VA-ECMO), a support system for patients with cardiopulmonary failure. VA-ECMO patients develop high rates of neurological complications, partially due to abnormal blood flow throughout the vasculature from the VA-ECMO system. To better understand these hemodynamic changes, it is important to resolve complex local flow parameters derived from three-dimensional (3D) fluid dynamics while also capturing the impact of VA-ECMO support throughout the systemic arterial system. As high-resolution 3D simulations of the arterial network remain computationally expensive and intractable for large studies, a validated, multiscale model is needed to compute both global effects and high-fidelity local hemodynamics. In this work, we developed and demonstrated a framework to model hemodynamics in VA-ECMO patients using coupled 3D and one-dimensional (1D) models (1D→3D). We demonstrated the ability of these multiscale models to simulate complex flow patterns in specific regions of interest while capturing bulk flow throughout the systemic arterial system. We compared 1D, 3D, and 1D→3D coupled models and found that multiscale models were able to sufficiently capture both global and local hemodynamics in the cerebral arteries and aorta in VA-ECMO patients. This study is the first to develop and compare 1D, 3D, and 1D→ 3D coupled models on the larger arterial system scale in VA-ECMO patients, with potential use for other large scale applications.

Physiology and coronary artery disease: emerging insights from computed tomography imaging based computational modeling

Improvements in spatial and temporal resolution now permit robust high quality characterization of presence, morphology and composition of coronary atherosclerosis in computed tomography (CT). These characteristics include high risk features such as large plaque volume, low CT attenuation, napkin-ring sign, spotty calcification and positive remodeling. Because of the high image quality, principles of patient-specific computational fluid dynamics modeling of blood flow through the coronary arteries can now be applied to CT and allow the calculation of local lesion-specific hemodynamics such as endothelial shear stress, fractional flow reserve and axial plaque stress. This review examines recent advances in coronary CT image-based computational modeling and discusses the opportunity to identify lesions at risk for rupture much earlier than today through the combination of anatomic and hemodynamic information.

Validation of Wall Shear Stress Assessment in Non-invasive Coronary CTA versus Invasive Imaging: A Patient-Specific Computational Study

Endothelial shear stress (ESS) identifies coronary plaques at high risk for progression and/or rupture leading to a future acute coronary syndrome. In this study an optimized methodology was developed to derive ESS, pressure drop and oscillatory shear index using computational fluid dynamics (CFD) in 3D models of coronary arteries derived from non-invasive coronary computed tomography angiography (CTA). These CTA-based ESS calculations were compared to the ESS calculations using the gold standard with fusion of invasive imaging and CTA. In 14 patients paired patient-specific CFD models based on invasive and non-invasive imaging of the left anterior descending (LAD) coronary arteries were created. Ten patients were used to optimize the methodology, and four patients to test this methodology. Time-averaged ESS (TAESS) was calculated for both coronary models applying patient-specific physiological data available at the time of imaging. For data analysis, each 3D reconstructed coronary artery was divided into 2 mm segments and each segment was subdivided into 8 arcs (45°).TAESS and other hemodynamic parameters were averaged per segment as well as per arc. Furthermore, the paired segment- and arc-averaged TAESS were categorized into patient-specific tertiles (low, medium and high). In the ten LADs, used for optimization of the methodology, we found high correlations between invasively-derived and non-invasively-derived TAESS averaged over segments (n = 263, r = 0.86) as well as arcs (n = 2104, r = 0.85, p < 0.001). The correlation was also strong in the four testing-patients with r = 0.95 (n = 117 segments, p = 0.001) and r = 0.93 (n = 936 arcs, p = 0.001).There was an overall high concordance of 78% of the three TAESS categories comparing both methodologies using the segment- and 76% for the arc-averages in the first ten patients. This concordance was lower in the four testing patients (64 and 64% in segment- and arc-averaged TAESS). Although the correlation and concordance were high for both patient groups, the absolute TAESS values averaged per segment and arc were overestimated using non-invasive vs. invasive imaging [testing patients: TAESS segment: 30.1(17.1-83.8) vs. 15.8(8.8-63.4) and TAESS arc: 29.4(16.2-74.7) vs 15.0(8.9-57.4) p < 0.001]. We showed that our methodology can accurately assess the TAESS distribution non-invasively from CTA and demonstrated a good correlation with TAESS calculated using IVUS/OCT 3D reconstructed models.

High precision invasive FFR, low-cost invasive iFR, or non-invasive CFR?: optimum assessment of coronary artery stenosis based on the patient-specific computational models

The objective of this paper is to apply computational fluid dynamic (CFD) as a complementary tool for clinical tests to not only predict the present and future status of left coronary artery stenosis but also to evaluate some clinical hypotheses. In order to assess the present status of the coronary artery stenosis severity, and thereby selecting the most appropriate type of treatment for each patient, fractional flow reserve (FFR), instantaneous wave free-ratio (iFR), and coronary flow reserve (CFR) are calculated. To examine FFR, iFR, and CFR results, the effect of geometric features of stenoses, including diameter reduction (%), lesion length (LL), and minimum lumen diameter (MLD), is studied on them. It is observed that FFR is a more conservative index than iFR and CFR to assess the severity of coronary stenosis. In addition, it is seen that FFR, iFR, and CFR decrease by increasing LL and decreasing MLD. Therefore, the morphological indices, LL/MLD and LL/MLD̂4, with the calculated conservative cut-off values equal to 5.5 and 3.6, are considered. Next, some controversial clinical hypotheses about the assessment of the severity of coronary stenosis are evaluated numerically. These include the examination of FFR, iFR, and CFR accuracies, investigating the effect of coronary hyperemia on iFR, as well as the reliability of the hybrid iFR-FFR decision-making strategy. The presented numerical model can also be used as a predictive tool to identify the atherosuseptible sites of arteries by calculating the time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

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.

Does the inflow velocity profile influence physiologically relevant flow patterns in computational hemodynamic models of left anterior descending coronary artery?

Patient-specific computational fluid dynamics is a powerful tool for investigating the hemodynamic risk in coronary arteries. Proper setting of flow boundary conditions in computational hemodynamic models of coronary arteries is one of the sources of uncertainty weakening the findings of in silico experiments, in consequence of the challenging task of obtaining in vivo 3D flow measurements within the clinical framework. Accordingly, in this study we evaluated the influence of assumptions on inflow velocity profile shape on coronary artery hemodynamics. To do that, (1) ten left anterior descending coronary artery (LAD) geometries were reconstructed from clinical angiography, and (2) eleven velocity profiles with realistic 3D features such as eccentricity and differently shaped (single- and double-vortex) secondary flows were generated analytically and imposed as inflow boundary conditions. Wall shear stress and helicity-based descriptors obtained prescribing the commonly used parabolic velocity profile were compared with those obtained with the other velocity profiles. Our findings indicated that the imposition of idealized velocity profiles as inflow boundary condition is acceptable as long the results of the proximal vessel segment are not considered, in LAD coronary arteries. As a pragmatic rule of thumb, a conservative estimation of the length of influence of the shape of the inflow velocity profile on LAD local hemodynamics can be given by the theoretical entrance length for cylindrical conduits in laminar flow conditions.

Analysis and manegement of laminar blood flow inside a cerebral blood vessel using a finite volume software program for biomedical engineering

BACKGROUND AND OBJECTIVE

Hemodynamic blood flow analysis in the cerebrovascular is has become one of the important research topics in the bio-mechanic in recent decades. The primary duty of the cerebral blood vessel is supplying Glucose and oxygen for the brain.

METHODS

In this investigation, the non-Newtonian blood flow in the cerebral blood vessels studied. For modeling the geometry of this problem, we used Magnetic Resonance Image (MRI) approach to take Digital Imaging and Communications in Medicine (DICOM) images and using an open-source software package to construct the geometry, which is a complicated one. The power-law indexes, heat flux, and Reynolds number range in the investigation are 0.6 ≤ n ≤ 0.8, 5 ≤ q ≤ 15Wm^-2 and 160≤Re≤310. Effects of Reynolds number, power-law indexes and heat fluxes are investigated.

RESULTS

We found that the pressure drop increase with increasing the Reynolds number and power-law index. The maximum Nusselt number in the cerebral blood vessels accrued in the running position of the body in n = 0.8. Also, the highest average wall shear stress occurs in maximum power-law indexes and Reynolds number.

CONCLUSION

By increasing the power-law index and Reynolds number, the wall shear stress increases.

Accelerating massively parallel hemodynamic models of coarctation of the aorta using neural networks

Comorbidities such as anemia or hypertension and physiological factors related to exertion can influence a patient’s hemodynamics and increase the severity of many cardiovascular diseases. Observing and quantifying associations between these factors and hemodynamics can be difficult due to the multitude of co-existing conditions and blood flow parameters in real patient data. Machine learning-driven, physics-based simulations provide a means to understand how potentially correlated conditions may affect a particular patient. Here, we use a combination of machine learning and massively parallel computing to predict the effects of physiological factors on hemodynamics in patients with coarctation of the aorta. We first validated blood flow simulations against in vitro measurements in 3D-printed phantoms representing the patient’s vasculature. We then investigated the effects of varying the degree of stenosis, blood flow rate, and viscosity on two diagnostic metrics – pressure gradient across the stenosis (ΔP) and wall shear stress (WSS) - by performing the largest simulation study to date of coarctation of the aorta (over 70 million compute hours). Using machine learning models trained on data from the simulations and validated on two independent datasets, we developed a framework to identify the minimal training set required to build a predictive model on a per-patient basis. We then used this model to accurately predict ΔP (mean absolute error within 1.18 mmHg) and WSS (mean absolute error within 0.99 Pa) for patients with this disease.

Coronary artery plaque growth: A two-way coupled shear stress-driven model

Atherosclerosis in coronary arteries can lead to plaque growth, stenosis formation, and blockage of the blood flow supplying the heart tissue. Several studies have shown that hemodynamics play an important role in the growth of coronary artery plaques. Specifically, low wall shear stress (WSS) appears to be the leading hemodynamic parameter promoting atherosclerotic plaque growth, which in turn influences the blood flow and WSS distribution. Therefore, a two-way coupled interaction exists between WSS and atherosclerosis growth. In this work, a computational framework was developed to study the coupling between WSS and plaque growth in coronary arteries. Computational fluid dynamics (CFD) was used to quantify WSS distribution. Surface mesh nodes were moved in the inward normal direction according to a growth model based on WSS. After each growth stage, the geometry was updated and the CFD simulation repeated to find updated WSS values for the next growth stage. One hundred twenty growth stages were simulated in an idealized tube and an image-based left anterior descending artery. An automated framework was developed using open-source software to couple CFD simulations with growth. Changes in plaque morphology and hemodynamic patterns during different growth stages are presented. The results show larger plaque growth towards the downstream segment of the plaque, agreeing with the reported clinical observations. The developed framework could be used to establish hemodynamic-driven growth models and study the interaction between these processes.