Formal analysis of the bubble phenomena of the second order relaxation scheme is presented when this method is applied to phase transition equations. The reason of oscillations is that the second order scheme can't find the discontinuities on the phase boundary. Based on this realization, a second order relaxation scheme is derived to eliminate it. This new method finds all components and characteristic discontinuities, thus the phase boundary is found exactly. The difference of presented new method and the well-known Nessyahu-Tadmor(NT) scheme is also studied. From the numerical experiment, the new derived scheme is shown much better to compute the phase transition equations.
展开▼
