We discuss high-resolution simulations of one-dimensional, plane-parallel shock waves with mean speeds between 150 and 240 km/s propagating into gas with Alfven velocities up to 40 km/s and outline the conditions under which these radiative shocks experience an oscillatory instability in the cooling length, shock velocity, and position of the shock front. We investigate two forms of postshock cooling: a truncated single power law and a more realistic piecewise power law. The degree of nonlinearity of the instability depends strongly on the cooling power law and the Alfven Mach number: for power-law indices \alpha < 0 typical magnetic field strengths may be insufficient either to stabilize the fundamental oscillatory mode or to prevent the oscillations from reaching nonlinear amplitudes.

High-resolution simulations of the cooling regions of spherically symmetric supernova remnants demonstrate a strong radiative instability. This instability, whose presence is dependent on the shock velocity, causes large-amplitude fluctuations in the shock velocity. The fluctuations begin almost immediately after the radiative phase begins (upon shell formation) if the shock velocity lies in the unstable range; they last until the shock slows to speeds less than approximately 130 km/s. We find that shock-velocity fluctuations from the reverberations of waves within the remnant are small compared to those due to the instability. Further, we find (in plane-parallel simulations) that advected inhomogeneities from the external medium do not interfere with the qualitative nature of the instability-driven fluctuations. Large-amplitude inhomogeneities may alter the phases of shock-velocity fluctuations, but do not substantially reduce their amplitudes.

We investigate asymptotic convergence in the~$\Delta x \!\rightarrow\! 0$ limit as a tool for determining whether numerical computations involving shocks are accurate. We use one-dimensional operator-split finite-difference schemes for hydrodynamics with a von Neumann artificial viscosity. An internal-energy scheme converges to demonstrably wrong solutions. We associate this failure with the presence of discontinuities in the limiting solution. Our extension of the Lax-Wendroff theorem guarantees that certain conservative, operator-split schemes converge to the correct continuum solution. For such a total-energy scheme applied to the formation of a single shock, convergence of a Cauchy error approaches the expected rate slowly. We relate this slowness to the effect of varying diffusion, due to varying linear artificial-viscous length, on small-amplitude waves. In an appendix we discuss the scaling of shock-transition regions with viscous lengths, and exhibit several difficulties for attempts to make extrapolations.