New Property Averaging Scheme for Volume of Fluid Method for Two-Phase Flows with Large Viscosity Ratios

TL;DR

This paper introduces a new property averaging scheme for the Volume-of-Fluid method, significantly improving the accuracy of two-phase flow simulations involving high-viscosity liquids, validated through multiple experimental comparisons.

Contribution

A novel property averaging technique for the VoF method that enhances accuracy in high-viscosity two-phase flow simulations, validated with experimental data.

Findings

Improved accuracy for high-viscosity liquids in VoF simulations.

Significant performance enhancement over the original VoF method.

Validated across diverse two-phase flow problems.

Abstract

To predict liquid-gas two-phase flow phenomena, accurate tracking and prediction of the evolving liquid-gas interface is required. Volume-of-Fluid or VoF method has been used in the literature for computationally modeling of such flows. In the VoF method, a single set of governing equations are solved for both phases along with an advection equation for the volume fraction. The properties in each computational cell are determined by a linear weighted average of the properties of the two fluids based on the phase fraction. While the method predicts water-air flows well, the predictions tend to deviate significantly from experimental data for liquids with high viscosity. A new property averaging technique is proposed in this paper, which is shown to provide accurate results for high viscosity liquids. Computational predictions using the open source VoF solver interFoam (available as a part of the OpenFOAM computational tool), and those obtained using the proposed method are compared with experimental data for multiple two-phase applications. Four different problems, viz., suspended droplet in air, jet breakup, drop impact on thin films, and air entrapment during drop interaction with liquid pool, are considered to extensively validate the new method. Experimental data are used to cover a range of surface tension and viscosities. For all cases, the modified VoF solver is observed to perform significantly better than original VoF method.