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

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.