An official website of the United States government
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
Effect of the Womersley number on transition to turbulence in pipe flow: An experimental study
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
- Award ID(s):
- 2335760
- PAR ID:
- 10566244
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 36
- Issue:
- 6
- ISSN:
- 1070-6631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)The instability and transition to turbulence and its evolution in pulsatile flows, which involve reverse flows and unsteady flow separations, is the primary focus of this experimental work. A piston driven by a programmable DC servo motor was used to set-up a water flow system and provide the pulsation characteristics. Time-resolved particle image velocimetry data were acquired in a refractive index matching set-up by using a continuous wave laser and a high-frame-rate digital camera. The position of the piston was continuously recorded by a laser proximity sensor. Five different experiments were carried out with Reynolds numbers in the range of 535–4825 and Womersley numbers from 11.91 to 23.82. The non-stationarity of the data was addressed by incorporating trend removal methods involving low- and high-pass filtering of the data, and using empirical mode decomposition together with the relevant Hilbert–Huang transform to determine the intrinsic mode functions. This latter method is more appropriate for nonlinear and non-stationary cases, for which traditional analysis involving classical Fourier decomposition is not directly applicable. It was found that transition to turbulence is a spontaneous event covering the whole near-wall region. The instantaneous vorticity profiles show the development of a large-scale ring-like attached wall vortical layer (WVL) with smaller vortices of higher frequencies than the pulsation frequency superimposed, which point to a shear layer Kelvin–Helmholtz (K–H) type of instability. Inflectional instability leads to flow separation and the formation of a major roll-up structure with the K–H vortices superimposed. This structure breaks down in the azimuthal direction into smaller turbulence patches with vortical content, which appears to be the prevailing structural content of the flow at each investigated Reynolds number ( Re ). At higher Re numbers, the strength and extent of the vortices are larger and substantial disturbances appear in the free stream region of the flow, which are typical of pipe flows at transitional Re numbers. Turbulence appears to be produced at the locations of maximum or minimum vorticity within the attached WVL, in the ridges between the K–H vortices around the separated WVL and the upstream side of the secondary vortex where the flow impinges on the wall. This wall turbulence breaks away into the middle section of the pipe, at approximately $$Re \ge 2200$$ , by strong eruptions of the K–H vortices.more » « less
-
The transition from laminar to turbulent flow is of great interest since it is one of the most difficult and unsolved problems in fluids engineering. The transition processes are significantly important because the transition has a huge impact on almost all systems that come in contact with a fluid flow by altering the mixing, transport, and drag properties of fluids even in simple pipe and channel flows. Generally, in most transportation systems, the transition to turbulence causes a significant increase in drag force, energy consumption, and, therefore, operating cost. Thus, understanding the underlying mechanisms of the laminar-to-turbulent transition can be a major benefit in many ways, especially economically. There have been substantial previous studies that focused on testing the stability of laminar flow and finding the critical amplitudes of disturbances necessary to trigger the transition in various wall-bounded systems, including circular pipes and square ducts. However, there is still no fundamental theory of transition to predict the onset of turbulence. In this study, we perform direct numerical simulations (DNS) of the transition flows from laminar to turbulence in a channel flow. Specifically, the effects of different magnitudes of perturbations on the onset of turbulence are investigated. The perturbation magnitudes vary from 0.001 (0.1%) to 0.05 (5%) of a typical turbulent velocity field, and the Reynolds number is from 5,000 to 40,000. Most importantly, the transition behavior in this study was found to be in good agreement with other reported studies performed for fluid flow in pipes and ducts. With the DNS results, a finite amplitude stability curve was obtained. The critical magnitude of perturbation required to cause transition was observed to be inversely proportional to the Reynolds number for the magnitude from 0.01 to 0.05. We also investigated the temporal behavior of the transition process, and it was found that the transition time or the time required to begin the transition process is inversely correlated with the Reynolds number only for the magnitude from 0.02 to 0.05, while different temporal behavior occurs for smaller perturbation magnitudes. In addition to the transition time, the transition dynamics were investigated by observing the time series of wall shear stress. At the onset of transition, the shear stress experiences an overshoot, then decreases toward sustained turbulence. As expected, the average values of the wall shear stress in turbulent flow increase with the Reynolds number. The change in the wall shear stress from laminar to overshoot was, of course, found to increase with the Reynolds number. More interestingly was the observed change in wall shear stress from the overshoot to turbulence. The change in magnitude appears to be almost insensitive to the Reynolds number and the perturbation magnitude. Because the change in wall shear stress is directly proportional to the pumping power, these observations could be extremely useful when determining the required pumping power in certain flow conditions. Furthermore, the stability curve and wall shear stress changes can be considered robust features for future applications, and ultimately interpreted as evidence of progress toward solving the unresolved fluids engineering problem.more » « less
-
In past experiments, simulations and theoretical analysis, rotation has been shown to dramatically effect the characteristics of turbulent flows, such as causing the mean velocity profile to appear laminar, leading to an overall drag reduction, as well as affecting the Reynolds stress tensor. The axially rotating pipe is an exemplary prototypical model problem that exhibits these complex turbulent flow physics. For this flow, the rotation of the pipe causes a region of turbulence suppression which is particularly sensitive to the rotation rate and Reynolds number. The physical mechanisms causing turbulence suppression are currently not well-understood, and a deeper understanding of these mechanisms is of great value for many practical examples involving swirling or rotating flows, such as swirl generators, wing-tip vortices, axial compressors, hurricanes, etc. In this work, Direct Numerical Simulations (DNS) of rotating turbulent pipe flows are conducted at moderate Reynolds numbers (Re=5300, 11,700, and 19,000) and rotation numbers of N=0 to 3. The main objectives of this work are to firstly quantify turbulence suppression for rotating turbulent pipe flows at different Reynolds numbers as well as study the effects of rotation on turbulence by analyzing the characteristics of the Reynolds stress tensor and the production and dissipation terms of the turbulence budgets.more » « less
-
The drag force on a spherical intruder in dense granular shear flows is studied using discrete element method simulations. Three regimes of the intruder dynamics are observed depending on the magnitude of the drag force (or the corresponding intruder velocity) and the flow inertial number: a fluctuation-dominated regime for small drag forces; a viscous regime for intermediate drag forces; and an inertial (cavity formation) regime for large drag forces. The transition from the viscous regime (linear force-velocity relation) to the inertial regime (quadratic force-velocity relation) depends further on the inertial number. Despite these distinct intruder dynamics, we find a quantitative similarity between the intruder drag in granular shear flows and the Stokesian drag on a sphere in a viscous fluid for intruder Reynolds numbers spanning five orders of magnitude. Beyond this first-order description, a modified Stokes drag model is developed that accounts for the secondary dependence of the drag coefficient on the inertial number and the intruder size and density ratios. When the drag model is coupled with a segregation force model for intruders in dense granular flows, it is possible to predict the velocity of gravity-driven segregation of an intruder particle in shear flow simulations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0210898
