skip to main content

This content will become publicly available on November 4, 2023

Title: Alternative formulation of weak magnetohydrodynamic turbulence theory

In a recent paper [P. H. Yoon and G. Choe, Phys. Plasmas 28, 082306 (2021)], the weak turbulence theory for incompressible magnetohydrodynamics is formulated by employing the method customarily applied in the context of kinetic weak plasma turbulence theory. Such an approach simplified certain mathematical procedures including achieving the closure relationship. The formulation in the above-cited paper starts from the equations of incompressible magnetohydrodynamic (MHD) theory expressed via Elsasser variables. The derivation of nonlinear wave kinetic equation therein is obtained via a truncated solution at the second-order of iteration following the standard practice. In the present paper, the weak MHD turbulence theory is alternatively formulated by employing the pristine form of incompressible MHD equation rather than that expressed in terms of Elsasser fields. The perturbative expansion of the nonlinear momentum equation is carried out up to the third-order iteration rather than imposing the truncation at the second order. It is found that while the resulting wave kinetic equation is identical to that obtained in the previous paper cited above, the third-order nonlinear correction plays an essential role for properly calculating derived quantities such as the total and residual energies.

 ;  ;  
Publication Date:
Journal Name:
Physics of Plasmas
Page Range or eLocation-ID:
Article No. 112303
American Institute of Physics
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A fundamental question in wave turbulence theory is to understand how the wave kinetic equation describes the long-time dynamics of its associated nonlinear dispersive equation. Formal derivations in the physics literature, dating back to the work of Peierls in 1928, suggest that such a kinetic description should hold (for well-prepared random data) at a large kinetic time scale $T_{\mathrm {kin}} \gg 1$ and in a limiting regime where the size L of the domain goes to infinity and the strength $\alpha $ of the nonlinearity goes to $0$ (weak nonlinearity). For the cubic nonlinear Schrödinger equation, $T_{\mathrm {kin}}=O\left (\alpha ^{-2}\right )$ and $\alpha $ is related to the conserved mass $\lambda $ of the solution via $\alpha =\lambda ^2 L^{-d}$ . In this paper, we study the rigorous justification of this monumental statement and show that the answer seems to depend on the particular scaling law in which the $(\alpha , L)$ limit is taken, in a spirit similar to how the Boltzmann–Grad scaling law is imposed in the derivation of Boltzmann’s equation. In particular, there appear to be two favourable scaling laws: when $\alpha $ approaches $0$ like $L^{-\varepsilon +}$ or like $L^{-1-\frac {\varepsilon }{2}+}$ (for arbitrary smallmore »$\varepsilon $ ), we exhibit the wave kinetic equation up to time scales $O(T_{\mathrm {kin}}L^{-\varepsilon })$ , by showing that the relevant Feynman-diagram expansions converge absolutely (as a sum over paired trees). For the other scaling laws, we justify the onset of the kinetic description at time scales $T_*\ll T_{\mathrm {kin}}$ and identify specific interactions that become very large for times beyond $T_*$ . In particular, the relevant tree expansion diverges absolutely there. In light of those interactions, extending the kinetic description beyond $T_*$ toward $T_{\mathrm {kin}}$ for such scaling laws seems to require new methods and ideas.« less
  2. Abstract Weakly turbulent processes that take place in plasmas are customarily formulated in terms of kinetic theory. However, owing to an inherent complexity associated with the problem, thus far the theory is fully developed largely for unmagnetized plasmas. In the present paper it is shown that a warm two fluid theory can successfully be employed in order to partially formulate the weak turbulence theory in spatially uniform plasma. Specifically, it is shown that the nonlinear wave-wave interaction, or decay processes, can be reproduced by the two-fluid formalism. The present finding shows that the same approach can in principle be extended to magnetized plasmas, which is a subject of future work.
  3. 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. Comparisonmore »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.« less
  4. Abstract Up to the third-order rogue wave solutions of the Sasa–Satsuma (SS) equation are derived based on the Hirota’s bilinear method and Kadomtsev–Petviashvili hierarchy reduction method. They are expressed explicitly by rational functions with both the numerator and denominator being the determinants of even order. Four types of intrinsic structures are recognized according to the number of zero-amplitude points. The first- and second-order rogue wave solutions agree with the solutions obtained so far by the Darboux transformation. In spite of the very complicated solution form compared with the ones of many other integrable equations, the third-order rogue waves exhibit two configurations: either a triangle or a distorted pentagon. Both the types and configurations of the third-order rogue waves are determined by different choices of free parameters. As the nonlinear Schrödinger equation is a limiting case of the SS equation, it is shown that the degeneration of the first-order rogue wave of the SS equation converges to the Peregrine soliton.
  5. We derive new kinetic and a porous medium equations from the nonlinear Schrödinger equation with random potentials. The kinetic equation has a very similar form compared to the four-wave turbulence kinetic equation in the wave turbulence theory. Moreover, we construct a class of self-similar solutions for the porous medium equation. These solutions spread with time, and this fact answers the “weak turbulence” question for the nonlinear Schrödinger equation with random potentials. We also derive Ohm’s law for the porous medium equation.