In this work a new solver is developed for the OpenFOAM (R) 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. (C) 2011 Published by Elsevier Ltd.     «

In this work a new solver is developed for the OpenFOAM (R) 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 a...     »