Multi-patient study for coronary vulnerable plaque model comparisons: 2D/3D and fluid–structure interaction simulations

Several image-based computational models have been used to perform mechanical analysis for atherosclerotic plaque progression and vulnerability investigations. However, differences of computational predictions from those models have not been quantified at multi-patient level. In vivo intravascular ultrasound (IVUS) coronary plaque data were acquired from seven patients. Seven 2D/3D models with/without circumferential shrink, cyclic bending and fluid–structure interactions (FSI) were constructed for the seven patients to perform model comparisons and quantify impact of 2D simplification, circumferential shrink, FSI and cyclic bending plaque wall stress/strain (PWS/PWSn) and flow shear stress (FSS) calculations. PWS/PWSn and FSS averages from seven patients (388 slices for 2D and 3D thin-layer models) were used for comparison. Compared to 2D models with shrink process, 2D models without shrink process overestimated PWS by 17.26%. PWS change at location with greatest curvature change from 3D FSI models with/without cyclic bending varied from 15.07% to 49.52% for the seven patients (average = 30.13%). Mean Max-FSS, Min-FSS and Ave-FSS from the flow-only models under maximum pressure condition were 4.02%, 11.29% and 5.45% higher than those from full FSI models with cycle bending, respectively. Mean PWS and PWSn differences between FSI and structure-only models were only 4.38% and 1.78%. Model differences had noticeable patient variations. FSI and flow-only model differences were greater for minimum FSS predictions, notable since low FSS is known to be related to plaque progression. Structure-only models could provide PWS/PWSn calculations as good approximations to FSI models for simplicity and time savings in calculation.

     

Fluid-structure interaction models based on patient-specific IVUS at baseline and follow-up for prediction of coronary plaque progression by morphological and biomechanical factors: A preliminary study

Plaque morphology and biomechanics are believed to be closely associated with plaque progression. In this paper, we test the hypothesis that integrating morphological and biomechanical risk factors would result in better predictive power for plaque progression prediction. A sample size of 374 intravascular ultrasound (IVUS) slices was obtained from 9 patients with IVUS follow-up data. 3D fluid-structure interaction models were constructed to obtain both structural stress/strain and fluid biomechanical conditions. Data for eight morphological and biomechanical risk factors were extracted for each slice. Plaque area increase (PAI) and wall thickness increase (WTI) were chosen as two measures for plaque progression. Progression measure and risk factors were fed to generalized linear mixed models and linear mixed-effect models to perform prediction and correlation analysis, respectively. All combinations of eight risk factors were exhausted to identify the optimal predictor(s) with highest prediction accuracy defined as sum of sensitivity and specificity. When using a single risk factor, plaque wall stress (PWS) at baseline was the best predictor for plaque progression (PAI and WTI). The optimal predictor among all possible combinations for PAI was PWS + PWSn + Lipid percent + Min cap thickness + Plaque Area (PA) + Plaque Burden (PB) (prediction accuracy = 1.5928) while Wall Thickness (WT) + Plaque Wall Strain (PWSn) + Plaque Area (PA) was the best for WTI (1.2589). This indicated that PAI was a more predictable measure than WTI. The combination including both morphological and biomechanical parameters had improved prediction accuracy, compared to predictions using only morphological features.     

A negative correlation between human carotid atherosclerotic plaque progression and plaque wall stress: in vivo MRI-based 2D/3D FSI models

It is well accepted that atherosclerosis initiation and progression correlate positively with low and oscillating flow wall shear stresses (FSS). However, this mechanism cannot explain why advanced plaques continue to grow under elevated FSS conditions. In vivo magnetic resonance imaging (MRI)-based 2D/3D multi-component models with fluid-structure interactions (FSI, 3D only) for human carotid atherosclerotic plaques were introduced to quantify correlations between plaque progression measured by wall thickness increase (WTI) and plaque wall (structure) stress (PWS) conditions. A histologically validated multi-contrast MRI protocol was used to acquire multi-year in vivo MRI images. Our results using 2D models (200-700 data points/patient) indicated that 18 out of 21 patients studied showed significant negative correlation between WTI and PWS at time 2 (T2). The 95% confidence interval for the Pearson correlation coefficient is (-0.443,-0.246), p<0.0001. Our 3D FSI model supported the 2D correlation results and further indicated that combining both plaque structure stress and flow shear stress gave better approximation results (PWS, T2: R(2)=0.279; FSS, T1: R(2)=0.276; combining both: R(2)=0.637). These pilot studies suggest that both lower PWS and lower FSS may contribute to continued plaque progression and should be taken into consideration in future investigations of diseases related to atherosclerosis.     

Advanced human carotid plaque progression correlates positively with flow shear stress using follow-up scan data: An in vivo MRI multi-patient 3D FSI study

Although it has been well-accepted that atherosclerosis initiation and early progression correlate negatively with flow wall shear stresses (FSS), increasing evidence suggests mechanisms governing advanced plaque progression are not well understood. Fourteen patients were scanned 2–4 times at 18 month intervals using a histologically validated multi-contrast magnetic resonance imaging (MRI) protocol to acquire carotid plaque progression data. Thirty-two scan pairs (baseline and follow-up scans) were formed with slices matched for model construction and analysis. 3D fluid–structure interaction (FSI) models were constructed and plaque wall stress (PWS) and flow shear stress (FSS) were obtained from all matching lumen data points (400–1000 per plaque; 100 points per matched slice) to quantify correlations with plaque progression measured by vessel wall thickness increase (WTI). Using FSS and PWS data from follow-up scan, 21 out of 32 scan pairs showed a significant positive correlation between WTI and FSS (positive/negative/no significance ratio=21/8/3), and 26 out of 32 scan pairs showed a significant negative correlation between WTI and PWS (positive/negative/no significance ratio=2/26/4). The mean FSS value of lipid core nodes (n=5294) from all 47 plaque models was 63.5 dyn/cm², which was 45% higher than that from all normal vessel nodes (n=27553, p<0.00001). The results from this intensive FSI study indicate that flow shear stress from follow-up scan correlates positively with advanced plaque progression which is different from what has been observed in plaque initiation and early-stage progression. It should be noted that the correlation results do not automatically lead to any causality conclusions.     

In vivo serial MRI-based models and statistical methods to quantify sensitivity and specificity of mechanical predictors for carotid plaque rupture: location and beyond

It has been hypothesized that mechanical risk factors may be used to predict future atherosclerotic plaque rupture. Truly predictive methods for plaque rupture and methods to identify the best predictor(s) from all the candidates are lacking in the literature. A novel combination of computational and statistical models based on serial magnetic resonance imaging (MRI) was introduced to quantify sensitivity and specificity of mechanical predictors to identify the best candidate for plaque rupture site prediction. Serial in vivo MRI data of carotid plaque from one patient was acquired with follow-up scan showing ulceration. 3D computational fluid-structure interaction (FSI) models using both baseline and follow-up data were constructed and plaque wall stress (PWS) and strain (PWSn) and flow maximum shear stress (FSS) were extracted from all 600 matched nodal points (100 points per matched slice, baseline matching follow-up) on the lumen surface for analysis. Each of the 600 points was marked "ulcer" or "nonulcer" using follow-up scan. Predictive statistical models for each of the seven combinations of PWS, PWSn, and FSS were trained using the follow-up data and applied to the baseline data to assess their sensitivity and specificity using the 600 data points for ulcer predictions. Sensitivity of prediction is defined as the proportion of the true positive outcomes that are predicted to be positive. Specificity of prediction is defined as the proportion of the true negative outcomes that are correctly predicted to be negative. Using probability 0.3 as a threshold to infer ulcer occurrence at the prediction stage, the combination of PWS and PWSn provided the best predictive accuracy with (sensitivity, specificity)?=?(0.97, 0.958). Sensitivity and specificity given by PWS, PWSn, and FSS individually were (0.788, 0.968), (0.515, 0.968), and (0.758, 0.928), respectively. The proposed computational-statistical process provides a novel method and a framework to assess the sensitivity and specificity of various risk indicators and offers the potential to identify the optimized predictor for plaque rupture using serial MRI with follow-up scan showing ulceration as the gold standard for method validation. While serial MRI data with actual rupture are hard to acquire, this single-case study suggests that combination of multiple predictors may provide potential improvement to existing plaque assessment schemes. With large-scale patient studies, this predictive modeling process may provide more solid ground for rupture predictor selection strategies and methods for image-based plaque vulnerability assessment.     

Morphological and Stress Vulnerability Indices for Human Coronary Plaques and Their Correlations with Cap Thickness and Lipid Percent: An IVUS-Based Fluid-Structure Interaction Multi-patient Study

Plaque vulnerability, defined as the likelihood that a plaque would rupture, is difficult to quantify due to lack of in vivo plaque rupture data. Morphological and stress-based plaque vulnerability indices were introduced as alternatives to obtain quantitative vulnerability assessment. Correlations between these indices and key plaque features were investigated. In vivo intravascular ultrasound (IVUS) data were acquired from 14 patients and IVUS-based 3D fluid-structure interaction (FSI) coronary plaque models with cyclic bending were constructed to obtain plaque wall stress/strain and flow shear stress for analysis. For the 617 slices from the 14 patients, lipid percentage, min cap thickness, critical plaque wall stress (CPWS), strain (CPWSn) and flow shear stress (CFSS) were recorded, and cap index, lipid index and morphological index were assigned to each slice using methods consistent with American Heart Association (AHA) plaque classification schemes. A stress index was introduced based on CPWS. Linear Mixed-Effects (LME) models were used to analyze the correlations between the mechanical and morphological indices and key morphological factors associated with plaque rupture. Our results indicated that for all 617 slices, CPWS correlated with min cap thickness, cap index, morphological index with r = -0.6414, 0.7852, and 0.7411 respectively (p<0.0001). The correlation between CPWS and lipid percentage, lipid index were weaker (r = 0.2445, r = 0.2338, p<0.0001). Stress index correlated with cap index, lipid index, morphological index positively with r = 0.8185, 0.3067, and 0.7715, respectively, all with p<0.0001. For all 617 slices, the stress index has 66.77% agreement with morphological index. Morphological and stress indices may serve as quantitative plaque vulnerability assessment supported by their strong correlations with morphological features associated with plaque rupture. Differences between the two indices may lead to better plaque assessment schemes when both indices were jointly used with further validations from clinical studies.     

Enhanced External Counterpulsation Treatment May Intervene The Advanced Atherosclerotic Plaque Progression by Inducing The Variations of Mechanical Factors: A 3D FSI Study Based on in vivo Animal Experiment

Growing evidences suggest that long-term enhanced external counterpulsation (EECP) treatment can inhibit the initiation of atherosclerotic lesion by improving the hemodynamic environment in aortas. However, whether this kind procedure will intervene the progression of advanced atherosclerotic plaque remains elusive and causes great concern in its clinical application presently. In the current paper, a pilot study combining animal experiment and numerical simulation was conducted to investigate the acute mechanical stress variations during EECP intervention, and then to assess the possible chronic effects. An experimentally induced hypercholesterolemic porcine model was developed and the basic hemodynamic measurement was performed in vivo before and during EECP treatment. Meanwhile, A 3D fluid-structure interaction (FSI) model of blood vessel with symmetric local stenosis was developed for the numerical calculation of some important mechanical factors. The results show that EECP augmented 12.21% of the plaque wall stress (PWS), 57.72% of the time average wall shear stress (AWSS) and 43.67% of the non-dimensional wall shear stress gradient (WSSGnd) at throat site of the stenosis. We suggest that long-term EECP treatment may intervene the advanced plaque progression by inducing the significant variations of some important mechanical factors, but its proper effects will need a further research combined follow-up observation in clinic.     

The influence of computational strategy on prediction of mechanical stress in carotid atherosclerotic plaques: Comparison of 2D structure-only, 3D structure-only,one-way and fully coupled fluid-structure interaction analyses

Background

Compositional and morphological features of carotid atherosclerotic plaques provide complementary information to luminal stenosis in predicting clinical presentations. However, they alone cannot predict cerebrovascular risk. Mechanical stress within the plaque induced by cyclical changes in blood pressure has potential to assess plaque vulnerability. Various modeling strategies have been employed to predict stress, including 2D and 3D structure-only, 3D one-way and fully coupled fluid-structure interaction (FSI) simulations. However, differences in stress predictions using different strategies have not been assessed.

Methods

Maximum principal stress (Stress-P1) within 8 human carotid atherosclerotic plaques was calculated based on geometry reconstructed from in vivo computerized tomography and high resolution, multi-sequence magnetic resonance images. Stress-P1 within the diseased region predicted by 2D and 3D structure-only, and 3D one-way FSI simulations were compared to 3D fully coupled FSI analysis.

Results

Compared to 3D fully coupled FSI, 2D structure-only simulation significantly overestimated stress level (94.1 kPa [65.2, 117.3] vs. 85.5 kPa [64.4, 113.6]; median [inter-quartile range], p=0.0004). However, when slices around the bifurcation region were excluded, stresses predicted by 2D structure-only simulations showed a good correlation (R²=0.69) with values obtained from 3D fully coupled FSI analysis. 3D structure-only model produced a small yet statistically significant stress overestimation compared to 3D fully coupled FSI (86.8 kPa [66.3, 115.8] vs. 85.5 kPa [64.4, 113.6]; p<0.0001). In contrast, one-way FSI underestimated stress compared to 3D fully coupled FSI (78.8 kPa [61.1, 100.4] vs. 85.5 kPa [64.4, 113.7]; p<0.0001).

Conclusions

A 3D structure-only model seems to be a computationally inexpensive yet reasonably accurate approximation for stress within carotid atherosclerotic plaques with mild to moderate luminal stenosis as compared to fully coupled FSI analysis.     

Higher critical plaque wall stress in patients who died of coronary artery disease compared with those who died of other causes: A 3D FSI study based on ex vivo MRI of coronary plaques

Mechanical forces play an important role in the rupture of vulnerable plaques. This process is often associated with cardiovascular syndromes, such as heart attack and stroke. In this study, magnetic resonance imaging (MRI)-based models were used to investigate the association between plaque wall stress (PWS) and coronary artery disease (CAD).     

Haemodynamic assessment of human coronary arteries is affected by degree of freedom of artery movement

Abnormal haemodynamic parameters are associated with atheroma plaque progression and instability in coronary arteries. Flow recirculation, shear stress and pressure gradient are understood to be important pathogenic mediators in coronary disease. The effect of freedom of coronary artery movement on these parameters is still unknown. Fluid–structure interaction (FSI) simulations were carried out in 25 coronary artery models derived from authentic human coronaries in order to investigate the effect of degree of freedom of movement of the coronary arteries on flow recirculation, wall shear stress (WSS) and wall pressure gradient (WPG). Each FSI model had distinctive supports placed upon it. The quantitative and qualitative differences in flow recirculation, maximum wall shear stress (MWSS), areas of low wall shear stress (ALWSS) and maximum wall pressure gradient (MWPG) for each model were determined. The results showed that greater freedom of movement was associated with lower MWSS, smaller ALWSS, smaller flow recirculation zones and lower MWPG. With increasing percentage diameter stenosis (%DS), the effect of degree of freedom on flow recirculation and WSS diminished. Freedom of movement is an important variable to be considered for computational modelling of human coronary arteries, especially in the setting of mild to moderate stenosis.

Abbreviations: 3D: Three-dimensional; 3DR: Three-dimensional Reconstruction; 3D-QCA: Three-dimensional quantitative coronary angiography; ALWSS: Areas of low wall shear stress; CAD: Coronary artery disease; CFD: Computational fluid dynamics; %DS: Diameter stenosis percentage; EPCS: End point of counter-rotating streamlines; FSI: Fluid–structure interaction; IVUS: Intravascular ultrasound; LAD: Left anterior descending; MWSS: Maximum wall shear stress; SST: Shear stress transport; TAWSS: Time-averaged wall shear stress; WSS: wall shear stress; WPG: Wall pressure gradient; MWPG: Maximum wall pressure gradient; FFR: Fractional flow reserve; iFR: Instantaneous wave-free ratio     

Cyclic Bending Contributes to High Stress in a Human Coronary Atherosclerotic Plaque and Rupture Risk: In Vitro Experimental Modeling and Ex Vivo MRI-Based Computational Modeling Approach

Many acute cardiovascular syndromes such as heart attack and stroke are caused by atherosclerotic plaque ruptures which often happen without warning. MRI-based models with fluid-structure interactions (FSI) have been introduced to perform flow and stress/strain analysis for atherosclerotic plaques and identify possible mechanical and morphological indices for accurate plaque vulnerability assessment. In this paper, cyclic bending was added to 3D FSI coronary plaque models for more accurate mechanical predictions. Curvature variation was prescribed using the data of a human left anterior descending (LAD) coronary artery. Five computational models were constructed based on ex vivo MRI human coronary plaque data to assess the effects of cyclic bending, pulsating pressure, plaque structure, and axial stretch on plaque stress/strain distributions. In vitro experiments using a hydrogel stenosis model with cyclical bending were performed to observe effect of cyclical bending on flow conditions. Our results indicate that cyclical bending may cause more than 100% or even up to more than 1000% increase in maximum principal stress values at locations where the plaque is bent most. Stress increase is higher when bending is coupled with axial stretch, non-smooth plaque structure, or resonant pressure conditions (zero phase angle shift). Effects of cyclic bending on flow behaviors are more modest (21.6% decrease in maximum velocity, 10.8% decrease in flow rate, maximum flow shear stress changes were < 5%). Computational FSI models including cyclic bending, plaque components and structure, axial stretch, accurate in vivo measurements of pressure, curvature, and material properties should lead to significant improvement on stress-based plaque mechanical analysis and more accurate coronary plaque vulnerability assessment.     

Measurement of the transverse strain tensor in the coronary arterial wall from clinical intravascular ultrasound images

Atherosclerotic plaque rupture is the major cause of acute coronary syndromes. Currently, there is no reliable diagnostic tool to predict plaque rupture. Knowledge of plaque mechanical properties based on local artery wall strain measurements would be useful for characterizing its composition and predicting its vulnerability. Due to cardiac motion, strain estimation in clinical intravascular ultrasound (IVUS) images is extremely challenging. A method is presented to estimate cross-sectional coronary artery wall strain in response to cardiac pulsatile pressure using clinically acquired IVUS images, which are acquired in continuous pullback mode. First, cardiac phase information is retrieved retrospectively from an IVUS image sequence using an image-based gating method, and image sub-sequences at systole and diastole are extracted. Then, images at branch sites are used as landmarks to align the two image sub-sequences. Finally, the paired images at each site are registered to measure the 2D strain tensor of the coronary artery cross-section. This method has been successfully applied to IVUS images of a left anterior descending (LAD) coronary artery acquired clinically during a standard procedure. Such complete strain information should be useful for identifying vulnerable plaque.     

Patient-Specific Artery Shrinkage and 3D Zero-Stress State in Multi-Component 3D FSI Models for Carotid Atherosclerotic Plaques Based on <i>In Vivo</i> MRI Data

Image-based computational models for atherosclerotic plaques have been developed to perform mechanical analysis to quantify critical flow and stress/strain conditions related to plaque rupture which often leads directly to heart attack or stroke. An important modeling issue is how to determine zero stress state from in vivo plaque geometries. This paper presents a method to quantify human carotid artery axial and inner circumferential shrinkages by using patient-specific ex vivo and in vivo MRI images. A shrink-stretch process based on patient-specific in vivo plaque morphology and shrinkage data was introduced to shrink the in vivo geometry first to find the zero-stress state (opening angle was ignored to reduce the complexity), and then stretch and pressurize to recover the in vivo plaque geometry with computed initial stress, strain, flow pressure and velocity conditions. Effects of the shrink-stretch process on plaque stress/strain distributions were demonstrated based on patient-specific data using 3D models with fluid-structure interactions (FSI). The average artery axial and inner circumferential shrinkages were 25% and 7.9%, respectively, based on a data set obtained from 10 patients. Maximum values of maximum principal stress and strain increased 349.8% and 249% respectively with 33% axial stretch. Influence of inner circumferential shrinkage (7.9%) was not very noticeable under 33% axial stretch, but became more noticeable under smaller axial stretch. Our results indicated that accurate knowledge of artery shrinkages and the shrink-stretch process will considerably improve the accuracy of computational predictions made based on results from those in vivo MRI-based FSI models.     

Stress-Based Plaque Vulnerability Index and Assessment for Carotid Atherosclerotic Plaques Using Patient-Specific Vessel Material Properties

Cardiovascular diseases are closely linked to atherosclerotic plaque development and rupture. Assessment of plaque vulnerability is of fundamental significance to cardiovascular research and disease diagnosis, prevention, treatment and management. Magnetic resonance image (MRI) data of carotid atherosclerotic plaques from 8 patients (5 male, 3 female; age: 62-83, mean=71) were acquired at the University of Washington (UW), Seattle by the Vascular Imaging Laboratory (VIL) with written informed consent obtained. Patient-specific vessel material properties were quantified using Cine MRI data for modeling use. 3D thin-layer models were used to obtain plaque stress and strain for plaque assessment. A stress-based plaque vulnerability index (SPVI) was proposed to combine mechanical analysis, plaque morphology and composition for more complete carotid plaque vulnerability assessment. The five intervals (unit: kPa) [0, 46.8), [46.8, 80), [80, 92), [92, 103), and [103, +∞) from in vivo material models were used for SPVI values of 0, 1, 2, 3 and 4, respectively. The optimized agreement rate was 85.19%. The use of patient-specific material properties in plaque models could potentially improve the accuracy of model stress/strain calculations. SPVI has the potential to improve the current image-based screening and plaque vulnerability assessment schemes.     

Cap inflammation leads to higher plaque cap strain and lower cap stress: An MRI-PET/CT-based FSI modeling approach

Plaque rupture may be triggered by extreme stress/strain conditions. Inflammation is also implicated and can be imaged using novel imaging techniques. The impact of cap inflammation on plaque stress/strain and flow shear stress were investigated. A patient-specific MRI-PET/CT-based modeling approach was used to develop 3D fluid-structure interaction models and investigate the impact of inflammation on plaque stress/strain conditions for better plaque assessment. 18FDG-PET/CT and MRI data were acquired from 4 male patients (average age: 66) to assess plaque characteristics and inflammation. Material stiffness for the fibrous cap was adjusted lower to reflect cap weakening causing by inflammation. Setting stiffness ratio (SR) to be 1.0 (fibrous tissue) for baseline, results for SR=0.5, 0.25, and 0.1 were obtained. Thin cap and hypertension were also considered. Combining results from the 4 patients, mean cap stress from 729 cap nodes was lowered by 25.2% as SR went from 1.0 to 0.1. Mean cap strain value for SR=0.1 was 0.313, 114% higher than that from SR=1.0 model. The thin cap SR=0.1 model had 40% mean cap stress decrease and 81% cap strain increase compared with SR=1.0 model. The hypertension SR=0.1 model had 19.5% cap stress decrease and 98.6% cap strain increase compared with SR=1.0 model. Differences of flow shear stress with 4 different SR values were limited (<10%). Cap inflammation may lead to large cap strain conditions when combined with thin cap and hypertension. Inflammation also led to lower cap stress. This shows the influence of inflammation on stress/strain calculations which are closely related to plaque assessment.     

3D Fluid-Structure Interaction Canine Heart Model with Patch to Quantify Mechanical Conditions for Optimal Myocardium Stem Cell Growth and Tissue Regeneration

Right ventricular (RV) dysfunction is a common cause of heart failure in patients with congenital heart defects and often leads to impaired functional capacity and premature death. Myocardial tissue regeneration techniques are being developed for the potential that viable myocardium may be regenerated to replace scar tissues in the heart or used as patch material in heart surgery. 3D computational RV/LV/Patch models with fluid-structure interactions (FSI) were constructed based on data from a healthy dog heart to obtain local fluid dynamics and structural stress/strain information and identify optimal conditions under which tissue regeneration techniques could achieve best outcome. RV/LV/Patch geometry and blood pressure data were obtained from a dog following established procedures. Four FSI models were used to quantify the influence of different patch materials (Dacron scaffold, treated pericardium) on local environment around the patch area, especially focusing on the thickness and stiffness of the patch. Our results indicated that changes in patch stiffness had little impact on the ejection fraction of the right ventricle because the total patch area was small. However, patch stiffness had huge impact on local RV maximum principal stress (Stress-P1) and strain (Strain-P1) around the patch area. Compared to the no-patch model, patch models had increased Stress-P1 and decreased Strain-P1 values in the patch area. Softer patches were associated with greater stress/strain variations. Thinner patch led to complex local flow environment which may have impact on myocytes seeding and RV remodeling. Our multi-physics RV/LV/Patch FSI model can serve as a useful tool to investigate cellular biology and tissue regeneration under localized flow and structural stress environment.     

Effect of a lipid pool on stress/strain distributions in stenotic arteries: 3-D fluid-structure interactions (FSI) models

Nonlinear 3-D models with fluid-structure interactions (FSI) based on in vitro experiments are introduced and solved by ADINA to perform flow and stress/strain analysis for stenotic arteries with lipid cores. Navier-Stokes equations are used as the governing equations for the fluid. Hyperelastic Mooney-Rivlin models are used for both the arteries and lipid cores. Our results indicate that critical plaque stress/strain conditions are affected considerably by stenosis severity, eccentricity, lipid pool size, shape and position, plaque cap thickness, axial stretch, pressure, and fluid-structure interactions, and may be used for possible plaque rupture predictions.     

3D reconstruction techniques of human coronary bifurcations for shear stress computations

Background

Heterogeneity in plaque composition in human coronary artery bifurcations is associated with blood flow induced shear stress. Shear stress is generally determined by combing 3D lumen data and computational fluid dynamics (CFD). We investigated two new procedures to generate 3D lumen reconstructions of coronary artery bifurcations for shear stress computations.

Methods

We imaged 10 patients with multislice computer tomography (MSCT) and intravascular ultrasound (IVUS). The 3D reconstruction of the main branch was based on the fusion of MSCT and IVUS. The proximal part of side branch was reconstructed using IVUS data or MSCT data, resulting in two different reconstructions of the bifurcation region. The distal part of the side branch was based on MSCT data alone. The reconstructed lumen was combined with CFD to determine the shear stress. Low and high shear stress regions were defined and shear stress patterns in the bifurcation regions were investigated.

Results

The 3D coronary bifurcations were successfully generated with both reconstruction procedures. The geometrical features of the bifurcation region for the two reconstruction procedures did not reveal appreciable differences. The shear stress maps showed a qualitative agreement, and the low and high shear stress regions were similar in size and average shear stress values were identical. The low and high shear stress regions showed an overlap of approximately 75%.

Conclusion

Reconstruction of the side branch with MSCT data alone is an adequate technique to study shear stress and wall thickness in the bifurcation region. The reconstruction procedure can be applied to further investigate the effect of shear stress on atherosclerosis in coronary bifurcations.     

The association of wall mechanics and morphology: a case study of abdominal aortic aneurysm growth

The purpose of this study is to evaluate the potential correlation between peak wall stress (PWS) and abdominal aortic aneurysm (AAA) morphology and how it relates to aneurysm rupture potential. Using in-house segmentation and meshing software, six 3-dimensional (3D) AAA models from a single patient followed for 28 months were generated for finite element analysis. For the AAA wall, both isotropic and anisotropic materials were used, while an isotropic material was used for the intraluminal thrombus (ILT). These models were also used to calculate 36 geometric indices characteristic of the aneurysm morphology. Using least squares regression, seven significant geometric features (p?相似文献   

Using in vivo Cine and 3D multi-contrast MRI to determine human atherosclerotic carotid artery material properties and circumferential shrinkage rate and their impact on stress/strain predictions

In vivo magnetic resonance image (MRI)-based computational models have been introduced to calculate atherosclerotic plaque stress and strain conditions for possible rupture predictions. However, patient-specific vessel material properties are lacking in those models, which affects the accuracy of their stress/strain predictions. A noninvasive approach of combining in vivo Cine MRI, multicontrast 3D MRI, and computational modeling was introduced to quantify patient-specific carotid artery material properties and the circumferential shrinkage rate between vessel in vivo and zero-pressure geometries. In vivo Cine and 3D multicontrast MRI carotid plaque data were acquired from 12 patients after informed consent. For each patient, one nearly-circular slice and an iterative procedure were used to quantify parameter values in the modified Mooney-Rivlin model for the vessel and the vessel circumferential shrinkage rate. A sample artery slice with and without a lipid core and three material parameter sets representing stiff, median, and soft materials from our patient data were used to demonstrate the effect of material stiffness and circumferential shrinkage process on stress/strain predictions. Parameter values of the Mooney-Rivlin models for the 12 patients were quantified. The effective Young's modulus (YM, unit: kPa) values varied from 137 (soft), 431 (median), to 1435 (stiff), and corresponding circumferential shrinkages were 32%, 12.6%, and 6%, respectively. Using the sample slice without the lipid core, the maximum plaque stress values (unit: kPa) from the soft and median materials were 153.3 and 96.2, which are 67.7% and 5% higher than that (91.4) from the stiff material, while the maximum plaque strain values from the soft and median materials were 0.71 and 0.293, which are about 700% and 230% higher than that (0.089) from the stiff material, respectively. Without circumferential shrinkages, the maximum plaque stress values (unit: kPa) from the soft, median, and stiff models were inflated to 330.7, 159.2, and 103.6, which were 116%, 65%, and 13% higher than those from models with proper shrinkage. The effective Young's modulus from the 12 human carotid arteries studied varied from 137 kPa to 1435 kPa. The vessel circumferential shrinkage to the zero-pressure condition varied from 6% to 32%. The inclusion of proper shrinkage in models based on in vivo geometry is necessary to avoid over-estimating the stresses and strains by up 100%. Material stiffness had a greater impact on strain (up to 700%) than on stress (up to 70%) predictions. Accurate patient-specific material properties and circumferential shrinkage could considerably improve the accuracy of in vivo MRI-based computational stress/strain predictions.