skip to main content


Title: Role of flow reversals in transition to turbulence and relaminarization of pulsatile flows
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
Award ID(s):
1803845
NSF-PAR ID:
10278476
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
917
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Drag for wall-bounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wall-normal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number$Re_\tau =1000$. The net transfer from the wall of spanwise vorticity created by downstream pressure drop is due to two large opposing fluxes, one which is ‘down-gradient’ or outward from the wall, where most vorticity concentrates, and the other which is ‘up-gradient’ or toward the wall and acting against strong viscous diffusion in the near-wall region. We present evidence that the up-gradient/down-gradient transport occurs by a mechanism of correlated inflow/outflow and spanwise vortex stretching/contraction that was proposed by Lighthill. This mechanism is essentially Lagrangian, but we explicate its relation to the Eulerian anti-symmetric vorticity flux tensor. As evidence for the mechanism, we study (i) statistical correlations of the wall-normal velocity and of wall-normal flux of spanwise vorticity, (ii) vorticity flux cospectra identifying eddies involved in nonlinear vorticity transport in the two opposing directions and (iii) visualizations of coherent vortex structures which contribute to the transport. The ‘D-type’ vortices contributing to down-gradient transport in the log layer are found to be attached, hairpin-type vortices. However, the ‘U-type’ vortices contributing to up-gradient transport are detached, wall-parallel, pancake-shaped vortices with strong spanwise vorticity, as expected by Lighthill's mechanism. We discuss modifications to the attached eddy model and implications for turbulent drag reduction.

     
    more » « less
  3. Rotating and swirling turbulence comprises an important class of flows, not only due to the complex physics that occur, but also due to their relevance to many engineering applications, such as combustion, cyclone separation, mixing, etc. In these types of flows, rotation strongly affects the characteristics and structure of turbulence. However, the underlying turbulent flow phenomena are complex and currently not well understood. The axially rotating pipe is an exemplary prototypical model problem that exhibits these complex turbulent flow physics. By examining the complex interaction of turbulent structures within rotating turbulent pipe flow, insight can be gained into the behavior of rotating flows relevant to engineering applications. Direct numerical simulations are conducted at a bulk Reynolds number up to Re_D = 19,000 with rotation numbers ranging from N = 0 to 3. Coherence analysis, including Proper Orthogonal Decomposition and Dynamic Mode Decomposition, are used to identify the relevant (highest energy) modes of the flow. Studying the influence of these modes on turbulent statistics (i.e. mean statistics, Reynolds stresses, turbulent kinetic energy, and turbulent kinetic energy budgets) allow for a deeper understanding of the effects of coherent turbulent flow structures in rotating flows. 
    more » « less
  4. One of the key factors in simulating realistic wall-bounded flows at high Reynolds numbers is the selection of an appropriate turbulence model for the steady Reynolds Averaged Navier–Stokes equations (RANS) equations. In this investigation, the performance of several turbulence models was explored for the simulation of steady, compressible, turbulent flow on complex geometries (concave and convex surface curvatures) and unstructured grids. The turbulence models considered were the Spalart–Allmaras model, the Wilcox k- ω model and the Menter shear stress transport (SST) model. The FLITE3D flow solver was employed, which utilizes a stabilized finite volume method with discontinuity capturing. A numerical benchmarking of the different models was performed for classical Computational Fluid Dynamic (CFD) cases, such as supersonic flow over an isothermal flat plate, transonic flow over the RAE2822 airfoil, the ONERA M6 wing and a generic F15 aircraft configuration. Validation was performed by means of available experimental data from the literature as well as high spatial/temporal resolution Direct Numerical Simulation (DNS). For attached or mildly separated flows, the performance of all turbulence models was consistent. However, the contrary was observed in separated flows with recirculation zones. Particularly, the Menter SST model showed the best compromise between accurately describing the physics of the flow and numerical stability. 
    more » « less
  5. The verification of whether small-scale turbulence is isotropic remains a grand challenge. The difficulty arises because the presence of small-scale anisotropy is tied to the dissipation tensor, whose components require the full three-dimensional information of the flow field in both high spatial and temporal resolution, a condition rarely satisfied in turbulence experiments, especially during field scale measurement of atmospheric turbulence. To circumvent this issue, an intermittency-anisotropy framework is proposed through which we successfully extract the features of small-scale anisotropy from single-point measurements of turbulent time series by exploiting the properties of small-scale intermittency. Specifically, this framework quantifies anisotropy by studying the contrasting effects of burstlike activities on the scalewise production of turbulence kinetic energy between the horizontal and vertical directions. The veracity of this approach is tested by applying it over a range of datasets covering an unprecedented range in the Reynolds numbers (Re≈10^3–10^6), sampling frequencies (10 kHz to 10 Hz), surface conditions (aerodynamically smooth surfaces to typical grasslands to forest canopies), and flow types (channel flows, boundary-layer flows, atmospheric flows, and flows over forest canopies). For these diverse datasets, the findings indicate that the effects of small-scale anisotropy persists up to the integral scales of the streamwise velocity fluctuations and there exists a universal relationship to predict this anisotropy from the two-component state of the Reynolds stress tensor. This relationship is important towards the development of next-generation closure models of wall turbulence by incorporating the effects of anisotropy at smaller scales of the flow. 
    more » « less