Piecewise parabolic method on a local stencil for magnetized supersonic turbulence simulation

DOI10.1016/J.JCP.2009.07.007zbMATH Open1391.76583arXiv0905.2960OpenAlexW2143638188MaRDI QIDQ733016FDOQ733016

Authors: Sergey D. Ustyugov, M. V. Popov, Michael L. Norman, Alexei G. Kritsuk 

Publication date: 15 October 2009

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local "fixes" for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.

Full work available at URL: https://arxiv.org/abs/0905.2960

Recommendations

• 10.1007/s11470-008-3011-1
• The piecewise parabolic method (PPM) for gas-dynamical simulations
• Pseudospectral simulation of compressible magnetohydrodynamic turbulence
• An unsplit Godunov method for ideal MHD via constrained transport
• A piecewise parabolic method for cosmological hydrodynamics

zbMATH Keywords

numerical methodsgas dynamicsmagnetohydrodynamicscompact stencilMHD turbulencePPM

Mathematics Subject Classification ID

Gas dynamics (general theory) (76N15) Magnetohydrodynamics and electrohydrodynamics (76W05) Turbulence (76F99) Supersonic flows (76J20)

Cites Work

• The numerical simulation of two-dimensional fluid flow with strong shocks
• Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy
• Accurate monotonicity-preserving schemes with Runge-Kutta time stepping
• Implicit Large Eddy Simulation
• The piecewise parabolic method (PPM) for gas-dynamical simulations
• A limiter for PPM that preserves accuracy at smooth extrema
• A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations
• The \(\nabla \cdot B=0\) constraint in shock-capturing magnetohydrodynamics codes
• Towards the ultimate conservative difference scheme. II: Monotonicity and conservation combined in a second-order scheme
• Title not available (Why is that?)
• Accurate monotonicity- and extrema-preserving methods through adaptive nonlinear hybridizations
• A contribution to the great Riemann solver debate
• An upwind differencing scheme for the equations of ideal magnetohydrodynamics
• The effect of nonzero \(\bigtriangledown\cdot B\) on the numerical solution of the magnetohydrodynamic equations
• An unsplit staggered mesh scheme for multidimensional magnetohydrodynamics
• An approximate Riemann solver for ideal magnetohydrodynamics
• An unsplit Godunov method for ideal MHD via constrained transport
• Development of low dissipative high order filter schemes for multiscale Navier-Stokes/MHD systems
• A multi-state HLL approximate Riemann solver for ideal magnetohydrodynamics
• An adaptive moving mesh method for two-dimensional ideal magnetohydrodynamics
• Extension of the piecewise parabolic method to multidimensional ideal magnetohydrodynamics
• A Third-Order Accurate Variation Nonexpansive Difference Scheme for Single Nonlinear Conservation Laws
• 10.1007/s11470-008-3011-1
• High resolution staggered mesh approach for nonlinear hyperbolic systems of conservation laws
• High resolution Godunov-type schemes with small stencils
• A high order staggered grid method for hyperbolic systems of conservation laws in one space dimension

Cited In (10)

• The piecewise parabolic method (PPM) for gas-dynamical simulations
• Modelling supernova-driven turbulence
• Numerical simulation of three-dimensional quasi-neutral gas flows based on smoothed magnetohydrodynamic equations
• 10.1007/s11470-008-3011-1
• Piecewise parabolic method on a local stencil in cylindrical coordinates for fluid dynamics simulations
• A new code for the numerical simulation of relativistic flows on supercomputers by means of a low-dissipation scheme
• Efficient implementation of ADER schemes for Euler and magnetohydrodynamical flows on structured meshes -- speed comparisons with Runge-Kutta methods
• Modern methods of mathematical modeling of the development of hydrodynamic instabilities and turbulent mixing
• The structure and statistics of interstellar turbulence
• Smoothed MHD equations for numerical simulations of ideal quasi-neutral gas dynamic flows

This page was built for publication: Piecewise parabolic method on a local stencil for magnetized supersonic turbulence simulation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q733016)