We use three-dimensional direct numerical simulations of homogeneous isotropic turbulence in a cubic domain to investigate the dynamics of heavy, chiral, finite-size inertial particles and their effects on the flow. Using an immersed-boundary method and a complex collision model, four-way coupled simulations have been performed, and the effects of particle-to-fluid density ratio, turbulence strength and particle volume fraction have been analysed. We find that freely falling particles on the one hand add energy to the turbulent flow but, on the other hand, they also enhance the flow dissipation: depending on the combination of flow parameters, the former or the latter mechanism prevails, thus yielding enhanced or weakened turbulence. Furthermore, particle chirality entails a preferential angular velocity which induces a net vorticity in the fluid phase. As turbulence strengthens, the energy introduced by the falling particles becomes less relevant and stronger velocity fluctuations alter the solid phase dynamics, making the effect of chirality irrelevant for the large-scale features of the flow. Moreover, comparing the time history of collision events for chiral particles and spheres (at the same volume fraction) suggests that the former tend to entangle, in contrast to the latter which rebound impulsively.

We present an experimental study of reactive control of turbulent jets, in which we target axisymmetric coherent structures, known to play a key role in the generation of sound. We first consider a forced jet, in which coherent structures are amplified above background levels, facilitating their detection, estimation and control. We then consider the more challenging case of an unforced jet. The linear control targets coherent structures in the region just downstream of the nozzle exit plane, where linear models are known to be appropriate for description of the lowest-order azimuthal modes of the turbulence. The control law is constructed in frequency space, based on empirically determined transfer functions. And the Wiener–Hopf formalism is used to enforce causality and to provide an optimal controller, as opposed to the sub-optimal control laws provided by simpler wave-cancellation methods. Significant improvements are demonstrated in the control of both forced and unforced jets. In the former case, order-of-magnitude reductions are achieved; and in the latter, turbulence levels are reduced by up to 60 %. The results open new perspectives for the control of turbulent flow at high Reynolds number.

The practical implementation of machine learning in flow control is limited due to its significant training expenses. In the present study the convolutional neural network (CNN) trained with the data of the restricted nonlinear (RNL) model is used to predict the normal velocity on a detection plane at $y^+=10$ in a turbulent channel flow, and the predicted velocity is used as wall blowing and suction for drag reduction. An active control test is carried out by using the well-trained CNN in direct numerical simulation (DNS). Substantial drag reduction rates up to 19 % and 16 % are obtained based on the spanwise and streamwise wall shear stresses, respectively. Furthermore, we explore the online control of wall turbulence by combining the RNL model with reinforcement learning (RL). The RL is constructed to determine the optimal wall blowing and suction based on its observation of the wall shear stresses without using the label data on the detection plane for training. The controlling and training processes are conducted synchronously in a RNL flow field. The control strategy discovered by RL has similar drag reduction rates with those obtained previously by the established method. Also, the training cost decreases by over thirty times at $Re_{\tau }=950$ compared with the DNS-RL model. The present results provide a perspective that combining the RNL model with machine learning control for drag reduction in wall turbulence can be effective and computationally economical. Also, this approach can be easily extended to flows at higher Reynolds numbers.

We optimize jet mixing using large eddy simulations (LES) at a Reynolds number of $3000$. Key methodological enablers consist of Bayesian optimization, a surrogate model enhanced by deep learning and persistent data topology for physical interpretation. The mixing performance is characterized by an equivalent jet radius ($R_{eq}$) derived from the streamwise velocity in a plane located $8$ diameters downstream. The optimization is performed in a 22-dimensional actuation space that comprises most known excitations. This search space parameterizes the distributed actuation imposed on the bulk flow and at the periphery of the nozzle in the streamwise and radial directions. The momentum flux measures the energy input of the actuation. The optimization quadruples the jet radius $R_{eq}$ with a $7$-armed blooming jet after around $570$ evaluations. The control input requires $2\,\%$ momentum flux of the main flow, which is one order of magnitude lower than an ad hoc dual-mode excitation. Intriguingly, a pronounced suboptimum in the search space is associated with a double-helix jet, a new flow pattern. This jet pattern results in a mixing improvement comparable to the blooming jet. A state-of-the-art Bayesian optimization converges towards this double-helix solution. The learning is accelerated and converges to another better optimum by including a deep-learning-enhanced surrogate model trained along the optimization. Persistent data topology extracts the global and many local minima in the actuation space. These minima can be identified with flow patterns beneficial to the mixing.

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. This provides a test condition in which the wall-normal temperature profile more realistically mimics that in an external flow, and 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.

This study demonstrates an experimental realization of turbulence control strategies previously explored by Choi et al. (J. Fluid Mech., vol. 262, 1994, pp. 75–110) through numerical simulations. To conduct the experiments, a deformable surface with a streamwise array of 16 independently controlled actuators was developed. A real-time particle image velocimetry (RT-PIV) system was also created for flow measurements. The objective of the control strategy was to target the sweep and ejection motions of the vortex shedding from a spherical cap placed in a laminar boundary layer. Reactive control strategies consisted of wall-normal surface deformations that opposed or complied with the wall-normal (v) or streamwise (u) velocity fluctuations obtained from the RT-PIV. The results showed two primary outcomes of the control approach. Firstly, it effectively hindered the advancement of sweep motions towards the wall. Secondly, it disrupted the periodic shedding of vortices. The v-control with opposing wall motions and u-control with compliant wall motions exhibited strong inhibition of sweep motions, while the v-control with compliant and u-control with opposing wall motions showed weaker inhibition. All reactive control cases resulted in the disruption of vortex shedding. In some instances, this disruption was accompanied by increased turbulent kinetic energy due to the generation of additional flow motions. However, the v-control with opposing wall motions significantly reduced the vortex-shedding energy while maintaining total turbulent kinetic energy close to or below that of the unforced flow. Overall, the experiments show the effectiveness of reactive control strategies in mitigating sweep motions and disrupting vortical structures, offering insights for developing reactive control strategies.

Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure–velocity coupling requires investigation to properly design a controller with such input data so that it can actuate upon the desired turbulent structures. A comprehensive database from direct numerical simulations (DNS) of turbulent channel flow is used for this purpose, spanning a Reynolds-number range $Re_\tau \approx 550\unicode{x2013}5200$. Spectral analysis reveals that the streamwise velocity is most strongly coupled to the linear term of the wall pressure, at a Reynolds-number invariant distance-from-the-wall scaling of $\lambda _x/y \approx 14$ (and $\lambda _x/y \approx 8$ for the wall-normal velocity). When extending the analysis to both homogeneous directions in $x$ and $y$, the peak coherence is centred at $\lambda _x/\lambda _z \approx 2$ and $\lambda _x/\lambda _z \approx 1$ for $p_w$ and $u$, and $p_w$ and $v$, respectively. A stronger coherence is retrieved when the quadratic term of the wall pressure is concerned, but there is only little evidence for a wall-attached-eddy type of scaling. An experimental dataset comprising simultaneous measurements of wall pressure and velocity complements the DNS-based findings at one value of $Re_\tau \approx 2$k, with ample evidence that the DNS-inferred correlations can be replicated with experimental pressure data subject to significant levels of (acoustic) facility noise. It is furthermore shown that velocity-state estimations can be achieved with good accuracy by including both the linear and quadratic terms of the wall pressure. An accuracy of up to 72 % in the binary state of the streamwise velocity fluctuations in the logarithmic region is achieved; this corresponds to a correlation coefficient of $\approx$0.6. This thus demonstrates that wall-pressure sensing for velocity-state estimation – e.g. for use in real-time control of wall-bounded turbulence – has merit in terms of its realization at a range of Reynolds numbers.

Eigenvectors of the observability and controllability Gramians represent responsive and receptive flow structures that enjoy a well-established connection to resolvent forcing and response modes. However, whereas resolvent modes have demonstrated great potential to guide sensor and actuator placement, observability and controllability modes have been leveraged exclusively in the context of model reduction via input and output projections. In this work, we introduce interpolatory, rather than orthogonal, input and output projections, that can be leveraged for sensor and actuator placement and open-loop control design. An interpolatory projector is an oblique projector with the property of preserving certain entries in the vector being projected. We review the connection between the resolvent operator and the Gramians, and present several numerical examples where we perform both orthogonal and interpolatory input and output projections onto the dominant forcing and response subspaces. Input projections are used to identify dynamically relevant disturbances, place sensors to measure disturbances, and place actuators for feedforward control in the linearized Ginzburg–Landau equation. Output projections are used to identify coherent structures and place sensors aiming at state reconstruction in the turbulent flow in a minimal channel at $Re_{\tau }=185$. The framework does not require data snapshots and relies only on knowledge of the steady or mean flow.

Turbulent shear flows are abundant in geophysical and astrophysical systems and in engineering-technology applications. They are often riddled with large-scale secondary flows that drastically modify the characteristics of the primary stream, preventing or enhancing mixing, mass and heat transfer. Using experiments and numerical simulations, we study the possibility of modifying these secondary flows by using superhydrophobic surface treatments that reduce the local shear. We focus on the canonical problem of Taylor–Couette flow, the flow between two coaxial and independently rotating cylinders, which has robust secondary structures called Taylor rolls that persist even at significant levels of turbulence. We generate these structures by rotating only the inner cylinder of the system, and show that an axially spaced superhydrophobic treatment can weaken the rolls through a mismatching surface heterogeneity, as long as the roll size can be fixed. The minimum hydrophobicity of the treatment required for this flow control is rationalized, and its effectiveness beyond the Reynolds numbers studied here is also discussed.

Turbulent drag reduction (DR) through streamwise travelling waves of the spanwise wall oscillation is investigated over a wide range of Reynolds numbers. Here, in Part 1, wall-resolved large-eddy simulations in a channel flow are conducted to examine how the frequency and wavenumber of the travelling wave influence the DR at friction Reynolds numbers $Re_\tau = 951$ and $4000$. The actuation parameter space is restricted to the inner-scaled actuation (ISA) pathway, where DR is achieved through direct attenuation of the near-wall scales. The level of turbulence attenuation, hence DR, is found to change with the near-wall Stokes layer protrusion height $\ell _{0.01}$. A range of frequencies is identified where the Stokes layer attenuates turbulence, lifting up the cycle of turbulence generation and thickening the viscous sublayer; in this range, the DR increases as $\ell _{0.01}$ increases up to $30$ viscous units. Outside this range, the strong Stokes shear strain enhances near-wall turbulence generation leading to a drop in DR with increasing $\ell _{0.01}$. We further find that, within our parameter and Reynolds number space, the ISA pathway has a power cost that always exceeds any DR savings. This motivates the study of the outer-scaled actuation pathway in Part 2, where DR is achieved through actuating the outer-scaled motions.

Direct numerical simulations of turbulent channel flow subjected to spanwise wall oscillations in the form of streamwise travelling waves (STW) were performed in an effort to elucidate the mechanism responsible for the observed drag reduction. We imposed large amplitudes to identify the proper effects of STW, while keeping the angular frequency and wavenumber fixed at a particular values. We primarily focus on the vorticity transport mechanism, to better understand the influence of STW actuation on the near-wall turbulence. We identify key terms appearing in the turbulent enstrophy transport equations that are directly linked to the STW actuation. The analysis reveals that the primary effect of the STW forcing is to attenuate the spanwise turbulent enstrophy at the wall, which is linked to the fluctuating wall shear stress. The suppression of the wall-normal turbulent enstrophy is deemed to be subordinate. To strengthen this point, we performed numerical experiments, where the streamwise fluctuating velocity, and consequently the spanwise vorticity, is artificially suppressed next to the wall. The anisotropic invariant maps show striking resemblance for large amplitude STW actuation and artificially forced cases. Detailed analysis of various structural features is provided, which includes the response of the near-wall streaks and shear layers of spanwise fluctuating velocity field. The quasistreamwise vortices, which play a key role in the regeneration mechanism, are shown to be pushed away from the wall, resulting in their weakened signature at the wall.

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.

The alteration of the near-wall flow field of a turbulent boundary layer flow subjected to spanwise travelling transversal surface waves at a friction Reynolds number $Re_\tau \approx 1525$ is investigated. The results of a spatial noise-assisted multivariate empirical mode decomposition reveal that this flow control method periodically induces near-wall large-scale bursts while simultaneously lowering the energetic content of small-scale features. The increasing occurrence of intense large-scale ejections in the near-wall region is of particular importance for reducing the wall-shear stress since these ejections balance large-scale sweeps originating from the outer layer. Thus, they corrupt the outer-layer impact on the near-wall dynamics and, consequently, the overall fluctuation intensity at the wall is attenuated. This disturbed top-down momentum exchange is highlighted by an inner–outer interaction analysis, which further reveals an increased bottom-up communication provoked by the large-scale ejections. Moreover, it is shown that the periodic secondary flow field induced by the actuation interferes with the quasi-streamwise vortices in the near-wall region. The velocity gradients of the secondary flow field deform the vortices’ cross-section into an elliptic shape, which yields an unstable vortex state resulting in vortex disintegration. In combination with the effect of the large-scale ejections, the reduced number of quasi-streamwise vortices compared with the undisturbed boundary layer flow results in a decreased wall-normal momentum exchange and the widening and weakening of near-wall streaks. This yields a reduced fluctuation intensity in the near-wall region that lowers the overall wall-shear stress level.

We investigate the Reynolds analogy over riblets, namely the analogy between the fractional increase in Stanton number $C_h$ and the fractional increase in the skin-friction coefficient $C_f$, relative to a smooth surface. We investigate the direct numerical simulation data of Endrikat et al. (Flow Turbul. Combust., vol. 107, 2021, pp. 1–29). The riblet groove shapes are isosceles triangles with tip angles $\alpha = {30}^{\circ }, {60}^{\circ }, {90}^{\circ }$, a trapezoid, a rectangle and a right triangle. The viscous-scaled riblet spacing varies between $s^+ \approx 10$ to $60$. The global Reynolds analogy is primarily influenced by Kelvin–Helmholtz rollers and secondary flows. Kelvin–Helmholtz rollers locally break the Reynolds analogy favourably, i.e. cause a locally larger fractional increase in $C_h$ than in $C_f$. These rollers induce negative wall shear stress patches which have no analogue in wall heat fluxes. Secondary flows at the riblets’ crests are associated with local unfavourable breaking of the Reynolds analogy, i.e. locally larger fractional increase in $C_f$ than in $C_h$. Only the triangular riblets with $\alpha = {30}^{\circ }$ trigger strong Kelvin–Helmholtz rollers without appreciable secondary flows. This riblet shape globally preserves the Reynolds analogy from $s^+ = 21$ to $33$. However, the other riblet shapes have weak or non-existent Kelvin–Helmholtz rollers, yet persistent secondary flows. These riblet shapes behave similarly to rough surfaces. They unfavourably break the global Reynolds analogy, and do so to a greater extent as $s^+$ increases.

A set of boundary conditions called the transpiration-resistance model (TRM) is investigated in altering near-wall turbulence. The TRM was proposed by Lacis et al. (J. Fluid Mech., vol. 884, 2020, p. A21) as a means of representing the net effect of surface micro-textures on their overlying bulk flows. It encompasses conventional Navier-slip boundary conditions relating the streamwise and spanwise velocities to their respective shears through the slip lengths $\ell _x$ and $\ell _z$. In addition, it features a transpiration condition accounting for the changes induced in the wall-normal velocity by expressing it in terms of variations of the wall-parallel velocity shears through the transpiration lengths $m_x$ and $m_z$. Greater levels of drag increase occur when more transpiration takes place at the boundary plane, with turbulent transpiration being predominately coupled to the spanwise shear component for canonical near-wall turbulence. The TRM reproduces the effect of a homogeneous and structured roughness up to $k^{+}\approx 18$, encompassing the regime of smooth-wall-like turbulence described using virtual origins (Luchini, 1996 Reducing the turbulent skin friction. In Computational Methods in Applied Sciences’ 96 (Paris, 9–13 Sept. 1996), pp. 465–470. Wiley; Ibrahim et al., J. Fluid Mech., vol. 915, 2021, p. A56) and slightly beyond it. The transpiration factor is defined as the product of the slip and transpiration lengths, i.e. $(m\ell )_{x,z}$. This factor contains the compound effect of the wall-parallel velocity occurring at the boundary plane and increased permeability, both of which lead to the transport of momentum in the wall-normal direction. A linear relation between the transpiration factor and the roughness function is observed for regularly textured surfaces in the transitionally rough regime of turbulence. This shows that such effective flow quantities can be suitable measures for characterizing rough surfaces in this flow regime.

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 billion 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.

A turbulent drag-reduction method employing synthetic jet sheets in a turbulent channel flow is investigated by direct numerical simulations. The jet sheets are wall-parallel and produced by periodic blowing and suction from pairs of thin slots aligned with the main streamwise flow. By varying the slot height and the jet-sheet angle with respect to the spanwise direction, drag-reduction margins between $10\,\%$ and $30\,\%$ are obtained for jet-sheet angles between $45 ^{\circ }$ and $75 ^{\circ }$, while a drag increase of almost 100 % is computed when the jet sheets are spanwise-oriented. When global skin-friction drag reduction occurs, the wall-shear stress near the jet-sheet exits increases during suction and decreases during blowing, while the velocity fluctuations weaken during suction and intensify during blowing. The global drag-reduction effect is produced by a finite counter flow induced by the nonlinear interaction between the jet-sheet flow and the main flow, although the turbulent intensity and Reynolds shear stresses increase. The power spent to generate the jet sheets is computed by modelling numerically the actuator underneath the channel flow as a piston oscillating sinusoidally along the spanwise direction in a round-shaped cavity from which the fluid is released into the channel through the cavity exits. A power balance leads to the computation of the efficiency of the actuator system, quantifying the portion of the piston power that is lost as internal power fluxes and heat transfer through the cavity walls. For the tested configurations, the power consumed by the piston to generate the jet sheets is larger than the power saved thanks to the drag reduction.

We present an experimental and numerical study of turbulent thermal convection in the presence of an effective horizontal buoyancy that generates extra shear at the boundary. Geometrical confinements are also applied by varying the streamwise and spanwise aspect ratios of the convection cell to condense the plumes. With these, we systematically explore the effects of plume and shear on heat transfer. It is found that a streamwise confinement results in increased plume coverage but decreased shear compared with spanwise confinement. The fact that streamwise confinement leads to a higher vertical heat transfer efficiency than the spanwise confined case suggests that the increase of plume coverage is the dominant effect responsible for the enhanced heat transfer. Our results highlight the potential applications of coherent structure manipulation in efficient passive heat transfer control and thermal engineering. We also analyse the energetics of the present system and derive the expression of mixing efficiency accordingly. The mixing efficiency is found to increase with both the buoyancy ratio and streamwise dimension.

We conduct a well-controlled model experiment for a wide variety of canopy flows. Examples of these include engineering flows such as wind flow, dispersion of scalars through and over urban areas, and the convective heat transfer in many heat exchangers, as well as natural canopies such as flows through terrestrial or aquatic vegetation. We aim to shed the light on fundamental flow and transport phenomena common to these applications. Specifically, the characteristics of mean flow and scalar concentration characteristics of a turbulent boundary layer flow impinging on a canopy, which comprises a cluster of tall obstacles (this can also be interpreted as a porous obstruction). The cluster is created with a group of cylinders of diameter $d$ and height $h$ arranged in a circular patch of diameter $D$. The solidity of the patch/obstruction is defined by $\phi$ (the total planar area covered by cylinders), which is systematically varied ($0.098 \leq \phi \leq 1$) by increasing the number of cylinders in a patch ($N_c$). A point source is placed at ground level upstream of the patch and its transport over and around the patch is examined. Time-averaged velocity and scalar fields, obtained from simultaneous planar particle image velocimetry-planar laser-induced fluorescence (PIV-PLIF) measurements, reveal that the characteristics of wake and flow above porous patches are heavily influenced by $\phi$. In particular, we observe that the horizontal and vertical extent of the wake and scalar concentration downstream of the patches decreases and increases with $\phi$, respectively. Here, the recirculation bubble is shifted closer to the trailing edge (TE) of the patches as $\phi$ increases, limiting the flow from convecting downstream, decreasing the scalar concentration and virtually ‘extending’ the patch in the streamwise direction. As the bubble forms in the TE, vertical bleeding increases and hence the concentration increases above the patch where the cylinders appear to ‘extend’ vertically towards the freestream.

Opposition flow control is a robust strategy that has been proved effective in turbulent wall-bounded flows. Its conventional set-up consists of measuring wall-normal velocity in the buffer layer and opposing it at the wall. This work explores the possibility of implementing this strategy with a detection plane in the logarithmic layer, where control could be feasible experimentally. We apply control on a channel flow at $Re_\tau = 932$, only on the eddies with relatively large wavelengths ($\lambda / h > 0.1$). Similarly to the buffer layer opposition control, our control strategy results in a virtual-wall effect for the wall-normal velocity, creating a minimum in its intensity. However, it also induces a large response in the streamwise velocity and Reynolds stresses near the wall, with a substantial drag increase. When the phase of the control lags with respect to the detection plane, spanwise-homogeneous rollers are observed near the channel wall. We show that they are a result of a linear instability. In contrast, when the control leads with respect to the detection plane, this instability is inactive and oblique waves are observed. Their wall-normal profiles can be predicted linearly as a response of the turbulent channel flow to a forcing with the advection velocity of the detection plane. The linearity, governing the flow, opens a possibility to affect large scales of the flow in a controlled manner, when enhanced turbulence intensity or mixing is desired.