Coupling of rigid body dynamics and moving particle semi-implicit method for simulating isothermal multi-phase fluid interactions

Abstract The moving particle semi-implicit (MPS) method does not require grids for simulating fluid motions. Therefore, the MPS method can easily handle a large deformation of fluid. However, the MPS method has some difficulties in simulating transfer of momentum caused by a physical collision between different fluids because fluid particles have no mass or volume and only have weights for interacting with other particles. To overcome this inherent defect of the MPS method, rigid body dynamics is explicitly coupled with the MPS method in this study. In the first step, the MPS calculation is performed with particles which are considered to have no mass or volume. In the second step, rigid body dynamics comes into the calculation and considers the particles to have a slightly lesser diameter than the initial distance between particles. Then, physical contacts between particles are simulated with the dynamic energy conserved while the incompressibility of fluids is effectively maintained. In the single fluid region, the coupled method deals with the behavior of the particles. For the interface of the different fluids, only rigid body dynamics is used to simulate the transfer of the momentum caused by physical collisions of fluids. Through this coupling of rigid body dynamics and the MPS method, the overall stability related with the incompressibility of a fluid is comparatively increased in the single-phase fluid simulation. For the calculation of the multi-phase fluids behavior, fluids interactions can be easily treated with improving stability of the MPS calculation. In this study, collapse of water column and the isothermal fuel–coolant interaction (FCI), in which a water jet is directed into a denser fluid pool, were simulated to validate the coupling method of the MPS method and rigid body dynamics.

[1]  Yoshiaki Oka,et al.  Two-dimensional simulation of drop deformation and breakup at around the critical Weber number , 2003 .

[2]  Yoshiaki Oka,et al.  Numerical analysis of fragmentation mechanisms in vapor explosions , 1999 .

[3]  David Baraff,et al.  Fast contact force computation for nonpenetrating rigid bodies , 1994, SIGGRAPH.

[4]  Yoshiaki Oka,et al.  Numerical analysis of boiling on high heat-flux and high subcooling condition using MPS-MAFL , 2002 .

[5]  Yoshiaki Oka,et al.  A hybrid particle-mesh method for viscous, incompressible, multiphase flows , 2005 .

[6]  David Baraff,et al.  Curved surfaces and coherence for non-penetrating rigid body simulation , 1990, SIGGRAPH.

[7]  Pixar Animation Studios,et al.  Physically Based Modeling , 2001 .

[8]  C. K. Thornhill,et al.  Part IV. An experimental study of the collapse of liquid columns on a rigid horizontal plane , 1952, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[9]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[10]  David Baraff,et al.  Analytical methods for dynamic simulation of non-penetrating rigid bodies , 1989, SIGGRAPH.

[11]  S. Koshizuka,et al.  Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid , 1996 .

[12]  Yoshiaki Oka,et al.  Numerical Analysis of Droplet Breakup Behavior using Particle Method , 2001 .

[13]  Seiichi KOSHIZUKA,et al.  Numerical Analysis of Jet Injection Behavior for Fuel-Coolant Interaction using Particle Method , 2001 .

[14]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .