3.1. Solenoidal and irrotational acceleration fluctuations

The Helmholtz decomposition states that a twice continously differentiable three- dimensional vector field $\boldsymbol {q}(\boldsymbol {x}, t)$ defined on a domain $V \subseteq \mathbb {R}^3$ can be expressed as the sum of an irrotational vector field $\boldsymbol {q}_I(\boldsymbol {x},t)$ and a solenoidal vector field $\boldsymbol {q}_{S}(\boldsymbol {x},t)$ (Helmholtz Reference Helmholtz1867; Stewart Reference Stewart2012; Bhatia et al. Reference Bhatia, Norgard, Pascucci and Bremer2013),

(3.1a,b)$$\begin{gather} \boldsymbol{q}_{I}(\boldsymbol{x},\,t) = -\boldsymbol{\nabla}_{x} \phi (\boldsymbol{x},\,t), \quad \boldsymbol{q}_{S}(\boldsymbol{x},\,t) = \boldsymbol{\nabla}_{x} \times \boldsymbol{B}(\boldsymbol{x},\,t), \end{gather}$$

where $\phi (\boldsymbol {x},t)$ is a scalar potential and $\boldsymbol {B}(\boldsymbol {x},t)$ is a vector potential. The Helmholtz decomposition and its interpretation can be applied to any vector field $\boldsymbol {q}(\boldsymbol {x},t)$ satisfying the above conditions, and Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) applied it to fluid accelerations and the incompressible NS equation.

The solenoidal and irrotational NS equations in homogeneous/periodic turbulence can be derived from the incompressible NS equation in Fourier space (see Appendix A). After transformation back to physical space, one obtains

(3.2)$$\begin{gather} \frac{\partial \boldsymbol{u}}{\partial t}+(\boldsymbol{u}\boldsymbol{\cdot} \boldsymbol{\nabla}_{\boldsymbol{x}}\boldsymbol{u})^{{\rm T}} = \nu\nabla_{\boldsymbol{x}}^2 \boldsymbol{u} + \boldsymbol{f}^{{\rm T}}, \end{gather}$$

(3.3)$$\begin{gather}(\boldsymbol{u}\boldsymbol{\cdot}\boldsymbol{\nabla}_{\boldsymbol{x}}\boldsymbol{u})^{L} =-\frac{1}{\rho}\boldsymbol{\nabla}_{\boldsymbol{x}} p+\boldsymbol{f}^{L}, \end{gather}$$

where the superscripts $L$ and $T$ denote fields obtained from longitudinal and transverse parts of respective Fourier vector fields (see Appendix A for precise definitions and Pope (Reference Pope2000); Stewart (Reference Stewart2012)), $p=p(\boldsymbol {x},t)$ is the pressure field and $\rho$ is the density. For any periodic vector field $\boldsymbol {q}$, $\boldsymbol {q}^{L}$ equals the irrotational field $\boldsymbol {q}_I$ and $\boldsymbol {q}^{{\rm T}}$ equals the solenoidal field $\boldsymbol {q}_S$ (see Appendix A and Stewart Reference Stewart2012). From (3.2)–(3.3), one arrives at (Tsinober et al. Reference Tsinober, Vedula and Yeung2001)

(3.4)$$\begin{gather} \frac{\partial \boldsymbol{u}}{\partial t}+(\boldsymbol{u}\boldsymbol{\cdot} \boldsymbol{\nabla}_{\boldsymbol{x}}\boldsymbol{u})_S = \nu\nabla_{\boldsymbol{x}}^2 \boldsymbol{u} + \boldsymbol{f}_S, \end{gather}$$

(3.5)$$\begin{gather}(\boldsymbol{u}\boldsymbol{\cdot}\boldsymbol{\nabla}_{\boldsymbol{x}}\boldsymbol{u})_I =-\frac{1}{\rho}\boldsymbol{\nabla}_{\boldsymbol{x}} p+\boldsymbol{f}_I, \end{gather}$$

which we refer to as Tsinober equations. Equation (3.4) contains only solenoidal vector fields and (3.5) contains only irrotational vector fields. Note that in the case of an incompressible periodic velocity field, the velocity field is solenoidal, i.e. $\boldsymbol {u}=\boldsymbol {u}_S$. This follows immediately from the scalar potential $\phi$ being the solution to $\nabla _{\boldsymbol {x}}^2\phi =0$ with periodic boundary conditions for $\boldsymbol {\nabla }_{\boldsymbol {x}}\phi$, yielding $\phi =const$.

In Appendix C we show that (3.4) is the integrated vorticity equation and that (3.5) is the integrated Poisson equation for pressure. The procedure presented in Appendix C for obtaining the Tsinober equations is also used in this same appendix to obtain generalized Tsinober equations for non-homogeneous/non-periodic turbulence with arbitrary boundary conditions.

Following the notation of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001), we define $\boldsymbol {a}_l \equiv \partial \boldsymbol {u}/\partial t$, $\boldsymbol {a}_c \equiv \boldsymbol {u}\boldsymbol{\cdot} \boldsymbol {\nabla }_{\boldsymbol {x}}\boldsymbol {u}$, $\boldsymbol {a} \equiv \boldsymbol {a}_l+\boldsymbol {a}_c$, $\boldsymbol {a}_p \equiv -1/\rho \boldsymbol {\nabla }_{\boldsymbol {x}}p$ and $\boldsymbol {a}_{\nu } \equiv \nu \nabla _{\boldsymbol {x}}^2 \boldsymbol {u}$. In such notation, (3.4)–(3.5) are $\boldsymbol {a}_l + \boldsymbol {a}_{c_S} = \boldsymbol {a}_{\nu } + \boldsymbol {f}_S$ and $\boldsymbol {a}_{c_I} = \boldsymbol {a}_p + \boldsymbol {f}_I$. Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) in fact wrote these equations for statistically homogeneous/periodic NS turbulence without body forces, i.e. with $\boldsymbol {f}=0$. In general, however, the body forcing can be considered, as in the present work, to be non-zero and typically incompressible, i.e. $\boldsymbol {f}_I = \boldsymbol {0}$ but $\boldsymbol {f}_{s}\not \equiv \boldsymbol {0}$, given that a compressible component of the forcing can be subsumed into the pressure field in incompressible turbulence. In body-forced statistically stationary homogeneous/periodic turbulence, the average forcing magnitude $\langle \boldsymbol {f}^2 \rangle$, where the brackets denote spatio-temporal averaging, tends to be small compared with $\langle \boldsymbol {a}_{\nu }^2 \rangle$ when the forcing is applied only to the largest scales (Vedula & Yeung Reference Vedula and Yeung1999). Given that $\langle \boldsymbol {f} \boldsymbol{\cdot} \boldsymbol {u}\rangle =\langle \epsilon \rangle$, where $\epsilon$ is the local turbulence dissipation rate, $\boldsymbol {f}^2$ can be quite small if $\boldsymbol {f}$ is not close to orthogonal to the velocity field. This is indeed the case with the negative damping and ABC forcings used in this study. In cases where $\boldsymbol {f}$ is close to orthogonal to the velocity field, which is conceivable in electromagnetic situations (Lorentz force), $\boldsymbol {f}^2$ needs to be large enough for $\langle \boldsymbol {f} \boldsymbol{\cdot} \boldsymbol {u}\rangle$ to balance $\langle \epsilon \rangle$. In this study we have not considered such forcings and some of our results might not be applicable to such situations. Our results for the forcings we used indicate that $\langle \boldsymbol {f}^2 \rangle$ is indeed much smaller than $\langle \boldsymbol {a}_{\nu }^2 \rangle$ (see results from our DNS in table 2) and the probability to find values of $\boldsymbol {f}^{2}$ large enough to be comparable to the other terms in the Tsinober equations is extremely small (see results from our DNS in figure 1 and table 3) where we see, in particular, that $|\boldsymbol {f}|>0.1 |\boldsymbol {a}_{c_S}|$ in $15.3\,\%$ and $6.3\,\%$ of the spatio-temporal domain for the two Reynolds numbers, respectively, the percentage being smaller for the higher Reynolds number. If we consider $|\boldsymbol {f}|>\sqrt {0.1} |\boldsymbol {a}_{c_S}| \approx 0.32 |\boldsymbol {a}_{c_S}|$, we see that this is only satisfied in $0.8\,\%$ and $0.3\,\%$ of the spatio-temporal domain, respectively. Furthermore, figure 1 and table 3 show that $\boldsymbol {f}$ is also typically much smaller than $\boldsymbol {a}_{\nu }$. We can therefore write $\boldsymbol {a}_l + \boldsymbol {a}_{c_S} \approx \boldsymbol {a}_{\nu }$, this being a good approximation in the majority of the flow for the majority of the time. With $\boldsymbol {a}_{c_I} = \boldsymbol {a}_p$ given that $\boldsymbol {f}_I = \boldsymbol {0}$, these two equations are very close to the way that Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) originally wrote them ($\boldsymbol {a}_l + \boldsymbol {a}_{c_S} = \boldsymbol {a}_{\nu }$ and $\boldsymbol {a}_{c_I} = \boldsymbol {a}_p$ for the ${\bf f}\equiv 0$ case) and we can therefore expect our DNS to retrieve the DNS results and conclusions of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001).

Figure 1. Probability density functions (PDFs) $P$ of NS acceleration and force magnitudes $\boldsymbol {q}^2$ for terms $\boldsymbol {q}$ listed at the top of (a). Here $P_{{max}}$ for the PDF of $\boldsymbol {q}^2$ denotes its maximum value. Results are shown for (a) $\langle {{Re}}_\lambda \rangle _t=112$, (b) $\langle {{Re}}_\lambda \rangle _t=174$.

Table 2. Normalised average magnitudes $\langle \boldsymbol {q}^2 \rangle / (3\langle \epsilon \rangle ^{3/2} \nu ^{-1/2})$ and $\langle \boldsymbol {q}^2 \rangle /\langle \boldsymbol {a}_{c}^2 \rangle$ for NS accelerations and forces $\boldsymbol {q}$ defined in the fourth paragraph of § 3.1 for our two $\langle {{Re}}_\lambda \rangle _t$. The accelerations and forces $\boldsymbol {q}$ are listed on the top row, $\boldsymbol {q}^2 \equiv q_i q_i$, $\epsilon$ denotes the viscous dissipation rate and $\langle \cdots \rangle$ denotes a spatio-temporal average.

Table 3. Probabilities of events $\boldsymbol {q}^2> \alpha \boldsymbol {p}^2$ for NS terms $(\boldsymbol {q},\boldsymbol {p})$ with $\alpha$ specified on the top row. The two probability values in the brackets for each $(\boldsymbol {q},\boldsymbol {p},\alpha )$ combination refer to $\langle {{Re}}_\lambda \rangle _t=112$ and $\langle {{Re}}_\lambda \rangle _t=174$, respectively.

The DNS of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) showed that $\boldsymbol {a}_{\nu }$ is typically negligible (i.e. in a statistical sense, not everywhere at any time in the flow) compared with all the other acceleration terms in the Tsinober equations, namely $\boldsymbol {a}_l$, $\boldsymbol {a}_{c_S}$, $\boldsymbol {a}_{c_I}$ and $\boldsymbol {a}_p$. This is confirmed by our DNS results in tables 2–3 and in figure 1 that are for similar Reynolds numbers to those of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) and where we report root-mean-square values, and probabilities of various acceleration terms. It is worth noting that $\boldsymbol {a}_{\nu }$ is not everywhere always negligible, at these Reynolds numbers at least. For example, $|\boldsymbol {a}_{\nu }|>0.1 |\boldsymbol {a}_{c_S}|$ in $44.1\,\%$ and $30.8\,\%$ of the space–time domain for our lower and higher Reynolds number, respectively; and if we consider $|\boldsymbol {a}_{\nu }|>0.32 |\boldsymbol {a}_{c_S}|$, this is satisfied in $6.8\,\%$ and $3.7\,\%$ of cases. Note that the DNS results of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) suggest that the viscous force typically decreases in magnitude compared with $\boldsymbol {a}_{c_S}$ as the Reynolds number increases and our results for our two Reynolds numbers agree with this trend. One may therefore expect the first of the two Tsinober equations for homogeneous/periodic turbulence with the kind of forcing we consider here to typically reduce to

(3.6)\begin{equation} \boldsymbol{a}_l + \boldsymbol{a}_{c_S} \approx 0 \end{equation}

at high enough Reynolds numbers, the approximation being valid in the sense that the neglected terms are significantly smaller than the retained ones in the majority of the flow for the majority of the time. There exist, however, some relatively rare spatio-temporal events where the neglegted viscous force and/or body force are significant (for example, as stated a few lines above, $|\boldsymbol {a}_{\nu }|$ is larger than $0.32 |\boldsymbol {a}_{c_S}|$ in $6.8\,\%$ and $3.7\,\%$ of all spatio-temporal events for our lower and higher Reynolds numbers, respectively) and where the right-hand side of (3.6) is therefore not zero. In fact, many of these relatively rare events can be expected to account for some or even much of the average turbulence dissipation that is a sum of squares of fluctuating velocity gradients. More generally, one cannot use (3.6) to derive statistics of fluctuating velocity gradients, as in Tang, Antonia & Djenidi (Reference Tang, Antonia and Djenidi2022) for example.

The second of the two Tsinober equations, namely

(3.7)\begin{equation} \boldsymbol{a}_{c_I} = \boldsymbol{a}_p, \end{equation}

is exact everywhere and at any time and we keep it as it is.

Equations (3.6)–(3.7) suggest similar magnitudes and strong alignment between $\boldsymbol {a}_l$ and $-\boldsymbol {a}_{c_S}$ and equal magnitudes as well as perfect alignment between $\boldsymbol {a}_{c_I}$ and $\boldsymbol {a}_p$. Such magnitudes and alignments were observed in the DNS of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) and are also strongly confirmed by our own DNS in table 4 ($\boldsymbol {a}_{c_S}$ and $\boldsymbol {a}_{c_I}$ are calculated on the basis of (A1a,b) in Appendix A and aliasing errors associated with nonlinear terms are removed by phase shifting and truncation Patterson & Orszag Reference Patterson and Orszag1971). As suggested by previous DNS and experimental results (e.g. Tsinober et al. Reference Tsinober, Vedula and Yeung2001; Chevillard et al. Reference Chevillard, Roux, Lévêque, Mordant, Pinton and Arnéodo2005; Yeung et al. Reference Yeung, Pope, Lamorgese and Donzis2006), and as also supported by our own DNS results in tables 2 and 4, $\boldsymbol {a}\approx \boldsymbol {a}_p$ and $\langle \boldsymbol {a}_{l}^2 \rangle /\langle \boldsymbol {a}^2 \rangle \sim \langle {{Re}}_{\lambda } \rangle _t^{1/2}$. In fact, the scaling $\langle \boldsymbol {a}_{l}^2 \rangle /\langle \boldsymbol {a}^2 \rangle \sim \langle {{Re}}_{\lambda } \rangle _t^{1/2}$ follows from the analysis of Tennekes (Reference Tennekes1975) who expressed the concept of passive sweeping by pointing out that ‘at high Reynolds number the dissipative eddies flow past an Eulerian observer in a time much shorter than the time scale which characterizes their own dynamics’. It then follows from (3.6)–(3.7), from $\langle \boldsymbol {a}_{l}^2 \rangle /\langle \boldsymbol {a}^2 \rangle \sim \langle {{Re}}_{\lambda } \rangle _t^{1/2}$ and from $\langle \boldsymbol {a}_{p}^2 \rangle \approx \langle \boldsymbol {a}^2 \rangle$ that $\langle \boldsymbol {a}_{c_S}^2 \rangle /\langle \boldsymbol {a}_{c_I}^2 \rangle \sim \langle {{Re}}_{\lambda } \rangle _t^{1/2}$ with increasing $\langle Re_\lambda \rangle _t$, i.e. $\boldsymbol {a}_c$ becomes increasingly solenoidal with increasing $\langle Re_\lambda \rangle _t$. In this way, the anti-alignment in (3.6) leads to an increasing anti-alignment tendency between $\boldsymbol {a}_l$ and $\boldsymbol {a}_c$ with increasing Reynolds number, which is consistent with the notion of passive sweeping of small eddies by large ones, i.e. the random Taylor hypothesis of Tennekes (Reference Tennekes1975). These observations and conclusions were all made by Tsinober et al. (Reference Tsinober, Vedula and Yeung2001). They are now confirmed by our DNS results in table 2 and this reiterates that they do not require a large Taylor length-based Reynolds number to emerge.

Table 4. NS average alignments $\langle \cos (\boldsymbol {q},\boldsymbol {p})\rangle$ for NS acceleration pairs $(\boldsymbol {q},\boldsymbol {p})$.

As a final point, it is a general property of isotropic random vector fields $\boldsymbol {q}$ that $\langle \boldsymbol {q}_I(\boldsymbol {x},t) \boldsymbol{\cdot} \boldsymbol {q}_S(\boldsymbol {x} + \boldsymbol {r},t)\rangle _x=0$ for any $\boldsymbol {r}$ (including $\boldsymbol {r}=0$), where $\langle \cdots \rangle _{\boldsymbol {x}}$ signifies a spatial average (Monin, Yaglom & Lumley Reference Monin, Yaglom and Lumley1975). Thus, $\langle \boldsymbol {a}_{c}^2\rangle = \langle \boldsymbol {a}_{c_I}^2\rangle + \langle \boldsymbol {a}_{c_S}^2\rangle$ if the small-scale turbulence is isotropic. Both our DNS and the DNS of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) confirm this equality. From this equality and from (3.6), $\langle \boldsymbol {a}_{c_S}^2 \rangle /\langle \boldsymbol {a}_{c_I}^2 \rangle \sim \langle {{Re}}_{\lambda } \rangle _t^{1/2}$, (3.7), $\boldsymbol {a} \approx \boldsymbol {a}_{p}$ and $\langle \boldsymbol {a}^{2}\rangle \gg \langle \boldsymbol {a}_{\nu }^{2}\rangle \gg \langle \boldsymbol {f}^{2}\rangle$, we have all in all

(3.8)\begin{equation} \langle \boldsymbol{a}_{c}^2 \rangle \ge \langle \boldsymbol{a}_{c_S}^2 \rangle \approx \langle \boldsymbol{a}_{l}^2 \rangle \gg \langle \boldsymbol{a}_{c_I}^2 \rangle = \langle\boldsymbol{a}_p^2 \rangle \approx \langle\boldsymbol{a}^2 \rangle \gg \langle \boldsymbol{a}_{\nu}^2 \rangle \gg \langle \boldsymbol{f}^2 \rangle \end{equation}

for large enough $\langle {{Re}}_{\lambda } \rangle _t$. The average magnitude ordering in (3.8) is confirmed in our DNS (see table 2) and the DNS of Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) even though the Reynolds numbers of these DNS are moderate and so the difference between $\langle \boldsymbol {a}_{c_I}^2 \rangle$ and $\langle \boldsymbol {a}_l^2\rangle$ is not so large. Tsinober's way to formulate sweeping is encapsulated in $\langle \boldsymbol {a}_{c_S}^2 \rangle \approx \langle \boldsymbol {a}_{l}^2 \rangle \gg \langle \boldsymbol {a}_{c_I}^2 \rangle = \langle \boldsymbol {a}_p^2 \rangle \approx \langle \boldsymbol {a}^2 \rangle$ and in the alignments implied by (3.6) and (3.7) that are also statistically confirmed by our DNS in table 4.

3.2. From one-point to two-point NS dynamics in periodic/homogeneous turbulence

The NSD equation at centroid $\boldsymbol {x}$ and separation vector $\boldsymbol {r}$ is derived by subtracting the NS equation at location $\boldsymbol {x}^+=\boldsymbol {x}+\boldsymbol {r}/2$ from the NS equation at location $\boldsymbol {x}^-=\boldsymbol {x}-\boldsymbol {r}/2$. Defining $\delta \boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t)\equiv \boldsymbol {q}(\boldsymbol {x}+\boldsymbol {r}/2,t)-\boldsymbol {q}(\boldsymbol {x}-\boldsymbol {r}/2,t)$ for any NS term $\boldsymbol {q}(\boldsymbol {x},t)$, the NSD equation (Hill Reference Hill2002) reads

(3.9)\begin{equation} \frac{\partial \delta \boldsymbol{u}}{\partial t}+\delta \boldsymbol{a}_c =-\frac{1}{\rho}\boldsymbol{\nabla}_{\boldsymbol{x}}\delta p +\delta \boldsymbol{a}_\nu +\delta \boldsymbol{f}. \end{equation}

The NSD equation governs the evolution of $\delta \boldsymbol {u}$, which can be thought of as pertaining to the momentum at scales smaller or comparable to $|\boldsymbol {r}|$. We derive the solenoidal NSD equation by subtracting (3.4) at $\boldsymbol {x}-\boldsymbol {r}/2$ from the same equation at $\boldsymbol {x}+\boldsymbol {r}/2$. The same operation is used to derive the irrotational NSD equation. The resulting equations read

(3.10)$$\begin{gather} \frac{\partial \delta \boldsymbol{u}}{\partial t}+\delta \boldsymbol{a}_{c_S} = \delta \boldsymbol{a}_\nu +\delta\boldsymbol{f}_S, \end{gather}$$

(3.11)$$\begin{gather}\delta \boldsymbol{a}_{c_I} =-\frac{1}{\rho}\boldsymbol{\nabla}_{\boldsymbol{x}}\delta p + \delta \boldsymbol{f}_I, \end{gather}$$

where $\delta \boldsymbol {a}_{c_S}(\boldsymbol {x},\boldsymbol {r},t) \equiv \boldsymbol {a}_{c_S}(\boldsymbol {x}+\boldsymbol {r}/2,t)- \boldsymbol {a}_{c_S}(\boldsymbol {x}-\boldsymbol {r}/2,t)$ and $\delta \boldsymbol {a}_{c_I}(\boldsymbol {x},\boldsymbol {r},t) \equiv \boldsymbol {a}_{c_I}(\boldsymbol {x}+\boldsymbol {r}/2,t) - \boldsymbol {a}_{c_I}(\boldsymbol {x}-\boldsymbol {r}/2,t)$ and note that these terms refer to solenoidal and irrotational terms in $\boldsymbol {x}$ space rather than $\boldsymbol {r}$ space. The forcings we consider have no irrotational part and so $\delta \boldsymbol {f}_I = 0$. At the moderate $\langle {{Re}}_\lambda \rangle _t$ of our DNS, the approximate (3.6) is valid in the sense explained in the text that accompanies it in the previous subsection, i.e. for a majority of spatio-temporal events. If the magnitude of the separation vector $\boldsymbol {r}$ is not too small for viscosity to matter directly nor too large for the forcing to be directly present, we may safely subtract (3.6) at $\boldsymbol {x}-\boldsymbol {r}/2$ from (3.6) at $\boldsymbol {x}+\boldsymbol {r}/2$ to obtain an approximation to (3.10) for a sufficiently high Reynolds number: this is the first of the two equations below where $\delta \boldsymbol {a}_l \equiv \partial \delta \boldsymbol {u}/\partial t$,

(3.12)$$\begin{gather} \delta \boldsymbol{a}_{l} + \delta \boldsymbol{a}_{c_S}\approx 0, \end{gather}$$

(3.13)$$\begin{gather}\delta \boldsymbol{a}_{c_I} =-\frac{1}{\rho}\boldsymbol{\nabla}_{\boldsymbol{x}}\delta p. \end{gather}$$

The second equation (3.13) follows directly from (3.11) with $\delta \boldsymbol {f}_I = 0$ without any restriction on either $\boldsymbol {r}$ or Reynolds number and is exact.

Like (3.6), (3.12) can be expected to be valid broadly except where and when $\delta \boldsymbol {a}_\nu +\delta \boldsymbol {f}_S$ is large enough not to be negligible. Figure 2 shows statistically converged estimations of exceedance probabilities of NSD viscous and external force terms that suggest that (3.12) is indeed a good approximation for most of space and time at the Reynolds numbers of our two DNS, at the very least for separation distances larger than $\langle \lambda \rangle _{t}$ and smaller than $\langle L\rangle _t$. With regards to the forcing, $\text {Prob}(|\delta \boldsymbol {f}|>0.32 |\delta \boldsymbol {a}_{c_S}|)$ is typically of the order of $1\,\%$, in particular for our higher Reynolds number. With regards to the viscous force, $\text {Prob}(|\delta \boldsymbol {a}_{\nu }|>0.32 |\delta \boldsymbol {a}_{c_S}|)$ is typically of the order of $5\,\%$ for $r\ge \langle \lambda \rangle _{t}$ and even less for our higher Reynolds number.

Figure 2. The NSD exceedance probabilities $\text {Prob}(\boldsymbol {q}^2>\alpha \boldsymbol {p}^2)$ for the NSD terms on top of (a) as a function of separation length $r_d=|\boldsymbol {r}|$. The legend entries read $(\boldsymbol {q},\alpha, \boldsymbol {p})$ for the NSD terms introduced in the first paragraph of § 3.2. Results are shown for (a) $\langle {{Re}}_\lambda \rangle _t = 112$, $\langle L\rangle _{t} = 3.5 \langle \lambda \rangle _{t}$; (b) $\langle {{Re}}_\lambda \rangle _t = 174$, $\langle L\rangle _{t} = 5.2 \langle \lambda \rangle _{t}$. The NSD terms are sampled at scale $r_d=|\boldsymbol {r}|$ at random orientations $\boldsymbol {r}$.

The link between nonlinearity and non-locality (via the pressure field) invoked in the two-point analysis of Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) has its root in (3.13) that parallels (3.7) and states that $\delta \boldsymbol {a}_{c_I}$ and $\delta \boldsymbol {a}_p$ are perfectly aligned and have the same magnitudes. Furthermore, similarly to the way that (3.6) supports the concept of sweeping of small turbulent eddies by large ones in the usual one-point sense, (3.12) suggests similar magnitudes for and strong alignment between $\delta \boldsymbol {a}_l$ and $-\delta \boldsymbol {a}_{c_S}$. A two-point concept of sweeping similar to that of Tennekes (Reference Tennekes1975) that relies on alignment between $\delta \boldsymbol {a}_l$ and $-\delta \boldsymbol {a}_{c}$ should also require that $\delta \boldsymbol {a}_{c}$ tends towards $\delta \boldsymbol {a}_{c_s}$ with increasing Reynolds number, i.e. that $\delta \boldsymbol {a}_c$ becomes increasingly solenoidal. We therefore seek to obtain inequalities and approximate equalities similar to (3.8). Note that (3.12) and (3.13) immediately imply $\langle \delta \boldsymbol {a}_{c_S}^2 \rangle \approx \langle \delta \boldsymbol {a}_{l}^2 \rangle$, $\langle \delta \boldsymbol {a}_{c_I}^2 \rangle = \langle \delta \boldsymbol {a}_p^2 \rangle$ and $\langle \delta \boldsymbol {a}_p^2 \rangle \approx \langle \delta \boldsymbol {a}^2 \rangle$. It therefore remains to argue that $\langle \delta \boldsymbol {a}_{c}^2 \rangle \ge \langle \delta \boldsymbol {a}_{c_S}^2 \rangle \gg \langle \delta \boldsymbol {a}_{c_I}^2 \rangle$, which is exactly what we need to complete the new concept of two-point sweeping.

We start from

(3.14)$$\begin{gather} \langle \delta \boldsymbol{q} \boldsymbol{\cdot} \delta \boldsymbol{q} \rangle (\boldsymbol{r}) = \langle \boldsymbol{q}^+ \boldsymbol{\cdot} \boldsymbol{q}^+ \rangle- \langle \boldsymbol{q}^+ \boldsymbol{\cdot} \boldsymbol{q}^- \rangle + \langle \boldsymbol{q}^- \boldsymbol{\cdot} \boldsymbol{q}^- \rangle- \langle \boldsymbol{q}^- \boldsymbol{\cdot} \boldsymbol{q}^+ \rangle, \end{gather}$$

(3.15)$$\begin{gather}= 2 [\langle \boldsymbol{q} \boldsymbol{\cdot} \boldsymbol{q} \rangle - \langle \boldsymbol{q}^+ \boldsymbol{\cdot} \boldsymbol{q}^- \rangle(\boldsymbol{r})], \end{gather}$$

where $\boldsymbol {q}^+ \equiv \boldsymbol {q}(\boldsymbol {x}+\boldsymbol {r}/2)$ and $\boldsymbol {q}^- \equiv \boldsymbol {q}(\boldsymbol {x}-\boldsymbol {r}/2)$ and where we used $\langle \boldsymbol {q}^+ \boldsymbol{\cdot} \boldsymbol {q}^+ \rangle =\langle \boldsymbol {q}^- \boldsymbol{\cdot} \boldsymbol {q}^- \rangle = \langle \boldsymbol {q} \boldsymbol{\cdot} \boldsymbol {q} \rangle$ because of statistical homogeneity/periodicity. Previous studies (Hill & Thoroddsen Reference Hill and Thoroddsen1997; Vedula & Yeung Reference Vedula and Yeung1999; Gulitski et al. Reference Gulitski, Kholmyansky, Kinzelbach, Lüthi, Tsinober and Yorish2007; Xu et al. Reference Xu, Ouellette, Vincenzi and Bodenschatz2007) demonstrated that fluid accelerations, pressure gradients and viscous forces have limited spatial correlations in terms of alignments at scales larger than $\langle \lambda \rangle _{t}$ for moderate and high $\langle {{Re}}_\lambda \rangle _t$. Thus, if we assume the two-point term to be negligible compared with the one-point term in (3.15) for scales $|\boldsymbol {r}|$ larger than $\langle \lambda \rangle _{t}$, we have that $\langle \delta \boldsymbol {q} \boldsymbol{\cdot} \delta \boldsymbol {q} \rangle (\boldsymbol {r})$ is approximately equal to $2 \langle \boldsymbol {q}\boldsymbol{\cdot} \boldsymbol {q} \rangle$ for $\vert \boldsymbol {r}\vert$ larger than $\langle \lambda \rangle _{t}$. From (3.8) we therefore obtain

(3.16)\begin{equation} \langle \delta \boldsymbol{a}_{c}^2 \rangle \ge \langle \delta\boldsymbol{a}_{c_S}^2 \rangle \approx \langle \delta \boldsymbol{a}_{l}^2 \rangle \gg \langle \delta\boldsymbol{a}_{c_I}^2 \rangle = \langle\delta\boldsymbol{a}_p^2 \rangle \approx \langle \delta\boldsymbol{a}^2 \rangle \gg \langle\delta \boldsymbol{a}_{\nu}^2 \rangle \gg \langle \delta\boldsymbol{f}^2 \rangle \end{equation}

for $\vert \boldsymbol {r}\vert$ larger than $\langle \lambda \rangle _{t}$, but $\langle \delta \boldsymbol {a}_{c}^2 \rangle \ge \langle \delta \boldsymbol {a}_{c_S}^2 \rangle$ and $\langle \delta \boldsymbol {a}_{c_I}^2 \rangle = \langle \delta \boldsymbol {a}_p^2 \rangle$ are in fact valid for any $\boldsymbol {r}$. Inequality $\langle \delta \boldsymbol {a}_{c}^2 \rangle \ge \langle \delta \boldsymbol {a}_{c_S}^2 \rangle$ follows from $\langle \delta \boldsymbol {a}_c^2 \rangle = \langle \delta \boldsymbol {a}_{c_I}^2 \rangle + \langle \delta \boldsymbol {a}_{c_S}^2 \rangle$ that itself follows from $\langle \boldsymbol {a}_{c_I}(\boldsymbol {x},t) \boldsymbol{\cdot} \boldsymbol {a}_{c_S}(\boldsymbol {x}+\boldsymbol {r},t)\rangle _x=0$ for any $\boldsymbol {r}$ if the turbulence is isotropic (Monin et al. Reference Monin, Yaglom and Lumley1975). Equality $\langle \delta \boldsymbol {a}_{c_I}^2 \rangle = \langle \delta \boldsymbol {a}_p^2\rangle$ follows directly from (3.13) that is exact and holds for any $\boldsymbol {r}$ and any Reynolds number. Of equalities/inequalities (3.16), the ones that we did not already directly derive from/with (3.12) and (3.13) are $\langle \delta \boldsymbol {a}_{c}^2 \rangle \ge \langle \delta \boldsymbol {a}_{c_S}^2 \rangle \gg \langle \delta \boldsymbol {a}_{c_I}^2 \rangle$ and $\langle \delta \boldsymbol {a}_{\nu }^2 \rangle \gg \langle \delta \boldsymbol {f}^2 \rangle$. The present way to formulate the new concept of two-point sweeping follows from Tsinober's way to formulate sweeping and is encapsulated in $\delta \boldsymbol {a}_{c_S}^2 \rangle \approx \langle \delta \boldsymbol {a}_{l}^2 \rangle \gg \langle \delta \boldsymbol {a}_{c_I}^2 \rangle = \langle \delta \boldsymbol {a}_p^2 \rangle \approx \langle \delta \boldsymbol {a}^2 \rangle$ and in the alignments implied by (3.12) and (3.13). We confirm (3.12), (3.13) and (3.16) with our DNS in the remainder of this subsection.

To test (3.16) with our DNS data in a manageable way, we calculate spatio-temporal averages of $\boldsymbol {r}$-orientation-averaged quantities

(3.17)\begin{equation} (\delta \boldsymbol{q} \boldsymbol{\cdot} \delta \boldsymbol{q} )^a(\boldsymbol{x},r_d,t) \equiv \frac{1}{{\rm \pi} r_d^2} \iiint_{|\boldsymbol{r}|=r_d} \delta \boldsymbol{q}(\boldsymbol{x},\boldsymbol{r},t) \boldsymbol{\cdot} \delta \boldsymbol{q}(\boldsymbol{x},\boldsymbol{r},t), {\rm d}\boldsymbol{r}, \end{equation}

which we plot in figure 3(a1,a2) as ratios of such quantities versus the two-point length $r_d$. In figure 3(a1,a2) we plot spatio-temporal averages of $\boldsymbol {r}$-orientation-averaged quantities (3.17) for various acceleration/force terms in the NSD and the Helmholtz decomposed NSD equations. A comparison of relative magnitudes in the plots of figure 3(a1,a2) with relative magnitudes in table 2 makes it clear that the results are consistent with (3.16) and $\langle \delta \boldsymbol {q} \boldsymbol{\cdot} \delta \boldsymbol {q} \rangle (\boldsymbol {r}) /\langle \boldsymbol {q}\boldsymbol{\cdot} \boldsymbol {q} \rangle$ close to $2$ for $r_{d} \geq \langle \lambda \rangle _t$ at both $\langle Re_\lambda \rangle _t$ to a good degree of accuracy ($\langle \delta \boldsymbol {q} \boldsymbol{\cdot} \delta \boldsymbol {q} \rangle (\boldsymbol {r}) /\langle \boldsymbol {q}\boldsymbol{\cdot} \boldsymbol {q} \rangle$ increases from 1.8 to 2.0 as $r_d$ grows from $\langle \lambda \rangle _t$ to $\langle L\rangle _t$). Note, in particular, that in figure 3(a1,b1) the average quantities corresponding to $\delta \boldsymbol {a}_l$ and $\delta \boldsymbol {a}_{c_S}$ overlap and those corresponding to $\delta \boldsymbol {a}_p$, $\delta \boldsymbol {a}$ and $\delta \boldsymbol {a}_{c_I}$ also overlap. At scales below $\langle \lambda \rangle _t$, the average relative magnitudes change slightly, but the NSD magnitude separations still abide by (3.16), the NSD analogue to (3.8), at all scales.

Figure 3. (a1,b1) Spatio-temporal averages of spherically averaged NSD magnitudes $(\delta \boldsymbol {q}^2)^a \equiv ({\rm \pi} r_d^2)^{-1} \iiint _{|\boldsymbol {r}|=r_d} \delta \boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t) \boldsymbol{\cdot} \delta \boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t),\ {\rm d}\boldsymbol {r}$ for NSD terms $\delta \boldsymbol {q}$ listed on top of the figures as a function of $r_d$: (a1) $\langle {{Re}}_\lambda \rangle _t=112$, (b1) $\langle {{Re}}_\lambda \rangle _t=174$. The magnitudes of the terms $\delta \boldsymbol {a}_l$ and $\delta \boldsymbol {a}_{c_S})$ overlap and the magnitudes of the terms $(\delta \boldsymbol {a}_p$, $\delta \boldsymbol {a}$ and $\delta \boldsymbol {a}_{c_I})$ also overlap. (a2,b2) Average NSD alignments between NSD terms $(\delta \boldsymbol {q},\delta \boldsymbol {w})$ listed on top of the figures as a function of $r_d$: (a2) $\langle {{Re}}_\lambda \rangle _t=112$, (b2) $\langle {{Re}}_\lambda \rangle _t=174$.

In figure 3(b1,b2) we use our DNS data to plot spatio-temporal averages of $\boldsymbol {r}$-orientation-averaged cosines of angles between various NSD terms $\delta \boldsymbol {q}$ and $\delta \boldsymbol {w}$ to test for average alignments as a function of $r_d$. These alignment results are of course in perfect agreement with (3.13) but they are also in good agreement with (3.12) and acceptable agreement with $\delta \boldsymbol {a} \approx \delta \boldsymbol {a}_{p}$ (the cosine of the angle between these two acceleration vectors is higher than $0.9$ for all $r_d$). They also show that we should not expect $\delta \boldsymbol {a}_l$ to be extremely well aligned with $-\delta \boldsymbol {a}_c$ at our moderate Reynolds numbers. This demonstrates the pertinence of the solenoidal–irrotational decomposition that has revealed very good alignments at our moderate Reynolds numbers for which there are significantly weaker alignments without this decomposition.

In conclusion, figure 3 provides strong support for (3.12), (3.13) and (3.16) that establish the two-point link between nonlinearity and non-locality, and also a concept of two-point sweeping.

3.3. Interscale transfer and physical space transport accelerations

The convective nonlinearity is responsible for nonlinear turbulence transport through space and nonlinear transfer through scales. We want to separate these two effects and, therefore, decompose the two-point nonlinear acceleration term $\delta \boldsymbol {a}_c$ into an interscale transfer acceleration $\boldsymbol {a}_{\varPi }$ and a physical space transport acceleration $\boldsymbol {a}_{\mathcal {T}}$ (Hill Reference Hill2002), i.e. $\delta \boldsymbol {a}_c = \boldsymbol {a}_{\varPi }+\boldsymbol {a}_\mathcal {T}$ with

(3.18a,b) \begin{equation} \boldsymbol{a}_{\mathcal{T}} (\boldsymbol{x},\boldsymbol{r},t) = \tfrac{1}{2}(\boldsymbol{u}^{+} + \boldsymbol{u}^{-})\boldsymbol{\cdot} \boldsymbol{\nabla}_{x}\delta{\boldsymbol{u}}, \quad \boldsymbol{a}_{\varPi}(\boldsymbol{x},\boldsymbol{r},t) = \delta{\boldsymbol{u}} \boldsymbol{\cdot} \boldsymbol{\nabla}_{r} \delta{\boldsymbol{u}} \end{equation}

With this decomposition of the nonlinear term, the NSD equation (3.9) reads

(3.19)\begin{equation} \frac{\partial \delta \boldsymbol{u}}{\partial t}+\boldsymbol{a}_{{\varPi}}+ \boldsymbol{a}_{\mathcal{T}} =-\frac{1}{\rho}\boldsymbol{\nabla}_{\boldsymbol{x}}\delta p + \delta \boldsymbol{a}_{\nu} + \delta \boldsymbol{f} . \end{equation}

We note relations $\boldsymbol {a}_{{\varPi }} = \delta \boldsymbol {a}_{C} + u^+_{j}({\partial }/{\partial x^-_{j}}) \boldsymbol {u}^- -u^-_{j}({\partial }/{\partial x^+_{j}}) \boldsymbol {u}^+$ and $\boldsymbol {a}_{\mathcal {T}}=\delta \boldsymbol {a}_{C} - u^+_{j}({\partial }/{\partial x^-_{j}}) \boldsymbol {u}^- + u^-_{j}({\partial }/{\partial x^+_{j}}) \boldsymbol {u}^+$ that can be easily used to show that $\langle \boldsymbol {a}_{{\varPi }}^{2}\rangle$ and $\langle \boldsymbol {a}_{\mathcal {T}}^{2}\rangle$ tend towards each other as the amplitude of the separation vector $\boldsymbol {r}$ grows above the integral length scale. We report DNS evidence of this tendency below in this paper.

We want to consider the effects of the interscale transfer and interspace transport terms in the solenoidal and irrotational NSD dynamics and we therefore need to break down the NSD equation (3.19) into two equations, one irrotational and one solenoidal. We therefore perform Helmholtz decompositions in centroid space $\boldsymbol {x}$ for a given separation $\boldsymbol {r}$ at time $t$, for example, $\delta \boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t)= \delta \boldsymbol {q}_{\bar {I}}(\boldsymbol {x},\boldsymbol {r},t)+\delta \boldsymbol {q}_{\bar {S}}(\boldsymbol {x},\boldsymbol {r},t)$, where $\delta \boldsymbol {q}_{\bar {I}}(\boldsymbol {x},\boldsymbol {r},t)$ and $\delta \boldsymbol {q}_{\bar {S}}(\boldsymbol {x},\boldsymbol {r},t)$ are, respectively, the irrotational and solenoidal parts in centroid space of $\delta \boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t)$. This decomposition in centroid space differs in general from the difference of the Helmholtz decomposed terms in the NS equations that gives (3.10) and (3.11), but in periodic/homogeneous turbulence $\delta \boldsymbol {q}_{I} = \delta \boldsymbol {q}_{{\bar {I}}}$ and $\delta \boldsymbol {q}_{S} = \delta \boldsymbol {q}_{{\bar {S}}}$ (see Appendix B). Furthermore, from $\delta \boldsymbol {a}_c = \boldsymbol {a}_{\varPi }+\boldsymbol {a}_\mathcal {T}$, immediately follow $\delta \boldsymbol {a}_{c_{\bar {S}}}=\boldsymbol {a}_{{\varPi }_{\bar {S}}}+\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\delta \boldsymbol {a}_{c_{\bar {I}}}=\boldsymbol {a}_{{\varPi }_{\bar {I}}}+\boldsymbol {a}_{\mathcal {T}_{\bar {I}}}$. Thus, we can rewrite the NSD solenoidal and irrotational equations (3.10) and (3.11) as

(3.20)$$\begin{gather} \boldsymbol{a}_{{\varPi}_{\bar{I}}}+\boldsymbol{a}_{\mathcal{T}_{\bar{I}}} = \delta \boldsymbol{a}_p, \end{gather}$$

(3.21)$$\begin{gather}\delta \boldsymbol{a}_l +\boldsymbol{a}_{{\varPi}_{\bar{S}}}+\boldsymbol{a}_{\mathcal{T}_{\bar{S}}} = \delta \boldsymbol{a}_{\nu}+ \delta \boldsymbol{f} , \end{gather}$$

in periodic/homogeneous turbulence.

We emphasize that the interscale transfer term $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ is related non-locally in space to two-point vortex stretching and compression terms governing the evolution of vorticity difference $\delta \boldsymbol {\omega }$. This follows from the fact that, as for the Tsinober equations, the NSD solenoidal equation is an integrated vorticity difference equation. We provide mathematical detail on the connection between $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ and $\delta \boldsymbol {\omega }$ in Appendix C. This relation between $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ and the vorticity difference dynamics provides an instantaneous connection between the interscale momentum dynamics and two-point vorticity stretching and compression dynamics.

Equation (3.20) can also be obtained by integrating the Poisson equation for $\delta p$ in centroid space similarly to (3.21) that, as already mentioned, can be obtained by integrating the vorticity difference equation in that same space. We use this approach in Appendix C to derive these equations for periodic/homogeneous turbulence but also their generalized form for non-homogeneous turbulence. By deriving the exact equations for $\boldsymbol {a}_{\mathcal {T}_{\bar {I}}}(\boldsymbol {x},\boldsymbol {r},t)$ and $\boldsymbol {a}_{{\varPi }_{\bar {I}}}(\boldsymbol {x},\boldsymbol {r},t)$ in Fourier centroid space we show in Appendix B that we have $\boldsymbol {a}_{\mathcal {T}_{\bar {I}}}(\boldsymbol {x},\boldsymbol {r},t)=\boldsymbol {a}_{{\varPi }_{\bar {I}}}(\boldsymbol {x},\boldsymbol {r},t)$ in periodic/homogeneous turbulence. This result combined with (3.20) yields

(3.22)\begin{equation} \boldsymbol{a}_{{\varPi}_{\bar{I}}} = \boldsymbol{a}_{\mathcal{T}_{\bar{I}}} = \tfrac{1}{2}\delta \boldsymbol{a}_p = \tfrac{1}{2}\delta \boldsymbol{a}_{c_I} \end{equation}

in periodic/homogeneous turbulence. In figure 4 we plot spatio-temporal averages of $\boldsymbol {r}$-orientation-averaged quantities (3.17) for various acceleration/force terms in the NSD and the Helmholtz decomposed NSD equations and in figure 5 we plot spatio-temporal averages of $\boldsymbol {r}$-orientation-averaged cosines of angles between various two-point acceleration terms in these equations. The overlapping magnitudes in figure 4 and the average alignments in figure 5 confirm (3.22), or rather validate our DNS given that (3.22) is exact.

Figure 4. Average magnitudes $\langle \delta \boldsymbol {q}^2 \rangle ^a$ of NSD terms present in the irrotational and solenoidal NSD equations (3.21) and (3.22) listed on top of (a). All values have been normalised with $\langle \delta \boldsymbol {a}_{c_S}^2 \rangle ^a$ at the largest considered separation $r_d$. The magnitudes of the terms $(\delta \boldsymbol {a}_l+\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\boldsymbol {a}_{{\varPi }_{\bar {S}}})$ overlap and the magnitudes of the terms $(1/2\delta \boldsymbol {a}_{c_I}$, $\boldsymbol {a}_{\mathcal {T}_{\bar {I}}}$ and $\boldsymbol {a}_{{\varPi }_{\bar {I}}})$ also overlap. Results are shown for (a) $\langle {{Re}}_\lambda \rangle _t=112$, (b) $\langle {{Re}}_\lambda \rangle _t=174$.

Figure 5. Average alignments of NSD terms $(\delta \boldsymbol {q},\delta \boldsymbol {w})$ listed on top of (a,b). The average alignments of $(\delta \boldsymbol {a}_p,\boldsymbol {a}_{\mathcal {T}_{\bar {I}}})$ and $(\delta \boldsymbol {a}_p,\boldsymbol {a}_{{\varPi }_{\bar {I}}})$ overlap: (a) $\langle {{Re}}_\lambda \rangle _t=112$, (b) $\langle {{Re}}_\lambda \rangle _t=174$.

The computational procedure to calculate the various $\boldsymbol {r}$-orientation-averaged terms in these figures is computationally expensive. To calculate the NSD irrotational and solenoidal parts of the interscale and interspace transport terms at a given time $t$ and separation $\boldsymbol {r}$, we use the pseudo-spectral algorithm of Patterson & Orszag (Reference Patterson and Orszag1971) with one phase shift and spherical truncation. We apply this algorithm to $\delta u_j$ and $\partial \delta u_i/\partial r_j$ for the interscale transfer and for $(u_j^++u_j^-)/2$ and $\partial \delta u_i/\partial x_j$ for the interspace transfer. Hence, we express these vectors/tensors in Fourier space (see (B13)–(B16) in Appendix B) and apply the pseudo-spectral method of Patterson & Orszag (Reference Patterson and Orszag1971) to calculate $\widehat {\boldsymbol {a}_{\mathcal {T}}}(\boldsymbol {k},\boldsymbol {r},t)$ and $\widehat {\boldsymbol {a}_{{\varPi }}}(\boldsymbol {k},\boldsymbol {r},t)$ without aliasing errors. We next decompose these fields to irrotational and solenoidal fields with the projection operator and inverse these fields to physical space to obtain $\boldsymbol {a}_{{\varPi }_{\bar {S}}}(\boldsymbol {x},\boldsymbol {r},t)$, $\boldsymbol {a}_{{\varPi }_{\bar {I}}}(\boldsymbol {x},\boldsymbol {r},t)$, $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}(\boldsymbol {x},\boldsymbol {r},t)$ and $\boldsymbol {a}_{\mathcal {T}_{\bar {I}}}(\boldsymbol {x},\boldsymbol {r},t)$. These fields can then be sampled over $\boldsymbol {x}$ to calculate, e.g. $\boldsymbol {a}_{{\varPi }_{\bar {S}}}^2(\boldsymbol {x},\boldsymbol {r},t)$ or KHMH terms such as $2\delta \boldsymbol {u} \boldsymbol{\cdot} \boldsymbol {a}_{{\varPi }_{\bar {S}}}(\boldsymbol {x},\boldsymbol {r},t)$ (see § 3.4). If we assume that the cost of a DNS time step is similar to the cost of the pseudo-spectral method to calculate the NS nonlinear term, the calculation of solenoidal and irrotational interspace and interscale transfers for one $t$ and one $\boldsymbol {r}$ has similar cost to one DNS time step. The total cost of the pseudo-spectral post-processing method is proportional to the total number $N_r$ of separation vectors $\boldsymbol {r}$ that we use in our spherical averaging across scales $r_d$ and to the total number $T_s/\Delta T$ of samples in time (see table 1). With a total number of separation vectors $N_r \sim 10^{3}\unicode{x2013}10^4$ and our $T_s/\Delta T$ values, the total cost of the pseudo-spectral post-processing method in terms of DNS time steps is at least one order of magnitude larger than the cost of the DNS itself. This high post-processing cost limits the possible number of grid points in this study. If we estimate the wall time of a $1024^{3}$ simulation to be approximately 10 days, the post-processing would require approximately three to four months.

The NSD solenoidal equation (3.21) describes a balance between the time derivative, solenoidal interscale transfer, solenoidal interspace transport, viscous and forcing terms. From the point we made in the sentence directly following (3.19), we expect $\langle \boldsymbol {a}_{\mathcal {T}_{\bar {S}}}^{2}\rangle$ and $\langle \boldsymbol {a}_{{\varPi }_{\bar {S}}}^{2}\rangle$ to tend to become equal to each other as the amplitude of $\boldsymbol {r}$ tends to values significantly larger than $\langle L \rangle _t$. Figure 4 confirms this trend for the second-order orientation-averaged moments of $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$. For brevity, in what follows we refer to such statistics as second-order magnitudes. With decreasing $r_d$, the magnitudes of $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ decrease relative to those of $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$. At all scales $r_{d}\geq \langle \lambda \rangle _t$ the second-order magnitudes of $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ are one order of magnitude larger than those of the viscous term $\delta \boldsymbol {a}_{\nu }$ and this separation is greater for the larger $\langle {{Re}}_\lambda \rangle _t$. The magnitudes of $\delta \boldsymbol {a}_{\nu }$ are themselves much larger than those of $\delta \boldsymbol {f}$ (not shown in figure 4 for not overloading the figure but see figure 3a1). These observations suggest that the solenoidal NSD equation (3.21) reduces to the approximate

(3.23)\begin{equation} \delta \boldsymbol{a}_l + \boldsymbol{a}_{\mathcal{T}_{\bar{S}}} \approx-\boldsymbol{a}_{{\varPi}_{\bar{S}}} , \end{equation}

where this equation is understood as typical in terms of second-order magnitudes: i.e. in most regions of the flow for the majority of the time, the removed terms are at least one order of magnitude smaller than the retained terms. (As for the NS dynamics, we do expect dynamically important regions localised in space and time where the dynamics differ from (3.23).) Figure 4 confirms (3.23) in a second-order sense and shows that the relatively rare spatio-temporal events that are neglected when writing (3.23) are indeed present as the second-order statistics do show a very small deviation from (3.23). An additional important observation to be made from figure 4 is that $\delta \boldsymbol {a}_{c_S}$ tends to become increasingly dominated by $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ rather than $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ as $r_d$ decreases.

Equation (3.23) is the same as (3.12), and similarly to figure 3 that provides support for (3.12), figures 4 and 5 provide strong support for (3.23), in particular for $r_{d} > \langle \lambda \rangle _{t}$. It is interesting to note that the average alignment between the left- and right-hand sides of (3.23) lies between 90 % and 100 % (typically 95 %) for $r_{d} > \langle \lambda \rangle _{t}$. Whilst this is strong support for approximate (3.23), the fact that the alignment is not 100 % is a reminder of the nature of the approximation, i.e. that relatively rare spatio-temporal events do exist where the viscous and/or driving forces are not negligible.

At length scales $r_{d} \le \langle \lambda \rangle _{t}$, the alignment between $\delta \boldsymbol {a}_l$ and $-\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ improves while the alignment between $\delta \boldsymbol {a}_l + \boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $-\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ worsens with decreasing $r_d$ (see figure 5), presumably because of direct dissipation and diffusion effects, so that $\delta \boldsymbol {a}_l + \boldsymbol {a}_{\mathcal {T}_{\bar {S}}} \approx 0$ becomes a better approximation than (3.23) at $r_d < 0.5 \langle \lambda \rangle _{t}$. This observation is consistent with our parallel observation that the magnitude of $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ increases while the magnitude of $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ decreases with decreasing $r_d$ and that $\delta \boldsymbol {a}_{c_S}$ in (3.12) tends to be dominated by $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ at the very smallest scales.

On the other end of the spectrum, i.e. as the length scale $r_d$ grows towards $\langle L \rangle _{t}$, the alignment between $\delta \boldsymbol {a}_l$ and $-\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ worsens while the alignment between $\delta \boldsymbol {a}_l$ and $-\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ improves (see figure 5), both reaching a comparable level of alignment/misalignment that contribute together to keep approximation (3.23) statistically well satisfied with 95 % alignment between $\delta \boldsymbol {a}_l + \boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $-\boldsymbol {a}_{{\varPi }_{\bar {S}}}$.

The strong anti-alignment between $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\delta \boldsymbol {a}_l$, increasingly so at smaller $r_d$ (see figure 5) expresses the sweeping of the two-point momentum difference $\delta \boldsymbol {u}$ at scales $r_d$ and smaller by the mainly large scale velocity $(\boldsymbol {u}^++\boldsymbol {u}^-)/2$. Note that this two-point sweeping differs from anti-alignment between $\delta \boldsymbol {a}_l$ and $\delta \boldsymbol {a}_c$ for two reasons. Firstly, by using the Helmholtz decomposition we have removed the pressure effect embodied in the $\boldsymbol {a}_{c_I}$ contribution to $\boldsymbol {a}_c$ that balances the pressure gradient. This was first understood in Tsinober et al. (Reference Tsinober, Vedula and Yeung2001) in a one-point setting and is here extended to a two-point setting. Secondly, $\delta \boldsymbol {a}_{c_S}$ is the sum of an interspace transport $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and an interscale transfer term $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ such that the interpretation of two-point sweeping as anti-alignment between $\boldsymbol {a}_{c_S}$ and $\boldsymbol {a}_l$ as sweeping cannot be exactly accurate. The advection of $\delta \boldsymbol {u}$ by the large-scale velocity is attributable to $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$, and figure 5 shows that the two-point sweeping anti-alignment between $\delta \boldsymbol {a}_l$ and $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ increases with decreasing $r_d$.

The sweeping anti-alignment between $\delta \boldsymbol {a}_l$ and $\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ is by no means perfect even if it reaches about 90 % accuracy at $r_{d} < \langle \lambda \rangle _{t}$, as is clear from the similar magnitudes and very strong alignment tendency between $\delta \boldsymbol {a}_l + \boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $-\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ at scales $|\boldsymbol {r}|\geq \langle \lambda \rangle _t$ (see figures 4 and 5). Note, in passing, that the Lagrangian solenoidal acceleration $\delta \boldsymbol {a}_l+\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ are both Galilean invariant. Equation (3.23) may be interpreted to mean that the Lagrangian solenoidal acceleration of $\delta \boldsymbol {u}$ (which is actually solenoidal) moving with the mainly large-scale velocity $(\boldsymbol {u}^++\boldsymbol {u}^-)/2$, namely $\delta \boldsymbol {a}_l+\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$, is evolving in time and space in response to $-\boldsymbol {a}_{{\varPi }_{\bar {S}}}$: when there is an influx of momentum from larger scales, there is an increase in $\delta \boldsymbol {a}_l+\boldsymbol {a}_{\mathcal {T}_{\bar {S}}}$ and $\delta \boldsymbol {u}$ and it vice versa.

3.4. From NSD dynamics to KHMH dynamics in homogeneous/periodic turbulence

The scale-by-scale evolution of $|\delta \boldsymbol {u}|^2$ locally in space and time is governed by a KHMH equation. This makes KHMH equations crucial tools for examining the turbulent energy cascade. The original KHMH equation and the new solenoidal and irrotational KHMH equations that we derive below are simply projections of the corresponding NSD equations onto $2 \delta \boldsymbol {u}$. Hence, KHMH dynamics depend on NSD dynamics and the various NSD terms’ alignment or non-alignment tendencies with $2 \delta \boldsymbol {u}$. In this subsection we present five KHMH results.

By contracting the NSD equation (3.9) with $2\delta \boldsymbol {u}$, one obtains the KHMH equation (Hill Reference Hill2002; Yasuda & Vassilicos Reference Yasuda and Vassilicos2018)

(3.24)\begin{align} & \frac{\partial}{\partial t}|\delta \boldsymbol{u}|^2+\frac{u_k^++u_k^-}{2} \frac{\partial}{\partial x_k}|\delta \boldsymbol{u}|^2 + \frac{\partial}{\partial r_k}\left(\delta u_k|\delta \boldsymbol{u}|^2\right) =-\frac{2}{\rho}\frac{\partial}{\partial x_k}\left(\delta u_k \delta p \right) + 2 \nu \frac{\partial^2}{\partial r_k^2}|\delta \boldsymbol{u}|^2 \nonumber\\ &\quad + \frac{\nu}{2}\frac{\partial^2}{\partial x_k^2}|\delta \boldsymbol{u}|^2 -\left[2 \nu\left(\frac{\partial u_i^+}{\partial x_k^+} \right)^2+ 2 \nu \left(\frac{\partial u_i^-}{\partial x_k^-} \right)^2 \right] +2\delta u_k \delta f_k, \end{align}

where no fluid velocity decomposition nor averaging operations have been used. In line with the naming convention of Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) this equation can be written as

(3.25)\begin{equation} \mathcal{A}_t+\mathcal{T}+{\varPi} = \mathcal{T}_p+\mathcal{D}_{r,\nu}+ \mathcal{D}_{x,\nu}-\mathcal{\epsilon}+\mathcal{I}, \end{equation}

where the first, second and third terms on the left-hand sides of (3.24) and (3.25) correspond to each other and so do the first, second, third, fourth and fifth terms on the right-hand sides. Pre-empting notation used further down in this paper, (3.25) is also written as $\mathcal {A} = \mathcal {T}_p + \mathcal {D} + \mathcal {I}$ or $\mathcal {A}_t + \mathcal {A}_c = \mathcal {T}_p + \mathcal {D} + \mathcal {I}$, where $\mathcal {A}_c \equiv \mathcal {T}+{\varPi }$, $\mathcal {A}\equiv \mathcal {A}_t+\mathcal {A}_c$ and $\mathcal {D}\equiv \mathcal {D}_{r,\nu }+\mathcal {D}_{x, \nu }-\mathcal {\epsilon }$.

To examine the KHMH dynamics in terms of irrotational and solenoidal dynamics, we contract the irrotational and solenoidal NSD equations with $2 \delta \boldsymbol {u}$ to derive what we refer to as irrotational and solenoidal KHMH equations. Each of the KHMH terms can be subdivided into a contribution from the NSD irrotational part and a contribution from the NSD solenoidal part of the respective term in the NSD equation. A solenoidal KHMH term corresponding to a $\delta \boldsymbol {q} (\boldsymbol {x},\boldsymbol {r},t)$ or $\boldsymbol {q} (\boldsymbol {x},\boldsymbol {r},t)$ term in (3.21) equals $\mathcal {Q}_{\bar {S}}=2\delta \boldsymbol {u} \boldsymbol{\cdot} \delta \boldsymbol {q}_{\bar {S}}$ or $\mathcal {Q}_{\bar {S}}=2\delta \boldsymbol {u} \boldsymbol{\cdot} \boldsymbol {q}_{\bar {S}}$, and an irrotational KHMH term corresponding to a $\delta \boldsymbol {q} (\boldsymbol {x},\boldsymbol {r},t)$ or $\boldsymbol {q} (\boldsymbol {x},\boldsymbol {r},t)$ term in (3.22) equals $\mathcal {Q}_{\bar {I}}=2\delta \boldsymbol {u} \boldsymbol{\cdot} \delta \boldsymbol {q}_{\bar {I}}$ or $\mathcal {Q}_{\bar {I}}=2\delta \boldsymbol {u} \boldsymbol{\cdot} \boldsymbol {q}_{\bar {I}}$. With $\mathcal {Q} = 2\delta \boldsymbol {u} \boldsymbol{\cdot} \delta \boldsymbol {q}$ or $\mathcal {Q} = 2\delta \boldsymbol {u} \boldsymbol{\cdot} \boldsymbol {q}$, we have $\mathcal {Q}= \mathcal {Q}_{\bar {I}} + \mathcal {Q}_{\bar {S}}$. The irrotational and solenoidal KHMH equations for periodic/homogeneous turbulence follow from (3.21) and (3.22), respectively, and read

(3.26)$$\begin{gather} \mathcal{A}_t +\mathcal{T}_{\bar{S}} + {\varPi}_{\bar{S}} = \mathcal{D}_{r,\nu} + \mathcal{D}_{x,\nu} - \mathcal{\epsilon} + \mathcal{I}, \end{gather}$$

(3.27)$$\begin{gather}{\varPi}_{\bar{I}} = \mathcal{T}_{\bar{I}} = \tfrac{1}{2}\mathcal{T}_p, \end{gather}$$

where use has been made of the fact that the velocity and velocity difference fields are solenoidal. These two equations are our first KHMH result.

Space-local changes in time of $|\delta \boldsymbol {u}|^2$, expressed via $\mathcal {A}_t$, are only due to solenoidal KHMH dynamics in (3.26) that include interspace transport, interscale transport, viscous and forcing effects. The irrotational KHMH equation (3.27) formulates how the imposition of incompressibility by the pressure field affects interspace and interscale dynamics and, in turn, energy cascade dynamics. Generalized solenoidal and irrotational KHMH equations also valid for non-periodic/non-homogeneous turbulence are given in Appendix B.

We first consider the spatio-temporal average of these equations in statistically steady forced periodic/homogeneous turbulence. As $\langle \mathcal {T}_p \rangle =0$, we obtain from (3.27), $\langle {\varPi }_{\bar {I}} \rangle = \langle \mathcal {T}_{\bar {I}} \rangle =0$. As $\langle \mathcal {T}_{\bar {S}}\rangle +\langle \mathcal {T}_{\bar {I}}\rangle =\langle \mathcal {T}\rangle =0$, we have $\langle \mathcal {T}_{\bar {S}}\rangle =0$, such that the spatio-temporal average of (3.26) reads

(3.28)\begin{equation} \langle {\varPi}\rangle = \langle {\varPi}_{\bar{S}}\rangle = \langle \mathcal{D}_{r,\nu} \rangle - \langle \mathcal{\epsilon} \rangle + \langle \mathcal{I} \rangle . \end{equation}

If an intermediate inertial subrange of scales $\vert \boldsymbol {r}\vert$ can be defined where viscous diffusion and forcing are negligible, (3.28) reduces to $\langle {\varPi }_{\bar {S}}\rangle \approx - \langle \mathcal {\epsilon } \rangle$ in that range. This theoretical conclusion (which is not part of our DNS study) is the backbone of the Kolmogorov (Reference Kolmogorov1941a,Reference Kolmogorovb,Reference Kolmogorovc) theory for a high Reynolds number statistically homogeneous stationary small-scale turbulence with the additional information that the part of the average interscale transfer rate involved in Kolmogorov's equilibrium balance is the solenoidal interscale transfer rate only. This is our second KHMH result. On average, there is a cascade of turbulence energy from large to small scales where the rate of interscale transfer is dominated by two-point vortex stretching (see Appendix C for the relation between the solenoidal interscale transfer and vortex stretching) and is equal to $- \langle \mathcal {\epsilon } \rangle$ independently of $\vert \boldsymbol {r}\vert$ over a range of scales where viscous diffusion and forcing are negligible.

In this paper we concentrate on the fluctuations around the average picture described by the scale-by-scale equilibrium (3.28) for any Reynolds number. If we subtract the spatio-temporal average solenoidal KHMH equation (3.28) from the solenoidal KHMH equation (3.26) and use the generic notation $\mathcal {Q}^{\prime}\equiv \mathcal {Q}-\langle \mathcal {Q}\rangle$, we attain the fluctuating solenoidal KHMH equation

(3.29)\begin{equation} \mathcal{A}_t +\mathcal{T}_{\bar{S}} +{\varPi}_{\bar{S}}^{\prime} = \mathcal{D}_{r,\nu}^{\prime} + \mathcal{D}_{x,\nu} - \mathcal{\epsilon}^{\prime} +\mathcal{I}^{\prime}. \end{equation}

This equation governs the fluctuations of the KHMH solenoidal dynamics around its spatio-temporal average. Clearly, if these non-equilibrium fluctuations are large relative to their average values, the average picture expressed by (3.28) is not characteristic of the interscale transfer dynamics. We now study the KHMH fluctuations in statistically stationary periodic/homogeneous turbulence on the basis of (3.27) and (3.29). Concerning (3.27), note that ${\varPi }_{\bar {I}}^{\prime}= {\varPi }_{\bar {I}}$, $\mathcal {T}_{\bar {I}}^{\prime}= \mathcal {T}_{\bar {I}}$ and $\mathcal {T}_{p}^{\prime} = \mathcal {T}_{\bar {p}}$.

We start by determining the relative fluctuation magnitudes of the spatio-temporal fluctuations of each term in the KHMH equations (3.27) and (3.29). These relative fluctuation magnitudes can emulate those of respective terms in the NSD equations under the following sufficient conditions: (i) the fluctuations are so intense that they dwarf averages, so that $\langle (\mathcal {Q}\prime)^2\rangle \approx \langle \mathcal {Q}^2\rangle$; (ii) the mean square of any KHMH term $\mathcal {Q} = 2\delta \boldsymbol {u} \boldsymbol{\cdot} \delta \boldsymbol {q}$ corresponding to a NSD term $\delta \boldsymbol {q} (\boldsymbol {x},\boldsymbol {r},t)$ (equivalently $\mathcal {Q} = 2\delta \boldsymbol {u} \boldsymbol{\cdot} \boldsymbol {q}$ corresponding to $\boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t)$) can be approximated as

(3.30)\begin{equation} {\langle \mathcal{Q}^2 \rangle(\boldsymbol{r}) \approx 4\langle |\delta \boldsymbol{u}|^2 \rangle \langle | \delta \boldsymbol{q}|^2 \rangle \langle \cos^2(\theta_{q}) \rangle,}\end{equation}

where the approximate equality results from a degree of decorrelation and $\theta _q$ is the angle between $\delta \boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t)$ (or $\boldsymbol {q}(\boldsymbol {x},\boldsymbol {r},t)$) and $\delta \boldsymbol {u}(\boldsymbol {x},\boldsymbol {r},t)$; (iii) $\langle \cos ^2(\theta _{q}) \rangle$ is not very sensitive to the choice of NSD term $\delta \boldsymbol {q}$ (or $\boldsymbol {q}$). Under these conditions, we get

(3.31)\begin{equation} {\frac{\langle (2 \delta \boldsymbol{u} \boldsymbol{\cdot} \delta \boldsymbol{q})^2\rangle (\boldsymbol{r})}{\langle (2\delta \boldsymbol{u} \boldsymbol{\cdot} \delta \boldsymbol{w})^2 \rangle (\boldsymbol{r})} \approx \frac{\langle |\delta \boldsymbol{u}|^2\rangle \langle | \delta \boldsymbol{q}|^2\rangle \langle \cos^2(\theta_q) \rangle (\boldsymbol{r})}{\langle|\delta \boldsymbol{u}|^2\rangle \langle |\delta \boldsymbol{w}|^2 \rangle \langle \cos^2(\theta_w)\rangle (\boldsymbol{r})} \approx \frac{\langle|\delta \boldsymbol{q}|^2\rangle (\boldsymbol{r})}{\langle|\delta \boldsymbol{w}|^2\rangle (\boldsymbol{r})}}, \end{equation}

which means that KHMH relative fluctuation magnitudes and NSD relative fluctuation magnitudes are approximately identical. The first approximate equality in (3.31) follows directly from (3.30) and the second approximate equality follows from hypothesis (iii) that $\cos ^2(\theta _q)$ and $\cos ^2(\theta _w)$ are about equal.

We test hypothesis (i) by comparing the plots in figure 6(a1,b1) with those in figure 6(a2,b2). Figure 6(a1,b1) shows average magnitudes of KHMH spatio-temporal fluctuations for terms with non-zero spatio-temporal averages. Comparing with figure 6(a2,b2), we find that $\langle (\mathcal {Q}\prime)^2\rangle ^a\approx \langle \mathcal {Q}^2\rangle ^a$, i.e. hypothesis (i), for all four terms plotted in figure 6(a1,b1) at all length scales $r_d$ considered. Note that this does not hold for $\mathcal {D}_{r,\nu }^{\prime}$ and $\mathcal {I}^{\prime}$ that are the only KHMH fluctuations such that $\sqrt {\langle (\mathcal {Q}\prime)^2\rangle ^a}/\langle \mathcal {\epsilon }\rangle ^a$ is smaller (in fact significantly smaller) than $1$ at all scales. Figure 6 also makes it clear that the magnitudes of the fluctuations of all other KHMH terms (solenoidal and irrotational) are much higher than those of the turbulence dissipation at all scales $r_d > 0.5 \langle \lambda \rangle _{t}$, and more so for the higher of the two Reynolds numbers. For scales $r_d \geq \langle \lambda \rangle _t$, the largest average fluctuating magnitudes are those of $\mathcal {A}_c^{\prime}$, followed closely by $\mathcal {A}_t$ and $\mathcal {T}_{\bar {S}}$. Next come the magnitudes of ${\varPi }_{\bar {S}}^{\prime}$ and $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$. Thereafter follow the irrotational terms ${\varPi }_{\bar {I}}=\mathcal {T}_{\bar {I}}$ $(=0.5 \mathcal {T}_{p})$ and finally the viscous, dissipative and forcing terms $\mathcal {D}^{\prime}, \mathcal {\epsilon }^{\prime}$ and $\mathcal {I}^{\prime}$ in that order. This order of fluctuations is our third KHMH result. An average description of the interscale turbulent energy transfer dynamics in terms of its spatio-temporal average cannot, therefore, be accurate. In order to characterize these dynamics, attention must be directed at most if not all KHMH term fluctuations, and in fact to much more than just the turbulence dissipation fluctuations given that they are among the weakest.

Figure 6. (a1,b1) The KHMH average square magnitudes $\langle \mathcal {Q}^2\rangle ^a$ and (a2,b2) KHMH average square fluctuating magnitudes $\langle (\mathcal {Q}^{\prime})^2\rangle ^a$, where $\mathcal {Q}^{\prime}=\mathcal {Q}-\langle \mathcal {Q}\rangle$, for the KHMH terms $\mathcal {Q}$ listed above the figures and introduced in the third and fourth paragraphs of § 3.4. All entries are normalised with $\langle \mathcal {\epsilon } \rangle ^a$ (see (3.24) and (3.25)). The following pairs of KHMH terms have overlapping magnitudes in (a2,b2): $\mathcal {A}_t$ and $\mathcal {A}_{c_S}$; $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ and ${\varPi }_{\bar {S}}$; $\mathcal {T}_{\bar {I}}$ and ${\varPi }_{\bar {I}}$. Results are shown for (a1,a2) $\langle {{Re}}_\lambda \rangle _t = 112$, (b1,b2) $\langle {{Re}}_\lambda \rangle _t=174$.

Next, we test hypothesis (ii) by testing the validity of (3.30) and hypothesis (iii) concerning approximately similar $\cos ^{2}(\theta _{q})$ behaviour for different KHMH terms. In figure 7(a1,b1) we plot ratios of the right-hand sides to left-hand sides of (3.30) and see that (3.30) is not valid, but that it is nevertheless about 65 % to 98 % accurate for $r_{d} \ge \langle \lambda \rangle _t$. Note that (3.30) might be sufficient but that it is by no means necessary for the leftmost and rightmost sides of (3.31) to approximately balance. In those cases where the variations between the ratios plotted in figure 7(a1,b1) are not too large and the assumption of approximately similar $\cos ^{2}(\theta _{q})$ for different KHMH terms more or less holds, the leftmost and rightmost sides of (3.31) can approximately balance.

Figure 7. Test of the assumptions (ii) and (iii) in the seventh paragraph of § 3.4 related to relations (3.30) and (3.31) connecting NSD and KHMH relative magnitudes. (a1,b1) Test of assumption (ii) by taking the ratio of the left-hand and right-hand sides of (3.30) for the KHMH terms $\mathcal {Q}$ listed above the figures. ($a2$,$b2$) Test of assumption (iii) used in (3.31) by comparing the behaviour of $\langle \cos ^2 (\theta _q)\rangle ^a$ for the various NSD terms listed above the figures. The black horizontal line $0.5$ corresponds to the value of $\langle \cos ^2 (\theta _q)\rangle$ if $\theta _q$ is uniformly distributed. Results are shown for (a1,a2) $\langle {{Re}}_\lambda \rangle _t = 112$, (b1,b2) $\langle {{Re}}_\lambda \rangle _t=174$.

Incidentally, figure 7(a2,b2) also shows that the angles $\theta _q$ are not random but that they are more likely to be small rather than large in an approximately similar way for all important NSD terms: $\cos ^2(\theta _{q})$ ranges between about 0.28 and 0.36 for all NSD terms (except the viscous acceleration difference and the viscous force difference) at all scales $r_d$. These values are much smaller than 0.5, the value that $\cos ^2(\theta _{q})$ would have taken if the angles $\theta _q$ were random. There is therefore an alignment tendency between $\delta \boldsymbol {u}$ and NSD terms that is similar for all the important NSD terms, thereby allowing the balance between the leftmost (ratio of KHMH terms) and the rightmost (ratio of NSD terms) sides of (3.31) to approximately hold as seen by comparing the plots (a1)-(b1) (mean square NSD terms) with the plots (a2)-(b2) (mean square KHMH terms) in figure 8. (Note that the viscous term is bounded from above, $\langle \mathcal {D}^2\rangle (\boldsymbol {r}) \leq 4 \langle |\delta \boldsymbol {u}|^2 |\delta \boldsymbol {a}_{\nu }|^2 \rangle$, which indicates limited magnitudes compared with the irrotational and dominant solenoidal terms because of the limited magnitude of $\langle \delta \boldsymbol {a}_{\nu }^2 \rangle$. The limited fluctuations of the viscous terms are clearly seen in figure 6.)

Figure 8. The NSD and KHMH relative average square magnitudes (which should be similar on the basis of (3.31)) for the terms listed above the figures: (a1) NSD and (a2) KHMH for $\langle {{Re}}_\lambda \rangle _t = 112$, (b1) NSD and (b2) KHMH for $\langle {{Re}}_\lambda \rangle _t=174$.

Figure 8 does indeed confirm the close correspondence between NSD and KHMH statistics that is a significant step further from the correspondence reported earlier in this paper between NS and NSD statistics. We can therefore use the approximate NSD relation (3.23) to deduce the following approximate KHMH relation:

(3.32)\begin{equation} \mathcal{A}_t + \mathcal{T}_{\bar{S}} + {\varPi}_{\bar{S}}^{\prime} \approx 0, \end{equation}

understood in the sense that it holds in the majority of the domain for the majority of the time but that there surely exist relatively rare events within the flow where this approximate KHMH relation is violated.

This approximate equation $\mathcal {A}_t + \mathcal {T}_{\bar {S}} + {\varPi }_{\bar {S}}^{\prime} \approx 0$ can be considered to be our fourth KHMH result. It is consistent with the order of fluctuation magnitudes in figure 8 that shows, in agreement with the NSD–KHMH correspondence just established, that the largest fluctuating magnitudes are those of $\mathcal {A}_c$, followed by the fluctuating magnitudes of $\mathcal {T}_{\bar {S}}$, $\mathcal {A}_t$ and $\mathcal {A}_{c_S}$ ($\mathcal {A}_{c_S} = \mathcal {T}_{\bar {S}} + {\varPi }_{\bar {S}}$). Note though that there is a crossover at about $r_{d} \approx 2 \langle \lambda \rangle _{t}$ for both Reynolds numbers considered here between the fluctuation magnitudes of $\mathcal {T}_{\bar {S}}$ and those of $\mathcal {A}_t$ and $\mathcal {A}_{c_S}$ that are about equal to each other in agreement with (3.32).

The fluctuation magnitudes of ${\varPi }_{\bar {S}}$ and ${\varPi }_{\bar {I}}$ are both smaller than those just mentioned, and those of ${\varPi }_{\bar {I}}$ are significantly smaller than those of ${\varPi }_{\bar {S}}$. Even smaller are the fluctuation magnitudes of $\mathcal {D}$ and $\mathcal {I}$, in that order. In agreement with (3.16), our third and fourth KHMH conclusions incorporate

(3.33)\begin{equation} \langle \mathcal{A}_{t}^{2} \rangle \approx \langle \mathcal{A}_{c_S}^{2} \rangle \gg \langle \mathcal{T}_{p}^{2}\rangle = 4 \langle {\varPi}_{\bar{I}}^{2}\rangle = 4 \langle \mathcal{T}_{\bar{I}}^{2}\rangle = \langle \mathcal{A}_{c_I}^{2}\rangle \gg \langle \mathcal{D}^{2}\rangle \gg \langle \mathcal{I}^{2}\rangle, \end{equation}

where $\mathcal {A}_{c_I} = \mathcal {T}_{\bar {I}} + {\varPi }_{\bar {I}}$.

An additional significant observation from figure 8 that we can count as our fifth KHMH result is that, as $r_{d}$ decreases towards about $0.5 \langle \lambda \rangle _{t}$, the fluctuation magnitude of $\mathcal {A}_{c_S} = \mathcal {T}_{\bar {S}} + {\varPi }_{\bar {S}}$ remains about constant but that of $\mathcal {T}_{\bar {S}}$ increases while that of ${\varPi }_{\bar {S}}$ decreases. (At scale $r_d$ smaller than $0.5 \langle \lambda \rangle _{t}$, the fluctuation magnitudes of both $\mathcal {A}_{c_S}$ and $\mathcal {T}_{\bar {S}}$ increase with diminishing $r_d$ whereas those of ${\varPi }_{\bar {S}}$ remain about constant.) The convective nonlinearity is increasingly of the spatial transport type and diminishingly of the interscale transfer type as the two-point separation length decreases.

We now consider correlations between different intermediate and large-scale fluctuating KHMH terms in light of (3.27) and (3.32).

4.1. Correlations

We start this section by assessing the existence or non-existence of local (in space and time) equilibrium between interscale transfer and dissipation at some intermediate scales. In figure 9 we plot correlations between various KHMH terms. In particular, this figure shows that the correlation coefficient between ${\varPi }_{\bar {S}}^{\prime}$ and $-\mathcal {\epsilon }^{\prime}$ lies well below $0.1$ for all scales $r_d \geq \langle \lambda \rangle _t$. The scatter plots of these quantities in figure 10 confirm the absence of a local relation between interscale transfer rate and dissipation rate. For example, for a given local/instantaneous dissipation fluctuation, the corresponding local/instantaneous interscale transfer rate fluctuation can be close to equally positive or negative. There is no local equilibrium between these quantities as they fluctuate at scales $r_d \geq \langle \lambda \rangle _t$. Such a correlation should of course not necessarily be expected. However, as $r_d$ decreases below $\langle \lambda \rangle _t$, the correlations between ${\varPi }_{\bar {S}}^{\prime}$ and either $-\mathcal {\epsilon }^{\prime}$ or $\mathcal {D}^{\prime}$ increase up to values between about $0.3$ and about $0.5$. This increased correlation may suggest a feeble tendency towards local/instantaneous equilibrium between interscale transfer rate and dissipation rate at scales $r_d < \langle \lambda \rangle _t$. However, these scales are strongly affected by direct viscous processes and can therefore not be inertial range scales.

Figure 9. Spherically averaged correlation coefficients between KHMH terms $(\mathcal {Q}_1, \mathcal {Q}_2)$ listed above the plots (a,b). They are plotted as functions of scale $r_d$. Results are shown for (a) $\langle {{Re}}_\lambda \rangle _t = 112$, (b) $\langle {{Re}}_\lambda \rangle _t = 174$.

Figure 10. Scatter plots of ${\varPi }_{\bar {S}}^{\prime}$ and $-\mathcal {\epsilon }^{\prime}$ at random orientations $\boldsymbol {r}$ with $r_d/\langle \lambda \rangle _t=(1.45,3.1)$ for (a,b), $\sigma _{{\varPi }_{\bar {S}}}$ is the standard deviation of ${\varPi }_{\bar {S}}$ and $\langle {{Re}}_\lambda \rangle _t=174$.

Following the question of local/instantaneous equilibrium, we now look for local/instantaneous sweeping. Figure 9 shows strong anti-correlation between $\mathcal {A}_t$ and $\mathcal {T}_{\bar {S}}$, increasingly so as $r_d$ decreases from large to small scales. Along with the fifth KHMH result at the end of the previous section (that the fluctuation magnitudes of $\mathcal {A}_t$ and $\mathcal {T}_{\bar {S}}$ become increasingly comparable as $r_d$ decreases), this anti-correlation tendency suggests a tendency towards $\mathcal {A}_t + \mathcal {T}_{\bar {S}} \approx 0$ at decreasing scales in agreement with the concept of two-point sweeping introduced in § 3.2. In other words, the sweeping of $|\delta \boldsymbol {u}|^2$ by the mainly large-scale advection velocity $(\boldsymbol {u}^++\boldsymbol {u}^-)/2$ becomes increasingly strong with decreasing $r_d$. The scatter plots of $\mathcal {A}_t$ and $\mathcal {T}_{\bar {S}}$ in figure 11 make this local/instantaneous two-point sweeping tendency with decreasing $r_d$ very evident, but also indicate that significant values of positive or negative ${\varPi }_{\bar {S}}$ can cause increasing deviations from $\mathcal {A}_t + \mathcal {T}_{\bar {S}}\approx 0$ as $r_d$ increases. Note that $\mathcal {A}_t+\mathcal {T}_{\bar {S}} +{\varPi }_{\bar {S}} \approx 0$ as indicated by the correlation coefficients in figure 9 between $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ and $-{\varPi }_{\bar {S}}$ (which exceed $0.95$ for $r_d\geq \langle \lambda \rangle _t$ at our Reynolds numbers) and by their overlapping fluctuation magnitudes in figure 6(a2,b2). The fluctuations of ${\varPi }_{\bar {S}}$ increase in magnitude as $r_d$ increases and so do high values of ${\varPi }_{\bar {S}}$ too. The scatter plots in figure 11 highlight how the $5\,\%$ most negative ${\varPi }_{\bar {S}}$ events (values of ${\varPi }_{\bar {S}}$ for which the probability that ${\varPi }_{\bar {S}}$ is smaller than a negative value ${\varPi }_{\bar {S}_{0.05}}$ is 0.05) and the $5\,\%$ most positive ${\varPi }_{\bar {S}}$ events (values of ${\varPi }_{\bar {S}}$ for which the probability that ${\varPi }_{\bar {S}}$ is larger than a positive value ${\varPi }_{\bar {S}_{0.95}}$ is also 0.05) cause significant deviations from ‘perfect sweeping’ $\mathcal {A}_t=-\mathcal {T}_{\bar {S}}$, increasingly so for increasing $r_d$, in agreement with $\mathcal {A}_t+\mathcal {T}_{\bar {S}} +{\varPi }_{\bar {S}} \approx 0$.

Figure 11. Scatter plots of $\mathcal {A}_t$ and $\mathcal {T}_{\bar {S}}$ at random orientations $\boldsymbol {r}$ normalised by $\sigma _{\mathcal {A}_t}$ and $\sigma _{\mathcal {T}_{\bar {S}}}$, their respective standard deviations. Here ${\varPi }_{\bar {S}_{0.05}}$ is the value of ${\varPi }_{\bar {S}}$ at the respective $r_d$ for which $5\,\%$ of the samples are more negative than ${\varPi }_{\bar {S}_{0.05}}$ and ${\varPi }_{\bar {S}_{0.95}}$ is the value of ${\varPi }_{\bar {S}}$ for which $95\,\%$ of the samples are more positive than ${\varPi }_{\bar {S}_{0.95}}$. The events ${\varPi }_{\bar {S}}<{\varPi }_{\bar {S}_{0.05}}$ and ${\varPi }_{\bar {S}}>{\varPi }_{\bar {S}_{0.95}}$ are marked in red and green, respectively, while the remaining events are marked in blue. The red line marks $\mathcal {A}_t = -\mathcal {T}_{\bar {S}} - \langle {\varPi }_{\bar {S}} | {\varPi }_{\bar {S}}<{\varPi }_{\bar {S}_{0.05}} \rangle$, where $\langle {\varPi }_{\bar {S}} | {\varPi }_{\bar {S}}<{\varPi }_{\bar {S}_{0.05}} \rangle$ is the average value of ${\varPi }_{\bar {S}}$ conditioned on ${\varPi }_{\bar {S}}<{\varPi }_{\bar {S}_{0.05}}$. The green line marks $\mathcal {A}_t = -\mathcal {T}_{\bar {S}} - \langle {\varPi }_{\bar {S}} | {\varPi }_{\bar {S}}>{\varPi }_{\bar {S}_{0.95}} \rangle$ and the blue line marks $\mathcal {A}_t = -\mathcal {T}_{\bar {S}}$ (with all terms appropriately normalised with $\sigma _{\mathcal {A}_t}$ and $\sigma _{\mathcal {T}_{\bar {S}}}$). Here $r_d/\langle \lambda \rangle _t=(0.12,1.45,3.1,5.2)$ for (a–d) and $\langle {{Re}}_\lambda \rangle _t=174$.

The scatter plots in figure 12 show that it is only in relatively rare circumstances that $\mathcal {A}_t+\mathcal {T}_{\bar {S}} +{\varPi }_{\bar {S}} \approx 0$ is significantly inaccurate for scales $r_d\geq \langle \lambda \rangle _t$. Similarly to NSD dynamics, $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ can be viewed as a Lagrangian time rate of change of $|\delta \boldsymbol {u}|^2$ moving with $(\boldsymbol {u}^++\boldsymbol {u}^-)/2$. As more than average $|\delta \boldsymbol {u}|^2$ is cascaded from larger to smaller scales at a particular location $({\varPi }_{\bar {S}}^{\prime}<0)$, $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ increases; and as more than average $|\delta \boldsymbol {u}|^2$ is inverse cascaded from smaller to larger scales $({\varPi }_{\bar {S}}^{\prime} > 0)$, $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ decreases. Here ${\varPi }_{\bar {S}}^{\prime}$ is to a large extent determined by $\boldsymbol {a}_{{\varPi }_{\bar {S}}}$ that, as we show in Appendix C, is a non-local function in space of the vortex stretching and compression dynamics determining the two-point vorticity difference $\delta \boldsymbol {\omega }$.

Figure 12. Scatter plots of $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ and ${\varPi }_{\bar {S}}^{\prime}$ at random orientations $\boldsymbol {r}$. The residual $-b \equiv \mathcal {A}_t+\mathcal {T}_{\bar {S}}+{\varPi }_{\bar {S}}^{\prime}$ and the values $b_{0.05}$ and $b_{0.95}$ are defined analogously as for ${\varPi }_{\bar {S}_{0.05}}$ and ${\varPi }_{\bar {S}_{0.95}}$ in the previous figure. The events $b< b_{0.05}$ and $b>b_{0.95}$ are marked in red and green, respectively, while the remaining events are marked in blue. The red line marks $\mathcal {A}_t + \mathcal {T}_{\bar {S}} = - {\varPi }_{\bar {S}}^{\prime} - \langle b | b< b_{0.05} \rangle$, the green line $\mathcal {A}_t + \mathcal {T}_{\bar {S}} = - {\varPi }_{\bar {S}}^{\prime} - \langle b | b>b_{0.95} \rangle$ and the blue line $\mathcal {A}_t + \mathcal {T}_{\bar {S}} = -{\varPi }_{\bar {S}}^{\prime}$ (with all terms appropriately normalised with $\sigma _{{\varPi }_{\bar {S}}}$). Here $r_d/\langle \lambda \rangle _t=(0.12,1.45,3.1,5.2)$ for $(a,b,c,d)$ and $\langle {{Re}}_\lambda \rangle _t=174$.

A fairly complete way to summarise the details of the balance $\mathcal {A}_t+\mathcal {T}_{\bar {S}} +{\varPi }_{\bar {S}} \approx 0$ at scales $r_d\geq \langle \lambda \rangle _t$ is by noting that, as $r_d$ decreases towards $\langle \lambda \rangle _t$, (i) the fluctuation magnitude of $\mathcal {T}_{\bar {S}}$ tends to become comparable to that of $\mathcal {A}_t$ while that of ${\varPi }_{\bar {S}}$ decreases by comparison, (ii) the correlation coefficient between $\mathcal {A}_t$ and $-\mathcal {T}_{\bar {S}}$ increases towards $0.9$, and also (iii) (not mentioned till now but evident in figure 9) the correlation coefficient between $\mathcal {A}_t$ and $-{\varPi }_{\bar {S}}$ decreases towards values below $0.2$.

4.2. Conditional correlations

At scales $r_d$ below $\langle \lambda \rangle _t$, the relation $\mathcal {A}_t+\mathcal {T}_{\bar {S}} +{\varPi }_{\bar {S}} \approx 0$ becomes less accurate as the correlation coefficient between $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ and $-{\varPi }_{\bar {S}}$ drops from $0.95$ to $0.7$ with decreasing $r_d$, reflecting the increase of correlation between $\mathcal {\epsilon }$ and $-{\varPi }_{\bar {S}}$ and the even higher increase towards values close to $0.5$ of the correlation coefficient between $\mathcal {D}$ and ${\varPi }_{\bar {S}}$. This increase of correlation appears to reflect the impact of relatively rare yet intense local/instantaneous occurrences of interscale transfer rate as shown in figure 13 where we plot correlations conditional on relatively rare interscale events where the magnitudes of the spherically averaged interscale transfer rates are higher than 95 % of all interscale transfer rates of same sign (positive for backward and negative for forward transfer) in our overall spatio-temporal sample. This impact is highest at scales smaller than $\langle \lambda \rangle _t$ where the correlation coefficient conditioned on intense forward or backward interscale transfer rate events of $\pm {\varPi }_{\bar {S}}$ and either $\mathcal {\epsilon }$ or $\mathcal {D}$ can be as high as $0.7$ ($+{\varPi }_{\bar {S}}$ in the case of backward events and $-{\varPi }_{\bar {S}}$ in the case of forward events that causes significantly higher correlations between $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ and either $-\mathcal {\epsilon }$ or $\mathcal {D}$ in the case of backward events than in the case of forward events as seen in figure 13). However, the impact of such relatively rare events is also manifest at scales larger than $\langle \lambda \rangle _t$ (see figure 13) where the conditioned correlation coefficient is significantly higher than the unconditioned one in figure 9. Interestingly, conditioning on these relatively rare events does not change the correlation coefficients of $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ with $-{\varPi }_{\bar {S}}^{\prime}$ except at scales $r_d$ smaller than $\langle \lambda \rangle _t$ where, consistently with the increased conditioned correlations between $-{\varPi }_{\bar {S}}$ and $\mathcal {D}$, they are smaller than the unconditional correlation coefficients of $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ with $-{\varPi }_{\bar {S}}^{\prime}$, particularly at relatively rare forward interscale events where this conditional correlation drops to values close to $0.3$ at scales well below $\langle \lambda \rangle _t$.

Figure 13. (a) Correlation coefficients among the $5\,\%$ strongest spherically averaged backward interscale transfer events ${\varPi }_{\bar {S}}^a>{\varPi }_{\bar {S}_{0.95}}^a$ for KHMH terms $(\mathcal {Q}_1^a,\mathcal {Q}_2^a)$ listed on top of the figure. (b) Correlation coefficients among the $5\,\%$ strongest spherically averaged forward interscale transfer events ${\varPi }_{\bar {S}}^a<{\varPi }_{\bar {S}_{0.05}}^a$ for KHMH terms $(\mathcal {Q}_1^a,\mathcal {Q}_2^a)$ listed on top of the figure. Here $\langle {{Re}}_\lambda \rangle _t=112$. (Corresponding plots for $\langle {{Re}}_\lambda \rangle _t=174$ are omitted because they are very similar.)

Given that our relatively rare intense interscale transfer rates can be the seat of some correlation between ${\varPi }_{\bar {S}}$ and either $-\mathcal {\epsilon }$ or $\mathcal {D}$ particularly for $r_d < \langle \lambda \rangle _t$, and given that $\mathcal {A}_t+\mathcal {T}_{\bar {S}} \approx 0$ is a good approximation at scales smaller than $\langle \lambda \rangle _t$, do we have approximate two-point sweeping and approximate equilibrium ${\varPi }_{\bar {S}} \approx \mathcal {D}$ if we condition on relatively rare forward or backward interscale transfer rate events? In fact the conditional correlations between $\mathcal {A}_t$ and $-\mathcal {T}_{\bar {S}}$ are very high (close to and above 0.95) at all scales (see figure 13), higher than the corresponding unconditional correlations. However, the conditional averages of $\mathcal {A}_t$ and $-\mathcal {T}_{\bar {S}}$ shown in figure 14 are also significantly different at all scales, implying that these strong conditional correlations do not actually amount to two-point sweeping at relatively rare forward and backward events. Furthermore, if we condition on high negative/positive values of ${\varPi }_{\bar {S}}$, the averages of both $\mathcal {A}_t$ and $\mathcal {T}_{\bar {S}}$ are positive/negative (figure 14), even though these conditional averages do tend to $0$ as $r_d$ tends to $0$. This has two implications.

1. (i) It implies that, even though $\mathcal {A}_t$ and $-\mathcal {T}_{\bar {S}}$ are very well correlated at these relatively rare events, $\mathcal {A}_t + \mathcal {T}_{\bar {S}}$ fluctuates around a constant $C$ where $C>0$ if we condition the fluctuations on relatively rare negative ${\varPi }_{\bar {S}}$ but $C<0$ if we condition them on relatively rare positive ${\varPi }_{\bar {S}}$ ($C=0$ if we do not condition). This amounts to a systematic deviation on the average from two-point sweeping even though the strong correlation between the high magnitude fluctuations of $\mathcal {A}_t$ and $-\mathcal {T}_{\bar {S}}$ point at a tendency towards sweeping that is frustrated by the presence of the comparatively low non-zero local ${\varPi }_{\bar {S}}$. Given (3.29), the presence of this non-zero constant $C$ (clearly non-zero for all scales, and non-zero but tending towards zero as $r_d$ tends to $0$ well below $\langle \lambda \rangle _t$) means that the equilibrium ${\varPi }_{\bar {S}} \approx \mathcal {D}$ for scales smaller than $\langle \lambda \rangle _t$ does not hold either, even at scales smaller than $\langle \lambda \rangle _t$ where the conditional correlation between ${\varPi }_{\bar {S}}$ and $\mathcal {D}$ is significant. In fact, figure 14 shows that the conditional averages of ${\varPi }_{\bar {S}}^{\prime}$ are much larger than those of both $\mathcal {D}^{\prime}$ and $-\mathcal {\epsilon }^{\prime}$; they are much closer to those of $\mathcal {A}_t + \mathcal {T}_{\bar {S}}$.

2. (ii) The second implication of the conditional signs of $\mathcal {T}_{\bar {S}}$ is the existence of a relation between the conditional average of solenoidal interspace transfer rate $\mathcal {T}_{\bar {S}}$ and the solenoidal interspace transfer rate ${\varPi }_{\bar {S}}$ on which the average is conditioned: when one is positive/negative, the other is negative/positive, and we also find that their absolute magnitudes increase together (see figure 15). This is an observation that may prove important in the future for both subgrid scale modelling and the detailed study of the very smallest scales of turbulence fluctuations.

Figure 14. (a) Spatio-temporal averages of KHMH terms $\mathcal {Q}^a$ conditioned on the $5\,\%$ strongest spherically averaged backward (a) and forward (b) interscale transfer events. The KHMH terms are listed above figure (a) and $\langle {{Re}}_\lambda \rangle _t=112$. (Corresponding plots for $\langle {{Re}}_\lambda \rangle _t=174$ are omitted because they are very similar.)

Figure 15. Spatio-temporal averages of $\mathcal {T}_{\bar {S}}^a$ across scales $r_d$ conditioned on ${\varPi }_{\bar {S}}^a$ being within a certain range of ${\varPi }_{\bar {S}}^a$ values and we consider $20$ such ranges of increasing values of ${\varPi }_{\bar {S}}^a$: the $5\,\%$ smallest/most negative ${\varPi }_{\bar {S}}^a$, the $5\,\%$ to $10\,\%$ smallest/most negative ${\varPi }_{\bar {S}}^a$ values and so on until the $5\,\%$ largest/most positive ${\varPi }_{\bar {S}}^a$ values. Results are shown for (a) $\langle {{Re}} \rangle _t=112$, (b) $\langle {{Re}} \rangle _t=174$.

In conclusion, ${\varPi }_{\bar {S}}$ does not fluctuate with neither $-\mathcal {\epsilon }$ nor $\mathcal {D}$. Instead, ${\varPi }_{\bar {S}}$ and $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ fluctuate together at all scales, in particular scales larger than $\langle \lambda \rangle _t$, and even at relatively rare interscale transfer events. At scales smaller than $\langle \lambda \rangle _t$, we have a general tendency towards two-point sweeping if we do not condition on particular events. At our relatively rare interscale transfer events this correlation tendency (now conditional) is in fact amplified but there is nevertheless a systematic average deviation from two-point sweeping consistent with the absence of equilibrium ${\varPi }_{\bar {S}} \approx \mathcal {D}$ at these events. Finally, a relation exists between interspace and interscale transfer rates because the average interspace transfer rate conditioned on positive/negative values of interscale transfer rate is negative/positive. It must be stressed, however, that this relation does not imply an anti-correlation between interscale and interspace transport rates. The unconditioned correlation coefficients between $-{\varPi }_{\bar {S}}$ and $\mathcal {T}_{\bar {S}}$ are around $0.2$ (see figure 9), and we checked that this $0.2$ correlation does not change significantly if we condition on relatively rare intense occurrences of interscale transfer rate.

5.1. Average values and probability density functions

The decomposition ${\varPi } = {\varPi }_{\bar {I}} +{\varPi }_{\bar {S}}$ helped us distinguish between the solenoidal vortex stretching/compression and the pressure-related aspects of the interscale transfer. As recently shown by Alves Portela et al. (Reference Alves Portela, Papadakis and Vassilicos2020), the interscale transfer rate ${\varPi }$ can also be decomposed in a way that brings out the fact that it has a direct inhomogeneity contribution to it. This last part of the present study is an examination of the decomposition introduced by Alves Portela et al. (Reference Alves Portela, Papadakis and Vassilicos2020) that is ${\varPi }={\varPi }_{I}+{\varPi }_H$, where

(5.1)$$\begin{gather} {\varPi}_I = \frac{1}{2} \delta u_i \frac{\partial}{\partial x_i }(u_k^+u_k^+-u_k^-u_k^-) , \end{gather}$$

(5.2)$$\begin{gather}{\varPi}_H =-2 \delta u_i \frac{\partial}{\partial r_i} (u_k^-u_k^+). \end{gather}$$

Here ${\varPi }_I$ can be locally/instantaneously non-zero only in the presence of a local/instantaneous inhomogeneity. However, it averages to zero, i.e. $\langle {\varPi }_I \rangle =0$, in periodic/statistically homogeneous turbulence. Note that ${\varPi }={\varPi }_{I}={\varPi }_{H}=0$ at $\boldsymbol {r}=\boldsymbol {0}$. With $\boldsymbol {r}$-orientation averaging, the decomposition ${\varPi }^{a}={\varPi }_{I}^{a}+{\varPi }_{H}^{a}$ is unique in the sense that any potentially suitable (e.g. such that it equals $0$ at $\boldsymbol {r}=\boldsymbol {0}$) $\boldsymbol {x}$-gradient term added to ${\varPi }_{I}$ vanishes after $\boldsymbol {r}$-orientation averaging (see Alves Portela et al. Reference Alves Portela, Papadakis and Vassilicos2020).

An equivalent expression for ${\varPi }_I$ that immediately reveals where the decomposition ${\varPi } = {\varPi }_I + {\varPi }_H$ comes from is ${\varPi }_I = \delta u_i ({\partial }/{\partial r_i})(u_k^+u_k^++u_k^-u_k^-)$. Given that the total interscale transfer rate is ${\varPi } = \delta u_i ({\partial }/{\partial r_i})(\delta u_k \delta u_k)$, the ${\varPi }_I$ part of the interscale transfer concerns the transferred energy differences coming mostly from differences between velocity amplitudes, i.e. local/instantaneous inhomogeneities of ‘turbulence intensity’ in the flow; the ${\varPi }_H$ part of the interscale transfer concerns transferred energy differences coming mostly from differences between velocity orientations. Consistently with its link to local/instantaneous non-homogeneity, ${\varPi }_I$ can be written in the form (5.1) making it clear that ${\varPi }_I$ is zero where and when fluctuating velocity magnitudes are locally uniform.

In comparing the decompositions ${\varPi } = {\varPi }_{\bar {S}} +{\varPi }_{\bar {I}}$ and ${\varPi }={\varPi }_{I}+{\varPi }_H$, it is worth noting that ${\varPi }_{I}= {\varPi }_{I_{\bar {I}}}$ given that ${\varPi }_{I_{\bar {S}}} = 0$ from its centroid gradient form (see (5.1)). It therefore follows that

(5.3)$$\begin{gather} {\varPi}_{\bar{S}} = {\varPi}_{H_{\bar{S}}} , \end{gather}$$

(5.4)$$\begin{gather}{\varPi}_{\bar{I}} = {\varPi}_{I} + {\varPi}_{H_{\bar{I}}}. \end{gather}$$

The inhomogeneity-based interscale transfer rate influences only the irrotational part of the total interscale transfer rate whereas ${\varPi }_{H}$ influences both the irrotational and the solenoidal parts. As $\langle {\varPi }_I \rangle =0$ and $\langle {\varPi }_{\bar {I}} \rangle =0$, it follows that $\langle {\varPi }_{H_{\bar {I}}}\rangle = 0$. More to the point, $\langle {\varPi }_{\bar {S}} \rangle$ equals $\langle {\varPi }_{H_{\bar {S}}}\rangle$ and so (3.28) reduces to

(5.5)\begin{equation} \langle {\varPi}\rangle = \langle{\varPi}_{H_{\bar{S}}}\rangle = \langle \mathcal{D}_{r,\nu} \rangle - \langle \mathcal{\epsilon} \rangle + \langle \mathcal{I} \rangle . \end{equation}

The part of the interscale transfer rate that is present in the average interscale transfer/cascade dynamics is in fact ${\varPi }_{H_{\bar {S}}}$.

Given that the average interscale transfer is controlled by ${\varPi }_{H_{\bar {S}}} = {\varPi }_{\bar {S}}$, it is worth asking whether the well-known negative skewness of the probability density function (PDF) of ${\varPi }^{a}$ (see, e.g. Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) and references therein) is also present in the PDF of ${\varPi }_{\bar {S}}^{a}$ or/and whether it is spread across different terms of our two interscale transfer rate decompositions. In figure 16 we plot the PDFs of ${\varPi }^{a}$ and of the different $\boldsymbol {r}$-orientation-averaged terms in the decompositions of ${\varPi }$ that we use. It is clear that the PDFs of ${\varPi }$ and ${\varPi }_{\bar {S}}$ are nearly identical whilst the PDFs of ${\varPi }_H$ are different though also negatively skewed. The PDFs of ${\varPi }_{H_{\bar {I}}}$, ${\varPi }_{\bar {I}}$ and ${\varPi }_{I}$ are not significantly skewed. In figure 17 we plot the skewnes factors of the various interscale transfer terms as well as some other KHMH terms. The inhomogeneity interscale transfer ${\varPi }_{I}$ has close to zero skewness across scales. Both ${\varPi }_{\bar {S}}$ and ${\varPi }_{H}$ are negatively skewed, the former more so than the latter. Given (5.3)–(5.4) and ${\varPi }_H= {\varPi }_{\bar {S}}+{\varPi }_{H_{\bar {I}}}$, this difference in skewness factors is due to the irrotational part of ${\varPi }_H$ that is not significantly skewed and reduces the skewness of ${\varPi }_H$ relative to that of ${\varPi }_{\bar {S}}$. All in all, the skewness towards forward rather than inverse interscale transfers is present in its homogeneous and solenoidal components but is absent in its non-homogeneous and irrotational parts.

Figure 16. (a–f) The PDFs of ${\varPi }$ decompositions $({\varPi },{\varPi }_{H_{\bar {I}}},{\varPi }_{\bar {S}},{\varPi }_{\bar {I}},{\varPi }_H,{\varPi }_I)$ at $\langle {{Re}}_\lambda \rangle _t=112$. Here $\sigma _{{\varPi }^a}$ denotes the standard deviation of ${\varPi }^a$ and $P_{{max}}$ denotes the maximum value of the PDF of ${\varPi }^a$. The inhomogeneity and homogeneity interscale transfer rates ${\varPi }_I$ and ${\varPi }_H$ are defined in (5.1)–(5.2) and the irrotational part of the homogeneity interscale transfer rate ${\varPi }_{H_{\bar {I}}}$ in (5.4).

Figure 17. Skewness factors for KHMH terms $\mathcal {Q}$ listed on top of (a): (a) $\langle {{Re}}_\lambda \rangle _t=112$, (b) $\langle {{Re}}_\lambda \rangle _t=174$.

Figure 17 also shows that $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ is slightly positively skewed with flatness factors of approximately $0.5$ at scales $r_d \geq \langle \lambda \rangle _t$ and close to 0 or below at scales below $\langle \lambda \rangle _t$. The skewness factor of $-{\varPi }_{\bar {S}}$ with which $\mathcal {A}_t+\mathcal {T}_{\bar {S}}$ is very well correlated (as we have seen in the previous section) is about the same at scales close to the integral scale but steadily increases to values well above $0.5$ as $r$ decreases, reaching nearly $6.0$ at scales close to $0.5\langle \lambda \rangle _t$. This is a concrete illustration of the fact already mentioned earlier in this paper that $\mathcal {A}_t+\mathcal {T}_{\bar {S}} \approx -{\varPi }_{\bar {S}}$ is a very good approximation for most locations and most times but not all. Given the very significantly increased correlation/anti-correlation of ${\varPi }_{\bar {S}}$ with both $\mathcal {D}$ and $\mathcal {\epsilon }$ at relatively intense forward/inverse interscale transfer events and with decreasing scale $r_d$, it is natural to expect the skewness factor of ${\varPi }_{\bar {S}}$ to veer towards the skewness factors of $\mathcal {D}$ and $-\mathcal {\epsilon }$ which, as can be seen in figure 17, are highly negative with values between $-3.0$ and $-7.0$.

5.2. Correlations

We now consider the local/instantaneous relations between the various interscale transfer terms in terms of correlation coefficients plotted in figure 18(a). First, note the very strong correlation between ${\varPi }$ and ${\varPi }_{\bar {S}}$ and the moderate correlation between ${\varPi }$ and ${\varPi }_{\bar {I}}$. Even though ${\varPi }$ and ${\varPi }_{\bar {S}}$ are highly correlated, we cannot ignore ${\varPi }_{\bar {I}}$ and cannot write ${\varPi } \approx {\varPi }_{\bar {S}}$. As seen earlier in the paper, we cannot ignore ${\varPi }_{\bar {I}}$ because it is the part of the interscale transfer that balances the pressure term, but we have also seen that the fluctuation magnitude of ${\varPi }_{\bar {S}}$ is significantly higher than the fluctuation magnitude of ${\varPi }_{\bar {I}}$. However, even if smaller, the fluctuation magnitude of ${\varPi }_{\bar {I}}$ is not negligible. There is no correlation between ${\varPi }_{\bar {S}}$ and ${\varPi }_{\bar {I}}$ (see figure 18b), and so ${\varPi }$ correlates with both ${\varPi }_{\bar {S}}$ (strongly) and ${\varPi }_{\bar {I}}$ (moderately) for different independent reasons. Here ${\varPi }$ feels the influence of solenoidal vortex stretching/compression via ${\varPi }_{\bar {S}}$ and the influence of pressure fluctuations via ${\varPi }_{\bar {I}}$, the former influencing ${\varPi }$ more than the latter.

Figure 18. Correlation coefficients between various ${\varPi }$ decompositions $(\mathcal {Q}_1,\mathcal {Q}_2)$ listed on top of the figures at $\langle {{Re}}_\lambda \rangle _t=112$. (Corresponding plots for $\langle {{Re}}_\lambda \rangle _t=174$ are omitted because they are very similar.)

Figure 18(a) also shows significantly smaller correlations between ${\varPi }$ and ${\varPi }_{H}$ than between ${\varPi }$ and ${\varPi }_{\bar {S}}$. This must be due to a decorrelating effect of ${\varPi }_{H_{\bar {I}}}$ as ${\varPi }_{H}={\varPi }_{\bar {S}}+{\varPi }_{H_{\bar {I}}}$. The correlations between ${\varPi }$ and ${\varPi }_{I}$ are even smaller at the smaller scales but at integral size scales these correlations are equal to those between ${\varPi }$ and ${\varPi }_{H}$ (figure 18a).

Figure 18(b) reveals a strong anti-correlation between ${\varPi }_I$ and ${\varPi }_{H}$ at the small scales and a weak one at the large scales. As the scales decrease, the interscale transfers of fluctuating velocity differences caused by local/instantaneous non-homogeneities and the interscale transfers of fluctuating velocity differences caused by orientation differences get progressively more anti-correlated. This anti-correlation tendency results in ${\varPi }_{H}$ and ${\varPi }_{I}$ having larger fluctuation magnitudes than ${\varPi }$ at smaller scales, in particular scales smaller than $\langle \lambda \rangle _t$ (verified with our DNS data but not shown here for economy of space).

The other significant correlations revealed in figure 18(b) are those between ${\varPi }_H$ and ${\varPi }_{\bar {S}}$ and those between ${\varPi }_I$ and ${\varPi }_{\bar {S}}$, particularly as $r_d$ increases from around/below $\langle \lambda \rangle _t$ to the integral length scale. These correlations relate to the very strong correlations between ${\varPi }$ and ${\varPi }_{\bar {S}}$ but are weaker. One can imagine that ${\varPi }_{\bar {S}}$ correlates with ${\varPi }_H$ sometimes and with ${\varPi }_{I}$ some other times, but not too often with both given that ${\varPi }_{I}$ and ${\varPi }_H$ tend to be anti-correlated, and that this happens in a way subjected to a continuously strong correlation between ${\varPi } = {\varPi }_{H} + {\varPi }_{I}$ and ${\varPi }_{\bar {S}}$.

We finally consider in figure 19 the average contributions of the various ${\varPi }$-decomposition terms conditional on relatively rare intense ${\varPi }$ events. We calculate averages conditioned on $5\,\%$ most negative (forward transfer) ${\varPi }_{\bar {S}}$ events (values of ${\varPi }_{\bar {S}}$ for which the probability that ${\varPi }_{\bar {S}}$ is smaller than a negative value ${\varPi }_{\bar {S}_{0.05}}$ is 0.05) and on $5\,\%$ most positive ${\varPi }_{\bar {S}}$ (inverse transfer) events (values of ${\varPi }_{\bar {S}}$ for which the probability that ${\varPi }_{\bar {S}}$ is larger than a positive value ${\varPi }_{\bar {S}_{0.95}}$ is also 0.05). All these averages tend to $0$ as $r_d$ tends to $0$ below $\langle \lambda \rangle _t$. The largest such conditional averages are those of ${\varPi }\prime$ followed by those of ${\varPi }_{\bar {S}}^{\prime}$. This is the forward-skewed part of the interscale transfer (in terms of PDFs) and it is dominant at both forward and backward intense interscale transfer events. The weakest such conditional averages are those of ${\varPi }_{\bar {I}}$ for all $r_d$ and both forward and inverse extreme interscale transfer events. This is consistent with our observation in § 3.4 that the unconditional fluctuation magnitude of ${\varPi }_{\bar {I}}$ is smaller than the unconditional fluctuation magnitudes of ${\varPi }$ followed by those of ${\varPi }_{\bar {S}}$.

Figure 19. Average values of ${\varPi }$ decompositions conditioned on (a) intense backward events, (b) intense forward events at $\langle {{Re}}_\lambda \rangle _t=112$. The top of (a) lists the ${\varPi }$ decompositions. (Corresponding plots for $\langle {{Re}}_\lambda \rangle _t=174$ are omitted because they are very similar.)

The most interesting point to note in figure 19, however, is the difference between conditional averages of ${\varPi }_{H}^{\prime}$ and ${\varPi }_{I}$ when conditioned on intense forward or intense inverse interscale transfer events. Whilst the conditional averages of these two quantities are about the same at intense inverse events, they differ substantially at forward transfer events where the conditional average of $-{\varPi }_{H}^{\prime}$ is substantially higher than the conditional average of $-{\varPi }_{I}$ except close to the integral length scale.