Stability analysis of a predictor/multi-corrector method for staggered-grid Lagrangian shock hydrodynamics

  • Authors:

  • E. Love;W. J. Rider;G. Scovazzi

  • Affiliations:

  • 1431 Computational Shock- and Multi-physics Department, Sandia National Laboratories, P.O. Box 5800, MS 1319, Albuquerque, NM 87185-1319, USA;1431 Computational Shock- and Multi-physics Department, Sandia National Laboratories, P.O. Box 5800, MS 1319, Albuquerque, NM 87185-1319, USA;1431 Computational Shock- and Multi-physics Department, Sandia National Laboratories, P.O. Box 5800, MS 1319, Albuquerque, NM 87185-1319, USA

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

Quantified Score

Hi-index 31.45

Visualization

Abstract

This article presents the complete von Neumann stability analysis of a predictor/multi-corrector scheme derived from an implicit mid-point time integrator often used in shock hydrodynamics computations in combination with staggered spatial discretizations. It is shown that only even iterates of the method yield stable computations, while the odd iterates are, in the most general case, unconditionally unstable. These findings are confirmed by, and illustrated with, a number of numerical computations. Dispersion error analysis is also presented.