POST ADAPTATION FOR A NUMERICAL SOLUTION OF THE SPHERICALLY-SYMMETRIC RIEMANN PROBLEM

A"folding" transformation is used to reduce a two point boundary value problem to one point boundary values (with double the number of functions). A transformation to quasi-Lagrangian coordinates is used to transform discontinuities to rest. After these transformations the general purpose numerical package POST (partial and ordinary differential equations solver in time and one space coordinates) can be successfully applied to a system of hydrodynamical equations, whose solution exhibits jumps and cusp-type discontinuities. Numerical results are presented for the spherically-symmetric shock problem.