Surface tension problems solved with the particle finite element method using large time-steps

Abstract

In previous works [S. R. Idelsohn, J. Marti, P. Becker, E. Oñate, Analysis of multifluid flows with large time steps using the particle finite element method, International Journal for Numerical Methods in Fluids 75 (9) (2014) 621–644. doi:10.1002/fld.3908. URL http://dx.doi.org/10.1002/fld.3908, Juan M. Gimenez and Leo M. González, An extended validation of the last generation of particle finite element method for free surface flows, J Comput Phys 284 (0) (2015) 186–205. doi:http://dx.doi.org/10.1016/j.jcp.2014.12.025. URL http://www.sciencedirect.com/science/article/pii/S0021999114008420 ], the authors have presented a highly efficient extension of the Particle Finite Element Method, called PFEM-2, to solve two-phase flows. The methodology which uses X-IVS [S. Idelsohn, N. Nigro, A. Limache, E. Oñate, Large time-step explicit integration method for solving problems with dominant convection, Comp Methods in Appl Mech Eng 217–220 (2012) 168–185.] to treat convection terms allowing large time-steps was validated for problems where the gravity forces and/or the inertial forces dominate the flow. Although that is the target range of problems to solve with PFEM-2, most of real problems that fall in these categories also includes other flow regimes in certain regions of the domain. Maybe the most common secondary regime is when the surface tension dominates, as an example when drops or bubbles are released from the main flow, and this feature must be taken into account in any complete numerical strategy.\ud \ud Attending to that, in this work the treatment of the surface tension to PFEM-2 is included. An implicit CSF methodology is employed together with a coupling between the marker function with a Level Set function to obtain a smooth representation of the normal of the interface which allows an accurate curvature calculation. Examples for curvature calculation and isolated bubbles and drops are presented where the accuracy and the computational efficiency are analyzed and contrasted with other numerical methodologies. Finally, a simulation of a jet atomization is analyzed. This case presents the above mentioned features: it is a inertia-dominant flow with a surface tension phenomena on drops and ligaments break up that can not be neglected.Peer ReviewedPostprint (author's final draft

A State of the Art Review of the Particle Finite Element Method (PFEM)

“…Gimenez et al [37] shows the potential of PFEM-2 to simulate industrial problems of large time duration. An application of the method to jet atomization simulation can be found in [39]. Becker and Idelsohn [5] shows the potential of PFEM-2 to simulate large scale landslides events.…”

A finite element level set method based on adaptive octree meshes for thermal free‐surface flows

“…Interface-capturing method has advantages for problems with large interface deformation, because it can automatically deal with topological changes. Interface-capturing method has been applied to challenging engineering problems in marine engineering, 16,17 jet atomization 18 and bubble dynamics. [19][20][21] In this paper, we utilize the level set interface-capturing method based on a large-scale finite element framework for modeling free-surface flows.…”

An immersogeometric formulation for free-surface flows with application to marine engineering problems

“…The free-surface topological changes can be automatically handled by solving an additional scalar partial differential equation. Interface capturing methods have been widely applied to a wide range of interfacial problems, including bubble dynamics [17][18][19], jet atomization [20], and free-surface flows [21,22].…”