An official website of the United States government
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
Magneto-Stokes flow in a shallow free-surface annulus
In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.
more »
« less
- Award ID(s):
- 2143939
- PAR ID:
- 10581132
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 996
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Mixed-layer depth (MLD) exhibits significant variability, which is important for atmosphere–ocean exchanges of heat and atmospheric gases. The origins of the mesoscale MLD variability in the Southern Ocean are studied here in an idealised regional ocean–atmosphere model (ROAM). The main conclusion from the analysis of the upper-ocean buoyancy budget is that, while the atmospheric forcing and oceanic vertical mixing, on average, induce the mesoscale variability of MLD, the three-dimensional oceanic advection of buoyancy counteracts and partially balances these atmosphere-induced vertical processes. The relative importance of advection changes with both season and average MLD. From January to May, when the mixed layer is shallow, the atmospheric forcing and oceanic mixing are the most important processes, with the advection playing a secondary role. From June to December, when the mixed layer is deep, both atmospheric forcing and oceanic advection are equally important in driving the MLD variability. Importantly, buoyancy advection by mesoscale ocean current anomalies can lead to both local shoaling and deepening of the mixed layer. The role of the atmospheric forcing is then directly addressed by two sensitivity experiments in which the mesoscale variability is removed from the atmosphere–ocean heat and momentum fluxes. The findings confirm that mesoscale atmospheric forcing predominantly controls MLD variability in summer and that intrinsic oceanic variability and surface forcing are equally important in winter. As a result, MLD variance increases when mesoscale anomalies in atmospheric fluxes are removed in winter, and oceanic advection becomes a dominant player in the buoyancy budget. This study highlights the importance of oceanic advection and intrinsic ocean dynamics in driving mesoscale MLD variability and underscores the importance of MLD in modulating the effects of advection on upper-ocean dynamics.more » « less
-
We present a combined experimental and theoretical study of the primary and secondary instabilities in a Kolmogorov-like flow. The experiment uses electromagnetic forcing with an approximately sinusoidal spatial profile to drive a quasi-two-dimensional (Q2D) shear flow in a thin layer of electrolyte suspended on a thin lubricating layer of a dielectric fluid. Theoretical analysis is based on a two-dimensional (2D) model (Suri et al. , Phys. Fluids , vol. 26 (5), 2014, 053601), derived from first principles by depth-averaging the full three-dimensional Navier–Stokes equations. As the strength of the forcing is increased, the Q2D flow in the experiment undergoes a series of bifurcations, which is compared with results from direct numerical simulations of the 2D model. The effects of confinement and the forcing profile are studied by performing simulations that assume spatial periodicity and strictly sinusoidal forcing, as well as simulations with realistic no-slip boundary conditions and an experimentally validated forcing profile. We find that only the simulation subject to physical no-slip boundary conditions and a realistic forcing profile provides close, quantitative agreement with the experiment. Our analysis offers additional validation of the 2D model as well as a demonstration of the importance of properly modelling the forcing and boundary conditions.more » « less
-
We propose a classification scheme for nocturnal atmospheric boundary layers and apply it to investigate the spatio‐temporal structure of air temperature and wind speed in a shallow valley during the Shallow Cold Pool Experiment. This field campaign was the first to collect spatially continuous temperature and wind information at high resolution (1 s, 0.25 m) using the distributed temperature sensing technique across a 220 m long transect at three heights (0.5, 1.0, 2.0 m). The night‐time classification scheme was motivated by a surface energy balance and used a combination of static stability, wind regime and longwave radiative forcing as quantities to determine physically meaningful boundary‐layer regimes. Out of all potential combinations of these three quantities, 14 night‐time classes contained observations, of which we selected three for detailed analysis and comparison. The three classes represent a transition from mechanical to radiative forcing. The first night class represents conditions with strong dynamic forcing caused by locally induced lee turbulence dominating near‐surface temperatures across the shallow valley. The second night class was a concurrence of enhanced dynamic mixing due to significant winds at the valley shoulders and cold‐air pooling at the bottom of the shallow valley as a result of strong radiative cooling. The third night class was characteristic of weak winds eliminating the impact of mechanical mixing but emphasizing the formation and pooling of cold air at the valley bottom. The proposed night‐time classification scheme was found to sort the experimental data into physically meaningful regimes of surface flow and transport. It is suitable to stratify short‐ and long‐term experimental data for ensemble averaging and to identify case studies.more » « less
-
Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities and rearranging flows, even on scales far above the particles’ gyroscales where kinetic effects are traditionally considered most important. Here, we study the influence of pressure anisotropy on turbulent plasmas threaded by a mean magnetic field (Alfvénic turbulence). Extending previous results that were concerned with Braginskii MHD, we consider a wide range of regimes and parameters using a simplified fluid model based on drift kinetics with heat fluxes calculated using a Landau-fluid closure. We show that viscous (pressure-anisotropy) heating dissipates between a quarter (in collisionless regimes) and half (in collisional regimes) of the turbulent-cascade power injected at large scales; this does not depend strongly on either plasma beta or the ion-to-electron temperature ratio. This will in turn influence the plasma's thermodynamics by regulating energy partition between different dissipation channels (e.g. electron and ion heat). Due to the pressure anisotropy's rapid dynamic feedback onto the flows that create it – an effect we term ‘magneto-immutability’ – the viscous heating is confined to a narrow range of scales near the forcing scale, supporting a nearly conservative, MHD-like inertial-range cascade, via which the rest of the energy is transferred to small scales. Despite the simplified model, our results – including the viscous heating rate, distributions and turbulent spectra – compare favourably with recent hybrid-kinetic simulations. This is promising for the more general use of extended-fluid (or even MHD) approaches to model weakly collisional plasmas such as the intracluster medium, hot accretion flows and the solar wind.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2024.674
