skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 8:00 PM ET on Friday, March 21 until 8:00 AM ET on Saturday, March 22 due to maintenance. We apologize for the inconvenience.


Title: Direct numerical simulation of a pulsatile flow in a stenotic channel using immersed boundary method
Abstract

A three‐dimensional direct numerical simulation model coupled with the immersed boundary method has been developed to simulate a pulsatile flow in a planar channel with single and double one‐sided semicircular constrictions. For relevance to blood flow in large arteries, simulations have been performed at Reynolds numbers of 750 and 1000. Flow physics and resultant wall shear stress (WSS)‐based hemodynamic parameters are presented. The instantaneous vortex dynamics, mean flow characteristics, and turbulent energy spectra are evaluated for flow physics. Subsequently, three WSS‐based parameters, namely the time‐averaged WSS, oscillatory shear index, and relative residence time, are calculated over the stenotic wall and correlated with flow physics to identify the regions prone to atherosclerotic plaque progression. Results show that the double stenotic channel leads to high‐intensity and broadband turbulent characteristics downstream, promoting critical values of the WSS‐based parameters in the post‐stenotic areas. In addition, the inter‐space area between two stenoses displays multiple strong recirculations, making this area highly prone to atherosclerosis progression. The effect of stenosis degree on the WSS‐based parameters is studied up to 60% degree. As the degree of occlusion is increased, larger regions are involved with the nonphysiological ranges of the WSS‐based parameters.

 
more » « less
Award ID(s):
1936065
PAR ID:
10361542
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Engineering Reports
Volume:
4
Issue:
1
ISSN:
2577-8196
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Spiral/helical forms of blood flow have been observed in large arteries of the cardiovascular system, but their benefits remain underappreciated. Spiral flow has been postulated to improve near‐wall washout, promoting anti‐atherothrombotic conditions. This research aims to study the washout characteristics of spiral flow, specifically, its ability to increase velocity and wall shear stress (WSS) in atherothrombotic‐prone regions. Using 1.2 cm diameter angled test‐conduits (45°, 90°, 135°) with known recirculation/stasis regions at the bend corners, spiral flow washout potential was evaluated in terms of low velocity and low WSS. Two sub‐studies were conducted: the first utilized a spiral flow‐inducing device to enable qualitative analysis of washout‐potential in both computational fluid dynamic (CFD) simulations and benchtop ultrasound visualization; the second used CFD to study the impact of several induced helical wavelengths on the conduit‐dependent recirculation/stasis zones. Physical models of the angled conduits and spiral flow‐inducer were 3D‐printed to facilitate ultrasound visualization. Compared to straight flow, spiral flow generated by the flow‐inducer significantly cleared the recirculation/stasis zones at the corners of the angled conduits. CFD simulations demonstrated that past a geometry‐dependent threshold, increased helical content improved washout, denoted by decreased regions of low velocity and low WSS. Overall, spiral flow markedly improved washout in difficult to reach areas in the angled conduits. This has several important clinical implications: spiral flow shows great promise in reducing blood‐transport‐related complications and can be used to enhance the performance of future medical devices (eg grafts, mechanical circulatory support devices, hemodialysis access ports).

     
    more » « less
  2. Direct numerical simulations are performed for incompressible, turbulent channel flow over a smooth wall and different sinusoidal wall roughness configurations at a constant $Re_\tau = 720$ . Sinusoidal walls are used to study the effects of well-defined geometric features of roughness-amplitude, $a$ , and wavelength, $\lambda$ , on the flow. The flow in the near-wall region is strongly influenced by both $a$ and $\lambda$ . Establishing appropriate scaling laws will aid in understanding the effects of roughness and identifying the relevant physical mechanisms. Using inner variables and the roughness function to scale the flow quantities provides support for Townsend's hypothesis, but inner scaling is unable to capture the flow physics in the near-wall region. We provide modified scaling relations considering the dynamics of the shear layer and its interaction with the roughness. Although not a particularly surprising observation, this study provides clear evidence of the dependence of flow features on both $a$ and $\lambda$ . With these relations, we are able to collapse and/or align peaks for some flow quantities and, thus, capture the effects of surface roughness on turbulent flows even in the near-wall region. The shear-layer scaling supports the hypothesis that the physical mechanisms responsible for turbulent kinetic energy production in turbulent flows over rough walls are greatly influenced by the shear layer and its interaction with the roughness elements. Finally, a semiempirical model is developed to predict the contribution of pressure and skin friction drag on the roughness element based purely on its geometric parameters and the corresponding shear-layer velocity scale. 
    more » « less
  3. Highly-resolved, direct numerical simulations of turbulent channel flows with sub- Kolmogorov grid resolution are performed to investigate the characteristics of wall-bounded turbulent flows in the presence of sinusoidal wall waviness. The wall waviness serves as a simplified model to study the effects of well-defined geometric parameters of roughness on the characteristics of wall-bounded turbulent flows. In this study, a two-dimensional wave profile with steepness ranging from 0.06 to 0.25 and wave amplitudes ranging from 9 to 36 wall units were considered. For the smooth and wavy-wall simulations, the Reynolds number based on the friction velocity was kept constant. To study the effects of wave amplitude and wavelength on turbulence, two-dimensional time and spanwise averaged distributions of the mean flow, turbulent kinetic energy, and Reynolds stresses as well as turbulent kinetic energy production and dissipation are examined. Furthermore, in order to provide a more direct comparison with the smooth-wall turbulent channel flow one-dimensional pro- files of these quantities are computed by averaging them over one wavelength of the wave profile. A strong effect of the wall-waviness and, in particular, the wave amplitude and wavelength on the characteristics of the turbulence was obtained. Wall waviness mainly affected the inner flow region while all recorded turbulent statistics collapsed in the outer flow region. Significant reductions in turbulent kinetic energy, production and dissipation were obtained with increasing wave amplitudes when reported in inner scale. While production is lower for all wavy wall cases considered here in comparison to the smooth wall, reducing the wavelength caused an increase in production and a decrease in dissipation. 
    more » « less
  4. null (Ed.)
    Recent direct numerical simulations (DNS) and experiments in turbulent channel flow have found intermittent low- and high-drag events in Newtonian fluid flows, at Reτ=uτh/ν between 70 and 100, where uτ, h and ν are the friction velocity, channel half-height and kinematic viscosity, respectively. These intervals of low-drag and high-drag have been termed “hibernating” and “hyperactive”, respectively, and in this paper, a further investigation of these intermittent events is conducted using experimental and numerical techniques. For experiments, simultaneous measurements of wall shear stress and velocity are carried out in a channel flow facility using hot-film anemometry (HFA) and laser Doppler velocimetry (LDV), respectively, for Reτ between 70 and 250. For numerical simulations, DNS of a channel flow is performed in an extended domain at Reτ = 70 and 85. These intermittent events are selected by carrying out conditional sampling of the wall shear stress data based on a combined threshold magnitude and time-duration criteria. The use of three different scalings (so-called outer, inner and mixed) for the time-duration criterion for the conditional events is explored. It is found that if the time-duration criterion is kept constant in inner units, the frequency of occurrence of these conditional events remain insensitive to Reynolds number. There exists an exponential distribution of frequency of occurrence of the conditional events with respect to their duration, implying a potentially memoryless process. An explanation for the presence of a spike (or dip) in the ensemble-averaged wall shear stress data before and after the low-drag (or high-drag) events is investigated. During the low-drag events, the conditionally-averaged streamwise velocities get closer to Virk’s maximum drag reduction (MDR) asymptote, near the wall, for all Reynolds numbers studied. Reynolds shear stress (RSS) characteristics during these conditional events are investigated for Reτ = 70 and 85. Except very close to the wall, the conditionally-averaged RSS is higher than the time-averaged value during the low-drag events. 
    more » « less
  5. Abstract

    The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC “footprints.” Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.

     
    more » « less