skip to main content

Title: Low- and High-Drag Intermittencies in Turbulent Channel Flows
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
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius$R$and axial length$10{\rm \pi} R$at friction Reynolds numbers up to$Re_\tau =5200$using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor$\lambda$overshoots by$2\,\%$and undershoots by$0.6\,\%$the Prandtl friction law at low and high$Re$ranges, respectively. In addition,$\lambda$in our results is slightly higher than in Pirozzoliet al.(J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichiet al.(Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to$y^+=250$, has not yet developed a plateau farther away from the wall in the pipes even for the$Re_\tau =5200$cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with$Re_\tau$– are lower in the pipe than in the channel, but the difference decreases with increasing$Re_\tau$. While the wall value is slightly lower in the channel than in the pipe at the same$Re_\tau$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a$k^{-1}$dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence.

    more » « less
  3. null (Ed.)
    Abstract Coherent structures are critical for controlling turbulent boundary layers due to their roles in momentum and heat transfer in the flow. Turbulent coherent structures can be detected by measuring wall shear stresses that are footprints of coherent structures. In this study, wall shear stress fluctuations were measured simultaneously in a zero pressure gradient turbulent boundary layer using two house-made wall shear stress probes aligned in the spanwise direction. The wall shear stress probe consisted of two hot-wires on the wall aligned in a V-shaped configuration for measuring streamwise and spanwise shear stresses, and their performance was validated in comparison with a direct numerical simulation result. Relationships between measured wall shear stress fluctuations and streamwise velocity fluctuations were analyzed using conditional sampling techniques. The peak detection method and the variable-interval time-averaging (VITA) method showed that quasi-streamwise vortices were inclined toward the streamwise direction. When events were simultaneously detected by the two probes, stronger fluctuations in streamwise velocity were detected, which suggests that stronger coherent structures were detected. In contrast to the former two methods, the hibernating event detection method detects events with lower wall shear stress fluctuations. The ensemble-averaged mean velocity profile of hibernating events was shifted upward compared to the law of the wall, which suggests low drag status of the coherent structures related with hibernating events. These methods suggest significant correlations between wall shear stress fluctuations and coherent structures, which could motivate flow control strategies to fully exploit these correlations. 
    more » « less
  4. Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications. 
    more » « less
  5. Turbulent boundary layers subject to severe acceleration or strong favorable pressure gradient (FPG) are of fundamental and technological importance. Scientifically, they elicit great interest from the points of view of scaling laws, the complex interaction between the outer and inner regions, and the quasi-laminarization phenomenon. Many flows of industrial and technological applications are subject to strong acceleration such as convergent ducts, turbines blades and nozzles. Our recent numerical predictions (J. Fluid Mech., vol. 775, pp. 189-200, 2015) of turbulent boundary layers subject to very strong FPG with high spatial/temporal resolution, i.e. Direct Numerical Simulation (DNS), have shown a meaningful weakening of the Reynolds shear stresses with an evident logarithmic behavior. In the present study, assessment of three different turbulence models (Shear Stress Transport, k-w and Spalart-Allmaras, henceforth SST, k-w and SA, respectively) in Reynolds-averaged Navier-Stokes (RANS) simulations is performed. The main objective is to evaluate the ability of popular turbulence models in capturing the characteristic features present during the quasi-laminarization phenomenon in highly accelerating turbulent boundary layers. Favorable pressure gradient is prescribed by a top converging surface (sink flow) with an approximately constant acceleration parameter of K = 4.0 x 10^(-6). Furthermore, the quasi-laminarization effect on the temperature field is also examined by solving the energy equation and assuming the temperature as a passive scalar. Validation of RANS results is carried out by means of a large DNS dataset. 
    more » « less