Fluid-structure interaction modeling of aneurysmal arteries under steady-state and pulsatile blood flow: a stability analysis

Tortuous aneurysmal arteries are often associated with a higher risk of rupture but the mechanism remains unclear. The goal of this study was to analyze the buckling and post-buckling behaviors of aneurysmal arteries under pulsatile flow. To accomplish this goal, we analyzed the buckling behavior of model carotid and abdominal aorta with aneurysms by utilizing fluid-structure interaction (FSI) method with realistic waveforms boundary conditions. FSI simulations were done under steady-state and pulsatile flow for normal (1.5) and reduced (1.3) axial stretch ratios to investigate the influence of aneurysm, pulsatile lumen pressure and axial tension on stability. Our results indicated that aneurysmal artery buckled at the critical buckling pressure and its deflection nonlinearly increased with increasing lumen pressure. Buckling elevates the peak stress (up to 118%). The maximum aneurysm wall stress at pulsatile FSI flow was (29%) higher than under static pressure at the peak lumen pressure of 130 mmHg. Buckling results show an increase in lumen shear stress at the inner side of the maximum deflection. Vortex flow was dramatically enlarged with increasing lumen pressure and artery diameter. Aneurysmal arteries are more susceptible than normal arteries to mechanical instability which causes high stresses in the aneurysm wall that could lead to aneurysm rupture.     

Effects of elastin degradation and surrounding matrix support on artery stability

Tortuous arteries are often associated with aging, hypertension, atherosclerosis, and degenerative vascular diseases, but the mechanisms are poorly understood. Our recent theoretical analysis suggested that mechanical instability (buckling) may lead to tortuous blood vessels. The objectives of this study were to determine the critical pressure of artery buckling and the effects of elastin degradation and surrounding matrix support on the mechanical stability of arteries. The mechanical properties and critical buckling pressures, at which arteries become unstable and deform into tortuous shapes, were determined for a group of five normal arteries using pressurized inflation and buckling tests. Another group of nine porcine arteries were treated with elastase (8 U/ml), and the mechanical stiffness and critical pressure were obtained before and after treatment. The effect of surrounding tissue support was simulated using a gelatin gel. The critical pressures of the five normal arteries were 9.52 kPa (SD 1.53) and 17.10 kPa (SD 5.11) at axial stretch ratios of 1.3 and 1.5, respectively, while model predicted critical pressures were 10.11 kPa (SD 3.12) and 17.86 kPa (SD 5.21), respectively. Elastase treatment significantly reduced the critical buckling pressure (P < 0.01). Arteries with surrounding matrix support buckled into multiple waves at a higher critical pressure. We concluded that artery buckling under luminal pressure can be predicted by a buckling equation. Elastin degradation weakens the arterial wall and reduces the critical pressure, which thus leads to tortuous vessels. These results shed light on the mechanisms of the development of tortuous vessels due to elastin deficiency.     

Mechanical instability of normal and aneurysmal arteries

Tortuous arteries associated with aneurysms have been observed in aged patients with atherosclerosis and hypertension. However, the underlying mechanism is poorly understood. The objective of this study was to determine the effect of aneurysms on arterial buckling instability and the effect of buckling on aneurysm wall stress. We investigated the mechanical buckling and post-buckling behavior of normal and aneurysmal carotid arteries and aorta’s using computational simulations and experimental measurements to elucidate the interrelationship between artery buckling and aneurysms. Buckling tests were done in porcine carotid arteries with small aneurysms created using elastase treatment. Parametric studies were done for model aneurysms with orthotropic nonlinear elastic walls using finite element simulations. Our results demonstrated that arteries buckled at a critical buckling pressure and the post-buckling deflection increased nonlinearly with increasing pressure. The presence of an aneurysm can reduce the critical buckling pressure of arteries, although the effect depends on the aneurysm’s dimensions. Buckled aneurysms demonstrated a higher peak wall stress compared to unbuckled aneurysms under the same lumen pressure. We conclude that aneurysmal arteries are vulnerable to mechanical buckling and mechanical buckling could lead to high stresses in the aneurysm wall. Buckling could be a possible mechanism for the development of tortuous aneurysmal arteries such as in the Loeys–Dietz syndrome.     

Artery buckling analysis using a four-fiber wall model

Artery bent buckling has been suggested as a possible mechanism that leads to artery tortuosity, which is associated with aging, hypertension, atherosclerosis, and other pathological conditions. It is necessary to understand the relationship between microscopic wall structural changes and macroscopic artery buckling behavior. To this end, the objectives of this study were to develop arterial buckling equations using a microstructure-based 4-fiber reinforced wall model, and to simulate the effects of vessel wall microstructural changes on artery buckling. Our results showed that the critical pressure increased nonlinearly with the axial stretch ratio, and the 4-fiber model predicted higher critical buckling pressures than what the Fung model predicted. The buckling equation using the 4-fiber model captured the experimentally observed reduction of critical pressure induced by elastin degradation and collagen fiber orientation changes in the arterial wall. These results improve our understanding of arterial stability and its relationship to microscopic wall remodeling, and the model provides a useful tool for further studies.     

Onset of pulsatile pressure causes transiently increased filtration through artery wall

Convective fluid motion through artery walls aids in the transvascular transport of macromolecules. Although many measurements of convective filtration have been reported, they were all obtained under constant transmural pressure. However, arterial pressure in vivo is pulsatile. Therefore, experiments were designed to compare filtration under steady and pulsatile pressure conditions. Rabbit carotid arteries were cannulated and excised from male New Zealand White rabbits anesthetized with pentobarbitol sodium (30 mg/kg i.v. administered). Hydraulic conductance was measured in cannulated excised rabbit carotid arteries at steady pressure. Next, pulsatile pressure trains were applied within the same vessels, and, simultaneously, arterial distension was monitored using Optical coherence tomography (OCT). For each pulse train, the volume of fluid lost through filtration was measured (subtracting volume change due to residual distension) and compared with that predicted from steady pressure measurements. At 60- and 80-mmHg baseline pressures, the experimental filtration volumes were significantly increased compared with those predicted for steady pressure (P < 0.05). OCT demonstrated that the excess fluid volume loss was significantly greater than the volume that would be lost through residual distension (P < 0.05). After 30 s, the magnitude of the excess of fluid loss was reduced. These results suggest that sudden onset of pulsatile pressure may cause changes in arterial interstitial hydration.     

The effect of collagenase on the critical buckling pressure of arteries

The stability of arteries is essential to normal arterial functions and loss of stability can lead to arterial tortuosity and kinking. Collagen is a main extracellular matrix component that modulates the mechanical properties of arteries and collagen degradation at pathological conditions weakens the mechanical strength of arteries. However, the effects of collagen degradation on the mechanical stability of arteries are unclear. The objective of this study was to investigate the effects of collagen degradation on the critical buckling pressure of arteries. Arterial specimens were subjected to pressurized inflation testing and fitted with nonlinear thick-walled cylindrical model equations to determine their stress strain relationships. The arteries were then tested for the critical buckling pressure at a set of axial stretch ratios. Then, arteries were divided into three groups and treated with Type III collagenase at three different concentrations (64, 128, and 400 U/ml). Mechanical properties and buckling pressures of the arteries were determined after collagenase treatment. Additionally, the theoretical buckling pressures were also determined using a buckling equation. Our results demonstrated that the buckling pressure of arteries was lower after collagenase treatment. The difference between pre- and post- treatment was statistically significant for the highest concentration of 400U/ml but not at the lower concentrations. The buckling equation was found to yield a fair estimation to the experimental critical pressure measurements. These results shed light on the role of matrix remodeling on the mechanical stability of arteries and developments of tortuous arteries.     

The Effect of Collagenase on the Critical Buckling Pressure of Arteries^*

The stability of arteries is essential to normal arterial functions and loss of stability can lead to arterial tortuosity and kinking. Collagen is a main extracellular matrix component that modulates the mechanical properties of arteries and collagen degradation at pathological conditions weakens the mechanical strength of arteries. However, the effects of collagen degradation on the mechanical stability of arteries are unclear. The objective of this study was to investigate the effects of collagen degradation on the critical buckling pressure of arteries. Arterial specimens were subjected to pressurized inflation testing and fitted with nonlinear thick-walled cylindrical model equations to determine their stress strain relationships. The arteries were then tested for the critical buckling pressure at a set of axial stretch ratios. Then, arteries were divided into three groups and treated with Type III collagenase at three different concentrations (64, 128, and 400U/ml). Mechanical properties and buckling pressures of the arteries were determined after collagenase treatment. Additionally, the theoretical buckling pressures were also determined using a buckling equation. Our results demonstrated that the buckling pressure of arteries was lower after collagenase treatment. The difference between pre- and post- treatment was statistically significant for the highest concentration of 400U/ml but not at the lower concentrations. The buckling equation was found to yield a fair estimation to the experimental critical pressure measurements. These results shed light on the role of matrix remodeling on the mechanical stability of arteries and developments of tortuous arteries.     

Nonlinear buckling of blood vessels: a theoretical study

Tortuosity and kinking often occur in arteries and veins but the underlying mechanisms are poorly understood. It has been suggested recently that long arteries may buckle and become tortuosity due to reduced axial tension or hypertensive pressure, but very few studies have been done to establish the biomechanical basis for artery buckling. Here we developed the arterial buckling equation using a nonlinear elastic thick-walled cylindrical model with residual stress. Our results demonstrated that arteries may buckle due to high blood pressure or low axial tension and that residual stress in the arteries increases the buckling pressure. These results are in general agreement with the previous linear elastic model. The buckling equation provides a useful tool for studying artery tortuosity and kinking.     

Steady flow and wall compression in stenotic arteries: a three-dimensional thick-wall model with fluid-wall interactions.

Severe stenosis may cause critical flow and wall mechanical conditions related to artery fatigue, artery compression, and plaque rupture, which leads directly to heart attack and stroke. The exact mechanism involved is not well understood. In this paper a nonlinear three-dimensional thick-wall model with fluid-wall interactions is introduced to simulate blood flow in carotid arteries with stenosis and to quantify physiological conditions under which wall compression or even collapse may occur. The mechanical properties of the tube wall were selected to match a thick-wall stenosis model made of PVA hydrogel. The experimentally measured nonlinear stress-strain relationship is implemented in the computational model using an incremental linear elasticity approach. The Navier-Stokes equations are used for the fluid model. An incremental boundary iteration method is used to handle the fluid-wall interactions. Our results indicate that severe stenosis causes considerable compressive stress in the tube wall and critical flow conditions such as negative pressure, high shear stress, and flow separation which may be related to artery compression, plaque cap rupture, platelet activation, and thrombus formation. The stress distribution has a very localized pattern and both maximum tensile stress (five times higher than normal average stress) and maximum compressive stress occur inside the stenotic section. Wall deformation, flow rates, and true severities of the stenosis under different pressure conditions are calculated and compared with experimental measurements and reasonable agreement is found.     

Carotid distensibility characterized via the isometric exercise pressor response

Distensibility of the large elastic arteries is a key index for cardiovascular health. Distensibility, usually estimated from resting values in humans, is not a static characteristic but a negative curvilinear function of pressure. We hypothesized that differences in vascular function with gender and age may only be recognized if distensibility is quantified over a range of pressures. We used isometric handgrip exercise to induce progressive increases in pressures and carotid diameters, thereby enhancing the characterization of distensibility. In 30 volunteers, evenly distributed by gender and age across the third to fifth decades of life, we derived pulsatile distensibility slopes as a function of arterial pressure for a dynamic distensibility index and compared it with a traditional static index at a reference pressure of 95 mmHg. We also assessed intima-media thickness (IMT). We found that women had greater distensibility slopes within each decade, despite comparable IMT. Furthermore, declines in distensibility slope with increasing age were correlated to increased IMT. The static distensibility index failed to show gender-related differences in distensibility but did show age-related differences. Our results indicate that gender- and age-related differences can be manifest even in young, healthy adults and may only be identified with techniques that assess carotid distensibility across a range of pressures.     

Modelling the circle of Willis to assess the effects of anatomical variations and occlusions on cerebral flows

Blood flow in the circle of Willis (CoW) is modelled using the 1-D equations of pressure and flow wave propagation in compliant vessels. The model starts at the left ventricle and includes the largest arteries that supply the CoW. Based on published physiological data, it is able to capture the main features of pulse wave propagation along the aorta, at the brachiocephalic bifurcation and throughout the cerebral arteries. The collateral ability of the complete CoW and its most frequent anatomical variations is studied in normal conditions and after occlusion of a carotid or vertebral artery (VA). Our results suggest that the system does not require collateral pathways through the communicating arteries to adequately perfuse the brain of normal subjects. The communicating arteries become important in cases of missing or occluded vessels, the anterior communicating artery (ACoA) being a more critical collateral pathway than the posterior communicating arteries (PCoAs) if an internal carotid artery (ICA) is occluded. Occlusions of the VAs proved to be far less critical than occlusions of the ICAs. The worst scenario in terms of reduction in the mean cerebral outflows is a CoW without the first segment of an anterior cerebral artery combined with an occlusion of the contralateral ICA. Furthermore, in patients without any severe occlusion of a carotid or VA, the direction of flow measured at the communicating arteries corresponds to the side of the CoW with an absent or occluded artery. Finally, we study the effect of partial occlusions of the communicating arteries on the cerebral flows, which again confirms that the ACoA is a more important collateral pathway than the PCoAs if an ICA is occluded.     

Twist buckling behavior of arteries

Arteries are often subjected to torsion due to body movement and surgical procedures. While it is essential that arteries remain stable and patent under twisting loads, the stability of arteries under torsion is poorly understood. The goal of this work was to experimentally investigate the buckling behavior of arteries under torsion and to determine the critical buckling torque, the critical buckling twist angle, and the buckling shape. Porcine common carotid arteries were slowly twisted in vitro until buckling occurred while subjected to a constant axial stretch ratio (1.1, 1.3, 1.5 (in vivo level) and 1.7) and lumen pressure (20, 40, 70 and 100 mmHg). Upon buckling, the arteries snapped to form a kink. For a group of six arteries, the axial stretch ratio significantly affected the critical buckling torque (\(p<0.002\)) and the critical buckling twist angle (\(p<0.001\)). Lumen pressure also significantly affected the critical buckling torque (\(p<0.001\)) but had no significant effect on the critical twist angle (\(p=0.067\)). Convex material constants for a Fung strain energy function were determined and fit well with the axial force, lumen pressure, and torque data measured pre-buckling. The material constants are valid for axial stretch ratios, lumen pressures, and rotation angles of 1.3–1.5, 20–100 mmHg, and 0–270\(^\circ \), respectively. The current study elucidates the buckling behavior of arteries under torsion and provides new insight into mechanical instability of blood vessels.     

Flow-induced control of arterial lumen

Blood flow velocity is a factor that affects the diameter of arteries. In order to investigate the flow-induced arterial dilatation, the outer diameter of the femoral, common carotid or renal arteries of anaesthetized cats was measured during perfusion of these arteries with blood or plasma-substituting solutions under conditions of stabilized perfusion pressure. It has been shown that, whatever the perfusate, blood or a substituent, an increase in flow makes the artery to dilate. Consequently, the flow-induced dilatation is not due to any blood-borne humoral factor. As an increase in the solution's viscosity causes dilatation even at constant flow-rate and pressure in the artery, the effect is to be ascribed to the ability of the vascular wall to perceive shear stress. As far as removal of endothelium eliminates the dilatation evoked by increasing flow or fluid viscosity, it may be concluded that the flow-induced dilatation is due to the sensitivity to shear stress of the endothelium.     

Proposition of an outflow boundary approach for carotid artery stenosis CFD simulation

The purpose of this study was to propose an innovative approach of setting outlet boundary conditions for the computational fluid dynamics (CFD) simulation of human common carotid arteries (CCAs) bifurcation based on the concept of energy loss minimisation at flow bifurcation. Comparisons between this new approach and previously reported boundary conditions were also made. The results showed that CFD simulation based on the proposed boundary conditions gave an accurate prediction of the critical stenosis ratio of carotid arteries (at around 65%). Other boundary conditions, such as the constant external pressure (P = 0) and constant outflow ratio, either overestimated or underestimated the critical stenosis ratio of carotid arteries. The patient-specific simulation results furthermore indicated that the calculated internal carotid artery flow ratio at CCA bifurcation (61%) coincided with the result obtained by clinical measurements through the use of Colour Doppler ultrasound.     

Artery buckling affects the mechanical stress in atherosclerotic plaques

Background

Tortuous arteries are often seen in patients with hypertension and atherosclerosis. While the mechanical stress in atherosclerotic plaque under lumen pressure has been studied extensively, the mechanical stability of atherosclerotic arteries and subsequent effect on the plaque stress remain unknown. To this end, we investigated the buckling and post-buckling behavior of model stenotic coronary arteries with symmetric and asymmetric plaque.

Methods

Buckling analysis for a model coronary artery with symmetric and asymmetric plaque was conducted using finite element analysis based on the dimensions and nonlinear anisotropic materials properties reported in the literature.

Results

Artery with asymmetric plaque had lower critical buckling pressure compared to the artery with symmetric plaque and control artery. Buckling increased the peak stress in the plaque and led to the development of a high stress concentration in artery with asymmetric plaque. Stiffer calcified tissue and severe stenosis increased the critical buckling pressure of the artery with asymmetric plaque.

Conclusions

Arteries with atherosclerotic plaques are prone to mechanical buckling which leads to a high stress concentration in the plaques that can possibly make the plaques prone to rupture.

     

Effects of pulsation frequency and endothelial integrity on enhanced arterial transmural filtration produced by pulsatile pressure

The role of the endothelium in regulating transmural fluid filtration into the artery wall under pulsatile pressure and the effects of changes in pulsatile frequency on filtration have received little attention. Previous experiments (Alberding JP, Baldwin AL, Barton JK, and Wiley E. Am J Physiol Heart Circ Physiol 286: H1827-H1835, 2004) demonstrated significantly increased filtration after initial onset of pulsatile pressure compared with that predicted by using parameters measured under steady pressure. To determine the role of the endothelium in this phenomenon, the following experiments were performed on five New Zealand White rabbits (anesthetized with 30 mg/kg pentobarbital sodium). One of each pair of carotid arteries was deendothelialized, and filtration measurements under steady and pulsatile pressure were compared with those made in intact vessels (Alberding JP, Baldwin AL, Barton JK, and Wiley E. Am J Physiol Heart Circ Physiol 286: H1827-H1835, 2004). To determine the effect of increasing pulsatile frequency on arterial filtration, transmural filtration was measured by using pulsatile pressure frequencies of 1 Hz, followed by 2 Hz, in another set of intact arteries (6 arteries and 3 animals). For deendothelialized vessels, the initial increase in filtration after onset of pulsatility was similar to that observed in intact vessels, but the subsequent reduction in filtration was less abrupt. When pulsatile frequency was increased from 1 to 2 Hz in intact arteries, an initial increase in filtration was observed, similar to that obtained after onset of pulsatile pressure subsequent to a steady pressure. The observed responses of arteries to pulsatile pressure, with and without endothelium, or undergoing a frequency change, suggest a dynamic role for the endothelium in regulating transvascular transport in vivo.     

Hemodynamics and the vascular endothelial cytoskeleton

Although there is considerable evidence to suggest that hemodynamics play an important role in vascular disease processes, the exact mechanisms are unknown. With this in mind, we have designed a pulsatile perfusion apparatus which reproducibly delivers pulsatile hemodynamics upon freshly excised canine carotid arteries in vitro. Quantifiable simulations included normotension with normal or lowered flow rates (120/80 mmHg, 120 and 40 ml/min), normotension with lowered or elevated transmural pressures (40-170 mmHg), and elevated pulse pressure (120 and 80 mmHg) with normal (150 ml/min) or elevated rates of flow (300 and 270 ml/min). Arterial biomechanical stresses and cellular behaviors were characterized biochemically and morphologically under all these stimulations which continued for 2-24 h. We found that increased pulse pressure alone had little effect on the total amount of radiolabeled [4-14C]cholesterol present within the medial compartment. However, normotension when coupled with altered transmural pressure yielded a three- to fourfold increase. Combinations of increased pulse pressure and flow potentiated cholesterol uptake by a factor of 10 when compared with normotension control values. Simulations that enhanced carotid arterial cholesterol uptake also influenced the endothelial cytoskeletal array of actin. Stress fibers were not present within the carotid endothelial cells of either the sham controls or the normotension and increased pulse pressure (normal flow) simulations. Endothelial cells lining carotids exposed to elevations in flow or those present within vessels perfused as per simulation b above assembled stress fibers (x = 4 and 10 per cell, respectively) within the time course of these studies. When endothelial cells were subjected to hemodynamic conditions that potentiated maximally cholesterol transport, no diffuse or stress fiber staining could be seen, but the cortical array of actin was intact. These results suggest that those biomechanical stresses that alter endothelial permeability and intimal integrity may do so via cytoskeletal actin signaling.     

A global mathematical model of the cerebral circulation in man

A mathematical model of the cerebral circulation has been formulated. It was based on non-linear equations of pulsatile fluid flow in distensible conduits and applied to a network simulating the entire cerebral vasculature, from the carotid and vertebral arteries to the sinuses and the jugular veins. The quasilinear hyperbolic system of equations was numerically solved using the two-step Lax-Wendroff scheme. The model's results were in good agreement with pressure and flow data recorded in humans during rest. The model was also applied to the study of autoregulation during arterial hypotension. A close relationship between cerebral blood flow (CBF) and capillary pressure was obtained. At arterial pressure of 80 mmHg, the vasodilation of the pial arteries was unable to maintain CBF at its control value. At the lower limit of autoregulation (60 mm Hg), CBF was maintained with a 25% increase of zero transmural pressure diameter of nearly the whole arterial network.     

A two-dimensional numerical analysis of unsteady flow in the carotid artery bifurcation. A comparison with three-dimensional in-vitro measurements and the influence of minor stenoses

In the present study a two-dimensional finite element model for incompressible Newtonian flow is applicated to the modelling of carotid artery flow. In earlier studies, the numerical model was validated experimentally for several flow configurations. In general the pulsatile flow is characterized by reversed flow regions at the non-divider side walls of both the internal and external carotid arteries. The unsteadiness of the flow is associated with rather complex spatial and temporal velocity distributions and leads to temporal variations of the location and length of the reversed flow regions. As a consequence, pronounced spatial and temporal variations in the wall shear stresses are found. At the non-divider side walls, wall shear stresses are relatively low and exhibits an oscillatory behaviour in space and time. At the divider side walls, wall shear stresses are relatively high and approximately follow the flow rate distribution in time. The aim of this study is not only to present two-dimensional calculations but also to compare the calculated two-dimensional velocity profiles with those from three-dimensional experiments. It is observed that in the common carotid artery and in the proximal parts of the internal and external carotid arteries, the two-dimensional numerical model provides valuable information with respect to the three-dimensional configuration. In the more distal parts of especially the internal carotid artery, deviations are found between the two-dimensional numerical and three-dimensional experimental model. These deviations can mainly be attributed to the neglect of the secondary velocity distribution in the two-dimensional model. In the two-dimensional numerical model the influence of a minor stenosis in the internal carotid artery is hardly distinguishable from a minor geometrical variation without stenosis. Full three-dimensional analyses of the influence of minor stenoses are needed to prove numerically whether in-vivo measurements of the axial velocity distribution are useful in the detection of minor stenoses.     

Turbulence detection in a stenosed artery bifurcation by numerical simulation of pulsatile blood flow using the low-Reynolds number turbulence model

The pulsatile blood flow in a partially blocked artery is significantly altered as the flow regime changes through the cardiac cycle. This paper reports on the application of a low-Reynolds turbulence model for computation of physiological pulsatile flow in a healthy and stenosed carotid artery bifurcation. The human carotid artery was chosen since it has received much attention because atherosclerotic lesions are frequently observed. The Wilcox low-Re k-omega turbulence model was used for the simulation since it has proven to be more accurate in describing transition from laminar to turbulent flow. Using the FIDAP finite element code a validation showed very good agreement between experimental and numerical results for a steady laminar to turbulent flow transition as reported in a previous publication by the same authors. Since no experimental or numerical results were available in the literature for a pulsatile and turbulent flow regime, a comparison between laminar and low-Re turbulent calculations was made to further validate the turbulence model. The results of this study showed a very good agreement for velocity profiles and wall shear stress values for this imposed pulsatile laminar flow regime. To explore further the medical aspect, the calculations showed that even in a healthy or non-stenosed artery, small instabilities could be found at least for a portion of the pulse cycle and in different sections. The 40% and 55% diameter reduction stenoses did not significantly change the turbulence characteristics. Further results showed that the presence of 75% stenoses changed the flow properties from laminar to turbulent flow for a good portion of the cardiac pulse. A full 3D simulation with this low-Re-turbulence model, coupled with Doppler ultrasound, can play a significant role in assessing the degree of stenosis for cardiac patients with mild conditions.