The Great East Japan Earthquake of March 2011 was the huge earthquake with a magnitude of 9.0 exceeding the imagination, ranging from the offshore to Iwate Prefecture to the offshore of Ibaraki Prefecture became the focal region. It is easy to imagine that damage increases to long-period earthquake ground motion as the magnitude increases. Damage to dangerous facilities due to long-period earthquake ground motion is caused by liquid surface fluctuation (sloshing) of oil, LNG and liquid storage tanks. In this paper, the interface capturing method deals with complex water surface large deformation problems. We use the density function method which is one of the effective methods of interface capturing method for numerical analysis of two phase fluid fields of gas and liquid. As a result, the proposed method shows a good agreement with the experimental result of the Martin and Moyce dam break problem, and this method improved the suppression of effortless motion. Also, under the calculation condition of the sloshing phenomenon, this method does not generate numerical diffusion. Therefore, this model is more robust than the conventional model.