Shock capturing with improved numerical viscosity in primitive Euler representation

A new and simple form of artificial viscosity is proposed to eliminate dire ctional dependence. Coupled with the CIP method [Comput. Phys. Commun. 66 ( 1991) 219, 233; J. Comput. Phys, 70 (1987) 355] which is known to be able t o treat all phases of matter together, the propagation of shock waves throu gh multiphase material has been accurately calculated without fictitious pr essure pulse at the phase boundary which is often seen in the conservative Euler scheme. Although the present formulation employs the primitive Euler scheme by the CIP, it has accurately captured weak and strong shock waves i n multidimensions ensuring the conservative property with the improved arti ficial viscosity of the Neumann-Richtmyer framework [J. Appl. Phys. 21 (195 0) 232]. (C) 1999 Elsevier Science B.V. All rights reserved.