PCIC: Cylindrical Volume Moment Calculation and Interface Reconstruction for Sub-Grid Modeling in Volume of Fluid Methods

TL;DR

This paper introduces PCIC, a novel cylindrical interface reconstruction method for sub-grid ligament modeling in VOF simulations, improving accuracy and efficiency in capturing small-scale flow features.

Contribution

The paper presents a new cylindrical reconstruction strategy, PCIC, enabling exact mass conservation and automatic ligament detection for enhanced sub-grid modeling in VOF methods.

Findings

PCIC achieves robust sub-grid ligament reconstruction.

Low-resolution PCIC simulations match high-resolution results.

Automatic detection allows dynamic switching between PLIC and PCIC.

Abstract

The accurate modeling of topology changes remains a significant challenge in geometric Volume of Fluid (VOF) simulations. When using traditional single-plane reconstruction (PLIC), fluid structures smaller than the mesh size cannot be resolved and spurious numerical breakup is triggered; this impacts important flow statistics such as drop size distributions. Recent advances have introduced paraboloid and two-plane reconstructions, which have improved high-curvature performance and enabled sub-grid film reconstructions, respectively. However, sub-grid ligament reconstructions have remained elusive. In this work, a novel cylindrical interface reconstruction strategy called PCIC is introduced for sub-grid ligament modeling. PCIC is facilitated by deriving the analytical volume moments of quadratic cylinders clipping polyhedra; this allows for exact mass conservation during volume moment transport. From the transported moments, a straight circular cylinder can be defined in the center cell of a 5x5x5 stencil. First, a quadratic principal curve is fitted to the normalized first-order moments in the stencil (the liquid barycenters), from which the cylinder's orientation and origin are approximated. The cylinder radius is then chosen to conserve volume. On-the-fly ligament detection is achieved using connected-component labeling and moments of inertia criteria, which allows for simulations to automatically choose between PLIC and PCIC in each interface cell at runtime. PCIC is demonstrated in multiphase flow test cases, where it exhibits robust reconstruction of sub-grid ligaments. This allows for relatively low-resolution PCIC simulations to provide comparable results to traditional high-resolution simulations.