Hostname: page-component-7b7567c8c-lpg2g Total loading time: 0 Render date: 2026-11-15T22:43:18.006Z Has data issue: false hasContentIssue false

Large-scale structures in a turbulent channel flow with a minimal streamwise flow unit

Published online by Cambridge University Press:  06 July 2018

Hiroyuki Abe*
Affiliation:
Japan Aerospace Exploration Agency, Tokyo 182-8522, Japan
Robert Anthony Antonia
Affiliation:
Discipline of Mechanical Engineering, School of Engineering, University of Newcastle, NSW 2308, Australia
Sadayoshi Toh
Affiliation:
Department of Physics and Astronomy, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
*
†Email address for correspondence: habe@chofu.jaxa.jp

Abstract

Direct numerical simulations are used to examine large-scale motions with a streamwise length $2\sim 4h$ ($h$ denotes the channel half-width) in the logarithmic and outer regions of a turbulent channel flow. We test a minimal ‘streamwise’ flow unit (Toh & Itano, J. Fluid Mech., vol. 524, 2005, pp. 249–262) (or MSU) for larger Kármán numbers ($h^{+}=395$ and 1020) than in the original work. This flow unit consists of a sufficiently long (${L_{x}}^{+}\approx 400$) streamwise domain to maintain near-wall turbulence (Jiménez & Moin, J. Fluid Mech., vol. 225, 1991, pp. 213–240) and a spanwise domain which is large enough to represent the spanwise behaviour of inner and outer structures correctly; as $h^{+}$ increases, the streamwise extent of the MSU domain decreases with respect to $h$. Particular attention is given to whether the spanwise organization of the large-scale structures may be represented properly in this simplified system at sufficiently large $h^{+}$ and how these structures are associated with the mean streamwise velocity $\overline{U}$. It is shown that, in the MSU, the large-scale structures become approximately two-dimensional at $h^{+}=1020$. In this case, the streamwise velocity fluctuation $u$ is energized, whereas the spanwise velocity fluctuation $w$ is weakened significantly. Indeed, there is a reduced energy redistribution arising from the impaired global nature of the pressure, which is linked to the reduced linear–nonlinear interaction in the Poisson equation (i.e. the rapid pressure). The logarithmic dependence of $\overline{ww}$ is also more evident due to the reduced large-scale spanwise meandering. On the other hand, the spanwise organization of the large-scale $u$ structures is essentially identical for the MSU and large streamwise domain (LSD). One discernible difference, relative to the LSD, is that the large-scale structures in the MSU are more energized in the outer region due to a reduced turbulent diffusion. In this region, there is a tight coupling between neighbouring structures, which yields antisymmetric pairs (with respect to centreline) of large-scale structures with a spanwise spacing of approximately $3h$; this is intrinsically identical with the outer energetic mode in the optimal transient growth of perturbations (del Álamo & Jiménez, J. Fluid Mech., vol. 561, 2006, pp. 329–358).

Information

Type
JFM Papers
Copyright
© 2018 Cambridge University Press 
Figure 0

Table 1.Domain size, grid points and spatial resolution for the DNS with the varying $L_{x}$.

Figure 1

Table 2.Domain size, grid points and spatial resolution of the DNS for the LSD.

Figure 2

Figure 1.Distributions of $\overline{U}$, $y(\unicode[STIX]{x2202}\overline{U}/\unicode[STIX]{x2202}y)$ and $-\overline{uv}$ normalized by wall units: (a,b) $\overline{U}^{+}$; (c) $y^{+}(\unicode[STIX]{x2202}\overline{U}^{+}/\unicode[STIX]{x2202}y^{+})$; (d) $-\overline{u^{+}v^{+}}$. In (a,b), the distribution of $\overline{U}^{+}$ is shown in the half and whole channels, respectively. In (d), the distribution of the total shear stress, $\unicode[STIX]{x1D70F}_{total}^{+}=\unicode[STIX]{x2202}\overline{U}^{+}/\unicode[STIX]{x2202}y^{+}-\overline{u^{+}v^{+}}$, is included. - - - - denotes the DNS data of del Álamo et al. (2004) for $h^{+}=934$ with the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$). ○ and ▫ refer to the experimental data of Wei & Willmarth (1989) for $h^{+}=1016$ and of Monty & Chong (2009) for $h^{+}=1040$.

Figure 3

Figure 2.Streamwise and spanwise two-point correlations of velocity fluctuations: (a) $R_{uu}(\unicode[STIX]{x0394}x)$, $R_{vv}(\unicode[STIX]{x0394}x)$ and $R_{ww}(\unicode[STIX]{x0394}x)$ for the MSU for $h^{+}=395$ and 1020 at $y^{+}=5$; (b) $R_{uu}(\unicode[STIX]{x0394}x)$, $R_{vv}(\unicode[STIX]{x0394}x)$ and $R_{ww}(\unicode[STIX]{x0394}x)$ for the MSU for $h^{+}=395$ and 1020 at $y/h=0.5$; (c) $R_{uu}(\unicode[STIX]{x0394}x)$ for the LSD for $h^{+}=1020$; (d) $R_{uu}(\unicode[STIX]{x0394}x)$ for the MSU for $h^{+}=1020$; (e) $R_{ww}(\unicode[STIX]{x0394}z)$ for the LSD for $h^{+}=1020$; (f) $R_{ww}(\unicode[STIX]{x0394}z)$ for the MSU for $h^{+}=1020$. In (c–f), the $y$ locations of the plotted data are $y^{+}=5$, 21, 61, 101 and $y/h=0.2$, 0.5, 1.

Figure 4

Figure 3.Distributions of normalized mean-square velocity fluctuations: (a,b) $\overline{uu}/U_{\unicode[STIX]{x1D70F}}^{2}$; (c,d) $\overline{vv}/U_{\unicode[STIX]{x1D70F}}^{2}$; (e,f) $\overline{ww}/U_{\unicode[STIX]{x1D70F}}^{2}$. The data for the LSD are plotted in (a,c,e) together with those for the VLSDs in the literature, whereas those for the MSU are given in (b,d,f). In (a,c,e), – ⋅ – ⋅ –, $h^{+}=180$; $\cdots \cdots$, $h^{+}=395$; - - - -, $h^{+}=640$; ——, $h^{+}=1020$; – – – –, $h^{+}=2003$ (Hoyas & Jiménez 2008); — ⋅ ⋅ —, $h^{+}=4079$ (Bernardini et al.2014); ○, Wei & Willmarth (1989) for $h^{+}=1016$. The bold dashed-lines in (a) and (e) come from relations of (3.3) and (3.4), respectively, with $A_{1}=2.1$, $A_{3}=0.85$, $B_{1}=1.26$, $B_{3}=0.45$, while those in (f) are from (3.4) with $A_{3}=0.15$, $B_{3}=0.45$. In (b,d,f), – ⋅ – ⋅ –, $h^{+}=137$ (Toh & Itano 2005); $\cdots \cdots$, $h^{+}=395$; ——, $h^{+}$$=1020$.

Figure 5

Figure 4.Contours in the $x$–$y$ plane of instantaneous $u$, $p$ and $w$ for $h^{+}=1020$: (a) $u$; (b) $p$; (c) $w$. In each panel, the top section shows the data for $L_{x}=12.8h$ (LSD), whereas the bottom section shows the data for $L_{x}=6.4h$, $L_{x}=3.2h$, $L_{x}=1.6h$, $L_{x}=0.8h$ and $L_{x}=0.4h$ (MSU), the order being from left to right. In the top panel, the black solid lines have inclinations of 10 (left line) and 20 (right line) degrees.

Figure 6

Figure 5.Effect of varying $L_{x}$ on premultiplied one-dimensional streamwise spectra of $u$, $k_{x}\unicode[STIX]{x1D719}_{uu}(k_{x})$, normalized by $U_{\unicode[STIX]{x1D70F}}^{2}$ for $h^{+}=1020$: (a) $y/h=0.1$; (b) $y/h=0.3$; (c) $y/h=0.5$; (d) $y/h=0.7$. – – – – denotes the DNS data of del Álamo et al. (2004) for $h^{+}=934$ with the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$). ○ and ▫ refer to the experimental data of Balakumar & Adrian (2007) for $h^{+}=960$ and of Monty & Chong (2009) for $h^{+}=1040$.

Figure 7

Figure 6.Effect of varying $L_{x}$ on streamwise two-point correlations of $u$, $R_{uu}(\unicode[STIX]{x0394}x)$, for $h^{+}=1020$: (a) $y/h=0.1$; (b) $y/h=0.3$; (c) $y/h=0.5$; (d) $y/h=0.7$.

Figure 8

Figure 7.Effect of varying $L_{x}$ on the mean-square velocity and pressure fluctuations for $h^{+}=1020$: (a) $\overline{uu}/U_{\unicode[STIX]{x1D70F}}^{2}$; (b) $\overline{ww}/U_{\unicode[STIX]{x1D70F}}^{2}$; (c) $\overline{vv}/U_{\unicode[STIX]{x1D70F}}^{2}$; (d) $\overline{pp}/\unicode[STIX]{x1D70F}_{w}^{2}$. - - - - denotes the DNS data of del Álamo et al. (2004) for $h^{+}=934$ with the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$). ▫ refers to the experimental data of Monty & Chong (2009) for $h^{+}=1040$. In (b,d), the bold dashed-lines are based on relations (3.4) with $A_{3}=0.85$, $B_{3}=0.45$ and (3.7), respectively.

Figure 9

Figure 8.Effect of varying $L_{x}$ on premultiplied one-dimensional streamwise and spanwise spectra of $w$ and $p$ for $h^{+}=1020$ at $y/h=0.5$: (a) $k_{x}\unicode[STIX]{x1D719}_{ww}(k_{x})/U_{\unicode[STIX]{x1D70F}}^{2}$; (b) $k_{z}\unicode[STIX]{x1D719}_{ww}(k_{z})/U_{\unicode[STIX]{x1D70F}}^{2}$; (c) $k_{x}\unicode[STIX]{x1D719}_{pp}(k_{x})/\unicode[STIX]{x1D70F}_{w}^{2}$; (d) $k_{z}\unicode[STIX]{x1D719}_{pp}(k_{z})/\unicode[STIX]{x1D70F}_{w}^{2}$. In (a,b), – – – – denotes the DNS data of del Álamo et al. (2004) for $h^{+}=934$ with the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$).

Figure 10

Figure 9.Effect of varying $L_{x}$ on budgets of the Reynolds normal stresses, energy dissipation rate, pressure strain and diffusion terms for $h^{+}=1020$: (a) budget of $\overline{uu}$; (b) budget of $\overline{vv}$; (c) budget of $\overline{ww}$; (d) $\overline{\unicode[STIX]{x1D700}}h/U_{\unicode[STIX]{x1D70F}}^{3}$; (e) $\unicode[STIX]{x1D719}_{11}^{+}$, $\unicode[STIX]{x1D719}_{22}^{+}$ and $\unicode[STIX]{x1D719}_{33}^{+}$; (f) $\overline{w^{+}w^{+}v^{+}}$ and $\overline{v^{+}p^{+}}$. In (a–c), the normalization has been made by $U_{\unicode[STIX]{x1D70F}}^{4}/\unicode[STIX]{x1D708}$. In (d), the bold dashed-line denotes the relation $\overline{\unicode[STIX]{x1D700}}h/U_{\unicode[STIX]{x1D70F}}^{3}=2.54/(y/h)-2.6$, which is equivalent to relation (3.9) or (3.10).

Figure 11

Figure 10.Contours in the $y$–$z$ plane of instantaneous streamwise, wall-normal and spanwise velocity fluctuations, $u$, $v$, $w$, for $h^{+}=1020$: (a,b) $u$; (c,d) $v$; (e,f) $w$. (a,c,e) Correspond to the LSD ($L_{x}=12.8h$), while (b,d,f) are for the MSU (${L_{x}}^{+}\approx 400$). In (a,b), lines denote positive values of the second invariant of $Q^{+}$ (line contour levels are from $5\times 10^{-4}$ to $5\times 10^{-3}$ with an increment of $5\times 10^{-5}$).

Figure 12

Figure 11.One-dimensional premultiplied spanwise spectra of $u$, $v$ and $w$ for $h^{+}=1020$: (a,b) $k_{z}\unicode[STIX]{x1D719}_{uu}(k_{z})/\overline{uu}$; (c,d) $k_{z}\unicode[STIX]{x1D719}_{vv}(k_{z})/\overline{vv}$; (e,f) $k_{z}\unicode[STIX]{x1D719}_{ww}(k_{z})/\overline{ww}$. (a,c,e) Correspond to the LSD ($L_{x}=12.8h$), while (b,d,f) are for the MSU (${L_{x}}^{+}\approx 400$). In each figure, the $y$ locations of the plotted data are $y^{+}=5$, 21, 61, 101 and $y/h=0.2$, 0.5, 1.

Figure 13

Figure 12.Most energetic spanwise scales, $\unicode[STIX]{x1D706}_{z,max}$, obtained from peaks of premultiplied spanwise spectra of $u$, $w$ and $p$ for the MSU and LSD $h^{+}=1020$: (a), $\unicode[STIX]{x1D706}_{z,max}/h$ for $u$ (b) $\unicode[STIX]{x1D706}_{z,max}^{+}$ for $u$; (c) $\unicode[STIX]{x1D706}_{z,max}/h$ for $w$; (d) $\unicode[STIX]{x1D706}_{z,max}/h$ for $p$. In (a,b), ▫ denotes the DNS data of del Álamo et al. (2004) for $h^{+}=934$ with the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$). In (a); the dashed, one-dotted and two-dotted lines denote $\unicode[STIX]{x1D706}_{z,max}/h=0.8$, 1.3 and 1.6, respectively; the ‘diamonds’ refer to the data of Tomkins & Adrian (2005) in a zero-pressure-gradient turbulent boundary layer for $\unicode[STIX]{x1D6FF}^{+}=426$. In (b–d), the dashed lines are linear fits to the data.

Figure 14

Figure 13.One-dimensional premultiplied spanwise spectra of $u$ for $h^{+}=1020$ for the LSD (the black line) and MSU (the red line) at $y/h=0.2$. The spectra are compared with those of del Álamo et al. (2004) and Hoyas & Jiménez (2008) for $h^{+}=934$ and 2003 at $y/h=0.16$ (the green line) and $y/h=0.2$ (the blue line), respectively, for the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$) in the present flow. Also included are the spectra of Tomkins & Adrian (2005) in a zero-pressure-gradient turbulent boundary layer for $\unicode[STIX]{x1D6FF}^{+}=2216$ at $y/\unicode[STIX]{x1D6FF}=0.2$ (circles) and those of Ahn et al. (2015) in a turbulent pipe for $R^{+}=3008$ at $y/R=0.1$ (triangles).

Figure 15

Figure 14.Time variation of premultiplied spanwise spectra of $u$ at $y^{+}=101$ for $h^{+}=1020$: (a) LSD ($L_{x}=12.8h$); (b) MSU (${L_{x}}^{+}\approx 400$). The spectra have been normalized by $U_{\unicode[STIX]{x1D70F}}^{2}$. In (a), —— and – ⋅ – ⋅ – correspond to the spectral peaks at $\unicode[STIX]{x1D706}_{z}/h=0.8$ and 1.3, respectively, for the LSD. In (b), – ⋅ – ⋅ – and — ⋅ ⋅ — correspond to the spectral peaks at $\unicode[STIX]{x1D706}_{z}/h=0.9$ and 1.6, respectively, for the MSU.

Figure 16

Figure 15.Comparison between x–z contours of $u$ for the LSD ($L_{x}=12.8h$) (a,b) and t–z contours of $u^{2D}$ for the MSU (c,d) for $h^{+}=1020$. The data at $y^{+}=15$ are plotted in (a,c), whereas those at $y/h=0.5$ are shown in (b,d). In (a,b), the instantaneous $u$ data in the LSD have been plotted at the two $y$ locations in the x–z plane. In (c,d), the time evolution data of $u^{2D}$, defined in (3.18), in the MSU have been plotted at the two $y$ locations as a function of the spanwise location.

Figure 17

Figure 16.Normalized one-dimensional spanwise spectra of velocity fluctuations: (a) $\unicode[STIX]{x1D719}_{uu}(k_{z})/U_{\unicode[STIX]{x1D70F}}\unicode[STIX]{x1D708}$, $\unicode[STIX]{x1D719}_{vv}(k_{z})/U_{\unicode[STIX]{x1D70F}}\unicode[STIX]{x1D708}$, $\unicode[STIX]{x1D719}_{ww}(k_{z})/U_{\unicode[STIX]{x1D70F}}\unicode[STIX]{x1D708}$ at $y^{+}=15$; (b) $\unicode[STIX]{x1D719}_{uu}(k_{z})/U_{\unicode[STIX]{x1D70F}}^{2}h$ at $y^{+}=15$ and $y/h=0.5$; (c) $\unicode[STIX]{x1D719}_{uu}(k_{z})/U_{\unicode[STIX]{x1D70F}}^{2}h$ for the LSD at $y^{+}=5$, 61 and $y/h=0.2$, 0.5, 0.8; (d) $\unicode[STIX]{x1D719}_{uu}(k_{z})/U_{\unicode[STIX]{x1D70F}}^{2}h$ for the MSU at $y^{+}=5$, 61 and $y/h=0.2$, 0.5, 0.8. In (a,b), the spectra for the MSU have been compared with those for the LSD. In (c,d), the spectra have been plotted in inner scales to highlight the low-wavenumber behaviour. In (a), – – – – denotes the DNS data of del Álamo et al. (2004) for $h^{+}=934$ with the VLSD ($L_{x}=8\unicode[STIX]{x03C0}h$).

Figure 18

Figure 17.Time variation of the large-scale structures in the MSU for $h^{+}=1020$: (a) $t^{+}=0$; (b) $t^{+}=600$; (c) $t^{+}=1200$; (d) $t^{+}=1800$: red, $u^{+}>1.5$; blue, $u^{+}<-1.5$. Spanwise scales below 0.5$h$ are filtered out with the use of a spectral sharp cutoff filter. The visualization domain size shown here is the whole computational domain for the MSU (i.e. $0.4h\times 2h\times 6.4h$ in the $x$, $y$ and $z$ directions, respectively).

Figure 19

Figure 18.Wall-normal two-point correlations for $h^{+}=1020$: (a) $R_{uu}(y_{r},y)$; (b) $R_{vv}(y_{r},y)$; (c) $R_{ww}(y_{r},y)$; (d) $R_{pp}(y_{r},y)$. In this figure, the averaging is based on a time period $t^{+}\approx 1300$ and 1800, for the LSD and MSU, respectively. The locations of $y_{r}$ are $y/h=0$, 0.2, 0.5 and 1. The data at $y_{r}=0$ have been obtained at the first grid point.

Figure 20

Figure 19.Contours in the $y$–$z$ plane of instantaneous total, slow rapid pressure for $h^{+}=1020$: (a,b) $p$; (c,d) $p_{s}$; (e,f) $p_{r}$. (a,c,e) correspond to the LSD ($L_{x}=12.8h$), while (b,d,f) are for the MSU (${L_{x}}^{+}\approx 400$). The visualized instantaneous field is the same as in figure 10.

Figure 21

Figure 20.Distributions of mean-square pressure fluctuations ($\overline{pp}/\unicode[STIX]{x1D70F}_{w}^{2}$, $\overline{p_{s}p_{s}}/\unicode[STIX]{x1D70F}_{w}^{2}$, $\overline{p_{r}p_{r}}/\unicode[STIX]{x1D70F}_{w}^{2}$) for $h^{+}=1020$: (a) LSD ($L_{x}=12.8h$); (b) MSU (${L_{x}}^{+}\approx 400$). In each figure, the relation of (3.8) is plotted as the two-dotted line.

Figure 22

Figure 21.Streamwise and spanwise power spectra of the total, slow and rapid pressure normalized by $\unicode[STIX]{x1D70F}_{w}^{2}$ and $h$, $\unicode[STIX]{x1D719}_{pp}(k_{x})/\unicode[STIX]{x1D70F}_{w}^{2}h$ and $\unicode[STIX]{x1D719}_{pp}(k_{z})/\unicode[STIX]{x1D70F}_{w}^{2}h$, for $h^{+}=1020$: (a) $\unicode[STIX]{x1D719}_{pp}(k_{x})/\unicode[STIX]{x1D70F}_{w}^{2}h$ at $y/h=0.5$; (b) $\unicode[STIX]{x1D719}_{pp}(k_{z})/\unicode[STIX]{x1D70F}_{w}^{2}h$ at $y/h=0.5$; (c) $\unicode[STIX]{x1D719}_{pp}(k_{x})/\unicode[STIX]{x1D70F}_{w}^{2}h$ for the slow pressure; (d) $\unicode[STIX]{x1D719}_{pp}(k_{x})/\unicode[STIX]{x1D70F}_{w}^{2}h$ for the rapid pressure. In (c,d), the $y$ stations of the spectra plotted are $y/h=0$, 0.2, 0.5, 0.7 and 1.

Figure 23

Figure 22.Distributions of the pressure strain terms for the total, slow and rapid pressure normalized by wall units for $h^{+}=1020$: (a) $\unicode[STIX]{x1D719}_{11}^{t}$, $\unicode[STIX]{x1D719}_{11}^{s}$, $\unicode[STIX]{x1D719}_{11}^{r}$; (b) $\unicode[STIX]{x1D719}_{22}^{t}$, $\unicode[STIX]{x1D719}_{22}^{s}$, $\unicode[STIX]{x1D719}_{22}^{r}$; (c) $\unicode[STIX]{x1D719}_{33}^{t}$, $\unicode[STIX]{x1D719}_{33}^{s}$, $\unicode[STIX]{x1D719}_{33}^{r}$.