Numerical simulation of viscoelastic two-phase flows using openFOAM®

In this work a new solver is developed for the OpenFOAM®OpenFOAM® CFD toolbox, which handles viscoelastic two-phase flows. A derivative of the volume-of-fluid (VoF) methodology is used to describe the interface. Established constitutive equations derived from kinetic theory, such as Oldroyd-B, Giesekus, FENE-P and FENE-CR, from network theory of concentrated solutions and melts, such as linear and exponential Phan-Thien–Tanner (PTT), and from reptation theory, such as Pom–Pom and XPP models, as well as multi-mode formulations are implemented in OpenFOAM. Validation of the numerical technique is performed by comparing detailed simulation predictions to data from several experimental studies, numerical studies and analytical models found in the literature. Two well-known viscoelastic free-surface effects, namely the Weissenberg and the Die Swell effect, are simulated. Furthermore, transient and steady-state droplet flow in shear and elongational flows is examined.

► Algorithm to simulate two-phase flows of viscoelastic fluids with surface tension.
► Capability of combining numerous different viscoelastic constitutive models.
► Validation and verification using several experimental, analytical and numerical studies.
► Qualitative and quantitative agreement with all reference data.