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

Title: Shock-capturing approach and nonevolutionary solutions in magnetohydrodynamics

Abstract

Shock-capturing methods have become an effective tool for the solution of hyperbolic partial differential equations. Both upwind and symmetric TVD schemes in the framework of the shock-capturing approach are thoroughly investigated and applied with great success to a number of complicated multidimensional gasdynamic problems. The extension of these schemes to magnetohydrodynamic (MHD) equations is not a simple task. First, the exact solution of the MHD Riemann problem is too multivariant to be used in regular calculations. On the other hand, the extensions of Roe`s approximate Riemann problem solvers for MHD equations in general case are nonunique and need further investigation. That is why, some simplified approaches should be constructed. In this work, the second order of accuracy in time and space high-resolution Lax-Friedrichs type scheme is suggested that gives a drastic simplification of the numerical algorithm comparing to the precise characteristic splitting of Jacobian matrices. The necessity is shown to solve the full set of MHD equations for modeling of multishocked flows, even when the problem is axisymmetric, to obtain evolutionary solutions. for the numerical example, the MHD Riemann problem is used with the initial data consisting of two constant states lying to the right and to the left frommore » the centerline of the computational domain. If the problem is solved as purely coplanar, a slow compound wave appears in the self-similar solution obtained by any shock-capturing scheme. If the full set of MHD equations is used and a small uniform tangential disturbance is added to the magnetic field vector, a rotational jump splits from the compound wave, and it degrades into a slow shock. The reconstruction process of the nonevolutionary compound wave into evolutionary shocks is investigated. Presented results should be taken into account in the development of shock-capturing methods for MHD flows. 18 refs., 17 figs.« less

Authors:
 [1];  [2];  [3]
  1. Moscow Univ. (Russian Federation)
  2. Steklov Institute of Mathematics, Moscow (Russian Federation)
  3. Institute for Problems in Mechanics, Moscow (Russian Federation)
Publication Date:
OSTI Identifier:
440757
Resource Type:
Journal Article
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 126; Journal Issue: 1; Other Information: PBD: Jun 1996
Country of Publication:
United States
Language:
English
Subject:
66 PHYSICS; 99 MATHEMATICS, COMPUTERS, INFORMATION SCIENCE, MANAGEMENT, LAW, MISCELLANEOUS; MAGNETOHYDRODYNAMICS; PLASMA FLUID EQUATIONS; NUMERICAL SOLUTION; PARTIAL DIFFERENTIAL EQUATIONS; SHOCK WAVES

Citation Formats

Barmin, A A, Kulikovskiy, A G, and Pogorelov, N V. Shock-capturing approach and nonevolutionary solutions in magnetohydrodynamics. United States: N. p., 1996. Web. doi:10.1006/jcph.1996.0121.
Barmin, A A, Kulikovskiy, A G, & Pogorelov, N V. Shock-capturing approach and nonevolutionary solutions in magnetohydrodynamics. United States. https://doi.org/10.1006/jcph.1996.0121
Barmin, A A, Kulikovskiy, A G, and Pogorelov, N V. 1996. "Shock-capturing approach and nonevolutionary solutions in magnetohydrodynamics". United States. https://doi.org/10.1006/jcph.1996.0121.
@article{osti_440757,
title = {Shock-capturing approach and nonevolutionary solutions in magnetohydrodynamics},
author = {Barmin, A A and Kulikovskiy, A G and Pogorelov, N V},
abstractNote = {Shock-capturing methods have become an effective tool for the solution of hyperbolic partial differential equations. Both upwind and symmetric TVD schemes in the framework of the shock-capturing approach are thoroughly investigated and applied with great success to a number of complicated multidimensional gasdynamic problems. The extension of these schemes to magnetohydrodynamic (MHD) equations is not a simple task. First, the exact solution of the MHD Riemann problem is too multivariant to be used in regular calculations. On the other hand, the extensions of Roe`s approximate Riemann problem solvers for MHD equations in general case are nonunique and need further investigation. That is why, some simplified approaches should be constructed. In this work, the second order of accuracy in time and space high-resolution Lax-Friedrichs type scheme is suggested that gives a drastic simplification of the numerical algorithm comparing to the precise characteristic splitting of Jacobian matrices. The necessity is shown to solve the full set of MHD equations for modeling of multishocked flows, even when the problem is axisymmetric, to obtain evolutionary solutions. for the numerical example, the MHD Riemann problem is used with the initial data consisting of two constant states lying to the right and to the left from the centerline of the computational domain. If the problem is solved as purely coplanar, a slow compound wave appears in the self-similar solution obtained by any shock-capturing scheme. If the full set of MHD equations is used and a small uniform tangential disturbance is added to the magnetic field vector, a rotational jump splits from the compound wave, and it degrades into a slow shock. The reconstruction process of the nonevolutionary compound wave into evolutionary shocks is investigated. Presented results should be taken into account in the development of shock-capturing methods for MHD flows. 18 refs., 17 figs.},
doi = {10.1006/jcph.1996.0121},
url = {https://www.osti.gov/biblio/440757}, journal = {Journal of Computational Physics},
number = 1,
volume = 126,
place = {United States},
year = {Sat Jun 01 00:00:00 EDT 1996},
month = {Sat Jun 01 00:00:00 EDT 1996}
}