Abstract:
Application of a method for dynamic adaptation to numerical solution of unsteady-state hyperbolic equations is considered. The basic for given method is a transition to arbitrary nonstationary coordinate system. This allows to formulate a unified difference model to determine both numerical solution and grid movement. Two different techniques for computation of shock waves by means of dynamic adaptation are examined by the example of numerical solution for a model problem which describes accelerating piston motion in gas. The first technique enables the discontinuity surface to show up as a region of high gradients which grid nodes are concentrated in. The second one permits to separate the discontinuity explicitly by means of Rankine–Hugoniot equations. Both of these techniques don't use artificial viscosity and allow to proceed computations when using grids with small number of nodes ($N\approx10-20$).