A new potential interface tension model for MPS method avoiding unphysical particle cohesion

Improvement of the capability to accurately track the interface between immiscible phases/components is of great importance to the numerical analysis with complex multi-component and multi-phase flow dynamics, which might be involved in severe accidents of light water reactors, such as Molten Core-Concrete Interaction (MCCI). In order to model the interface tension between immiscible phases/components, Zhu et, al., precedingly developed a potential interface tension model in MPS method based on Kondo et al.‘s earlier model for surface tension and succeeded to stably model the interface tension qualitatively by introducing different attractive forces for different pairs of particle types. However, it has been known that pressure is overestimated in Zhu et, al.‘s potential interface tension model. Moreover, this study has shown that the fluid flow can be numerically hindered with Zhu et al.‘s potential interface tension model. These issues of Zhu et al.‘s potential interface tension model is caused by unphysical particle cohesion due to the inter-particle attractive forces determined by the surface tension coefficients, which are required even when the free surface is not considered. Meanwhile, in the framework of SPH method, Zhou et, al., proposed to model interface tension by introducing the repulsive force between different phases. Hence, this study aimed to develop a new potential interface tension model for MPS method which can overcome the conventional drawbacks of pressure overestimation or unphysical flow hindrance, by incorporating the repulsive force between immiscible fluid phases referring to Zhou et, al.‘s idea within SPH method. The interface tension in the present method is modeled as the repulsive force derived from the “repulsive potential” and is considered only near the interface to avoid the unphysical particle cohesion due to the inter-particle attractive forces. The amplitude coefficient in the repulsive potential is determined by the methodology similar to that proposed by Kondo, et, al. With the newly developed interface tension model, the common issue of numerical flow hindrance that occurred in the previous interface tension models in MPS method has been resolved. The validation of the newly developed interface tension model shows good agreement with theoretical evaluations of Laplace pressure and the droplet deformation ratio in a shear flow. Moreover, the new interface tension model has been validated against a 3-dimensional immiscible fluid stratification test, in which the new interface tension model shows good agreement with the test results and improvement of the pressure profile evaluations compared to that by Zhu et al.‘s model. Further extension of the present model so as to incorporate surface tension and wall wetting is one of the most important issues which require future development.