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.

A kind of numerical simulation algorithms for instabilities of fluid interfaces,Ghost Fluid Method, is discussed.The setting of fluid variables in ghost regions is directly based on multiple dimensional cases.Moreover,non-splitting high resolution schemes are used to solve the fluid dynamic equations. At last, the algorithm is applied to the simulation of interface instabilities of fluid dyanmics,and satisfied numerical results are obtained.

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.

Computation of three dimensional fluid interface instability with random perturbation is performed on a 8 CPU personal parallel computer system.The parallel computation platform is a message passage interface (MPI).A second order TVD scheme with Ghost method is applied with a fully parallel algorithm to the 3D Euler equations.Level Set function is used to track the motion of a fluid interface in an Eulerian framework.Buffer zone and data communication are discussed.The number of computation meshes is about 1 million.Bubble evolution with random perturbation in Rayleigh Taylor instability is obtained by numerical simulation.

It uses a Level Set function to track the motion of a fluid interface in an Eulerian framework.In addition,the use of Ghost cells (actually Ghost nodes in the difference framework) can keep the density profile from smearing out,while still preserving the scheme robust 2-nd easy to program with and accuracy TVD scheme in spacial and two step Runge Kutta methods.In contrast, the references[6][7] allow the density profile smearing out.This method significantly improves the resolution at the contact discontinuity.