Quantifying wall turbulence via a symmetry approach. Part 2. Reynolds stresses

• Published in: Journal of fluid mechanics Vol. 850; pp. 401 - 438
• Main Authors: Chen, Xi, Hussain, Fazle, She, Zhen-Su
• Format: Journal Article
• Language: English
• Published: Cambridge, UK Cambridge University Press 10.09.2018
• Subjects: Computational fluid dynamics; Computer simulation; Data; Direct numerical simulation; Domain names; Educational institutions; Fluid flow; Formulae; Formulas (mathematics); Hypotheses; JFM Papers; Mathematical models; Navier-Stokes equations; Pipes; Reynolds number; Reynolds stress; Reynolds stresses; Scaling; Shear flow; Shear stress; Symmetry; Turbulence; Turbulent flow; Turbulent shear flow; Velocity distribution; China; turbulent boundary layers; pipe flow boundary layer; turbulence theory
• ISSN: 0022-1120; 1469-7645
• DOI: 10.1017/jfm.2018.405

We present new scaling expressions, including high-Reynolds-number ( $Re$ ) predictions, for all Reynolds stress components in the entire flow domain of turbulent channel and pipe flows. In Part 1 (She et al., J. Fluid Mech., vol. 827, 2017, pp. 322–356), based on the dilation symmetry of the mean Navier–Stokes equation a four-layer formula of the Reynolds shear stress length $\ell _{12}$ – and hence also the entire mean velocity profile (MVP) – was obtained. Here, random dilations on the second-order balance equations for all the Reynolds stresses (shear stress $-\overline{u^{\prime }v^{\prime }}$ , and normal stresses $\overline{u^{\prime }u^{\prime }}$ , $\overline{v^{\prime }v^{\prime }}$ , $\overline{w^{\prime }w^{\prime }}$ ) are analysed layer by layer, and similar four-layer formulae of the corresponding stress length functions $\ell _{11}$ , $\ell _{22}$ , $\ell _{33}$ (hence the three turbulence intensities) are obtained for turbulent channel and pipe flows. In particular, direct numerical simulation (DNS) data are shown to agree well with the four-layer formulae for $\ell _{12}$ and $\ell _{22}$ – which have the celebrated linear scalings in the logarithmic layer, i.e. $\ell _{12}\approx \unicode[STIX]{x1D705}y$ and $\ell _{22}\approx \unicode[STIX]{x1D705}_{22}y$ . However, data show an invariant peak location for $\overline{w^{\prime }w^{\prime }}$ , which theoretically leads to an anomalous scaling in $\ell _{33}$ in the log layer only, namely $\ell _{33}\propto y^{1-\unicode[STIX]{x1D6FE}}$ with $\unicode[STIX]{x1D6FE}\approx 0.07$ . Furthermore, another mesolayer modification of $\ell _{11}$ yields the experimentally observed location and magnitude of the outer peak of $\overline{u^{\prime }u^{\prime }}$ . The resulting $-\overline{u^{\prime }v^{\prime }}$ , $\overline{u^{\prime }u^{\prime }}$ , $\overline{v^{\prime }v^{\prime }}$ and $\overline{w^{\prime }w^{\prime }}$ are all in good agreement with DNS and experimental data in the entire flow domain. Our additional results include: (1) the maximum turbulent production is located at $y^{+}\approx 12$ ; (2) the location of peak value $-\overline{u^{\prime }v^{\prime }}_{p}$ has a scaling transition from $5.7Re_{\unicode[STIX]{x1D70F}}^{1/3}$ to $1.5Re_{\unicode[STIX]{x1D70F}}^{1/2}$ at $Re_{\unicode[STIX]{x1D70F}}\approx 3000$ , with a $1+\overline{u^{\prime }v^{\prime }}_{p}^{+}$ scaling transition from $8.5Re_{\unicode[STIX]{x1D70F}}^{-2/3}$ to $3.0Re_{\unicode[STIX]{x1D70F}}^{-1/2}$ ( $Re_{\unicode[STIX]{x1D70F}}$ the friction Reynolds number); (3) the peak value $\overline{w^{\prime }w^{\prime }}_{p}^{+}\approx 0.84Re_{\unicode[STIX]{x1D70F}}^{0.14}(1-48/Re_{\unicode[STIX]{x1D70F}})$ ; (4) the outer peak of $\overline{u^{\prime }u^{\prime }}$ emerges above $Re_{\unicode[STIX]{x1D70F}}\approx 10^{4}$ with its location scaling as $1.1Re_{\unicode[STIX]{x1D70F}}^{1/2}$ and its magnitude scaling as $2.8Re_{\unicode[STIX]{x1D70F}}^{0.09}$ ; (5) an alternative derivation of the log law of Townsend (1976, The Structure of Turbulent Shear Flow, Cambridge University Press), namely, $\overline{u^{\prime }u^{\prime }}^{+}\approx -1.25\ln y+1.63$ and $\overline{w^{\prime }w^{\prime }}^{+}\approx -0.41\ln y+1.00$ in the bulk.