Finite-Element Modeling of the Hemodynamics of Stented Aneurysms

2004 ◽  
Vol 126 (3) ◽  
pp. 382-387 ◽  
Author(s):  
Gordan R. Stuhne and ◽  
David A. Steinman

Background. Computational fluid dynamics (CFD) simulations are used to analyze the wall shear stress distribution and flow streamlines near the throat of a stented basilar side-wall aneurysm. Previous studies of stented aneurysm flows used low mesh resolution, did not include mesh convergence analyses, and depended upon conformal meshing techniques that apply only to very artificial stent geometries. Method of Approach. We utilize general-purpose computer assisted design and unstructured mesh generation tools that apply in principle to stents and vasculature of arbitrary complexity. A mesh convergence analysis for stented steady flow is performed, varying node spacing near the stent. Physiologically realistic pulsatile simulations are then performed using the converged mesh. Results. Artifact-free resolution of the wall shear stress field on stent wires requires a node spacing of approximately 1/3 wire radius. Large-scale flow features tied to the velocity field are, however, captured at coarser resolution (nodes spaced by about one wire radius or more). Conclusions. Results are consistent with previous work, but our methods yield more detailed insights into the complex flow dynamics. However, routine applications of CFD to anatomically realistic cases still depend upon further development of dedicated algorithms, most crucially to handle geometry definition and mesh generation for complicated stent deployments.

2001 ◽  
Vol 124 (1) ◽  
pp. 44-51 ◽  
Author(s):  
Francis Loth ◽  
Steven A. Jones ◽  
Christopher K. Zarins ◽  
Don P. Giddens ◽  
Raja F. Nassar ◽  
...  

Background : Intimal hyperplastic thickening (IHT) is a frequent cause of prosthetic bypass graft failure. Induction and progression of IHT is thought to involve a number of mechanisms related to variation in the flow field, injury and the prosthetic nature of the conduit. This study was designed to examine the relative contribution of wall shear stress and injury to the induction of IHT at defined regions of experimental end-to-side prosthetic anastomoses. Methods and Results: The distribution of IHT was determined at the distal end-to-side anastomosis of seven canine Iliofemoral PTFE grafts after 12 weeks of implantation. An upscaled transparent model was constructed using the in vivo anastomotic geometry, and wall shear stress was determined at 24 axial locations from laser Doppler anemometry measurements of the near wall velocity under conditions of pulsatile flow similar to that present in vivo. The distribution of IHT at the end-to-side PTFE graft was determined using computer assisted morphometry. IHT involving the native artery ranged from 0.0±0.1 mm to 0.05±0.03 mm. A greater amount of IHT was found on the graft hood (PTFE) and ranged from 0.09±0.06 to 0.24±0.06 mm. Nonlinear multivariable logistic analysis was used to model IHT as a function of the reciprocal of wall shear stress, distance from the suture line, and vascular conduit type (i.e. PTFE versus host artery). Vascular conduit type and distance from the suture line independently contributed to IHT. An inverse correlation between wall shear stress and IHT was found only for those regions located on the juxta-anastomotic PTFE graft. Conclusions: The data are consistent with a model of intimal thickening in which the intimal hyperplastic pannus migrating from the suture line was enhanced by reduced levels of wall shear stress at the PTFE graft/host artery interface. Such hemodynamic modulation of injury induced IHT was absent at the neighboring artery wall.


Author(s):  
Takuya Sawada ◽  
Osamu Terashima ◽  
Yasuhiko Sakai ◽  
Kouji Nagata ◽  
Mitsuhiro Shikida ◽  
...  

The objective of this study is to establish a technique for accurately measuring the wall shear stress in turbulent flows using a micro-fabricated hot-film sensor. Previously, we developed a hot-film sensor with a flexible polyimide-film substrate. This sensor can be attached to curved walls and be used in various situations. Furthermore, the sensor has a 20-μm-wide, heated thin metal film. However, the temporal resolution of this hot-film sensor is not very high owing to its substrate’s high heat capacity. Consequently, its performance is inadequate for measuring the wall shear stress “fluctuations” in turbulent flows. Therefore, we have developed another type of hot-film sensor in which the substrate is replaced with silicon, and a cavity has been introduced under the hot-film for reducing heat loss from the sensor and achieving high temporal resolution. Furthermore, for improving the sensor’s spatial resolution, the width of the hot-film is decreased to 10 μm. The structure of the hot-film’s pattern and the flow-detection mechanism are similar to those of the previous sensor. Experimental results show that new hot-film sensor works as expected and has better temporal resolution than the previous hot-film sensor. As future work, we will measure the wall shear stress for a turbulent wall-jet and discuss the relationship between a large-scale coherent vortex structure and wall shear stress based on data obtained using the new hot-film sensor.


1980 ◽  
Vol 102 (4) ◽  
pp. 456-461 ◽  
Author(s):  
J. H. G. Howard ◽  
S. V. Patankar ◽  
R. M. Bordynuik

A parabolic numerical analysis procedure has been used to predict the flow in a straight, radial rotating channel of rectangular cross-section, chosen as a simple model of an impeller passage. A two equation turbulence model was employed, with alternative modifications, to include the influence of Coriolis force on turbulent kinetic energy. Alternative Coriolis force terms were evaluated by comparisons in a high-aspect-ratio duct with measured velocity, wall shear stress and turbulent viscosity. Secondary velocity predictions were checked with data from a low-aspect-ratio duct where the Coriolis modification of turbulence was found less influential than the secondary flow in the modification of side wall shear stress.


1994 ◽  
Vol 259 ◽  
pp. 191-218 ◽  
Author(s):  
Anthony Wietrzak ◽  
Richard M. Lueptow

Instantaneous streamwise fluctuations of the wall shear stress have been measured using a hot-element probe in a thick axisymmetric turbulent boundary layer on a cylinder aligned parallel to the flow. The measurements were made at a momentum-thickness Reynolds number Rθ = 3050 and a ratio of boundary-layer thickness to cylinder radius of δ/a = 5.7. The ratio of the r.m.s. of the fluctuation to the mean value of the wall shear stress, $\tau_{rms}/\bar{\tau}$, is about 0.32, a value slightly lower than that for recent measurements for flow over a flat plate. The probability density function of the wall shear stress is similar to that for planar wall-bounded flows within experimental error. The power spectral density of the wall shear stress shows that a cylindrical boundary layer contains less energy at lower frequencies and more energy at higher frequencies than other wall-bounded flows. Analysis of simultaneous measurement of the streamwise wall shear stress and the streamwise velocity using VITA and peak detection suggests that transverse curvature has little effect on the near-wall burst–sweep cycle compared to planar wall-bounded flows. The angle of inclination of the structures is similar to that measured for large-scale structures in planar wall-bounded flows. However, measurements of the cross-correlation between the shear stress and the velocity suggest the existence of smaller structures yawed to the axis of the cylinder. The coherence between shear stress and velocity shows a low frequency associated with the inclined structures and a higher frequency associated with the yawed structures. The yawed structures could have an arrowhead or half-arrowhead shape and may be associated with fluid from the outer flow washing over the cylinder.


2004 ◽  
Vol 126 (5) ◽  
pp. 835-843 ◽  
Author(s):  
Hiroyuki Abe ◽  
Hiroshi Kawamura ◽  
Haecheon Choi

Direct numerical simulation of a fully developed turbulent channel flow has been carried out at three Reynolds numbers, 180, 395, and 640, based on the friction velocity and the channel half width, in order to investigate very large-scale structures and their effects on the wall shear-stress fluctuations. It is shown that very large-scale structures exist in the outer layer and that they certainly contribute to inner layer structures at high Reynolds number. Moreover, it is revealed that very large-scale structures exist even in the wall shear-stress fluctuations at high Reynolds number, which are essentially associated with the very large-scale structures in the outer layer.


2007 ◽  
Vol 293 (2) ◽  
pp. H1031-H1037 ◽  
Author(s):  
A. D. Augst ◽  
B. Ariff ◽  
S. A. G. McG. Thom ◽  
X. Y. Xu ◽  
A. D. Hughes

Background: Previous clinical studies have observed relationships between increased intima-media thickness (IMT) in the carotid artery, elevated blood pressure, and low wall shear stress (WSS) calculated from the Poiseuille equation. This study used numerical methods to more accurately determine WSS in the carotid artery and to investigate possible determinants of increased IMT. Methods: IMT [common carotid artery (CCA) and bulb], CCA flow velocity, brachial systolic (SBP) and diastolic blood pressure (DBP), and carotid systolic pressure (cSBP) were measured in 14 healthy subjects (aged 44 ± 16 yr). Flow patterns in the carotid bifurcation were determined by computational fluid dynamics (CFD) based on three-dimensional ultrasound geometry. Instantaneous and time-averaged wall shear stress (WSSav), oscillatory shear index (OSI), and wall shear stress angle gradients (WSSAG) were calculated. Results: IMT was positively related to SBP, DBP, cSBP, and WSSAG and inversely related to WSSav in the CCA. In the bulb, IMT was positively related to SBP and cSBP but was not significantly related to WSSav or WSSAG. IMT was unrelated to OSI in both the CCA and the bulb. Conclusion: Increased carotid artery IMT in healthy subjects with no evidence of focal plaques is primarily a response to elevated pressure.


1986 ◽  
Vol 108 (2) ◽  
pp. 166-172 ◽  
Author(s):  
M. Abdelghany ◽  
R. Eichhorn

Hot film probe measurements of the distribution of the wall shear stress were made for axial flow along a rectangular 3 × 6 array of rods with a pitch to diameter ratio, P/D = 4/3, and a wall to diameter ratio, W/D = 2/3. Measurements were performed on rods at several locations and on two duct side walls at a position 62 hydraulic diameters from the entrance. Local shear stress maxima occur near the largest subchannel flow areas with the lowest maximum local shear stress on rods nearest the sidewalls. Maximum to the minimum shear stress ratio on an individual rod is largest for the corner rod. Side wall maximum local shear stress occurs in the first wall subchannel. Overall friction factors calculated from the wall shear stress measurements agree with those calculated from pressure drop data.


2014 ◽  
Vol 743 ◽  
pp. 202-248 ◽  
Author(s):  
Sébastien Deck ◽  
Nicolas Renard ◽  
Romain Laraufie ◽  
Pierre-Élie Weiss

AbstractA numerical investigation of the mean wall shear stress properties on a spatially developing turbulent boundary layer over a smooth flat plate was carried out by means of a zonal detached eddy simulation (ZDES) technique for the Reynolds number range $3060\leq Re_{\theta }\leq 13\, 650$. Some asymptotic trends of global parameters are suggested. Consistently with previous findings, the calculation confirms the occurrence of very large-scale motions approximately $5\delta $ to $6 \delta $ long which are meandering with a lateral amplitude of $0.3 \delta $ and which maintain a footprint in the near-wall region. It is shown that these large scales carry a significant amount of Reynolds shear stress and their influence on the skin friction, denoted $C_{f,2}$, is revisited through the FIK identity by Fukagata, Iwamoto & Kasagi (Phys. Fluids, vol. 14, 2002, p. L73). It is argued that $C_{f,2}$ is the relevant parameter to characterize the high-Reynolds-number turbulent skin friction since the term describing the spatial heterogeneity of the boundary layer also characterizes the total shear stress variations across the boundary layer. The behaviour of the latter term seems to follow some remarkable self-similarity trends towards high Reynolds numbers. A spectral analysis of the weighted Reynolds stress with respect to the distance to the wall and to the wavelength is provided for the first time to our knowledge and allows us to analyse the influence of the largest scales on the skin friction. It is shown that structures with a streamwise wavelength $\lambda _x >\delta $ contribute to more than $60\, \%$ of $C_{f,2}$, and that those larger than $\lambda _x >2\delta $ still represent approximately $45\, \%$ of $C_{f,2}$.


2019 ◽  
Author(s):  
Haifeng Wang ◽  
Timm Krüger ◽  
Fathollah Varnik

AbstractBlood flow in an artery is a fluid-structure interaction problem. It is widely accepted that aneurysm formation, enlargement and failure are associated with wall shear stress (WSS) which is exerted by flowing blood on the aneurysmal wall. To date, the combined effect of aneurysm size and wall elasticity on intra-aneurysm (IA) flow characteristics, particularly in the case of side-wall aneurysms, is poorly understood. Here we propose a model of three-dimensional viscous flow in a compliant artery containing an aneurysm by employing the immersed boundary-lattice Boltzmann-finite element method. This model allows to adequately account for the elastic deformation of both the blood vessel and aneurysm walls. Using this model, we perform a detailed investigation of the flow through aneurysm under different conditions with a focus on the parameters which may influence the wall shear stress. Most importantly, it is shown in this work that the use of flow velocity as a proxy for wall shear stress is well justified only in those sections of the vessel which are close to the ideal cylindrical geometry. Within the aneurysm domain, however, the correlation between wall shear stress and flow velocity is largely lost due to the complexity of the geometry and the resulting flow pattern. Moreover, the correlations weaken further with the phase shift between flow velocity and transmural pressure. These findings have important implications for medical applications since wall shear stress is believed to play a crucial role in aneurysm rupture.


2018 ◽  
Vol 24 (3) ◽  
pp. 288-296 ◽  
Author(s):  
Gerald J Riccardello ◽  
Abhinav R Changa ◽  
Fawaz Al-Mufti ◽  
I Paul Singh ◽  
Chirag Gandhi ◽  
...  

Objective The natural history intracranial aneurysms (IA) remains poorly understood despite significant morbidity and mortality associated with IA rupture. Hemodynamic impingement resulting in elevations in wall shear stress and wall shear stress gradient (WSSG) has been shown to induce aneurysmal remodeling at arterial bifurcations. We investigate the hemodynamic environment specific to side-wall pre-aneurysmal vasculature. We hypothesize that fluid impingement and secondary flow patterns play a role in side-wall aneurysm initiation. Methods Eight side-wall internal carotid artery aneurysms from the Aneurisk repository were identified. Pre-aneurysmal vasculature was algorithmically reconstructed. Blood flow was simulated with computational fluid dynamic simulations. An indicator of isolated fluid impingement energy was developed by insetting the vessel surface and calculating the impinging component of the fluid dynamic pressure. Results Isolated fluid impingement was found to be elevated in the area of aneurysm initiation in 8/8 cases. The underlying fluid flow for each area of initiation was found to harbor secondary flow patterns known as Dean’s vortices, the result of changes in momentum imparted by bends in the internal carotid artery (ICA). Conclusion Isolated fluid impingement and secondary flow patterns may play a major role in the initiation of side-wall aneurysm initiation. We are unable to determine if this role is through direct or indirect mechanisms but hypothesize that elevations in isolated fluid impingement mark areas of cerebral vasculature that are at risk for aneurysm initiation. Thus, this indicator provides vascular locations to focus future study of side-wall aneurysm initiation.


Sign in / Sign up

Export Citation Format

Share Document