Fluid-structure interaction study on human arterial plaque with patient specific geometry and boundary conditions

Atherosclerotic plaque rupture has been extensively considered as the leading cause of death in western countries. It is believed that high stresses within plaque can be an important factor on triggering the rupture of the plaque. Stress analysis in the coronary and carotid arteries with plaque have been developed by many researchers from 2D to 3-D models, from structure analysis only to the Fluid-Structure Interaction (FSI) models[1].

Stress analysis on human arterial plaques by fluid structure interactions - Multi-case study

Atherosclerotic plaque rupture has been extensively considered as the leading cause of death in the world. It is believed that high stress within plaque can be an important factor which can trigger the rupture of the plaque. High resolution multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components (arterial wall, lipids, and fibrous cap) to be visualized in vivo [1]. The patient specific finite element model can be generated from the image data to perform stress analysis and provide critical information on understanding plaque rupture mechanisms [2]. The present work is to apply the procedure to a total of 14 patients (S1 ∼ S14), to study the stress distributions on carotid artery plaque reconstructed from multi-spectral magnetic resonance images, and the possible relationships between stress and plaque burdens.

Study of reproducibility of human arterial plaque reconstruction and its effects on stress analysis based on multispectral in vivo magnetic resonance imaging

Purpose: To quantify the uncertainties of carotid plaque morphology reconstruction based on patient-specific multispectral in vivo magnetic resonance imaging (MRI) and their impacts on the plaque stress analysis. Materials and Methods: In this study, three independent investigators were invited to reconstruct the carotid bifurcation with plaque based on MR images from two subjects to study the geometry reconstruction reproducibility. Finite element stress analyses were performed on the carotid bifurcations, as well as the models with artificially modified plaque geometries to mimic the image segmentation uncertainties, to study the impacts of the uncertainties to the stress prediction. Results: Plaque reconstruction reproducibility was generally high in the study. The uncertainties among interobservers are around one or the subpixel level. It also shows that the predicted stress is relatively less sensitive to the arterial wall segmentation uncertainties, and more affected by the accuracy of lipid region definition. For a model with lipid core region artificially increased by adding one pixel on the lipid region boundary, it will significantly increase the maximum Von Mises Stress in fibrous cap (>100%) compared with the baseline model for all subjects. Conclusion: The current in vivo MRI in the carotid plaque could provide useful and reliable information for plaque morphology. The accuracy of stress analysis based on plaque geometry is subject to MRI quality. The improved resolution/quality in plaque imaging with newly developed MRI protocols would generate more realistic stress predictions.

The impact of wall shear stress and pressure drop on the stability of the atherosclerotic plaque

Rupture of vulnerable atheromatous plaque in the carotid and coronary arteries often leads to stroke and heart attack respectively. The mechanism of blood flow and plaque rupture in stenotic arteries is still not fully understood. A three dimensional rigid wall model was solved under steady state conditions and unsteady conditions by assuming a time-varying inlet velocity profile to investigate the relative importance of axial forces and pressure drops in arteries with asymmetric stenosis. Flow-structure interactions were investigated for the same geometry and the results were compared with those retrieved with the corresponding 2D cross-section structural models. The Navier-Stokes equations were used as the governing equations for the fluid. The tube wall was assumed hyperelastic, homogeneous, isotropic and incompressible. The analysis showed that the three dimensional behavior of velocity, pressure and wall shear stress is in general very different from that predicted by cross-section models. Pressure drop across the stenosis was found to be much higher than shear stress. Therefore, pressure may be the more important mechanical trigger for plaque rupture other than shear stress, although shear stress is closely related to plaque formation and progression.

Does PGA external stenting reduce compliance mismatch in venous grafts?

Background: Autogenous vein grafting is widely used in regular bypassing procedures. Due to its mismatch with the host artery in both mechanical property and geometry, the graft often over expands under high arterial blood pressure and forms a step-depth where eddy flow develops, thus causing restenosis, fibrous graft wall, etc. External stents, such as sheaths being used to cuff the graft, have been introduced to eliminate these mismatches and increase the patency. Although histological and immunochemical studies have shown some positive effects of the external stent, the mechanical mismatch under the protection of an external stent remains poorly analyzed. Methods: In this study, the jugular veins taken from hypercholesterolemic rabbits were transplanted into the carotid arteries, and non-woven polyglycolic acid (PGA) fabric was used to fabricate the external stents to study the effect of the biodegradable external stent. Eight weeks after the operation, the grafts were harvested to perform mechanical tests and histological examinations. An arc tangent function was suggested to describe the relationship between pressure and cross-sectional area to analyse the compliance of the graft. Results: The results from the mechanical tests indicated that grafts either with or without external stents displayed large compliance in the low-pressure range and were almost inextensible in the high-pressure range. This was very different from the behavior of the arteries or veins in vivo. The data from histological tests showed that, with external stents, collagen fibers were more compact, whilst those in the graft without protection were looser and thicker. No elastic fiber was found in either kind of grafts. Furthermore, grafts without protection were over-expanded which resulted in much bigger cross-sectional areas. Conclusion: The PGA external extent contributes little to the reduction of the mechanical mismatch between the graft and its host artery while remodeling develops. For the geometric mismatch, it reduces the cross-section area, therefore matching with the host artery much better. Although there are some positive effects, conclusively the PGA is not an ideal material for external stent.

Structural analysis and magnetic resonance imaging predict plaque vulnerability: A study comparing symptomatic and asymptomatic individuals

Background: More than half of all cerebral ischemic events are the result of rupture of extracranial plaques. The clinical determination of carotid plaque vulnerability is currently based solely on luminal stenosis; however, it has been increasingly suggested that plaque morphology and biomechanical stress should also be considered. We used finite element analysis based on in vivo magnetic resonance imaging (MRI) to simulate the stress distributions within plaques of asymptomatic and symptomatic individuals. Methods: Thirty nonconsecutive subjects (15 symptomatic and 15 asymptomatic) underwent high-resolution multisequence in vivo MRI of the carotid bifurcation. Stress analysis was performed based on the geometry derived from in vivo MRI of the carotid artery at the point of maximal stenosis. The finite element analysis model considered plaque components to be hyperelastic. The peak stresses within the plaques of symptomatic and asymptomatic individuals were compared. Results: High stress concentrations were found at the shoulder regions of symptomatic plaques, and the maximal stresses predicted in this group were significantly higher than those in the asymptomatic group (508.2 ± 193.1 vs 269.6 ± 107.9 kPa; P = .004). Conclusions: Maximal predicted plaque stresses in symptomatic patients were higher than those predicted in asymptomatic patients by finite element analysis, suggesting the possibility that plaques with higher stresses may be more prone to be symptomatic and rupture. If further validated by large-scale longitudinal studies, biomechanical stress analysis based on high resolution in vivo MRI could potentially act as a useful tool for risk assessment of carotid atheroma. It may help in the identification of patients with asymptomatic carotid atheroma at greatest risk of developing symptoms or mild-to-moderate symptomatic stenoses, which currently fall outside current clinical guidelines for intervention.

Does calcium deposition play a role in the stability of atheroma? Location may be the key

Background: Rupture of vulnerable atheromatous plaque in the carotid and coronary arteries often leads to stroke and heart attack respectively. The role of calcium deposition and its contribution to plaque stability is controversial. This study uses both an idealized and a patient-specific model to evaluate the effect of a calcium deposit on the stress distribution within an atheromatous plaque. Methods: Using a finite-element method, structural analysis was performed on an idealized plaque model and the location of a calcium deposit within it was varied. In addition to the idealized model, in vivo high-resolution MR imaging was performed on 3 patients with carotid atheroma and stress distributions were generated. The individual plaques were chosen as they had calcium at varying locations with respect to the lumen and the fibrous cap. Results: The predicted maximum stress was increased by 47.5% when the calcium deposit was located in the thin fibrous cap in the model when compared with that in a model without a deposit. The result of adding a calcium deposit either to the lipid core or remote from the lumen resulted in almost no increase in maximal stress. Conclusion: Calcification at the thin fibrous cap may result in high stress concentrations, ultimately increasing the risk of plaque rupture. Assessing the location of calcification may, in the future, aid in the risk stratification of patients with carotid stenosis.

Simulation of the interaction between blood flow and atherosclerotic plaque

It has been well accepted that over 50% of cerebral ischemic events are the result of rupture of vulnerable carotid atheroma and subsequent thrombosis. Such strokes are potentially preventable by carotid interventions. Selection of patients for intervention is currently based on the severity of carotid luminal stenosis. It has been, however, widely accepted that luminal stenosis alone may not be an adequate predictor of risk. To evaluate the effects of degree of luminal stenosis and plaque morphology on plaque stability, we used a coupled nonlinear time-dependent model with flow-plaque interaction simulation to perform flow and stress/strain analysis for stenotic artery with a plaque. The Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the fluid. The Ogden strain energy function was used for both the fibrous cap and the lipid pool. The plaque Principal stresses and flow conditions were calculated for every case when varying the fibrous cap thickness from 0.1 to 2mm and the degree of luminal stenosis from 10% to 90%. Severe stenosis led to high flow velocities and high shear stresses, but a low or even negative pressure at the throat of the stenosis. Higher degree of stenosis and thinner fibrous cap led to larger plaque stresses, and a 50% decrease of fibrous cap thickness resulted in a 200% increase of maximum stress. This model suggests that fibrous cap thickness is critically related to plaque vulnerability and that, even within presence of moderate stenosis, may play an important role in the future risk stratification of those patients when identified in vivo using high resolution MR imaging.