Spectral energy cascade and decay in nonlinear acoustic waves

Abstract

We present a numerical and theoretical investigation of nonlinear spectral energy cascade of decaying finite-amplitude planar acoustic waves in a single-component ideal gas at standard temperature and pressure (STP). We analyze various one-dimensional canonical flow configurations: a propagating traveling wave (TW), a standing wave (SW), and randomly initialized Acoustic Wave Turbulence (AWT). Due to nonlinear wave propagation, energy at the large scales cascades down to smaller scales dominated by viscous dissipation, analogous to hydrodynamic turbulence. We use shock-resolved mesh-adaptive direct numerical simulations (DNS) of the fully compressible onedimensional Navier-Stokes equations to simulate the spectral energy cascade in nonlinear acoustic waves. The simulation parameter space for the TW, SW, and AWT cases spans three orders of magnitude in initial wave pressure amplitude and dynamic viscosity, thus covering a wide range of both spectral energy cascade and the viscous dissipation rates. The shock waves formed as a result of energy cascade are weak (M < 1.4), and hence we neglect thermodynamic non-equilibrium effects such as molecular vibrational relaxation in the current study. We also derive a new set of nonlinear acoustics equations truncated to second order and the corresponding perturbation energy corollary yielding the expression for a new perturbation energy norm E (2) . Its spatial average, <E (2) > satisfies the definition of a Lyapunov function, correctly capturing the inviscid (or lossless) broadening of spectral energy in the initial stages of evolution -analogous to the evolution of kinetic energy during the hydrodynamic break down of three-dimensional coherent vorticity -resulting in the formation of smaller scales. Upon saturation of the spectral energy cascade i.e. fully broadened energy spectrum, the onset of viscous losses causes a monotonic decay of <E (2) > in time. In this regime, the DNS results yield <E (2) > ∼ t −2 for TW and SW, and <E (2) > ∼ t −2/3 for AWT initialized with white noise. Using the perturbation energy corollary, we derive analytical expressions for the energy, energy flux, and dissipation rate in the wavenumber space. These yield the definitions of characteristic length scales such as the integral length scale (characteristic initial energy containing scale) and the Kolmogorov length scale η (shock thickness scale), analogous to K41 theory of hydrodynamic turbulence (A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 30 , 9 (1941)). Finally, we show that the fully developed energy spectrum of the nonlinear acoustic waves scales as E k k 2 −2/3 1/3 ∼ C f (kη), with C ≈ 0.075 constant for TW and SW but decaying in time for AWT. arXiv:1809.02202v1 [physics.flu-dyn]

Geophysical turbulence dominated by inertia–gravity waves

“…At smaller scales, the effects of rotation would be negligible and strongly nonlinear dynamics would dominate. Given that the shallow water equations are similar to the compressible gas dynamics equations (Whitham 2011), shocks in RSW models at small scales may be expected to generate a wave spectrum close to at these small scales, as is the scenario in compressible flows (Kuznetsov 2004; Falkovich & Kritsuk 2017; Gupta & Scalo 2018; Murray & Bustamante 2018).…”

Sub-filter-scale shear stress analysis in hypersonic turbulent Couette flow

“…The evaluation of the energy decay gives an estimation of the quality of the resolution in each element and allows for adapting the intensity of the subgrid dissipation locally. Recently, Sousa & Scalo (2022 b ) introduced the quasi-spectral viscosity (QSV) method, which is capable of unifying shock capturing and SFS modelling under a LES mathematical framework based on the concept that both hydrodynamic turbulence and shock formation are characterized by the energy cascade from large to small scales due to nonlinear interactions (Frisch 1995; Gupta & Scalo 2018), and they should be treated in a similar fashion. The QSV approach was also developed to be applicable to unstructured grids, via a block-spectral Legendre spectral decomposition (Sousa & Scalo 2022 a ).…”

High‐resolution numerical simulations of electrophoresis using the Fourier pseudo‐spectral method

“…Moreover, physical problems that are not naturally periodic, such as initial‐value problems, can also be solved on a periodic domain. For these reasons, pseudo‐spectral methods have found applications in various disciplines including meteorology [31], geophysics [32], fluid dynamics [33], and acoustics [34]. All electrophoresis techniques involve wave‐type phenomena and can be formulated as initial value problems [35–37].…”