The spectral model of Perry, Henbest & Chong (1986) predicts that the integral length-scale varies very slowly with distance to the wall in the intermediate layer. The only way for the integral length scales variation to be more realistic while keepi
ng with the Townsend-Perry attached eddy spectrum is to add a new wavenumber range to the model at wavenumbers smaller than that spectrum. This necessary addition can also account for the high Reynolds number outer peak of the turbulent kinetic energy in the intermediate layer. An analytic expression is obtained for this outer peak in agreement with extremely high Reynolds number data by Hultmark, Vallikivi, Bailey & Smits (2012, 2013). The finding of Dallas, Vassilicos & Hewitt (2009) that it is the eddy turnover time and not the mean flow gradient which scales with distance to the wall and skin friction velocity in the intermediate layer implies, when combined with Townsends (1976) production-dissipation balance, that the mean flow gradient has an outer peak at the same location as the turbulent kinetic energy. This is seen in the data of Hultmark, Vallikivi, Bailey & Smits (2012, 2013). The same approach also predicts that the mean flow gradient has a logarithmic decay at distances to the wall larger than the position of the outer peak, a qualitative prediction which the aforementioned data also support.
We focus in this paper on the effect of the resolution of Direct Numerical Simulations (DNS) on the spatio-temporal development of the turbulence downstream of a single square grid. The aims of this study are to validate our numerical approach by com
paring experimental and numerical one-point statistics downstream of a single square grid and then investigate how the resolution is impacting the dynamics of the flow. In particular, using the Q-R diagram, we focus on the interaction between the strain-rate and rotation tensors, the symmetric and skew-symmetric parts of the velocity gradient tensor respectively. We first show good agreement between our simulations and hot-wire experiment for one-point statistics on the centreline of the single square grid. Then, by analysing the shape of the Q-R diagram for various streamwise locations, we evaluate the ability of under-resolved DNS to capture the main features of the turbulence downstream of the single square grid.
We investigate non-equilibrium turbulence where the non-dimensionalised dissipation coefficient $C_{varepsilon}$ scales as $C_{varepsilon} sim Re_{M}^{m}/Re_{ell}^{n}$ with $mapprox 1 approx n$ ($Re_M$ and $Re_{ell}$ are global/inlet and local Reynol
ds numbers respectively) by measuring the downstream evolution of the scale-by-scale energy transfer, dissipation, advection, production and transport in the lee of a square-mesh grid and compare with a region of equilibrium turbulence (i.e. where $C_{varepsilon}approx mathrm{constant}$). These are the main terms of the inhomogeneous, anisotropic version of the von K{a}rm{a}n-Howarth-Monin equation. It is shown in the grid-generated turbulence studied here that, even in the presence of non-negligible turbulence production and transport, production and transport are large-scale phenomena that do not contribute to the scale-by-scale balance for scales smaller than about a third of the integral-length scale, $ell$, and therefore do not affect the energy transfer to the small-scales. In both the non-equilibrium and the equilibrium decay regions, the peak of the scale-by-scale energy transfer scales as $(overline{u^2})^{3/2}/ell$ ($overline{u^2}$ is the variance of the longitudinal fluctuating velocity). In the non-equilibrium case this scaling implies an imbalance between the energy transfer to the small scales and the dissipation. This imbalance is reflected on the small-scale advection which becomes larger in proportion to the maximum energy transfer as the turbulence decays whereas it stays proportionally constant in the further downstream equilibrium region where $C_{varepsilon} approx mathrm{constant}$ even though $Re_{ell}$ is lower.
The previously reported non-equilibrium dissipation law is investigated in turbulent flows generated by various regular and fractal square grids. The flows are documented in terms of various turbulent profiles which reveal their differences. In spite
of significant inhomogeneity and anisotropy differences, the new non-equilibrium dissipation law is observed in all these flows. Various transverse and longitudinal integral scales are measured and used to define the dissipation coefficient $C_{varepsilon}$. It is found that the new non-equilibrium dissipation law is not an artefact of a particular choice of the integral scale and that the usual equilibrium dissipation law can actually coexist with the non-equilibrium law in different regions of the same flow.
