skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence

Abstract

Direct numerical simulations of low Mach number compressible three-dimensional magnetohydrodynamic (CMHD3D) turbulence in the presence of a strong mean magnetic field are compared with simulations of reduced magnetohydrodynamics (RMHD). Periodic boundary conditions in the three spatial coordinates are considered. Different sets of initial conditions are chosen to explore the applicability of RMHD and to study how close the solution remains to the full compressible MHD solution as both freely evolve in time. In a first set, the initial state is prepared to satisfy the conditions assumed in the derivation of RMHD, namely, a strong mean magnetic field and plane-polarized fluctuations, varying weakly along the mean magnetic field. In those circumstances, simulations show that RMHD and CMHD3D evolve almost indistinguishably from one another. When some of the conditions are relaxed the agreement worsens but RMHD remains fairly close to CMHD3D, especially when the mean magnetic field is large enough. Moreover, the well-known spectral anisotropy effect promotes the dynamical attainment of the conditions for RMHD applicability. Global quantities (mean energies, mean-square current, and vorticity) and energy spectra from the two solutions are compared and point-to-point separation estimations are computed. The specific results shown here give support to the use of RMHD asmore » a valid approximation of compressible MHD with a mean magnetic field under certain but quite practical conditions.« less

Authors:
; ;  [1];  [2]
  1. Bartol Research Institute, University of Delaware, Newark, Delaware 19716 (United States)
  2. (New Zealand)
Publication Date:
OSTI Identifier:
20782333
Resource Type:
Journal Article
Resource Relation:
Journal Name: Physics of Plasmas; Journal Volume: 12; Journal Issue: 11; Other Information: DOI: 10.1063/1.2128573; (c) 2005 American Institute of Physics; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; ANISOTROPY; APPROXIMATIONS; BOUNDARY CONDITIONS; CHARGED-PARTICLE TRANSPORT; COMPARATIVE EVALUATIONS; ELECTRIC CURRENTS; ENERGY SPECTRA; FLUCTUATIONS; MACH NUMBER; MAGNETIC FIELDS; MAGNETOHYDRODYNAMICS; MATHEMATICAL SOLUTIONS; PERIODICITY; PLASMA; PLASMA SIMULATION; THREE-DIMENSIONAL CALCULATIONS; TURBULENCE; VORTICES

Citation Formats

Dmitruk, Pablo, Matthaeus, William H., Oughton, Sean, and Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton. Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence. United States: N. p., 2005. Web. doi:10.1063/1.2128573.
Dmitruk, Pablo, Matthaeus, William H., Oughton, Sean, & Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton. Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence. United States. doi:10.1063/1.2128573.
Dmitruk, Pablo, Matthaeus, William H., Oughton, Sean, and Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton. Tue . "Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence". United States. doi:10.1063/1.2128573.
@article{osti_20782333,
title = {Direct comparisons of compressible magnetohydrodynamics and reduced magnetohydrodynamics turbulence},
author = {Dmitruk, Pablo and Matthaeus, William H. and Oughton, Sean and Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton},
abstractNote = {Direct numerical simulations of low Mach number compressible three-dimensional magnetohydrodynamic (CMHD3D) turbulence in the presence of a strong mean magnetic field are compared with simulations of reduced magnetohydrodynamics (RMHD). Periodic boundary conditions in the three spatial coordinates are considered. Different sets of initial conditions are chosen to explore the applicability of RMHD and to study how close the solution remains to the full compressible MHD solution as both freely evolve in time. In a first set, the initial state is prepared to satisfy the conditions assumed in the derivation of RMHD, namely, a strong mean magnetic field and plane-polarized fluctuations, varying weakly along the mean magnetic field. In those circumstances, simulations show that RMHD and CMHD3D evolve almost indistinguishably from one another. When some of the conditions are relaxed the agreement worsens but RMHD remains fairly close to CMHD3D, especially when the mean magnetic field is large enough. Moreover, the well-known spectral anisotropy effect promotes the dynamical attainment of the conditions for RMHD applicability. Global quantities (mean energies, mean-square current, and vorticity) and energy spectra from the two solutions are compared and point-to-point separation estimations are computed. The specific results shown here give support to the use of RMHD as a valid approximation of compressible MHD with a mean magnetic field under certain but quite practical conditions.},
doi = {10.1063/1.2128573},
journal = {Physics of Plasmas},
number = 11,
volume = 12,
place = {United States},
year = {Tue Nov 15 00:00:00 EST 2005},
month = {Tue Nov 15 00:00:00 EST 2005}
}
  • In this work we numerically test a model of Hall magnetohydrodynamics in the presence of a strong mean magnetic field: the reduced Hall magnetohydrodynamic model (RHMHD) derived by [Gomez et al., Phys. Plasmas 15, 102303 (2008)] with the addition of weak compressible effects. The main advantage of this model lies in the reduction of computational cost. Nevertheless, up until now the degree of agreement with the original Hall MHD system and the range of validity in a regime of turbulence were not established. In this work direct numerical simulations of three-dimensional Hall MHD turbulence in the presence of a strongmore » mean magnetic field are compared with simulations of the weak compressible RHMHD model. The results show that the degree of agreement is very high (when the different assumptions of RHMHD, such as spectral anisotropy, are satisfied). Nevertheless, when the initial conditions are isotropic but the mean magnetic field is maintained strong, the results differ at the beginning but asymptotically reach a good agreement at relatively short times. We also found evidence that the compressibility still plays a role in the dynamics of these systems, and the weak compressible RHMHD model is able to capture these effects. In conclusion the weak compressible RHMHD model is a valid approximation of the Hall MHD turbulence in the relevant physical context.« less
  • We numerically study two-dimensional quantum turbulence with a Gross-Pitaevskii model. With the energy initially accumulated at large scale, quantum turbulence with many quantized vortex points is generated. Due to the lack of enstrophy conservation in this model, direct energy cascade with a Kolmogorov-Obukhov energy spectrum E(k){proportional_to}k{sup -5/3} is observed, which is quite different from two-dimensional incompressible classical turbulence in the decaying case. A positive value for the energy flux guarantees a direct energy cascade in the inertial range (from large to small scales). After almost all the energy at the large scale cascades to the small scale, the compressible kineticmore » energy realizes the thermodynamic equilibrium state without quantized vortices.« less
  • The recently developed Gas Kinetic Method (GKM) for computing fluid flow is enhanced with advanced reconstruction (interpolation) schemes to enable direct simulations of highly compressible transition and turbulence fields. Variants of Weighted Essentially Non-Oscillatory (WENO) reconstruction schemes of different orders of accuracy are implemented and examined along with more elementary van Leer method. The competing schemes are evaluated for their accuracy, efficiency and numerical stability. The computed results are compared against the Rapid Distortion Theory for the case of compressible shear turbulence and 'pressure-released' Burgers solution. In the case of decaying isotropic turbulence, the efficacy of the reconstruction schemes ismore » evaluated by comparison against a 'gold standard' high-resolution simulation. The capabilities of the reconstruction schemes to capture linear, non-linear, pressure-released and viscous flow physics as well as solenoidal and dilatational features of the flow fields are established in isolation and combination. The most suitable WENO variant for integration with GKM is identified. Another important outcome of the study is the finding that temperature-interpolation is superior to energy-interpolation in avoiding negative temperatures arising due to the Gibbs phenomenon. Overall, this work advances the applicability of kinetic theory based GKM to a wider range of high Mach number flow physics.« less
  • We present direct numerical simulations of inhomogeneous reduced magnetohydrodynamic (RMHD) turbulence between the Sun and the Alfvén critical point. These are the first such simulations that take into account the solar-wind outflow velocity and the radial inhomogeneity of the background solar wind without approximating the nonlinear terms in the governing equations. RMHD turbulence is driven by outward-propagating Alfvén waves (z {sup +} fluctuations) launched from the Sun, which undergo partial non-WKB reflection to produce sunward-propagating Alfvén waves (z {sup –} fluctuations). We present 10 simulations with different values of the correlation time τ{sub c{sub sun}{sup +}} and perpendicular correlation lengthmore » L{sub ☉} of outward-propagating Alfvén waves at the coronal base. We find that between 15% and 33% of the z {sup +} energy launched into the corona dissipates between the coronal base and Alfvén critical point. Between 33% and 40% of this input energy goes into work on the solar-wind outflow, and between 22% and 36% escapes as z {sup +} fluctuations through the simulation boundary at r = r{sub A}. The z {sup ±} power spectra scale like k{sub perpendicular}{sup -α{sup ±}}, where k is the wavenumber in the plane perpendicular to B{sub 0}. In our simulation with the smallest value of τ{sub c{sub sun}{sup +}} (∼2 minutes) and largest value of L{sub ☉} (2 × 10{sup 4} km), we find that α{sup +} decreases approximately linearly with increasing ln (r), reaching a value of 1.3 at r = 11.1 R{sub ☉}. Our simulations with larger values of τ{sub c{sub sun}{sup +}} exhibit alignment between the contours of constant φ{sup +}, φ{sup –}, Ω{sub 0}{sup +}, and Ω{sub 0}{sup -}, where φ{sup ±} are the Elsässer potentials and Ω{sub 0}{sup ±} are the outer-scale parallel Elsässer vorticities.« less