A multiphase MPS solver for modeling multi-fluid interaction with free surface and its application in oil spill

Abstract Numerical simulation of free surface multi-fluid flow is a challenging problem owning to the interaction between deformed interface and free surface. In this paper, a multiphase particle solver (free surface multiphase moving particle semi-implicit method, FS-MMPS) is developed to predict the early spreading flow of spilled oil where exist the oil–water interface and air–oil/air–water free surface. First, a multiphase virtual particle model is proposed to substitute the inaccurate free surface boundary condition for multiphase flow in conventional MPS methods. Specifically, virtual particles of different liquid phases are compensated outside free surface so that the pressure of free surface particles can be solved from pressure Poisson equation, thereby improving the accuracy of multi-fluid interaction at free surface. Meanwhile, a pressure gradient model based on the coupling of Taylor series expansion and dynamic specification of particle stabilizing term (PST) is proposed to simultaneously enhance accuracy and depress instability caused by multiphase virtual particles. Experiments of early spreading of thick oil slicks and continuous oil spill from a damaged tank are conducted for validation and demonstration of the accuracy and stability enhancements in 2D. Finally, the effects of traveling wave and continuous spilling on oil spreading are investigated to show the capability of the proposed method. The spreading of thick oil slicks is accelerated when wave height increases or wave length decreases, and the largest increase of spreading rate accelerated by traveling wave is 14%.

[1]  Damien Violeau,et al.  Numerical modelling of boom and oil spill with SPH , 2007 .

[2]  Alan A. Allen,et al.  Advantages and disadvantages of burning spilled oil , 1993 .

[3]  H. M. Cekirge,et al.  Empirical studies of the spreading of oil spills , 1984 .

[4]  J. Monaghan.  Simulating Free Surface Flows with SPH , 1994 .

[5]  G. Dilts.  MOVING-LEAST-SQUARES-PARTICLE HYDRODYNAMICS-I. CONSISTENCY AND STABILITY , 1999 .

[6]  Mehmet Yildiz,et al.  Numerical simulation of single droplet dynamics in three-phase flows using ISPH , 2013, Comput. Math. Appl..

[7]  L. Lucy.  A numerical approach to the testing of the fission hypothesis. , 1977 .

[8]  Tiangang Zhang,et al.  Improvement of boundary conditions for non-planar boundaries represented by polygons with an initial particle arrangement technique , 2016 .

[9]  Abbas Khayyer,et al.  Modified Moving Particle Semi-implicit methods for the prediction of 2D wave impact pressure , 2009 .

[10]  Seiichi Koshizuka,et al.  Improvement of stability in moving particle semi‐implicit method , 2011 .

[11]  J. Monaghan,et al.  Smoothed particle hydrodynamics: Theory and application to non-spherical stars , 1977 .

[12]  Pavel Tkalich,et al.  The Numerical Investigation on Oil Slick Behavior behind the Oil Boom , 2014 .

[13]  Hitoshi Gotoh,et al.  Current achievements and future perspectives for projection-based particle methods with applications in ocean engineering , 2016 .

[14]  Hitoshi Gotoh,et al.  Enhancement of stability and accuracy of the moving particle semi-implicit method , 2011, J. Comput. Phys..

[15]  A. Colagrossi,et al.  Viscous bubbly flows simulation with an interface SPH model , 2013 .

[16]  Guoqian Chen,et al.  Numerical analysis of a lock-release oil slick , 2010 .

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

[18]  David Le Touzé,et al.  An Hamiltonian interface SPH formulation for multi-fluid and free surface flows , 2009, J. Comput. Phys..

[19]  Abbas Khayyer,et al.  Space potential particles to enhance the stability of projection-based particle methods , 2015 .

[20]  Seiichi Koshizuka,et al.  Current Achievements and Future Perspectives on Particle Simulation Technologies for Fluid Dynamics and Heat Transfer , 2011 .

[21]  Tao Jiang,et al.  Numerical analysis of the impact of two droplets with a liquid film using an incompressible SPH method , 2014 .

[22]  Torstein K. Fannelop,et al.  DYNAMICS OF OIL SLICKS , 1972 .

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

[24]  Kohei Murotani,et al.  Improved pressure calculation for the moving particle semi-implicit method , 2015, CPM 2015.

[25]  Yee-Chung Jin,et al.  MPS mesh-free particle method for multiphase flows , 2012 .

[26]  Fabian Denner,et al.  Fully-Coupled Balanced-Force VOF Framework for Arbitrary Meshes with Least-Squares Curvature Evaluation from Volume Fractions , 2014 .

[27]  Hitoshi Gotoh,et al.  Key issues in the particle method for computation of wave breaking , 2006 .

[28]  Hitoshi Gotoh,et al.  Enhancement of performance and stability of MPS mesh-free particle method for multiphase flows characterized by high density ratios , 2013, J. Comput. Phys..

[29]  J. Fay.  Physical Processes in the Spread of Oil on a Water Surface , 1971 .

[30]  Moo-Hyun Kim,et al.  Step-by-step improvement of MPS method in simulating violent free-surface motions and impact-loads , 2011 .

[31]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[32]  Gary A. Dilts,et al.  Moving least‐squares particle hydrodynamics II: conservation and boundaries , 2000 .

[33]  Jong-Chun Park,et al.  Dynamic Coupling Between Ship Motion and Three-Layer-Liquid Separator by Using Moving Particle Simulation , 2014 .

[34]  Moubin Liu,et al.  Numerical modeling of oil spill containment by boom using SPH , 2013 .

[35]  Joseph V. Mullin,et al.  Introduction/Overview to In Situ Burning of Oil Spills , 2003 .

[36]  P. G. Jeffery,et al.  Large-Scale Experiments on the Spreading of Oil at Sea and its Disappearance by Natural Factors , 1973 .

[37]  Nikolaus A. Adams,et al.  An incompressible multi-phase SPH method , 2007, J. Comput. Phys..

[38]  Zhen Chen,et al.  An SPH model for multiphase flows with complex interfaces and large density differences , 2015, J. Comput. Phys..

[39]  Seiichi Koshizuka,et al.  Least squares moving particle semi-implicit method , 2014, Computational Particle Mechanics.

[40]  Nikolaus A. Adams,et al.  A constant-density approach for incompressible multi-phase SPH , 2009, J. Comput. Phys..

[41]  Xiaoyang Xu.  An improved SPH approach for simulating 3D dam-break flows with breaking waves , 2016 .

[42]  Guang Xi,et al.  Improving stability of MPS method by a computational scheme based on conceptual particles , 2014 .

[43]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[44]  Y. Oka,et al.  Experiments and MPS analysis of stratification behavior of two immiscible fluids , 2013 .

[45]  S. Koshizuka,et al.  International Journal for Numerical Methods in Fluids Numerical Analysis of Breaking Waves Using the Moving Particle Semi-implicit Method , 2022 .

[46]  Chan Ghee Koh,et al.  A new particle method for simulation of incompressible free surface flow problems , 2012 .

[47]  J. Monaghan.  Smoothed particle hydrodynamics , 2005 .

[48]  S. J. Lind,et al.  Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves , 2012, J. Comput. Phys..

[49]  David P. Hoult,et al.  Oil Spreading on the Sea , 1972 .

[50]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[51]  Moo-Hyun Kim,et al.  Numerical prediction of oil amount leaked from a damaged tank using two-dimensional moving particle simulation method , 2013 .

[52]  Kenji Fukuda,et al.  Simulation of the Rayleigh-Taylor Instability with the MPS Method , 2004 .

[53]  Inertia–gravity spreading of oil on water , 2000 .

[54]  Gretar Tryggvason,et al.  Direct numerical simulations of gas/liquid multiphase flows , 2011 .

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

[56]  Seiichi Koshizuka,et al.  Stable multiphase moving particle semi-implicit method for incompressible interfacial flow , 2017 .

[57]  Furen Ming,et al.  An SPH modeling of bubble rising and coalescing in three dimensions , 2015 .

[58]  Bin Chen,et al.  A contoured continuum surface force model for particle methods , 2015, J. Comput. Phys..

[59]  Nikolaus A. Adams,et al.  A multi-phase SPH method for macroscopic and mesoscopic flows , 2006, J. Comput. Phys..

[60]  Kamil Szewc,et al.  Simulations of single bubbles rising through viscous liquids using Smoothed Particle Hydrodynamics , 2013 .

[61]  Abbas Khayyer,et al.  A higher order Laplacian model for enhancement and stabilization of pressure calculation by the MPS method , 2010 .

[62]  Kh. M. Gamzaev.  Modeling the spread of an oil slick on the sea surface , 2009 .

[63]  Masayuki Tanaka,et al.  Stabilization and smoothing of pressure in MPS method by Quasi-Compressibility , 2010, J. Comput. Phys..

[64]  Yuichi Yamamoto,et al.  Analysis of Flows around a BWR Spacer by the Two-Fluid Particle Interaction Method , 2002 .

[65]  Abdulrauf Rasheed Adebayo,et al.  A sustainable approach to controlling oil spills. , 2012, Journal of environmental management.

[66]  Guangtao Duan,et al.  Large Eddy Simulation by particle method coupled with Sub-Particle-Scale model and application to mixing layer flow , 2015 .