Multiphase Flows Simulation with the Particle Finite Element Method and its Comparison with Eulerian Alternatives

The latest version of the Particle-Finite Element Method (PFEM), which incorporates the novel explicit integration strategy named eXplicit Integration of Velocity and Acceleration following Streamlines (X-IVAS), has proven to be fast and accurate to solve homogeneous flows, mainly thanks to the possibility of using large time-steps. In this work the extension of this strategy to solve multiphase flows is presented, where the calculation of the interface evolution is of fundamental importance. In Eulerian formulations, one of the most used strategies to determine the interface position is the advection of an indicator function. This approach is followed, by example, in the Volume of Fluid (VoF) technique, which can add limiters as a method of guaranteeing boundedness and/or sharpness of phasefractions. On the other hand, Lagrangian frames use typically marker particles. In the case of PFEM, the same set of particles transported for flow calculation allows to carry a marker function to determinate the interface position without any extra cost. In order to compare the accuracy of PFEM interface evolution strategy with the Eulerian one implemented in the widely used OpenFOAM®suite, several classical tests are presented. In addition to capture the sharpen interface evolution, PFEM algorithm includes the use of an enriched finite element space to avoid spurious solutions due to the discontinuity of pressure gradients, which also requires some changes in the streamline integration strategy. The description of those improvements are also presented in this work. Finally the PFEM algorithm is tested for a number of problems involving free surface flows with different ratio between the densities and viscosities of the fluids involved. The accuracy of the results are compared with reference ones, focusing in the capability of using large time steps in contrast with Eulerian solvers, including those with implicit phase-fractions advection treatment, represented by the OpenFOAM®suite.

[1]  Santiago Márquez Damián,et al.  An extended mixture model for the simultaneous treatment of short and long scale interfaces , 2013 .

[2]  Akira Umemura,et al.  Surface instability and primary atomization characteristics of straight liquid jet sprays , 2011 .

[3]  Heinz Pitsch,et al.  An accurate conservative level set/ghost fluid method for simulating turbulent atomization , 2008, J. Comput. Phys..

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

[5]  Sergio Idelsohn,et al.  A fast and accurate method to solve the incompressible Navier‐Stokes equations , 2013 .

[6]  A. Umemura,et al.  Simulation of liquid jet primary breakup: Dynamics of ligament and droplet formation , 2010 .

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

[8]  Joseph J Monaghan,et al.  An introduction to SPH , 1987 .

[9]  Francis H. Harlow,et al.  The Particle-In-Cell Method for Hydrodynamic Calculations , 1957 .

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

[11]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

[12]  Eugenio Oñate,et al.  Large time-step explicit integration method for solving problems with dominant convection , 2012 .

[13]  Leo M. González,et al.  Two phase analysis of sloshing in a rectangular container with Volume of Fluid (VOF) methods , 2013 .

[14]  G. Tryggvason.  Numerical simulations of the Rayleigh-Taylor instability , 1988 .

[15]  P. Koumoutsakos,et al.  A Lagrangian particle level set method. , 2005 .

[16]  E. Sergent,et al.  Vers une méthodologie de couplage entre la simulation des grandes échelles et les modèles statistiques , 2002 .

[17]  Eugenio Oñate,et al.  The particle finite element method: a powerful tool to solve incompressible flows with free‐surfaces and breaking waves , 2004 .

[18]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[19]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .