Simulation of viscoelastic two-phase flows with insoluble surfactants

Venkatesan, Jagannath and Padmanabhan, Adhithya and Ganesan, Sashikumaar (2019) Simulation of viscoelastic two-phase flows with insoluble surfactants. In: JOURNAL OF NON-NEWTONIAN FLUID MECHANICS, 267 . pp. 61-77.

PDF Jou_Non_New_Flu_Mec_267_61-77_2019.pdf - Published Version Restricted to Registered users only Download (11MB) | Request a copy

Abstract

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.

Item Type:	Journal Article
Publication:	JOURNAL OF NON-NEWTONIAN FLUID MECHANICS
Publisher:	ELSEVIER SCIENCE BV
Additional Information:	copyright for this article belongs to JOURNAL OF NON-NEWTONIAN FLUID MECHANICS
Keywords:	Viscoelastic fluids; Giesekus model; Insoluble surfactants; Finite elements; ALE approach; Local projection stabilization
Department/Centre:	Division of Interdisciplinary Sciences > Computational and Data Sciences
Date Deposited:	24 Jun 2019 18:00
Last Modified:	24 Jun 2019 18:00
URI:	http://eprints.iisc.ac.in/id/eprint/63071

Actions (login required)

View Item