MLPG_R method for modelling 2D flows of two immiscible fluids

This is a first attempt to develop the Meshless Local Petrov-Galerkin method with Rankine source solution (MLPG_R method) to simulate multiphase flows. In this paper, we do not only further develop the MLPG_R method to model two-phase flows but also propose two new techniques to tackle the associated challenges. The first technique is to form an equation for pressure on the explicitly identified interface between different phases by considering the continuity of the pressure and the discontinuity of the pressure gradient (i.e. the ratio of pressure gradient to fluid density), the latter reflecting the fact that the normal velocity is continuous across the interface. The second technique is about solving the algebraic equation for pressure, which gives reasonable solution not only for the cases with low density ratio but also for the cases with very high density ratio, such as more than 1000. The numerical tests show that the results of the newly developed two-phase MLPG_R method agree well with analytical solutions and experimental data in the cases studied. The numerical results also demonstrate that the newly developed method has a second-order convergent rate in the cases for sloshing motion with small amplitudes.

[1]  Yong Zhao,et al.  A novel coupled level set and volume of fluid method for sharp interface capturing on 3D tetrahedral grids , 2010, J. Comput. Phys..

[2]  Q. W. Ma,et al.  Incompressible SPH method based on Rankine source solution for violent water wave simulation , 2014, J. Comput. Phys..

[3]  J. Simpson,et al.  Gravity currents produced by instantaneous releases of a heavy fluid in a rectangular channel , 1983, Journal of Fluid Mechanics.

[4]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[5]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[6]  M. Rudman.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

[7]  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..

[8]  Qingwei Ma,et al.  A review on approaches to solving Poisson’s equation in projection-based meshless methods for modelling strongly nonlinear water waves , 2016 .

[9]  Ryan J. Lowe,et al.  The non-Boussinesq lock exchange problem , 2000 .

[10]  Juntao Zhou,et al.  MLPG_R Method for Numerical Simulation of 2D Breaking Waves , 2009 .

[11]  Ashkan Rafiee,et al.  A simple SPH algorithm for multi‐fluid flow with high density ratios , 2013 .

[12]  Rex Britter,et al.  The propagation of intrusion fronts of high density ratios , 1993, Journal of Fluid Mechanics.

[13]  J. Monaghan,et al.  Sedimentation in homogeneous and inhomogeneous fluids using SPH , 2015 .

[14]  H. Schwaiger.  An implicit corrected SPH formulation for thermal diffusion with linear free surface boundary conditions , 2008 .

[15]  Abbas Khayyer,et al.  Corrected higher order Laplacian for enhancement of pressure calculation by projection-based particle methods with applications in ocean engineering , 2015 .

[16]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

[17]  Flekkoy,et al.  Foundations of dissipative particle dynamics , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[18]  Qingwei Ma.  MLPG Method Based on Rankine Source Solution for Simulating Nonlinear Water Waves , 2005 .

[19]  K. Szewc,et al.  Spurious interface fragmentation in multiphase SPH , 2015 .

[20]  Benedict D. Rogers,et al.  Numerical predictions of water–air wave slam using incompressible–compressible smoothed particle hydrodynamics , 2015 .

[21]  J. J. Keller,et al.  On the hydraulic lock-exchange problem , 1991 .

[22]  J. López,et al.  On the reinitialization procedure in a narrow‐band locally refined level set method for interfacial flows , 2005 .

[23]  A. Colagrossi,et al.  Numerical simulation of interfacial flows by smoothed particle hydrodynamics , 2003 .

[24]  Benedict D. Rogers,et al.  Incompressible-compressible flows with a transient discontinuous interface using smoothed particle hydrodynamics (SPH) , 2016, J. Comput. Phys..

[25]  Ryan J. Lowe,et al.  The non-Boussinesq lock-exchange problem. Part 1. Theory and experiments , 2005, Journal of Fluid Mechanics.

[26]  Rui Xu,et al.  Comparisons of weakly compressible and truly incompressible algorithms for the SPH mesh free particle method , 2008, J. Comput. Phys..

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

[28]  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..

[29]  Qingwei Ma,et al.  Numerical simulation of interaction between wind and 2D freak waves , 2010 .

[30]  T. Ménard,et al.  Coupling level set/VOF/ghost fluid methods: Validation and application to 3D simulation of the primary break-up of a liquid jet , 2007 .

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

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

[33]  Masashi Kashiwagi,et al.  Numerical simulation of violent sloshing by a CIP-based method , 2006 .

[34]  Hitoshi Gotoh,et al.  Lagrangian Two-Phase Flow Model of the Settling Behavior of Fine Sediment Dumped into Water , 2001 .

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

[36]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

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

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

[39]  Songdong Shao,et al.  Incompressible smoothed particle hydrodynamics simulation of multifluid flows , 2012 .

[40]  Q. W. Ma,et al.  Improved MLPG_R method for simulating 2D interaction between violent waves and elastic structures , 2012, J. Comput. Phys..

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

[42]  Qingwei Ma,et al.  A New Interface Identification Technique Based on Absolute Density Gradient for Violent Flows , 2018 .

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

[44]  J. Monaghan,et al.  SPH simulation of multi-phase flow , 1995 .