An official website of the United States government
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
A statistically stationary minimal flow unit for self-similar Rayleigh–Taylor turbulence in the mode-coupling limit
We propose a computational framework for simulating the self-similar regime of turbulent Rayleigh–Taylor (RT) mixing layers in a statistically stationary manner. By leveraging the anticipated self-similar behaviour of RT mixing layers, a transformation of the vertical coordinate and velocities is applied to the Navier–Stokes equations (NSE), yielding modified equations that resemble the original NSE but include two sets of additional terms. Solving these equations, a statistically stationary RT (SRT) flow is achieved. Unlike temporally growing Rayleigh–Taylor (TRT) flow, SRT flow is independent of initial conditions and can be simulated over infinite simulation time without escalating resolution requirements, hence guaranteeing statistical convergence. Direct numerical simulations (DNS) are performed at an Atwood number of 0.5 and unity Schmidt number. By varying the ratio of the mixing layer height to the domain width, a minimal flow unit of aspect ratio 1.5 is found to approximate TRT turbulence in the self-similar mode-coupling regime. The SRT minimal flow unit has one-sixteenth the number of grid points required by the equivalent TRT simulation of the same Reynolds number and grid resolution. The resultant flow corresponds to a theoretical limit where self-similarity is observed in all fields and across the entire spatial domain – a late-time state that existing experiments and DNS of TRT flow have difficulties attaining. Simulations of the SRT minimal flow unit span TRT-equivalent Reynolds numbers (based on mixing layer height) ranging from 500 to 10 800. The SRT results are validated against TRT data from this study as well as from Cabot & Cook (Nat. Phys., vol. 2, 2006, pp. 562–568).
more »
« less
- Award ID(s):
- 2422513
- PAR ID:
- 10631937
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 1002
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A numerical investigation of an asymptotically reduced model for quasigeostrophic Rayleigh-Bénard convection is conducted in which the depth-averaged flows are numerically suppressed by modifying the governing equations. At the largest accessible values of the Rayleigh number Ra, the Reynolds number and Nusselt number show evidence of approaching the diffusion-free scalings of Re ∼ RaE/Pr and Nu ∼ Pr−1/2Ra3/2E2, respectively, where E is the Ekman number and Pr is the Prandtl number. For large Ra, the presence of depth-invariant flows, such as large-scale vortices, yield heat and momentum transport scalings that exceed those of the diffusion-free scaling laws. The Taylor microscale does not vary significantly with increasing Ra, whereas the integral length scale grows weakly. The computed length scales remain O(1) with respect to the linearly unstable critical wave number; we therefore conclude that these scales remain viscously controlled. We do not find a point-wise Coriolis-inertia-Archimedean (CIA) force balance in the turbulent regime; interior dynamics are instead dominated by horizontal advection (inertia), vortex stretching (Coriolis) and the vertical pressure gradient. A secondary, subdominant balance between the Archimedean buoyancy force and the viscous force occurs in the interior and the ratio of the root mean square (rms) of these two forces is found to approach unity with increasing Ra. This secondary balance is attributed to the turbulent fluid interior acting as the dominant control on the heat transport. These findings indicate that a pointwise CIA balance does not occur in the high Rayleigh number regime of quasigeostrophic convection in the plane layer geometry. Instead, simulations are characterized by what may be termed a nonlocal CIA balance in which the buoyancy force is dominant within the thermal boundary layers and is spatially separated from the interior Coriolis and inertial forces.more » « less
-
Convection is a ubiquitous process driving geophysical/astrophysical fluid flows, which are typically strongly constrained by planetary rotation on large scales. A celebrated model of such flows, rapidly rotating Rayleigh-Bénard convection, has been extensively studied in direct numerical simulations (DNS) and laboratory experiments, but the parameter values attainable by state-of-the-art methods are limited to moderately rapid rotation (Ekman numbers Ek≳10−8), while realistic geophysical/astrophysical Ek are significantly smaller. Asymptotically reduced equations of motion, the nonhydrostatic quasi-geostrophic equations (NHQGE), describing the flow evolution in the limit Ek→0, do not apply at finite rotation rates. The geophysical/astrophysical regime of small but finite Ek therefore remains currently inaccessible. Here, we introduce a new, numerically advantageous formulation of the Navier-Stokes-Boussinesq equations informed by the scalings valid for Ek→0, the \textit{Rescaled Rapidly Rotating incompressible Navier-Stokes Equations} (RRRiNSE). We solve the RRRiNSE using a spectral quasi-inverse method resulting in a sparse, fast algorithm to perform efficient DNS in this previously unattainable parameter regime. We validate our results against the literature across a range of Ek and demonstrate that the algorithmic approaches taken remain accurate and numerically stable at Ek as low as 10−15. Like the NHQGE, the RRRiNSE derive their efficiency from adequate conditioning, eliminating spurious growing modes that otherwise induce numerical instabilities at small Ek. We show that the time derivative of the mean temperature is inconsequential for accurately determining the Nusselt number in the stationary state, significantly reducing the required simulation time, and demonstrate that full DNS using RRRiNSE agree with the NHQGE at very small Ek.more » « less
-
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
-
We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work $$\varLambda$$ , due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work $$\varPi$$ , due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that $$\varLambda$$ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where $$\varPi$$ takes over. In three dimensions, $$\varPi$$ continues a persistent cascade to smaller scales, whereas in two dimensions $$\varPi$$ rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger $$\alpha$$ values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work $$\varLambda$$ .more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2024.1104
