Abstract:Since viscoelastic two-phase flows arise in various industrial and natural processes, developing accurate and efficient software for their detailed numerical simulation is a highly relevant and challenging research task. We present a geometrical unstructured Volume-of-Fluid (VOF) method for handling two-phase flows with viscoelastic liquid phase, where the latter is modeled via generic rate-type constitutive equations and a one-field description is derived by conditional volume averaging of the local instantaneous bulk equations and interface jump conditions. The method builds on the plicRDF-isoAdvector geometrical VOF solver that is extended and combined with the modular framework DeboRheo for viscoelastic computational fluid dynamics (CFD). A piecewise-linear geometrical interface reconstruction technique on general unstructured meshes is employed for discretizing the viscoelastic stresses across the fluid interface. DeboRheo facilitates a flexible combination of different rheological models with appropriate stabilization methods to address the high Weissenberg number problem.

What problem does this paper attempt to address?

The paper attempts to address the significant challenges faced in direct numerical simulation (DNS) of viscoelastic two-phase fluids. These challenges mainly arise from the abrupt changes in physical properties and rheological behavior at the interface of the two-phase fluids, as well as the additional numerical difficulties brought by the strong elastic components of viscoelastic fluids at high Weissenberg numbers (Wi). Specifically: 1. **Abrupt changes in physical properties and rheological behavior**: The physical properties (such as density, viscosity, etc.) and rheological behavior (such as elastic modulus, etc.) of two-phase fluids undergo abrupt changes at the interface, which brings great difficulty to numerical simulation. 2. **High Weissenberg Number Problem (HWNP)**: When the Weissenberg number is high, the elastic effects of viscoelastic fluids increase significantly, leading to instability in numerical simulations. The Weissenberg number is a dimensionless parameter that describes the strength of fluid elasticity; the higher its value, the stronger the elastic effects of the fluid. 3. **Interface deformation and breakup**: In shear flow, droplets may undergo severe deformation or even breakup, which poses higher requirements for the stability and computational efficiency of numerical simulations. To address these challenges, the authors propose an unstructured finite volume method based on the geometric Volume of Fluid (VOF) approach, specifically designed to handle two-phase fluids containing viscoelastic phases. This method integrates the DeboRheo framework, allowing flexible combination of different rheological models and employing appropriate stabilization techniques to solve the high Weissenberg number problem. Additionally, this method can use unstructured grids in complex geometric domains, improving the overall accuracy of the solution. Through this method, the authors aim to: - Establish an extended geometric Volume of Fluid method suitable for simulating viscoelastic fluids. - Validate the numerical stability and robustness of this method at high Weissenberg numbers. - Investigate the deformation behavior of a single droplet in shear flow, including 2D and 3D cases, covering scenarios of viscoelastic droplets in a Newtonian fluid matrix and Newtonian droplets in a viscoelastic matrix.

An unstructured geometrical un-split VOF method for viscoelastic two-phase flows

Development of a Volume of Fluid Method for Computing Interfacial Incompressible Fluid Flows

A Volume of Fluid Method Algorithm for Simulation of Surface Tension Dominant Two-Phase Flows

A numerical stabilization framework for viscoelastic fluid flow using the finite volume method on general unstructured meshes

voFoam - A geometrical Volume of Fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM

High-Precision Reconstruction Of Gas-Liquid Interface In Plic-Vof Framework On Unstructured Mesh

Inconsistencies in Unstructured Geometric Volume-of-Fluid Methods for Two-Phase Flows with High Density Ratios

Numerical simulation of non-isothermal viscoelastic flows at high Weissenberg numbers using a finite volume method on general unstructured meshes

Unstructured un-split geometrical Volume-of-Fluid methods -- A review

A coupled VOF/embedded boundary method to model two-phase flows on arbitrary solid surfaces

Numerical simulation of 3D-unsteady viscoelastic free surface flows by improved smoothed particle hydrodynamics method

Application of Bounded and Compressed VOF Method to Interfacial Flow

Numerical Simulation of Viscoelastic Fluid Flows Using a Robust FVM Framework on Triangular Grid

Numerical Simulation of the Non-isothermal Viscoelastic Flow Past a Confined Cylinder

A Lattice Boltzmann Method and Asynchronous Model Coupling for Viscoelastic Fluids

A Three Dimensional Coupled VOF and Level Set (VOSET) Method with and Without Phase Change on General Curvilinear Grids

A Volume of Fluid Method for Three Dimensional Direct Numerical Simulations of Immiscible Droplet Collisions

A Geometric VOF Method for Interface Flow Simulations

A Coupled Volume-of-fluid and Level Set Method Based on Multi-Dimensional Advection for Unstructured Triangular Meshes

Computation of Deformable Interface Two‐Phase Flows: A Semi‐Lagrangian Finite Element Approach