Evaluation and Improvements to Interfacial Curvature Predictions in interFoam

Abstract

Improvements to the interfacial curvature of interFoam based on (i) the smoothing of the liquid fraction field and (ii) the creation of a signed distance function (ϕ-based) are implemented. While previous work in this area has focused on evaluating spurious currents and similar configurations, the tests implemented in this work are more applicable to sprays and hydrodynamic breakup problems. For the ϕ-based method, a dual approach is developed based on a geometric reconstruction of the interface at interfacial cells and the solution of the Hamilton-Jacobi equation away from these cells. The more promising results are from this method, where the lack of convergence of Laplace pressure predictions existing in the standard version of interFoam is fixed, resulting in second-order convergence. Similar but less drastic improvements are observed for other exercises consisting of the oscillation of a droplet, a 2-phase Orr–Sommerfeld problem, the Rayleigh–Plateau instability, and the retraction of a liquid column. It is only when the dynamics are either entirely governed by surface tension or are heavily influenced by it that we see the need to substitute the standard interFoam curvature approach with a more accurate scheme. For more realistic problems, which naturally include more complicated dynamics, the difference between the standard approach and the ϕ-based approach is minimal.

Breakup of planar liquid sheets injected at high speed in a quiescent gas environment

“…Capturing the full nonlinear dynamics of the perturbation process leading to the sheet breakup is done in the present investigation using a VoF methodology (Tryggvason, Scardovelli & Zaleski 2011). The particular flavour of the VoF method employed in the present study consists of an algebraic VoF solver phibasedFoam (Agarwal, Ananth & Trujillo 2022). This solver is a slightly modified version of the interFoam code, which is part of the OpenFoam open source distribution of continuum mechanics solvers (OpenFoam foundation, https://openfoam.org/).…”