On locally embedded two-scale solution for wall-bounded turbulent flows

Recent findings on wall-bounded turbulence have prompted a new impetus for modelling development to capture and resolve the Reynolds-number-dependent influence of outer flow on near-wall turbulence in terms of the ‘foot-printing’ of the large-scale coherent structures and the scale-interaction associated ‘modulation’. We develop a two-scale method to couple a locally embedded near-wall fine-mesh direct numerical simulation (DNS) block with a global coarser mesh domain. The influence of the large-scale structures on the local fine-mesh block is captured by a scale-dependent coarse–fine domain interface treatment. The coarse-mesh resolved disturbances are directly exchanged across the interface, while only the fine-mesh resolved fluctuations around the coarse-mesh resolved variables are subject to periodic conditions in the streamwise and spanwise directions. The global near-wall coarse-mesh region outside the local fine-mesh block is governed by the augmented flow governing equations with forcing source terms generated by upscaling the space–time-averaged fine-mesh solution. The validity and effectiveness of the method are examined for canonical incompressible channel flows at several Reynolds numbers. The mean statistics and energy spectra are in good agreement with the corresponding full DNS data. The results clearly illustrate the ‘foot-printing’ and ‘modulation’ in the local fine-mesh block. Noteworthy also is that neither spectral-gap nor scale-separation is assumed, and a smooth overlap between the global-domain and the local-domain energy spectra is observed. It is shown that the mesh-count scaling with Reynolds number is potentially reduced from $O(R{e^2})$ for the conventional fully wall-resolved large-eddy simulation (LES) to $O(Re)$ for the present locally embedded two-scale LES.

The logarithmic variance of streamwise velocity and conundrum in wall turbulence

The logarithmic dependence of streamwise turbulence intensity has been observed repeatedly in recent experimental and direct numerical simulation data. However, its spectral counterpart, a well-developed $k^{-1}$ spectrum ( $k$ is the spatial wavenumber in a wall-parallel direction), has not been convincingly observed from the same data. In the present study, we revisit the spectrum-based attached eddy model of Perry and co-workers, who proposed the emergence of a $k^{-1}$ spectrum in the inviscid limit, for small but finite $z/\delta$ and for finite Reynolds numbers ( $z$ is the wall-normal coordinate, and $\delta$ is the outer length scale). In the upper logarithmic layer (or inertial sublayer), a reexamination reveals that the intensity of the spectrum must vary with the wall-normal location at order of $z/\delta$ , consistent with the early observation argued with ‘incomplete similarity’. The streamwise turbulence intensity is subsequently calculated, demonstrating that the existence of a well-developed $k^{-1}$ spectrum is not a necessary condition for the approximate logarithmic wall-normal dependence of turbulence intensity – a more general condition is the existence of a premultiplied power-spectral intensity of $O(1)$ for $O(1/\delta ) < k < O(1/z)$ . Furthermore, it is shown that the Townsend–Perry constant must be weakly dependent on the Reynolds number. Finally, the analysis is semi-empirically extended to the lower logarithmic layer (or mesolayer), and a near-wall correction for the turbulence intensity is subsequently proposed. All the predictions of the proposed model and the related analyses/assumptions are validated with high-fidelity experimental data (Samie et al., J. Fluid Mech., vol. 851, 2018, pp. 391–415).

Direct numerical simulation of turbulent duct flow with inclined or V-shaped ribs mounted on one wall

In this research, highly disturbed turbulent flow of distinct three-dimensional characteristics in a square duct with inclined or V-shaped ribs mounted on one wall is investigated using direct numerical simulation. The turbulence field is highly sensitive to not only the rib geometry but also the boundary layers developed over the side and top walls. In a cross-stream plane secondary flows appear as large longitudinal vortices in both inclined and V-shaped rib cases due to the confinement of four sidewalls of the square duct. However, owing to the difference in the pattern of cross-stream secondary flow motions, the flow physics is significantly different in these two ribbed duct cases. It is observed that the mean flow structures in the cross-stream directions are asymmetrical in the inclined rib case but symmetrical in the V-shaped rib case, causing substantial differences in the momentum transfer across the spanwise direction. The impacts of rib geometry on near-wall turbulence structures are investigated using vortex identifiers, joint probability density functions between the streamwise and vertical velocity fluctuations, statistical moments of different orders, spatial two-point autocorrelations and velocity spectra. It is found that near the leeward and windward rib faces, the mean inclination angle of turbulence structures in the V-shaped rib case is greater than that of the inclined rib case, which subsequently enhances momentum transport between the ribbed bottom wall and the smooth top wall.

Asymptotics of streamwise Reynolds stress in wall turbulence

The scaling of different features of streamwise normal stress profiles $\langle uu\rangle ^+(y^+)$ in turbulent wall-bounded flows is the subject of a long-running debate. Particular points of contention are the scaling of the ‘inner’ and ‘outer’ peaks of $\langle uu\rangle ^+$ at $y^+\approxeq ~15$ and $y^+ ={O}(10^3)$ , respectively, their infinite Reynolds number limit, and the rate of logarithmic decay in the outer part of the flow. Inspired by the thought-provoking paper of Chen & Sreenivasan (J. Fluid Mech., vol. 908, 2021, p. R3), two terms of an inner asymptotic expansion of $\langle uu\rangle ^+$ in the small parameter $Re_{\tau }^{-1/4}$ are constructed from a set of direct numerical simulations (DNS) of channel flow. This inner expansion is for the first time matched through an overlap layer to an outer expansion, which not only fits the same set of channel DNS within 1.5 % of the peak stress, but also provides a good match of laboratory data in pipes and the near-wall part of boundary layers, up to the highest $Re_{\tau }$ values of $10^5$ . The salient features of the new composite expansion are first, an inner $\langle uu\rangle ^+$ peak, which saturates at 11.3 and decreases as $Re_{\tau }^{-1/4}$ . This inner peak is followed by a short ‘wall log law’ with a slope that becomes positive for $Re_{\tau }$ beyond ${O}(10^4)$ , leading up to an outer peak, followed by the logarithmic overlap layer with a negative slope going continuously to zero for $Re_{\tau }\to \infty$ .

Collapse of transitional wall turbulence captured using a rare events algorithm

This text presents one of the first successful applications of a rare events sampling method for the study of multistability in a turbulent flow without stochastic energy injection. The trajectories of collapse of turbulence in plane Couette flow, and their probability and rate of occurrence are systematically computed using adaptive multilevel splitting (AMS). The AMS computations are performed in a system of size $L_x\times L_z=24\times 18$ at Reynolds number $R=370$ with an acceleration by a factor ${O}(10)$ with respect to direct numerical simulations (DNS) and in a system of size $L_x\times L_z=36\times 27$ at Reynolds number $R=377$ with an acceleration by a factor ${O}(10^3)$ . The AMS results are validated by a comparison with DNS in the smaller system. Visualisations indicate that turbulence collapses because the self-sustaining process of turbulence fails locally. The streamwise vortices decay first in streamwise elongated holes, leaving streamwise invariant streamwise velocity tubes that experience viscous decay. These holes then extend in the spanwise direction. The examination of more than a thousand trajectories in the $(E_{k,x}=\int u_x^2/2\,\textrm {d}^3\boldsymbol {x},E_{k,y-z}=\int (u_y^2/2+u_z^2/2)\,\textrm {d}^3\boldsymbol {x})$ plane in the smaller system confirms the faster decay of streamwise vortices and shows concentration of trajectories. This hints at an instanton phenomenology in the large size limit. The computation of turning point states, beyond which laminarisation is certain, confirms the hole formation scenario and shows that it is more pronounced in larger systems. Finally, the examination of non-reactive trajectories indicates that both the vortices and the streaks reform concomitantly when the laminar holes close.

Evolution of turbulent kinetic energy during the entire sandstorm process

An adaptive segmented stationary method for non-stationary signal is proposed to reveal the turbulent kinetic energy evolution during the entire sandstorm process observed at the Qingtu Lake Observation Array. Sandstorm which is a common natural disaster is mechanically characterized by a particle-laden two-phase flow experiencing wall turbulence, with an extremely high Reynolds number and significant turbulent kinetic energy. Turbulence energy transfer is important to the understanding of sandstorm dynamics. This study indicates that large-/very-large-scale coherent structures originally exist in the rising stage of sandstorms with a streamwise kinetic energy of 75 % rather than gradually forming. In addition to carrying a substantial portion of energy, the very-large-scale-motions are active structures with strong nonlinear energy transfer. These structures gain energy from strong nonlinear interaction. As sandstorm evolves, these large structures are gradually broken by quadratic phase coupling, with the energy fraction reducing to 40 % in the declining stage. The nonlinear process in the steady and declining stages weakens and maintains a balanced budget of energy. The systematic bispectrum results provide a new perspective for further insight of sandstorms.