It presents numerical computation results of the refraction of a plane shock wave at the different interfaces. The unsteady, two-dimensional, compressible, Euler equations are solved numerically assuming stiffened gas equation of state. The numerical method used is PPM (Piecewise Parabolic Method) method for multi-component. On basis of solving the Riemann problem using two-shock approximation the interface can be determined by computing the volume-fraction in terms of the high-order PPM method in the whole computational field. This method is also used to compute some case of the inter-action when the shock wave meets the interfaces of the different fluids. It also gives the evolvement of interfaces due to baroclinic effects on the interface, especially in the cases of strong shock striking at the interfaces.

Two-dimensional compressible flow field for the interaction of shocks with cylinder interface is directly simulated by using Ghost Fluid method (Ghost) and γ-model method, respectively, with the same discrete order in space and time. Numerical results are compared with the experimental results. They are almost alike in the beginning time and show the right position of interface, right strength and velocity of shocks. With the time developing, the effect of large numerical dissipation of γ-model method becomes greater and greater while the low numerical dissipation of Ghost method makes it efficient to simulate the moving interface well. Comparison with experiments proves that Ghost method is better than γ-model method in simulating the problem of interface instability.

An approach combining the fifth Upwind Compact Difference Scheme and Group Velocity Control is applied to simulate directly the interaction of shocks and the density stratified interface of light-gas cylinder.Some research about the generation and development of vorticity and the distortion of edge of cylinder is done.The instability of cylinder density-interface is analyzed.Main character of the linear developing phase is that the gas of higher density goes into the one of lower denstiy at a constant velocity and the spike structure is formed and the character of the nonlinear developing phase is that a pair of vortex structures rotating adversely are formed on both sides of the top of the spike.The numerical solution accords with the experiment.