Axisymmetric VOF Solver

A Navier-Stokes solver for interfacial two-phase flows under axisymmetric conditions.

Dependencies

• C Compiler
• GNU Make
• OpenMP (optional)

Quick Start

Check out the readme of the other library.

Methodology

We consider two incompressible and immiscible Newtonian liquids (e.g., water and air) which are separated by free and deformable surface. The liquids are confined inside polar coordinates under axisymmetric condition. Their dynamics are governed by the incompressibility constraint, the mass conservation, and the momentum balance. The whole equations are solved on the general orthogonal coordinate system (see the documentation of the other library).

The following numerical techniques are adopted to solve the equations.

• Second-order accurate central-finite-difference scheme for spatial discretization
• Predictor-corrector method (SMAC / fractional-step methods) for enforcing incompressibility constraint
• Volume-of-fluid (THINC) method for interface capturing
• Continuum surface-tension force model for describing surface-tension force
• Approximating scheme of the variable-coefficient Poisson equation proposed by Dodd and Ferrante to facilitate the fast Poisson solver
• Energy-consistent discretization for stable integration

The governing equations are briefly documented here. For more details, check-out the following pages.

• Simple TC Solver
• Simple Bubbly Flow Solver