https://doi.org/10.1140/epjs/s11734-025-01758-x
Regular Article
A residual-based non-orthogonality correction for force-balanced unstructured volume-of-fluid methods
1
Mathematical Modeling and Analysis, Technische Universität Darmstadt, Darmstadt, Germany
2
Robert Bosch GmbH, Corperate Research, Robert-Bosch-Campus 1, 71272, Renningen, Germany
Received:
12
March
2025
Accepted:
22
June
2025
Published online:
4
August
2025
The collocated unstructured finite volume method is widely used to simulate flows in geometrically complex flow systems. It supports arbitrary finite volume shapes, significantly simplifying the generation of boundary-fitted domain discretizations (meshes). However, fitting the mesh to the domain boundary introduces the so-called non-orthogonality discretization error, severely impacting the stability and accuracy of two-phase flow simulations where the fluid interface meets the domain boundary. We show that the non-orthogonality error introduces a significant imbalance of forces at fluid interfaces when the one-field formulation of the Navier–Stokes equations is discretized by the collocated unstructured volume-of-fluid method on boundary-fitted non-orthogonal meshes. This force imbalance error leads to large errors when the flow is driven by the surface tension or gravitational force, manifesting as artificial (parasitic) velocities. We deterministically recover the force balance using a residual-based control of the non-orthogonality defect correction. Our approach for applying the non-orthogonality correction offers three key advantages for simulations of wetting processes on geometrically complex surfaces. First, contrary to conventional non-orthogonal correction methods, in our approach, the number of non-orthogonal corrections is no longer a user-defined parameter, to be determined by trial and error. Second, our approach significantly increases accuracy and simultaneously speeds up the solution algorithm for the one-field formulation of the Navier–Stokes equations for two-phase flows. Finally, our approach can be applied with minimal modifications to existing methods. The canonical stationary droplet and stationary water column verification cases and a complex droplet wetting a spherical obstacle verify and validate the proposed method.
Research was conducted at Technische Universität Darmstadt, Darmstadt, Germany.
© The Author(s) 2025
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
