2.1. Triple decomposition of the velocity field

An incompressible turbulent flow, varying in space $\boldsymbol {x}$ and time $t$, is typically described via its mean and fluctuating velocity and pressure fields, defined after the classic Reynolds decomposition. Provided the flow exhibits well-defined non-stochastic (e.g. periodic) features, the fluctuating field can be further decomposed into a coherent and a stochastic part. Therefore, the velocity field reads

(2.1)\begin{equation} \boldsymbol{u} = \boldsymbol{U} + \underbrace{ \tilde{ {\boldsymbol{u}} } + \boldsymbol{u}''}_{\boldsymbol{u}'}, \end{equation}

where $\boldsymbol {U}$, $\boldsymbol {u}'$, $\tilde { {\boldsymbol {u}} }$ and $\boldsymbol {u}''$ indicate the mean, fluctuating, coherent and stochastic parts of the velocity field $\boldsymbol {u}$. The mean velocity $\boldsymbol {U}$ is defined as $\boldsymbol {U} \equiv \left\langle {\boldsymbol {u}} \right \rangle$, with the operator $\left\langle {\cdot } \right \rangle$ indicating ensemble averaging, which under the ergodic hypothesis becomes equivalent to averaging over homogeneous directions and time (if the flow is statistically stationary). For a single realisation without homogeneous directions, the mean is simply a temporal average:

(2.2)\begin{equation} \boldsymbol{U}(\boldsymbol{x}) \equiv \lim_{\tau \to +\infty} \frac{1}{\tau} \int_0^\tau \boldsymbol{u}(\boldsymbol{x},t) \,{\rm d}t . \end{equation}

Considering a periodic motion with period $T$ and phase $\varphi \in (0,2{\rm \pi} ]$, the overbar $\overline {\cdot}$ denotes the phase average operator over an integer number $N$ of periods. Like $\left\langle {\cdot } \right \rangle$, it includes averaging over the homogeneous directions. Considering again a single realisation without homogeneous directions, $\overline {\cdot}$ is defined as

(2.3)\begin{equation} \overline{{\boldsymbol{u}}}(\boldsymbol{x},\varphi) \equiv \lim_{N\to +\infty} \frac{1}{N} \sum_{n = 0}^{N-1} \boldsymbol{u}\left(\boldsymbol{x}, \left( \frac{\varphi}{2 {\rm \pi}} + n \right) T \right). \end{equation}

The coherent field $\tilde { {\boldsymbol {u}} }$ is thus defined as

(2.4)\begin{equation} \tilde{ {\boldsymbol{u}} } (\boldsymbol{x},\varphi) = \overline{{\boldsymbol{u}}}(\boldsymbol{x},\varphi) - \boldsymbol{U}(\boldsymbol{x}) , \end{equation}

and the stochastic vector field $\boldsymbol {u}''$ is defined after the triple decomposition (2.1) as $\boldsymbol {u}''=\boldsymbol {u}-\boldsymbol {U}- \tilde { {\boldsymbol {u}} }$. An analogous triple decomposition is used to decompose the pressure field $p=P+ \tilde { {p} } +p''$, with $\tilde { {p} } +p''=p'$.

2.2. The AGKE

Before presenting the $\varphi$AGKE, the standard AGKE based on the Reynolds decomposition are recalled. Full details of their derivation from the incompressible Navier–Stokes equations are provided by Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020).

Exact budget equations can be written for the components of the second-order structure function tensor $\big\langle {\delta u_i \delta u_j}\big\rangle $, where $\delta u_i=u_i(\boldsymbol {X}+\boldsymbol {r}/2,t)-u_i(\boldsymbol {X}-\boldsymbol {r}/2,t)$ is the $i$th component of the velocity difference between two points $\boldsymbol {x}_1$ and $\boldsymbol {x}_2$, identified by their midpoint $\boldsymbol {X} = (\boldsymbol {x}_1+\boldsymbol {x}_2)/2$ and their separation vector $\boldsymbol {r} = (\boldsymbol {x}_2-\boldsymbol {x}_1)$, as shown by the sketch in figure 1. The Reynolds decomposition leads to budget equations for $\delta U_i \delta U_j$ and $\big\langle {\delta u'_i \delta u'_j} \big\rangle $. In general, the time-independent tensor $\delta U_i \delta U_j$ depends upon six independent variables, i.e. the six coordinates of $\boldsymbol {X}$ and $\boldsymbol {r}$. The tensor $\big\langle {\delta u'_i \delta u'_j} \big\rangle$ additionally features time $t$ as an independent variable if the process is not statistically stationary (e.g. periodic), and is related to the Reynolds stresses $\big\langle {u'_i u'_j} \big\rangle$ and to the spatial correlation tensor ${\mathsf{R}}_{ij}$ (Davidson et al. Reference Davidson, Nickels and Krogstad2006; Agostini & Leschziner Reference Agostini and Leschziner2017) as

(2.5)\begin{equation} \big\langle {\delta u'_i \delta u'_j} \big\rangle(\boldsymbol{X},\boldsymbol{r},t) = {\mathsf{V}}_{ij}(\boldsymbol{X},\boldsymbol{r},t) - {\mathsf{R}}_{ij}(\boldsymbol{X},\boldsymbol{r},t) - {\mathsf{R}}_{ij}(\boldsymbol{X},-\boldsymbol{r},t), \end{equation}

where

(2.6)\begin{equation} {\mathsf{V}}_{ij}(\boldsymbol{X},\boldsymbol{r},t)=\big\langle {u'_i u'_j} \big\rangle \left( \boldsymbol{X} + \frac{\boldsymbol{r}}{2}, t \right) + \big\langle {u'_i u'_j} \big\rangle \left( \boldsymbol{X} - \frac{\boldsymbol{r}}{2}, t \right) \end{equation}

is the sum of the single-point Reynolds stresses evaluated at the two points $\boldsymbol {X} \pm \boldsymbol {r}/2$ and

(2.7)\begin{equation} {\mathsf{R}}_{ij}(\boldsymbol{X},\boldsymbol{r},t) = \left\langle { u'_i \left(\boldsymbol{X}+\frac{\boldsymbol{r}}{2},t \right) u'_j \left(\boldsymbol{X}-\frac{\boldsymbol{r}}{2},t \right)} \right \rangle \end{equation}

is the two-point spatial correlation function.

Figure 1. Sketch of two points $\boldsymbol {x}_1$ and $\boldsymbol {x}_2$ involved in the definition of the second-order structure function tensor. Here, $\boldsymbol {X}=( \boldsymbol {x}_1 + \boldsymbol {x}_2 )/2$ and $\boldsymbol {r}=\boldsymbol {x}_2-\boldsymbol {x}_1$ indicate their mid-point and separation vector, respectively, and $\delta \boldsymbol {u} = \boldsymbol {u}_2 - \boldsymbol {u}_1$ is the velocity increment between the two points.

The budget equations for the components of the mean second-order structure function tensor $\delta U_i \delta U_j$ are presented here for the first time; they were not reported by Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020), and the tensor has received little attention so far, owing to its irrelevance in homogeneous isotropic turbulence, where there is no mean flow. The mean AGKE are written compactly as

(2.8)\begin{equation} \frac{\partial \varPhi^m_{k,ij}}{\partial r_k} + \frac{\partial \varPsi^m_{k,ij}}{\partial X_k} = \varXi^m_{ij} , \end{equation}

where the repeated index $k$ implies summation. The following notation is adopted. Uppercase letters (e.g. $\varPhi$, $\varPsi$ and $\varXi$) are used to denote time-averaged quantities and lowercase letters (e.g. $\phi$, $\psi$ and $\xi$) are used for phase-dependent quantities. Furthermore, superscripts $m$, $f$, $c$ and $s$ are used to label terms in the budget equations for the mean structure function tensor $\delta U_i \delta U_j$, the fluctuating structure function tensor $\big\langle {\delta u'_i \delta u'_j} \big\rangle $, the coherent structure function tensor $\overline {{\delta \tilde{u} _i \delta \tilde{u} _j}} = \delta \tilde{u} _i \delta \tilde{u} _j$ and the stochastic structure function tensor $\overline {{\delta u''_i \delta u''_j}}$.

The fluxes $\varPhi ^m_{k,ij}$ and $\varPsi ^m_{k,ij}$ are the mean scale- and physical-space fluxes, i.e.

(2.9)\begin{equation} \varPhi^m_{k,ij} =\underbrace{ \delta U_k \delta U_i \delta U_j }_{ \text{Mean transport} } + \underbrace{ \delta U_j \big\langle { \delta u'_k \delta u'_i } \big\rangle + \delta U_i \big\langle { \delta u'_k \delta u'_j } \big\rangle }_{ \text{Fluctuating transport} } \underbrace{ - 2 \nu \frac{ \partial \delta U_i \delta U_j}{ \partial r_k} }_{ \text{Viscous diffusion} } \quad k = 1,2,3 \end{equation}

and

(2.10)$$\begin{align} \varPsi^m_{k,ij} &=\underbrace{ U_k^* \delta U_i \delta U_j }_{ \text{Mean transport} } + \underbrace{ \delta U_j \big\langle { u'^*_k \delta u'_i } \big\rangle + \delta U_i \big\langle { u'^*_k \delta u'_j } \big\rangle }_{ \text{Fluctuating transport} } \underbrace{+ \frac{1}{\rho} \delta P \delta U_j \delta_{ki} + \frac{1}{\rho} \delta P \delta U_i \delta_{kj} }_{ \text{Pressure transport} } \nonumber\\ &\quad \underbrace{ - \frac{\nu}{2} \frac{ \partial \delta U_i \delta U_j}{ \partial X_k} }_{ \text{Viscous diffusion} } \quad k = 1,2,3, \end{align}$$

where $\delta _{ij}$ is the Kronecker delta, $\nu$ is the kinematic viscosity and the asterisk superscript $(\cdot )^{\ast }$ indicates the arithmetic average of a quantity between the two points $\boldsymbol {X} \pm \boldsymbol {r}/2$. The term $\varXi ^m_{ij}$, instead, is the mean source and reads

(2.11) $$\begin{align} \varXi^m_{ij} &=- \underbrace{\left[- \big\langle {u^{\prime *}_k \delta u'_j } \big\rangle \delta \left( \frac{ \partial U_i }{ \partial x_k } \right) -\left\langle {u^{\prime *}_k \delta u'_i} \right \rangle \delta \left( \frac{ \partial U_j }{ \partial x_k } \right) -\big\langle {\delta u'_k \delta u'_j } \big\rangle \left( \frac{\partial U_i}{\partial x_k} \right)^*- \left\langle {\delta u'_k \delta u'_i } \right \rangle \left( \frac{\partial U_j}{\partial x_k} \right)^* \right]}_{\text{Mean-fluctuating production} \ ({\mathsf{P}}^{mf}_{ij})}\nonumber\\ &\quad \underbrace{ + \frac{1}{\rho} \delta P \frac{\partial \delta U_i}{\partial X_j} + \frac{1}{\rho} \delta P \frac{\partial \delta U_j}{\partial X_i} }_{\text{Pressure strain} \ (\varPi^m_{ij})} \underbrace{- 4 \epsilon_{ij}^{m*} }_{\text{Dissipation} \ ({\mathsf{D}}^m_{ij})} + \underbrace{\delta U_j \delta F_i + \delta U_i \delta F_j}_{\text{Forcing interaction} \ ({\mathsf{F}}^m_{ij})}. \end{align}$$

The standard AGKE, presented by Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020), pertain to increments of the fluctuating velocity field, and describe the production, transport, redistribution and dissipation of each component, in physical space $\boldsymbol {X}$ and in the space of scales $\boldsymbol {r}$. They can be written compactly as

(2.12)\begin{equation} \frac{\partial \big\langle {\delta u'_i \delta u'_j} \big\rangle }{\partial t} + \frac{\partial \varPhi^f_{k,ij}}{\partial r_k} + \frac{\partial \varPsi^f_{k,ij}}{\partial X_k} = \varXi^f_{ij} . \end{equation}

The scale-space fluxes $\varPhi ^f_{k,ij}$ and physical-space fluxes $\varPsi ^f_{k,ij}$ are defined as

(2.13)\begin{equation} \varPhi^f_{k,ij} =\underbrace{\big\langle {\delta U_k \delta u'_i \delta u'_j } \big\rangle }_{\text{Mean transport}} + \underbrace{\big\langle {\delta u'_k \delta u'_i \delta u'_j} \big\rangle }_{\text{Fluctuating transport}} \underbrace{ -2 \nu \frac{\partial}{\partial r_k} \big\langle {\delta u'_i \delta u'_j} \big\rangle}_{\text{Viscous diffusion}} \quad k=1,2,3 \end{equation}

and

(2.14)$$\begin{align} \varPsi^f_{k,ij} &= \underbrace{ \big\langle { U_k^{{\ast}} \delta u'_i \delta u'_j} \big\rangle }_{\text{Mean transport}}+ \underbrace{ \big\langle {u_k'^{{\ast}} \delta u'_i \delta u'_j} \big\rangle }_{\text{Fluctuating transport}} + \underbrace{\frac{1}{\rho} \big\langle {\delta p' \delta u'_i} \big\rangle \delta_{kj} + \frac{1}{\rho} \big\langle {\delta p' \delta u'_j} \big\rangle \delta_{ki}}_{\text{Pressure transport}}\nonumber\\ &\quad\underbrace{- \frac{\nu}{2} \frac{\partial}{\partial X_k} \big\langle {\delta u'_i \delta u'_j} \big\rangle}_{\text{Viscous diffusion}} \quad k=1,2,3. \end{align}$$

The term $\varXi ^f_{ij}$ in (2.12) is the source for $\big\langle { \delta u'_i \delta u'_j} \big\rangle$ and reads

(2.15)$$\begin{gather} \varXi^f_{ij} = \underbrace{-\big\langle {u_k'^{{\ast}} \delta u'_j} \big\rangle \delta \left( \frac{ \partial U_i}{\partial x_k} \right) - \big\langle {u_k'^{{\ast}} \delta u'_i} \big\rangle \delta \left(\frac{\partial U_j}{ \partial x_k}\right) -\big\langle {\delta u'_k \delta u'_j} \big\rangle \left(\frac{ \partial U_i}{ \partial x_k}\right)^{{\ast}} - \big\langle {\delta u'_k \delta u'_i} \big\rangle \left( \frac{\partial U_j}{ \partial x_k} \right)^{{\ast}}}_{\text{Mean-fluctuating production } ({\mathsf{P}}^{mf}_{ij})} \nonumber\\ + \underbrace{\frac{1}{\rho}\left\langle {\delta p' \frac{ \partial \delta u'_i}{ \partial X_j}} \right \rangle + \frac{1}{\rho}\left\langle {\delta p' \frac{\partial \delta u'_j}{ \partial X_i}} \right \rangle}_{\text{Pressure strain } (\varPi^f_{ij})} \underbrace{-4 \epsilon_{ij}^{f \ast} }_{\text{Dissipation } ({\mathsf{D}}^f_{ij})} + \underbrace{\big\langle {\delta u'_j \delta f'_i} \big\rangle + \big\langle {\delta u'_i \delta f'_j} \big\rangle}_{\text{Forcing interaction } ({\mathsf{F}}^f_{ij})}, \end{gather}$$

in which $\epsilon _{ij}^f$ is the pseudo-dissipation tensor $\big\langle {\partial u_i'/\partial x_k \partial u_j' / \partial x_k} \big \rangle $. The source term $\varXi ^f_{ij}$ identifies scales and positions with a net sink ($\varXi ^f_{ij}<0$) or a net source ($\varXi ^f_{ij} > 0$) for each component of the Reynolds stresses. The separation of $\varXi ^f_{ij}$ into its constituent terms provides insight into mean-fluctuating production ${\mathsf{P}}^{mf}_{ij}$ (which also appears in (2.11) with opposite sign), redistribution $\varPi ^f_{ij}$, dissipation ${\mathsf{D}}^f_{ij}$ and interaction with external fluctuating volume forces ${\mathsf{F}}^f_{ij}$ of turbulent stresses among scales and positions (note that the forcing interaction term was missing in the original AGKE formulated by Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020)). The flux vectors describe the various transfer processes, and their field lines visualise how fluctuations are transferred among scales and positions, via direct and inverse cascades. It should be recalled that, as stressed by Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020), when interpreting AGKE results to extract structural turbulence information, local peaks of the structure functions always need to be connected to local maxima/minima of the correlation functions whenever a separation along an inhomogeneous direction is involved.

2.3. The phase-aware AGKE ($\varphi$AGKE)

By using the triple decomposition (2.1), the phase-averaged fluctuating structure function tensor $\overline {{\delta u'_i \delta u'_j}}(\boldsymbol {X},\boldsymbol {r},\varphi )$ can be separated into its coherent and stochastic parts, i.e.

(2.16)\begin{equation} \overline{{\delta u_i' \delta u_j'}}(\boldsymbol{X},\boldsymbol{r},\varphi) = \overline{{\delta \tilde{ {u} } _i \delta \tilde{ {u} } _j}} (\boldsymbol{X},\boldsymbol{r},\varphi) + \overline{{\delta u''_i \delta u''_j}} (\boldsymbol{X},\boldsymbol{r},\varphi); \end{equation}

note that $\overline {{\delta \tilde{u} _i \delta \tilde { {u_j} } }} \equiv \delta \tilde{u} _i \delta \tilde { {u_j} }$ owing to the definition of the phase-average operator. Two budget equations, called $\varphi$AGKE, can be written for $\delta \tilde{u} _i \delta \tilde{u} _j$ and $\overline {{\delta u''_i \delta u''_j}}$, which include, unlike the standard AGKE, the interplay among the mean, coherent and stochastic fields at each phase $\varphi$. These new equations extend in a significant way the work of Thiesset et al. (Reference Thiesset, Danaila and Antonia2014) and Alves Portela et al. (Reference Alves Portela, Papadakis and Vassilicos2020), which considered the budget equations for $\left\langle {\delta \tilde{u} _i \delta \tilde{u} _i} \right \rangle(\boldsymbol {X},\boldsymbol {r})$ and $\left\langle {\delta u''_i \delta u''_i} \right \rangle(\boldsymbol {X},\boldsymbol {r})$. They applied the triple decomposition to the trace $\left\langle {\delta u'_i \delta u'_i} \right \rangle$ of the second-order structure function tensor, instead of considering the whole tensor. The major difference, though, is that the dependence on the phase $\varphi$ of the coherent motion (or external forcing) was lost, because of the use of the $\left\langle {\cdot } \right \rangle$ operator. On the contrary, the $\varphi$AGKE retain full phase information.

The step-by-step derivation of the $\varphi$AGKE from the incompressible Navier–Stokes equations is described in Appendix A. At each phase $\varphi$, they link the phase variation of each component of the coherent and stochastic structure function tensors, at a given scale $\boldsymbol {r}$ and position $\boldsymbol {X}$, to the unbalance among inter-component redistribution, scale-space transport, dissipation and mean–coherent–stochastic interaction. The last term is obviously absent in the classic AGKE.

The equations for the coherent and stochastic parts can be compactly written as

(2.17)\begin{equation} \frac{2{\rm \pi}}{T} \frac{\partial \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }{ \partial \varphi} + \frac{\partial \phi^c_{k,ij}}{\partial r_k} + \frac{\partial \psi^c_{k,ij}}{\partial X_k} = \xi^c_{ij} + \zeta^c_{ij} \end{equation}

and

(2.18)\begin{equation} \frac{2{\rm \pi}}{T} \frac{\partial \overline{{\delta u''_i \delta u''_j}}}{\partial \varphi} + \frac{\partial \phi^s_{k,ij}}{\partial r_k} + \frac{\partial \psi^s_{k,ij}}{\partial X_k} = \xi^s_{ij} , \end{equation}

where, as above, the repeated index $k$ implies summation.

The first term in each of (2.17) and (2.18) represents the phase variation of the coherent and stochastic components of the structure function tensor. The coherent and stochastic scale fluxes $\phi ^c_{k,ij}$ and $\phi ^s_{k,ij}$, i.e. the fluxes of $\delta \tilde{u} _i \delta \tilde{u} _j$ and $\overline {{\delta u''_i \delta u''_j}}$ in the space of scales, are defined as

(2.19)\begin{equation} \phi^c_{k,ij}=\underbrace{ \delta U_k \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }_{\text{Mean transport}} + \underbrace{ \delta \tilde{ {u} } _k \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }_{ \text{Coherent transport} } + \underbrace{\overline{{ \delta u''_k \delta u''_i }} \delta \tilde{ {u} } _j + \overline{{ \delta u''_k \delta u''_j }} \delta \tilde{ {u} } _i }_{ \text{Stochastic transport} } \underbrace{-2 \nu \frac{ \partial \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }{\partial r_k} }_{ \text{Viscous diffusion} } \quad k = 1,2,3 \end{equation}

and

(2.20)\begin{gather} \phi^s_{k,ij} =\underbrace{ \delta U_k \overline{{ \delta u''_i \delta u''_j}} }_{ \text{Mean transport} } + \underbrace{ \delta \tilde{ {u} } _k \overline{{ \delta u''_i \delta u''_j }} }_{ \text{Coherent transport} } + \underbrace{ \overline{{\delta u''_k \delta u''_i \delta u''_j}} }_{ \text{Stochastic transport} } \underbrace{ - 2 \nu \frac{ \partial \overline{{ \delta u''_i \delta u''_j}}}{ \partial r_k} }_{ \text{Viscous diffusion} } \quad k = 1,2,3. \end{gather}

The coherent and stochastic spatial flux terms $\psi ^c_{k,ij}$ and $\psi ^s_{k,ij}$, i.e. the fluxes of $\delta \tilde{u} _i \delta \tilde{u} _j$ and $\overline {{\delta u''_i \delta u''_j}}$ in physical space, are defined as

(2.21)$$\begin{align} \psi^c_{k,ij} &=\underbrace{ U_k^* \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }_{ \text{Mean transport} } + \underbrace{ \tilde{ {u} } _k^* \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }_{ \text{Coherent transport} } + \underbrace{ \overline{{ u^{\prime \prime * }_k \delta u^{\prime \prime}_i }} \delta \tilde{ {u} } _j + \overline{{ u^{\prime \prime *}_k \delta u''_j }} \delta \tilde{ {u} } _i }_{ \text{Stochastic transport} } + \underbrace{ \frac{1}{\rho} \delta \tilde{ {p} } \delta \tilde{ {u} } _i \delta_{kj} }_{ \text{Pressure transport} } \nonumber\\ &\quad \underbrace{+ \frac{1}{\rho} \delta \tilde{ {p} } \delta \tilde{ {u} } _j \delta_{ki} }_{ \text{Pressure transport} } \underbrace{ - \frac{\nu}{2} \frac{\partial \delta \tilde{ {u} } _i \delta \tilde{ {u} } _j }{ \partial X_k } }_{ \text{Viscous diffusion} } \quad k = 1,2,3, \end{align}$$

(2.22)$$\begin{align}\psi^s_{k,ij} &=\underbrace{ U_k^* \overline{{ \delta u''_i \delta u''_j }} }_{ \text{Mean transport} } + \underbrace{ \tilde{ {u} } _k^* \overline{{\delta u''_i \delta u''_j}} }_{ \text{Coherent transport} } + \underbrace{ \overline{{ u_k^{\prime \prime * } \delta u''_i \delta u''_j }} }_{ \text{Stochastic transport} } + \underbrace{ \frac{1}{\rho}\overline{{ \delta p'' \delta u''_i}} \delta_{kj} + \frac{1}{\rho}\overline{{ \delta p'' \delta u''_j}} \delta_{ki} }_{ \text{Pressure transport} }\nonumber\\ &\quad \underbrace{ - \frac{\nu}{2} \frac{\partial \overline{{\delta u''_i \delta u''_j}}}{\partial X_k } }_{ \text{Viscous diffusion} }\quad k = 1,2,3. \end{align}$$

The differences with the fluxes (2.13) and (2.14) appearing in the standard AGKE are worth noticing. Two new terms appear here to account for the effect of the coherent field upon transport in the stochastic field, labelled as coherent transport in (2.20) and (2.22). Vice versa, how the stochastic field affects transport in the coherent field is reflected by the stochastic transport term in (2.19) and (2.21).

The coherent and stochastic source terms $\xi ^c_{ij}$ and $\xi ^s_{ij}$ denote the scale-space net production of $\delta \tilde{u} _i \delta \tilde{u} _j$ and $\overline {{\delta u''_i \delta u''_j}}$. They can be either positive or negative, and read

(2.23) \begin{align} \xi^c_{ij} & =\underbrace{- \delta \tilde{ {u} } _j \delta \tilde{ {u} } _k \left( \frac{\partial U_i}{\partial x_k} \right)^* - \delta \tilde{ {u} } _i \delta \tilde{ {u} } _k \left( \frac{\partial U_j}{\partial x_k} \right)^* - \delta \tilde{ {u} } _j \tilde{ {u} } _k^* \delta \left( \frac{ \partial U_i }{ \partial x_k } \right) - \delta \tilde{ {u} } _i \tilde{ {u} } _k^* \delta \left( \frac{ \partial U_j }{ \partial x_k } \right)}_{\text{Mean-coherent production} \ ({\mathsf{p}}^{mc}_{ij})}\nonumber\\ & \quad - \underbrace{\left[- \overline{{ \delta u''_j \delta u''_k}} \left( \frac{\partial \tilde{ {u} } _i}{\partial x_k} \right)^* - \overline{{ \delta u''_i \delta u''_k}} \left( \frac{\partial \tilde{ {u} } _j}{\partial x_k} \right)^* - \overline{{ \delta u''_j u^{\prime \prime *}_k }} \delta \left( \frac{ \partial \tilde{ {u} } _i }{\partial x_k} \right) - \overline{{ \delta u''_i u^{\prime \prime *}_k }} \delta \left( \frac{ \partial \tilde{ {u} } _j }{\partial x_k} \right) \right] }_{\text{Coherent-stochastic production} \ ({\mathsf{p}}^{cs}_{ij})} \nonumber\\ &\quad \underbrace{+ \frac{1}{\rho} \delta \tilde{ {p} } \frac{\partial \delta \tilde{ {u} } _i}{\partial X_j} + \frac{1} {\rho} \delta \tilde{ {p} } \frac{\partial \delta \tilde{ {u} } _j}{\partial X_i} }_{\text{Pressure strain} \ ({\rm \pi}^c_{ij})} \underbrace{- 4 \epsilon_{ij}^{c \ast} }_{\text{Dissipation} ({\mathsf{d}}^c_{ij})} + \underbrace{\delta \tilde{ {u} } _j \delta \tilde{ {f} } _i + \delta \tilde{ {u} } _i \delta \tilde{ {f} } _j}_{\text{Forcing interaction} ({\mathsf{f}}^c_{ij})}, \end{align}

(2.24)\begin{align} \xi^s_{ij} & =\underbrace{- \overline{{ \delta u''_j \delta u''_k }} \left( \frac{\partial U_i}{\partial x_k} \right)^* - \overline{{ \delta u''_i \delta u''_k }} \left( \frac{\partial U_j}{\partial x_k} \right)^* - \overline{{ \delta u^{\prime \prime }_j u^{''* }_k}} \delta \left( \frac{ \partial U_i }{ \partial x_k } \right) - \overline{{ \delta u''_i u^{\prime \prime *}_k }} \delta \left( \frac{ \partial U_j }{ \partial x_k } \right)}_{\text{Mean-stochastic production} ({\mathsf{p}}^{ms}_{ij})} \nonumber\\ &\quad +\underbrace{\left[ - \overline{{ \delta u''_j \delta u''_k}} \left( \frac{\partial \tilde{ {u} } _i}{\partial x_k} \right)^* - \overline{{ \delta u''_i \delta u''_k}} \left( \frac{\partial \tilde{ {u} } _j}{\partial x_k} \right)^* - \overline{{ \delta u''_j u^{\prime \prime *}_k }} \delta \left( \frac{ \partial \tilde{ {u} } _i }{\partial x_k} \right) - \overline{{ \delta u''_i u^{\prime \prime *}_k }} \delta \left( \frac{ \partial \tilde{ {u} } _j }{\partial x_k} \right) \right] }_{\text{Coherent-stochastic production} \ ({\mathsf{p}}^{cs}_{ij})} \nonumber\\ &\quad \underbrace{+ \frac{1}{\rho} \overline{{ \delta p'' \frac{\partial \delta u''_i}{\partial X_j}}} + \frac{1}{\rho} \overline{{ \delta p'' \frac{\partial \delta u''_j}{\partial X_i} }}}_{\text{Pressure strain} \ ({\rm \pi}^s_{ij})} \underbrace{- 4 \epsilon_{ij}^{s \ast}}_{\text{Dissipation} ({\mathsf{d}}^s_{ij})} + \underbrace{\overline{{\delta u''_j \delta f''_i}} + \overline{{\delta u''_i \delta f''_j}}}_{\text{Forcing interaction} ({\mathsf{f}}^s_{ij})}. \end{align}

Among the terms appearing in the source, the mean–coherent and mean–stochastic productions ${\mathsf{p}}^{mc}_{ij}$ and ${\mathsf{p}}^{ms}_{ij}$ indicate the scales and positions where the mean flow feeds, or drains energy from, the coherent and stochastic fields: they are not positive definite, and therefore can be either sources or sinks. They both contribute to the mean-fluctuating production ${\mathsf{P}}^{mf}_{ij}$ in (2.11), as ${\mathsf{P}}^{mf}_{ij}=\big\langle {{\mathsf{p}}^{mc}_{ij}} \big\rangle + \big\langle {{\mathsf{p}}^{ms}_{ij}} \big\rangle $. The coherent–stochastic production ${\mathsf{p}}^{cs}_{ij}$ indicates the exchange of stresses between the coherent and stochastic fields, and appears in the budgets for $\delta \tilde{u} _i \delta \tilde{u} _j$ and $\overline {{\delta u''_i \delta u''_j}}$ with opposite sign. Terms ${\mathsf{d}}^c_{ij}$ and ${\mathsf{d}}^s_{ij}$ denote viscous dissipation, and the pressure–strain terms ${\rm \pi} ^c_{ij}$ and ${\rm \pi} ^s_{ij }$ describe the interplay between pressure and velocity fields. Pressure–strain terms involve neither production nor dissipation of energy, and no cross-talk between coherent and fluctuating fields. Overall, among the source terms, the productions ${\mathsf{p}}^{mc}_{ij}$, ${\mathsf{p}}^{ms}_{ij}$ and ${\mathsf{p}}^{cs}_{ij}$ are the only ones that connect the mean, coherent and fluctuating budgets, and are essential to ascertain how the mean, stochastic and coherent fields force each other. The forcing interactions ${\mathsf{f}}^c_{ij}$ and ${\mathsf{f}}^s_{ij}$ represent the power injected into the system by the interaction of a coherent and stochastic external volume forcing with the coherent and stochastic flow fields, respectively. Finally, in (2.17) for $\delta \tilde{u} _i \delta \tilde{u} _j$ a new term $\zeta ^c_{ij}$ appears on the right-hand side. It describes the inter-phase interaction driven by the coherent flow field, and is defined as

(2.25)\begin{align} \zeta^c_{ij} & = \frac{\partial }{\partial r_k} \left[ \left\langle {\delta \tilde{ {u} } _i \delta \tilde{ {u} } _k} \right \rangle \delta \tilde{ {u} } _j + \left\langle {\delta \tilde{ {u} } _j \delta \tilde{ {u} } _k} \right \rangle \delta \tilde{ {u} } _i \right] + \frac{\partial}{\partial X_k} \left[ \left\langle { \tilde{ {u} } _k^* \delta \tilde{ {u} } _i } \right \rangle \delta \tilde{ {u} } _j + \left\langle { \tilde{ {u} } _k^* \delta \tilde{ {u} } _j } \right \rangle \delta \tilde{ {u} } _i \right] \nonumber\\ &\quad + \frac{\partial }{\partial r_k} \left[ \left\langle {\delta u''_i \delta u''_k} \right \rangle \delta \tilde{ {u} } _j + \big\langle {\delta u''_j \delta u''_k} \big\rangle \delta \tilde{ {u} } _i \right] + \frac{\partial}{\partial X_k} \left[ \left\langle {u''^*_k \delta u''_i } \right \rangle \delta \tilde{ {u} } _j + \big\langle {u''^*_k \delta u''_j } \big\rangle \delta \tilde{ {u} } _i \right] \nonumber\\ &\quad - \left\langle {\delta \tilde{ {u} } _i \delta \tilde{ {u} } _k } \right \rangle \left( \frac{\partial \tilde{ {u} } _j}{\partial x_k} \right)^* - \left\langle {\delta \tilde{ {u} } _j \delta \tilde{ {u} } _k } \right \rangle \left( \frac{\partial \tilde{ {u} } _i}{\partial x_k} \right)^* - \left\langle {\delta \tilde{ {u} } _i \tilde{ {u} } _k^* } \right \rangle \delta \left( \frac{\partial \tilde{ {u} } _j}{\partial x_k} \right) - \left\langle {\delta \tilde{ {u} } _j \tilde{ {u} } _k^* } \right \rangle \delta \left( \frac{\partial \tilde{ {u} } _i}{\partial x_k} \right)\nonumber\\ &\quad - \left\langle {\delta u''_i \delta u''_k } \right \rangle \left( \frac{\partial \tilde{ {u} } _j}{\partial x_k} \right)^* - \big\langle {\delta u''_j \delta u''_k } \big\rangle \left( \frac{\partial \tilde{ {u} } _i}{\partial x_k} \right)^* - \left\langle {\delta u''_i u''^*_k } \right \rangle \delta \left( \frac{\partial \tilde{ {u} } _j}{\partial x_k} \right) - \big\langle {\delta u''_j u''^*_k } \big\rangle \delta \left( \frac{\partial \tilde{ {u} } _i}{\partial x_k} \right). \end{align}

The terms in the last two rows above resemble a production term, and indicate the production of $\delta \tilde{u} _i \delta \tilde{u} _j$ due to the correlation of each phase with all the others.

By averaging equations (2.17) and (2.18) over the phases, the budget equations for $\left\langle {\delta \tilde{u} _i \delta \tilde{u} _j} \right \rangle(\boldsymbol {X},\boldsymbol {r})$ and $\big\langle {\delta u''_i \delta u''_j} \big\rangle(\boldsymbol {X},\boldsymbol {r})$ are obtained. In doing this, the inter-phase contributions vanish, since by definition they have zero average. The sum of the equations for the three diagonal components of $\left\langle {\delta \tilde{u} _i \delta \tilde{u} _j} \right \rangle$ and $\big\langle {\delta u''_i \delta u''_j} \big\rangle$ yields the GKE used by Alves Portela et al. (Reference Alves Portela, Papadakis and Vassilicos2020). If the equations for $\left\langle {\delta \tilde{u} _i \delta \tilde{u} _j} \right \rangle$ and $\big\langle {\delta u''_i \delta u''_j} \big\rangle$ are added together, the standard AGKE for the fluctuating field $\big\langle {\delta u'_i \delta u'_j} \big\rangle$ are recovered.

3.1. Database and computational details

The $\varphi$AGKE terms are computed from two datasets obtained by direct numerical simulations. They are described by Gallorini, Quadrio & Gatti (Reference Gallorini, Quadrio and Gatti2022), where the interested reader can find full details.

The simulations are run under a constant pressure gradient (Quadrio, Frohnapfel & Hasegawa Reference Quadrio, Frohnapfel and Hasegawa2016), with a friction Reynolds number of $Re_\tau = u_\tau h / \nu =200$, where $h$ is the channel half-height. A constant pressure gradient provides a unique value of $u_\tau$ with/without drag reduction, thus avoiding ambiguities in viscous scaling. The size of the computational domain is $(L_x,L_y,L_z) = (4{\rm \pi} h, 2h, 2{\rm \pi} h)$ in the streamwise, wall-normal and spanwise directions. The number of Fourier modes is $N_x= N_z= 256$ in the two homogeneous (streamwise and spanwise) directions, further increased by a factor of $3/2$ to remove aliasing error. In the wall-normal direction, a hyperbolic tangent distribution of $N_y = 192$ points provides a finer grid near the wall. The spatial resolution is $\Delta x^+ = 6.6$ and $\Delta z^+ = 3.3$ by considering the extra modes, while $\Delta y^+$ varies from $\Delta y^+ \approx 0.5$ close to the wall to $\Delta y^+ \approx 3.7$ at the centreline.

A first simulation of a plane channel with fixed walls is run as a reference, followed by two others in which wall oscillation according to (3.1) is enforced. The oscillation amplitude is fixed at $A^+=7$: a rather small value, which keeps the energy cost of the actuation limited, and might even provide a small net energy saving at optimal periods. As in Agostini & Leschziner (Reference Agostini and Leschziner2014), we consider two control periods, namely $T^+=75$ and $T^+=250$. The value $T^+=75$ is nearly optimal, and yields drag reduction (defined here as a percentage decrease of the friction coefficient, determined by the increase in bulk velocity) of 25.2 %. The value $T^+=250$ is suboptimal, and yields only 13.2 % drag reduction. These figures are in agreement with existing information (see e.g. Gatti & Quadrio Reference Gatti and Quadrio2016).

Simulations are started from an uncontrolled turbulent flow field. During the initial, transient phase, the solution is advanced by setting the Courant–Friedrichs–Lewy number at $1$. After the transient, however, the time step is set to a fixed value, in order to synchronise data saving with predetermined control phases. The value of the time step is thus chosen as an integer submultiple of the forcing period that keeps the maximum Courant–Friedrichs–Lewy number below unity: it is $\Delta t^+ = 0.0938$ for the smaller period and $\Delta t^+ = 0.0781$ for the longer period. After the transient, $376$ complete velocity fields are saved, so that $47$ control periods are stored for later analysis, each of them divided in 8 equally spaced phases.

The $\varphi$AGKE terms are computed from the database with a post-processing code derived with modifications from that described by Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020). It employs the same important numerical optimisations described in Gatti et al. (Reference Gatti, Remigi, Chiarini, Cimarelli and Quadrio2019), which include the computation of correlations pseudo-spectrally whenever possible. The code, written in the CPL computer programming language (Luchini Reference Luchini2020, Reference Luchini2021), has been validated by checking that the sum of each term of the budget of coherent and stochastic fields equals the corresponding term of $\big\langle {\delta u'_i \delta u'_j} \big\rangle$ within round-off. Statistical convergence of the results is verified by ensuring that the residuals of the budgets are negligible compared with the values of the production, pressure–strain and dissipation.

3.2. The $\varphi$AGKE tailored to channel flow with oscillating walls

The general form (2.17) and (2.18) of the $\varphi$AGKE can be simplified for the problem at hand. Since $x$ and $z$ are homogeneous, in an indefinite plane channel the $\varphi$AGKE depend on five independent variables: the three components of the separation vector $(r_x,r_y,r_z)$, the wall-normal component of the midpoint $Y$ and the phase $\varphi$. Note that the finite distance between the two walls implies the constraint $r_y < 2 Y$.

In an indefinite channel flow, the $x$ direction aligns with the mean flow, hence $\boldsymbol {U}(y)=(U(y),0,0)$, and the wall-parallel derivatives of the mean velocity are zero. Moreover, in the specific case of the oscillating wall, the coherent velocity field is independent on $x$ and $z$, as the wall control law (3.1) is a function of time only, so that $\partial \tilde{u} _i / \partial x = \partial \tilde{u} _i / \partial z = 0$. Therefore, incompressibility and no penetration at the wall dictate that the wall-normal component of the coherent field is null everywhere, i.e. $\tilde { {v} } (y,t)=0$. The streamwise coherent velocity $\tilde{u}$, instead, does not vanish, albeit it is known to be extremely small: (Yakeno et al. Reference Yakeno, Hasegawa and Kasagi2014) report it to be two orders of magnitude smaller than the spanwise coherent velocity $\tilde { {w} }$. The non-zero components of the $\delta \tilde{u} _i \delta \tilde{u} _j$ tensor are $\delta \tilde{u} \delta \tilde{u}$, $\delta \tilde { {w} } \delta \tilde { {w} }$ and $\delta \tilde{u} \delta \tilde { {w} }$.

The specialised form of the $\varphi$AGKE for the channel flow with oscillating walls is reported in Appendix B. It can be observed that the mean–coherent production ${\mathsf{p}}^{mc}_{ij}$ is zero: in this particular case, there is no exchange of stresses between the mean and coherent fields, as the coherent field interacts directly with the external forcing and with the stochastic field only. However, this term does appear in other flows, and for example is important for the flow past a bluff body (Alves Portela et al. Reference Alves Portela, Papadakis and Vassilicos2020), where the mean flow supports the coherent vortex shedding, which in turn supports the stochastic fluctuations. In the budget for the stochastic part, the productions ${\mathsf{p}}^{ms}_{ij}$ and ${\mathsf{p}}^{cs}_{ij}$ represent the two avenues for the stochastic field to interact with the mean and coherent fields, involving distinct components of $\overline {{\delta u''_i \delta u''_j}}$. The mean–stochastic production ${\mathsf{p}}^{ms}_{ij}$ is non-zero only for $\overline {{\delta u'' \delta u''}}$ and for the off-diagonal components $\overline {{\delta u'' \delta v''}}$ and $\overline {{\delta u'' \delta w''}}$. In contrast, the coherent–stochastic production contributes to all the elements of $\overline {{\delta u''_i \delta u''_j}}$ except for $\overline {{\delta v'' \delta v''}}$, being ${\mathsf{p}}^{cs}_{22}=0$.

The flow symmetries and the type of forcing make only certain paths available for energy exchanges. This is represented graphically in figure 3, which shows an ‘energy circle’ (Quadrio Reference Quadrio2011) to describe energy exchanges among the mean, coherent and stochastic fields after spatial and temporal integration. In the following, thanks to the $\varphi$AGKE, these global energy exchanges and redistributions are expanded and described in space and among scales, with a phase-by-phase analysis.

Figure 3. Sketch of the energy exchanges between mean, coherent and stochastic fields for the turbulent channel flow modified by spanwise-oscillating walls. Blue (red) arrows indicate energy entering (leaving) the system. The blue arrows $P_m$ and $\tilde { {w} }$ represent the pumping energy required to move the flow, and the energy introduced by the moving walls.

4.1. Near-wall structures

4.1.1. Description at a fixed phase

Figure 4 shows the diagonal components of $\overline {{\delta u_i'' \delta u_j''}}$ in the $r_y=r_x=0$ plane for the uncontrolled channel (figure 4a–c), $T^+=75$ (figure 4d–f) and $T^+=250$ (figure 4g–i). For the two controlled cases, only phase $\varphi _4$ is shown, but the discussion that follows is qualitatively valid for all phases.

Figure 4. Diagonal components of the stochastic tensor $\overline {{\delta u_i'' \delta u_j''}}^+$ at $\varphi _4$ in the ($r_z^+, Y^+$) plane. (a–c) Uncontrolled case with $A=0$, (d–f) $T^+=75$ and (g–i) $T^+=250$. The contour is set at 95 % of each maximum. The coordinates of the maximum, marked with a cross, can be read on the axes.

The local maxima of $\overline {{\delta u'' \delta u''}}$ and $\overline {{\delta v'' \delta v''}}$, hereafter denoted with subscript $m$, are the statistical trace of the structures of the near-wall cycle. In the $r_x=r_y=0$ space, indeed, they indicate a negative peak of the streamwise and vertical stochastic correlation functions ${\mathsf{R}}_{11}$ and ${\mathsf{R}}_{22}$; see (2.7). The coordinates $Y^+ \approx 14\unicode{x2013}18$ and $r_z^+ \approx 55\unicode{x2013}65$ of $\overline {{\delta u'' \delta u''}}_m$ in the $(r_z^+,Y^+)$ plane indicate the characteristic wall distance and spanwise spacing of low- and high-speed streaks. The coordinates $Y^+ \approx 43\unicode{x2013}55$ and $r_z^+ \approx 49\unicode{x2013}59$ of $\overline {{\delta v'' \delta v''}}_m$ indicate the characteristic wall distance and spanwise size of the quasi-streamwise vortices, which induce at their spanwise sides regions of vertical velocity with negative correlation.

Figure 4 shows that the oscillating wall leaves $\overline {{\delta u'' \delta u''}}$ and $\overline {{\delta v'' \delta v''}}$ almost unchanged, indicating that the size and strength of the near-wall structures only marginally depend on the amount of drag reduction.

This is consistent with the constant pressure gradient driving strategy, which forces the same level of wall friction; the large changes observed by various authors under different driving strategies simply derive trivially from the different friction, as discussed by Frohnapfel, Hasegawa & Quadrio (Reference Frohnapfel, Hasegawa and Quadrio2012). However, the velocity streaks are slightly moved away from the wall: an upward shift of $\overline {{\delta u'' \delta u''}}_m$ can be seen in figure 4. The previous observation is confirmed by numerical data: the maximum moves from ${Y^+=14.1}$ in the reference case to $Y^+=17.8$ for $T^+=75$ and to ${Y^+=14.7}$ for $T^+=250$ (at phase $\varphi _4$). Both shifts are upwards, and the $T^+=75$ case with larger drag reduction has a larger shift. The quasi-streamwise vortices react differently from control: $\overline {{\delta v'' \delta v''}}_m$ moves from $Y^+=53$ in the reference case to $Y^+=55$ for $T^+=75$ and to $Y^+=43$ for $T^+=250$. These contrasting trends are consistent with the wall-normal displacement found by Gallorini et al. (Reference Gallorini, Quadrio and Gatti2022) for conditionally averaged quasi-streamwise vortices, but are extracted from the present analysis without the need for an (inevitably subjective) procedure for conditional structure extraction.

In the canonical channel flow, the map of $\overline {{\delta w'' \delta w''}}$ embeds information of the quasi-streamwise vortices only when the $r_y \neq 0$ space is considered, which contains the peak $\left\langle {\delta w' \delta w'} \right \rangle_m$ (Gatti et al. Reference Gatti, Chiarini, Cimarelli and Quadrio2020). Indeed, the quasi-streamwise vortices induce negatively correlated regions of $w''$ fluctuations at their vertical sides only, and the $r_y$ coordinate of the maximum indicates their characteristic wall-normal size. In the controlled cases, however, a local peak of $\overline {{\delta w'' \delta w''}}$ appears in the $r_x=r_y=0$ (figure 4) and $r_z=r_y=0$ (not shown) planes. Interestingly, the local peak is particularly evident for $T^+=250$, extending for $r_z^+ \approx 50\unicode{x2013}100$, $r_x^+ \approx 85\unicode{x2013}270$ and $Y^+ \approx 13\unicode{x2013}25$, but it is hardly visible for $T^+=75$, where the $w''$ fluctuations are weaker. The next subsection, which examines how these quantities vary with $\varphi$, shows that this derives from a combination of the streaks tilting in the $x\unicode{x2013}z$ plane and from the interaction of the quasi-streamwise vortices with the coherent spanwise shear.

4.1.2. Evolution during the cycle

Figure 5 shows the phase evolution of $\overline {{\delta u'' \delta u''}}$, $\overline {{\delta v'' \delta v''}}$ and $\overline {{\delta w'' \delta w''}}$ in the $r_x=r_y=0$ plane, to describe how the organisation of the near-wall stochastic fluctuations changes during the oscillation cycle, i.e. the very type of information that the $\varphi$AGKE are designed to provide. Only the suboptimal $T^+=250$ is considered, as the large period emphasises the phase dependence; moreover, only one half of the forcing period is shown (from $\varphi _1$ to $\varphi _4$), because of temporal symmetry. Extra quantitative information is provided by figure 6, which plots the phase evolution of the maxima $\overline {{\delta u'' \delta u''}}_m$, $\overline {{\delta v'' \delta v''}}_m$ and $\overline {{\delta w'' \delta w''}}_m$.

Figure 5. Diagonal components of the stochastic tensor $\overline {{\delta u_i'' \delta u_j''}}^+$ in the $(r_z^+,Y^+)$ plane, at phases (a–c) $\varphi _1$, (d–f) $\varphi _2$, (g–i) $\varphi _3$ and (j–l) $\varphi _4$, for the period $T^+=250$. For each component, the white contour is set at the 95 % of the smallest peak over the phases (i.e. at $\varphi _2$ for $\overline {{\delta u'' \delta u''}}$ and $\overline {{\delta w'' \delta w''}}$, and at $\varphi _3$ for $\overline {{\delta v'' \delta v''}}$).

Figure 6. (a–c) Phase variation of the maxima $\overline {{\delta u_i'' \delta u_i''}}_m^+$ in the $(r_z^+,Y^+)$ plane.

The streamwise velocity streaks cyclically strengthen and weaken under the action of the alternating Stokes layer. The maximum $\overline {{\delta u'' \delta u''}}_m$ assumes its lowest value at $\varphi _2$, and then grows to reach the highest value at $\varphi _4$, with an intra-cycle variation of 27 %. The quasi-streamwise vortices, instead, show a much smaller phase dependence: the intra-cycle variation of $\overline {{\delta v'' \delta v''}}$ is 8 % only. This is not surprising, since the quasi-streamwise vortices reside at larger wall distances, where the intensity of the Stokes layer is lower; at ${y^+ = 14}$, the average position of the streaks, the maximum $\tilde { {w} } ^+$ is $1.15$, while at $y^+=50$, representative wall-normal distance of the vortices, it is only $0.2$. A different wall distance for streaks and vortices also implies a phase shift; in fact the intensity of $\overline {{\delta v'' \delta v''}}$ is minimum at $\varphi _3$ and maximum at $\varphi _1$, whereas $\overline {{\delta u'' \delta u''}}$ and $\overline {{\delta w'' \delta w''}}$ are minimum at $\varphi _2$ and maximum at $\varphi _4$. This is consistent with the early observation (Baron & Quadrio Reference Baron and Quadrio1996) that streaks and quasi-streamwise vortices are displaced by the spanwise Stokes layer differently.

From figure 6, one notices that the phase evolution of $\overline {{\delta w'' \delta w''}}_m$ resembles that of $\overline {{\delta u'' \delta u''}}_m$, thus suggesting that part of the stochastic $w''$ fluctuations derives from a redistribution of the streamwise fluctuations. The near-wall structures are tilted in the $x\unicode{x2013}z$ plane and follow the shear vector $(\text {d}U/\text {d}y,0,\partial \tilde { {w} } /\partial y)$ (Yakeno et al. Reference Yakeno, Hasegawa and Kasagi2014). The tilting causes the streamwise high- and low-speed streaks to reorient, thus contributing via pressure–strain redistribution (see below § 4.3) to the spanwise stochastic fluctuations. When the tilting angle is positive (negative), the low- and high-speed streaks contribute to respectively positive (negative) and negative (positive) $w''$. This produces regions of $w''$ fluctuations that correlate negatively for scales $r_x$ and $r_z$ and position $Y$ compatible with the position of $\overline {{\delta w'' \delta w''}}_m$ observed in figure 5. This is shown with a sketch in figure 7, and confirmed with a phase-by-phase conditional average of events extracted from the present database in Appendix C. The picture is also consistent with the lower $\overline {{\delta w'' \delta w''}}_m$ observed in figure 4 for $T^+=75$: for periods close to the optimum, the oscillation is too fast for the streaks to align with the shear vector (Touber & Leschziner Reference Touber and Leschziner2012), and this redistribution mechanism becomes weaker.

Figure 7. Sketch of the contribution of $u''$ and $w''$ for positively (a) and negatively (b) tilted low-speed (blue) and high-speed (red) streaks induced by a positively rotating quasi-streamwise vortex (white).

Similar information is usually extracted (Yakeno et al. Reference Yakeno, Hasegawa and Kasagi2014) from phase-locked conditional averages. However, such statistics are unavoidably arbitrary to some degree; e.g. ‘short’ structures have to be excluded from averaging, and one needs to predetermine a specific wall distance for the eduction procedure. Here we obtain information that is equivalent to conditional averaging, but via a statistical analysis that is free from assumptions and hypotheses.

For example, the scales $r_{z,m}$ and $r_{x,m}$ identified by $\overline {{\delta w'' \delta w''}}_m$ can be used to track the phase evolution of the tilting angle $\theta$ of the flow structures during the cycle:

(4.1)\begin{equation} |\theta(\varphi)| = \tan^{-1} \left( \frac{r_{z,m}(\varphi)}{r_{x,m}(\varphi)} \right) . \end{equation}

Similarly, the evolution of the wall-normal position $Y_m$ of $\overline {{\delta w'' \delta w''}}_m$ (or, equivalently, of $\overline {{\delta u'' \delta u''}}_m$) quantifies the vertical displacement of the streaks during the cycle.

Figure 8 compares $|\theta |$ with the shear angle $\theta _s$ evaluated at $Y_m$, i.e.

(4.2)\begin{equation} \theta_s = \tan^{-1} \left( \frac{ \partial \tilde{ {w} } /\partial y}{\text{d}U /\text{d}y} \right), \end{equation}

that is conventionally used to describe the tilting angle of the near-wall structures (Yakeno et al. Reference Yakeno, Hasegawa and Kasagi2014; Gallorini et al. Reference Gallorini, Quadrio and Gatti2022). The two quantities $\theta$ and $\theta _s$ are quantitatively similar and present the same phase dependence, with a nearly constant difference of approximately $8^\circ$. Figure 8(b) also shows that when the tilting angle of the streaks is maximum, their distance from the wall is minimum (and vice versa). This implies that a higher coherent spanwise velocity yields a larger tilting.

Figure 8. (a) Evolution of the tilt angle of the wall streaks during the cycle. Comparison between present results (blue symbols) and the shear angle introduced by Yakeno et al. (Reference Yakeno, Hasegawa and Kasagi2014) (red symbols). (b) Wall-normal position of the structures, educed from the wall-normal position $Y_m^+$ of $\overline {{\delta w'' \delta w''}}_m$.

Part of the wall-parallel modulation of $\overline {{\delta w'' \delta w''}}$ induced by the wall oscillation derives from the interaction of the quasi-streamwise vortices with the coherent spanwise shear. When the coherent shear $\partial \tilde { {w} } /\partial y$ is positive, the quasi-streamwise vortices move low-spanwise-velocity fluid upwards, and high-spanwise-velocity fluid downwards. The opposite happens when $\partial \tilde { {w} } /\partial y<0$. This creates two regions with spanwise velocity of opposite sign at the vortex sides, resulting in negative ${\mathsf{R}}_{33}$ correlation and a positive peak of $\overline {{\delta w'' \delta w''}}$ at their characteristic spanwise separation. This process, quantified by the coherent–stochastic production ${\mathsf{p}}^{cs}_{33}$ (see § 4.2), resembles the ejections and sweeps typical of the near-wall cycle, where the mean streamwise shear is involved; its description is similar to the explanation provided by Agostini et al. (Reference Agostini, Touber and Leschziner2014) for the non-zero $\big\langle {v''w''} \big\rangle $. Once again, our interpretation is supported by the velocity field induced by the ensemble-averaged quasi-streamwise vortex, computed at various phases and shown in Appendix C.

4.2. Interaction of the mean, coherent and fluctuating fields

The energy exchanges of the mean field with the stochastic and coherent fields are described by the two mean production terms ${\mathsf{p}}^{mc}_{ij}$ and ${\mathsf{p}}^{ms}_{ij}$. However, as shown in figure 3, for the present problem ${\mathsf{p}}^{mc}_{ij}=0$, and the mean field interacts directly with the stochastic field only, by feeding (or draining from) streamwise fluctuations. Moreover, energy is exchanged between the coherent and stochastic fields via the coherent–stochastic production ${\mathsf{p}}^{cs}_{ij}$, which involves only $\overline {{\delta u'' \delta u''}}$ and $\overline {{\delta w'' \delta w''}}$ among the diagonal components of the $\overline {{\delta u''_i \delta u''_j}}$ tensor.

Figure 9 shows how the mean–stochastic production ${\mathsf{p}}^{ms}_{11}$ varies with $\varphi$ for $T^+=75$ (figure 9a,c,e,g) and $T^+=250$ (figure 9b,d,f,h) in the $r_x=r_y=0$ plane, where the production terms are maxima. Here, ${\mathsf{p}}^{ms}_{11}$ reduces to

(4.3)\begin{equation} {\mathsf{p}}^{ms}_{11} =- 2 \overline{{\delta u'' \delta v''}} \left( \frac{\text{d}U}{\text{d}y} \right). \end{equation}

The mean–stochastic production is positive everywhere, with a peak in the range $r_{z,m}^+ = 36\unicode{x2013}42$ and $Y_m^+= 13\unicode{x2013}17$ for $T^+=75$ and $r_{z,m}^+ =36\unicode{x2013}39$ and $Y_m^+ =12\unicode{x2013}14$ for $T^+=250$. Hence, the interaction of the near-wall cycle ($\overline {{\delta u'' \delta v''}}$) with the mean shear ($\text {d}U/\text {d}y$) invariably moves energy from the mean field towards the stochastic streamwise fluctuations. Note that the smaller $Y^+$ for $T^+=250$ is consistent with the reduced thickening of the viscous sublayer for suboptimal periods. The production intensity is largest at $\varphi _1$ and lowest at $\varphi _3$ for $T^+=75$, whereas it is largest at $\varphi _3$ and lowest at $\varphi _1$ for $T^+=250$. Since $\text {d}U/\text {d}y$ is phase-independent, this can only descend from $\overline {{\delta u'' \delta v''}}$, which includes the phase evolution of the streaks and of the quasi-streamwise vortices (see § 4.1.2).

Figure 9. Mean–stochastic production ${{\mathsf{p}}^{ms}_{11}}^+$ in the $(r_z^+,Y^+)$ plane for $T^+=75$ (a,c,e,g) and $T^+=250$ (b,d,f,h), at phases (a,b) $\varphi _1$, (c,d) $\varphi _2$, (e,f) $\varphi _3$ and (g,h) $\varphi _4$. The contour line is set at 95 % of the smallest maximum over the phases. The coordinates of the maximum, marked with a cross, can be read on the axes.

Figures 10 for $T^+=75$ and 11 for $T^+=250$ show how ${\mathsf{p}}^{cs}_{11}$ and ${\mathsf{p}}^{cs}_{33}$ change with $\varphi$. Like for ${\mathsf{p}}^{ms}_{11}$, the expressions for ${\mathsf{p}}^{cs}_{11}$ and ${\mathsf{p}}^{cs}_{33}$ simplify in the $r_x=r_y=0$ plane where their maxima occur, i.e.

(4.4a,b)\begin{equation} {\mathsf{p}}^{cs}_{11} =- 2 \overline{{\delta u'' \delta v''}} \left( \frac{ \partial \tilde{ {u} } }{\partial y} \right) \quad \text{and} \quad {\mathsf{p}}^{cs}_{33} =- 2 \overline{{\delta v'' \delta w''}} \left( \frac{ \partial \tilde{ {w} } }{\partial y} \right). \end{equation}

Unlike ${\mathsf{p}}^{ms}_{11}$, however, these productions can take either sign. Their maps show evident horizontal stripes of alternating sign, from the wall up to $Y^+ \approx 40$: hence, at a given phase the coherent field feeds the stochastic field at certain wall distances, but extracts energy from it at others. It is worth noting that, although ${\mathsf{p}}^{cs}_{11}$ and ${\mathsf{p}}^{cs}_{33}$ at a given phase are both positive and negative, after averaging over the phases $\left\langle {{\mathsf{p}}^{cs}_{11}} \right \rangle$ almost vanishes and $\left\langle {{\mathsf{p}}^{cs}_{33}} \right \rangle$ is positive everywhere. This is not entirely new, and confirms the single-point analysis by Agostini et al. (Reference Agostini, Touber and Leschziner2014) (see their figure 14); however, scale information is added here so that this mechanism can be related to the structures of the flow. At every phase, the positive/negative peaks of ${\mathsf{p}}^{cs}_{11}$ and ${\mathsf{p}}^{cs}_{33}$ occur at $r_z^+ \approx 25\unicode{x2013}50$, a spanwise separation which points to the structures of the near-wall cycle.

Figure 10. Coherent–stochastic production ${{\mathsf{p}}^{cs}_{11}}^+$ (a,c,e,g) and ${{\mathsf{p}}^{cs}_{33}}^+$ (b,d,f,h) in the $(r_z^+,Y^+)$ plane for $T^+=75$, at phases (a,b) $\varphi _1$, (c,d) $\varphi _2$, (e,f) $\varphi _3$ and (g,h) $\varphi _4$. The thin contour line is set at 95 % of the smallest (positive and negative) maximum over the phases; the thick black contour line is ${\mathsf{p}}^{cs}_{ii}=0$. The coordinates of the maximum, marked with a cross, can be read on the axes.

Figure 11. As in figure 10, but for $T^+=250$.

The intensity of ${\mathsf{p}}^{cs}_{11}$ and ${\mathsf{p}}^{cs}_{33}$ at the two periods is comparable, at all scales and positions. However, for ${\mathsf{p}}^{cs}_{11}$ the contribution of the shear stresses is dominant, whereas the opposite occurs for ${\mathsf{p}}^{cs}_{33}$, where the coherent spanwise shear dominates. Indeed, $\partial \tilde { {w} } /\partial y$ is two orders of magnitude larger than $\partial \tilde{u} /\partial y$, and $\overline {{\delta v'' \delta w''}}$ is two orders of magnitude smaller than $\overline {{\delta u'' \delta v''}}$. Note, moreover, that for both control periods ${\mathsf{p}}^{ms}_{11} \gg {\mathsf{p}}^{cs}_{11}$, meaning that the streamwise stochastic fluctuations are predominantly fed by the mean field.

The alternating positive/negative stripes for ${\mathsf{p}}^{cs}_{11}$ and ${\mathsf{p}}^{cs}_{33}$ are due to the change of sign of $\partial \tilde{u} /\partial y$ and $\partial \tilde { {w} } /\partial y$ with $y$. For ${\mathsf{p}}^{cs}_{33}$, the changing sign of the shear is also indirectly responsible for the alternating positive/negative $\overline {{\delta v'' \delta w''}}$, due to the quasi-streamwise vortex–shear interaction described in § 4.1.2. In contrast, for ${\mathsf{p}}^{cs}_{11}$, $\overline {{\delta u'' \delta v''}}$ is entirely due to the interaction of the near-wall structures with the mean shear $\text {d}U/\text {d}y$, which overwhelms $\partial \tilde{u} /\partial y$ everywhere.

Comparing figures 10 and 11 highlights that the slower oscillation introduces substantial differences in the coherent–stochastic energy exchange. The positive/negative maxima of ${\mathsf{p}}^{cs}_{11}$ increase, and their positions move towards larger $r_z$ and larger $Y$, but the effect of $T^+$ on ${\mathsf{p}}^{cs}_{33}$ is even more evident. At $T^+=250$, the stripes of negative ${\mathsf{p}}^{cs}_{33}$ weaken, while those with ${\mathsf{p}}^{cs}_{33}>0$ strengthen; overall, the spanwise contribution to the energy flowing from the coherent to the stochastic field becomes larger. A larger oscillating period implies a larger thickness of the Stokes layer, proportional to $\sqrt {\nu T}$, thus stretching outwards the coherent spanwise shear and, as a consequence, the scale-space map of $\overline {{\delta v'' \delta w''}}$, yielding an overall increase of the positive ${\mathsf{p}}^{cs}_{33}$. At $\varphi _2$ and $\varphi _3$, for example, $\partial \tilde { {w} } /\partial y$ is negative close to the wall and changes sign only at $y^+ \approx 30\unicode{x2013}50$ for $T^+=250$ (see figure 2), while it changes sign already at $y^+ \approx 13\unicode{x2013}18$ for $T^+=75$ (not shown). For $T^+=250$ this results in a large increase of the near-wall positive ${\mathsf{p}}^{cs}_{33}$, as highlighted by the dark red colour in figure 11. Due to the negative $\partial \tilde { {w} } /\partial y$, indeed, the quasi-streamwise vortices induce on their sides positive/negative $v''$ and convect upwards/downwards high/low spanwise velocity $w''$, thus yielding positive $\overline {{\delta v'' \delta w''}}$ and an intense energy exchange from the coherent to the stochastic field. The scale-space information of this exchange mechanism is highlighted by the positive peak of ${\mathsf{p}}^{cs}_{33}$ placed at $(r_z^+,Y^+) \approx (38,9)$ for the considered $\varphi _2$ and $\varphi _3$ phases.

4.3. Pressure–strain redistribution

As seen schematically in figure 3, the pressure–strain action partially redistributes the streamwise energy $\overline {{\delta u'' \delta u''}}$ drained from the mean flow towards the cross-stream fluctuations $\overline {{\delta v'' \delta v''}}$ and $\overline {{\delta w'' \delta w''}}$. Figure 12(a,c,e) shows that ${\rm \pi} _{11}^s<0$, ${\rm \pi} _{22}^s>0$ and ${\rm \pi} _{33}^s>0$ at almost all scales and positions: only in a very thin region close to the wall are ${\rm \pi} _{11}^s>0$, ${\rm \pi} _{22}^s<0$ and ${\rm \pi} _{33}^s>0$, according to the reorientation of vertical fluctuations into wall-parallel ones because of the impermeable wall (Mansour, Kim & Moin Reference Mansour, Kim and Moin1988). The peaks of ${\rm \pi} _{11}^s$, ${\rm \pi} _{22}^s$ and ${\rm \pi} _{33}^s$ in the ($r_z,Y$) plane have $Y_m^+ \approx 11\unicode{x2013}27$ and $r_{z,m}^+ \approx 30\unicode{x2013}52$, indicating that the energy redistribution is dominated by the near-wall cycle.

Figure 12. (a,c,e) Pressure–strain redistribution ${\rm \pi} _{ii}^{s+}$ from $\overline {{\delta u'' \delta u''}}$ towards $\overline {{\delta v'' \delta v''}}$ and $\overline {{\delta w'' \delta w''}}$ at phase $\varphi _4$ for $T^+=250$; the coordinates of the maximum, marked with a cross, can be read on the axes. (b,d,f) Phase variation of their maxima in the $(r_z^+,Y^+)$ plane, with a horizontal solid line indicating the value of the uncontrolled flow.

It is known (Touber & Leschziner Reference Touber and Leschziner2012; Yakeno et al. Reference Yakeno, Hasegawa and Kasagi2014) that the spanwise oscillation of the wall enhances the energy redistribution, mainly towards spanwise fluctuations. Compared with the uncontrolled case, the negative peak of ${\rm \pi} _{11}^s$ increases by 23 %–67 % for $T^+=75$ and by 36 %–77 % for $T^+=250$, while the positive peak of ${\rm \pi} _{22}^s$ decreases by 2 %–11 % for $T^+=75$ and increases by 4 %–29 % for $T^+=250$. The positive peak of ${\rm \pi} _{33}^s$, instead, has the largest variation, with an increase of 30 %–53 % for $T^+=75$ and 40 %–87 % for $T^+=250$.

The phase evolution of the pressure-mediated energy redistribution is described in figure 12(b,d,f) for the $T^+=250$ case, by considering the maxima of the diagonal components of ${\rm \pi} _{ij}^s$. Only their values are plotted, since their position remains nearly constant at $(Y^+,r_z^+) \approx (20,52)$ for ${\rm \pi} ^s_{11,m}$, $\approx (27,30)$ for ${\rm \pi} ^s_{22,m}$ and $\approx (12,46)$ for ${\rm \pi} ^s_{33,m}$. Like $\overline {{\delta v'' \delta v''}}_m$, ${\rm \pi} _{22,m}^s$ is the component with the smallest intra-cycle variation, with a 21 % excursion during the cycle compared with 30 % and 35 % for ${\rm \pi} _{11,m}^s$ and ${\rm \pi} _{33,m}^s$. In fact, the largest energy redistribution towards $\overline {{\delta v'' \delta v''}}$ occurs quite far from the wall, where the influence of the Stokes layer is weak. The phase dependence of ${\rm \pi} _{11,m}^s$ is qualitatively different from that of the others. The redistribution of $\overline {{\delta u'' \delta u''}}$ towards the cross-stream components is maximum at $\varphi _3$ and minimum at $\varphi _1$, following the absolute value of ${\rm \pi} _{11,m}^s$. In contrast, ${\rm \pi} _{22,m}^s$ and ${\rm \pi} _{33,m}^s$ are minima at $\varphi _2$ and maxima at $\varphi _4$ (this is not inconsistent with the incompressibility constraint ${\rm \pi} _{11}^s+ {\rm \pi}_{22}^s + {\rm \pi}_{33}^s = 0$, since the three maxima occur at different scales and positions). As already mentioned in § 4.1.2, ${\rm \pi} _{33,m}^s$ and $\overline {{\delta w'' \delta w''}}_m$ have the same phase dependence, confirming that the tilting of the near-wall structures is accompanied by a redistribution of the streamwise fluctuations towards the spanwise ones.

4.4. Transfers of the spanwise stresses

A peculiarity of the present flow is the direct connection between the Stokes layer and the stochastic stresses, described by the coherent–stochastic production ${\mathsf{P}}_{11}^{cs}$ and ${\mathsf{P}}_{33}^{cs}$ shown in figure 3. It is therefore interesting to examine the variable-phase scale-space transfers of the stochastic stresses by looking at their fluxes in the scale and physical spaces. In this analysis, we only consider the transfer of spanwise stresses $\overline {{\delta w'' \delta w''}}$, since for the streamwise stresses ${\mathsf{p}}^{cs}_{11}$ is negligible compared with ${\mathsf{p}}^{ms}_{11}$. Moreover, only the $T^+=250$ case is considered, as the one where the effect of the Stokes layer on the $w''$ field is larger. For simplicity, the analysis is restricted to the $r_x=r_y=0$ subspace, where the budget of $\overline {{\delta w'' \delta w''}}$ can be rewritten by moving to the right-hand side the off-plane flux divergence terms $\partial \phi _{x,33}^s/\partial r_x$, $\partial \phi _{y,33}^s/\partial r_y$ and the phase evolution term, as follows:

(4.5)\begin{equation} \frac{ \partial \phi_{z,33}^s }{ \partial r_z } + \frac{ \partial \psi_{33}^s }{ \partial Y } = \underbrace{ {\mathsf{p}}^{cs}_{33} + {\rm \pi}_{33}^s + {\mathsf{d}}_{33}^s}_{\xi^s_{33}} - \frac{ \partial \phi_{x,33}^s }{ \partial r_x } - \frac{ \partial \phi_{y,33}^s }{ \partial r_y } - \omega \frac{\partial \overline{{\delta w'' \delta w''}}}{\partial \varphi} . \end{equation}

In this way, the left-hand side features the divergence of the in-plane flux vector, which provides information on the energetic relevance of the fluxes with its intensity, and shows their direction via its field lines. Moreover, the off-plane fluxes (i.e. the last three terms in the equation above) are always very small, and the in-plane divergence approximates well the full source term $\xi _{33}^s$ everywhere (Gatti et al. Reference Gatti, Chiarini, Cimarelli and Quadrio2020). This descends from a combination of the symmetries of the plane channel flow system and of the approximate alignment of the dominant vortical structures with the streamwise direction. Hence, the scale-space properties of the source term $\xi _{33}^s$ approximate well those of the divergence of the in-plane flux.

Figure 13 plots the map of $\xi _{33}^s= {\mathsf{p}}^{cs}_{33} +{\rm \pi} _{33}^s + {\mathsf{d}}_{33}^s$ for the uncontrolled case (where ${\mathsf{p}}^{cs}_{33}=0$) and the controlled case at $T^+=250$ for $\varphi _1$, $\varphi _2$, $\varphi _3$ and $\varphi _4$, with the field lines of the in-plane flux coloured with its divergence. In the uncontrolled case, a region with $\xi ^s_{33}>0$ extends for $5 \lessapprox Y^+ \lessapprox 100$ and for $r_z^+ \gtrapprox 15$, at scales and positions where the pressure–strain dominates over dissipation. When control is active, instead, $\xi _{33}^s$ receives the additional contribution from coherent–stochastic production, and the values of $\xi ^s_{33}$ are generally larger. Two regions with $\xi ^s_{33}>0$ exist. One is close to the wall at $Y^+ \approx 10\unicode{x2013}20$, and extends for all scales $r_z^+ \gtrapprox 15$, with a peak at $r_z^+ \approx 40$. A second, connected region involves larger wall distances and scales, in the $40 \lessapprox r_z^+ \lessapprox 200$ range. It is clearly visible in figure 14, where the ratio ${\rm \pi} _{33}^s/( \xi _{33}^s - {\mathsf{d}}_{33}^s )$ is plotted to determine the main contribution to these positive sources at the different phases. When ${\rm \pi} _{33}^s/( \xi _{33}^s - {\mathsf{d}}_{33}^s )> 0.5$, ${\rm \pi} _{33}^s> {\mathsf{p}}^{cs}_{33}$ meaning that the pressure–strain is the largest contribution to the positive source. When ${\rm \pi} _{33}^s/( \xi _{33}^s - {\mathsf{d}}_{33}^s )< 0.5$, instead, the main contributor is the coherent–stochastic production ${\mathsf{p}}^{cs}_{33}$. Figure 14 shows that ${\mathsf{p}}^{cs}_{33}$ and ${\rm \pi} ^s_{33}$ both contribute to the near-wall source, but their relative importance changes with the phase. For $\varphi _2$ and $\varphi _3$, ${\mathsf{p}}^{cs}_{33}$ is the main contributor to the intense source peak. For $\varphi _1$ and $\varphi _4$, instead, ${\mathsf{p}}^{cs}_{33}$ weakens (see figure 11): now the (weaker) source is mainly fed by the pressure–strain. The source at larger $Y$, instead, is dominated by the pressure–strain at all phases; this is reasonable, as for $y^+ > 30$ the Stokes layer and consequently the coherent–stochastic production are weak.

Figure 13. Source $\xi _{33}^{s+}$ in the $(r_z^+,Y^+)$ plane, with field lines of the in-plane flux vector coloured with its divergence for (a) the uncontrolled case and for the controlled case at $T^+=250$ at phases (b) $\varphi _1$, (c) $\varphi _2$, (d) $\varphi _3$ and (e) $\varphi _4$. The thin contour line marks the zero level. Dots (coloured according to figure 2) indicate the singularity point for the near-wall source, and the black vertical line marks the cut-off spanwise scale $r_{z,min}^+$ (see text).

Figure 14. Region with positive source in the $(r_z^+,Y^+)$ plane at phases (a) $\varphi _1$, (b) $\varphi _2$, (c) $\varphi _3$ and (d) $\varphi _4$ for $T^+=250$. The colour scale is for the ratio ${\rm \pi} _{33}^s/( \xi _{33}^s - {\mathsf{d}}_{33}^s )$ and is centred at 0.5: red means ${\rm \pi} _{33}^s>{\mathsf{p}}^{cs}_{33}$ and blue means ${\rm \pi} _{33}^s<{\mathsf{p}}^{cs}_{33}$.

As for the sinks, figure 13 shows three of them: viscous effects dominate the very near-wall region ($Y \rightarrow 0$), the bulk flow ($Y \rightarrow h$) and the smallest scales ($r_z \rightarrow 0$). Extension and intensity of these sinks change with $\varphi$, according to the evolution of ${\mathsf{p}}^{cs}_{33}$, ${\rm \pi} _{33}^s$ and ${\mathsf{d}}_{33}^s$. A cut-off scale $r_{z,min}$ (Chiarini et al. Reference Chiarini, Gatti, Cimarelli and Quadrio2022a) can also be plotted to quantify the minimal scale where (spanwise) energy is always dissipated, regardless of the wall distance.

The field lines of $\overline {{\delta w'' \delta w''}}$ drawn in figure 13 originate from a singularity point, i.e. a point near the source peak where the direction of the fluxes is undefined. Here the lines are energised by the intense positive source and transfer $\overline {{\delta w'' \delta w''}}$ towards the sinks. Three types of lines are recognised, depending on where they vanish, and reflect the three sinks described above. Overall, these fluxes indicate the coexistence of ascending/descending and direct/inverse energy transfers, as described by Cimarelli et al. (Reference Cimarelli, De Angelis and Casciola2013, Reference Cimarelli, De Angelis, Jimenez and Casciola2016) and Chiarini et al. (Reference Chiarini, Mauriello, Gatti and Quadrio2022b) in the context of Poiseuille and Couette turbulent flows.

The three line types possess the same topology in the uncontrolled and controlled cases. For the latter, though, the amount of spanwise energy withdrawn from the sources and released to the sinks changes with $\varphi$. An estimate of this change is provided by the phase evolution of the positive peak of the two-dimensional divergence of the flux vector. Its value is maximum at $\varphi _3$ where it is $3.36$, $1.56$ and $1.29$ times larger than at $\varphi _1$, $\varphi _2$ and $\varphi _4$, respectively. This is consistent with the phase evolution of the positive peak of ${\mathsf{p}}^{cs}_{33}$ visualised in figure 11. Moreover, the singularity point lies in the source region dominated by ${\mathsf{p}}^{cs}_{33}$, and its $r_z$ position moves with $\varphi$ following the peak of ${\mathsf{p}}^{cs}_{33}$, being $r_z^+=24$, $33$, $40$ and $45$ for $\varphi _1$, $\varphi _2$, $\varphi _3$ and $\varphi _4$; for the uncontrolled case it is $r_z^+=26$.

We therefore conclude that, at least for the $T^+=250$ case discussed here, the phase dependence of the transfers of $\overline {{\delta w'' \delta w''}}$ is governed by the ${\mathsf{p}}^{cs}_{33}$ contribution to $\xi _{33}^s$ rather than by ${\rm \pi} _{33}^s$. At all phases, the largest part of $\overline {{\delta w'' \delta w''}}$ withdrawn by the source is released in the near-wall region; a relatively smaller part goes to the smallest scales, and a minimal part goes towards the channel centre, where the turbulent activity is low. By comparing the negative peaks of the divergence of the in-plane flux vector at the wall and at the smallest scales, it is established that in the uncontrolled case the amount of $\overline {{\delta w'' \delta w''}}$ released at $Y \rightarrow 0$ is $5.67$ times larger than that released at $r_z \rightarrow 0$. The oscillating wall alters the relative importance of the fluxes: the amount of $\overline {{\delta w'' \delta w''}}$ released at $Y \rightarrow 0$ is significantly less, being $2.62$, $3.85$, $2.46$ and $2.41$ times larger than that released at $r_z \rightarrow 0$ at phases $\varphi _1$, $\varphi _2$, $\varphi _3$ and $\varphi _4$, respectively.