The mechanisms driving the transition to turbulence in pulsatile flows are not well understood. Prior studies in this domain have noted the dynamics of this flow regime to depend on the mean Reynolds number, pulsation frequency (i.e., Womersley number), and inflow pulsatile waveform shape. Conflicting findings, particularly regarding the role of the Womersley number on the critical Reynolds number and the development of turbulence, have been reported. The discord has primarily been observed for flows, with Womersley numbers ranging from 4 to 12. Hence, in this work, we use particle image velocimetry to explore the effects of the Womersley number within this 4–12 range on the dynamics of the pulsatile transition. Eighteen test cases were captured using six mean Reynolds numbers (range 800–4200) and five Womersley numbers. Turbulent kinetic energy, turbulence intensity (TI), and phase lag were computed. Our results indicated that the critical Reynolds number was roughly independent of the Womersley number. At high Womersley numbers, the TI trend maintained lower pulsatility, and the flow was observed to mimic a steady transitional flow regime. A plateau of the TI-velocity and TI-acceleration phase lag was observed at a Womersley number of 8, highlighting that this may be the critical value where further increases to the Womersley number do not alter the transition dynamics. Furthermore, this suggests that the phase lag may provide a universal indicator of the specific influence of the Womersley number on transition for a given flow. Overall, these findings elucidate critical details regarding the role of the Womersley number in the transition to turbulence.
more »
« less
This content will become publicly available on July 27, 2026
An Experimental Study of Turbulence Transition Thresholds in Physiologically Relevant Helical Flows
Abstract Helical blood flow, characterized by its spiral motion, is a crucial physiological phenomenon observed in various circulatory structures, including the heart, aorta, vessel bifurcations, umbilical cord, and respiratory system. Despite its importance, a comprehensive understanding of helical flow dynamics is limited. This study aims to investigate the transition from laminar to turbulent flow in helical tubes under both steady and pulsatile conditions. An experimental setup was developed, featuring a closed-loop system with a pulsatile flow generator and steady pump, producing sinusoidal waveforms with adjustable frequencies and amplitudes. Five helical tube models, varying in curvature radius and torsion, were fabricated with a fixed inner diameter of 15 mm. High-frequency pressure transducers, ultrasonic flow sensors, and Laser Doppler Velocimetry (LDV) were employed to visualize flow fields and measure turbulence kinetic energy (TKE). Results indicate that helical tube geometry significantly impacts turbulence transitions. Specifically, larger curvature radii stabilize the flow and reduce turbulence downstream, while smaller radii lead to earlier transitions to turbulence. Under steady flow conditions, the critical Reynolds number (Re) for turbulence onset was found to be around 2300 upstream, similar to straight pipes, but significantly higher turbulence levels were observed downstream. Pulsatile flow with high pulsatility indices (PI > 3) near transitional Re markedly increases turbulence, particularly downstream, with large bursts of turbulence occurring during the deceleration phase. These findings enhance the understanding of helical flow dynamics and have implications for biomedical device design, diagnostics, and treatments for circulatory disorders.
more »
« less
- Award ID(s):
- 2342517
- PAR ID:
- 10654575
- Publisher / Repository:
- American Society of Mechanical Engineers
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Low-inertia pulsatile flows in highly distensible viscoelastic vessels exist in many biological and engineering systems. However, many existing works focus on inertial pulsatile flows in vessels with small deformations. As such, here we study the dynamics of a viscoelastic tube at large deformation conveying low-Reynolds-number oscillatory flow using a fully coupled fluid–structure interaction computational model. We focus on a detailed study of the effect of wall (solid) viscosity and oscillation frequency on tube deformation, flow rate, phase shift and hysteresis, as well as the underlying flow physics. We find that the general behaviour is dominated by an elastic flow surge during inflation and a squeezing effect during deflation. When increasing the oscillation frequency, the maximum inlet flow rate increases and tube distention decreases, whereas increasing solid viscosity causes both to decrease. As the oscillation frequency approaches either$$0$$(quasi-steady inflation cycle) or$$\infty$$(steady flow), the behaviours of tubes with different solid viscosities converge. Our results suggest that deformation and flow rate are most affected in the intermediate range of solid viscosity and oscillation frequency. Phase shifts of deformation and flow rate with respect to the imposed pressure are analysed. We predict that the phase shifts vary throughout the oscillation; while the deformation always lags the imposed pressure, the flow rate may either lead or lag depending on the parameter values. As such, the flow rate shows hysteresis behaviour that traces either a clockwise or counterclockwise curve, or a mix of both, in the pressure–flow rate space. This directional change in hysteresis is fully characterised here in the appropriate parameter space. Furthermore, the hysteresis direction is shown to be predicted by the signs of the flow rate phase shifts at the crest and trough of the oscillation. A distinct change in the tube dynamics is also observed at high solid viscosity which leads to global or ‘whole-tube’ motion that is absent in purely elastic tubes.more » « less
-
Fritts et al. (J. Fluid Mech., vol. xx, 2022, xx) describe a direct numerical simulation of interacting Kelvin–Helmholtz instability (KHI) billows arising due to initial billow cores that exhibit variable phases along their axes. Such KHI exhibit strong ‘tube and knot’ dynamics identified in early laboratory studies by Thorpe ( Geophys. Astrophys. Fluid Dyn. , vol. 34, 1985, pp. 175–199). Thorpe ( Q.J.R. Meteorol. Soc. , vol. 128, 2002, pp. 1529–1542) noted that these dynamics may be prevalent in the atmosphere, and they were recently identified in atmospheric observations at high altitudes. Tube and knot dynamics were found by Fritts et al. ( J. Fluid. Mech. , 2022) to drive stronger and faster turbulence transitions than secondary instabilities of individual KH billows. Results presented here reveal that KHI tube and knot dynamics also yield energy dissipation rates $$\sim$$ 2–4 times larger as turbulence arises and that remain $$\sim$$ 2–3 times larger to later stages of the flow evolution, compared with those of secondary convective instabilities (CI) and secondary KHI accompanying KH billows without tube and knot influences. Elevated energy dissipation rates occur due to turbulence transitions by tube and knot dynamics arising on much larger scales than secondary CI and KHI where initial KH billows are misaligned. Tube and knot dynamics also excite large-scale Kelvin ‘twist waves’ that cause vortex tube and billow core fragmentation, more energetic cascades of similar interactions to smaller scales and account for the strongest energy dissipation events accompanying such KH billow evolutions.more » « less
-
Accurately identifying the onset of transitional and turbulent flow within any pipe flow environment is of great interest. Most often, the critical Reynolds number (Re) is used to pinpoint the onset of turbulence. However, the critical Re is known to be highly variable, depending on the specifics of the flow system. Thus, for flows (e.g., blood flows), where only one realization (i.e., one mean Re) exists, the presence of transitional and turbulent flow behaviors cannot be accurately determined. In this work, we aim to address this by evaluating the extent to which instantaneous time-frequency (TF)-based analysis of the fluctuating velocity field can be used to evaluate the onset of transitional and turbulent flow regimes. Because current TF analysis methods are not suitable for this, we propose a novel “wavelet-Hilbert time-frequency” (WHTF) method, which we validate herein. Using the WHTF method, we analyzed the instantaneous dominant frequency of three planar particle image velocimetry-captured pipe flows, which included one steady and two pulsatile with Womersley numbers of 4 and 12. For each case, data were captured at Re's spanning 800–4500. The instantaneous dominant frequency analysis of these flows revealed that the magnitude, size, and coherence of two-dimensional spatial frequency structures were uniquely different across flow regimes. Specifically, the transitional regime maintained the most coherent, but lowest magnitude frequency structures, while the laminar regime had the highest magnitude, lowest coherence, and smallest frequency structures. Overall, this study demonstrates the efficacy of TF-based metrics for characterizing the progression of transition and turbulent flow development.more » « less
-
The interactions between fluid flow and structural components of collapsible tubes are representative of those in several physiological systems. Although extensively studied, there exists a lack of characterization of the three-dimensionality in the structural deformations of the tube and its influence on the flow field. This experimental study investigates the spatio-temporal relationship between 3D tube geometry and the downstream flow field under conditions of fully open, closed, and slamming-type oscillating regimes. A methodology is implemented to simultaneously measure three-dimensional surface deformations in a collapsible tube and the corresponding downstream flow field. Stereophotogrammetry was used to measure tube deformations, and simultaneous flow field measurements included pressure and planar Particle Image Velocimetry (PIV) data downstream of the collapsible tube. The results indicate that the location of the largest collapse in the tube occurs close to the downstream end. In the oscillating regime, sections of the tube downstream of the largest mean collapse experience the largest oscillations in the entire tube that are completely coherent and in phase. At a certain streamwise distance upstream of the largest collapse, a switch in the direction of oscillations occurs with respect to those downstream. Physically, when the tube experiences constriction downstream of the location of the largest mean collapse, this causes the accumulation of fluid and build-up of pressure in the upstream regions and an expansion of these sections. Fluctuations in the downstream flow field are significantly influenced by tube fluctuations along the minor axes. The fluctuations in the downstream flowfield are influenced by the propagation of disturbances due to oscillations in tube geometry, through the advection of fluid through the tube. Further, the manifestation of the LU-type pressure fluctuations is found to be due to the variation in the propagation speed of the disturbances during the different stages within a period of oscillation of the tube.more » « less
An official website of the United States government
