An official website of the United States government
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
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
- 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
-
-
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
-
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
-
Nguyen, Thao Vicky; Ethier, C Ross (Ed.)Abstract Cell-laden, scaffold-based tissue engineering methods have been successfully utilized for the treatment of bone fractures and diseases, caused by factors such as trauma, tumors, congenital anomalies, and aging. In such methods, the rate of scaffold biodegradation, transport of nutrients and growth factors, as well as removal of cell metabolic wastes at the site of injury are critical fluid-dynamics factors, affecting cell proliferation and ultimately tissue regeneration. Therefore, there is a critical need to identify the underlying material transport mechanisms and factors associated with cell-seeded, scaffold-based bone tissue engineering. The overarching goal of this study is to contribute to patient-specific, clinical treatment of bone pathology. The overall objective of the work is to establish computational fluid dynamics (CFD) models: (i) to identify the consequential mechanisms behind internal and external material transport through/over porous bone scaffolds designed based on the principles of triply periodic minimal surfaces (TPMS) and (ii) to identify TPMS designs with optimal geometry and flow characteristics for the treatment of bone fractures in clinical practice. In this study, advanced CFD models were established based on ten TPMS scaffold designs for (i) single-unit internal flow analysis, (ii) single-unit external flow analysis, and (iii) cubic, full-scaffold external flow analysis, where the geometry of each design was parametrically created. The influence of several design parameters, such as surface representation iteration, wall thickness, and pore size on geometry accuracy as well as computation time, was investigated in order to obtain computationally efficient and accurate CFD models. The fluid properties (such as density and dynamic viscosity) as well as the boundary conditions (such as no-slip condition, inlet flow velocity, and pressure outlet) of the CFD models were set based on clinical/research values reported in the literature, according to the fundamentals of internal and external Newtonian flow modeling. The main fluid characteristics influential in bone regeneration, including flow velocity, flow pressure, and wall shear stress (WSS), were analyzed to observe material transport internally through and externally over the TPMS scaffold designs. Regarding the single-unit internal flow analysis, it was observed that P.W. Hybrid and Neovius designs had the highest level of not only flow pressure but also WSS. This can be attributed to their relatively flat surfaces when compared to the rest of the TPMS designs. Schwarz primitive (P) appeared to have the lowest level of flow pressure and WSS (desirable for development of bone tissues) due to its relatively open channels allowing for more effortless fluid transport. An analysis of streamline velocity exhibited an increase in velocity togther with a depiction of potential turbulent motion along the curved sections of the TPMS designs. Regarding the single-unit external flow analysis, it was observed that Neovius and Diamond yielded the highest level of flow pressure and WSS, respectively, while Schwarz primitive (P) similarly had a relatively low level of flow pressure and WSS suitable for bone regeneration. Besides, pressure buildup was observed within the inner channels of almost all the TPMS designs due to flow resistance and the intrinsic interaction between the fluid flow and the scaffold walls. Regarding the cubic (full-scaffold) external flow analysis, the Diamond and Schwarz gyroid (G) designs appeared to have a relatively high level of both flow pressure and WSS, while Schwarz primitive (P) similarly yielded a low level of flow pressure and WSS. Overall, the outcomes of this study pave the way for optimal design and fabrication of complex, bone-like tissues with desired material transport properties for cell-laden, scaffold-based treatment of bone fractures.more » « less
-
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
