Temporal variations of wall shear stress parameters in intracranial aneurysms—importance of patient-specific inflow waveforms for CFD calculations

2010 ◽  
Vol 152 (8) ◽  
pp. 1391-1398 ◽  
Author(s):  
Christof Karmonik ◽  
Christopher Yen ◽  
Orlando Diaz ◽  
Richard Klucznik ◽  
Robert G. Grossman ◽  
...  
2009 ◽  
Vol 7 (47) ◽  
pp. 967-988 ◽  
Author(s):  
H. Baek ◽  
M. V. Jayaraman ◽  
P. D. Richardson ◽  
G. E. Karniadakis

We investigate the flow dynamics and oscillatory behaviour of wall shear stress (WSS) vectors in intracranial aneurysms using high resolution numerical simulations. We analyse three representative patient-specific internal carotid arteries laden with aneurysms of different characteristics: (i) a wide-necked saccular aneurysm, (ii) a narrower-necked saccular aneurysm, and (iii) a case with two adjacent saccular aneurysms. Our simulations show that the pulsatile flow in aneurysms can be subject to a hydrodynamic instability during the decelerating systolic phase resulting in a high-frequency oscillation in the range of 20–50 Hz, even when the blood flow rate in the parent vessel is as low as 150 and 250 ml min −1 for cases (iii) and (i), respectively. The flow returns to its original laminar pulsatile state near the end of diastole. When the aneurysmal flow becomes unstable, both the magnitude and the directions of WSS vectors fluctuate at the aforementioned high frequencies. In particular, the WSS vectors around the flow impingement region exhibit significant spatio-temporal changes in direction as well as in magnitude.


2016 ◽  
Vol 9 (4) ◽  
pp. 376-380 ◽  
Author(s):  
Waleed Brinjikji ◽  
Bong Jae Chung ◽  
Carlos Jimenez ◽  
Christopher Putman ◽  
David F Kallmes ◽  
...  

BackgroundWhile clinical and angiographic risk factors for intracranial aneurysm instability are well established, it is reasonable to postulate that intra-aneurysmal hemodynamics also have a role in aneurysm instability.ObjectiveTo identify hemodynamic characteristics that differ between radiologically unstable and stable unruptured intracranial aneurysms.Materials and methods12 pairs of unruptured intracranial aneurysms with a 3D rotational angiographic set of images and followed up longitudinally without treatment were studied. Each pair consisted of one stable aneurysm (no change on serial imaging) and one unstable aneurysm (demonstrated growth of at least 1 mm diameter or ruptured during follow-up) of matching size (within 10%) and locations. Patient-specific computational fluid dynamics models were created and run under pulsatile flow conditions. Relevant hemodynamic and geometric variables were calculated and compared between groups using the paired Wilcoxon test.ResultsThe area of the aneurysm under low wall shear stress (low shear stress area (LSA)) was 2.26 times larger in unstable aneurysms than in stable aneurysms (p=0.0499). The mean aneurysm vorticity was smaller by a factor of 0.57 in unstable aneurysms compared with stable aneurysms (p=0.0499). No statistically significant differences in geometric variables or shape indices were found.ConclusionsThis pilot study suggests there may be hemodynamic differences between unstable and stable unruptured cerebral aneurysms. In particular, the area under low wall shear stress was larger in unstable aneurysms. These findings should be considered tentative until confirmed by future larger studies.


BMC Neurology ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Mikołaj Zimny ◽  
Edyta Kawlewska ◽  
Anna Hebda ◽  
Wojciech Wolański ◽  
Piotr Ładziński ◽  
...  

Abstract Background Previously published computational fluid dynamics (CFD) studies regarding intracranial aneurysm (IA) formation present conflicting results. Our study analysed the involvement of the combination of high wall shear stress (WSS) and a positive WSS gradient (WSSG) in IA formation. Methods We designed a case-control study with a selection of 38 patients with an unruptured middle cerebral artery (MCA) aneurysm and 39 non-aneurysmal controls to determine the involvement of WSS, oscillatory shear index (OSI), the WSSG and its absolute value (absWSSG) in aneurysm formation based on patient-specific CFD simulations using velocity profiles obtained from transcranial colour-coded sonography. Results Among the analysed parameters, only the WSSG had significantly higher values compared to the controls (11.05 vs − 14.76 [Pa/mm], P = 0.020). The WSS, absWSSG and OSI values were not significantly different between the analysed groups. Logistic regression analysis identified WSS and WSSG as significant co-predictors for MCA aneurysm formation, but only the WSSG turned out to be a significant independent prognosticator (OR: 1.009; 95% CI: 1.001–1.017; P = 0.025). Significantly more patients (23/38) in the case group had haemodynamic regions of high WSS combined with a positive WSSG near the bifurcation apex, while in the control group, high WSS was usually accompanied by a negative WSSG (14/39). From the analysis of the ROC curve for WSSG, the area under the curve (AUC) was 0.654, with the optimal cut-off value −0.37 Pa/mm. The largest AUC was recognised for combined WSS and WSSG (AUC = 0.671). Our data confirmed that aneurysms tend to form near the bifurcation apices in regions of high WSS values accompanied by positive WSSG. Conclusions The development of IAs is determined by an independent effect of haemodynamic factors. High WSS impacts MCA aneurysm formation, while a positive WSSG mainly promotes this process.


2019 ◽  
Vol 11 (10) ◽  
pp. 999-1003 ◽  
Author(s):  
Michael R Levitt ◽  
Christian Mandrycky ◽  
Ashley Abel ◽  
Cory M Kelly ◽  
Samuel Levy ◽  
...  

ObjectivesTo study the correlation between wall shear stress and endothelial cell expression in a patient-specific, three-dimensional (3D)-printed model of a cerebral aneurysm.Materials and methodsA 3D-printed model of a cerebral aneurysm was created from a patient’s angiogram. After populating the model with human endothelial cells, it was exposed to media under flow for 24 hours. Endothelial cell morphology was characterized in five regions of the 3D-printed model using confocal microscopy. Endothelial cells were then harvested from distinct regions of the 3D-printed model for mRNA collection and gene analysis via quantitative polymerase chain reaction (qPCR.) Cell morphology and mRNA measurement were correlated with computational fluid dynamics simulations.ResultsThe model was successfully populated with endothelial cells, which survived under flow for 24 hours. Endothelial morphology showed alignment with flow in the proximal and distal parent vessel and aneurysm neck, but disorganization in the aneurysm dome. Genetic analysis of endothelial mRNA expression in the aneurysm dome and distal parent vessel was compared with the proximal parent vessels. ADAMTS-1 and NOS3 were downregulated in the aneurysm dome, while GJA4 was upregulated in the distal parent vessel. Disorganized morphology and decreased ADAMTS-1 and NOS3 expression correlated with areas of substantially lower wall shear stress and wall shear stress gradient in computational fluid dynamics simulations.ConclusionsCreating 3D-printed models of patient-specific cerebral aneurysms populated with human endothelial cells is feasible. Analysis of these cells after exposure to flow demonstrates differences in both cell morphology and genetic expression, which correlate with areas of differential hemodynamic stress.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yu Wang ◽  
Junwei Wang ◽  
Jing Peng ◽  
Mingming Huo ◽  
Zhiqiang Yang ◽  
...  

Patients with heart failure (HF) or undergoing cardiogenic shock and percutaneous coronary intervention require short-term cardiac support. Short-term cardiac support using a left ventricular assist device (LVAD) alters the pressure and flows of the vasculature by enhancing perfusion and improving the hemodynamic performance for the HF patients. However, due to the position of the inflow and outflow of the LVAD, the local hemodynamics within the aorta is altered with the LVAD support. Specifically, blood velocity, wall shear stress, and pressure difference are altered within the aorta. In this study, computational fluid dynamics (CFD) was used to elucidate the effects of a short-term LVAD for hemodynamic performance in a patient-specific aorta model. The three-dimensional (3D) geometric models of a patient-specific aorta and a short-term LVAD, Impella CP, were created. Velocity, wall shear stress, and pressure difference in the patient-specific aorta model with the Impella CP assistance were calculated and compared with the baseline values of the aorta without Impella CP support. Impella CP support augmented cardiac output, blood velocity, wall shear stress, and pressure difference in the aorta. The proposed CFD study could analyze the quantitative changes in the important hemodynamic parameters while considering the effects of Impella CP, and provide a scientific basis for further predicting and assessing the effects of these hemodynamic signals on the aorta.


Author(s):  
Navid Freidoonimehr ◽  
Rey Chin ◽  
Anthony C. Zander ◽  
Maziar Arjomandi

Abstract Temporal variations of the coronary arteries during a cardiac cycle are defined as the superposition of the changes in the position, curvature, and torsion of the coronary artery axis markers and the variations in the lumen cross-sectional shape due to the distensible wall motion induced by the pulse pressure and contraction of the myocardium in a cardiac cycle. This review discusses whether the modelling the temporal variations of the coronary arteries is needed for the investigation of the hemodynamics specifically in time critical applications such as a clinical environment. The numerical modellings in the literature which model or disregard the temporal variations of the coronary arteries on the hemodynamic parameters are discussed. The results in the literature show that neglecting the effects of temporal geometric variations is expected to result in about 5\% deviation of the time-averaged pressure drop and wall shear stress values and also about 20\% deviation of the temporal variations of hemodynamic parameters, such as time-dependent wall shear stress and oscillatory shear index. This review study can be considered as a guide for the future studies to outline the conditions in which temporal variations of the coronary arteries can be neglected, while providing a reliable estimation of hemodynamic parameters.


PLoS ONE ◽  
2019 ◽  
Vol 14 (6) ◽  
pp. e0217271 ◽  
Author(s):  
A. M. Moerman ◽  
K. Dilba ◽  
S. Korteland ◽  
D. H. J. Poot ◽  
S. Klein ◽  
...  

2015 ◽  
Vol 8 (8) ◽  
pp. 808-812 ◽  
Author(s):  
Ying Zhang ◽  
Linkai Jing ◽  
Jian Liu ◽  
Chuanhui Li ◽  
Jixing Fan ◽  
...  

ObjectiveTo identify clinical, morphological, and hemodynamic independent characteristic factors that discriminate posterior communicating artery (PCoA) aneurysm rupture status.Methods173 patients with single PCoA aneurysms (108 ruptured, 65 unruptured) between January 2012 and June 2014 were retrospectively collected. Patient-specific models based on their three-dimensional digital subtraction angiography images were constructed and analyzed by a computational fluid dynamic method. All variables were analyzed by univariate analysis and multivariate logistic regression analysis.ResultsTwo clinical factors (younger age and atherosclerosis), three morphological factors (higher aspect ratio, bifurcation type, and irregular shape), and six hemodynamic factors (lower mean and minimum wall shear stress, higher oscillatory shear index, a greater portion of area under low wall shear stress, unstable and complex flow pattern) were significantly associated with PCoA aneurysm rupture. Independent factors characterizing the rupture status were identified as age (OR 0.956, p=0.015), irregular shape (OR 6.709, p<0.001), and minimum wall shear stress (OR 0.001, p=0.038).ConclusionsWe combined clinical, morphological, and hemodynamic characteristics analysis and found the three strongest independent factors for PCoA aneurysm rupture were younger age, irregular shape, and low minimum wall shear stress. This may be useful for guiding risk assessments and subsequent treatment decisions for PCoA aneurysms.


1994 ◽  
Vol 116 (3) ◽  
pp. 294-301 ◽  
Author(s):  
D. A. Steinman ◽  
C. Ross Ethier

The development of intimal hyperplasia at the distal anastomosis is the major cause of long-term bypass graft failure. To evaluate the suspected role of hemodynamic factors in the pathogenesis of distal intimal hyperplasia, an understanding of anastomotic flow patterns is essential. Due to the complexity of arterial flow, model studies typically make simplifying assumptions, such as treating the artery and graft walls as rigid. In the present study this restriction is relaxed to consider the effects of vessel wall distensibility on anastomotic flow patterns. Flow was simulated in an idealized 2-D distensible end-to-side anastomosis model, using parameters appropriate for the distal circulation and assuming a purely elastic artery wall. A novel numerical approach was developed in which the wall velocities are solved simultaneously with the fluid and pressure fields, while the wall displacements are treated via an iterative update. Both the rigid and distensible cases indicated the presence of elevated temporal variations and low average magnitudes of wall shear stress at sites known to be susceptible to the development of intimal hyperplasia. At these same sites, large spatial gradients of wall shear stress were also noted. Comparison between distensible-walled and corresponding rigid-walled simulations showed moderate changes in wall shear stress at isolated locations, primarily the bed, toe and heel. For example, in the case of a distensible geometry and a physiologic pressure waveform, the heel experienced a 38 percent increase in cycle-averaged shear stress, with a corresponding 15 percent reduction in shear stress variability, both relative to the corresponding values in the rigid-walled case. However, other than at these isolated locations, only minor changes in overall wall shear stress patterns were observed. While the physiological implications of such changes in wall shear stress are not known, it is suspected that the effects of wall distensibility are less pronounced than those brought about by changes in arterial geometry and flow conditions.


2012 ◽  
Vol 134 (9) ◽  
Author(s):  
Matthew D. Ford ◽  
Ugo Piomelli

Cerebral aneurysms are a common cause of death and disability. Of all the cardiovascular diseases, aneurysms are perhaps the most strongly linked with the local fluid mechanic environment. Aside from early in vivo clinical work that hinted at the possibility of high-frequency intra-aneurysmal velocity oscillations, flow in cerebral aneurysms is most often assumed to be laminar. This work investigates, through the use of numerical simulations, the potential for disturbed flow to exist in the terminal aneurysm of the basilar bifurcation. The nature of the disturbed flow is explored using a series of four idealized basilar tip models, and the results supported by four patient specific terminal basilar tip aneurysms. All four idealized models demonstrated instability in the inflow jet through high frequency fluctuations in the velocity and the pressure at approximately 120 Hz. The instability arises through a breakdown of the inflow jet, which begins to oscillate upon entering the aneurysm. The wall shear stress undergoes similar high-frequency oscillations in both magnitude and direction. The neck and dome regions of the aneurysm present 180 deg changes in the direction of the wall shear stress, due to the formation of small recirculation zones near the shear layer of the jet (at the frequency of the inflow jet oscillation) and the oscillation of the impingement zone on the dome of the aneurysm, respectively. Similar results were observed in the patient-specific models, which showed high frequency fluctuations at approximately 112 Hz in two of the four models and oscillations in the magnitude and direction of the wall shear stress. These results demonstrate that there is potential for disturbed laminar unsteady flow in the terminal aneurysm of the basilar bifurcation. The instabilities appear similar to the first instability mode of a free round jet.


Sign in / Sign up

Export Citation Format

Share Document