Showing 1–44 of 44 results for author: Quadrio, M

Turbulent skin-friction drag reduction via spanwise forcing at high Reynolds number

Authors: Davide Gatti, Maurizio Quadrio, Alessandro Chiarini, Federica Gattere, Sergio Pirozzoli

Abstract: We address the Reynolds-number dependence of the turbulent skin-friction drag reduction induced by streamwise-travelling waves of spanwise wall oscillations. The study relies on direct numerical simulations of drag-reduced flows in a plane open channel at friction Reynolds numbers in the range $1000 \le Re_τ\le 6000$, which is the widest range considered so far in simulations with spanwise forcing… ▽ More We address the Reynolds-number dependence of the turbulent skin-friction drag reduction induced by streamwise-travelling waves of spanwise wall oscillations. The study relies on direct numerical simulations of drag-reduced flows in a plane open channel at friction Reynolds numbers in the range $1000 \le Re_τ\le 6000$, which is the widest range considered so far in simulations with spanwise forcing. Our results corroborate the validity of the predictive model proposed by Gatti & Quadrio, J. Fluid Mech. (2016): regardless of the control parameters, the drag reduction decreases monotonically with $Re$, at a rate that depends on the drag reduction itself and on the skin-friction of the uncontrolled flow. We do not find evidence in support of the results of Marusic et al., Nat. Comm. (2021), which instead report by experiments an increase of the drag reduction with $Re$ in turbulent boundary layers, for control parameters that target low-frequency, outer-scaled motions. Possible explanations for this discrepancy are provided, including obvious differences between open channel flows and boundary layers, and possible limitations of laboratory experiments. △ Less

Submitted 11 September, 2024; originally announced September 2024.

An adjoint-based approach for the surgical correction of nasal septal deviations

Authors: Marcello Macellari, Andrea Schillaci, Umberto Tanzini, Matteo Trimarchi, Maurizio Quadrio

Abstract: Deviations of the septal wall are widespread anatomic anomalies of the human nose; they vary significantly in shape and location, and often cause the obstruction of the nasal airways. When severe, septal deviations need to be surgically corrected by ear-nose-throat (ENT) specialists. Septoplasty, however, has a low success rate, owing to the lack of suitable standardized clinical tools for assessi… ▽ More Deviations of the septal wall are widespread anatomic anomalies of the human nose; they vary significantly in shape and location, and often cause the obstruction of the nasal airways. When severe, septal deviations need to be surgically corrected by ear-nose-throat (ENT) specialists. Septoplasty, however, has a low success rate, owing to the lack of suitable standardized clinical tools for assessing type and severity of obstructions, and for surgery planning. Moreover, the restoration of a perfectly straight septal wall is often impossible and possibly unnecessary. This paper introduces a procedure, based on advanced patient-specific Computational Fluid Dynamics (CFD) simulations, to support ENT surgeons in septoplasty planning. The method hinges upon the theory of adjoint-based optimization, and minimizes a cost function that indirectly accounts for viscous losses. A sensitivity map is computed on the mucosal wall to provide the surgeon with a simple quantification of how much tissue removal at each location would contribute to easing the obstruction. The optimization procedure is applied to three representative nasal anatomies, reconstructed from CT scans of patients affected by complex septal deviations. The computed sensitivity consistently identifies all the anomalies correctly. Virtual surgery, i.e. morphing of the anatomies according to the computed sensitivity, confirms that the characteristics of the nasal airflow improve significantly after small anatomy changes derived from adjoint-based optimization. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Spatial discretization effects in spanwise forcing for turbulent drag reduction

Authors: Emanuele Gallorini, Maurizio Quadrio

Abstract: Wall-based spanwise forcing has been experimentally used with success by Auteri et al. (Phys. Fluids vol. 22, 2010, 115103) to obtain large reductions of turbulent skin-friction drag and considerable energy savings in a pipe flow. The spatial distribution of the azimuthal wall velocity used in the experiment was not continuous, but piecewise-constant. The present study is a numerical replica of th… ▽ More Wall-based spanwise forcing has been experimentally used with success by Auteri et al. (Phys. Fluids vol. 22, 2010, 115103) to obtain large reductions of turbulent skin-friction drag and considerable energy savings in a pipe flow. The spatial distribution of the azimuthal wall velocity used in the experiment was not continuous, but piecewise-constant. The present study is a numerical replica of the experiment, based on a set of direct numerical simulations (DNS); its goal is the identification of the effects of spatially discrete forcing, as opposed to the idealized sinusoidal forcing considered in the majority of numerical studies. Regardless of the discretization, with DNS the maximum drag reduction is found to be larger: the flow easily reaches complete relaminarization, whereas the experiment was capped at 33\% drag reduction. However, the key result stems from the observation that, for the piecewise-constant forcing, the apparent irregularities of the experimental data appear in the simulation data too. They derive from the rich harmonic content of the discontinuous travelling wave, which alters the drag reduction of the sinusoidal forcing. A detailed understanding of the contribution of each harmonic reveals that, whenever e.g. technological limitations constrain one to work far from the optimal forcing parameters, a discrete forcing may perform very differently from the corresponding ideal sinusoid, and in principle can outperform it. However, care should be exercised in comparison, as discrete and continuous forcing have different energy requirements. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comparing flow-based and anatomy-based features in the data-driven study of nasal pathologies

Authors: Andrea Schillaci, Kazuto Hasegawa, Carlotta Pipolo, Giacomo Boracchi, Maurizio Quadrio

Abstract: In several problems involving fluid flows, Computational Fluid Dynamics (CFD) provides detailed quantitative information, and often allows the designer to successfully optimize the system, by minimizing a cost function. Sometimes, however, one cannot improve the system with CFD alone, because a suitable cost function is not readily available: one notable example is diagnosis in medicine. The field… ▽ More In several problems involving fluid flows, Computational Fluid Dynamics (CFD) provides detailed quantitative information, and often allows the designer to successfully optimize the system, by minimizing a cost function. Sometimes, however, one cannot improve the system with CFD alone, because a suitable cost function is not readily available: one notable example is diagnosis in medicine. The field of interest considered here is rhinology: a correct air flow is key for the functioning of the human nose, yet the notion of a functionally normal nose is not available, and a cost function cannot be written. An alternative and attractive pathway to diagnosis and surgery planning is offered by data-driven methods. In this work, we consider the machine-learning study of nasal pathologies caused by anatomic malformations, with the aim of understanding whether fluid dynamic features, available after a CFD analysis, are more effective than purely geometric features in the training of a neural network for regression. Our experiments are carried out on an extremely simplified anatomic model and a correspondingly simple CFD approach; nevertheless, they demonstrate that flow-based features perform better than geometry-based ones, and allow the training of a neural network with fewer inputs, a crucial advantage in fields like medicine. △ Less

Submitted 18 December, 2023; originally announced December 2023.

Comments: Submitted to Flow

Turbulent drag reduction with streamwise travelling waves in the compressible regime

Authors: Federica Gattere, Massimo Zanolini, Davide Gatti, Matteo Bernardini, Maurizio Quadrio

Abstract: The ability of streamwise-travelling waves of spanwise velocity to reduce the turbulent skin friction drag is assessed in the compressible regime. Direct numerical simulations are carried out to compare drag reduction in subsonic, transonic and supersonic channel flows. Compressibility improves the benefits of the travelling waves, in a way that depends on the control parameters: drag reduction be… ▽ More The ability of streamwise-travelling waves of spanwise velocity to reduce the turbulent skin friction drag is assessed in the compressible regime. Direct numerical simulations are carried out to compare drag reduction in subsonic, transonic and supersonic channel flows. Compressibility improves the benefits of the travelling waves, in a way that depends on the control parameters: drag reduction becomes larger than the incompressible one for small frequencies and wavenumbers. However, the improvement depends on the specific procedure employed for comparison. When the Mach number is varied and, at the same time, wall friction is changed by the control, the bulk temperature in the flow can either evolve freely in time until the aerodynamic heating balances the heat flux at the walls, or be constrained such that a fixed percentage of kinetic energy is transformed into thermal energy. Physical arguments suggest that, in the present context, the latter approach should be preferred. Not only it provides a test condition in which the wall-normal temperature profile more realistically mimics that in an external flow, but also leads to a much better scaling of the results, over both the Mach number and the control parameters. Under this comparison, drag reduction is only marginally improved by compressibility. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Submitted to J. Fluid Mech., 22 pages

Linear stability of Poiseuille flow over a steady spanwise Stokes layer

Authors: Daniele Massaro, Fulvio Martinelli, Peter J. Schmid, Maurizio Quadrio

Abstract: The temporal linear stability of plane Poiseuille flow modified by spanwise forcing applied at the walls is considered. The forcing consists of a stationary streamwise distribution of spanwise velocity that generates a steady transversal Stokes layer, known to reduce skin-friction drag in a turbulent flow with little energetic cost. A large numerical study is carried out, where the effects of both… ▽ More The temporal linear stability of plane Poiseuille flow modified by spanwise forcing applied at the walls is considered. The forcing consists of a stationary streamwise distribution of spanwise velocity that generates a steady transversal Stokes layer, known to reduce skin-friction drag in a turbulent flow with little energetic cost. A large numerical study is carried out, where the effects of both the physical and the discretization parameters are thoroughly explored, for three representative subcritical values of the Reynolds number Re. Results show that the spanwise Stokes layer significantly affects the linear stability of the system. For example, at Re=2000 the wall forcing is found to more than double the negative real part of the least-stable eigenvalue, and to decrease by nearly a factor of four the maximum transient growth of perturbation energy. These observations are Re-dependent and further improve at higher $Re$. Comments on the physical implications of the obtained results are provided, suggesting that spanwise forcing might be effective to obtain at the same time a delayed transition to turbulence and a reduced turbulent friction. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Under consideration by PRF

On the comparison between pre- and post-surgery nasal anatomies via computational fluid dynamics

Authors: Eric Segalerba, Gabriele Dini Ciacci, Maurizio Quadrio, Jan O. Pralits

Abstract: Nasal breathing difficulties (NBD) are widespread and difficult to diagnose; the failure rate of their surgical corrections is high. Computational Fluid Dynamics (CFD) enables diagnosis of NBD and surgery planning, by comparing a pre-operative (pre-op) situation with the outcome of virtual surgery (post-op). An equivalent comparison is involved when considering distinct anatomies in the search for… ▽ More Nasal breathing difficulties (NBD) are widespread and difficult to diagnose; the failure rate of their surgical corrections is high. Computational Fluid Dynamics (CFD) enables diagnosis of NBD and surgery planning, by comparing a pre-operative (pre-op) situation with the outcome of virtual surgery (post-op). An equivalent comparison is involved when considering distinct anatomies in the search for the functionally normal nose. Currently, this comparison is carried out in more than one way, under the implicit assumption that results are unchanged, which reflects our limited understanding of the driver of the respiratory function. The study describes how to set up a meaningful comparison. A pre-op anatomy, derived via segmentation from a CT scan, is compared with a post-op anatomy obtained via virtual surgery. State-of-the-art numerical simulations for a steady inspiration carry out the comparison under three types of global constraints, derived from the field of turbulent flow control: a constant pressure drop (CPG) between external ambient and throat, a constant flow rate (CFR) through the airways and a constant power input (CPI) from the lungs can be enforced. A significant difference in the quantities of interest is observed depending on the type of comparison. Global quantities (flow rate, pressure drop, nasal resistance) as well as local ones are affected. The type of flow forcing affects the outcome of the comparison between pre-op and post-op anatomies. Among the three available options, we argue that CPG is the least adequate. Arguments favouring either CFR or CPI are presented. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Curvature effects on the structure of near-wall turbulence

Authors: Davide Selvatici, Maurizio Quadrio, Alessandro Chiarini

Abstract: The interaction between near-wall turbulence and wall curvature is described for the incompressible flow in a plane channel with a small concave-convex-concave bump on the bottom wall, with height comparable to the wall-normal location of the main turbulent structures. The analysis starts from a database generated by a direct numerical simulation and hinges upon the anisotropic generalised Kolmogo… ▽ More The interaction between near-wall turbulence and wall curvature is described for the incompressible flow in a plane channel with a small concave-convex-concave bump on the bottom wall, with height comparable to the wall-normal location of the main turbulent structures. The analysis starts from a database generated by a direct numerical simulation and hinges upon the anisotropic generalised Kolmogorov equations, i.e. the exact budget equations for the second-order structure function tensor. The influence of the bump on the wall cycle and on the energy production, redistribution and transfers is described in the physical and scale spaces. Over the upstream side of the bump, the energy drained from the mean flow to sustain the streamwise fluctuations decreases, and the streaks of high and low streamwise velocity weaken and are stretched spanwise. After the bump tip, instead, the production of streamwise fluctuations grows and the streaks intensify, progressively recovering their characteristic spanwise scale. The wall-normal fluctuations, and thus the quasi-streamwise vortices, are sustained by the mean flow over the upstream side of the bump, while energy flows from the vertical fluctuations to the mean field over the downstream side. On the concave portion of the upstream side, the near-wall fluctuations form structures of spanwise velocity which are consistent with Taylor-Görtler vortices at an early stage of development. Their evolution is described by analysing the scale-space pressure-strain term. A schematic description of the bump flow is presented, in which various regions are identified according to the signs of curvature and streamwise pressure gradient. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: Under consideration by Journal of Fluid Mechanics

On-off pumping for drag reduction in a turbulent channel flow

Authors: Giulio Foggi Rota, Alessandro Monti, Marco E. Rosti, Maurizio Quadrio

Abstract: We show that the energy required by a turbulent flow to displace a given amount of fluid through a straight duct in a given time interval can be reduced by modulating in time the pumping power. The control strategy is hybrid: it is passive, as it requires neither a control system nor control energy, but it manipulates how pumping energy is delivered to the system (as in active techniques) to incre… ▽ More We show that the energy required by a turbulent flow to displace a given amount of fluid through a straight duct in a given time interval can be reduced by modulating in time the pumping power. The control strategy is hybrid: it is passive, as it requires neither a control system nor control energy, but it manipulates how pumping energy is delivered to the system (as in active techniques) to increase the pumping efficiency. Our control employs a temporally periodic pumping pattern, where a short and intense acceleration (in which the pumping system is on) followed by a longer deceleration (in which the pumping system is off) makes the flow alternately visit a quasi-laminar and a turbulent state. The computational study is for a plane channel flow, and employs direct numerical simulations, which present specific computational challenges, for example the highly varying instantaneous value of the Reynolds number, and the importance of discretisation effects. Particular care is devoted to a meaningful definition of drag reduction in the present context. The ability of the forcing to yield significant savings is demonstrated. Since only a small portion of the parameter space is investigated, the best performance of the control technique remains to be assessed. △ Less

Submitted 29 April, 2023; originally announced May 2023.

Comments: Submitted to J. Fluid Mech

Journal ref: Journal of Fluid Mechanics 966, A12 (2023)

Structure function tensor equations with triple decomposition

Authors: Federica Gattere, Alessandro Chiarini, Emanuele Gallorini, Maurizio Quadrio

Abstract: Exact budget equations are derived for the coherent and stochastic contributions to the second-order structure function tensor. They extend the anisotropic generalised Kolmogorov equations (AGKE) by considering the coherent and stochastic parts of the Reynolds stress tensor, and are useful for the statistical description of turbulent flows with periodic or quasi-periodic features, like e.g. the al… ▽ More Exact budget equations are derived for the coherent and stochastic contributions to the second-order structure function tensor. They extend the anisotropic generalised Kolmogorov equations (AGKE) by considering the coherent and stochastic parts of the Reynolds stress tensor, and are useful for the statistical description of turbulent flows with periodic or quasi-periodic features, like e.g. the alternate shedding after a bluff body. While the original AGKE describe production, transport, inter-component redistribution and dissipation of the Reynolds stresses in the combined space of scales and positions, the new equations, called $\varphi$AGKE, contain the phase $\varphi$ as an additional independent variable, and describe the interplay among the mean, coherent and stochastic fields at the various phases. The newly derived $\varphi$AGKE are then applied to a case where an exactly periodic external forcing drives the flow: a turbulent plane channel flow modified by harmonic spanwise oscillations of the wall to reduce drag. The phase-by-phase action of the oscillating transversal Stokes layer generated by the forcing on the near-wall turbulent structures is observed, and a detailed description of the scale-space interaction among mean, coherent and stochastic fields is provided thanks to the $\varphi$AGKE. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Saving Energy in Turbulent Flows with Unsteady Pumping

Authors: Giulio Foggi Rota, Alessandro Monti, Marco E. Rosti, Maurizio Quadrio

Abstract: Viscous dissipation causes significant energy losses in fluid flows; in ducts, laminar flows provide the minimum resistance to the motion, whereas turbulent currents substantially increase the friction at the wall and the energy requirement for pumping. Great effort is currently being devoted to find new ways of reducing the energy losses induced by turbulence. Here we propose a simple and novel d… ▽ More Viscous dissipation causes significant energy losses in fluid flows; in ducts, laminar flows provide the minimum resistance to the motion, whereas turbulent currents substantially increase the friction at the wall and the energy requirement for pumping. Great effort is currently being devoted to find new ways of reducing the energy losses induced by turbulence. Here we propose a simple and novel drag-reduction technique which achieves substantial energy savings in internal flows, even after the energy cost of the control is accounted for. Our approach consists in driving the flow with a pulsed pumping, unlike the common practice of a constant pumping. We alternate "pump on" phases where the flow accelerates, and "pump off" phases where the flow decays freely. The flow cyclically enters a quasi-laminar state during the acceleration, and transitions to a more classic turbulent state during the deceleration. Our numerical results demonstrate that important net energy savings can be achieved by simply modulating the power injection in the system over time. The physical understanding we have achieved can help the industry in reducing the waste of energy, making it cheaper and preserving the environment by reducing harmful emissions. △ Less

Submitted 19 January, 2023; v1 submitted 17 December, 2022; originally announced December 2022.

Comments: 12 pages, 3 figures

Journal ref: Scientific Reports 13, 1299 (2023)

Coherent near-wall structures and drag reduction by spanwise forcing

Authors: Emanuele Gallorini, Maurizio Quadrio, Davide Gatti

Abstract: The effect of streamwise-traveling waves of spanwise wall velocity (StTW) on the quasi-streamwise vortices (QSV) populating the near-wall region of turbulent channels is studied via a conditional averaging technique applied to flow snapshots obtained via Direct Numerical Simulation. The analysis by A. Yakeno, Y. Hasegawa, N. Kasagi, "Modification of quasi-streamwise vortical structure in a drag-re… ▽ More The effect of streamwise-traveling waves of spanwise wall velocity (StTW) on the quasi-streamwise vortices (QSV) populating the near-wall region of turbulent channels is studied via a conditional averaging technique applied to flow snapshots obtained via Direct Numerical Simulation. The analysis by A. Yakeno, Y. Hasegawa, N. Kasagi, "Modification of quasi-streamwise vortical structure in a drag-reduced turbulent channel flow with spanwise wall oscillation", Phys. Fluids 26, 085109 (2014), where the special case of spatially uniform wall oscillation (OW) was considered, is extended to the general case of StTW, which yield both reduction and increase of turbulent skin-friction drag. StTW are found to significantly impact the wall-normal distribution of the vortex population. The conditionally averaged velocity field around the vortices shows that the contributions of the QSV to the quadrant Reynolds shear stresses change significantly during the control cycle. On the one hand, as for OW, the suppression of Q2 events (with upwelling of low-speed fluid away from the wall) dominates the drag-reduction process. On the other, the enhancement of Q2 and also Q4 events (with downwelling of high-speed fluid toward the wall) is related to drag increase. Based on the link identified between the phase changes of the Reynolds stresses and the principal directions of the rate-of-strain tensor induced by the StTW, a predictive correlation for drag reduction by StTW is proposed which uses physically significant parameters to overcome the shortcomings of existing models. △ Less

Submitted 21 August, 2022; originally announced August 2022.

Comments: 17 pages, 13 figures

Structure of turbulence in the flow around a rectangular cylinder

Authors: Alessandro Chiarini, Davide Gatti, Andrea Cimarelli, Maurizio Quadrio

Abstract: The separating and reattaching turbulent flow past a rectangular cylinder is studied to describe how small and large scales contribute to the sustaining mechanism of the velocity fluctuations. The work is based on the Anisotropic Generalised Kolmogorov Equations (AGKE), exact budget equations for the second-order structure function tensor in the space of scales and in the physical space. Scale-spa… ▽ More The separating and reattaching turbulent flow past a rectangular cylinder is studied to describe how small and large scales contribute to the sustaining mechanism of the velocity fluctuations. The work is based on the Anisotropic Generalised Kolmogorov Equations (AGKE), exact budget equations for the second-order structure function tensor in the space of scales and in the physical space. Scale-space energy fluxes show that forward and reverse energy transfers simultaneously occur in the flow, with interesting modelling implications. Over the longitudinal cylinder side, the Kelvin-Helmholtz instability of the leading-edge shear layer generates large spanwise rolls, which get stretched into hairpin-like vortices and eventually break down into smaller streamwise vortices. Independent sources of velocity fluctuations act at different scales. The flow dynamics is dominated by pressure-strain: the flow impingement on the cylinder surface in the reattachment zone produces spanwise velocity fluctuations very close to the wall and, at larger wall distances, reorients them to feed streamwise-aligned vortices. In the near wake, large von Kàrmàn-like vortices are shed from the trailing edge and coexist with smaller turbulent structures, each with its own independent production mechanism. At the trailing edge, the sudden disappearance of the wall changes the structure of turbulence: streamwise vortices progressively vanish, while spanwise structures close to the wall are suddenly turned into vertical fluctuations by the pressure strain. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Dimples for skin-friction drag reduction: status and perspectives

Authors: Federica Gattere, Alessandro Chiarini, Maurizio Quadrio

Abstract: Dimples are small concavities imprinted on a flat surface, known to affect heat transfer and also flow separation and aerodynamic drag on bluff bodies when acting as a standard roughness. Recently, dimples have been proposed as a roughness pattern that is capable to reduce the turbulent drag of a flat plate, by providing a reduction of skin friction that compensates the dimple-induced pressure dra… ▽ More Dimples are small concavities imprinted on a flat surface, known to affect heat transfer and also flow separation and aerodynamic drag on bluff bodies when acting as a standard roughness. Recently, dimples have been proposed as a roughness pattern that is capable to reduce the turbulent drag of a flat plate, by providing a reduction of skin friction that compensates the dimple-induced pressure drag, and leads to a global benefit. The question whether dimples do actually work to reduce friction drag is still unsettled. In this paper, we provide a comprehensive review of the available information, touching upon the many parameters that characterize the problem. A number of reasons that contribute to explaining the contrasting literature information are discussed. We also provide guidelines for future studies, by highlighting key methodological steps required for a meaningful comparison between a flat and dimpled surface in view of drag reduction. △ Less

Submitted 10 July, 2022; originally announced July 2022.

Drag reduction on a transonic airfoil

Authors: Maurizio Quadrio, Alessandro Chiarini, Jacopo Banchetti, Davide Gatti, Antonio Memmolo, Sergio Pirozzoli

Abstract: Flow control for turbulent skin-friction drag reduction is applied to a transonic airfoil to improve its aerodynamic performance. The study is based on direct numerical simulations (with up to 1.8 billions cells) of the compressible turbulent flow around a supercritical airfoil, at Reynolds and Mach numbers of $Re_\infty= 3 \times 10^5$ and $M_\infty =0.7$. Control via spanwise forcing is applied… ▽ More Flow control for turbulent skin-friction drag reduction is applied to a transonic airfoil to improve its aerodynamic performance. The study is based on direct numerical simulations (with up to 1.8 billions cells) of the compressible turbulent flow around a supercritical airfoil, at Reynolds and Mach numbers of $Re_\infty= 3 \times 10^5$ and $M_\infty =0.7$. Control via spanwise forcing is applied over a fraction of the suction side of the airfoil. Besides locally reducing friction, the control modifies the shock wave and significantly improves the aerodynamic efficiency of the airfoil by increasing lift and decreasing drag. Hence, the airfoil can achieve the required lift at a lower angle of attack and with a lower drag. Estimates at the aircraft level indicate that substantial savings are possible; when control is active, its energy cost becomes negligible thanks to the small application area. We suggest that skin-friction drag reduction should be considered not only as a goal, but also as a tool to improve the global aerodynamics of complex flows. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Ascending-descending and direct-inverse cascades of Reynolds stresses in turbulent Couette flow

Authors: Alessandro Chiarini, Mariadebora Mauriello, Davide Gatti, Maurizio Quadrio

Abstract: The interaction between small- and large-scale structures, and the coexisting bottom-up and top-down processes are studied in a turbulent plane Couette flow, where space-filling longitudinal rolls appear at relatively low values of the Reynolds number $Re$. A DNS database at $Re_τ=101$ is built to replicate the highest $Re$ considered in a recent experimental work by Kawata and Alfredsson (Phys. R… ▽ More The interaction between small- and large-scale structures, and the coexisting bottom-up and top-down processes are studied in a turbulent plane Couette flow, where space-filling longitudinal rolls appear at relatively low values of the Reynolds number $Re$. A DNS database at $Re_τ=101$ is built to replicate the highest $Re$ considered in a recent experimental work by Kawata and Alfredsson (Phys. Rev. Lett., vol.120, 2018, 244501). Our study is based on the exact budget equations for the second-order structure function tensor $\langle δu_i δu_j \rangle$, i.e. the Anisotropic Generalized Kolmogorov Equations (AGKE). The AGKE study production, redistribution, transport and dissipation of every Reynolds stress tensor component, considering simultaneously the physical space and the space of scales, and properly define the concept of scale in the inhomogeneous wall-normal direction. We show how the large-scale energy-containing motions are involved in the production and redistribution of the turbulent fluctuations. Both bottom-up and top-down interactions occur, and the same is true for direct and inverse cascading. The wall-parallel components $\langle δu δu \rangle$ and $\langle δw δw \rangle$ show that the both small and large near-wall scales feed the large scales away from the wall. The wall-normal component $\langle δv δv \rangle$ is different, and shows a dominant top-down dynamics, being produced via pressure-strain redistribution away from the wall and transferred towards near-wall larger scales via an inverse cascade. The off-diagonal component shows a top-down interaction, with both direct and inverse cascade, albeit the latter takes place within a limited range of scales. △ Less

Submitted 23 September, 2021; originally announced September 2021.

The importance of corner sharpness in the BARC test case: a numerical study

Authors: Alessandro Chiarini, Maurizio Quadrio

Abstract: The BARC flow is studied via Direct Numerical Simulation at a relatively low turbulent Reynolds number, with focus on the geometrical representation of the leading-edge (LE) corners. The study contributes to further our understanding of the discrepancies between existing numerical and experimental BARC data. In a first part, rounded LE corners with small curvature radii are considered. Results sho… ▽ More The BARC flow is studied via Direct Numerical Simulation at a relatively low turbulent Reynolds number, with focus on the geometrical representation of the leading-edge (LE) corners. The study contributes to further our understanding of the discrepancies between existing numerical and experimental BARC data. In a first part, rounded LE corners with small curvature radii are considered. Results show that a small amount of rounding does not lead to abrupt changes of the mean fields, but that the effects increase with the curvature radius. The shear layer separates from the rounded LE at a lower angle, which reduces the size of the main recirculating region over the cylinder side. In contrast, the longitudinal size of the recirculating region behind the trailing edge (TE) increases, as the TE shear layer is accelerated. The effect of the curvature radii on the turbulent kinetic energy and on its production, dissipation and transport are addressed. The present results should be contrasted with the recent work of Rocchio et al, JWEIA 2020, who found via implicit Large-Eddy Simulations at larger Reynolds numbers than even a small curvature radius leads to significant changes of the mean flow. In a second part, the LE corners are fully sharp and the exact analytical solution of the Stokes problem in the neighborhood of the corner is used to locally restore the solution accuracy degraded by the singularity. Changes in the mean flow reveal that the analytical correction leads to streamlines that better follow the corners. The flow separates from the LE with a lower angle, resulting in a slightly smaller recirculating region. The corner-correction approach is valuable in general, and is expected to help developing high-quality numerical simulations at the high Reynolds numbers typical of the experiments with reasonable meshing requirements. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Global energy budgets in turbulent Couette and Poiseuille flows

Authors: Andrea Andreolli, Maurizio Quadrio, Davide Gatti

Abstract: Turbulent plane Poiseuille and Couette flows share the same geometry, but produce their flow rate owing to different external drivers, pressure gradient and shear respectively. By looking at integral energy fluxes, we pose and answer the question of which flow performs better at creating flow rate. We define a flow {\em efficiency}, that quantifies the fraction of power used to produce flow rate i… ▽ More Turbulent plane Poiseuille and Couette flows share the same geometry, but produce their flow rate owing to different external drivers, pressure gradient and shear respectively. By looking at integral energy fluxes, we pose and answer the question of which flow performs better at creating flow rate. We define a flow {\em efficiency}, that quantifies the fraction of power used to produce flow rate instead of being wasted as a turbulent overhead; {\em effectiveness}, instead, describes the amount of flow rate produced by a given power. The work by Gatti \emph{et al.} (\emph{J. Fluid Mech.} vol.857, 2018, pp. 345--373), where the constant power input (CPI) concept was developed to compare turbulent Poiseuille flows with drag reduction, is here extended to compare different flows. By decomposing the mean velocity field into a laminar contribution and a deviation, analytical expressions are derived which are the energy-flux equivalents of the FIK identity. These concepts are applied to literature data supplemented by a new set of direct numerical simulations, to find that Couette flows are less efficient but more effective than Poiseuille ones. The reason is traced to the more effective laminar component of Couette flows, which compensates for their higher turbulent activity. It is also observed that, when the fluctuating fields of the two flows are fed with the same total power fraction, Couette flows dissipate a smaller percentage of it via turbulent dissipation. A decomposition of the fluctuating field into large and small scales explains this feature: Couette flows develop stronger large-scale structures, which alter the mean flow while contributing less significantly to dissipation. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Journal ref: J. Fluid Mech. 924 (2021) A25

Structure function tensor equations in inhomogeneous turbulence

Authors: Davide Gatti, Alessandro Chiarini, Andrea Cimarelli, Maurizio Quadrio

Abstract: Exact budget equations for the second-order structure function tensor $\langle δu_i δu_j \rangle$ are used to study the two-point statistics of velocity fluctuations in inhomogeneous turbulence. The Anisotropic Generalized Kolmogorov Equations (AGKE) describe the production, transport, redistribution and dissipation of every Reynolds stress component occurring simultaneously among different scales… ▽ More Exact budget equations for the second-order structure function tensor $\langle δu_i δu_j \rangle$ are used to study the two-point statistics of velocity fluctuations in inhomogeneous turbulence. The Anisotropic Generalized Kolmogorov Equations (AGKE) describe the production, transport, redistribution and dissipation of every Reynolds stress component occurring simultaneously among different scales and in space, i.e. along directions of statistical inhomogeneity. The AGKE are effective to study the inter-component and multi-scale processes of turbulence. In contrast to more classic approaches, such as those based on the spectral decomposition of the velocity field, the AGKE provide a natural definition of scales in the inhomogeneous directions, and describe fluxes across such scales too. Compared to the Generalized Kolmogorov Equation, which is recovered as their half trace, the AGKE can describe inter-component energy transfers occurring via the pressure-strain term and contain also budget equations for the off-diagonal components of $\langle δu_i δu_j \rangle$. The non-trivial physical interpretation of the AGKE terms is demonstrated with three examples. First, the near-wall cycle of a turbulent channel flow at $Re_τ=200$ is considered. The off-diagonal component $\langle -δu δv \rangle$, which can not be interpreted in terms of scale energy, is discussed in detail. Wall-normal scales in the outer turbulence cycle are then discussed by applying the AGKE to channel flows at $Re_τ=500$ and $1000$. In a third example, the AGKE are computed for a separating and reattaching flow. The process of spanwise-vortex formation in the reverse boundary layer within the separation bubble is discussed for the first time. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Journal ref: J. Fluid Mech. 898 (2020) A5

Turbulent drag reduction over curved walls

Authors: Jacopo Banchetti, Paolo Luchini, Maurizio Quadrio

Abstract: This work studies the effects of skin-friction drag reduction in a turbulent flow over a curved wall, with a view to understanding the relationship between the reduction of friction and changes to the total aerodynamic drag. Direct numerical simulations (DNS) are carried out for an incompressible turbulent flow in a channel where one wall has a small bump; two bump geometries are considered, that… ▽ More This work studies the effects of skin-friction drag reduction in a turbulent flow over a curved wall, with a view to understanding the relationship between the reduction of friction and changes to the total aerodynamic drag. Direct numerical simulations (DNS) are carried out for an incompressible turbulent flow in a channel where one wall has a small bump; two bump geometries are considered, that produce mildly separated and attached flows. Friction drag reduction is achieved by applying streamwise-travelling waves of spanwise velocity (StTW). The local friction reduction produced by the StTW is found to vary along the curved wall, leading to a global friction reduction that, for the cases studied, is up to 10\% larger than that obtained in the plane-wall case. Moreover, the modified skin friction induces non-negligible changes of pressure drag, which is favorably affected by StTW and globally reduces by up to 10\%. The net power saving, accounting for the power required to create the StTW, is positive and, for the cases studied, is one half larger than the net saving of the planar case. The study suggests that reducing friction at the surface of a body of complex shape induces further effects, a simplistic evaluation of which might lead to underestimating the total drag reduction. △ Less

Submitted 25 April, 2020; originally announced April 2020.

Journal ref: J. Fluid Mech. 896 (2020) A10

An Efficient Numerical Method for the Generalised Kolmogorov Equation

Authors: Davide Gatti, Alberto Remigi, Alessandro Chiarini, Andrea Cimarelli, Maurizio Quadrio

Abstract: An efficient algorithm for computing the terms appearing in the Generalised Kolmogorov Equation (GKE) written for the indefinite plane channel flow is presented. The algorithm, which features three distinct strategies for parallel computing, is designed such that CPU and memory requirements are kept to a minimum, so that high-Re wall-bounded flows can be afforded. Computational efficiency is mainl… ▽ More An efficient algorithm for computing the terms appearing in the Generalised Kolmogorov Equation (GKE) written for the indefinite plane channel flow is presented. The algorithm, which features three distinct strategies for parallel computing, is designed such that CPU and memory requirements are kept to a minimum, so that high-Re wall-bounded flows can be afforded. Computational efficiency is mainly achieved by leveraging the Parseval's theorem for the two homogeneous directions available in the plane channel geometry. A speedup of 3-4 orders of magnitude, depending on the problem size, is reported in comparison to a key implementation used in the literature. Validation of the code is demonstrated by computing the residual of the GKE, and example results are presented for channel flows at Re_τ=200 and Re_τ=1000, where for the first time they are observed in the whole four-dimensional domain. It is shown that the space and scale properties of the scale-energy fluxes change for increasing values of the Reynolds number. Among all scale-energy fluxes, the wall-normal flux is found to show the richest behaviour for increasing streamwise scales. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Journal ref: Journal of Turbulence 20(8); 457--480; 2019

Global energy fluxes in fully-developed turbulent channels with flow control

Authors: Davide Gatti, Andrea Cimarelli, Yosuke Hasegawa, Bettina Frohnapfel, Maurizio Quadrio

Abstract: This paper addresses the integral energy fluxes in natural and controlled turbulent channel flows, where active skin-friction drag reduction techniques allow a more efficient use of the available power. We study whether the increased efficiency shows any general trend in how energy is dissipated by the mean velocity field (mean dissipation) and by the fluctuating velocity field (turbulent dissipat… ▽ More This paper addresses the integral energy fluxes in natural and controlled turbulent channel flows, where active skin-friction drag reduction techniques allow a more efficient use of the available power. We study whether the increased efficiency shows any general trend in how energy is dissipated by the mean velocity field (mean dissipation) and by the fluctuating velocity field (turbulent dissipation). Direct Numerical Simulations (DNS) of different control strategies are performed at Constant Power Input (CPI), so that at statistical equilibrium each flow (either uncontrolled or controlled by different means) has the same power input, hence the same global energy flux and, by definition, the same total energy dissipation rate. The simulations reveal that changes in mean and turbulent energy dissipation rates can be of either sign in a successfully controlled flow. A quantitative description of these changes is made possible by a new decomposition of the total dissipation, stemming from an extended Reynolds decomposition, where the mean velocity is split into a laminar component and a deviation from it. Thanks to the analytical expressions of the laminar quantities, exact relationships are derived that link the achieved flow rate increase and all energy fluxes in the flow system with two wall-normal integrals of the Reynolds shear stress and the Reynolds number. The dependence of the energy fluxes on the Reynolds number is elucidated with a simple model in which the control-dependent changes of the Reynolds shear stress are accounted for via a modification of the mean velocity profile. The physical meaning of the energy fluxes stemming from the new decomposition unveils their inter-relations and connection to flow control, so that a clear target for flow control can be identified. △ Less

Submitted 3 September, 2018; originally announced September 2018.

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

Authors: Davide Gatti, Maurizio Quadrio

Abstract: 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 create… ▽ More 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_τ=200$ and $Re_τ=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. △ Less

Submitted 13 July, 2016; v1 submitted 26 December, 2015; originally announced December 2015.

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

Authors: Maurizio Quadrio, Bettina Frohnapfel, Yosuke Hasegawa

Abstract: 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… ▽ More 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. △ Less

Submitted 16 September, 2015; originally announced September 2015.

Modal and nonmodal stability analysis of electrohydrodynamic flow with and without cross-flow

Authors: Mengqi Zhang, Fulvio Martinelli, Jian Wu, Peter J. Schmid, Maurizio Quadrio

Abstract: We report the results of a complete modal and nonmodal linear stability analysis of the electrohydrodynamic flow (EHD) for the problem of electroconvection in the strong injection region. Convective cells are formed by Coulomb force in an insulating liquid residing between two plane electrodes subject to unipolar injection. Besides pure electroconvection, we also consider the case where a cross-fl… ▽ More We report the results of a complete modal and nonmodal linear stability analysis of the electrohydrodynamic flow (EHD) for the problem of electroconvection in the strong injection region. Convective cells are formed by Coulomb force in an insulating liquid residing between two plane electrodes subject to unipolar injection. Besides pure electroconvection, we also consider the case where a cross-flow is present, generated by a streamwise pressure gradient, in the form of a laminar Poiseuille flow. The effect of charge diffusion, often neglected in previous linear stability analyses, is included in the present study and a transient growth analysis, rarely considered in EHD, is carried out. In the case without cross-flow, a non-zero charge diffusion leads to a lower linear stability threshold and thus to a more unstable low. The transient growth, though enhanced by increasing charge diffusion, remains small and hence cannot fully account for the discrepancy of the linear stability threshold between theoretical and experimental results. When a cross-flow is present, increasing the strength of the electric field in the high-$Re$ Poiseuille flow yields a more unstable flow in both modal and nonmodal stability analyses. Even though the energy analysis and the input-output analysis both indicate that the energy growth directly related to the electric field is small, the electric effect enhances the lift-up mechanism. The symmetry of channel flow with respect to the centerline is broken due to the additional electric field acting in the wall-normal direction. As a result, the centers of the streamwise rolls are shifted towards the injector electrode, and the optimal spanwise wavenumber achieving maximum transient energy growth increases with the strength of the electric field. △ Less

Submitted 2 April, 2015; originally announced April 2015.

Direct numerical simulation of turbulent channel flow over porous walls

Authors: Marco E. Rosti, Luca Cortelezzi, Maurizio Quadrio

Abstract: 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 dimens… ▽ More 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_τ=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. △ Less

Submitted 23 September, 2015; v1 submitted 27 October, 2014; originally announced October 2014.

Comments: Revised version, with additional data and more in-depth analysis

Numerical simulation of turbulent duct flows with constant power input

Authors: Yosuke Hasegawa, Maurizio Quadrio, Bettina Frohnapfel

Abstract: 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… ▽ More 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. △ Less

Submitted 7 May, 2014; originally announced May 2014.

Comments: Accepted for publication in J. Fluid Mech

Prediction of turbulence control for arbitrary periodic spanwise wall movement

Authors: Andrea Cimarelli, Bettina Frohnapfel, Yosuke Hasegawa, Elisabetta De Angelis, Maurizio Quadrio

Abstract: In order to generalize the well-known spanwise-oscillating-wall technique for drag reduction, non-sinusoidal oscillations of a solid wall are considered as a means to alter the skin-friction drag in a turbulent channel flow. A series of Direct Numerical Simulations is conducted to evaluate the control performance of nine different temporal waveforms, in addition to the usual sinusoid, systematical… ▽ More In order to generalize the well-known spanwise-oscillating-wall technique for drag reduction, non-sinusoidal oscillations of a solid wall are considered as a means to alter the skin-friction drag in a turbulent channel flow. A series of Direct Numerical Simulations is conducted to evaluate the control performance of nine different temporal waveforms, in addition to the usual sinusoid, systematically changing the wave amplitude and the period for each waveform. The turbulent average spanwise motion is found to coincide with the laminar Stokes solution that is constructed, for the generic waveform, through harmonic superposition. This allows us to define and compute, for each waveform, a new penetration depth of the Stokes layer which correlates with the amount of turbulent drag reduction, and eventually to predict both turbulent drag reduction and net energy saving rate for arbitrary waveforms. Among the waveforms considered, the maximum net energy saving rate is obtained by the sinusoidal wave at its optimal amplitude and period. However, the sinusoid is not the best waveform at every point in the parameter space. Our predictive tool offers simple guidelines to design waveforms that outperform the sinusoid for given (suboptimal) amplitude and period of oscillation. This is potentially interesting in view of applications, where physical limitations often preclude the actuator to reach its optimal operating conditions. △ Less

Submitted 19 June, 2013; originally announced June 2013.

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

Authors: Marco Belan, Maurizio Quadrio

Abstract: 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 p… ▽ More 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. △ Less

Submitted 22 January, 2013; originally announced January 2013.

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

Authors: Davide Gatti, Maurizio Quadrio

Abstract: 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… ▽ More 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_τ=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. △ Less

Submitted 30 December, 2013; v1 submitted 25 December, 2012; originally announced December 2012.

Comments: Revised version: change of title, revised intro, small improvements to figures and text

Money versus Time: Evaluation of Flow Control in Terms of Energy Consumption and Convenience

Authors: Bettina Frohnapfel, Yosuke Hasegawa, Maurizio Quadrio

Abstract: Flow control with the goal of reducing the skin friction drag on the fluid-solid interface is an active fundamental research area, motivated by its potential for significant energy savings and reduced emissions in the transport sector. Customarily, the performance of drag reduction techniques in internal flows is evaluated under two alternative flow conditions, i.e. at constant mass flow rate or c… ▽ More Flow control with the goal of reducing the skin friction drag on the fluid-solid interface is an active fundamental research area, motivated by its potential for significant energy savings and reduced emissions in the transport sector. Customarily, the performance of drag reduction techniques in internal flows is evaluated under two alternative flow conditions, i.e. at constant mass flow rate or constant pressure gradient. Successful control leads to reduction of drag and pumping power within the former approach, whereas the latter leads to an increase of the mass flow rate and pumping power. In practical applications, however, money and time define the flow control challenge: a compromise between the energy expenditure (money) and the corresponding convenience (flow rate) achieved with that amount of energy has to be reached so as to accomplish a goal which in general depends on the specific application. Based on this idea, we derive two dimensionless parameters which quantify the total energy consumption and the required time (convenience) for transporting a given volume of fluid through a given duct. Performances of existing drag reduction strategies as well as the influence of wall roughness are re-evaluated within the present framework; how to achieve the (application-dependent) optimum balance between energy consumption and convenience is addressed. It is also shown that these considerations can be extended to external flows. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Changes in turbulent dissipation in a channel flow with oscillating walls

Authors: Pierre Ricco, Claudio Ottonelli, Yosuke Hasegawa, Maurizio Quadrio

Abstract: Harmonic oscillations of the walls of a turbulent plane channel flow are studied by direct numerical simulations to improve our understanding of the physical mechanism for skin-friction drag reduction. The simulations are carried out at constant pressure gradient in order to define an unambiguous inner scaling: in this case, drag reduction manifests itself as an increase of mass flow rate. Energy… ▽ More Harmonic oscillations of the walls of a turbulent plane channel flow are studied by direct numerical simulations to improve our understanding of the physical mechanism for skin-friction drag reduction. The simulations are carried out at constant pressure gradient in order to define an unambiguous inner scaling: in this case, drag reduction manifests itself as an increase of mass flow rate. Energy and enstrophy balances, carried out to emphasize the role of the oscillating spanwise shear layer, show that the viscous dissipations of the mean flow and of the turbulent fluctuations increase with the mass flow rate, and the relative importance of the latter decreases. We then focus on the turbulent enstrophy: through an analysis of the temporal evolution from the beginning of the wall motion, the dominant, oscillation-related term in the turbulent enstrophy is shown to cause the turbulent dissipation to be enhanced in absolute terms, before the slow drift towards the new quasi-equilibrium condition. This mechanism is found to be responsible for the increase in mass flow rate. We finally show that the time-average volume integral of the dominant term relates linearly to the drag reduction. △ Less

Submitted 16 February, 2012; originally announced February 2012.

Linear feedback control of transient energy growth and control performance limitations in subcritical plane Poiseuille flow

Authors: F. Martinelli, M. Quadrio, J. McKernan, J. F. Whidborne

Abstract: Suppression of the transient energy growth in subcritical plane Poiseuille flow via feedback control is addressed. It is assumed that the time derivative of any of the velocity components can be imposed at the walls as control input, and that full-state information is available. We show that it is impossible to design a linear state-feedback controller that leads to a closed-loop flow system witho… ▽ More Suppression of the transient energy growth in subcritical plane Poiseuille flow via feedback control is addressed. It is assumed that the time derivative of any of the velocity components can be imposed at the walls as control input, and that full-state information is available. We show that it is impossible to design a linear state-feedback controller that leads to a closed-loop flow system without transient energy growth. In a subsequent step, full-state feedback controllers -- directly targeting the transient growth mechanism -- are designed, using a procedure based on a Linear Matrix Inequalities approach. The performance of such controllers is analyzed first in the linear case, where comparison to previously proposed linear-quadratic optimal controllers is made; further, transition thresholds are evaluated via Direct Numerical Simulations of the controlled three-dimensional Poiseuille flow against different initial conditions of physical interest, employing different velocity components as wall actuation. The present controllers are effective in increasing the transition thresholds in closed loop, with varying degree of performance depending on the initial condition and the actuation component employed. △ Less

Submitted 13 January, 2011; originally announced January 2011.

DNS-based measurement of the mean impulse response of homogeneous isotropic turbulence

Authors: Marco Carini, Maurizio Quadrio

Abstract: A technique for measuring the mean impulse response function of stationary homogeneous isotropic turbulence is proposed. Such measurement is carried out here on the basis of Direct Numerical Simulation (DNS). A zero-mean white-noise volume forcing is used to probe the turbulent flow, and the response function is obtained by accumulating the space-time correlation between the white forcing and the… ▽ More A technique for measuring the mean impulse response function of stationary homogeneous isotropic turbulence is proposed. Such measurement is carried out here on the basis of Direct Numerical Simulation (DNS). A zero-mean white-noise volume forcing is used to probe the turbulent flow, and the response function is obtained by accumulating the space-time correlation between the white forcing and the velocity field. This technique to measure the turbulent response in a DNS numerical experiment is a new research tool in that field of spectral closures where the linear response concept is invoked either by resorting to renormalized perturbations theories or by introducing the well-known Fluctuation-Dissipation Relation (FDR). Though the results obtained in the present work are limited to relatively low values of the Reynolds number, a preliminary analysis is possible. Both the characteristic form and the time scaling properties of the response function are investigated in the universal subrange of dissipative wavenumbers; a comparison with the response approximation given by the FDR is proposed through the independent DNS measurement of the correlation function. Very good agreement is found between the measured response and Kraichnan's description of random energy-range advection effects. △ Less

Submitted 10 November, 2010; originally announced November 2010.

The laminar generalized Stokes layer and turbulent drag reduction

Authors: Maurizio Quadrio, Pierre Ricco

Abstract: 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 e… ▽ More 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. △ Less

Submitted 24 August, 2010; originally announced August 2010.

Comments: Accepted for publication on J. Fluid Mech

Experimental assessment of drag reduction by traveling waves in a turbulent pipe flow

Authors: Franco Auteri, Arturo Baron, Marco Belan, Gabriele Campanardi, Maurizio Quadrio

Abstract: We experimentally assess the capabilities of an active, open-loop technique for drag reduction in turbulent wall flows recently introduced by Quadrio et al. [J. Fluid Mech., v.627, 161, (2009)]. The technique consists in generating streamwise-modulated waves of spanwise velocity at the wall, that travel in the streamwise direction. A proof-of-principle experiment has been devised to measure the… ▽ More We experimentally assess the capabilities of an active, open-loop technique for drag reduction in turbulent wall flows recently introduced by Quadrio et al. [J. Fluid Mech., v.627, 161, (2009)]. The technique consists in generating streamwise-modulated waves of spanwise velocity at the wall, that travel in the streamwise direction. A proof-of-principle experiment has been devised to measure the reduction of turbulent friction in a pipe flow, in which the wall is subdivided into thin slabs that rotate independently in the azimuthal direction. Different speeds of nearby slabs provide, although in a discrete setting, the desired streamwise variation of transverse velocity. Our experiment confirms the available DNS results, and in particular demonstrates the possibility of achieving large reductions of friction in the turbulent regime. Reductions up to 33% are obtained for slowly forward-traveling waves; backward-traveling waves invariably yield drag reduction, whereas a substantial drop of drag reduction occurs for waves traveling forward with a phase speed comparable to the convection speed of near-wall turbulent structures. A Fourier analysis is employed to show that the first harmonics introduced by the discrete spatial waveform that approximates the sinusoidal wave are responsible for significant effects that are indeed observed in the experimental measurements. Practical issues related to the physical implementation of this control scheme and its energetic efficiency are briefly discussed. △ Less

Submitted 16 August, 2010; originally announced August 2010.

Comments: Article accepted by Phys. Fluids. After it is published, it will be found at http://pof.aip.org

Turbulent superfluid profiles in a counterflow channel

Authors: L. Galantucci, C. F. Barenghi, M. Sciacca, M. Quadrio, P. Luchini

Abstract: We have developed a two-dimensional model of quantised vortices in helium II moving under the influence of applied normal fluid and superfluid in a counterflow channel. We predict superfluid and vortex-line density profiles which could be experimentally tested using recently developed visualization techniques. We have developed a two-dimensional model of quantised vortices in helium II moving under the influence of applied normal fluid and superfluid in a counterflow channel. We predict superfluid and vortex-line density profiles which could be experimentally tested using recently developed visualization techniques. △ Less

Submitted 16 June, 2010; originally announced June 2010.

Comments: 3 double figures, 9 pages

Very fine near-wall structures in scalar mixing

Authors: Luca Galantucci, Maurizio Quadrio

Abstract: Passive scalar dynamics in wall-bounded turbulence is studied via Direct Numerical Simulations of plane channel flow, for a friction Reynolds number $Re_* = 160$ and a Schmidt number $Sc=1$. Peculiar to the present research is that the spatial resolution reaches far beyond what has been employed in similar past studies. Our aim is to examine the statistics of the most dissipative events across the… ▽ More Passive scalar dynamics in wall-bounded turbulence is studied via Direct Numerical Simulations of plane channel flow, for a friction Reynolds number $Re_* = 160$ and a Schmidt number $Sc=1$. Peculiar to the present research is that the spatial resolution reaches far beyond what has been employed in similar past studies. Our aim is to examine the statistics of the most dissipative events across the various layers of the channel flow, and to compare them to the homogeneous isotropic case, where the recent studies by Schumacher et al (2005) and Watanabe and Gotoh (2007) have described a range of scalar micro-scales that require extremely high spatial resolution to be properly resolved. Resolution effects are observed on integral-scale quantities such as the mean profiles of the scalar dissipation and its variance. By examining probability distributions, it is found that the finest resolution is essential for correctly computing small-scale statistics in the near-wall region of the channel. As expected, this high-resolution requirement extends outwards to the channel centerline, where the behavior of isotropic turbulence is recovered. However, high-intensity scalar dissipation events are overemphasized by marginal resolution near the wall, while they are underemphasized in the central region. △ Less

Submitted 13 January, 2011; v1 submitted 20 May, 2010; originally announced May 2010.

Comments: Revised version: a word was missing from the title

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

Authors: Claudio Viotti, Maurizio Quadrio, Paolo Luchini

Abstract: 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 f… ▽ More 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. △ Less

Submitted 31 October, 2009; originally announced November 2009.

Streamwise-traveling waves of spanwise wall velocity for turbulent drag reduction

Authors: M. Quadrio, P. Ricco, C. Viotti

Abstract: Waves of spanwise velocity imposed at the walls of a plane turbulent channel flow are studied by Direct Numerical Simulations. We consider sinusoidal waves of spanwise velocity which vary in time and are modulated in space along the streamwise direction. The phase speed may be null, positive or negative, so that the waves may be either stationary or traveling forward or backward in the direction… ▽ More Waves of spanwise velocity imposed at the walls of a plane turbulent channel flow are studied by Direct Numerical Simulations. We consider sinusoidal waves of spanwise velocity which vary in time and are modulated in space along the streamwise direction. The phase speed may be null, positive or negative, so that the waves may be either stationary or traveling forward or backward in the direction of the mean flow. Such a forcing includes as particular cases two known techniques for reducing friction drag: the oscillating wall technique (a traveling wave with infinite phase speed) and the recently proposed steady distribution of spanwise velocity (a wave with zero phase speed). The traveling waves alter the friction drag significantly. Waves which slowly travel forward produce a large reduction of drag, that can relaminarize the flow at low values of the Reynolds number. Faster waves yield a totally different outcome, i.e. drag increase. Even faster waves produce a drag reduction effect again. Backward-traveling waves instead lead to drag reduction at any speed. The traveling waves, when they reduce drag, operate in similar fashion to the oscillating wall, with an improved energetic efficiency. Drag increase is observed when the waves travel at a speed comparable with that of the convecting near-wall turbulence structures. A diagram illustrating the different flow behaviors is presented. △ Less

Submitted 8 January, 2009; originally announced January 2009.

Wall-oscillation conditions for drag reduction in turbulent channel flow

Authors: Pierre Ricco, Maurizio Quadrio

Abstract: The drag reduction properties of a turbulent channel flow modified by spanwise sinusoidal oscillations of the walls are investigated by direct numerical simulations. The work is based on the linear relation between the drag reduction and the parameter $S$, function of the maximum wall velocity and the period of the oscillation. This quantity has been found through physical arguments pertaining t… ▽ More The drag reduction properties of a turbulent channel flow modified by spanwise sinusoidal oscillations of the walls are investigated by direct numerical simulations. The work is based on the linear relation between the drag reduction and the parameter $S$, function of the maximum wall velocity and the period of the oscillation. This quantity has been found through physical arguments pertaining to the action of the oscillating Stokes layer on the near-wall turbulence dynamics. The predictive potential of the scaling parameter is exploited to gain insight into the drag-reducing effects of the oscillating wall technique. The period of oscillation which guarantees the maximum drag reduction for a given maximum wall displacement is studied for the first time. The issue of the minimum intensity of wall forcing required to produce a non-zero drag reduction effect and the dependence of the drag reduction on the Reynolds number are also addressed. The drag reduction data available in the literature are compared with the prediction given by the scaling parameter, thus attaining a comprehensive view of the state of the art. △ Less

Submitted 1 February, 2008; originally announced February 2008.

Comments: Accepted for publication in the International Journal of Heat and Fluid Flow

Direct Numerical Simulation of turbulent Taylor-Couette flow

Authors: Davide Pirro, Maurizio Quadrio

Abstract: The direct numerical simulation (DNS) of the Taylor--Couette flow in the fully turbulent regime is described. The numerical method extends the work by Quadrio & Luchini (Eur. J. Mech. B / Fluids, v.21, pp.413--427, 2002), and is based on a parallel computer code which uses mixed spatial discretization (spectral schemes in the homogeneous directions, and fourth-order, compact explicit finite-diff… ▽ More The direct numerical simulation (DNS) of the Taylor--Couette flow in the fully turbulent regime is described. The numerical method extends the work by Quadrio & Luchini (Eur. J. Mech. B / Fluids, v.21, pp.413--427, 2002), and is based on a parallel computer code which uses mixed spatial discretization (spectral schemes in the homogeneous directions, and fourth-order, compact explicit finite-difference schemes in the radial direction). A DNS is carried out to simulate for the first time the turbulent Taylor--Couette flow in the turbulent regime. Statistical quantities are computed to complement the existing experimental information, with a view to compare it to planar, pressure-driven turbulent flow at the same value of the Reynolds number. The main source for differences in flow statistics between plane and curved-wall flows is attributed to the presence of large-scale rotating structures generated by curvature effects. △ Less

Submitted 9 November, 2007; originally announced November 2007.

Comments: To appear in European Journal of Mechanics B / Fluids

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

Authors: Paolo Luchini, Maurizio Quadrio, Simone Zuccher

Abstract: 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… ▽ More 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. △ Less

Submitted 21 November, 2006; originally announced November 2006.

Comments: Accepted for publication in Physics of Fluids

A low-cost parallel implementation of direct numerical simulation of wall turbulence

Authors: Paolo Luchini, Maurizio Quadrio

Abstract: A numerical method for the direct numerical simulation of incompressible wall turbulence in rectangular and cylindrical geometries is presented. The distinctive feature resides in its design being targeted towards an efficient distributed-memory parallel computing on commodity hardware. The adopted discretization is spectral in the two homogeneous directions; fourth-order accurate, compact finit… ▽ More A numerical method for the direct numerical simulation of incompressible wall turbulence in rectangular and cylindrical geometries is presented. The distinctive feature resides in its design being targeted towards an efficient distributed-memory parallel computing on commodity hardware. The adopted discretization is spectral in the two homogeneous directions; fourth-order accurate, compact finite-difference schemes over a variable-spacing mesh in the wall-normal direction are key to our parallel implementation. The parallel algorithm is designed in such a way as to minimize data exchange among the computing machines, and in particular to avoid taking a global transpose of the data during the pseudo-spectral evaluation of the non-linear terms. The computing machines can then be connected to each other through low-cost network devices. The code is optimized for memory requirements, which can moreover be subdivided among the computing nodes. The layout of a simple, dedicated and optimized computing system based on commodity hardware is described. The performance of the numerical method on this computing system is evaluated and compared with that of other codes described in the literature, as well as with that of the same code implementing a commonly employed strategy for the pseudo-spectral calculation. △ Less

Submitted 18 June, 2005; originally announced June 2005.

Comments: To be published in J. Comp. Physics