A method to simulate linear stability of impulsively accelerated density interfaces in ideal-MHD and gas dynamics

  • Authors:
  • Ravi Samtaney

  • Affiliations:
  • Princeton Plasma Physics Laboratory, Princeton University, P.O. Box 451, Princeton, NJ 08543-0451, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2009

Quantified Score

Hi-index 31.45

Visualization

Abstract

We present a numerical method to solve the linear stability of impulsively accelerated density interfaces in two dimensions such as those arising in the Richtmyer-Meshkov instability. The method uses an Eulerian approach, and is based on an upwind method to compute the temporally evolving base state and a flux vector splitting method for the perturbations. The method is applicable to either gas dynamics or magnetohydrodynamics. Numerical examples are presented for cases in which a hydrodynamic shock interacts with a single or double density interface, and a doubly shocked single density interface. Convergence tests show that the method is spatially second-order accurate for smooth flows, and between first and second-order accurate for flows with shocks.