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

Comparative analysis of RADAR vs. conventional techniques for AVF maturation in patients with blood viscosity and vessel elasticity-related diseases through fluid-structure interaction modeling: Anemia, hypertension, and diabetes

PURPOSE

This study aims to compare two surgical techniques, the standard Vein-to-Artery and the newer Artery-to-Vein (Radial Artery Deviation And Reimplantation; RADAR), for enhancing the success of Arterio-Venous Fistula maturation in end-stage renal disease patients. The impact of diseases like anemia, diabetes, hypertension, and chronic kidney disease were considered. The goals are to advance Arterio-Venous Fistula (AVF) surgery, improve patient outcomes, and contribute to evidence-based surgical guidelines.

METHODS

Fluid-structure interaction modeling was employed to investigate how hemodynamic and mechanical stresses impact arteriovenous fistula maturation, with a particular focus on the role of wall shear stress in determining maturation outcomes. The critical threshold for vessel injury was identified as wall shear stress values exceeding 35 N/m2, while stenosis formation was projected to occur at levels below 1 N/m2. This work introduced a novel approach by considering disease-related factors, including blood viscosity (anemia), and vessel elasticity (diabetes, hypertension, and chronic kidney diseases), which directly influence hemodynamics and the generation of wall shear stress. Furthermore, the model was designed to incorporate varying thicknesses and elasticities for both the vein and artery, accurately representing authentic vascular anatomy.

RESULTS

The RADAR technique has demonstrated superior performance compared to the standard technique by providing appropriate wall shear stress in critical regions and minimizing the risk of wall damage. Its use of a thicker vessel also reduces the risk of vessel injury, making it particularly effective for patients with Chronic Kidney Disease (CKD), hypertension, anemia, and diabetes, ensuring optimal blood flow and fewer complications. However, there are minor concerns about stenosis formation in hypertension and anemia cases, which could be mitigated by adjusting the anastomosis angle to be lower than 30°.

CONCLUSION

Diabetes and hypertension have significant physiological effects that increase the risks associated with arteriovenous fistula maturation. The anemic condition resulting from CKD may help reduce vessel injury but raises concerns about potential stenosis formation. Despite these co-morbidities, the RADAR technique has demonstrated its ability to induce more favorable hemodynamic changes, promoting arteriovenous fistula maturation.

Pathophysiology, cellular and molecular mechanisms of large and small vessel diseases

Cerebrovascular disease (CVD) is the second most common cause of cognitive impairment and dementia in aged population. CVD presents in a myriad number of clinical ways based on the functional location of pathology. While primary clinical emphasis has been placed on motor, speech and visual deficits, vascular cognitive decline is a vastly under recognized and devastating condition afflicting millions of Americans. CVD, a disease of the blood vessels that supply blood to brain involves an integration between small and large vessels. Cerebral large vessel diseases (LVD) are associated with atherosclerosis, artery-to-artery embolism, intracardiac embolism and a large vessel stroke leading to substantial functional disability. Cerebral small vessel disease (SVD) is critically involved in stroke, brain hemorrhages, cognitive decline and functional loss in elderly patients. An evolving understanding of cellular and molecular mechanisms emphasizes that inflammatory vascular changes contribute to systemic pathologic conditions of the central nervous systems (CNS), with specific clinical presentations including, cognitive decline. Advances in an understanding of pathophysiology of disease processes and therapeutic interventions may help improve outcomes. This review will focus on large and small vessels diseases and their relationship to vascular cognitive decline, atherosclerosis, stroke, and inflammatory neurodegeneration. We will also emphasize the molecular and cellular mechanisms, as well as genetic and epigenetic factors associated with LVD and SVD.

Inducing Vascular Grammars for Anomaly Classification in Brain Angiograms

As machine learning is used to make strides in medical diagnostics, few methods provide heuristics from which human doctors can learn directly. This work introduces a method for leveraging human observable structures, such as macroscale vascular formations, for producing assessments of medical conditions with relatively few training cases, and uncovering patterns that are potential diagnostic aids. The approach draws on shape grammars, a rule-based technique, pioneered in design and architecture, and accelerated through a recursive subgraph mining algorithm. The distribution of rule instances in the data from which they are induced is then used as an intermediary representation enabling common classification and anomaly detection approaches to identify indicative rules with relatively small data sets. The method is applied to seven-tesla time-of-flight angiography MRI (n = 54) of human brain vasculature. The data were segmented and induced to generate representative grammar rules. Ensembles of rules were isolated to implicate vascular conditions reliably. This application demonstrates the power of automated structured intermediary representations for assessing nuanced biological form relationships, and the strength of shape grammars, in particular for identifying indicative patterns in complex vascular networks.

In Vitro Validation of Regional Circumferential Strain Assessment in a Phantom Aortic Model Using Cine Displacement Encoding With Stimulated Echoes MRI

BACKGROUND

A novel application of cine Displacement ENcoding with Stimulated Echoes Magnetic Resonance Imaging (DENSE MRI) has recently been described to assess regional heterogeneities in circumferential strain around the aortic wall in vivo; however, validation is first required for successful clinical translation.

PURPOSE

To validate the quantification of regional circumferential strain around the wall of an aortic phantom using DENSE MRI.

STUDY TYPE

In vitro phantom study.

POPULATION

Three polyvinyl alcohol aortic phantoms with eight axially oriented nitinol wires embedded evenly around the walls.

FIELD STRENGTH/SEQUENCE

3 T; gradient-echo aortic DENSE MRI with spiral cine readout, gradient-echo phase-contrast MRI (PCMR) with Cartesian cine readout.

ASSESSMENT

Phantoms were connected to a pulsatile flow loop and peak DENSE-derived regional circumferential Green strains at 16 equally spaced sectors around the wall were assessed according to previously published algorithms. "True" regional circumferential strains were calculated by manually tracking displacements of the nitinol wires by two independent observers. Normalized circumferential strains (NCS) were calculated by dividing regional strains by the mean strain. Finally, DENSE-derived regional strain was corrected by multiplying regional DENSE NCS by the mean strain calculated from the diameter change on the PCMR.

STATISTICAL TESTS

One-sample t-test, Paired-sample t-test, and analysis of variance with Bonferroni correction, coefficient of variation (CoV), Bland-Altman analysis; P < 0.05 was considered statistically significant.

RESULTS

Aortic DENSE MRI significantly overestimated circumferential strain compared to the wire-tracking method (mean difference and SD 0.030 ± 0.014, CoV 0.31). However, NCS demonstrated good agreement between DENSE and wire-tracking data (mean difference 0.000 ± 0.172, CoV 0.15). After correcting the DENSE-derived regional strain, the mean difference in regional circumferential strain between DENSE and wire-tracking was significantly reduced to 0.006 ± 0.008, and the CoV was reduced to 0.18.

DATA CONCLUSION

For aortic phantoms with mild spatial heterogeneity in circumferential strain, the previously published aortic DENSE MRI technique successfully assessed the regional NCS distribution but overestimated the mean strain. This overestimation is correctable by computing a more accurate mean circumferential strain using a separate cine scan.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Numerical simulation of pulsatile blood flow characteristics in a multi stenosed coronary artery

BACKGROUND

Coronary artery disease is reported as one of the most common sources of death all over the world. The presence of stenosis (plaque) in the coronary arteries results in the restriction of blood supply, which leads to myocardial infarction.

OBJECTIVE

The aim of this study was to investigate the effect of multi stenosis on hemodynamics parameters in idealized coronary artery models with varying degrees of stenosis and interspace distance between the stenosis.

METHODS

A finite volume-based software package (Ansys CFX version 17.2) was employed to model the blood flow. The hemodynamic stenosis parameters of blood, such as the pressure, velocity, and wall shear stress were obtained.

RESULTS

The computed results showed that the pressure drop is maximum across the 90% area stenosis (AS). The pressure drop is increased as the distance between the proximal and distal stenosis is decreased across the proximal stenosis for the model P70_D70 during the systolic period of the cardiac cycle. A recirculation zone is formed behind the stenosis and is restricted by the occurrence of distal stenosis as the interspacing distance decreases, which could lead to further progression of stenosis in the flow-disturbed area. The wall shear stress was found to increase as the distance between the proximal and distal stenosis is increased across the distal stenosis. The maximum wall shear stress was found at 90% AS.

CONCLUSIONS

In the clinical diagnosis, an overestimation of distal stenosis severity could be possible. Furthermore, the low wall shear stress zone in between the proximal and distal stenosis may help atherosclerotic growth or merge adjacent stenosis.

Fluid-Structure Interaction Analyses of Biological Systems Using Smoothed-Particle Hydrodynamics

Due to the inherent complexity of biological applications that more often than not include fluids and structures interacting together, the development of computational fluid-structure interaction models is necessary to achieve a quantitative understanding of their structure and function in both health and disease. The functions of biological structures usually include their interactions with the surrounding fluids. Hence, we contend that the use of fluid-structure interaction models in computational studies of biological systems is practical, if not necessary. The ultimate goal is to develop computational models to predict human biological processes. These models are meant to guide us through the multitude of possible diseases affecting our organs and lead to more effective methods for disease diagnosis, risk stratification, and therapy. This review paper summarizes computational models that use smoothed-particle hydrodynamics to simulate the fluid-structure interactions in complex biological systems.

Using intravascular ultrasound image-based fluid-structure interaction models and machine learning methods to predict human coronary plaque vulnerability change

Plaque vulnerability prediction is of great importance in cardiovascular research. In vivo follow-up intravascular ultrasound (IVUS) coronary plaque data were acquired from nine patients to construct fluid-structure interaction models to obtain plaque biomechanical conditions. Morphological plaque vulnerability index (MPVI) was defined to measure plaque vulnerability. The generalized linear mixed regression model (GLMM), support vector machine (SVM) and random forest (RF) were introduced to predict MPVI change (ΔMPVI = MPVI_follow-up‒MPVI_baseline) using ten risk factors at baseline. The combination of mean wall thickness, lumen area, plaque area, critical plaque wall stress, and MPVI was the best predictor using RF with the highest prediction accuracy 91.47%, compared to 90.78% from SVM, and 85.56% from GLMM. Machine learning method (RF) improved the prediction accuracy by 5.91% over that from GLMM. MPVI was the best single risk factor using both GLMM (82.09%) and RF (78.53%) while plaque area was the best using SVM (81.29%).

Effect of Axial Stretch on Lumen Collapse of Arteries

The stability of the arteries under in vivo pressure and axial tension loads is essential to normal arterial function, and lumen collapse due to buckling can hinder the blood flow. The objective of this study was to develop the lumen buckling equation for nonlinear anisotropic thick-walled arteries to determine the effect of axial tension. The theoretical equation was developed using exponential Fung strain function, and the effects of axial tension and residual stress on the critical buckling pressure were illustrated for porcine coronary arteries. The buckling behavior was also simulated using finite-element analysis. Our results demonstrated that lumen collapse of arteries could occur when the transmural pressure is negative and exceeded a critical value. This value depends upon the axial stretch ratio and material properties of the arterial wall. Axial tensions show a biphasic effect on the critical buckling pressure. The lumen aspect ratio of arteries increases nonlinearly with increasing external pressure beyond the critical value as the lumen collapses. These results enhance our understanding of artery lumen collapse behavior.

Nonlinear Dynamics of Dacron Aortic Prostheses Conveying Pulsatile Flow

This study addresses the dynamic response to pulsatile physiological blood flow and pressure of a woven Dacron graft currently used in thoracic aortic surgery. The model of the prosthesis assumes a cylindrical orthotropic shell described by means of nonlinear Novozhilov shell theory. The blood flow is modeled as Newtonian pulsatile flow, and unsteady viscous effects are included. Coupled fluid–structure Lagrange equations for open systems with wave propagation subject to pulsatile flow are applied. Physiological waveforms of blood pressure and velocity are approximated with the first eight harmonics of the corresponding Fourier series. Time responses of the prosthetic wall radial displacement are considered for two physiological conditions: at rest (60 bpm) and at high heart rate (180 bpm). While the response at 60 bpm reproduces the behavior of the pulsatile pressure, higher harmonics frequency contributions are observed at 180 bpm altering the shape of the time response. Frequency-responses show resonance peaks for heart rates between 130 bpm and 200 bpm due to higher harmonics of the pulsatile flow excitation. These resonant peaks correspond to unwanted high-frequency radial oscillations of the vessel wall that can compromise the long-term functioning of the prosthesis in case of significant physical activity. Thanks to this study, the dynamic response of Dacron prostheses to pulsatile flow can be understood as well as some possible complications in case of significant physical activity.

The relationship between coronary artery distensibility and fractional flow reserve

Discordance between angiography-based anatomical assessment of coronary stenosis severity and fractional flow reserve (FFR) has been attributed to several factors including lesion length and irregularity, and the myocardial territory supplied by the target vessel. We sought to examine if coronary arterial distensibility is an independent contributor to this discordance. There were two parts to this study. The first consisted of “in silico” models of 26 human coronary arteries. Computational fluid dynamics-derived FFR was calculated for fully rigid, partially distensible and fully distensible models of the 26 arteries. The second part of the study consisted of 104 patients who underwent coronary angiography and FFR measurement. Distensibility at the lesion site (Distensibility_MLA) and for the reference vessel (Distensibility_Ref) was determined by analysing three-dimensional angiography images during end-systole and end-diastole. Computational fluid dynamics-derived FFR was 0.67±0.19, 0.70±0.18 and 0.75±0.17 (P<0.001) in the fully rigid, partially distensible and fully distensible models respectively. FFR correlated with both Distensibility_MLA (r = 0.36, P<0.001) and Distensibility_Ref (r = 0.44, P<0.001). Two-way ANCOVA analysis revealed that Distensibility_MLA (F (1, 100) = 4.17, p = 0.031) and percentage diameter stenosis (F (1, 100) = 60.30, p < 0.01) were both independent predictors of FFR. Coronary arterial distensibility is a novel, independent determinant of FFR, and an important factor contributing to the discordance between anatomical and functional assessment of stenosis severity.

MaterialCloning: Acquiring Elasticity Parameters from Images for Medical Applications

We present a practical approach for automatically estimating the material properties of soft bodies from two sets of images, taken before and after deformation. We reconstruct 3D geometry from the given sets of multiple-view images; we use a coupled simulation-optimization-identification framework to deform one soft body at its original, non-deformed state to match the deformed geometry of the same object in its deformed state. For shape correspondence, we use a distance-based error metric to compare the estimated deformation fields against the actual deformation field from the reconstructed geometry. The optimal set of material parameters is thereby determined by minimizing the error metric function. This method can simultaneously recover the elasticity parameters of multiple types of soft bodies using Finite Element Method-based simulation (of either linear or nonlinear materials undergoing large deformation) and particle-swarm optimization methods. We demonstrate this approach on real-time interaction with virtual organs in patient-specific surgical simulation, using parameters acquired from low-resolution medical images. We also highlight the results on physics-based animation of virtual objects using sketches from an artist's conception.

Reconfiguration of the Carotid Artery after Angioplasty and Stenting: A Case Report and Review of the Literature

Severe carotid stenosis or occlusion may cause insufficient blood flow and lead to distal artery wall collapse and extensive lumen contraction. Whether this 'adaptive narrowing' can restitute after carotid recanalization is unclear. We report a patient with global ischemia due to occlusions of bilateral carotid and right vertebral arteries. The occluded left carotid was recanalized successfully with angioplasty and stenting. The adaptively narrowed distal carotid did not restitute immediately but regained its morphology 1 week after the procedure. Carotid adaptive narrow distal occlusion or stenosis may not regain its original morphology immediately but several days after recanalization. This knowledge is instructive for treating occlusive carotid diseases.

Assessing abdominal aorta narrowing using computational fluid dynamics

This paper investigates the effect of developing arterial blockage at the abdominal aorta on the blood pressure waves at an externally accessible location suitable for invasive measurements such as the brachial and the femoral arteries. Arterial blockages are created surgically within the abdominal aorta of healthy Wistar rats to create narrowing resemblance conditions. Blood pressure is measured using a catheter inserted into the right femoral artery. Measurements are taken at the baseline healthy condition as well as at four different severities (20, 50, 80 and 100 %) of arterial blockage. In vivo and in vitro measurements of the lumen diameter and wall thickness are taken using magnetic resonance imaging and microscopic techniques, respectively. These data are used to validate a 3D computational fluid dynamics model which is developed to generalize the outcomes of this work and to determine the arterial stress and strain under the blockage conditions. This work indicates that an arterial blockage in excess of 20 % of the lumen diameter significantly influences the pressure wave and reduces the systolic blood pressure at the right femoral artery. High wall shear stresses and low circumferential strains are also generated at the blockage site.

Stability of carotid artery under steady-state and pulsatile blood flow: a fluid-structure interaction study

It has been shown that arteries may buckle into tortuous shapes under lumen pressure, which in turn could alter blood flow. However, the mechanisms of artery instability under pulsatile flow have not been fully understood. The objective of this study was to simulate the buckling and post-buckling behaviors of the carotid artery under pulsatile flow using a fully coupled fluid-structure interaction (FSI) method. The artery wall was modeled as a nonlinear material with a two-fiber strain-energy function. FSI simulations were performed under steady-state flow and pulsatile flow conditions with a prescribed flow velocity profile at the inlet and different pressures at the outlet to determine the critical buckling pressure. Simulations were performed for normal (160 ml/min) and high (350 ml/min) flow rates and normal (1.5) and reduced (1.3) axial stretch ratios to determine the effects of flow rate and axial tension on stability. The results showed that an artery buckled when the lumen pressure exceeded a critical value. The critical mean buckling pressure at pulsatile flow was 17-23% smaller than at steady-state flow. For both steady-state and pulsatile flow, the high flow rate had very little effect (&lt;5%) on the critical buckling pressure. The fluid and wall stresses were drastically altered at the location with maximum deflection. The maximum lumen shear stress occurred at the inner side of the bend and maximum tensile wall stresses occurred at the outer side. These findings improve our understanding of artery instability in vivo.

A numerical study of the effect of varied blood pressure on the stability of carotid atherosclerotic plaque

Background

Blood pressure (BP) is associated with early atherosclerosis and plaque rupture because the BP variability can significantly affect the blood flow velocity and shear stress over the plaque. However, the mechanical response of BP variability to the plaque remains unclear. Therefore, we investigated the correlation between different maximum systolic blood pressure (SBP) and the stress distribution on plaque, as well as the stress over the plaque and blood velocity around the plaque using different BP variations, which are the BP variability in different phases during one cardiac cycle and beat-to-beat BP variability.

Method

We established a two-dimensional artery model with stenosis at the degree of 62.5%. Eight combinations of pulsatile pressure gradients between the inflow and outflow were implemented at the model. Three levels of fibrous cap thickness were taken into consideration to investigate the additional effect on the BP variability. Wall shear stress and stress/strain distribution over the plaque were derived as well as the oscillation shear index (OSI) to analyze the impact of the changing rate of BP.

Result

The stresses at diastole were 2.5% ± 1.8% lower than that at systole under the same pressure drop during one cycle. It was also found that elevated SBP might cause the immediate increment of stress in the present cycle (292% ± 72.3%), but slight reduction in the successive cycle (0.48% ± 0.4%).

Conclusion

The stress/strain distribution over the plaque is sensitive to the BP variability during one cardiac cycle, and the beat-to-beat BP variability could cause considerable impact on the progression of atherosclerosis in long-term.

Effect of postural changes on normal and stenosed common carotid artery using FSI

Gravity associated with postural changes has a strong bearing on haemodynamics of blood flow in arteries. Its effect on stenosed cases has not been widely investigated. In the present study, variation observed in blood flow during postural changes is investigated for different conditions like standing, sleeping and head-down position. A fluid structure interaction study is carried out for idealized normal and 75% eccentric and concentric stenosed common carotid normal artery. The results clearly indicate the effects of altered gravity on flow conditions. It was found to be very significant during head-down position and demonstrated very high arterial blood pressure in stenosed common carotid when compared with normal carotid.

Artery buckling: new phenotypes, models, and applications

Arteries are under significant mechanical loads from blood pressure, flow, tissue tethering, and body movement. It is critical that arteries remain patent and stable under these loads. This review summarizes the common forms of buckling that occur in blood vessels including cross-sectional collapse, longitudinal twist buckling, and bent buckling. The phenomena, model analyses, experimental measurements, effects on blood flow, and clinical relevance are discussed. It is concluded that mechanical buckling is an important issue for vasculature, in addition to wall stiffness and strength, and requires further studies to address the challenges. Studies of vessel buckling not only enrich vascular biomechanics but also have important clinical applications.

A simplified murine intimal hyperplasia model founded on a focal carotid stenosis

Murine models offer a powerful tool for unraveling the mechanisms of intimal hyperplasia and vascular remodeling, although their technical complexity increases experimental variability and limits widespread application. We describe a simple and clinically relevant mouse model of arterial intimal hyperplasia and remodeling. Focal left carotid artery (LCA) stenosis was created by placing 9-0 nylon suture around the artery using an external 35-gauge mandrel needle (middle or distal location), which was then removed. The effect of adjunctive diet-induced obesity was defined. Flowmetry, wall strain analyses, biomicroscopy, and histology were completed. LCA blood flow sharply decreased by ∼85%, followed by a responsive right carotid artery increase of ∼71%. Circumferential strain decreased by ∼2.1% proximal to the stenosis in both dietary groups. At 28 days, morphologic adaptations included proximal LCA intimal hyperplasia, which was exacerbated by diet-induced obesity. The proximal and distal LCA underwent outward and negative inward remodeling, respectively, in the mid-focal stenosis (remodeling indexes, 1.10 and 0.53). A simple, defined common carotid focal stenosis yields reproducible murine intimal hyperplasia and substantial differentials in arterial wall adaptations. This model offers a tool for investigating mechanisms of hemodynamically driven intimal hyperplasia and arterial wall remodeling.

CABG MODELS FLOW SIMULATION STUDY ON THE EFFECTS OF VALVE REMNANTS IN THE VENOUS GRAFT

Venous valves and sinuses are frequently observed in vein grafts in the coronary artery bypass grafts (CABG). However, from the biomedical engineering viewpoint, vein grafts are always assumed as smooth tubes in the existing simulations, and no effort has been made to investigate the effects of jaggedness of the graft inner wall due to the valve cusps remnants and valve sinus (in case of valve-stripped saphenous vein (SV) grafts) on the blood flow patterns and hemodynamic parameters (HPs). In this paper, the effects of the inner surface irregularities of a vein graft on the blood flow is investigated in the graft as well as in the distal anastomotic region, with a more realistic geometry of valve-stripped SV, by means of numerical simulation of pulsatile, Newtonian blood flow. The simulation results demonstrate that the valve remnants and sinuses cause disturbances in the flow field within the graft (due to vortices formation within the valve sinuses) and undesirable distribution of HPs, which can result in early atherosclerotic lesion development in the graft.

Effect of branchings on blood flow in the system of human coronary arteries

In this work, we investigate the behavior of the pulsatile blood flow in the system of human coronary arteries. Blood is modeled as an incompressible non-Newtonian fluid. The transient phenomena of blood flow through the coronary system are simulated by solving the three dimensional unsteady state Navier-Stokes equations and continuity equation. Distributions of velocity, pressure and wall shear stresses are determined in the system under pulsatile conditions on the boundaries. Effect of branching vessel on the flow problem is investigated. The numerical results show that blood pressure in the system with branching vessels of coronary arteries is lower than the one in the system with no branch. The magnitude of wall shear stresses rises at the bifurcation.

The effects of stenosis severity on the hemodynamic parameters-assessment of the correlation between stress phase angle and wall shear stress

To study the effects of increase in the degree of stenosis severity and subsequent complexity of hemodynamic patterns on hemodynamic parameters, experimental investigations and numerical simulations were performed. The correlations between the large negative Stress Phase Angle (SPA), the low mean Wall Shear Stress (WSS) and high Oscillatory Shear Index (OSI) were investigated at the distal shoulder and post-stenotic regions as the outcomes of elevated stenosis severity. Models included non-Newtonian fluid flow in stenotic arteries with 30-80% symmetrical stenoses. To study the interactions between pulsatile WSS and pulsatile wall circumferential stress (WCS) acting on endothelial cells, SPA as the phase difference between WSS and WCS waves was used. Moreover, the distribution of SPA on the lumen axis was compared to the distributions of the mean WSS and OSI that have been regarded until now as the determinants of atherosclerosis-prone regions. Results indicate that an increase in stenosis severity, not only affects the mean WSS, mean WCS and pulse amplitudes, but also influences the phase difference between them. The SPA is large negative on the distal shoulder and post-stenotic areas where atherosclerotic plaque develops. The increasing stenosis severity and the subsequent increasing complexity of hemodynamic patterns affect the correlation between any of the low mean WSS and high OSI with large negative SPA, such that it not only leads to create and develop some regions where the correlation between any of the low mean WSS and high OSI with large negative SPA is well but also leads to create and develop other regions where such correlations fail.

Coupled fluid-wall modelling of steady flow in stenotic carotid arteries

Arterial stenoses may cause critical blood flow and wall conditions leading to clinical complications. In this paper computational models of stenotic carotid arteries are proposed and the vessel wall collapse phenomenon is studied. The models are based on fluid-structure interactions (FSI) between blood and the arterial walls. Coupled finite element and computational fluid dynamics methods are used to simultaneously solve for stress and displacement in the solid, and for pressure, velocity and shear stress in the fluid domain. Results show high wall shear stress at the stenosis throat and low (negative) values accompanied by disturbed flow patterns downstream of the stenosis. The wall circumferential stress varies abruptly from tensile to compressive along the stenosis with high stress concentration on the plaque shoulders showing regions of possible plaque rupture. Wall compression and collapse are observed for severe cases. Post-stenotic collapse of the arterial wall occurs for stenotic severity as low as 50%, with the assumption that a given amount of blood flow needs to pass the stenotic artery; whereas if constant pressure drop should be maintained across a constriction, then collapse happens at severity of 75% and above. The former assumption is based on the requirement of adequate blood supply to the downstream organs/tissue, while the latter stems from the fact that the pumping mechanism of the body has a limited capacity in regulating blood pressure, in case a stenosis appears in the vasculature.

Introduction to the biomechanics of carotid plaque pathogenesis and rupture: review of the clinical evidence

The management of patients with asymptomatic carotid disease is currently under debate and new methods are warranted for better risk stratification. The role of the biomechanical properties of the atherosclerotic arterial wall together with the effect of different stress types in plaque destabilisation has only been recently investigated. PubMed and Scopus databases were reviewed. There is preliminary clinical evidence demonstrating that the analysis of the combined effect of the various types of biomechanical stress acting on the carotid plaque may help us to identify the vulnerable plaque. At present, MRI and two-dimensional ultrasound are combined with fluid-structure interaction techniques to produce maps of the stress variation within the carotid wall, with increased cost and complexity. Stress wall analysis can be a useful tool for carotid plaque evaluation; however, further research and a multidisciplinary approach are deemed as necessary for further development in this direction.

A review of carotid atherosclerosis and vascular cognitive decline: a new understanding of the keys to symptomology

This review encourages the reader to consider cerebral vascular disease beyond the traditional clinical end points of major motor and speech strokes and to consider the possible impact of embolic cerebral vascular disease on vascular cognitive decline. This article examines the issue of "silent" strokes in the relationship between the structural stability of atherosclerotic carotid plaque and the development of nonmotor symptomatology, including cognitive decline. It addresses the question of the role of carotid emboli in silent stroke and their cognitive sequelae. In a study of endarterectomy patients, we relate plaque elasticity and its development of mechanical strain features and thinning of stabilizing fibrous cap at the point of these mechanical strain features. The possibility that microemboli from such mechanically unstable carotid plaques could contribute to silent strokes led to a study of cognitive function in such patients. A linear relationship between the process of mechanically unstable areas of carotid plaques and cognitive decline suggests a contributory role for such a process in silent strokes.

Characterization of the highly nonlinear and anisotropic vascular tissues from experimental inflation data: a validation study toward the use of clinical data for in-vivo modeling and analysis

We study whether an inverse modeling approach is applicable for characterizing vascular tissue subjected to various levels of internal pressure and axial stretch that approximate in-vivo conditions. To compensate for the limitation of axial-displacement/pressure/diameter data typical of clinical data, which does not provide information about axial force, we propose to constrain the ratio of axial to circumferential elastic moduli to a typical range. Vessel wall constitutive behavior is modeled with a transversely isotropic hyperelastic equation that accounts for dispersed collagen fibers. A single-layer and a bi-layer approximation to vessel ultrastructure are examined, as is the possibility of obtaining the fiber orientation as part of the optimization. Characterization is validated against independent pipette-aspiration biaxial data on the same samples. It was found that the single-layer model based on homogeneous wall assumption could not reproduce the validation data. In contrast, the constrained bi-layer model was in excellent agreement with both types of experimental data. Due to covariance, estimations of fiber angle were slightly outside of the normal range, which can be resolved by predefining the angles to normal values. Our approach is relatively invariant to a constant or a variable axial response. We believe that it is suitable for in-vivo characterization.

Direct numerical simulation of transitional flow in a stenosed carotid bifurcation

The blood flow dynamics of a stenosed, subject-specific, carotid bifurcation were numerically simulated using the spectral element method. Pulsatile inlet conditions were based on in vivo color Doppler ultrasound measurements of blood velocity. The results demonstrated the transitional or weakly turbulent state of the blood flow, which featured rapid velocity and pressure fluctuations in the post-stenotic region of the internal carotid artery (ICA) during systole and laminar flow during diastole. High-frequency vortex shedding was greatest downstream of the stenosis during the deceleration phase of systole. Velocity fluctuations had a frequency within the audible range of 100-300Hz. Instantaneous wall shear stress (WSS) within the stenosis was relatively high during systole ( approximately 25-45Pa) compared to that in a healthy carotid. In addition, high spatial gradients of WSS were present due to flow separation on the inner wall. Oscillatory flow reversal and low pressure were observed distal to the stenosis in the ICA. This study predicts the complex flow field, the turbulence levels and the distribution of the biomechanical stresses present in vivo within a stenosed carotid artery.

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.

Influence of microcalcifications on vulnerable plaque mechanics using FSI modeling

Sudden heart attacks remain one of the primary causes of premature death in the developed world. Asymptomatic vulnerable plaques that rupture are believed to prompt such fatal heart attacks and strokes. The role of microcalcifications in the vulnerable plaque rupture mechanics is still debated. Recent studies suggest the microcalcifications increase the plaque vulnerability. In this manuscript we present a numerical study of the role of microcalcifications in plaque vulnerability in an eccentric stenosis model using a transient fluid-structure interaction (FSI) analysis. Two cases are being compared (i) in the absence of a microcalcification (ii) with a microcalcification spot fully embedded in the fibrous cap. Critical plaque stress/strain conditions were affected considerably by the presence of a calcified spot, and were dependent on the timing (phase) during the flow cycle. The vulnerable plaque with the embedded calcification spot presented higher wall stress concentration region in the fibrous cap a bit upstream to the calcified spot, with stress propagating to the deformable parts of the structure around the calcified spot. Following previous studies, this finding supports the hypothesis that microcalcifications increase the plaque vulnerability. Further studies in which the effect of additional microcalcifications and parametric studies of critical plaque cap thickness based on plaque properties and thickness, will help to establish the mechanism by which microcalcifications weaken the plaque and may lead to its rupture.

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.

Wall shear stress gradient analysis within an idealized stenosis using non-Newtonian flow

OBJECTIVE

The endothelium is functionally regulated by the magnitude and spatiotemporal gradients of wall shear stress (WSS). Although flow separation and reversal occur beyond high-grade stenoses, little is known of the WSS pattern within clinically relevant mild to moderate stenoses.

METHODS

An axisymmetric geometry with 25, 50, and 75% stenosis criteria (quantified in accordance with the North American Symptomatic Carotid Endarterectomy Trial) was used to generate a high-resolution, hybrid, tetrahedral-hexahedral computational mesh with boundary-layer enrichment to improve near-wall shear stress gradient (WSSG) computation. Time-dependent computational fluid dynamic analysis was performed using a non-Newtonian Carreau-Yasuda model of blood to yield the shear-dependent viscosity.

RESULTS

Transition to secondary flow patterns was demonstrated in stenoses of 25, 50, and 75%. A focal region with near-wall flow reversal and retrograde WSS was identified within the stenosis itself and was found to migrate cyclically during the cardiac pulse. A zone of zero WSS and divergent WSSG that shifts in toward the throat with increasing stenotic severity was identified. Focal zones of high WSSG with converging and/or diverging direction were uncovered within the stenosis itself, as were expected changes in the distal poststenotic region. These zones of divergent WSSG shift over a substantial length of the stenosis during the course of the cardiac cycle.

CONCLUSION

Luminal WSS demonstrates dynamic direction reversal and high spatial gradients within the distal stenosis throat of even clinically moderate lesions. These findings shed light on the complex vessel wall hemodynamics within clinical stenoses and reveal a mechanical microenvironment that is conducive to perpetual endothelial functional dysregulation and stenosis progression.

PREDICTION OF COMPLEX FLOW PATTERNS IN INTRACRANIAL ATHEROSCLEROTIC DISEASE USING COMPUTATIONAL FLUID DYNAMICS

OBJECTIVE

Although carotid and vertebral intracranial atherosclerotic disease (ICAD) can lead to both hemodynamic insufficiency and thromboembolism, its fluid dynamic properties remain undefined because of its intricate features and complex three-dimensional geometry. We used computational fluid dynamic (CFD) analysis to model the hemodynamics of symptomatic ICAD lesions.

METHODS

Nine ICAD lesions (six carotid, two vertebral, one middle cerebral) underwent high-resolution catheter-based digital rotational angiography. The reconstructed three-dimensional volumes of the target lesions were segmented and used to generate hybrid computational meshes. Dynamic pulsatile CFD analysis was performed using a non-Newtonian shear-dependent model of blood's viscosity.

RESULTS

CFD results revealed complex flow patterns within ICAD lesions with midstenotic shear rates of greater than 19,000/s, sufficiently high to induce high-shear platelet activation. Vorticity and helicity within the stenoses were followed by sudden deceleration with formation of vortex cores. Pressure gradients were significant mostly at greater than 75% stenosis with a mean time-averaged drop of 27.2 ±17.8 mmHg. Unlike the smoothly-varying helicity imparted by the three-dimensional anatomy of the intracranial circulation, poststenotic regions of ICAD lesions showed significant and rapidly fluctuating helicity and vorticity patterns, which may contribute to the propagation of platelets activated by the high shear region within the stenosis throat. Stent angioplasty restored the hemodynamic profile of ICAD lesions to within contralateral controls.

CONCLUSION

Patient-based symptomatic ICAD lesions studied using CFD analysis appear to harbor a hemodynamically pathological environment that favors the activation, aggregation and distal embolization of platelets and is reversed by endovascular stent angioplasty.

A biomechanical model of artery buckling

The stability of arteries under blood pressure load is essential to the maintenance of normal arterial function and the loss of stability can lead to tortuosity and kinking that are associated with significant clinical complications. However, mechanical analysis of arterial bent buckling is lacking. To address this issue, this paper presents a biomechanical model of arterial buckling. Using an elastic cylindrical arterial model, the mechanical equations for arterial buckling were developed and the critical buckling pressure was found to be a function of the wall stiffness (Young's modulus), arterial radius, length, wall thickness, and the axial strain. Both the model equations and experimental results demonstrated that the critical pressure is related to the axial strain. Arteries may buckle and become tortuous due to reduced (subphysiological) axial strain, hypertensive pressure, and a weakened wall. These results are in accordance with, and provide a possible explanation to the clinical observations that hypertension and aging are the risk factors for arterial tortuosity and kinking. The current model is also applicable to veins and ureters.

Numerical analysis of pulsatile blood flow and vessel wall mechanics in different degrees of stenoses

Hemodynamics factors and biomechanical forces play key roles in atherogenesis, plaque development and final rupture. In this paper, we investigated the flow field and stress field for different degrees of stenoses under physiological conditions. Disease is modelled as axisymmetric cosine shape stenoses with varying diameter reductions of 30%, 50% and 70%, respectively. A simulation model which incorporates fluid-structure interaction, a turbulence model and realistic boundary conditions has been developed. The results show that wall motion is constrained at the throat by 60% for the 30% stenosis and 85% for the 50% stenosis; while for the 70% stenosis, wall motion at the throat is negligible through the whole cycle. Peak velocity at the throat varies from 1.47 m/s in the 30% stenosis to 3.2m/s in the 70% stenosis against a value of 0.78 m/s in healthy arteries. Peak wall shear stress values greater than 100 Pa were found for > or =50% stenoses, which in vivo could lead to endothelial stripping. Maximum circumferential stress was found at the shoulders of plaques. The results from this investigation suggest that severe stenoses inhibit wall motion, resulting in higher blood velocities and higher peak wall shear stress, and localization of hoop stress. These factors may contribute to further development and rupture of plaques.

Multiphysics Simulation of Blood Flow and LDL Transport in a Porohyperelastic Arterial Wall Model

Atherosclerosis localizes at a bend and∕or bifurcation of an artery, and low density lipoproteins (LDL) accumulate in the intima. Hemodynamic factors are known to affect this localization and LDL accumulation, but the details of the process remain unknown. It is thought that the LDL concentration will be affected by the filtration flow, and that the velocity of this flow will be affected by deformation of the arterial wall. Thus, a coupled model of a blood flow and a deformable arterial wall with filtration flow would be invaluable for simulation of the flow field and concentration field in sequence. However, this type of highly coupled interaction analysis has not yet been attempted. Therefore, we performed a coupled analysis of an artery with multiple bends in sequence. First, based on the theory of porous media, we modeled a deformable arterial wall using a porohyperelastic model (PHEM) that was able to express both the filtration flow and the viscoelastic behavior of the living tissue, and simulated a blood flow field in the arterial lumen, a filtration flow field and a displacement field in the arterial wall using a fluid-structure interaction (FSI) program code by the finite element method (FEM). Next, based on the obtained results, we further simulated LDL transport using a mass transfer analysis code by the FEM. We analyzed the PHEM in comparison with a rigid model. For the blood flow, stagnation was observed downward of the bends. The direction of the filtration flow was only from the lumen to the wall for the rigid model, while filtration flows from both the wall to the lumen and the lumen to the wall were observed for the PHEM. The LDL concentration was high at the lumen∕wall interface for both the PHEM and rigid model, and reached its maximum value at the stagnation area. For the PHEM, the maximum LDL concentration in the wall in the radial direction was observed at the position of 3% wall thickness from the lumen∕wall interface, while for the rigid model, it was observed just at the lumen∕wall interface. In addition, the peak LDL accumulation area of the PHEM moved about according to the pulsatile flow. These results demonstrate that the blood flow, arterial wall deformation, and filtration flow all affect the LDL concentration, and that LDL accumulation is due to stagnation and the presence of filtration flow. Thus, FSI analysis is indispensable.

Simulations of congenital septal defect closure and reactivity testing in patient-specific models of the pediatric pulmonary vasculature: A 3D numerical study with fluid-structure interaction

Clinical imaging methods are highly effective in the diagnosis of vascular pathologies, but they do not currently provide enough detail to shed light on the cause or progression of such diseases, and would be hard pressed to foresee the outcome of surgical interventions. Greater detail of and prediction capabilities for vascular hemodynamics and arterial mechanics are obtained here through the coupling of clinical imaging methods with computational techniques. Three-dimensional, patient-specific geometric reconstructions of the pediatric proximal pulmonary vasculature were obtained from x-ray angiogram images and meshed for use with commercial computational software. Two such models from hypertensive patients, one with multiple septal defects, the other who underwent vascular reactivity testing, were each completed with two sets of suitable fluid and structural initial and boundary conditions and used to obtain detailed transient simulations of artery wall motion and hemodynamics in both clinically measured and predicted configurations. The simulation of septal defect closure, in which input flow and proximal vascular stiffness were decreased, exhibited substantial decreases in proximal velocity, wall shear stress (WSS), and pressure in the post-op state. The simulation of vascular reactivity, in which distal vascular resistance and proximal vascular stiffness were decreased, displayed negligible changes in velocity and WSS but a significant drop in proximal pressure in the reactive state. This new patient-specific technique provides much greater detail regarding the function of the pulmonary circuit than can be obtained with current medical imaging methods alone, and holds promise for enabling surgical planning.

Quantifying effects of plaque structure and material properties on stress distributions in human atherosclerotic plaques using 3D FSI models

BACKGROUND

Atherosclerotic plaques may rupture without warning and cause acute cardiovascular syndromes such as heart attack and stroke. Methods to assess plaque vulnerability noninvasively and predict possible plaque rupture are urgently needed.

METHOD

MRI-based three-dimensional unsteady models for human atherosclerotic plaques with multi-component plaque structure and fluid-structure interactions are introduced to perform mechanical analysis for human atherosclerotic plaques.

RESULTS

Stress variations on critical sites such as a thin cap in the plaque can be 300% higher than that at other normal sites. Large calcification block considerably changes stress/strain distributions. Stiffness variations of plaque components (50% reduction or 100% increase) may affect maximal stress values by 20-50%. Plaque cap erosion causes almost no change on maximal stress level at the cap, but leads to 50% increase in maximal strain value.

CONCLUSIONS

Effects caused by atherosclerotic plaque structure, cap thickness and erosion, material properties, and pulsating pressure conditions on stress/strain distributions in the plaque are quantified by extensive computational case studies and parameter evaluations. Computational mechanical analysis has good potential to improve accuracy of plaque vulnerability assessment.

Local maximal stress hypothesis and computational plaque vulnerability index for atherosclerotic plaque assessment

It is believed that atherosclerotic plaque rupture may be related to maximal stress conditions in the plaque. More careful examination of stress distributions in plaques reveals that it may be the local stress/strain behaviors at critical sites such as very thin plaque cap and locations with plaque cap weakness that are more closely related to plaque rupture risk. A "local maximal stress hypothesis" and a stress-based computational plaque vulnerability index (CPVI) are proposed to assess plaque vulnerability. A critical site selection (CSS) method is proposed to identify critical sites in the plaque and critical stress conditions which are be used to determine CPVI values. Our initial results based on 34 2D MRI slices from 14 human coronary plaque samples indicate that CPVI plaque assessment has an 85% agreement rate (91% if the square root of stress values is used) with assessment given by histopathological analysis. Large-scale and long-term patient studies are needed to further validate our findings for more accurate quantitative plaque vulnerability assessment.

Quantitative assessment of coronary artery plaque vulnerability by high-resolution magnetic resonance imaging and computational biomechanics: a pilot study ex vivo

The risk of atherosclerotic plaque disruption is thought to be closely related to plaque composition and rupture triggers such as external mechanical forces. The purpose of this study was to integrate MR imaging and computational techniques for the quantification of plaque vulnerability with morphologic data and biomechanical stress/strain distributions that were all based on high-resolution MR images of coronary artery plaque specimens ex vivo. Twenty-two coronary artery plaque specimens were selectively collected from 10 cadavers. Multislice T(2)-weighted spin echo images were acquired with a resolution of 100 x 100 microm(2). Histopathological images were used as the gold standard for the identification of plaque components and vulnerability. Plaque components were classified on MR images, and the stress/strain components were calculated with a two-dimensional computational model with fluid-structure interactions. As expected, vulnerable plaques appeared to result from a large lipid pool, a thin fibrous cap, and some critical stress/strain conditions. An empiric vulnerability marker was derived and was closely related to the vulnerability score that was determined through pathologic examination. Noninvasive quantification of the MR contrast and mechanical properties of plaque may provide a comprehensive biomarker for the assessment of vulnerability of atherosclerotic plaques.

Computational fluid dynamics modeling of the upper airway of children with obstructive sleep apnea syndrome in steady flow

Computational fluid dynamic (CFD) analysis was used to model the effect of airway geometry on internal pressure in the upper airway of three children with obstructive sleep apnea syndrome (OSAS), and three controls. Model geometry was reconstructed from magnetic resonance images obtained during quiet tidal breathing, meshed with an unstructured grid, and solved at normative peak resting flow. The unsteady Reynolds-averaged Navier-Stokes equations were solved with steady flow boundary conditions in inspiration and expiration, using a two-equation low-Reynolds number turbulence model. Model results were validated using an in-vitro scale model, unsteady flow simulation, and reported nasal resistance measurements in children. Pharynx pressure drop strongly correlated to airway area restriction. Inspiratory pressure drop was primarily proportional to the square of flow, consistent with pressure losses due to convective acceleration caused by area restriction. On inspiration, in OSAS pressure drop occurred primarily between the choanae and the region where the adenoids overlap the tonsils (overlap region) due to airway narrowing, rather than in the nasal passages; in controls the majority of pressure drop was in the nasal passages. On expiration, in OSAS the majority of pressure drop occurred between the oropharynx (posterior to the tongue) and overlap region, and local minimum pressure in the overlap region was near atmospheric due to pressure recovery in the anterior nasopharynx. The results suggest that pharyngeal airway shape in children with OSAS significantly affects internal pressure distribution compared to nasal resistance. The model may also help explain regional dynamic airway narrowing during expiration.