Abstract
A method of treating a sharp discontinuity with a density function is proposed. The surface of the density function is captured within a grid cell throughout the calculation even when the discontinuity surface is largely distorted. This description is made possible by the CIP (cubic-interpolated propagation) method combined with variable transformation. This scheme can be used for detecting a sharp boundary in compressible fluid as well as incompressible fluid. Lower order schemes with the same procedure cannot reproduce the result. This scheme is applied to the interaction of a shock wave with a liquid drop. The liquid mass is conserved within an error of 0.15%.