Volume of Fluid Method
Also known as: VOF, VoF, volume fraction method
The Volume of Fluid (VOF) method is an Eulerian interface-tracking technique for multiphase flow simulations, developed by Hirt and Nichols in 1981. Instead of explicitly tracking the interface between phases, VOF advects a scalar field (the volume fraction) that represents the fractional volume of one phase in each grid cell. This approach elegantly handles topological changes (splashing, merging, breaking) without explicit interface reconstruction, making it ideal for complex multiphase flows on fixed Eulerian grids.
Read the full method
Sign in with a free account to read this section.
Method map
The neighbourhood of related methods — select a node to explore.
When to use it
VOF is ideal for multiphase flows involving free surfaces, bubbles, droplets, and complex mixing on structured or unstructured grids. Use VOF when topological changes (breaking, coalescence) are significant and explicit interface tracking is impractical. VOF works well with standard Navier-Stokes solvers on existing Eulerian grids. Prefer VOF for industrially relevant problems; it is mature, robust, and available in commercial CFD packages. Avoid VOF if extreme interface sharpness is critical (use Level Set instead) or if the interface geometry is extremely complex and cannot be accurately reconstructed.
Strengths & limitations
- Handles topological changes naturally without explicit interface tracking or merging algorithms
- Conservative by construction: volume fractions integrate to conserved phase volumes
- Computationally efficient on structured grids; integrates readily with RANS simulations
- Works with arbitrary numbers of phases; generalizes to multiphase flows naturally
- Well-established industrial tool with mature implementations in ANSYS, OpenFOAM, Star-CCM+
- Interface is smeared over several grid cells; accuracy depends on interface reconstruction scheme
- Smearing increases with time due to advection diffusion; longer simulations lose interface sharpness
- Interface reconstruction (PLIC, ELIAS) adds complexity and computational cost
- Surface tension calculation from volume fractions can be noisy; requires smoothing
- Handling of wetting and contact angle is model-dependent and often inadequate
Frequently asked
How does VOF differ from Level Set method?
VOF advects a volume fraction (conserved, bounded 0-1), resulting in conservative phase volumes but smeared interfaces. Level Set tracks a signed distance function (not conserved), yielding sharper interfaces but requiring mass correction. VOF is more robust; Level Set is more accurate. Hybrid methods (coupled VOF-LS) combine benefits of both.
What is interface reconstruction and why is it necessary?
Volume fractions define phase amounts per cell but not interface position. Interface reconstruction (PLIC, ELIAS) geometrically estimates the interface location from F values, allowing accurate curvature computation and surface tension application. Without it, interfaces become thick and diffuse over time. PLIC assumes the interface is a straight line (2D) or plane (3D) per cell.
How do I reduce interface diffusion in VOF simulations?
Use high-order, compressive discretization schemes for volume fraction advection (CICSAM, MULES). These schemes limit numerical diffusion while avoiding oscillations. Ensure sufficient grid resolution at the interface: typically 3-5 cells across the interface. For long simulations, periodically apply interface sharpening to restore crisp interfaces.
How is surface tension computed in VOF?
Curvature κ is estimated from the interface geometry (reconstructed from volume fractions) using finite differences. The surface tension force f_s = σ κ n̂δ (where σ is surface tension, n̂ is normal, δ is interface thickness) is applied to the momentum equation. Noisy curvature estimates are smoothed to prevent spurious oscillations.
Can VOF simulate contact angles and wetting?
Yes, by specifying contact angle boundary conditions at walls. The interface-wall intersection angle is enforced through wall-adjacent volume fractions or explicit contact angle models. However, accuracy is limited; thin films and capillary-dominated wetting are better handled with slip models or precursor films.
Sources
- Hirt, C. W., & Nichols, B. D. (1981). Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of Computational Physics, 39(1), 201-225. DOI: 10.1016/0021-9991(81)90145-5 ↗
- Youngs, D. L. (1982). Time-dependent multi-material flow with large fluid distortion. In Numerical Methods for Fluid Dynamics (pp. 273-285). Academic Press. link ↗
- Ubbink, O., & Issa, R. I. (1999). A method for capturing sharp fluid interfaces on arbitrary meshes. Journal of Computational Physics, 153(1), 26-50. DOI: 10.1006/jcph.1999.6276 ↗
How to cite this page
ScholarGate. (2026, June 3). Volume of Fluid Method. ScholarGate. https://scholargate.app/en/fluid-dynamics/volume-of-fluid
Which method?
Set this method beside its closest kin and read them side by side — the library lays the books on the table; the choice is yours.
- Eulerian-Lagrangian ModelFluid Dynamics↔ compare
- Large Eddy SimulationFluid Dynamics↔ compare
- Lattice Boltzmann MethodFluid Dynamics↔ compare
- Level Set MethodFluid Dynamics↔ compare
- Smoothed Particle HydrodynamicsFluid Dynamics↔ compare