Simulation of viscoelastic two-phase flows with insoluble surfactants

A stabilized finite element scheme is developed for computations of buoyancy driven 3D-axisymmetric viscoelastic two-phase flows with insoluble surfactants. The numerical scheme solves the Navier-Stokes equations for the fluid flow, Giesekus constitutive equation for the effects of viscoelasticity and simultaneously an evolution equation for the surfactant concentration on the interface. The interface is tracked by the coupled arbitrary Lagrangian-Eulerian (ALE) and Lagrangian approach. The interface-resolved moving meshes allow accurate incorporation of the interfacial tension force, Marangoni forces and the jumps in the material properties. Further, the tangential gradient operator technique is used to handle the curvature approximation in a semi-implicit manner. An one-level Local Projection Stabilization (LPS), which is based on an enriched approximation space and a discontinuous projection space, where both spaces are defined on a same mesh is used to stabilize the model equations. The stabilized numerical scheme allows us to use isoparametric second order conforming finite elements enriched with cubic bubble functions for velocity and viscoelastic stress, second order finite elements for surfactant concentration and discontinuous first order finite element for pressure. A number of computations are performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column with insoluble surfactants on the interface. The influence of the Marangoni number, initial surfactant concentration and Peclet number on the dynamics of the rising drop are analyzed. The numerical study shows that a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape without insoluble surfactants. The presence of insoluble surfactants forces the drop to rise slowly but the drop at the tail end is pulled up more. However, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape without insoluble surfactants. The presence of surfactants pulls the tail end of the drop up slightly and makes the tail flatter with/without small undulations depending on the magnitude of the surfactant concentrations.