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.     

Is arterial wall-strain stiffening an additional process responsible for atherosclerosis in coronary bifurcations?: an in vivo study based on dynamic CT and MRI

Coronary bifurcations represent specific regions of the arterial tree that are susceptible to atherosclerotic lesions. While the effects of vessel compliance, curvature, pulsatile blood flow, and cardiac motion on coronary endothelial shear stress have been widely explored, the effects of myocardial contraction on arterial wall stress/strain (WS/S) and vessel stiffness distributions remain unclear. Local increase of vessel stiffness resulting from wall-strain stiffening phenomenon (a local process due to the nonlinear mechanical properties of the arterial wall) may be critical in the development of atherosclerotic lesions. Therefore, the aim of this study was to quantify WS/S and stiffness in coronary bifurcations and to investigate correlations with plaque sites. Anatomic coronary geometry and cardiac motion were generated based on both computed tomography and MRI examinations of eight patients with minimal coronary disease. Computational structural analyses using the finite element method were subsequently performed, and spatial luminal arterial wall stretch (LW(Stretch)) and stiffness (LW(Stiff)) distributions in the left main coronary bifurcations were calculated. Our results show that all plaque sites were concomitantly subject to high LW(Stretch) and high LW(Stiff), with mean amplitudes of 34.7 ± 1.6% and 442.4 ± 113.0 kPa, respectively. The mean LW(Stiff) amplitude was found slightly greater at the plaque sites on the left main coronary artery (mean value: 482.2 ± 88.1 kPa) compared with those computed on the left anterior descending and left circumflex coronary arteries (416.3 ± 61.5 and 428.7 ± 181.8 kPa, respectively). These findings suggest that local wall stiffness plays a role in the initiation of atherosclerotic lesions.     

Biaxial elastic material properties of porcine coronary media and adventitia

The importance of mechanical stresses and strains has become well recognized in vascular physiology and pathology. To compute the stress and strain on the various components of the vessel wall, we must know the constitutive equations for the different layers of the vessel wall. The objective of the present study is to determine the constitutive equation of the coronary artery treated as a two-layer composite: intima-media and adventitial layers. Twelve hearts were obtained from a local slaughterhouse, and the right coronary artery and left anterior descending artery were dissected free from the myocardium. The vessel wall was initially mechanically tested biaxially (inflation and axial extension) as a whole (intact wall) and subsequently as intima-media or adventitial layer. A Fung-type exponential strain energy function was used to curve fit the experimental data for the intact wall and individual layers for the right coronary artery and left anterior descending artery. Two methods were used for the determination of material constants, including the Marquardt-Levenberg nonlinear least squares method and the genetic algorithm method. Our results show that there were no statistically significant differences in the material constants obtained from the two methods and that either set of elastic constants results in good fit of the data. Furthermore, at an in vivo value of axial stretch ratio, we find that the stiffness is as follows: intima-media > intact > adventitia. These results underscore the composite nature of coronary arteries with different material properties in each layer. The present results are necessary for analysis of coronary artery mechanics and to provide a fundamental understanding of vessel physiology.     

Site specific inelasticity of arterial tissue

Understanding the mechanical behaviour of arterial tissue is vital to the development and analysis of medical devices targeting diseased vessels. During angioplasty and stenting, stress softening and permanent deformation of the vessel wall occur during implantation of the device, however little data exists on the inelastic behaviour of cardiovascular tissue and how this varies through the arterial tree. The aim of this study was to characterise the magnitude of stress softening and inelastic deformations due to loading throughout the arterial tree and to investigate the anisotropic inelastic behaviour of the tissue. Cyclic compression tests were used to investigate the differences in inelastic behaviour for carotid, aorta, femoral and coronary arteries harvested from 3-4 month old female pigs, while the anisotropic behaviour of aortic and carotid tissue was determined using cyclic tensile tests in the longitudinal and circumferential directions. The differences in inelastic behaviour were correlated to the ratio of collagen to elastin content of the arteries. It was found that larger inelastic deformations occurred in muscular arteries (coronary), which had a higher collagen to elastin ratio than elastic arteries (aorta), where the smallest inelastic deformations were observed. Lower magnitude inelastic deformations were observed in the circumferential tensile direction than in the longitudinal tensile direction or due to radial compression. This may be as a result of non-collagenous components in the artery becoming more easily damaged than the collagen fibres during loading. Stress softening was also found to be dependent on artery type. In the future, computational models should consider such site dependant, anisotropic inelastic behaviour in order to better predict the outcomes of interventional procedures such as angioplasty and stenting.     

Geometry guided data averaging enables the interpretation of shear stress related plaque development in human coronary arteries

The average low shear stress (SS) is known to determine predilection sites of atherosclerotic plaques. However, as plaques encroach into the lumen and thereby increase SS, interpretation of patient-specific data obtained at one moment in time regarding the influence of SS in the generation of atherosclerosis is not straightforward. This study aims to compare two methods of data analysis for the aid of data interpretation: (a) point-wise analysis of the raw data, (b) global analysis: to assess the history related natural SS distribution in coronary arteries by averaging the data in the axial vessel direction. Normal to mildly diseased human coronary arteries were investigated applying a combination of 3-D reconstruction technique and computational fluid dynamics (CFD). Point-wise analysis relating local wall thickness to local SS showed in only 4% of the cases an inverse relationship. In contrast, averaging the data in the axial vessel direction, showed in 38% a significant inverse relation between wall thickness and SS, resulting in an average negative slope of -0.70+/-0.46 mm/Pa. These data suggest that using a geometry guided way of data averaging may reveal history related effects of SS, which in part explains localization of atherosclerotic plaques.     

Experimental determination of circumferential properties of fresh carotid artery plaques

Carotid endarterectomy (CEA) is currently accepted as the gold standard for interventional revascularisation of diseased arteries belonging to the carotid bifurcation. Despite the proven efficacy of CEA, great interest has been generated in carotid angioplasty and stenting (CAS) as an alternative to open surgical therapy. CAS is less invasive compared with CEA, and has the potential to successfully treat lesions close to the aortic arch or distal internal carotid artery (ICA). Following promising results from two recent trials (CREST; Carotid revascularisation endarterectomy versus stenting trial, and ICSS; International carotid stenting study) it is envisaged that there will be a greater uptake in carotid stenting, especially amongst the group who do not qualify for open surgical repair, thus creating pressure to develop computational models that describe a multitude of plaque models in the carotid arteries and their reaction to the deployment of such interventional devices. Pertinent analyses will require fresh human atherosclerotic plaque material characteristics for different disease types. This study analysed atherosclerotic plaque characteristics from 18 patients tested on site, post-surgical revascularisation through endarterectomy, with 4 tissue samples being excluded from tensile testing based on large width-length ratios. According to their mechanical behaviour, atherosclerotic plaques were separated into 3 grades of stiffness. Individual and group material coefficients were then generated analytically using the Yeoh strain energy function. The ultimate tensile strength (UTS) of each sample was also recorded, showing large variation across the 14 atherosclerotic samples tested. Experimental Green strains at rupture varied from 0.299 to 0.588 and the Cauchy stress observed in the experiments was between 0.131 and 0.779 MPa. It is expected that this data may be used in future design optimisation of next generation interventional medical devices for the treatment and revascularisation of diseased arteries of the carotid bifurcation.     

Blood flow patterns in the proximal human coronary arteries: relationship to atherosclerotic plaque occurrence

Atherosclerotic plaques in human coronary arteries are focal manifestations of systemic disease, and biomechanical factors have been hypothesized to contribute to plaque genesis and localization. We developed a computational fluid dynamics (CFD) model of the ascending aorta and proximal sections of the right and left coronary arteries of a normal human subject using computed tomography (CT) and magnetic resonance imaging (MRI) and determined the pulsatile flow field. Results demonstrate that flow patterns in the ascending aorta contribute to a pro-atherosclerotic flow environment, specifically through localization of low and oscillatory wall shear stress in the neighborhood of coronary orifices. Furthermore, these patterns differ in their spatial distribution between right and left coronary arteries. Entrance effects of aortic flow diminish within two vessel diameters. We examined relationships between spatial distributions of wall shear stress and reports of plaque occurrence in the literature. Results indicate low wall shear stress is co-located with increased incidence of lesions, and higher wall shear stresses are associated with lesion-resistant areas. This investigation does not consider plaque progression or advanced lesions, inasmuch as the CFD model was developed from a normal individual and the clinical data used for comparisons were obtained from autopsy specimens of subjects who died from non-cardiovascular causes. The data reported are consistent with the hypothesis that low wall shear stress is associated with the localization of atherosclerotic lesions, and the results demonstrate the importance of aortic flow on flow patterns in the proximal segments of the coronary arteries.     

An experimental study on the ultimate strength of the adventitia and media of human atherosclerotic carotid arteries in circumferential and axial directions

Atherosclerotic plaque may rupture without warning causing heart attack or stroke. Knowledge of the ultimate strength of human atherosclerotic tissues is essential for understanding the rupture mechanism and predicting cardiovascular events. Despite its great importance, experimental data on ultimate strength of human atherosclerotic carotid artery remains very sparse. This study determined the uniaxial tensile strength of human carotid artery sections containing type II and III lesions (AHA classifications). Axial and circumferential oriented adventitia, media and intact specimens (total=73) were prepared from 6 arteries. The ultimate strength in uniaxial tension was taken as the peak stress recorded when the specimen showed the first evidence of failure and the extensibility was taken as the stretch ratio at failure. The mean adventitia strength values calculated using the first Piola–Kirchoff stress were 1996±867 and 1802±703 kPa in the axial and circumferential directions respectively, while the corresponding values for the media sections were 519±270 and 1230±533 kPa. The intact specimens showed ultimate strengths similar to media in circumferential direction but were twice as strong as the media in the axial direction. Results also indicated that adventitia, media and intact specimens exhibited similar extensibility at failure, in both the axial and circumferential directions (stretch ratio 1.50±0.22). These measurements of the material strength limits for human atherosclerotic carotid arteries could be useful in improving computational models that assess plaque vulnerability.     

A theoretical and non-destructive experimental approach for direct inclusion of measured collagen orientation and recruitment into mechanical models of the artery wall

Gradual collagen recruitment has been hypothesized as the underlying mechanism for the mechanical stiffening with increasing stress in arteries. In this work, we investigated this hypothesis in eight rabbit carotid arteries by directly measuring the distribution of collagen recruitment stretch under increasing circumferential loading using a custom uniaxial (UA) extension device combined with a multi-photon microscope (MPM). This approach allowed simultaneous mechanical testing and imaging of collagen fibers without traditional destructive fixation methods. Fiber recruitment was quantified from 3D rendered MPM images, and fiber orientation was measured in projected stacks of images. Collagen recruitment was observed to initiate at a finite strain, corresponding to a sharp increase in the measured mechanical stiffness, confirming the previous hypothesis and motivating the development of a new constitutive model to capture this response. Previous constitutive equations for the arterial wall have modeled the collagen contribution with either abrupt recruitment at zero strain, abrupt recruitment at finite strain or as gradual recruitment beginning at infinitesimal strain. Based on our experimental data, a new combined constitutive model was presented in which fiber recruitment begins at a finite strain with activation stretch represented by a probability distribution function. By directly including this recruitment data, the collagen contribution was modeled using a simple Neo-Hookean equation. As a result, only two phenomenological material constants were required from the fit to the stress stretch data. Three other models for the arterial wall were then compared with these results. The approach taken here was successful in combining stress-strain analysis with simultaneous microstructural imaging of collagen recruitment and orientation, providing a new approach by which underlying fiber architecture may be quantified and included in constitutive equations.     

Tissue Doppler imaging of carotid plaque wall motion: a pilot study

Background

Studies suggest the physical and mechanical properties of vessel walls and plaque may be of clinical value in the diagnosis and treatment of cardiovascular atherosclerotic disease. The purpose of this pilot study was to investigate the potential clinical application of ultrasound Tissue Doppler Imaging (TDI) of Arterial Wall Motion (AWM) and to quantify simple wall motion indices in normal and diseased carotid arteries.

Methods

224 normal and diseased carotid arteries (0–100% stenoses) were imaged in 126 patients (age 25–88 years, mean 68 ± 11). Longitudinal sections of the carotid bifurcation were imaged using a Philips HDI5000 scanner and L12-5 probe under optimized TDI settings. Temporal and spatial AWMs were analyzed to evaluate the vessel wall displacements and spatial gradients at peak systole averaged over 5 cardiac cycles.

Results

AWM data were successfully extracted in 91% of cases. Within the carotid bifurcation/plaque region, the maximum wall dilation at peak systole ranged from -100 to 750 microns, mean 335 ± 138 microns. Maximum wall dilation spatial gradients ranged 0–0.49, mean 0.14 ± 0.08. The AWM parameters showed a wide variation and had poor correlation with stenoses severity. Case studies illustrated a variety of pertinent qualitative and quantitative wall motion features related to the biophysics of arterial disease.

Conclusion

Our clinical experience, using a challenging but realistic imaging protocol, suggests the use of simple quantitative AWM measures may have limitations due to high variability. Despite this, pertinent features of AWM in normal and diseased arteries demonstrate the potential clinical benefit of the biomechanical information provided by TDI.     

Atherosclerosis and flow in carotid arteries with authentic geometries

The influence of blood flow on the depositions and development of atherosclerotic lesions have been observed and described since the 19th century. Observations have shown that depositions correlate with regions of low wall shear stress. However, the exact correlations between depositions, vessel geometry and flow parameters are not yet known. The purpose of this study was the quantification of atherosclerosis risk factors in carotid bifurcation. This artery has attracted particular interest because lesions are often found in this bifurcation. Post mortem, the arteries are excised and vessel casts are produced. Afterwards, the arteries are analyzed morphometrically. The vessel casts are used for the assessment of some geometrical parameters. 31 carotid bifurcations were analyzed in this study. Eight vessel casts were digitized and rendered three-dimensional mathematical models of the arteries. These data were imported by the computational fluid dynamics program FLUENT. Further, the blood flow was reconstructed in a computer model based on the individual vessel geometry. The flow parameters, such as velocity, pressure and wall shear stress were computed. At the same time the geometrical parameters and wall alterations are known. This permits the comparison of the anatomical shape and its flow with the distribution and level of the wall alterations.     

Estrogen modulates the mechanical homeostasis of mouse arterial vessels through nitric oxide

We have recently shown that estrogen causes vessel dilation through receptor-mediated stimulation of nitric oxide (NO) production. Here, we hypothesize that estrogen modulates the mechanical homeostasis in the blood vessel wall through NO production. The mechanical properties of female ovariectomized (ovx) mice, female mice lacking the gene for endothelial NO synthase (eNOS(-/-)), and control female and male mice were studied to test the hypothesis. The femoral and carotid arteries and aorta were cannulated in situ and mechanically distended. The stress, strain, elastic modulus, and wall thickness of vessels in ovx and eNOS(-/-) mice, as well as intact female and male mice, were determined. Western blot and immunohistochemistry were used to assess eNOS protein expression in the aorta. Moreover, NO by-products of the femoral and carotid artery were determined by measuring the levels of nitrite and nitrate. Our results show that ovariectomy and eNOS(-/-) significantly decrease the strain in all arteries. Furthermore, the eNOS protein was significantly reduced in ovx mice. Finally, the NO metabolites were significantly decreased both in ovx and eNOS(-/-) mice. We found statistically significant correlations between the structural (wall thickness), mechanical (stress, strain, and elastic modulus), and biochemical parameters (NO by-products). These novel results connect NO to the structural and mechanical properties of the vessel wall. Hence, the effect of endogenous estrogen on the arterial mechanical properties is mediated by the regulation of NO derived from eNOS.     

Assessment of endothelial function of large, medium, and small vessels: a unified myograph

Endothelial dysfunction precedes the development of morphological atherosclerotic changes and can also contribute to lesion development in cardiovascular diseases. Currently, there is a lack of a single method to determine endothelial function of the entire range of vessel dimensions from aorta to arterioles. Here we assessed endothelial function of a large range of size arteries using a unified isovolumic myograph method. The method maintains a constant volume of fluid in the lumen of the vessel during contraction and relaxation, which are characterized by an increase and a decrease of pressure, respectively. Segments of six aortas, six common femoral arteries, and six mesenteric arteries from rats; six carotid arteries from mice; and six coronary and carotid arteries from pigs were used. The endothelium-dependent dose-response vasorelaxation was determined with endothelium-dependent vasodilators while arterial preconstriction was induced with vasoconstrictors at a submaximal dose. The circumferential midtension during vascular reactivity varied from 43.1 ± 7.9 to 2.59 ± 0.46 mN/mm (from large to small arteries), whereas the circumferential midstress showed a much smaller variation from 217 ± 23.5 to 123 ± 15.3 kPa (in the same range of vessels). We also found that overinflation and axial overelongation compromised endothelium-dependent vasorelaxation to underscore the significance of vessel preload. In conclusion, an isovolumic myograph was used to unify arterial vasoreactivity from large to small arteries and shows the uniformity of wall stress and %tension throughout the range of vessel sizes.     

Blood flow and vessel mechanics in a physiologically realistic model of a human carotid arterial bifurcation

The pulsatile flow in an anatomically realistic compliant human carotid bifurcation was simulated numerically. Pressure and mass flow waveforms in the carotid arteries were obtained from an individual subject using non-invasive techniques. The geometry of the computational model was reconstructed from magnetic resonance angiograms. Maps of time-average wall shear stress, contours of velocity in the flow field as well as wall movement and tensile stress on the arterial wall are all presented. Inconsistent with previous findings from idealised geometry models, flow in the carotid sinus is dominated by a strong helical flow accompanied by a single secondary vortex motion. This type of flow is induced primarily by the asymmetry and curvature of the in vivo geometry. Flow simulations have been carried out under the rigid wall assumption and for the compliant wall, respectively. Comparison of the results demonstrates the quantitative influence of the vessel wall motion. Generally there is a reduction in the magnitude of wall shear stress, with its degree depending on location and phase of the cardiac cycle. The region of slow or reversed flow was greater, in both spatial and temporal terms in the compliant model, but the global characteristics of the flow and stress patterns remain unchanged. The analysis of mechanical stresses on the vessel surface shows a complicated stress field. Stress concentration occurs at both the anterior and posterior aspects of the proximal internal bulb. These are also regions of low wall shear stress. The comparison of computed and measured wall movement generally shows good agreement.     

Effect of cigarette smoking on nitric oxide, structural, and mechanical properties of mouse arteries

Cigarette smoking (CS) is a major risk factor for vascular disease. The aim of this study was to quantitatively assess the influence of CS on mouse arteries. We studied the effect of short-term (6 wk) and long-term (16 wk) CS exposure on structural and mechanical properties of coronary arteries compared with that of control mice. We also examined the reversibility of the deleterious effects of CS on structural [e.g., wall thickness (WT)], mechanical (e.g., stiffness), and biochemical [e.g., nitric oxide (NO) by-products] properties with the cessation of CS. The left and right coronary arteries were cannulated in situ and mechanically distended. The stress, strain, elastic modulus, and WT of coronary arteries were determined. Western blot analysis was used to analyze endothelial NO synthase (eNOS) in the femoral and carotid arteries of the same mice, and NO by-products were determined by measuring the levels of nitrite. Our results show that the mean arterial pressure was increased by CS. Furthermore, CS significantly increased the elastic modulus, decreased stress and strain, and increased the WT and WT-to-radius ratio compared with those of control mice. The reduction of eNOS protein expression was found only after long-term CS exposure. Moreover, the NO metabolite was markedly decreased in CS mice after short- and long-term exposure of CS. These findings suggest that 16 wk of CS exposure can cause an irreversible deterioration of structural and elastic properties of mouse coronary arteries. The decrease in endothelium-derived NO in CS mice was seen to significantly correlate with the remodeling of arterial wall.     

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     

Determination of layer-specific mechanical properties of human coronary arteries with nonatherosclerotic intimal thickening and related constitutive modeling

At autopsy, 13 nonstenotic human left anterior descending coronary arteries [71.5 +/- 7.3 (mean +/- SD) yr old] were harvested, and related anamnesis was documented. Preconditioned prepared strips (n = 78) of segments from the midregion of the left anterior descending coronary artery from the individual layers in axial and circumferential directions were subjected to cyclic quasi-static uniaxial tension tests, and ultimate tensile stresses and stretches were documented. The ratio of outer diameter to total wall thickness was 0.189 +/- 0.014; ratios of adventitia, media, and intima thickness to total wall thickness were 0.4 +/- 0.03, 0.36 +/- 0.03, and 0.27 +/- 0.02, respectively; axial in situ stretch of 1.044 +/- 0.06 decreased with age. Stress-stretch responses for the individual tissues showed pronounced mechanical heterogeneity. The intima is the stiffest layer over the whole deformation domain, whereas the media in the longitudinal direction is the softest. All specimens exhibited small hysteresis and anisotropic and strong nonlinear behavior in both loading directions. The media and intima showed similar ultimate tensile stresses, which are on average three times smaller than ultimate tensile stresses in the adventitia (1,430 +/- 604 kPa circumferential and 1,300 +/- 692 kPa longitudinal). The ultimate tensile stretches are similar for all tissue layers. A recently proposed constitutive model was extended and used to represent the deformation behavior for each tissue type over the entire loading range. The study showed the need to model nonstenotic human coronary arteries with nonatherosclerotic intimal thickening as a composite structure composed of three solid mechanically relevant layers with different mechanical properties. The intima showed significant thickness, load-bearing capacity, and mechanical strength compared with the media and adventitia.     

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

Mounting evidence suggests that the pulsatile character of blood pressure and flow within large arteries plays a particularly important role as a mechano-biological stimulus for wall growth and remodeling. Nevertheless, understanding better the highly coupled interactions between evolving wall geometry, structure, and properties and the hemodynamics will require significantly more experimental data. Computational fluid–solid-growth models promise to aid in the design and interpretation of such experiments and to identify candidate mechanobiological mechanisms for the observed arterial adaptations. Motivated by recent aortic coarctation models in animals, we used a computational fluid–solid interaction model to study possible local and systemic effects on the hemodynamics within the thoracic aorta and coronary, carotid, and cerebral arteries due to a distal aortic coarctation and subsequent spatial variations in wall adaptation. In particular, we studied an initial stage of acute cardiac compensation (i.e., maintenance of cardiac output) followed by early arterial wall remodeling (i.e., spatially varying wall thickening and stiffening). Results suggested, for example, that while coarctation increased both the mean and pulse pressure in the proximal vessels, the locations nearest to the coarctation experienced the greatest changes in pulse pressure. In addition, after introducing a spatially varying wall adaptation, pressure, left ventricular work, and wave speed all increased. Finally, vessel wall strain similarly experienced spatial variations consistent with the degree of vascular wall adaptation.     

An efficient two-stage approach for image-based FSI analysis of atherosclerotic arteries

Patient-specific biomechanical modeling of atherosclerotic arteries has the potential to aid clinicians in characterizing lesions and determining optimal treatment plans. To attain high levels of accuracy, recent models use medical imaging data to determine plaque component boundaries in three dimensions, and fluid–structure interaction is used to capture mechanical loading of the diseased vessel. As the plaque components and vessel wall are often highly complex in shape, constructing a suitable structured computational mesh is very challenging and can require a great deal of time. Models based on unstructured computational meshes require relatively less time to construct and are capable of accurately representing plaque components in three dimensions. These models unfortunately require additional computational resources and computing time for accurate and meaningful results. A two-stage modeling strategy based on unstructured computational meshes is proposed to achieve a reasonable balance between meshing difficulty and computational resource and time demand. In this method, a coarsegrained simulation of the full arterial domain is used to guide and constrain a fine-scale simulation of a smaller region of interest within the full domain. Results for a patient-specific carotid bifurcation model demonstrate that the two-stage approach can afford a large savings in both time for mesh generation and time and resources needed for computation. The effects of solid and fluid domain truncation were explored, and were shown to minimally affect accuracy of the stress fields predicted with the two-stage approach.