Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-28T01:59:34.544Z Has data issue: false hasContentIssue false

On the simultaneous cascades of energy, helicity, and enstrophy in incompressible homogeneous turbulence

Published online by Cambridge University Press:  24 November 2022

P. Baj
Affiliation:
Department of Energy and Process Engineering, Norwegian University of Science and Technology, Trondheim 7491, Norway Faculty of Power and Aeronautical Engineering, Warsaw University of Technology, Warsaw 00-661, Poland
F. Alves Portela
Affiliation:
Univ. Lille, CNRS, ONERA, Arts et Metiers Institute of Technology, Centrale Lille, UMR 9014 – LMFL – Laboratoire de Mécanique des Fluides de Lille – Kampé de Fériet, F-59000 Lille, France Department of Aeronautical and Astronautical Engineering, University of Southampton, Southampton SO17 1BJ, UK
D.W. Carter*
Affiliation:
Department of Aeronautical and Astronautical Engineering, University of Southampton, Southampton SO17 1BJ, UK
*
Email address for correspondence: [email protected]

Abstract

We characterise the incompressible turbulence cascade in terms of the concurrent inter-scale and inter-space exchanges of the scale-by-scale energy, helicity and enstrophy. The governing equations for the scale-by-scale helicity and enstrophy are derived in a similar fashion to that of the second order structure function following Hill (J. Fluid Mech., vol. 468, 2002, pp. 317–326). We examine the instantaneous dynamics, applying these equations to forced periodic turbulence and a von Kármán flow focusing on scales in the dissipative range $r=2.5\eta$, the near-dissipative range $r=0.5\lambda$ and the onset inertial range $r=\lambda$ (where $\eta$ and $\lambda$ are the Kolmogorov and Taylor length scales, respectively). The signature of the random sweeping effect is observed in all three individual budgets and between the energy and enstrophy transfers. As in the energy cascade, the anti correlation of the pressure transport and non-linear transfer is identified also in the helicity cascade. Owing to its lack of positive definiteness, the helicity transfers are found to be decorrelated from the others. However a connection between the energy cascade and helicity is identified kinematically. This connection reveals the large-scale sweeping motions are a key element in the overall energy cascade and underpins previous observations of large-scale intermittency. Taken together, this work extends a classic framework to gain novel insight on turbulence dynamics that underlay the statistically steady state, and demonstrates how transfers are interconnected.

Type
JFM Papers
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2022. Published by Cambridge University Press.

1. Introduction

Turbulence is an inherently multi-scale phenomenon arising from the nonlinear nature of the motion of fluids. In general, a turbulent flow is populated by a hierarchy of eddies, the largest of which appear through interactions between the fluid and its boundaries (or some explicit forcing) and the smallest being conditioned by viscous forces (Falkovich Reference Falkovich2009), with energy flowing, on average, from the former to the latter. This statistical picture of turbulence is described by the so-called Richardson–Kolmogorov cascade (Richardson Reference Richardson1920; Kolmogorov Reference Kolmogorov1941a,Reference Kolmogorovb). However, this cascade is a result of dynamic interactions between eddies of various sizes (Pelz et al. Reference Pelz, Yakhot, Orszag, Shtilman and Levich1985). The distinction between the instantaneous interactions in the turbulence and the statistically stationary cascade that they give rise to are a well known problem in turbulence modelling. For example, in large eddy simulations (LES) the directionality of the cascade at any given instant cannot be predicted from statistically stationary models (Germano et al. Reference Germano, Piomelli, Moin and Cabot1991; Alexakis & Chibbaro Reference Alexakis and Chibbaro2020). Recent work by Goto & Vassilicos (Reference Goto and Vassilicos2016) and Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) highlights how viewing the cascade from a statistically stationary viewpoint overlooks important dynamics that determine how such a stationary state itself is established. The present work seeks to expand on that perspective by analysing the concurrent scale-by-scale budgets of energy, helicity and enstrophy by deriving the equations governing those cascades directly from the Navier–Stokes equations without any averaging or statistical assumptions. We then apply these equations to two distinct homogeneous turbulent flows: forced periodic turbulence and a von Kármán flow.

The kinetic energy $q^2=({\boldsymbol u}\boldsymbol {\cdot }{\boldsymbol u})/2$ is a central quantity in characterising turbulent flows and it has been, historically, at the core of theories of turbulence. Away from any boundaries, $q^2$ can only be destroyed by the turbulent dissipation $\varepsilon$ arising from small scale motions in the fluid (i.e. in the absence of viscosity, $q^2$ is a conserved quantity of the Navier–Stokes equations). Indeed, in the classical theory of Kolmogorov (Reference Kolmogorov1941b), where the details of the largest scales are dispensed with, the turbulence cascade is formulated in terms of a budget for $\delta q^2$, representing, to a large extent (Davidson & Pearson Reference Davidson and Pearson2005), the energy density at a given scale $r$. This energy cascade is governed solely by the nonlinear transfer of $\delta q^2$ from large to small scales which is in equilibrium with $\varepsilon$. Under the assumptions of homogeneity and isotropy, Kolmogorov closed the relationship between $q^2$ and $\varepsilon$ through the $4/5$-ths law leading to the inertial-scale distribution $\delta q^2 \sim r^{2/3}$ (the $5/3$-rds law for the energy spectrum), practically ubiquitous with turbulence (Kraichnan Reference Kraichnan1974).

In fact, $\varepsilon$ arises due to the finest structures in the turbulence that are characterised by strong vorticity ${\boldsymbol \omega }=\boldsymbol {\nabla } \times {\boldsymbol u}$ and therefore large enstrophy $\omega ^2={\boldsymbol \omega }\boldsymbol {\cdot }{\boldsymbol \omega }$. While dynamically distinct (see discussion of Carbone & Bragg Reference Carbone and Bragg2020), it can be shown that the dissipation and enstrophy are closely related (Raynal Reference Raynal1996) upon integration over all space (Tsinober Reference Tsinober2001; Tennekes & Lumley Reference Tennekes and Lumley2018). Enstrophy exhibits its own cascade across scales but, due to the presence of vortex stretching, this cascade is not inviscidly conservative. Indeed, it is known since the work of Fjørtoft (Reference Fjørtoft1953) that in two dimensions (i.e. in the absence of vortex stretching), enstrophy cascades from large to small scales, effectively reversing the direction of the energy cascade (Eyink Reference Eyink1996; Kraichnan Reference Kraichnan1967). Although purely two-dimensional flows are seldom realised in nature, such an idealised cascade is an important element in our understanding of a variety of turbulent flows (see reviews by Boffetta & Ecke Reference Boffetta and Ecke2012; Falkovich et al. Reference Falkovich, Boffetta, Shats and Lanotte2017) as would be the case of geostrophic flows (Lindborg Reference Lindborg2007; Lindborg & Nordmark Reference Lindborg and Nordmark2022), strongly rotating turbulence (Deusebio et al. Reference Deusebio, Boffetta, Lindborg and Musacchio2014; van Kan & Alexakis Reference van Kan and Alexakis2022) but also in flows that exhibit coherent structures (Dascaliuc & Grujić Reference Dascaliuc and Grujić2013). Recently, Bos (Reference Bos2021) analysed three-dimensional turbulence, where vortex stretching is suppressed altogether; they found that enstrophy is preserved but energy is not and showed how this modified system displayed a dual direct cascade of both energy (no longer inviscidly conserved) and enstrophy.

Vorticity also features in the helicity $h = {\boldsymbol u}\boldsymbol {\cdot }{\boldsymbol \omega }$, which is another inviscid invariant of the Navier–Stokes equations and acts as a measure of breakage in mirror-symmetry, or parity-invariance, within the fluid (Moreau Reference Moreau1960; Moffatt Reference Moffatt1969). Topologically, it describes the degree of knottedness of vortex tubes (Moffatt & Tsinober Reference Moffatt and Tsinober1992), a concept which is quite useful in the study of superfluid turbulence, where very large Reynolds numbers can be achieved (see e.g. Kleckner, Kauffman & Irvine Reference Kleckner, Kauffman and Irvine2016; Kivotides & Leonard Reference Kivotides and Leonard2021). As noted by Brissaud et al. (Reference Brissaud, Frisch, Léorat, Lesieur and Mazure1973), because $h$ is not positive definite (as are $q^2$ and $\omega ^2$), the possibility of a dual direct cascade exists in tandem with that of one analogous to the two-dimensional scenario described above (see also Kraichnan Reference Kraichnan1973). It is established that the most physically sound scenario is that of a dual (direct) cascade (see e.g. Chen, Chen & Eyink Reference Chen, Chen and Eyink2003). The bulk of the work on helicity cascades has been carried out in Fourier space (see Scott & Wang Reference Scott and Wang2005; Alexakis Reference Alexakis2017; Alexakis & Biferale Reference Alexakis and Biferale2018; Pouquet et al. Reference Pouquet, Rosenberg, Stawarz and Marino2019), where the nonlinearity of the Navier–Stokes equations is more easily described by triad interactions of wavenumbers (as motivated by Waleffe Reference Waleffe1992) at the sacrifice of requiring homogeneity. For flows in which some degree of anisotropy or inhomogeneity is present, non-zero helicity can give rise to the spontaneous formation of large-scale coherence in the flow (see Yokoi & Yoshizawa Reference Yokoi and Yoshizawa1993, and references therein). Their effect on the cascade is likely non-trivial.

Recognising that the Richardson–Kolmogorov cascade is insufficient for a general description of turbulence, Hill (Reference Hill2002) derived a budget for $\delta q^2$ that reduces to that of Kolmogorov (under the assumptions of homogeneity, isotropy and stationarity) but is instead obtained directly from the Navier–Stokes equations without requiring any information regarding the structure of the flow or the Reynolds number. This equation, often referred to as the generalised Kolmogorov equation (Mollicone et al. Reference Mollicone, Battista, Gualtieri and Casciola2018; Gatti et al. Reference Gatti, Chiarini, Cimarelli and Quadrio2020) or the Kármán–Howarth–Monin-Hill equation (Alves Portela, Papadakis & Vassilicos Reference Alves Portela, Papadakis and Vassilicos2017; Yasuda & Vassilicos Reference Yasuda and Vassilicos2018), characterises the energy cascade in terms of physical and scale space exchanges of $\delta q^2$ (as well as any sources or sinks) such that the nonlinear transfer of $\delta q^2$ is but one of the mechanisms in balance with $\varepsilon$, which always acts as a sink of $\delta q^2$. In recent years, this equation has been extensively used to characterise the effects of inhomogeneity and anisotropy on the energy cascade (Gomes-Fernandes, Ganapathisubramani & Vassilicos Reference Gomes-Fernandes, Ganapathisubramani and Vassilicos2015; Cimarelli et al. Reference Cimarelli, De Angelis, Jimenez and Casciola2016; Knutsen et al. Reference Knutsen, Baj, Lawson, Bodenschatz, Dawson and Worth2020; Zimmerman et al. Reference Zimmerman, Antonia, Djenidi, Philip and Klewicki2022) and has been extended to variable density and compressible flows (Lai, Charonko & Prestridge Reference Lai, Charonko and Prestridge2018; Arun et al. Reference Arun, Sameen, Srinivasan and Girimaji2021).

It is clear that a full characterisation of the turbulence cascade must involve not only the turbulent kinetic energy but also quantities such as enstrophy and helicity, as they bear relation to dissipation (a core element of the energy cascade) and capture the presence of coherence, intermittency and other such phenomena that are known to break with the classical picture of turbulence. Following this observation, the objectives of the present study are threefold:

  1. (i) To provide a generalised framework for analysing the transfers of energy, enstrophy and helicity through two-point equations for these quantities, following Hill (Reference Hill2002).

  2. (ii) To combine this framework with the correlation-based analysis of the energy transfers (Yasuda & Vassilicos Reference Yasuda and Vassilicos2018), examining how the different mechanisms involved in the scale-by-space energy, helicity and enstrophy budgets are related.

  3. (iii) To leverage the two-point framework towards a deeper understanding of the connection between the instantaneous energy cascade and helicity.

We begin with the formal derivation of the scale-space equations in § 2. A brief description and characterisation of the data sets on which the scale-space equations will be applied is given in § 3. The correlations between the individual terms of each of the three transfer budgets is shown in § 4. In § 5 we discuss the role played by helicity in the energy cascade from a two-point scale-by-scale perspective. We conclude with a summary of our results and suggestions for future work in § 6.

2. Formulation

The single-point quantities of interest that will be cast into the scale-space framework are the energy $q^2$, helicity $h$ and enstrophy $\omega ^2$. The formulation begins with the familiar single-point incompressible Navier–Stokes equations,

(2.1)\begin{equation} \frac{\partial}{\partial t} u_i + u_k \frac{\partial u_i}{\partial x_k} ={-}\frac{1}{\rho}\frac{\partial p}{\partial x_i} + \nu \frac{\partial^2 u_i}{\partial x_k \partial x_k}, \end{equation}

together with the incompressibility condition $\boldsymbol {\nabla } \boldsymbol {\cdot } \boldsymbol {u} = 0$, where $u_i(\boldsymbol {x},t)$ and $p(\boldsymbol {x},t)$ are the instantaneous velocity and pressure fields, respectively, with $\nu$ the kinematic viscosity and $\rho$ the density. In addition to (2.1), the vorticity form of the Navier–Stokes equations is invoked taking the curl:

(2.2)\begin{equation} \frac{\partial}{\partial t} \omega_i + u_k\frac{\partial \omega_i}{\partial x_k} = \omega_k \frac{\partial u_i}{\partial x_k} + \nu \frac{\partial^2 \omega_i}{\partial x_k \partial x_k}, \end{equation}

where $\boldsymbol {\omega } = \boldsymbol {\nabla } \times \boldsymbol {u}$ is the vorticity with the kinematic property $\boldsymbol {\nabla } \boldsymbol {\cdot } \boldsymbol {\omega } = 0$.

2.1. Generalised scale-to-scale transfers

In this section we outline the general method for obtaining the scale-space energy, helicity and enstrophy budget equations. These originate from studies focused on the energy and can be found in the classic work of Von Kármán & Howarth (Reference Von Kármán and Howarth1938) and later made general in order to account for anisotropy (Monin & Yaglom Reference Monin and Yaglom1975) and inhomogeneity (Hill Reference Hill2002). As such the evolution equations are commonly referred to as the generalised Kármán–Howarth–Monin (KHM) equations or the Kármán–Howarth–Monin–Hill (KHMH) equations. For the present work, these equations will be used to describe the scale-space dynamics of three specific structure functions: the energy structure function $\delta q^2= \delta u_i \delta u_i$, the helicity structure function $\delta h = \delta u_i \delta \omega _i$ and the enstrophy structure function $\delta \omega ^2 = \delta \omega _i \delta \omega _i$. Here, $\delta u_i$ and $\delta \omega _i$ represent, respectively, velocity and vorticity increments, taken as the difference of that quantity at two points $x_i$ and $x_i'$, i.e. $\delta u_i = u_i-u_i'$ and $\delta \omega _i = \omega _i-\omega _i'$. Henceforth, primed variables denote belonging to a set of spatial points $x_i'$ independent of $x_i$. A detailed derivation and interpretation of the KHMH equations in the scale-space coordinates system is given in Hill (Reference Hill2002) and Marati, Casciola & Piva (Reference Marati, Casciola and Piva2004). In the present study we abstain from using the Reynolds decomposition, instead following a similar approach to Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) where the instantaneous dynamics of each term at various scales can be analysed.

The derivation of the transfer equations begins with the equations for velocity and vorticity, (2.1) and (2.2). There are three common steps in deriving the budgets for $\delta q^2$, $\delta h$ and $\delta \omega ^2$:

  1. (i) Take the difference of (2.1) and (2.2) evaluated on two independent coordinates $x_k$ and $x_k'$. These constitute the two-point equations for the velocity and vorticity differences ($\delta u_i = u_i-u_i'$ and $\delta \omega _i = \omega _i-\omega _i'$, respectively)

    1. (a) For $\delta q^2$: Multiply the equation for $\delta u_i$ by $2\delta u_i$.

    2. (b) For $\delta h$: Add the product between $\delta \omega _i$ and the equation for $\delta u_i$ to the product between $\delta u_i$ and the equation for $\delta \omega _i$.

    3. (c) For $\delta \omega ^2$: Multiply the equation for $\delta \omega _i$ by 2$\delta \omega _i$.

  2. (ii) Express all quantities as two-point differences ($\delta$) or sums ($\varSigma$). The strain rate tensor, that appears in the equation for $\delta \omega ^2$, is expressed as $S_{ik} = ({\partial u_i}/{\partial x_k} + {\partial u_k}/{\partial x_i} )/2$ with $\delta S_{ik}=S_{ik}-S_{ik}'$ and $\varSigma S_{ik}=S_{ik}+S_{ik}'$.

  3. (iii) Change the coordinate system from $(x_k,x_k')$ to $(X_k,r_k)$, where $X_k \equiv (x_k+x_k')/2$ represents physical space (in the sense that it represents a centroid position) and $r_k \equiv x_k'-x_k$ represents the space of scales (as it relates to a distance vector).

The resulting equations describe the various budgets of the scale-space energy, helicity and enstrophy structure functions. For the energy structure function $\delta q^2$ one obtains:

(2.3)\begin{align} \frac{\partial}{\partial t}\delta q^2 + \frac{\partial \delta u_k\delta q^2 }{\partial r_k} &=-\frac{\partial}{\partial X_k}\left(\frac{\varSigma u_k\delta q^2}{2}\right)- \frac{2}{\rho}\frac{\partial \delta u_k \delta p}{\partial X_k}\nonumber\\ &\quad +\nu\left[2\frac{\partial^2}{\partial r_k^2}+\frac{1}{2}\frac{\partial^2}{\partial X_k^2}\right]\delta q^2 -2\nu\left[\left(\frac{\partial u_i}{\partial x_k}\right)^2+\left(\frac{\partial u_i'}{\partial x_k'}\right)^2 \right] . \end{align}

In shorthand notation:

(2.4)\begin{equation} \mathcal{A}_t^{q^2} + \varPi^{q^2} ={-}\mathcal{T}^{q^2}-\mathcal{T}_p^{q^2}+\mathcal{D}_{\nu}^{q^2}-\mathcal{E}^{q^2}, \end{equation}

where the superscript denotes the energy $q^2$, and the terms of (2.4) reflect the terms presented in (2.3) sequentially. For the helicity structure function $\delta h$ one obtains:

(2.5)\begin{align} &\frac{\partial}{\partial t}\delta h + \frac{\partial}{\partial r_k}\left[\delta u_k\delta h -\frac{1}{2}\delta \omega_k \delta q^2 \right] ={-}\frac{\partial}{\partial X_k}\left[\frac{\varSigma u_k\delta h}{2}-\frac{\varSigma \omega_k\delta q^2}{4}\right]\nonumber\\ &\quad - \frac{1}{\rho}\frac{\partial \delta \omega_k \delta p}{\partial X_k} +\nu\left[2\frac{\partial^2}{\partial r_k^2}+\frac{1}{2}\frac{\partial^2}{\partial X_k^2}\right]\delta h -2\nu\left[\left(\frac{\partial \omega_i}{\partial x_k}\frac{\partial u_i}{\partial x_k}\right)+\left(\frac{\partial \omega_i'}{\partial x_k'}\frac{\partial u_i'}{\partial x_k'}\right) \right] . \end{align}

In shorthand notation:

(2.6)\begin{equation} \mathcal{A}_t^{h} + \varPi^{h} ={-}\mathcal{T}^{h}-\mathcal{T}_p^{h}+\mathcal{D}_{\nu}^{h}-\mathcal{E}^{h}, \end{equation}

where the superscript denotes the helicity $h$ for the respective terms of (2.5). Finally, for the enstrophy $\delta \omega ^2$ one obtains:

(2.7)\begin{align} &\frac{\partial}{\partial t}\delta \omega^2 + \frac{\partial \delta u_k\delta \omega^2}{\partial r_k} =-\frac{\partial}{\partial X_k}\left(\frac{\varSigma u_k\delta \omega^2}{2}\right)+ (\delta \omega_k \delta \omega_i \varSigma S_{ik} + \varSigma \omega_k \delta \omega_i \delta S_{ik})\nonumber\\ &\quad + \nu\left[2\frac{\partial^2}{\partial r_k^2}+\frac{1}{2}\frac{\partial^2}{\partial X_k^2}\right]\delta \omega^2 -2\nu\left[\left(\frac{\partial \omega_i}{\partial x_k}\right)^2+\left(\frac{\partial \omega_i'}{\partial x_k'}\right)^2 \right]. \end{align}

In shorthand notation:

(2.8)\begin{equation} \mathcal{A}_t^{\omega^2} + \varPi^{\omega^2} ={-}\mathcal{T}^{\omega^2}+\mathcal{G}_{S}^{\omega^2}+\mathcal{D}_{\nu}^{\omega^2}-\mathcal{E}^{\omega^2}, \end{equation}

where the superscript denotes the enstrophy $\omega ^2$.

Notice that these three budgets are similarly structured, where:

  1. (i) $\mathcal {A}_t$ is an unsteady term and represents the temporal increase or decrease of scale energy/helicity/enstrophy at each instant.

  2. (ii) $\varPi$ represents the nonlinear exchanges in scale space of energy/helicity/enstrophy.

  3. (iii) $\mathcal {T}$ is the nonlinear turbulent transport in physical space.

  4. (iv) $\mathcal {T}_p$ results from the interaction of the pressure and velocity/vorticity fields to produce a pressure transport that acts to transport energy/helicity at a particular scale.

  5. (v) $\mathcal {G}_S$ is a generation term in scale space resulting from the coupling between the rate-of-strain and the enstrophy.

  6. (vi) $\mathcal {D}_{\nu }$ represents the viscous diffusion of energy/helicity/enstrophy both in scale and physical space.

  7. (vii) $\mathcal {E}$ represents the two-point average dissipation rate. This can be seen, for example for the scale-to-scale energy transfers, dividing equation (2.3) by 4 on both sides such that $\mathcal {E}=(\varepsilon + \varepsilon ')/2$, where $\varepsilon = \nu (\partial u_i/\partial x_k)^2$.

We refer to Marati et al. (Reference Marati, Casciola and Piva2004), Danaila et al. (Reference Danaila, Krawczynski, Thiesset and Renou2012), Valente & Vassilicos (Reference Valente and Vassilicos2015), Alves Portela et al. (Reference Alves Portela, Papadakis and Vassilicos2017), Mollicone et al. (Reference Mollicone, Battista, Gualtieri and Casciola2018) and Gatti et al. (Reference Gatti, Chiarini, Cimarelli and Quadrio2020), who give a more detailed interpretation of the various terms in the scale-space framework.

There are several useful observations that follow immediately from deriving the evolution equations for the scale energy, helicity and enstrophy. These equations resemble their one-point counterparts except that they account for nonlinear exchanges across scales through $\mathcal {\varPi }$. Closer inspection reveals several differences amongst the three budgets. It is clear that transfers of helicity (2.5) not only arise from the interactions of the scale to scale helicity $\delta h$ with the velocity increment $\delta u_k$, but also from the interaction of the scale-to-scale energy $\delta q^2$ with the vorticity increment $\delta \omega _k$. This is consistent with the two channels of inter-scale transfer identified in the study of Yan et al. (Reference Yan, Li, Yu, Wang and Chen2020), who highlight distinct helicity transfers through the combined action of vortex twisting and vortex stretching arising naturally from the vorticity equation. The pressure term only has an explicit role in the scale energy and helicity budgets. For the scale enstrophy budget, the pressure transport is effectively ‘curled out’; however a non-conservative generation $\mathcal {G}_S$ term emerges whose role is dynamically distinct compared to the pressure transport.

3. Experimental and numerical datasets

To exemplify the scale-space framework, two homogeneous turbulence data sets, one experimental and one numerical, are selected. We opt to limit the present scope to homogeneous turbulence for its relative simplicity, but remark that the instantaneous framework is equally applicable to inhomogeneous flows.

3.1. Von-Kármán mixing tank

The experimental data set considered in this work was gathered in the Göttingen Turbulence Facility #3. The experimental rig consisted of a steel cylinder of diameter and height equal to $48\ {\rm cm}$ and $58\ {\rm cm}$ respectively. It featured $8$ axial baffles, attached to its inside wall, and two counter-rotating impellers of diameter $25\ {\rm cm}$ driving the flow inside the cylinder. The rotation frequency was set to $0.2\ {\rm Hz}$. Water was used as the working fluid (kinematic viscosity of $\nu = 0.98\ {\rm mm}^2\ {\rm s}^{-1}$). The flow was seeded with PMMA microspheres with mean diameter $6.0\ {\rm m}$ and specific gravity $\rho _{seeding} /\rho = 1.22$ (the resultant particle Stokes number was estimated as $6\times 10^{-5} \ll 1$). A high-speed Nd:YAG pulse laser was used as the source of illumination. Two high-speed Phantom v640 cameras were used to record the flow. A more in-depth description can be found in Knutsen et al. (Reference Knutsen, Baj, Lawson, Bodenschatz, Dawson and Worth2020).

The flow inside a cubic domain of $8.5\times 8.5\times 8.5\ {\rm mm}^3$ located at the centre of the tank (where the mean flow vanishes) was measured using scanning particle image velocimetry (scanning PIV). This measurement technique is extensively discussed in Lawson & Dawson (Reference Lawson and Dawson2014), and here we only include a brief summary of its principles. In the first step, two cameras acquire multiple stereo PIV images, with a very small time separation, as the laser sheet traverses the measurement volume. A volumetric snapshot of the scattered light intensity is subsequently reconstructed from the recorded series of images. Ultimately, cross-correlation of the reconstructed snapshots is performed, yielding a single, volumetric velocity field.

Short, time-resolved sequences of velocity snapshots were captured during the experiment (a correction proposed by Wang et al. Reference Wang, Gao, Wei, Li and Wang2017, was used to reduce the residual divergence of the data). Each sequence consisted of six snapshots with the time separation below one-tenth of the Kolmogorov time scale. A vast collection of $2\times 10^5$ such sequences was gathered throughout the experiment. They were considered independent samples of the velocity field as the time separation between sequences was of the order of the impeller revolution period. A Lagrangian filtering, similar to that proposed by Novara & Scarano (Reference Novara and Scarano2013), was applied to the data in the post-processing phase. Each sequence was used to advect artificial tracers, whose initial positions coincided with the measurement grid points, forward and backward in time. The resultant traces were utilised to evaluate the filtered velocity and Lagrangian acceleration vectors. The approach described by Lawson & Dawson (Reference Lawson and Dawson2015) was subsequently employed to reconstruct pressure fields. The discretised momentum equation was rearranged to form an over-determined set of linear equations used to solve for pressure in the least-square sense.

The experimental data was stored as cubes containing fields for all components of the velocity, material derivatives and pressure field (each cube corresponded to one time-resolved measurement sequence). A fourth-order central difference scheme was used to evaluate spatial derivatives.

3.2. Direct numerical simulation

A direct numerical simulation (DNS) data set of forced homogeneous and isotropic turbulence (HIT) provided by the Johns Hopkins University (JHU) turbulence database was used to provide numerical comparison (Li et al. Reference Li, Perlman, Wan, Yang, Meneveau, Burns, Chen, Szalay and Eyink2008). The DNS is triply periodic on a cubic domain of size length $2{\rm \pi}$ with Reynolds number $Re_{\lambda }=433$ on a regular grid of $1024^3$ points. In time intervals of one second, 256 spatially independent sub-cubes of side length $5\lambda$ ($0.4L$) were extracted from the database at full resolution for a total of 2560 sub-cubes over 10 time units (or $5T_L$, where $T_L=u_{rms}/L$ is the large eddy turnover time scale), for which third-order statistics (such as the structure functions reported in § 3.4.2) were found to satisfactorily converge.

The DNS data were stored similarly to the experimental data, where each cube contained fields for all components of the velocity, the material derivatives (accounting for the numerical forcing) and the pressure field. Fourth-order differentiation in physical space was applied directly to the fields as necessary to calculate gradient quantities used in the analysis. As the JHU database enforces the divergence-free condition via spectral methods, a small but non-negligible divergence residual was found for the gradients in physical space. This error was found to propagate to higher-order derivatives, impacting the energy, helicity and enstrophy budget residuals at approximately 2 %, 5 %, and 15 %, respectively (Appendix A). Though this error is non-negligible, it is not believed to overshadow the correlation-based conclusions drawn in this study that are primarily rooted in phase information. For the Von-Kármán (VK) data, the error propagation of gradients in physical space is more significant as a result of the limited sub-pixel accuracy. This is discussed in more detail in Appendix A but is particularly problematic for enstrophy dissipation (fourth order gradients). As a result, conclusions based on such high-order gradients in the VK data are avoided.

3.3. Characterisation of the data sets

In this section we characterise the data sets to elucidate the statistical nature of the scale-space energy, helicity and enstrophy quantities before investigating their transfers in § 4. Basic statistical quantities for both data sets are reported in table 1.

Table 1. Parameters from the Von-Kármán mixing tank scanning PIV experiment (Knutsen et al. Reference Knutsen, Baj, Lawson, Bodenschatz, Dawson and Worth2020) and the Johns Hopkins University DNS of homogeneous and isotropic turbulence (Li et al. Reference Li, Perlman, Wan, Yang, Meneveau, Burns, Chen, Szalay and Eyink2008). Note the symbols $\Delta x$ and $\Delta t$ denote the grid spacing and temporal spacing of the data sets.

The probability density functions (p.d.f.s) of the dissipation rates of the energy, helicity and enstrophy are shown in figure 1. These distributions are of central importance towards comprehending the scale-by-scale cascades, as both energy and helicity can only be destroyed by viscous dissipation. The wide tails of these p.d.f.s extend to extreme values resulting in large flatness (Van Atta & Antonia Reference Van Atta and Antonia1980) and are indicative of small-scale intermittency, characteristic of turbulent flows (Sreenivasan & Antonia Reference Sreenivasan and Antonia1997; Laval, Dubrulle & Nazarenko Reference Laval, Dubrulle and Nazarenko2001). The dissipation p.d.f.s agree well between data sets for the energy and helicity but deviate at extreme values for the dissipation of enstrophy. This is due to the unavoidable PIV error propagation (see Appendix A).

Figure 1. Premulitplied p.d.f.s of the dissipation of energy (a), helicity (b) and enstrophy (c). The inset of (b) highlights the range of helicity dissipation between ${\pm }5$ standard deviations.

3.4. Ensemble-averaged structure functions

In this section, in addition to the classical energy structure function, the structure functions associated with the helicity and enstrophy are probed; these are central to core theories of turbulence (Zhou Reference Zhou2021). This analysis provides further context on the extension of the scale-space framework beyond that of the scale-space energy budget.

3.4.1. Energy, helicity and enstrophy structure functions

The energy, helicity and enstrophy structure functions are presented in figure 2. Here, and for all scale-space quantities, the overbar denotes both ensemble averaging and averages over all scale-space orientations at a particular scale $r=\|{\boldsymbol r}\|$. We note that here (and elsewhere in this study) the error bars correspond to statistical errors obtained using the convergence of the moments of the p.d.f.s (using a 95 % confidence interval). Though bias errors arising from experimental artefacts and discretisation are also of concern, they do not dominate the uncertainty. This is supported by the small residual in the energy, helicity and enstrophy budgets (Appendix A).

Figure 2. Normalised mean energy structure function $\overline {\delta q^2}$ (a), helicity structure function $\bar {\delta h}$ (b) and enstrophy structure function $\overline {\delta \omega ^2}$ (c).

The energy structure functions in figure 2(a) are compensated by the classical inertial range scaling $(r\overline {\epsilon ^{q^2}})^{2/3}$ (Kolmogorov Reference Kolmogorov1941a). The universal constant for the analogous longitudinal second-order structure function is $C_2\approx 2.0\pm 0.2$ (Sreenivasan Reference Sreenivasan1995; Pope Reference Pope2001); however, the energy structure functions feature a sum over two transverse components (that scale as $\tfrac {4}{3}C_2$) and one longitudinal component. Therefore, in the inertial range, the energy structure function is expected to plateau at a value of approximately $7.3\pm 0.7$. For the VK data, this value is achieved only at the largest available separation. This is likely due to the limited range of separations but may also be attributed to the Reynolds number ($Re_{\lambda }=199$) for which the inertial range is expected to be limited in breadth (Ishihara, Gotoh & Kaneda Reference Ishihara, Gotoh and Kaneda2009). For the JHU data the plateau is more evident; however, it is also considerably limited due to the sizes of the extracted sub-cubes of side length $5\lambda$.

The helicity structure functions are shown in figure 2(b). As the helicity is not a sign-definite quantity, the structure functions were found to converge very slowly due to their characteristically high variation (Kurien, Taylor & Matsumoto Reference Kurien, Taylor and Matsumoto2004). The helicity structure functions are normalised using the inertial range scaling in physical space as $\bar {\delta h} /\overline {\epsilon ^{h}}(r^2/\overline {\epsilon ^{q^2}})^{1/3}$, outlined by Brissaud et al. (Reference Brissaud, Frisch, Léorat, Lesieur and Mazure1973). A developing plateau is observed for the JHU data for $r>\lambda$, in agreement with the energy structure function. Also consistent with its energy structure function and limited range of separations, the VK data does not show a clear plateau for $\bar {\delta h}$. The overlap within uncertainty between the two data sets indicates satisfactory collapse, though curiously this occurs mostly for $r$ in the near-dissipative range. The lack of large separations $r$ in the VK data prevents from evaluating the collapse of the data sets for $r>\lambda$ (in the inertial range).

Despite the lack of explicit helical forcing in the DNS, the scale-space helicity is non-zero for both data sets. This is consistent with the helicity spectrum computed over a wide range of simulations (Chen et al. Reference Chen, Chen and Eyink2003; Mininni, Alexakis & Pouquet Reference Mininni, Alexakis and Pouquet2006) with and without helical forcing (Alexakis Reference Alexakis2017). This can be seen directly through an expansion of $\bar {\delta h}$, i.e. for homogeneous turbulence, $\overline {(u_i-u_i')(\omega _i-\omega _i')}=2\bar {h}-2\overline {u_i\omega _i'}$. This reveals that even when small-scale mirror symmetry holds ($\bar {h}=0$), the non-zero scale space helicity (and spectrum of helicity) arises from the coherence in the velocity–vorticity correlation (Levich & Shtilman Reference Levich and Shtilman1988). This was shown explicitly for isotropic turbulence by Deusebio et al. (Reference Deusebio, Boffetta, Lindborg and Musacchio2014). The present results confirm this non-zero correlation as well as the tendency of the turbulence cascade to restore small-scale mirror symmetry with decreasing scale (Kraichnan Reference Kraichnan1973; Chen et al. Reference Chen, Chen and Eyink2003).

The enstrophy structure functions are presented in figure 2(c), normalised using the energy dissipation rate and the Kolmogorov length scale $\eta =(\nu ^3/\overline {\epsilon ^{q^2}})^{1/4}$ as $(\overline {\epsilon ^{q^2}}/\eta ^2)^{2/3}$. A normalisation using $\overline {\epsilon ^{\omega ^2}}$ was tested but found to give unsatisfactory agreement (likely due to error propagation, see Appendix A). The close relationship between the enstrophy and the dissipative small scales (Jiménez et al. Reference Jiménez, Wray, Saffman and Rogallo1993) implies a constant scaling in the inertial range. Consistent with the energy and helicity structure functions, the enstrophy structure functions reach a maximum only for the largest separations and more conclusively for the JHU data. This reaffirms enstrophy as a quantity that is confined to small scales (Jiménez et al. Reference Jiménez, Wray, Saffman and Rogallo1993; Davidson, Morishita & Kaneda Reference Davidson, Morishita and Kaneda2008; Ishihara, Kaneda & Hunt Reference Ishihara, Kaneda and Hunt2013; Elsinga et al. Reference Elsinga, Ishihara, Goudar, Da Silva and Hunt2017).

3.4.2. Nonlinear flux energy, helicity and enstrophy structure functions

In this section the nonlinear flux (NLF) structure functions, of which the divergence is taken in the nonlinear inter-scale transfers $\varPi$, are explored. A small departure to justify this terminology will be taken here. These structure functions are commonly referred to as ‘third-order structure functions’, but as the origin of ‘third order’ is in the statistical moment of the increment of a single quantity (i.e. $\overline {\delta u^3}$), a different terminology is adopted for improved generality. As noted by Hill (Reference Hill2002), making use of Gauss’ theorem when integrating the nonlinear energy transfer term of (2.3) inside a ball $\mathcal {V}_R$ (defined in the space of scales as $\{\boldsymbol {r}\in \mathbb {R}^{3} : \lvert \boldsymbol {r} \rvert \le R \}$) gives $\unicode{x222D} \varPi ^{q^2} \,\mathrm {d}\mathcal {V}_R = \unicode{x222F} _{\partial \mathcal {V}_R}\overline {\delta u_k\delta q^2} n_k \,\mathrm {d}\mathcal {S}$ with $n_k$ the outward normal vector and $d\mathcal {S}$ the surface of the spherical shell (the boundary of $\mathcal {V}_R$). The orientation average over the spherical shell is identically zero, leaving only a flux in the radial direction. This motivates referring to these quantities as NLF structure functions in the present context. In addition to the central role of NLF structure functions in classical turbulence theory, i.e. the 4/5ths (K41) and 2/15ths laws (Chkhetiani Reference Chkhetiani1996; L'vov, Podivilov & Procaccia Reference L'vov, Podivilov and Procaccia1997), their physical significance is well documented in the context of a spherical scale-space coordinate system (Gomes-Fernandes et al. Reference Gomes-Fernandes, Ganapathisubramani and Vassilicos2015; Valente & Vassilicos Reference Valente and Vassilicos2015; Alves Portela et al. Reference Alves Portela, Papadakis and Vassilicos2017). The NLF structure functions for the experimental and numerical data sets are presented in figure 3.

Figure 3. Normalised orientation-averaged NLF structure functions corresponding to the nonlinear scale-to-scale flux of energy $\overline {\delta q^2}$ (a), helicity $\bar {\delta h}$ via advection (b) and vortex stretching (c), and enstrophy $\overline {\delta \omega ^2}$ (d).

The orientation-averaged NLF energy structure function is shown in figure 3(a), normalised using the classical 4/3rds law (Hill Reference Hill2002). The subscript $r$ is used to denote the radial flux as described above. For $r>\lambda$ the plateau seen in the JHU data suggests classical behavior of $\overline {\delta u_r \delta q^2}$. Consistent with figure 2, the VK data does not show a clear plateau developing for $\overline {\delta u_r \delta q^2}/\tfrac {4}{3}\overline {\varepsilon ^{q^2}}r$ but does reach value close to one at $r\approx \lambda$ (the largest available separation). This may imply the nonlinear inter-scale transfer is dominant at $r=\lambda$ (Yasuda & Vassilicos Reference Yasuda and Vassilicos2018), but it must be noted that such an apparent balance is possible when some inhomogeneous effects are at play (Alves Portela et al. Reference Alves Portela, Papadakis and Vassilicos2017).

The orientation-averaged NLF helicity structure functions are shown in figures 3(b,c), normalised using the product of $r$ and dissipation rate of helicity $\overline {\epsilon ^h}$. This normalisation is chosen to compare the two distinct mechanisms associated with $\mathcal {\varPi }^h$ in (2.5). The first, $\overline {\delta u_k \delta h}$, originates from the advective term and the second, $\overline {\delta \omega _k \delta q^2}$, from the vortex stretching term of (2.2) (Yan et al. Reference Yan, Li, Yu, Wang and Chen2020). It is again seen that, due to its non-positive definiteness and inherently high variation, the nonlinear flux of helicity is very slow to converge, leading to large uncertainty bars. Despite this, the fluxes are found to be non-zero, and the JHU data appears to exhibit a plateau (albeit within large uncertainty bars).

Finally, the normalised NLF enstrophy structure functions can be seen in figure 3(d). Similarly to the enstrophy structure function, a collapse of the data sets is observed in the near-dissipation range between $r/\eta =10$ and $r/\lambda =1$. The flux is seen to decrease as $r^{-3}$ ($r^{-2}$ in the figure due to compensating $-\overline {\delta u_r\delta \omega ^2}$ by $r$) into the inertial range. This power law decrease is consistent with the results of Davidson et al. (Reference Davidson, Morishita and Kaneda2008) (in their case, a scale-by-scale enstrophy flux is defined that is not motivated by the present spherical scale-space coordinate system. This necessitates compensating $-\overline {\delta u_r\delta \omega ^2}$ by $r$ for a one-to-one comparison).

4. Correlations of instantaneous cascades

4.1. Correlations within budgets

Normalised correlations ($-$1: perfectly anti-correlated, 0: uncorrelated, and 1: perfectly correlated) were tabulated as outlined in Appendix A. We shall focus on the separations $r=2.5\eta$, $0.5\lambda$ and $\lambda$ that correspond to the dissipative and near-dissipative ranges and the onset of inertial range, respectively. The latter is of particular interest, as it is the scale for which Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) identified the largest variations in energy transfers (standard deviation) in periodic box turbulence with $Re_{\lambda }=178$, indicative of strong dynamics at that scale.

The correlations of the terms in the scale energy, helicity and enstrophy budgets are presented in figures 4, 5 and 7, respectively. Starting with the energy in figure 4, the correlations of terms are in good agreement across the three considered scales between both the JHU and VK data. At the smallest $r$, the high correlation between the viscous diffusion and dissipation seen in figure 4(a,d) reflects the average balance between $\mathcal {D}_{\nu }^{q^2}$ and $\mathcal {E}^{q^2}$, as noted by (Valente & Vassilicos Reference Valente and Vassilicos2015).

Figure 4. Correlation coefficients between terms of (2.4) for $\delta q^2$ at various separations; results based on VK (ac) and JHU (df) datasets.

Figure 5. Correlation coefficients between terms of (2.6) for $\delta h$ at various separations; results based on VK (ac) and JHU (df) datasets.

The most robust correlation of figure 4 lies between the unsteady transport $\mathcal {A}_t^{q^2}$ and the turbulent transport $\mathcal {T}^{q^2}$. This is taken to be the signature of the random sweeping effect identified by Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018), who found that random sweeping indeed extends from single- to two-point quantities. (Strictly speaking, the random sweeping effect is reflected in the correlation between $\mathcal {A}_t^{q^2}$ and $(\mathcal {T}^{q^2}-\varPi ^{q^2})$. However, in the present analysis, we restrict to correlations between individual terms but remark that the signature of random sweeping is consistent with the results of Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018).) As $\mathcal {T}$ is associated with transport in physical space, its high correlation with $\mathcal {A}_t^{q^2}$ is likely a result of the known preferential anti-alignment between the unsteady acceleration $({\partial }/{\partial t}) u_i$ and the convective acceleration $u_k ({\partial }/{\partial x_k}) u_i$ (Tsinober Reference Tsinober2001).

Turning attention to the nonlinear transfer $\varPi ^{q^2}$, an anti-correlation with the pressure transport term $\mathcal {T}_p^{q^2}$ is seen to persist across all scales investigated here. As the pressure (and resulting pressure transport) results from a volumetric integration of the velocity field over the entire flow domain, the non-local influence of the baffles of the VK tank and the numerical forcing of the JHU data could non-trivially impact the role of the pressure transport. Despite this, a significant anti-correlation between $\mathcal {T}_p^{q^2}$ and $\varPi ^{q^2}$ persists between both data sets, indicating a dynamical link between those quantities.

As seen in figure 5, the correlations between the terms involved in the scale to scale helicity budget show some similarities with those observed for the energy budget. At the smallest scales, the same positive correlation between $\mathcal {D}_{\nu }^h$ and $\mathcal {E}^h$ is seen just as in the energy, reflecting the destruction of helicity by viscous forces. The signature of the random sweeping effect is again identified from the correlation between $\mathcal {A}_t^{h}$ and $\mathcal {T}^{h}$ at all $r$. In contrast to the energy transfers, however, the pressure transport $\mathcal {T}_p^{h}$ appears to have an increasing correlation with the nonlinear transfer $\varPi ^h$ as $r$ decreases.

Before this result is discussed for the helicity cascade in more depth, it is necessary to adopt the framework presented by Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) for the energy budget. As in their work, the correlation of $\mathcal {T}_p^{q^2}$ and $|\delta \boldsymbol {u}|\, |\delta \boldsymbol {\nabla } p |$ is confirmed to be negligible in the present study. On the other hand, defining the cosine of the angle between $\delta \boldsymbol u$ and $-\delta {\boldsymbol {\nabla }} p$ as

(4.1)\begin{equation} \cos{\phi_u} = (\widehat{\delta\boldsymbol{u}}\boldsymbol{\cdot}\hat{\boldsymbol{r}})(-\widehat{\delta\boldsymbol{\nabla} p}\boldsymbol{\cdot}\hat{\boldsymbol{r}}) + (\hat{\delta\boldsymbol{u}}\times\hat{\boldsymbol{r}})\boldsymbol{\cdot}(-\widehat{\delta\boldsymbol{\nabla} p}\times\hat{\boldsymbol{r}}), \end{equation}

where $\hat {\boldsymbol {\cdot }}$ denotes unit norm, we find that the correlation of $\mathcal {T}_p^{q^2}$ with $\cos {\phi _u}$ is substantial. This is shown in figure 6(a) for the JHU data, but similar results were found for the VK data (these are omitted for brevity).

Figure 6. Mean value of the alignment between unit vectors of the velocity (a) and vorticity (b) increments via ((4.1)) and ((4.2)). The correlations of the nonlinear helicity transfer mechanisms $\varPi ^{h,u}$ and $\varPi ^{h,\omega }$ (see text for definitions) with the pressure transport of helicity $\mathcal {T}_p^h$ are shown in (c). All results are from the JHU data.

From both the dot product and cross product contributions of (4.1), Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) found that the mean was small but positive. This implies that in the averaged picture, the convergence events (i.e. ${\delta \boldsymbol {u}}\boldsymbol {\cdot }\hat {\boldsymbol {r}}<0$) occur simultaneously with compressing forces (i.e. $-\widehat {\delta \boldsymbol {\nabla } p}\boldsymbol {\cdot }\hat {\boldsymbol {r}}<0$), and vice-versa: divergence events are coupled to expanding pressure forces. Importantly, one may conclude based on the identity between the (scale-space) volume integral $\unicode{x222D} _{\mathcal {V}_R} \varPi ^{q^2} \,\mathrm {d}\mathcal {V}_R$ and surface integral $\unicode{x222F} _{\partial \mathcal {V}_R} \delta \boldsymbol {u}\boldsymbol {\cdot }\hat {\boldsymbol {r}}\delta q^2\,\mathrm {d}\mathcal {S}$ that the convergence events facilitate the forward $\delta q^2$ cascade, whereas divergence events facilitate the inverse cascade (Yasuda & Vassilicos Reference Yasuda and Vassilicos2018). And further, it can be inferred from the positive value of mean $\cos {\phi _u}$ that the forward cascading of $\delta q^2$ coincides also with compressing forces.

Returning to the correlation of $\varPi ^h$ with $\mathcal {T}_p^h$ in the present study, it is apparent that the role of the pressure transport is comparatively more complex. This is because $\varPi ^h$ features two distinct NLF structure functions arising from the vortex twisting and stretching mechanisms (Yan et al. Reference Yan, Li, Yu, Wang and Chen2020). The mechanisms are separated here such that $\varPi ^h=\varPi ^{h,u}-\varPi ^{h,\omega }$, where $\varPi ^{h,u}=({\partial }/{\partial r_k})\delta u_k \delta h$ (vortex twisting) and $\varPi ^{h,\omega }=\tfrac {1}{2}({\partial }/{\partial r_k})\delta \omega _k \delta q^2$ (vortex stretching). As a near-zero correlation of $\mathcal {T}_p^h$ with $|\delta \boldsymbol {\omega }| |\delta \boldsymbol {\nabla } p |$ was found, an analogous analysis to Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) may be attempted using the angle defined via

(4.2)\begin{equation} \cos{\phi_\omega} = (\widehat{\delta\boldsymbol{\omega}}\boldsymbol{\cdot}\hat{\boldsymbol{r}})(-\widehat{\delta\boldsymbol{\nabla} p}\boldsymbol{\cdot}\hat{\boldsymbol{r}}) + (\widehat{\delta\boldsymbol{\omega}}\times\hat{\boldsymbol{r}})\boldsymbol{\cdot}(-\widehat{\delta\boldsymbol{\nabla} p}\times\hat{\boldsymbol{r}}). \end{equation}

The alignment for both the dot product and cross product contributions were similarly found to be small but positive, as can be seen in figure 6(b). The meaning of this alignment, however, is not trivial in the case of helicity. Such a correlation is expected to vanish in mirror-symmetric flow for lack of a preference for any specific sense of rotation (imposed by the sign of $\widehat {\delta \boldsymbol {\omega }}\boldsymbol {\cdot }\hat {\boldsymbol {r}}$). The non-zero average of $\cos {\phi _\omega }$, increasing with $r$, is reminiscent of the scale-to-scale helicity itself (figure 2b).

To reconcile the meaning of this result, we shall clarify the direction of the helicity cascade in the following. In contrast to $\delta q^2$, the terms ‘forward cascade’ and ‘inverse cascade’ are not well-established in the context of helicity (Alexakis Reference Alexakis2017), due to the lack of positive definiteness of $\delta h$. To address this, an auxiliary quantity $H_R$ is introduced as a reference for the sense of rotation as

(4.3)\begin{equation} H_R = \unicode{x222D}_{\mathcal{V}_R} \underbrace{(\delta\boldsymbol{u}\boldsymbol{\cdot}\hat{\boldsymbol{r}})(\delta\boldsymbol{\omega}\boldsymbol{\cdot}\hat{\boldsymbol{r}})}_{\delta h_R} \,\mathrm{d}\mathcal{V}_R. \end{equation}

The forward and backward directions of the helicity cascade is then defined based on the sign of $H_R$ (i.e. ${H_R}/{\lvert H_R\rvert })$. The helicity cascading events causing $({H_R}/{\lvert H_R\rvert })\unicode{x222D} _{\mathcal {V}_R} \varPi ^{h} \,\mathrm {d}\mathcal {V}_R < 0$ are deemed to cascade forward, while $({H_R}/{\lvert H_R\rvert })\unicode{x222D} _{\mathcal {V}_R} \varPi ^{h} \mathrm {d}\mathcal {V}_R> 0$ implies inverse cascading. Similarly to $\delta q^2$, analogous relations invoking the divergence theorem may be drawn for the nonlinear transfer terms:

(4.4)\begin{equation} \unicode{x222D}_{\mathcal{V}_R} \varPi^{h,u}\,\mathrm{d}\mathcal{V}_R = \unicode{x222F}_{\partial\mathcal{V}_R} \delta \boldsymbol{u}\boldsymbol{\cdot}\hat{\boldsymbol{r}}\delta h\,\mathrm{d}\mathcal{S}\quad {\rm and}\quad \unicode{x222D}_{\mathcal{V}_R} \varPi^{h,\omega}\,\mathrm{d}\mathcal{V}_R = \frac{1}{2}\unicode{x222F}_{\partial\mathcal{V}_R} \delta \boldsymbol{\omega}\boldsymbol{\cdot}\hat{\boldsymbol{r}}\delta q^2\,\mathrm{d}\mathcal{S}. \end{equation}

It follows that events of $({\delta h H_R}/{\lvert \delta h\rvert \lvert H_R\rvert })\delta \boldsymbol {u}\boldsymbol {\cdot }\hat {\boldsymbol {r}} < 0$ and $({H_R}/{\lvert H_R\rvert })\delta \boldsymbol {\omega }\boldsymbol {\cdot }\hat {\boldsymbol {r}} > 0$ facilitate the forward helicity cascade and vice-versa for the inverse cascade.

It is important to note that the two mechanisms (i.e. stretching and twisting) cascade helicity in the same direction (either upscale or downscale) only when $\delta h \delta h_r <0$. Additionally, a positive mean value of $\cos {\phi _\omega }$ (as in figure 6b) suggests that when the signs of $H_R$ and $\delta \boldsymbol {\omega }\boldsymbol {\cdot }\hat {\boldsymbol {r}}$ are in agreement, the downscale helicity cascade couples with expanding pressure forces and the upscale helicity cascade associates with compressing pressure forces (in the averaged sense).

The correlations of the individual vortex twisting and stretching mechanisms with $\mathcal {T}_p^h$ in figure 6(c) shows a clear distinction between mechanisms. In the former, the correlation is relatively steady within the probed separation range, similar to what is seen in the context of $\delta q^2$. On the other hand, the stretching mechanism exhibits strong interactions with pressure forces (stronger than the twisting) at dissipative scales, which diminish almost entirely by the time $r$ reaches $\lambda$.

Finally, we turn out attention to the correlations of the transfer terms of the scale-to-scale enstrophy are presented in figure 7. As explained in Appendix A, noise propagation in the VK data renders correlations of $\mathcal {A}_t^{\omega ^2}$ and $\mathcal {E}^{\omega ^2}$ particularly erroneous. The results of the VK data are therefore omitted.

Figure 7. Correlation coefficients between terms of (2.8) for $\delta \omega ^2$ at various separations of JHU data (ac).

A striking similarity with the energy and helicity budgets is immediately apparent through the anti-correlation of the unsteady transport $\mathcal {A}_t^{\omega ^2}$ and turbulent transport $\mathcal {T}^{\omega ^2}$. This confirms the random sweeping effect persists across all budgets and across all scales. The generation term $\mathcal {G}_s^{\omega ^2}$ exhibits a positive correlation with $\mathcal {E}^{\omega ^2}$ and a negative correlation with $\mathcal {D}_{\nu }^{\omega ^2}$. This implies a concurrent enstrophy generation and dissipation mechanism (Davidson et al. Reference Davidson, Morishita and Kaneda2008). This simultaneous generation and dissipation is tied to a reduction (during high dissipative events) or increase (during low dissipative events) in the diffusion of enstrophy. This is seen to persist across all scales considered. A small positive correlation between $\mathcal {D}_{\nu }^{\omega ^2}$ and $\mathcal {E}^{\omega ^2}$, akin to the energy and helicity cascades, is seen only for $r=2.5\eta$.

4.2. Correlations between budgets

We now explore possible connections between budget equations, i.e. correlations between the mechanisms involved in the budget of a given scale-by-scale structure function and those of the others, as shown in figure 8. The correlations between both the scale energy and enstrophy transfers and the helicity transfers were found to be universally of order $10^{-2}$ and are thus omitted. At first glance, this appears to contradict a wealth of research that has identified clear causal relationships between the three quantities (Bershadskii et al. Reference Bershadskii, Kit, Tsinober and Vaisburd1994; Biferale, Musacchio & Toschi Reference Biferale, Musacchio and Toschi2013; Alexakis Reference Alexakis2017; Bos Reference Bos2021), even analytically (via Lagrangian closure theory, see Inagaki Reference Inagaki2021). As discussed by Tsinober (Reference Tsinober2001), near-zero correlations are necessary, but not sufficient, to determine that two quantities are unrelated. For example, the unsteady term ${\partial u_i}/{\partial t}$ (2.1) is almost entirely decorrelated with the material acceleration ${{\rm D} u_i}/{{\rm D}t}$ simply due to the underlying anti-correlation of ${\partial u_i}/{\partial t}$ with $u_k ({\partial u_i}/{\partial x_k})$.

Figure 8. Correlation coefficients between terms of (2.4) and (2.8) (respectively transfer terms of $\delta q^2$ and $\delta \omega ^2$ on the horizontal and vertical axis) at various separations of JHU data (ac).

In the present case, the root of the decorrelation is likely to be the fact that helicity is not positive-definite (Kurien et al. Reference Kurien, Taylor and Matsumoto2004). Alternatively, a geometric approach allows to separate the transfer quantities into the product of a (positive-definite) norm with an appropriately defined angle (Yasuda & Vassilicos Reference Yasuda and Vassilicos2018). This avenue for further insight is explored in § 5.

Let us now return to figure 8, which shows correlations between energy and enstrophy budgets. The most robust correlation is that between dissipation rate of energy $\mathcal {E}^{q^2}$ and the dissipation rate of enstrophy $\mathcal {E}^{\omega ^2}$, with a magnitude exceeding 0.5 across all considered scales. Figure 8 also shows that the energy dissipation $\mathcal {E}^{q^2}$ and enstrophy diffusion $\mathcal {D}_{\nu }^{\omega ^2}$ are negatively correlated, whereas the energy dissipation $\mathcal {E}^{q^2}$ and enstrophy generation $\mathcal {G}_s^{\omega ^2}$ (for $r=0.5\lambda$ and $r=\lambda$) are positively correlated. Together, these observations support the classical picture that the generation of scale enstrophy at small scales (Siggia Reference Siggia1981) is largely responsible for the dissipation of energy (as well as enstrophy generation, Siggia Reference Siggia1981). Finally, the signature of the random sweeping effect can be clearly seen amongst the correlations between $\mathcal {A}_t$ and $\mathcal {T}$ in figure 8, particularly at the smallest scales. Together, these mechanisms point to the interconnected nature of the simultaneous turbulence exchanges.

5. Influence of scale-space helicity on the energy cascade

5.1. Application of Lamb decomposition

In the previous section, it was found that the correlation between transfers of energy and transfers of helicity was zero. There is however evidence that the magnitude of the helicity itself has an impact on the local transfer of energy. Regions of high helicity are thought to suppress nonlinear energy transfer and dissipation (Pelz et al. Reference Pelz, Yakhot, Orszag, Shtilman and Levich1985; Stepanov et al. Reference Stepanov, Golbraikh, Frick and Shestakov2015), though some results suggest otherwise (Zhou et al. Reference Zhou, Nagata, Sakai, Ito and Hayase2016). To probe more closely the connection between the helicity and the energy cascade in the present study, the Lamb decomposition is invoked on the nonlinear term of (2.1):

(5.1)\begin{equation} u_k \frac{\partial}{\partial x_k}u_i = \frac{1}{2}\frac{\partial}{\partial x_i}(u_k u_k)+\epsilon_{jki}\omega_j u_k, \end{equation}

where $\epsilon _{jki}\omega _j u_k={\boldsymbol \omega }\times {\boldsymbol u}$ is known as the Lamb vector. This quantity is central to the nonlinear turbulence cascade as maximal helicity $h={\boldsymbol \omega }\boldsymbol {\cdot }{\boldsymbol u}$, i.e. the alignment of ${\boldsymbol u}$ and ${\boldsymbol \omega }$, corresponds to a suppression of nonlinearity ${\boldsymbol \omega }\times {\boldsymbol u}$ in the Navier–Stokes equation.

A kinematic connection can be gleaned from (2.1), invoking (5.1). Similarly to the derivation of the KHMH energy (§ 2) but focusing on the nonlinear terms exclusively, it follows that

(5.2)\begin{align} &2\delta u_i \left(u_k \frac{\partial}{\partial x_k}u_i-u_k' \frac{\partial}{\partial x_k'}u_i'\right) ={-}2\varSigma u_k \epsilon_{jik}\delta \omega_j\delta u_i \nonumber\\ &\quad + \frac{\partial}{\partial r_k}(\delta q^2\delta u_k)-\frac{1}{2}\frac{\partial}{\partial X_k}(\varSigma u_k\delta q^2)+2\frac{\partial}{\partial r_i}(\delta u_i \varSigma q^2). \end{align}

After simplification, this leads to the triple product relation

(5.3)\begin{equation} \varSigma {\boldsymbol u}\boldsymbol{\cdot}(\delta{\boldsymbol \omega}\times\delta{\boldsymbol u})=\varPi^{q^2>}-\mathcal{T}^{q^2<}, \end{equation}

where $\varPi ^{q^2>}=({\partial }/{\partial r_i})(\delta u_i \varSigma q^2)$ is the nonlinear scale-to-scale transfer for the energy sum $\varSigma q^2=\varSigma u_i \varSigma u_i$, and $\mathcal {T}^{q^2<}=\tfrac {1}{2}({\partial }/{\partial X_k})(\varSigma u_k\delta q^2)$ is the usual nonlinear energy transport from (2.4). Hereafter the superscripts < and > are used for shorthand to refer to transfer of the increment and the sum quantities, respectively. A detailed derivation of this result is provided in Appendix B.

An immediate curiosity stems from the presence of the energy sum nonlinear transfer $\varPi ^{q^2>}$. The dynamical significance of the energy (or velocity) sum is a matter of some debate but has been linked to so-called ‘large-scale intermittency’ in previous studies (Sreenivasan & Antonia Reference Sreenivasan and Antonia1997; Blum et al. Reference Blum, Kunwar, Johnson and Voth2010; Chien, Blum & Voth Reference Chien, Blum and Voth2013; Carter & Coletti Reference Carter and Coletti2018). To summarise, instances when the energy sum is large and within the range of correlated scales (large-scale sweeping motions) have been shown to increase the content of the energy structure functions relative to the mean at each scale. This implies a connection between large and small scales in turbulence (Mininni et al. Reference Mininni, Alexakis and Pouquet2006; Hosokawa Reference Hosokawa2007). Equation (5.3) indicates that these sweeping motions play a central role in nonlinear energy transport. It shows that sufficiently large $\varSigma {\boldsymbol u}$ can overcome the alignment of $\delta {\boldsymbol u}$ and $\delta {\boldsymbol \omega }$ (as long as they not perfectly aligned). This could explain some conflicting results where strong dissipation is found even when these quantities are aligned and helicity is maximal (Zhou et al. Reference Zhou, Nagata, Sakai, Ito and Hayase2016).

Particularly for turbulence theory, the kinematic connection revealed by (5.2) provides novel insight into the mechanics of the energy cascade. By substituting for $\mathcal {T}^{q^2<}$ in (2.4) using (5.2) and invoking steady, homogeneous turbulence within an independent intermediate range of scales (Kolmogorov Reference Kolmogorov1941a,Reference Kolmogorovb), it is found upon ensemble averaging that

(5.4)\begin{equation} \overline{\varSigma {\boldsymbol u}\boldsymbol{\cdot}(\delta{\boldsymbol \omega}\times\delta{\boldsymbol u})}=4\overline{\epsilon^{q^2}}. \end{equation}

By virtue of homogeneity, the sum $\overline {\varPi ^{q^2>}}+\overline {\varPi ^{q^2<}}=0$ and therefore no longer features in the balance with the dissipation rate of energy. This can be shown analytically and was verified for both data sets in the present study (not shown for brevity). This demonstrates the high relevance of the triple product on the left-hand side of (5.4). It is not only a feature of the instantaneous kinematics but also central to the overall dynamics. It is perhaps unwelcome to see the velocity sum play a central role in the mean rate of energy dissipation in the inertial range, as this challenges the possibility that there exists a range of scales independent of the large-scale effects (or of the forcing). Such an observation has been previously identified in terms of triadic interactions (Yeung & Brasseur Reference Yeung and Brasseur1991), but to the authors’ knowledge this is the first such demonstration in real space. This provides direct insight on the role of large-scale intermittency (via $\varSigma {\boldsymbol u}$) on the inertial range in the cascade of energy (Sreenivasan & Antonia Reference Sreenivasan and Antonia1997).

Equation (5.4) highlights a geometric perspective of the energy dissipation, where the mean projection of the velocity sum onto the cross product of the velocity and vorticity increment balances the mean dissipation rate. By virtue of the fact that $\|\delta \boldsymbol {u}\times \delta \boldsymbol {\omega }\|^2=\|\delta \boldsymbol {u}\|^2\|\delta \boldsymbol {\omega }\|^2 - \|\delta h\|^2$, there is an implicit connection between the dissipation of energy and the scale-by-scale helicity. Although this does not necessarily connect the transfers, or dynamics, of the two quantities together, it does show that the magnitude of the scale-space helicity influences the dissipation of energy (but also, critically, depends on the alignment with sweeping motions).

Having identified the statistical significance of the triple product, we turn our attention to the instantaneous behaviour. Motivated by the kinematics of (5.3), the joint p.d.f.s of the triple product with the nonlinear energy sum transfer $\varPi ^{q^2>}$ and the nonlinear energy transport $\mathcal {T}^{q^2<}$ were computed for the VK and JHU data sets and plotted in figure 9. There is a clear positive and negative correlation between the triple product and $\varPi ^{q^2>}$ and $\mathcal {T}^{q^2<}$, respectively, with very close agreement in terms of the slope ($\pm 0.5$) of the major axis of the ellipse formed by the joint p.d.f.s. In fact, (5.2) requires the slopes be complementary. The slope, and therefore the underlying correlation, is robust across all scales considered and for both the experimental and numerical data sets. The increase in intermittency (i.e. the flatness of the joint p.d.f.s) with increasing scale is reminiscent of large-scale intermittency (Sreenivasan & Antonia Reference Sreenivasan and Antonia1997).

Figure 9. Joint p.d.f.s of the triple product of (5.3) with the nonlinear energy sum transfer $\varPi ^{q^2>}$ (a,c,e) and nonlinear energy transport $\mathcal {T}^{q^2<}$ (b,df) for $r=2.5\eta$ (a,b), $r=0.5\lambda$ (c,d) and $r=\lambda$ (ef). Contour values range from $10^{-12}$ to $10^{-8}$ (counting from the outermost line).

It is of interest that the triple product reveals a similar correlation to both $\varPi ^{q^2>}$ and $\mathcal {T}^{q^2<}$ separately. This originates with the velocity sum $\varSigma {\boldsymbol u}$, featuring both in the triple product and via $\mathcal {T}^{q^2<}=\varSigma {\boldsymbol u}\boldsymbol {\cdot }({\partial }/{\partial {\boldsymbol X}})({\delta q^2}/{2})$. This establishes an immediate and scale-independent connection between the two quantities. The remaining correlation between the triple product and $\varPi ^{q^2>}$ follows as a direct consequence of (5.3). Despite their closely tied kinematics, there are important dynamic distinctions amongst the quantities involved in terms of their statistical behaviour. Under the assumption of homogeneity, the ensemble average $\overline {\mathcal {T}^{q^2<}}=0$, whereas $\overline {\varPi ^{q^2>}}=-\overline {\varPi ^{q^2<}}$, and $\overline {\varSigma {\boldsymbol u}\boldsymbol {\cdot }(\delta {\boldsymbol \omega }\times \delta {\boldsymbol u})}=4\overline {\epsilon ^{q^2}}$. This motivates the interpretation of $\varPi ^{q^2>}$ as a transfer of energy to larger scales whose average is a complement to $\overline {\varPi ^{q^2<}}$.

The establishment of an overall balance between the triple product and the dissipation of energy in the inertial range (5.4) motivates further insight into the alignment of $\varSigma {\boldsymbol u}$ and $(\delta {\boldsymbol \omega }\times \delta {\boldsymbol u})$. This is quantified using the angle defined by

(5.5)\begin{equation} \cos\theta=\frac{\varSigma {\boldsymbol u}\boldsymbol{\cdot}(\delta{\boldsymbol \omega}\times\delta{\boldsymbol u})}{\|\varSigma {\boldsymbol u}\|\,\|\delta {\boldsymbol u}\times\delta{\boldsymbol \omega}\|}. \end{equation}

The p.d.f. is presented for both data sets in figure 10. Note that an extra separation at $r=4.7\lambda$ (the maximum separation attainable using the sub-cubes extracted from the JHU database) is presented for further insight. As required by (5.4) to maintain dissipation downscale, the alignment is found to be biased towards positive values of $\cos \theta$ for all separations such that $\overline {\cos \theta }>0$. The distributions appear approximately quadratic in nature, with an inflection occurring somewhere between $2.5\eta$ and $0.5\lambda$.

Figure 10. P.d.f. of $\cos \theta$ characterising the alignment between the velocity sum and the cross product (5.3) for the VK data (a) and JHU data (b), and an extrapolation of the p.d.f. at extreme values of $\cos \theta$ for JHU data (c) ($L$ marks the integral length scale).

As the separation increases, the probability of parallel ($\cos \theta =1$) and anti-parallel ($\cos \theta =-1$) alignment is increasingly likely. Most evident from the JHU distributions (figure 10b), the increase in likelihood across separations is more rapid for anti-parallel alignment than parallel alignment. A linear regression for these points with increasing scale indicates they will be equal in probability at $r\approx 13\lambda$, corresponding almost exactly to the integral scale of the DNS data set (figure 10c). This implies that at $r=L$, the p.d.f. will be symmetric about zero and, due to the parabolic shape, with equal likelihood of parallel and anti-parallel alignment, a minimum at $\cos \theta =0$ and $\overline {\cos \theta }=0$. A similar extrapolation for the experimental data was not conclusive due to the limited number of available separations.

6. Conclusions

A scale-space framework, already well established for generalisation of the turbulent energy transfers via the KHMH equation (Hill Reference Hill2002; Valente & Vassilicos Reference Valente and Vassilicos2015), has been extended analogously to the scale-space helicity and enstrophy together for the first time. This has provided a novel framework for insight on the instantaneous turbulence dynamics that are responsible for establishing the statistically steady state. To exemplify this framework, an analysis of the various inter-scale and inter-space transfers of energy, helicity, and enstrophy in two homogeneous turbulent flows focusing on scales up to the Taylor microscale has been presented.

The correlation-based analysis has revealed that the random sweeping effect is indeed present within all three considered budgets. This was identified through the significant anti-correlation of the unsteady and physical-space transport terms ($\mathcal {A}_t$ and $\mathcal {T}$, respectively) within the individual budgets. A significant correlation of these transfers between budgets was also identified via the energy and enstrophy, suggesting that random sweeping also occurs simultaneously across budgets. In addition, the energy and enstrophy transfers were also found to correlate through the generation/dissipation of enstrophy and dissipation of energy, consistent with previous studies (Davidson et al. Reference Davidson, Morishita and Kaneda2008). This again points to the connection of mechanisms across budgets.

The anti-correlation between the pressure transport and nonlinear inter-scale transfer terms ($\varPi$ and $\mathcal {T}_p$, respectively), already reported by Yasuda & Vassilicos (Reference Yasuda and Vassilicos2018) in the context of the energy cascade, was also observed in the helicity cascade. In contrast to the energy cascade, this anti-correlation was strongest at the dissipative scales, where both vortex-stretching and twisting mechanisms correlate equally. Overall (i.e. on average) it was found that the downscale helicity cascade couples with expanding pressure forces, and the upscale helicity cascade associates with compressing pressure forces.

As a direct consequence of the non-positive definiteness of helicity, it was found that the terms governing the transfers of helicity were entirely decorrelated from those governing the energy and enstrophy. Despite this, it was shown that the scale-to-scale helicity itself influences the energy transfers. On the basis of kinematic arguments, it was demonstrated that the alignment between the two-point velocity sum and the cross product of the velocity and vorticity increments is connected to the overall inter-scale transfer of energy (5.4). This is at odds with the theoretical independence of small-scale dynamics from large-scale motions but is consistent with so-called large-scale intermittency phenomena observed in numerous studies (e.g. Sreenivasan & Antonia Reference Sreenivasan and Antonia1997; Blum et al. Reference Blum, Kunwar, Johnson and Voth2010; Carter & Coletti Reference Carter and Coletti2018). It is further seen that, to maintain the dissipation of energy, this alignment is non-zero and positive on average.

Taken together, the present work illustrates how the the extended scale-space framework (introduced here in its most general form) can be used to probe the various mechanisms involved in the concurrent exchanges of energy, enstrophy and helicity. In this first iteration, this framework was applied to homogeneous turbulence, allowing us to identify links and similarities between the various transfers. Furthermore, by dissecting the term associated with inter-scale energy exchanges, we find that the scale helicity regulates the imbalance between inter-scale exchanges of large-scale energy ($\varSigma q^2$) and the inter-space exchanges of small-scale energy ($\delta q^2$) through $\varSigma {\boldsymbol u}\boldsymbol {\cdot }(\delta {\boldsymbol \omega }\times \delta {\boldsymbol u})$.

Future applications of the present framework to other canonical problems, such as a Taylor–Green vortex, but also (and perhaps more importantly) to flows where the turbulence is markedly inhomogeneous, such as wall-bounded and free-shear flows, are warranted. A detailed analysis of the scale-space triple product, particularly as it pertains to the large scale sweeping motions (e.g. a conditional analysis for when $u_k$ and $u_k'$ are aligned/anti-aligned) is also low-hanging fruit. Finally, the use of physics-driven instantaneous correlations (as in the present study) as parameters in the loss functions of neural networks for sub-grid stress models in LES is a promising avenue for future research.

Acknowledgements

The authors wish to acknowledge J. Dawson for kindly sharing experimental data. D.W.C. wishes to thank F. Coletti, A. Alexakis and B. Ganapathisubramani for support in early stages. This work is the product of exchanges of ideas at the Eleventh International Symposium on Turbulence and Shear Flow Phenomena in Southampton, UK.

Declaration of interests

The authors report no conflict of interest.

Data availability

Data for all figures are available from the University of Southampton repository at https://doi.org/10.5258/SOTON/D2420.

Appendix A. Energy, enstrophy and helicity transfer budgets

The extensive size of the data sets makes tabulating all data points at each location in scale space into a single storage variable prohibitively expensive. For this reason, a limited number of separations were chosen in order to store the relevant information in confined memory.

For each volumetric cube of Cartesian data and for each separation, spheres corresponding to scale-space radius $r$ were positioned at all possible centroid locations on the grid. As the number of possible centroid locations decreases with increasing sphere radius, the number of samples used for calculating statistics decreases as $r$ increases (Camussi et al. Reference Camussi, Baudet, Benzi and Ciliberto1996), leading to increased uncertainty at larger separations (e.g. see the error bars of figure 3). Statistics (i.e. the first and second moments of the quantities of interest) were calculated by averaging over the sphere. As the grid points did not provide a uniform probing of the sphere, a weighting procedure was adopted to account for the non-uniformity. The number of surface points on the sphere was chosen to be $N_s=96$, with no significant impact on the results tested with a variety of $N_s$ between 24 and 288. The correlation between two quantities, $\zeta$ and $\xi$, is obtained as

(A1)\begin{equation} {\rm corr}(\zeta,\xi)=\frac{\overline{\zeta\xi}-\bar{\zeta}\bar{\xi}}{\sqrt{\overline{\zeta^{2}}-{\bar{\zeta}}^{2}}\sqrt{\overline{\xi^{2}}-\bar{\xi}^{2}}}. \end{equation}

By normalising by the product of the square root of the variance of both quantities, values of 1, 0 and $-1$ correspond to perfectly correlated, uncorrelated and anti-correlated quantities, respectively.

The histogram for each sphere at each separation and for all quantities of interest was tabulated and stored in 2500 linearly spaced bins. Similarly, the joint histogram for each pair of investigated quantities was stored in a $500\times 500$ bin matrix. The large number of bins allowed to capture the most extreme values of the highly intermittent quantities. As these extreme values are not known a priori, an adaptive histogram approach was used in which the bins and bin counts were adjusted when a value exceeding the current limits of the bins was encountered. Once tabulated, the histograms were normalised into p.d.f.s (or joint p.d.f.s).

To validate the adopted approach, the budget equations representing the balance of the various inter-scale transfers are presented for the VK and JHU data in figure 11. Similar results are obtained for both data sets. The energy budgets exhibit the classical behaviour: they are dominated by terms $\varPi ^{q^2}$ and $\mathcal {E}^{q^2}$ at larger separations, while $\mathcal {D}_{\nu }^{q^2}$ and $\mathcal {E}^{q^2}$ are the most significant at the smallest separation. The accuracy of the method and overall statistics is also supported by the small residuals in each case. However, it must be noted that the particularly small residuals associated with the VK data are largely artificial. The post-processing routine applied to the scanning PIV data yielded the acceleration and pressure fields in addition to the measured velocity fields. This extra information was resolved by minimization of the instantaneous momentum equation residual (i.e. a data assimilation approach). Therefore, small values of the budgets residuals are enforced in the case of VK data.

Figure 11. Budgets of terms of (2.4) (a,b), (2.6) (c,d) and (2.8) (ef) evaluated for VK (a,c,e) and JHU (b,df) datasets (res. stands for residuals of the respective budgets).

Although there is a generally good match between budgets based on the two data sets, this is certainly not the case when it comes to the unsteady term of the enstrophy budget $\mathcal {A}_t^{\omega ^2}$ and enstrophy dissipation rate $\mathcal {E}^{\omega ^2}$. The former term ought to vanish in a statistically stationary flow, such as the two considered in this manuscript, yet it dominates the VK enstrophy budget (for the JHU case $\mathcal {A}_t^{\omega ^2}$ is indeed small and of the order of the residual). The dissipation term $\mathcal {E}^{\omega ^2}$, on the other hand, is twice as large in the VK case than in the JHU data. This mismatch is due to bias error propagation via discrete differentiation. This effect is expected to be the strongest for the enstrophy transfer terms, owing to the high order derivatives involved (acting as noise amplifiers). For example, the measured value of $\mathcal {E}^{\omega ^2}$ is proportional to $(\partial \omega _i/\partial x_j)^2_{meas.} = (\partial \omega _i/\partial x_j)^2_{true}+2e(\partial \omega _i/\partial x_j)_{true}+e^2$, where $e$ represents error. The non-negative term $e^2$ is the dominating source of bias error (a similar example in the context of the energy dissipation rate is given by Tanaka & Eaton Reference Tanaka and Eaton2007). Since the VK data is the noisier of the two sets, it follows that the measured value of $\mathcal {E}^{\omega ^2}$ is larger, in agreement with figure 11.

Taking into account the above considerations, any conclusions regarding terms $\mathcal {A}_t^{\omega ^2}$ and $\mathcal {E}^{\omega ^2}$ based on the VK data are questionable. This ambiguity is avoided in the manuscript by discarding the VK data and relying purely on the JHU data for the enstrophy transfers.

Appendix B. Derivation of triple product relation

We start by introducing shorthand notation for the KHMH equation for the energy structure function and the energy sum, respectively:

(B1)\begin{align} &\overbrace{\frac{\partial}{\partial t}\delta q^2 + \frac{\partial \delta u_k\delta q^2 }{\partial r_k} ={-}\frac{\partial}{\partial X_k}\left(\frac{\varSigma u_k\delta q^2}{2}\right)- \frac{2}{\rho}\frac{\partial \delta u_k \delta p}{\partial X_k}}^{\mathcal{A}_t^{q^2<}+\varPi^{q^2<}={-}\mathcal{T}^{q^2<}-\mathcal{T}_p^{q^2<}+\mathcal{D}_{\nu}^{q^2<}-\mathcal{E}^{q^2}} \nonumber\\ &\quad +\nu\left[2\frac{\partial^2}{\partial r_k^2}+\frac{1}{2}\frac{\partial^2}{\partial X_k^2}\right]\delta q^2 -2\nu\left[\left(\frac{\partial u_i}{\partial x_k}\right)^2+\left(\frac{\partial u_i'}{\partial x_k'}\right)^2 \right], \end{align}
(B2)\begin{align} &\overbrace{\frac{\partial}{\partial t}\varSigma q^2 + \frac{\partial \delta u_k\varSigma q^2 }{\partial r_k} ={-}\frac{\partial}{\partial X_k}\left(\frac{\varSigma u_k\varSigma q^2}{2}\right)- \frac{2}{\rho}\frac{\partial \varSigma u_k \delta p}{\partial X_k}}^{\mathcal{A}_t^{q^2>}+\varPi^{q^2>}={-}\mathcal{T}^{q^2>}-\mathcal{T}_p^{q^2>}+\mathcal{D}_{\nu}^{q^2>}-\mathcal{E}^{q^2}} \nonumber\\ &\quad+\nu\left[2\frac{\partial^2}{\partial r_k^2}+\frac{1}{2}\frac{\partial^2}{\partial X_k^2}\right]\varSigma q^2 -2\nu\left[\left(\frac{\partial u_i}{\partial x_k}\right)^2+\left(\frac{\partial u_i'}{\partial x_k'}\right)^2 \right]. \end{align}

Simplification starts with the Lamb decomposition applied to the nonlinear terms. This can begin with the nonlinear terms of either (B1) or (B2) and leads to the same result:

(B3)\begin{align} &\overbrace{2\delta u_i\left(u_j\frac{\partial}{\partial x_j}u_i-u_j'\frac{\partial}{\partial x_j'}u_i'\right)}^{\mathcal{T}^{q^2<}+\varPi^{q^2<}}\nonumber\\ &\quad =2\delta u_i \left(\epsilon_{kji}\omega_k u_j - \epsilon_{kji}\omega_k'u_j' + \frac{1}{2}\frac{\partial}{\partial x_i}(u_j u_j)-\frac{1}{2}\frac{\partial}{\partial x_i'}(u_j' u_j')\right) \nonumber\\ &\quad =2\delta u_i\epsilon_{kji}\left(\delta \omega_k \varSigma u_j + \omega_k'u_j - \omega_k u_j'\right) + 2\delta u_i\left( \frac{1}{2}\frac{\partial}{\partial x_i}(u_j u_j)-\frac{1}{2}\frac{\partial}{\partial x_i'}(u_j' u_j')\right) \nonumber\\ &\quad ={-}2\varSigma u_j \epsilon_{kij}\delta \omega_k \delta u_i + 2\delta u_i \epsilon_{kji}(\omega_k'u_j - \omega_k u_j') \nonumber\\ &\qquad + 2\delta u_i\left( \frac{1}{2}\frac{\partial}{\partial x_i}(u_j u_j)-\frac{1}{2}\frac{\partial}{\partial x_i'}(u_j' u_j')\right) \nonumber\\ &\quad ={-}2\varSigma u_j \epsilon_{kij}\delta \omega_k \delta u_i + 2\delta u_i\epsilon_{kji}\epsilon_{mlk}\left(u_j \frac{\partial}{\partial x_m'}u_l'-u_j'\frac{\partial}{\partial x_m}u_l\right) \nonumber\\ &\qquad + 2\delta u_i\left( \frac{1}{2}\frac{\partial}{\partial x_i}(u_j u_j)-\frac{1}{2}\frac{\partial}{\partial x_i'}(u_j' u_j')\right) \nonumber\\ &\quad =-2\varSigma u_j \epsilon_{kij}\delta \omega_k \delta u_i + 2\delta u_i(\delta_{jm}\delta_{il}-\delta_{jl}\delta_{im})\left(u_j \frac{\partial}{\partial x_m'}u_l'-u_j'\frac{\partial}{\partial x_m}u_l\right) \nonumber\\ &\qquad + 2\delta u_i\left( \frac{1}{2}\frac{\partial}{\partial x_i}(u_j u_j)-\frac{1}{2}\frac{\partial}{\partial x_i'}(u_j' u_j')\right) \nonumber\\ &\quad =-2\varSigma u_j \epsilon_{kij}\delta \omega_k \delta u_i + 2\delta u_i\left(u_j\frac{\partial}{\partial x_j'}u_i'-u_j'\frac{\partial}{\partial x_j}u_i\right) \nonumber\\ &\qquad + 2\delta u_i\left({-}u_j\frac{\partial}{\partial x_i'}u_j' + u_j'\frac{\partial}{\partial x_i}u_j+ \frac{1}{2}\frac{\partial}{\partial x_i}(u_j u_j)-\frac{1}{2}\frac{\partial}{\partial x_i'}(u_j' u_j')\right) \nonumber\\ &\quad =-2\varSigma u_j \epsilon_{kij}\delta \omega_k \delta u_i + 2\delta u_i\left({-}u_j\frac{\partial}{\partial x_j'}\delta u_i -u_j'\frac{\partial}{\partial x_j}\delta u_i\right) + \delta u_i\left( \frac{\partial}{\partial x_i}\varSigma q^2 - \frac{\partial}{\partial x_i'}\varSigma q^2\right) \nonumber\\ &\quad = \overbrace{-2\varSigma u_j\epsilon_{kij}\delta \omega_k\delta u_i}^{{-}2\varSigma{\boldsymbol u}\boldsymbol{\cdot}(\delta{\boldsymbol \omega}\times\delta{\boldsymbol u})} + \overbrace{\frac{\partial}{\partial r_j}(\delta u_j\delta q^2)}^{+\varPi^{q^2<}} -\overbrace{\frac{1}{2}\frac{\partial}{\partial X_j}(\varSigma u_j \delta q^2)}^{-\mathcal{T}^{q^2<}}+\overbrace{2\frac{\partial}{\partial r_j}(\delta u_j \varSigma q^2)}^{{+}2\varPi^{q^2>}} \end{align}

The final expression in shorthand notation is

(B4)\begin{equation} \mathcal{T}^{q^2<}+\varPi^{q^2<} ={-}2\varSigma{\boldsymbol u}\boldsymbol{\cdot}(\delta{\boldsymbol \omega}\times\delta{\boldsymbol u})+\varPi^{q^2<}-\mathcal{T}^{q^2<} +2\varPi^{q^2>}. \end{equation}

Subtracting $\varPi ^{q^2<}$ from both sides, dividing by two and rearranging, we arrive at (5.3).

References

REFERENCES

Alexakis, A. 2017 Helically decomposed turbulence. J.Fluid Mech. 812, 752770.CrossRefGoogle Scholar
Alexakis, A. & Biferale, L. 2018 Cascades and transitions in turbulent flows. Phys. Rep. 767, 1101.CrossRefGoogle Scholar
Alexakis, A. & Chibbaro, S. 2020 Local energy flux of turbulent flows. Phys. Rev. Fluids 5, 094604.CrossRefGoogle Scholar
Alves Portela, F., Papadakis, G. & Vassilicos, J.C. 2017 The turbulence cascade in the near wake of a square prism. J.Fluid Mech. 825, 315352.CrossRefGoogle Scholar
Arun, S., Sameen, A., Srinivasan, B. & Girimaji, S.S. 2021 Scale-space energy density function transport equation for compressible inhomogeneous turbulent flows. J.Fluid Mech. 920, A31.CrossRefGoogle Scholar
Bershadskii, A., Kit, E., Tsinober, A. & Vaisburd, H. 1994 Strongly localized events of energy, dissipation, enstrophy and enstrophy generation in turbulent flows. Fluid Dyn. Res. 14 (2), 71.CrossRefGoogle Scholar
Biferale, L., Musacchio, S. & Toschi, F. 2013 Split energy-helicity cascades in three-dimensional homogeneous and isotropic turbulence. J. Fluid Mech. 730, 309327.CrossRefGoogle Scholar
Blum, D.B., Kunwar, S.B., Johnson, J. & Voth, G.A. 2010 Effects of nonuniversal large scales on conditional structure functions in turbulence. Phys. Fluids 22 (1), 015107.CrossRefGoogle Scholar
Boffetta, G. & Ecke, R.E. 2012 Two-dimensional turbulence. Annu. Rev. Fluid Mech. 44, 427451.CrossRefGoogle Scholar
Bos, W.J.T. 2021 Three-dimensional turbulence without vortex stretching. J.Fluid Mech. 915, A21.CrossRefGoogle Scholar
Brissaud, A., Frisch, U., Léorat, J., Lesieur, M. & Mazure, A. 1973 Helicity cascades in fully developed isotropic turbulence. Phys. Fluids 16 (8), 13661367.CrossRefGoogle Scholar
Camussi, R., Baudet, C., Benzi, R. & Ciliberto, S. 1996 Statistical uncertainty in the analysis of structure functions in turbulence. Phys. Rev. E 54 (4), R3098.CrossRefGoogle ScholarPubMed
Carbone, M. & Bragg, A.D. 2020 Is vortex stretching the main cause of the turbulent energy cascade? J.Fluid Mech. 883, R2.CrossRefGoogle Scholar
Carter, D.W. & Coletti, F. 2018 Small-scale structure and energy transfer in homogeneous turbulence. J.Fluid Mech. 854, 505543.CrossRefGoogle Scholar
Chen, Q., Chen, S. & Eyink, G.L. 2003 The joint cascade of energy and helicity in three-dimensional turbulence. Phys. Fluids 15 (2), 361374.CrossRefGoogle Scholar
Chien, C.-C., Blum, D.B. & Voth, G.A. 2013 Effects of fluctuating energy input on the small scales in turbulence. J.Fluid Mech. 737, 527551.CrossRefGoogle Scholar
Chkhetiani, O.G. 1996 On the third moments in helical turbulence. J.Expl Theor. Phys. Lett. 63 (10), 808812.CrossRefGoogle Scholar
Cimarelli, A., De Angelis, E., Jimenez, J. & Casciola, C.M. 2016 Cascades and wall-normal fluxes in turbulent channel flows. J.Fluid Mech. 796, 417436.CrossRefGoogle Scholar
Danaila, L., Krawczynski, J.F., Thiesset, F. & Renou, B. 2012 Yaglom-like equation in axisymmetric anisotropic turbulence. Physica D 241 (3), 216223.CrossRefGoogle Scholar
Dascaliuc, R. & Grujić, Z. 2013 Coherent vortex structures and 3d enstrophy cascade. Commun. Math. Phys. 317 (2), 547561.CrossRefGoogle Scholar
Davidson, P.A., Morishita, K. & Kaneda, Y. 2008 On the generation and flux of enstrophy in isotropic turbulence. J.Turbul. 9, N42.CrossRefGoogle Scholar
Davidson, P.A. & Pearson, B.R. 2005 Identifying turbulent energy distributions in real, rather than fourier, space. Phys. Rev. Lett. 95 (21), 214501.CrossRefGoogle ScholarPubMed
Deusebio, E., Boffetta, G., Lindborg, E. & Musacchio, S. 2014 Dimensional transition in rotating turbulence. Phys. Rev. E 90 (2), 023005.CrossRefGoogle ScholarPubMed
Elsinga, G.E., Ishihara, T., Goudar, M.V., Da Silva, C.B. & Hunt, J.C.R. 2017 The scaling of straining motions in homogeneous isotropic turbulence. J.Fluid Mech. 829, 3164.CrossRefGoogle Scholar
Eyink, G.L. 1996 Exact results on stationary turbulence in 2d: consequences of vorticity conservation. Physica D 91 (1–2), 97142.CrossRefGoogle Scholar
Falkovich, G. 2009 Symmetries of the turbulent state. J.Phys. A 42 (12), 123001.CrossRefGoogle Scholar
Falkovich, G., Boffetta, G., Shats, M. & Lanotte, A.S. 2017 Introduction to focus issue: two-dimensional turbulence. Phys. Fluids 29 (11), 110901.CrossRefGoogle Scholar
Fjørtoft, R. 1953 On the changes in the spectral distribution of kinetic energy for twodimensional, nondivergent flow. Tellus 5 (3), 225230.CrossRefGoogle Scholar
Gatti, D., Chiarini, A., Cimarelli, A. & Quadrio, M. 2020 Structure function tensor equations in inhomogeneous turbulence. J.Fluid Mech. 898, A5.CrossRefGoogle Scholar
Germano, M., Piomelli, U., Moin, P. & Cabot, W.H. 1991 A dynamic subgrid-scale eddy viscosity model. Phys. Fluids A 3 (7), 17601765.CrossRefGoogle Scholar
Gomes-Fernandes, R., Ganapathisubramani, B. & Vassilicos, J.C. 2015 The energy cascade in near-field non-homogeneous non-isotropic turbulence. J.Fluid Mech. 771, 676705.CrossRefGoogle Scholar
Goto, S. & Vassilicos, J.C. 2016 Unsteady turbulence cascades. Phys. Rev. E 94 (5), 053108.CrossRefGoogle ScholarPubMed
Hill, R.J. 2002 Exact second-order structure-function relationships. J.Fluid Mech. 468, 317326.CrossRefGoogle Scholar
Hosokawa, I. 2007 A paradox concerning the refined similarity hypothesis of kolmogorov for isotropic turbulence. Prog. Theor. Phys. 118 (1), 169173.CrossRefGoogle Scholar
Inagaki, K. 2021 Scale-similar structures of homogeneous isotropic non-mirror-symmetric turbulence based on the lagrangian closure theory. J.Fluid Mech. 926, A14.CrossRefGoogle Scholar
Ishihara, T., Gotoh, T. & Kaneda, Y. 2009 Study of high–Reynolds number isotropic turbulence by direct numerical simulation. Annu. Rev. Fluid Mech. 41, 165180.CrossRefGoogle Scholar
Ishihara, T., Kaneda, Y. & Hunt, J.C.R. 2013 Thin shear layers in high Reynolds number turbulence–DNS results. Flow Turbul. Combust. 91 (4), 895929.CrossRefGoogle Scholar
Jiménez, J., Wray, A.A., Saffman, P.G. & Rogallo, R.S. 1993 The structure of intense vorticity in isotropic turbulence. J.Fluid Mech. 255, 6590.CrossRefGoogle Scholar
van Kan, A. & Alexakis, A. 2022 Energy cascades in rapidly rotating and stratified turbulence within elongated domains. J.Fluid Mech. 933, A11.CrossRefGoogle Scholar
Kivotides, D. & Leonard, A. 2021 Helicity spectra and topological dynamics of vortex links at high Reynolds numbers. J.Fluid Mech. 911, A25.CrossRefGoogle Scholar
Kleckner, D., Kauffman, L.H. & Irvine, W.T.M. 2016 How superfluid vortex knots untie. Nat. Phys. 12 (7), 650655.CrossRefGoogle Scholar
Knutsen, A.N., Baj, P., Lawson, J.M., Bodenschatz, E., Dawson, J.R. & Worth, N.A. 2020 The inter-scale energy budget in a von kármán mixing flow. J.Fluid Mech. 895, A11.CrossRefGoogle Scholar
Kolmogorov, A.N. 1941 a Dissipation of energy in locally isotropic turbulence. In Dokl. Akad. Nauk SSSR, vol. 32, pp. 16–18.Google Scholar
Kolmogorov, A.N. 1941 b The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers. In Dokl. Akad. Nauk SSSR, vol. 30, pp. 299–303.Google Scholar
Kraichnan, R.H. 1967 Inertial ranges in two-dimensional turbulence. Phys. Fluids 10 (7), 14171423.CrossRefGoogle Scholar
Kraichnan, R.H. 1973 Helical turbulence and absolute equilibrium. J.Fluid Mech. 59 (4), 745752.CrossRefGoogle Scholar
Kraichnan, R.H. 1974 On kolmogorov's inertial-range theories. J.Fluid Mech. 62 (2), 305330.CrossRefGoogle Scholar
Kurien, S., Taylor, M.A. & Matsumoto, T. 2004 Isotropic third-order statistics in turbulence with helicity: the 2/15-law. J. Fluid Mech. 515, 8797.CrossRefGoogle Scholar
Lai, C.C.K., Charonko, J.J. & Prestridge, K. 2018 A Kármán–Howarth–Monin equation for variable-density turbulence. J.Fluid Mech. 843, 382418.CrossRefGoogle Scholar
Laval, J.P., Dubrulle, B. & Nazarenko, S. 2001 Nonlocality and intermittency in three-dimensional turbulence. Phys. Fluids 13 (7), 19952012.CrossRefGoogle Scholar
Lawson, J.M. & Dawson, J.R. 2014 A scanning PIV method for fine-scale turbulence measurements. Exp. Fluids 55 (12), 1857.CrossRefGoogle Scholar
Lawson, J.M. & Dawson, J.R. 2015 On velocity gradient dynamics and turbulent structure. J.Fluid Mech. 780, 6098.CrossRefGoogle Scholar
Levich, E. & Shtilman, L. 1988 Coherence and large fluctuations of helicity in homogeneous turbulence. Phys. Lett. A 126 (4), 243248.CrossRefGoogle Scholar
Li, Y., Perlman, E., Wan, M., Yang, Y., Meneveau, C., Burns, R., Chen, S., Szalay, A. & Eyink, G. 2008 A public turbulence database cluster and applications to study lagrangian evolution of velocity increments in turbulence. J.Turbul. 9, N31.CrossRefGoogle Scholar
Lindborg, E. 2007 Third-order structure function relations for quasi-geostrophic turbulence. J.Fluid Mech. 572, 255260.CrossRefGoogle Scholar
Lindborg, E. & Nordmark, A. 2022 Two-dimensional turbulence on a sphere. J.Fluid Mech. 933, A60.CrossRefGoogle Scholar
L'vov, V.S., Podivilov, E. & Procaccia, I. 1997 Exact result for the 3rd order correlations of velocity in turbulence with helicity. arXiv:chao-dyn/9705016.Google Scholar
Marati, N., Casciola, C.M. & Piva, R. 2004 Energy cascade and spatial fluxes in wall turbulence. J.Fluid Mech. 521, 191215.CrossRefGoogle Scholar
Mininni, P.D., Alexakis, A. & Pouquet, A. 2006 Large-scale flow effects, energy transfer, and self-similarity on turbulence. Phys. Rev. E 74 (1), 016303.CrossRefGoogle ScholarPubMed
Moffatt, H.K. 1969 The degree of knottedness of tangled vortex lines. J.Fluid Mech. 35 (1), 117129.CrossRefGoogle Scholar
Moffatt, H.K. & Tsinober, A. 1992 Helicity in laminar and turbulent flow. Annu. Rev. Fluid Mech. 24 (1), 281312.CrossRefGoogle Scholar
Mollicone, J.-P., Battista, F., Gualtieri, P. & Casciola, C.M. 2018 Turbulence dynamics in separated flows: the generalised kolmogorov equation for inhomogeneous anisotropic conditions. J.Fluid Mech. 841, 10121039.CrossRefGoogle Scholar
Monin, A.S. & Yaglom, A.M. 1975 Statistical Fluid Mechanics (ed. J. Lumley). MIT Press.Google Scholar
Moreau, J.J. 1960 Constantes d'un îlot tourbillonnaire en fluide parfait barotrope. C. R. Hebd. Seances Acad. Sci. 252, 28102812.Google Scholar
Novara, M. & Scarano, F. 2013 A particle-tracking approach for accurate material derivative measurements with tomographic PIV. Exp. Fluids 54 (8), 1584.CrossRefGoogle Scholar
Pelz, R.B., Yakhot, V., Orszag, S.A., Shtilman, L. & Levich, E. 1985 Velocity-vorticity patterns in turbulent flow. Phys. Rev. Lett. 54 (23), 2505.CrossRefGoogle ScholarPubMed
Pope, S.B. 2001 Turbulent Flows, Cambridge University Press.CrossRefGoogle Scholar
Pouquet, A., Rosenberg, D., Stawarz, J.E. & Marino, R. 2019 Helicity dynamics, inverse, and bidirectional cascades in fluid and magnetohydrodynamic turbulence: a brief review. Earth Space Sci. 6 (3), 351369.CrossRefGoogle Scholar
Raynal, F. 1996 Exact relation between spatial mean enstrophy and dissipation in confined incompressible flows. Phys. Fluids 8 (8), 22422244.CrossRefGoogle Scholar
Richardson, L.F. 1920 The supply of energy from and to atmospheric eddies. Proc. R. Soc. Lond. A 97 (686), 354373.Google Scholar
Scott, R.B. & Wang, F. 2005 Direct evidence of an oceanic inverse kinetic energy cascade from satellite altimetry. J.Phys. Oceanogr. 35 (9), 16501666.CrossRefGoogle Scholar
Siggia, E.D. 1981 Numerical study of small-scale intermittency in three-dimensional turbulence. J.Fluid Mech. 107, 375406.CrossRefGoogle Scholar
Sreenivasan, K.R. 1995 On the universality of the kolmogorov constant. Phys. Fluids 7 (11), 27782784.CrossRefGoogle Scholar
Sreenivasan, K.R. & Antonia, R.A. 1997 The phenomenology of small-scale turbulence. Annu. Rev. Fluid Mech. 29 (1), 435472.CrossRefGoogle Scholar
Stepanov, R., Golbraikh, E., Frick, P. & Shestakov, A. 2015 Hindered energy cascade in highly helical isotropic turbulence. Phys. Rev. Lett. 115 (23), 234501.CrossRefGoogle ScholarPubMed
Tanaka, T. & Eaton, J.K. 2007 A correction method for measuring turbulence kinetic energy dissipation rate by PIV. Exp. Fluids 42 (6), 893902.CrossRefGoogle Scholar
Tennekes, H. & Lumley, J.L. 2018 A First Course in Turbulence. MIT Press.Google Scholar
Tsinober, A. 2001 An Informal Introduction to Turbulence, vol. 63. Springer Science & Business Media.Google Scholar
Valente, P.C. & Vassilicos, J.C. 2015 The energy cascade in grid-generated non-equilibrium decaying turbulence. Phys. Fluids 27 (4), 045103.CrossRefGoogle Scholar
Van Atta, C.W. & Antonia, R.A. 1980 Reynolds number dependence of skewness and flatness factors of turbulent velocity derivatives. Phys. Fluids 23 (2), 252257.CrossRefGoogle Scholar
Von Kármán, T. & Howarth, L. 1938 On the statistical theory of isotropic turbulence. Proc. R. Soc. Lond. A 164 (917), 192215.CrossRefGoogle Scholar
Waleffe, F. 1992 The nature of triad interactions in homogeneous turbulence. Phys. Fluids A 4 (2), 350363.CrossRefGoogle Scholar
Wang, C.Y., Gao, Q., Wei, R.J., Li, T. & Wang, J.J. 2017 Weighted divergence correction scheme and its fast implementation. Exp. Fluids 58 (5), 44.CrossRefGoogle Scholar
Yan, Z., Li, X., Yu, C., Wang, J. & Chen, S. 2020 Dual channels of helicity cascade in turbulent flows. J.Fluid Mech. 894, R2.CrossRefGoogle Scholar
Yasuda, T. & Vassilicos, J.C. 2018 Spatio-temporal intermittency of the turbulent energy cascade. J.Fluid Mech. 853, 235252.CrossRefGoogle Scholar
Yeung, P.K. & Brasseur, J.G. 1991 The response of isotropic turbulence to isotropic and anisotropic forcing at the large scales. Phys. Fluids A 3 (5), 884897.CrossRefGoogle Scholar
Yokoi, N. & Yoshizawa, A. 1993 Statistical analysis of the effects of helicity in inhomogeneous turbulence. Phys. Fluids A 5 (2), 464477.CrossRefGoogle Scholar
Zhou, Y. 2021 Turbulence theories and statistical closure approaches. Phys. Rep. 935, 1117.CrossRefGoogle Scholar
Zhou, Y., Nagata, K., Sakai, Y., Ito, Y. & Hayase, T. 2016 Spatial evolution of the helical behavior and the 2/3 power-law in single-square-grid-generated turbulence. Fluid Dyn. Res. 48 (2), 021404.CrossRefGoogle Scholar
Zimmerman, S.J., Antonia, R.A., Djenidi, L., Philip, J. & Klewicki, J.C. 2022 Approach to the 4/3 law for turbulent pipe and channel flows examined through a reformulated scale-by-scale energy budget. J.Fluid Mech. 931, A28.CrossRefGoogle Scholar
Figure 0

Table 1. Parameters from the Von-Kármán mixing tank scanning PIV experiment (Knutsen et al.2020) and the Johns Hopkins University DNS of homogeneous and isotropic turbulence (Li et al.2008). Note the symbols $\Delta x$ and $\Delta t$ denote the grid spacing and temporal spacing of the data sets.

Figure 1

Figure 1. Premulitplied p.d.f.s of the dissipation of energy (a), helicity (b) and enstrophy (c). The inset of (b) highlights the range of helicity dissipation between ${\pm }5$ standard deviations.

Figure 2

Figure 2. Normalised mean energy structure function $\overline {\delta q^2}$ (a), helicity structure function $\bar {\delta h}$ (b) and enstrophy structure function $\overline {\delta \omega ^2}$ (c).

Figure 3

Figure 3. Normalised orientation-averaged NLF structure functions corresponding to the nonlinear scale-to-scale flux of energy $\overline {\delta q^2}$ (a), helicity $\bar {\delta h}$ via advection (b) and vortex stretching (c), and enstrophy $\overline {\delta \omega ^2}$ (d).

Figure 4

Figure 4. Correlation coefficients between terms of (2.4) for $\delta q^2$ at various separations; results based on VK (ac) and JHU (df) datasets.

Figure 5

Figure 5. Correlation coefficients between terms of (2.6) for $\delta h$ at various separations; results based on VK (ac) and JHU (df) datasets.

Figure 6

Figure 6. Mean value of the alignment between unit vectors of the velocity (a) and vorticity (b) increments via ((4.1)) and ((4.2)). The correlations of the nonlinear helicity transfer mechanisms $\varPi ^{h,u}$ and $\varPi ^{h,\omega }$ (see text for definitions) with the pressure transport of helicity $\mathcal {T}_p^h$ are shown in (c). All results are from the JHU data.

Figure 7

Figure 7. Correlation coefficients between terms of (2.8) for $\delta \omega ^2$ at various separations of JHU data (ac).

Figure 8

Figure 8. Correlation coefficients between terms of (2.4) and (2.8) (respectively transfer terms of $\delta q^2$ and $\delta \omega ^2$ on the horizontal and vertical axis) at various separations of JHU data (ac).

Figure 9

Figure 9. Joint p.d.f.s of the triple product of (5.3) with the nonlinear energy sum transfer $\varPi ^{q^2>}$ (a,c,e) and nonlinear energy transport $\mathcal {T}^{q^2<}$ (b,df) for $r=2.5\eta$ (a,b), $r=0.5\lambda$ (c,d) and $r=\lambda$ (ef). Contour values range from $10^{-12}$ to $10^{-8}$ (counting from the outermost line).

Figure 10

Figure 10. P.d.f. of $\cos \theta$ characterising the alignment between the velocity sum and the cross product (5.3) for the VK data (a) and JHU data (b), and an extrapolation of the p.d.f. at extreme values of $\cos \theta$ for JHU data (c) ($L$ marks the integral length scale).

Figure 11

Figure 11. Budgets of terms of (2.4) (a,b), (2.6) (c,d) and (2.8) (ef) evaluated for VK (a,c,e) and JHU (b,df) datasets (res. stands for residuals of the respective budgets).