Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains ( $${>}300\delta _i$$ , with $$\delta _i$$ the inflow boundary-layer thickness) at three different Mach numbers, $2.5$ , $4.9$ and $10.9$ , with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number $$Re_\tau \approx 1200$$ ) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region.
more »
« less
Prograde vortices, internal shear layers and the Taylor microscale in high-Reynolds-number turbulent boundary layers
The statistical properties of prograde spanwise vortex cores and internal shear layers (ISLs) are evaluated for a series of high-Reynolds-number turbulent boundary layers. The considered flows span a wide range of both Reynolds number and surface roughness. In each case, the largest spanwise vortex cores in the outer layer of the boundary layer have size comparable to the Taylor microscale $$\lambda _T$$ , and the azimuthal velocity of these large vortex cores is governed by the friction velocity $${u_\tau }$$ . The same scaling parameters describe the average thickness and velocity difference across the ISLs. The results demonstrate the importance of the local large-eddy turnover time in determining the strain rate confining the size of the vortex cores and shear layers. The relevance of the turnover time, and more generally the Taylor microscale, can be explained by a stretching mechanism involving the mutual interaction of coherent velocity structures such as uniform momentum zones with the evolving shear layers separating the structures.
more »
« less
- Award ID(s):
- 2031312
- PAR ID:
- 10273225
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 920
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Since Taylor’s seminal paper, the existence of large-scale quasi-axisymmetric structures has been a matter of interest when studying Taylor–Couette flow. In this article, we probe their formation in the highly turbulent regime by conducting a series of numerical simulations at a fixed Reynolds number Re s = 3.6 × 10 4 while varying the Coriolis parameter to analyse the flow characteristics as the structures arise and dissipate. We show how the Coriolis force induces a one-way coupling between the radial and azimuthal velocity fields inside the boundary layer, but in the bulk, there is a two-way coupling that causes competing effects. We discuss how this complicates the analogy of narrow-gap Taylor–Couette to other convective flows. We then compare these statistics with a similar shear flow without no-slip boundary layers, showing how this double coupling causes very different effects. We finish by reflecting on the possible origins of turbulent Taylor rolls. This article is part of the theme issue ‘Taylor–Couette and related flows on the centennial of Taylor’s seminal Philosophical Transactions paper (part 1)’.more » « less
-
Direct Numerical Simulation (DNS) of compressible spatially-developing turbulent boundary layers (SDTBL) is performed at a Mach number of 2.5 and low/high Reynolds numbers over isothermal Zero-Pressure Gradient (ZPG) flat plates. Turbulent inflow information is generated via a dynamic rescaling-recycling approach (J. Fluid Mech., 670, pp. 581-605, 2011), which avoids the use of empirical correlations in the computation of inlet turbulent scales. The range of the low Reynolds number case is approximately 400-800, based on the momentum thickness, freestream velocity and wall viscosity. DNS at higher Reynolds numbers (~3,000, about four-fold larger) is also carried out with the purpose of analyzing the effect of Reynolds number on the transport phenomena in the supersonic regime. Additionally, low/high order flow statistics are compared with DNS of an incompressible isothermal ZPG boundary layer at similar low Reynolds numbers and the temperature regarded as a passive scalar. Peaks of turbulence intensities move closer to the wall as the Reynolds number increases in the supersonic flat plate. Furthermore, Reynolds shear stresses depict a much larger "plateau" (constant shear layer) at the highest Reynolds number considered in present study.more » « less
-
In this work, we introduce a scalable and efficient GPU-accelerated methodology for volumetric particle advection and finite-time Lyapunov exponent (FTLE) calculation, focusing on the analysis of Lagrangian coherent structures (LCS) in large-scale direct numerical simulation (DNS) datasets across incompressible, supersonic, and hypersonic flow regimes. LCS play a significant role in turbulent boundary layer analysis, and our proposed methodology offers valuable insights into their behavior in various flow conditions. Our novel owning-cell locator method enables efficient constant-time cell search, and the algorithm draws inspiration from classical search algorithms and modern multi-level approaches in numerical linear algebra. The proposed method is implemented for both multi-core CPUs and Nvidia GPUs, demonstrating strong scaling up to 32,768 CPU cores and up to 62 Nvidia V100 GPUs. By decoupling particle advection from other problems, we achieve modularity and extensibility, resulting in consistent parallel efficiency across different architectures. Our methodology was applied to calculate and visualize the FTLE on four turbulent boundary layers at different Reynolds and Mach numbers, revealing that coherent structures grow more isotropic proportional to the Mach number, and their inclination angle varies along the streamwise direction. We also observed increased anisotropy and FTLE organization at lower Reynolds numbers, with structures retaining coherency along both spanwise and streamwise directions. Additionally, we demonstrated the impact of lower temporal frequency sampling by upscaling with an efficient linear upsampler, preserving general trends with only 10% of the required storage. In summary, we present a particle search scheme for particle advection workloads in the context of visualizing LCS via FTLE that exhibits strong scaling performance and efficiency at scale. Our proposed algorithm is applicable across various domains, requiring efficient search algorithms in large, structured domains. While this article focuses on the methodology and its application to LCS, an in-depth study of the physics and compressibility effects in LCS candidates will be explored in a future publication.more » « less
-
In this work, we introduce a scalable and efficient GPU-accelerated methodology for volumetric particle advection and finite-time Lyapunov exponent (FTLE) calculation, focusing on the analysis of Lagrangian Coherent Structures (LCS) in large-scale Direct Numerical Simulation (DNS) datasets across incompressible, supersonic, and hypersonic flow regimes. LCS play a significant role in turbulent boundary layer analysis, and our proposed methodology offers valuable insights into their behavior in various flow conditions. Our novel owning-cell locator method enables efficient, constant-time cell search, and the algorithm draws inspiration from classical search algorithms and modern multi-level approaches in numerical linear algebra. The proposed method is implemented for both multi-core CPUs and Nvidia GPUs, demonstrating strong scaling up to 32,768 CPU cores and up to 62 Nvidia V100 GPUs. By decoupling particle advection from other problems, we achieve modularity and extensibility, resulting in consistent parallel efficiency across different architectures. Our methodology was applied to calculate and visualize the FTLE on four turbulent boundary layers at different Reynolds and Mach numbers, revealing that coherent structures grow more isotropic proportional to the Mach number, and their inclination angle varies along the streamwise direction. We also observed increased anisotropy and FTLE organization at lower Reynolds numbers, with structures retaining coherency along both spanwise and streamwise directions. Additionally, we demonstrated the impact of lower temporal frequency sampling by upscaling with an efficient linear upsampler, preserving general trends with only 10% of the required storage. In summary, we present a particle search scheme for particle advection workloads in the context of visualizing LCS via FTLE that exhibits strong scaling performance and efficiency at scale. Our proposed algorithm is applicable across various domains requiring efficient search algorithms in large structured domains. While this manuscript focuses on the methodology and its application to LCS, an in-depth study of the physics and compressibility effects in LCS candidates will be explored in a future publication.more » « less
An official website of the United States government

