No Arabic 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). We use shock-resolved mesh-adaptive direct numerical simulations (DNS) of the fully compressible one-dimensional Navier-Stokes equations to simulate the spectral energy cascade in nonlinear acoustic waves. 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.
We begin with the theoretical study of spectral energy cascade due to the propagation of high amplitude sound in the absence of thermal sources. To this end, a first-principles-based system of governing equations, correct up to second order in perturbation variables is derived. The exact energy corollary of such second-order system of equations is then formulated and used to elucidate the spectral energy dynamics of nonlinear acoustic waves. We then extend this analysis to thermoacoustically unstable waves -- i.e. amplified as a result of thermoacoustic instability. We drive such instability up until the generation of shock waves. We further study the nonlinear wave propagation in geometrically complex case of waves induced by the spark plasma between the electrodes. This case adds the geometrical complexity of a curved, three-dimensional shock, yielding vorticity production due to baroclinic torque. Finally, detonation waves are simulated by using a low-order approach, in a periodic setup subjected to high pressure inlet and exhaust of combustible gaseous mixture. An order adaptive fully compressible and unstructured Navier Stokes solver is currently under development to enable higher fidelity studies of both the spark plasma and detonation wave problem in the future.
Stratified turbulence shows scale- and direction-dependent anisotropy and the coexistence of weak turbulence of internal gravity waves and strong turbulence of eddies. Straightforward application of standard analyses developed in isotropic turbulence sometimes masks important aspects of the anisotropic turbulence. To capture detailed structures of the energy distribution in the wave-number space, it is indispensable to examine the energy distribution with non-integrated spectra by fixing the codimensional wave-number component or in the two-dimensional domain spanned by both the horizontal and vertical wave numbers. Indices which separate the range of the anisotropic weak-wave turbulence in the wave-number space are proposed based on the decomposed energies. In addition, the dominance of the waves in the range is also verified by the small frequency deviation from the linear dispersion relation. In the wave-dominant range, the linear wave periods given by the linear dispersion relation are smaller than approximately one third of the eddy-turnover time. The linear wave periods reflect the anisotropy of the system, while the isotropic Brunt-Vaisala period is used to evaluate the Ozmidov wave number, which is necessarily isotropic. It is found that the time scales in consideration of the anisotropy of the flow field must be appropriately selected to obtain the critical wave number separating the weak-wave turbulence.
The 4/5-law of turbulence, which characterizes the energy cascade from large to small-sized eddies at high Reynolds numbers in classical fluids, is verified experimentally in a superfluid 4He wind tunnel, operated down to 1.56 K and up to R_lambda ~ 1640. The result is corroborated by high-resolution simulations of Landau-Tiszas two-fluid model down to 1.15 K, corresponding to a residual normal fluid concentration below 3 % but with a lower Reynolds number of order R_lambda ~ 100. Although the Karman-Howarth equation (including a viscous term) is not valid emph{a priori} in a superfluid, it is found that it provides an empirical description of the deviation from the ideal 4/5-law at small scales and allows us to identify an effective viscosity for the superfluid, whose value matches the kinematic viscosity of the normal fluid regardless of its concentration.
The conventional approach to the turbulent energy cascade, based on Richardson-Kolmogorov phenomenology, ignores the topology of emerging vortices, which is related to the helicity of the turbulent flow. It is generally believed that helicity can play a significant role in turbulent systems, e.g., supporting the generation of large-scale magnetic fields, but its impact on the energy cascade to small scales has never been observed. We suggest for the first time a generalized phenomenology for isotropic turbulence with an arbitrary spectral distribution of the helicity. We discuss various scenarios of direct turbulent cascades with new helicity effect, which can be interpreted as a hindering of the spectral energy transfer. Therefore the energy is accumulated and redistributed so that the efficiency of non-linear interactions will be sufficient to provide a constant energy flux. We confirm our phenomenology by high Reynolds number numerical simulations based on a shell model of helical turbulence. The energy in our model is injected at a certain large scale only, whereas the source of helicity is distributed over all scales. In particular, we found that the helical bottleneck effect can appear in the inertial interval of the energy spectrum.
Nonlinear dynamics of surface gravity waves trapped by an opposing jet current is studied analytically and numerically. For wave fields narrowband in frequency but not necessarily with narrow angular distributions the developed asymptotic weakly nonlinear theory based on the modal approach of (V. Shrira, A. Slunyaev, J. Fluid. Mech, 738, 65, 2014) leads to the one-dimensional modified nonlinear Schr{o}dinger equation of self-focusing type for a single mode. Its solutions such as envelope solitons and breathers are considered to be prototypes of rogue waves; these solutions, in contrast to waves in the absence of currents, are robust with respect to transverse perturbations, which suggests potentially higher probability of rogue waves. Robustness of the long-lived analytical solutions in form of the modulated trapped waves and solitary wave groups is verified by direct numerical simulations of potential Euler equations.