A perturbative model for predicting the high-Reynolds-number behaviour of the streamwise travelling waves technique in turbulent drag reduction

The background of this work is the problem of reducing the aerodynamic turbulent friction drag, which is an important source of energy waste in innumerable technological fields. We develop a theoretical framework aimed at predicting the behaviour of existing drag reduction techniques when used at the large values of Re which are typical of applications. We focus on one recently proposed and very promising technique, which consists in creating at the wall streamwise-travelling waves of spanwise velocity. A perturbation analysis of the Navier-Stokes equations that govern the fluid motion is carried out, for the simplest wall-bounded flow geometry, i.e. the plane channel flow. The streamwise base flow is perturbed by the spanwise time-varying base flow induced by the travelling waves. An asymptotic expansion is then carried out with respect to the velocity amplitude of the travelling wave. The analysis, although based on several assumptions, leads to predictions of drag reduction that agree well with the measurements available in literature and mostly computed through DNS of the full Navier-Stokes equations. New DNS data are produced on purpose in this work to validate our method further. The method is then applied to predict the drag-reducing performance of the streamwise-travelling waves at increasing Re, where comparison data are not available. The current belief, based on a Re-range of about one decade only above the transitional value, that drag reduction obtained at low Re is deemed to decrease as Re is increased is fully confirmed by our results. From a quantitative standpoint, however, our outlook based on several decades of increase in Re is much less pessimistic than other existing estimates, and motivates further, more accurate studies on the present subject

Reynolds-dependence of turbulent skin-friction drag reduction induced by spanwise forcing

This paper examines how increasing the value of the Reynolds number $Re$ affects the ability of spanwise-forcing techniques to yield turbulent skin-friction drag reduction. The considered forcing is based on the streamwise-travelling waves of spanwise wall velocity (Quadrio {\em et al. J. Fluid Mech.}, vol. 627, 2009, pp. 161--178). The study builds upon an extensive drag-reduction database created with Direct Numerical Simulation of a turbulent channel flow for two, 5-fold separated values of $Re$, namely $Re_\tau=200$ and $Re_\tau=1000$. The sheer size of the database, which for the first time systematically addresses the amplitude of the forcing, allows a comprehensive view of the drag-reducing characteristics of the travelling waves, and enables a detailed description of the changes occurring when $Re$ increases. The effect of using a viscous scaling based on the friction velocity of either the non-controlled flow or the drag-reduced flow is described. In analogy with other wall-based drag reduction techniques, like for example riblets, the performance of the travelling waves is well described by a vertical shift of the logarithmic portion of the mean streamwise velocity profile. Except when $Re$ is very low, this shift remains constant with $Re$, at odds with the percentage reduction of the friction coefficient, which is known to present a mild, logarithmic decline. Our new data agree with the available literature, which is however mostly based on low-$Re$ information and hence predicts a quick drop of maximum drag reduction with $Re$. The present study supports a more optimistic scenario, where for an airplane at flight Reynolds numbers a drag reduction of nearly 30\% would still be possible thanks to the travelling waves

Performance losses of drag-reducing spanwise forcing at moderate values of the Reynolds number

A fundamental problem in the field of turbulent skin-friction drag reduction is to determine the performance of the available control techniques at high values of the Reynolds number $Re$. We consider active, predetermined strategies based on spanwise forcing (oscillating wall and streamwise-traveling waves applied to a plane channel flow), and explore via Direct Numerical Simulations (DNS) up to $Re_\tau=2100$ the rate at which their performance deteriorates as $Re$ is increased. To be able to carry out a comprehensive parameter study, we limit the computational cost of the simulations by adjusting the size of the computational domain in the homogeneous directions, compromising between faster computations and the increased need of time-averaging the fluctuating space-mean wall shear-stress. Our results, corroborated by a few full-scale DNS, suggest a scenario where drag reduction degrades with $Re$ at a rate that varies according to the parameters of the wall forcing. In agreement with already available information, keeping them at their low-$Re$ optimal value produces a relatively quick decrease of drag reduction. However, at higher $Re$ the optimal parameters shift towards other regions of the parameter space, and these regions turn out to be much less sensitive to $Re$. Once this shift is accounted for, drag reduction decreases with $Re$ at a markedly slower rate. If the slightly favorable trend of the energy required to create the forcing is considered, a chance emerges for positive net energy savings also at large values of the Reynolds number.Comment: Revised version: change of title, revised intro, small improvements to figures and tex

The phase-locked mean impulse response of a turbulent channel flow

We describe the first DNS-based measurement of the complete mean response of a turbulent channel flow to small external disturbances. Space-time impulsive perturbations are applied at one channel wall, and the linear response describes their mean effect on the flow field as a function of spatial and temporal separations. The turbulent response is shown to differ from the response a laminar flow with the turbulent mean velocity profile as base flow.Comment: Accepted for publication in Physics of Fluid

The laminar generalized Stokes layer and turbulent drag reduction

This paper considers plane channel flow modified by waves of spanwise velocity applied at the wall and travelling along the streamwise direction. Laminar and turbulent regimes for the streamwise flow are both studied. When the streamwise flow is laminar, it is unaffected by the spanwise flow induced by the waves. This flow is a thin, unsteady and streamwise-modulated boundary layer that can be expressed in terms of the Airy function of the first kind. We name it the generalized Stokes layer because it reduces to the classical oscillating Stokes layer in the limit of infinite wave speed. When the streamwise flow is turbulent, the laminar generalized Stokes layer solution describes well the space-averaged turbulent spanwise flow, provided that the phase speed of the waves is sufficiently different from the turbulent convection velocity, and that the time scale of the forcing is smaller than the life time of the near-wall turbulent structures. Under these conditions, the drag reduction is found to scale with the Stokes layer thickness, which renders the laminar solution instrumental for the analysis of the turbulent flow. A classification of the turbulent flow regimes induced by the waves is presented by comparing parameters related to the forcing conditions with the space and time scales of the turbulent flow.Comment: Accepted for publication on J. Fluid Mec

Numerical simulation of turbulent duct flows with constant power input

The numerical simulation of a flow through a duct requires an externally specified forcing that makes the fluid flow against viscous friction. To this aim, it is customary to enforce a constant value for either the flow rate (CFR) or the pressure gradient (CPG). When comparing a laminar duct flow before and after a geometrical modification that induces a change of the viscous drag, both approaches (CFR and CPG) lead to a change of the power input across the comparison. Similarly, when carrying out the (DNS and LES) numerical simulation of unsteady turbulent flows, the power input is not constant over time. Carrying out a simulation at constant power input (CPI) is thus a further physically sound option, that becomes particularly appealing in the context of flow control, where a comparison between control-on and control-off conditions has to be made. We describe how to carry out a CPI simulation, and start with defining a new power-related Reynolds number, whose velocity scale is the bulk flow that can be attained with a given pumping power in the laminar regime. Under the CPI condition, we derive a relation that is equivalent to the Fukagata--Iwamoto--Kasagi relation valid for CFR (and to its extension valid for CPG), that presents the additional advantage of natively including the required control power. The implementation of the CPI approach is then exemplified in the standard case of a plane turbulent channel flow, and then further applied to a flow control case, where the spanwise-oscillating wall is used for skin friction drag reduction. For this low-Reynolds number flow, using 90% of the available power for the pumping system and the remaining 10% for the control system is found to be the optimum share that yields the largest increase of the flow rate above the reference case, where 100% of the power goes to the pump.Comment: Accepted for publication in J. Fluid Mec

Does the choice of the forcing term affect flow statistics in DNS of turbulent channel flow?

We seek possible statistical consequences of the way a forcing term is added to the Navier--Stokes equations in the Direct Numerical Simulation (DNS) of incompressible channel flow. Simulations driven by constant flow rate, constant pressure gradient and constant power input are used to build large databases, and in particular to store the complete temporal trace of the wall-shear stress for later analysis. As these approaches correspond to different dynamical systems, it can in principle be envisaged that these differences are reflect by certain statistics of the turbulent flow field. The instantaneous realizations of the flow in the various simulations are obviously different, but, as expected, the usual one-point, one-time statistics do not show any appreciable difference. However, the PDF for the fluctuations of the streamwise component of wall friction reveals that the simulation with constant flow rate presents lower probabilities for extreme events of large positive friction. The low probability value of such events explains their negligible contribution to the commonly computed statistics; however, the very existence of a difference in the PDF demonstrates that the forcing term is not entirely uninfluential. Other statistics for wall-based quantities (the two components of friction and pressure) are examined; in particular spatio-temporal autocorrelations show small differences at large temporal separations, where unfortunately the residual statistical uncertainty is still of the same order of the observed difference. Hence we suggest that the specific choice of the forcing term does not produce important statistical consequences, unless one is interested in the strongest events of high wall friction, that are underestimated by a simulation run at constant flow rate

Direct numerical simulation of turbulent channel flow over porous walls

We perform direct numerical simulations (DNS) of a turbulent channel flow over porous walls. In the fluid region the flow is governed by the incompressible Navier--Stokes (NS) equations, while in the porous layers the Volume-Averaged Navier--Stokes (VANS) equations are used, which are obtained by volume-averaging the microscopic flow field over a small volume that is larger than the typical dimensions of the pores. In this way the porous medium has a continuum description, and can be specified without the need of a detailed knowledge of the pore microstructure by indipendently assigning permeability and porosity. At the interface between the porous material and the fluid region, momentum-transfer conditions are applied, in which an available coefficient related to the unknown structure of the interface can be used as an error estimate. To set up the numerical problem, the velocity-vorticity formulation of the coupled NS and VANS equations is derived and implemented in a pseudo-spectral DNS solver. Most of the simulations are carried out at $Re_\tau=180$ and consider low-permeability materials; a parameter study is used to describe the role played by permeability, porosity, thickness of the porous material, and the coefficient of the momentum-transfer interface conditions. Among them permeability, even when very small, is shown to play a major role in determining the response of the channel flow to the permeable wall. Turbulence statistics and instantaneous flow fields, in comparative form to the flow over a smooth impermeable wall, are used to understand the main changes introduced by the porous material. A simulations at higher Reynolds number is used to illustrate the main scaling quantities.Comment: Revised version, with additional data and more in-depth analysi

Streamwise oscillation of spanwise velocity at the wall of a channel for turbulent drag reduction

Steady forcing at the wall of a channel flow is studied via DNS to assess its ability of yielding reductions of turbulent friction drag. The wall forcing consists of a stationary distribution of spanwise velocity that alternates in the streamwise direction. The idea behind the forcing builds upon the existing technique of the spanwise wall oscillation, and exploits the convective nature of the flow to achieve an unsteady interaction with turbulence. The analysis takes advantage of the equivalent laminar flow, that is solved analytically to show that the energetic cost of the forcing is unaffected by turbulence. In a turbulent flow, the alternate forcing is found to behave similarly to the oscillating wall; in particular an optimal wavelength is found that yields a maximal reduction of turbulent drag. The energetic performance is significantly improved, with more than 50% of maximum friction saving at large intensities of the forcing, and a net energetic saving of 23% for smaller intensities. Such a steady, wall-based forcing may pave the way to passively interacting with the turbulent flow to achieve drag reduction through a suitable distribution of roughness, designed to excite a selected streamwise wavelength