An official website of the United States government
Abstract Several new methods are proposed that can diagnose the interscale transfer (or spectral flux) of kinetic energy (KE) and other properties in oceanic and broader geophysical systems, using integrals of advective structure functions and Bessel functions (herein “Bessel methods”). The utility of the Bessel methods is evaluated using simulations of anisotropic flow within two-dimensional (2D), surface quasigeostrophic (SQG), and two-layer QG systems. The Bessel methods diagnose various spectral fluxes within all of these systems, even under strong anisotropy and complex dynamics (e.g., multiple cascaded variables, coincident and opposing spectral fluxes, and nonstationary systems). In 2D turbulence, the Bessel methods capture the inverse KE cascade at large scales and the downscale enstrophy cascade (and associated downscale energy flux) at small scales. In SQG turbulence, the Bessel methods capture the downscale buoyancy variance cascade and the coincident upscale wavenumber-dependent KE flux. In QG turbulence, the Bessel methods capture the upscale kinetic energy flux. It is shown that these Bessel methods can be applied to data with limited extent or resolution, provided the scales of interest are captured by the range of separation distances. The Bessel methods are shown to have several advantages over other flux-estimation methods, including the ability to diagnose downscale energy cascades and to identify sharp transition scales. Analogous Bessel methods are also discussed for third-order structure functions, along with some caveats due to boundary terms. Significance StatementBig ocean eddies play an important role in Earth’s energy cycle by moving energy to both larger and smaller scales, but it is difficult to measure these “eddy energy fluxes” from oceanic observations. We develop a new method to estimate eddy energy fluxes that utilizes spatial differences between pairs of points and can be applied to various ocean data. This new method accurately diagnoses key eddy energy flux properties, as we demonstrate using idealized numerical simulations of various large-scale ocean systems.
more »
« less
Advective structure functions in anisotropic two-dimensional turbulence
In inertial-range turbulence, structure functions can diagnose transfer or dissipation rates of energy and enstrophy, which are difficult to calculate directly in flows with complex geometry or sparse sampling. However, existing relations between third-order structure functions and these rates only apply under isotropic conditions. We propose new relations to diagnose energy and enstrophy dissipation rates in anisotropic two-dimensional (2-D) turbulence. These relations use second-order advective structure functions that depend on spatial increments of vorticity, velocity, and their advection. Numerical simulations of forced-dissipative anisotropic 2-D turbulence are used to compare new and existing relations against model-diagnosed dissipation rates of energy and enstrophy. These simulations permit a dual cascade where forcing is applied at an intermediate scale, energy is dissipated at large scales, and enstrophy is dissipated at small scales. New relations to estimate energy and enstrophy dissipation rates show improvement over existing methods through increased accuracy, insensitivity to sampling direction, and lower temporal and spatial variability. These benefits of advective structure functions are present under weakly anisotropic conditions, and increase with the flow anisotropy as third-order structure functions become increasingly inappropriate. Several of the structure functions also show promise for diagnosing the forcing scale of 2-D turbulence. Velocity-based advective structure functions show particular promise as they can diagnose both enstrophy and energy cascade rates, and are robust to changes in the effective resolution of local derivatives. Some existing and future datasets that are amenable to advective structure function analysis are discussed.
more »
« less
- PAR ID:
- 10248590
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 916
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We derive and test a new heuristic theory for third-order structure functions that resolves the forcing scale in the scenario of simultaneous spectral energy transfer to both small and large scales, which can occur naturally, for example, in rotating stratified turbulence or magnetohydrodynamical (MHD) turbulence. The theory has three parameters – namely the upscale/downscale energy transfer rates and the forcing scale – and it includes the classic inertial-range theories as local limits. When applied to measured data, our global-in-scale theory can deduce the energy transfer rates using the full range of data, therefore it has broader applications compared with the local theories, especially in situations where the data is imperfect. In addition, because of the resolution of forcing scales, the new theory can detect the scales of energy input, which was impossible before. We test our new theory with a two-dimensional simulation of MHD turbulence.more » « less
-
Abstract Exact laws for evaluating cascade rates, tracing back to the Kolmogorov “4/5” law, have been extended to many systems of interest including magnetohydrodynamics (MHD), and compressible flows of the magnetofluid and ordinary fluid types. It is understood that implementations may be limited by the quantity of available data and by the lack of turbulence symmetry. Assessment of the accuracy and feasibility of such third-order (or Yaglom) relations is most effectively accomplished by examining the von Kármán–Howarth equation in increment form, a framework from which the third-order laws are derived as asymptotic approximations. Using this approach, we examine the context of third-order laws for incompressible MHD in some detail. The simplest versions rely on the assumption of isotropy and the presence of a well-defined inertial range, while related procedures generalize the same idea to arbitrary rotational symmetries. Conditions for obtaining correct and accurate values of the dissipation rate from these laws based on several sampling and fitting strategies are investigated using results from simulations. The questions we address are of particular relevance to sampling of solar wind turbulence by one or more spacecraft.more » « less
-
Abstract The cross‐shore transformation of breaking‐wave roller momentum and energy on observed barred surfzone bathymetry is investigated with a two‐phase Reynolds Averaged Navier Stokes model driven with measured incident waves. Modeled wave spectra, wave heights, and wave‐driven increases in the mean water level (setup) agree well with field observations along transects extending from 5‐m water depth to the shoreline. Consistent with prior results the roller forcing contributes 50%–60% to the setup, whereas the advective terms contribute ∼20%, with the contribution of bottom stress largest (up to 20%) for shallow sandbar crest depths. The model simulations suggest that an energy‐flux balance between wave dissipation, roller energy, and roller dissipation is accurate. However, as little as 70% of the modeled wave energy ultimately dissipated by breaking was first transferred from the wave to the roller. Furthermore, of the energy transferred to the roller, 15%–25% is dissipated by turbulence in the water column below the roller, with the majority of energy dissipated in the aerated region or near the roller‐surface interface. The contributions of turbulence to the momentum balance are sensitive to the parameterized turbulent anisotropy, which observations suggest increases with increasing turbulence intensity. Here, modeled turbulent kinetic energy dissipation decreases with increasing depth of the sandbar crest, possibly reflecting a change from plunging (on the steeper offshore slope of the bar) to spilling breakers (over the flatter bar crest and trough). Thus, using a variable roller front slope in the roller‐wave energy flux balance may account for these variations in breaker type.more » « less
-
We present an idealized study of rotating stratified wave turbulence in a two-dimensional vertical slice model of the Boussinesq equations, focusing on the peculiar case of equal Coriolis and buoyancy frequencies. In this case the fully nonlinear fluid dynamics can be shown to be isotropic in the vertical plane, which allows the classical methods of isotropic turbulence to be applied. Contrary to ordinary two-dimensional turbulence, here a robust downscale flux of total energy is observed in numerical simulations that span the full parameter regime between Ozmidov and forcing scales. Notably, this robust downscale flux of the total energy does not hold separately for its various kinetic and potential components, which can exhibit both upscale and downscale fluxes, depending on the parameter regime. Using a suitable extension of the classical Kármán–Howarth–Monin equation, exact expressions that link third-order structure functions and the spectral energy flux are derived and tested against numerical results. These expressions make obvious that even though the total energy is robustly transferred downscale, the third-order structure functions are sign indefinite, which illustrates that the sign and the form of measured third-order structure functions are both crucially important in determining the direction of the spectral energy transfer.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.247
