Horizontal density layers are commonly observed in stratified turbulence. Recent work (e.g. Taylor & Zhou, J. Fluid Mech. , vol. 823, 2017, R5) has reinvigorated interest in the Phillips instability (PI), by which density layers form via negative diffusion if the turbulent buoyancy flux weakens as stratification increases. Theoretical understanding of PI is incomplete, in part because it remains unclear whether and by what mechanism the flux-gradient relationship for a given example of turbulence has the required negative-diffusion property. Furthermore, the difficulty of analysing the flux-gradient relation in evolving turbulence obscures the operating mechanism when layering is observed. These considerations motivate the search for an example of PI that can be analysed clearly. Here PI is shown to occur in two-dimensional Boussinesq sheared stratified turbulence maintained by stochastic excitation. PI is analysed using the second-order S3T closure of statistical state dynamics, in which the dynamics is written directly for statistical variables of the turbulence. The predictions of S3T are verified using nonlinear simulations. This analysis provides theoretical underpinning of PI based on the fundamental equations of motion that complements previous analyses based on phenomenological models of turbulence.
more »
« less
Statistical state dynamics of vertically sheared horizontal flows in two-dimensional stratified turbulence
Simulations of strongly stratified turbulence often exhibit coherent large-scale structures called vertically sheared horizontal flows (VSHFs). VSHFs emerge in both two-dimensional (2D) and three-dimensional (3D) stratified turbulence with similar vertical structure. The mechanism responsible for VSHF formation is not fully understood. In this work, the formation and equilibration of VSHFs in a 2D Boussinesq model of stratified turbulence is studied using statistical state dynamics (SSD). In SSD, equations of motion are expressed directly in the statistical variables of the turbulent state. Restriction to 2D turbulence facilitates application of an analytically and computationally attractive implementation of SSD referred to as S3T, in which the SSD is expressed by coupling the equation for the horizontal mean structure with the equation for the ensemble mean perturbation covariance. This second-order SSD produces accurate statistics, through second order, when compared with fully nonlinear simulations. In particular, S3T captures the spontaneous emergence of the VSHF and associated density layers seen in simulations of turbulence maintained by homogeneous large-scale stochastic excitation. An advantage of the S3T system is that the VSHF formation mechanism, which is wave–mean flow interaction between the emergent VSHF and the stochastically excited large-scale gravity waves, is analytically understood in the S3T system. Comparison with fully nonlinear simulations verifies that S3T solutions accurately predict the scale selection, dependence on stochastic excitation strength, and nonlinear equilibrium structure of the VSHF. These results constitute a theory for VSHF formation applicable to interpreting simulations and observations of geophysical examples of turbulent jets such as the ocean’s equatorial deep jets.
more »
« less
- Award ID(s):
- 1640989
- PAR ID:
- 10088442
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 854
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 544 to 590
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Geophysical turbulence has a wide range of spatiotemporal scales that requires a multiscale prediction model for efficient and fast simulations. Stochastic parameterization is a class of multiscale methods that approximates the large-scale behaviors of the turbulent system without relying on scale separation. In the stochastic parameterization of unresolved subgrid-scale dynamics, there are several modeling parameters to be determined by tuning or fitting to data. We propose a strategy to estimate the modeling parameters in the stochastic parameterization of geostrophic turbulent systems. The main idea of the proposed approach is to generate data in a spatiotemporally local domain and use physical/statistical information to estimate the modeling parameters. In particular, we focus on the estimation of modeling parameters in the stochastic superparameterization, a variant of the stochastic parameterization framework, for an idealized model of synoptic-scale turbulence in the atmosphere and oceans. The test regimes considered in this study include strong and moderate turbulence with complicated patterns of waves, jets, and vortices.more » « less
-
Understanding and predicting turbulent flow phenomena remain a challenge for both theory and applications. The nonlinear and nonlocal character of small-scale turbulence can be comprehensively described in terms of the velocity gradients, which determine fundamental quantities like dissipation, enstrophy, and the small-scale topology of turbulence. The dynamical equation for the velocity gradient succinctly encapsulates the nonlinear physics of turbulence; it offers an intuitive description of a host of turbulence phenomena and enables establishing connections between turbulent dynamics, statistics, and flow structure. The consideration of filtered velocity gradients enriches this view to express the multiscale aspects of nonlinearity and flow structure in a formulation directly applicable to large-eddy simulations. Driven by theoretical advances together with growing computational and experimental capabilities, recent activities in this area have elucidated key aspects of turbulence physics and advanced modeling capabilities.more » « less
-
We investigate how turbulence in liquid affects the rising speed of gas bubbles within the inertial range. Experimentally, we employ stereoscopic tracking of bubbles rising through water turbulence created by the convergence of turbulent jets and characterized with particle image velocimetry performed throughout the measurement volume. We use the spatially varying, time-averaged mean water velocity field to consider the physically relevant bubble slip velocity relative to the mean flow. Over a range of bubble sizes within the inertial range, we find that the bubble mean rise velocity $$\left \langle v_z \right \rangle$$ decreases with the intensity of the turbulence as characterized by its root-mean-square fluctuation velocity, $u'$ . Non-dimensionalized by the quiescent rise velocity $$v_{q}$$ , the average rise speed follows $$\left \langle v_z \right \rangle /v_{q}\propto 1/{\textit {Fr}}$$ at high $${\textit {Fr}}$$ , where $${\textit {Fr}}=u'/\sqrt {dg}$$ is a Froude number comparing the intensity of the turbulence to the bubble buoyancy, with $$d$$ the bubble diameter and $$g$$ the acceleration due to gravity. We complement these results by performing numerical integration of the Maxey–Riley equation for a point bubble experiencing nonlinear drag in three-dimensional, homogeneous and isotropic turbulence. These simulations reproduce the slowdown observed experimentally, and show that the mean magnitude of the slip velocity is proportional to the large-scale fluctuations of the flow velocity. Combining the numerical estimate of the slip velocity magnitude with a simple theoretical model, we show that the scaling $$\left \langle v_z \right \rangle /v_{q}\propto 1/{\textit {Fr}}$$ originates from a combination of the nonlinear drag and the nearly isotropic behaviour of the slip velocity at large $${\textit {Fr}}$$ that drastically reduces the mean rise speed.more » « less
-
Abstract Dispersion processes in the ocean surface boundary layer (OSBL) determine marine material distributions such as those of plankton and pollutants. Sheared velocities drive shear dispersion, which is traditionally assumed to be due to mean horizontal currents that decrease from the surface. However, OSBL turbulence supports along-wind jets; located in near-surface convergence and downwelling regions, such turbulent jets contain strong local shear. Through wind-driven idealized and large-eddy simulation (LES) models of the OSBL, this study examines the role of turbulent along-wind jets in dispersing material. In the idealized model, turbulent jets are generated by prescribed cellular flow with surface convergence and associated downwelling regions. Numeric and analytic model solutions reveal that horizontal jets substantially contribute to along-wind dispersion for sufficiently strong cellular flows and exceed contributions due to vertical mean shear for buoyant surface-trapped material. However, surface convergence regions also accumulate surface-trapped material, reducing shear dispersion by jets. Turbulence resolving LES results of a coastal depth-limited ocean agree qualitatively with the idealized model and reveal long-lived coherent jet structures that are necessary for effective jet dispersion. These coastal results indicate substantial jet contributions to along-wind dispersion. However, jet dispersion is likely less effective in the open ocean because jets are shorter lived, less organized, and distorted due to spiraling Ekman currents.more » « less
An official website of the United States government

