Simulation of viscoelastic free-surface flows with the Particle Finite Element Method

Viscoelastic fluids are central in numerous applications from polymer manufacturing to the pharmaceutical industry and biological research. However, since analytical solutions are generally not available or too complex, it is common practice to study free-surface viscoelastic flows through numerical simulation techniques. This work proposes the use of the so-called particle finite element method (PFEM), a Lagrangian approach combining standard FEM techniques with a remeshing strategy. The PFEM is able to efficiently handle mesh distortion and to accurately track the free-surface evolution. Therefore, it is exploited in this work to deal with large displacements problems in the context of nonlinear viscoelasticity. An implementation of the Oldroyd-B constitutive model in the PFEM framework is here presented including details regarding how to deal with the transfer of the internal variables during remeshing events. Additionally, an innovative approach to impose unilateral Dirichlet boundary conditions ensuring optimal mass conservation is presented. The implementation is verified with two free-surface highly viscous benchmark flows: the impacting drop and the jet buckling problems. The results show perfect agreement with those obtained with other numerical techniques. The proposed framework opens the way for using PFEM in various applications, ranging from polymer extrusion to more sophisticated scenarios involving viscoelastic and viscoelasto-plastic constitutive laws.