Stochastic modal velocity field in rough-wall turbulence

• Published in: Journal of fluid mechanics 2024-11, Vol.999, Article A56
• Main Authors: Ehsani, Roozbeh, Heisel, Michael, Puccioni, Matteo, Hong, Jiarong, Iungo, Valerio, Voller, Vaughan, Guala, Michele
• Format: Article
• Language: English
• Subjects: Atmospheric boundary layer; Autocorrelation; Autocorrelation functions; Boundary layers; Datasets; Flow simulation; Fluid dynamics; Fluid flow; Height; JFM Papers; Momentum; Pressure dependence; Pressure gradients; Reynolds number; Scaling; Shear layers; Spatial discrimination; Spatial resolution; Statistical analysis; Stochasticity; Thickness; Turbulence; Turbulent boundary layer; Turbulent flow; Velocity; Velocity distribution; Velocity profiles; Vortices; Wind tunnels
• ISSN: 0022-1120; 1469-7645
• DOI: 10.1017/jfm.2024.933

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the $-u^{\prime }w^{\prime }$ co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.