A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries

This paper presents a novel, fully geometric, and conservative 3D Volume-of-Fluid method that combines a modified advection scheme for mixed cells with a height-function-based contact angle imposition to accurately simulate moving contact lines with hysteresis on complex embedded boundaries while overcoming severe time-step limitations.

The Gist

Imagine you are trying to simulate how a drop of water behaves when it lands on a surface. If that surface is a perfectly flat, smooth table, it's relatively easy to predict where the water goes. But what if the surface is a bumpy rock, a cone, a wire mesh, or a plate with a hole in it? And what if the water doesn't just slide; it gets "stuck" in some spots due to friction or chemical differences (a phenomenon called hysteresis)?

This paper presents a new, super-accurate computer program designed to solve exactly that problem in 3D. The authors, Chong-Sen Huang and colleagues, have built a digital toolkit that can track the edge of a water droplet (the "contact line") as it moves over incredibly complex, bumpy, or irregular shapes without losing any "water" in the simulation.

Here is a breakdown of their breakthrough using simple analogies:

1. The Problem: The "Tiny Puddle" Dilemma

In computer simulations, the world is divided into a grid of tiny 3D boxes (like a giant Rubik's cube).

• The Issue: When a solid object (like a rock) cuts through these boxes, it creates "mixed cells"—boxes that are partly rock and partly water.
• The Trap: If the rock cuts the box at a very sharp angle, the water part becomes a tiny, sliver-like puddle. Traditional computer methods get confused by these tiny puddles. They either lose the water (mass conservation error) or force the computer to take incredibly tiny steps in time to avoid crashing, making the simulation painfully slow.
• The Analogy: Imagine trying to pour water from a bucket into a cup, but the cup has a weirdly shaped hole in the bottom. If you pour too fast, the water spills out the sides (error). If you pour too slowly to be safe, it takes forever. The authors found a way to pour at a normal speed without spilling a single drop.

2. Solution A: The "Smart Redistribution" (Fixing the Flow)

The authors developed a new way to move the water through the grid.

• The Old Way: It was like trying to push a heavy box through a narrow, jagged doorway. If the box got stuck, the whole process halted.
• The New Way: They introduced a Redistribution Strategy. Think of it like a team of people passing a bucket of water. If one person (a cell) is too small to hold the water they are supposed to receive, they don't stop the line. Instead, they immediately pass the "overflow" to their neighbors.
• The Result: This removes the "time-step" bottleneck. The simulation can run fast and smooth, even with jagged rocks, because the water is never lost; it's just temporarily shared with neighbors until it finds a permanent home.

3. Solution B: The "Parabolic Crystal Ball" (Fixing the Angle)

When a drop of water hits a surface, it forms a specific angle (the contact angle). On a flat table, this is easy to calculate. But on a curved or bumpy surface, the angle changes constantly.

• The Old Way: Previous methods tried to guess the shape of the water surface by drawing a straight line (linear fitting) through the data points.
• The Analogy: Imagine trying to guess the curve of a rainbow by drawing a straight stick through it. It works okay if the rainbow is flat, but if it's a big arc, the stick misses the mark.
• The New Way: The authors use a Paraboloid Fitting method. Instead of a straight stick, they use a flexible, curved template (like a parabola or a bowl shape) that fits perfectly over the data points.
• The "Pre-Fitting" Trick: Before they fit the final curve, they do a "pre-fit" to check which data points are reliable. It's like a detective checking alibis before building a case. This ensures that even on a jagged rock, the computer knows exactly where the water edge is and what angle it should make.

4. Solution C: The "Memory Effect" (Hysteresis)

In the real world, water doesn't always slide smoothly. Sometimes it gets "stuck" (pinned) on a rough spot until the force pushing it is strong enough to break it free. This is called hysteresis.

• The Innovation: The new method includes a "memory" for the droplet. It remembers the "receding angle" (how much it wants to shrink) and the "advancing angle" (how much it wants to grow).
• The Result: The simulation can now show a droplet sliding down a bumpy hill, getting stuck on a bump, and then suddenly snapping forward when the slope gets steep enough. This matches real-life physics much better than older models.

Why Does This Matter?

This isn't just about water droplets. This technology is crucial for:

• Inkjet Printing: Ensuring ink lands exactly where it's supposed to on curved screens or uneven paper.
• Microfluidics: Designing tiny lab-on-a-chip devices where fluids move through microscopic, winding channels.
• Oil Recovery: Understanding how oil moves through the complex, rocky pores of the earth.
• Self-Cleaning Surfaces: Designing surfaces that shed water efficiently, even if they are bumpy.

The Bottom Line

The authors have built a 3D digital microscope for fluids. It is "conservative" (it never loses a drop of water), "sharp" (it sees the edge of the water clearly, not blurry), and "robust" (it doesn't crash when the surface gets weird). By combining a smart way to move water through tiny gaps with a curved, intelligent way to calculate angles, they have solved a decades-old headache in computer physics, allowing us to simulate wetting on any shape imaginable.

Technical Summary

Here is a detailed technical summary of the paper "A 3D sharp and conservative VOF method for modeling the contact line dynamics with hysteresis on complex boundaries".

1. Problem Statement

The paper addresses the numerical challenges associated with simulating moving contact line dynamics (the intersection of liquid, gas, and solid phases) on arbitrarily complex 3D solid boundaries. While significant progress has been made in simulating these phenomena on regular (flat, aligned) boundaries, extending these methods to complex, embedded geometries (e.g., curved surfaces, porous media, or irregular substrates) introduces two fundamental difficulties:

1. Mass Conservation in Mixed Cells: Standard geometric Volume-of-Fluid (VOF) advection schemes fail to maintain strict local mass conservation in "mixed cells" (cells intersected by both the solid boundary and the liquid/gas interface). This leads to non-physical fluid accumulation or depletion, particularly when small cut cells are present.
2. Accurate Contact Angle Imposition: Enforcing the correct contact angle condition on curved or irregular solid surfaces is difficult. Existing methods often rely on linear extrapolation or simple ghost-cell approaches, which introduce significant errors in curvature estimation and contact line motion, especially for extreme contact angles ($\theta < 45^\circ$ or $\theta > 135^\circ$) or when the solid boundary is not aligned with the grid.

2. Methodology

The authors propose a unified numerical framework implemented within the open-source Basilisk solver, combining a Geometric Volume-of-Fluid (VOF) method with an Embedded Boundary Method (EBM).

A. Conservative 3D Interface Advection

To handle the advection of the liquid/gas interface in mixed cells:

• Flood-Fill Reconstruction: A 3D flood-fill algorithm is developed to reconstruct the piecewise linear interface (PLIC) within irregular convex polyhedra formed by the intersection of the solid boundary and the computational cell. This involves rotating the polyhedron, subdividing it into slabs, and solving for the interface height iteratively.
• Flux Correction: A modified advection scheme is introduced to correct flux imbalances caused by the solid boundary. The effective advection width is adjusted ($u_n^*$) to compensate for the volume blocked by the solid, ensuring strict local volume conservation.
• Redistribution Strategy: To overcome the severe time-step constraints imposed by small cut cells (where the fluid volume fraction is tiny), a flux redistribution scheme is introduced. Instead of removing small cells (which distorts geometry in 3D) or merging them (which complicates geometry), excess or deficit fluid volume is redistributed to downstream neighbors. This completely removes the CFL constraint associated with small cut cells while preserving global and local volume conservation.

B. 3D Contact Angle Imposition (Paraboloid Fitting)

To accurately enforce contact angles on complex surfaces:

• Height Function (HF) Extension: The standard HF method is extended to handle mixed cells. A pre-fitting paraboloid is constructed using valid HF points from the fluid region to approximate the local interface profile.
• Constrained Paraboloid Fitting: A final paraboloid is fitted using a nonlinear system that includes:
  1. HF values from the fluid region.
  2. Four auxiliary points derived from the pre-fitted surface.
  3. The Contact Angle Constraint: The system explicitly enforces that the angle between the reconstructed interface normal and the solid normal matches the prescribed contact angle $\theta$.
• Vertical Height Function (VHF): For degenerate cases where the solid surface is nearly parallel to the grid plane (making standard "horizontal" HFs invalid), a Vertical HF method is introduced. This constructs the interface based on the contact angle-dependent normal vector, allowing accurate curvature estimation even for steep or nearly horizontal surfaces.
• Hysteresis Modeling: A bisection method is employed within a hysteresis window ($[\theta_{rec}, \theta_{adv}]$) to determine the dynamic contact angle. The contact line is allowed to move only when the local angle exceeds the advancing or receding limits, ensuring physically realistic pinning/depinning behavior.

3. Key Contributions

1. First Fully Geometric & Conservative 3D VOF Scheme for Complex Boundaries: The paper presents the first method capable of strictly conserving mass in 3D mixed cells containing liquid, gas, and solid phases, eliminating the fluid accumulation issues seen in previous sharp-interface approaches.
2. Time-Step Unrestricted Advection: The introduction of the redistribution strategy allows simulations to proceed without the restrictive time-step limitations typically caused by small cut cells in complex 3D geometries, significantly improving computational efficiency.
3. Robust 3D Contact Angle Algorithm: The novel paraboloid-fitting technique, combined with the Vertical HF method, provides high-fidelity contact angle enforcement on arbitrary 3D surfaces. It overcomes the errors of linear fitting methods, particularly for extreme contact angles and curved substrates.
4. Integration of Hysteresis: The framework successfully incorporates contact angle hysteresis into the 3D geometric VOF framework, enabling the simulation of realistic wetting dynamics where the contact line pins and depins based on local geometry and flow history.

4. Results and Validation

The method was validated through a series of challenging benchmark tests:

• Volume Conservation: In a rotating concentric sphere model, the proposed scheme achieved machine-precision local volume conservation ($L_\infty \approx 10^{-15}$), whereas traditional schemes showed significant errors ($L_\infty \approx 10^{-2}$).
• Droplet Spreading on Flat/Inclined Plates: The paraboloid fitting method produced equilibrium contact radii matching theoretical predictions with high accuracy, regardless of the solid surface's embedding depth or orientation. It maintained perfect circular symmetry, unlike linear fitting methods which produced distorted shapes.
• Droplet Spreading on Spherical Surfaces: The method accurately captured droplet equilibrium shapes on curved substrates with relative radius errors $< 0.6\%$ and volume loss $< 10^{-6}$.
• Shear-Driven Motion with Hysteresis: Simulations of droplets sliding under shear flow on flat and sine-wavy substrates successfully reproduced pinning/depinning dynamics and steady-state morphologies consistent with experimental data and previous literature.
• Complex Geometries: The method successfully simulated:
  • Droplet impact on plates with round and sharp orifices (capturing pinning and partial entrapment).
  • Spontaneous climbing of droplets on conical surfaces driven by curvature gradients.
  • Droplet spreading on mesh-shaped surfaces with high local curvature and abrupt topological changes.

5. Significance

This work represents a major advancement in computational multiphase fluid dynamics. By resolving the dual challenges of strict mass conservation and accurate contact angle enforcement on complex 3D geometries, the proposed method enables high-fidelity simulations of wetting phenomena in realistic engineering applications. These include inkjet printing, microfluidics, water electrolysis, and coating processes where fluid interfaces interact with intricate solid structures. The removal of the small-cut-cell time-step constraint makes these simulations computationally feasible for long-duration or large-scale problems, bridging the gap between theoretical sharp-interface models and practical engineering needs.