Surface tension problems solved with the Particle Finite Element Method using large time-steps

In previous works [1,2], 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 [3] 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. 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