Free Surface Flow
Free Surface Flow: Theoretical Foundations
Overview
Professor, what is free surface flow? Is it different from the VOF method?
Free surface flow refers to the general category of flows where an interface between liquid and gas (the free surface) exists. This includes dam breaks, sloshing, ship wave-making, filling processes, and surface waves. The VOF method is a representative technique for solving it, but there are also other solution methods such as the Level Set method, SPH (Smoothed Particle Hydrodynamics) method, and Phase-Field method.
What physics is unique to free surface flow?
Surface tension and contact angle become important at small scales. At large scales, the Froude number $Fr$, which is the ratio of gravity to inertia, is the governing parameter.
CSF Model for Surface Tension
How is surface tension modeled?
The CSF (Continuum Surface Force) model by Brackbill et al. (1992) is the most widely used. It calculates it as a volume force from the interface curvature $\kappa$.
I've heard that curvature calculation accuracy is problematic.
The biggest challenge of CSF is parasitic currents (spurious currents). Unphysical velocity fields arise from discretization errors at the interface. This is particularly severe in flows with low capillary number $Ca = \mu U / \sigma$ (surface tension dominated).
Comparison of Interface Tracking Methods
| Method | Interface Representation | Mass Conservation | Interface Sharpness | Computational Cost |
|---|---|---|---|---|
| VOF Method | Volume fraction $\alpha$ | Exact | Scheme dependent | Low ~ Medium |
| Level Set Method | Signed distance function $\phi$ | Non-conservative (requires reinitialization) | Sharp | Medium |
| CLSVOF | VOF + Level Set | Good | Sharp | Medium ~ High |
| Phase-Field | Order parameter | Conservative | Diffusive | High |
| SPH | Particles | Conservative | Particle resolution dependent | High |
Which method should I choose?
In industry, the VOF method is the most widely used due to its balance of cost and accuracy. For cases where interface curvature accuracy is critical (e.g., microfluidics), CLSVOF or Phase-Field methods are advantageous. For flows involving large deformation and splashing, the SPH method is also an option.
The Mechanics of Free SurfacesโWhere Bernoulli Meets Capillarity
The free surface (gas-liquid interface) is a place where the richest physical phenomena in fluid mechanics are condensed. Gravity, surface tension, viscosity, and inertia all compete, and the physics governing them is determined by the combination of dimensionless numbers (Bond number, Weber number, Capillary number, Froude number) that represent their ratios. When the Bond number Bo = ฯgL^2/ฯ exceeds 1, gravity overwhelms surface tension; conversely, when it is small, capillary phenomena dominate. In microfluidic devices, Bo << 1, and the free surface converges to a spherical shape, while for large ship waves, Bo >> 1, and surface tension can be ignored. This scale diversity spanning over 20 orders of magnitude in Bo number makes free surface CFD a technology required across all industries.
Computational Methods for Free Surface Flow
Details of Numerical Methods
Please tell me the key numerical points for free surface flow.
Interface Courant number management and mesh design are most important. In the VOF method, the interface Courant number $Co_\alpha < 0.25$ (max 0.5) must be maintained.
Open Channel Flow
For river and channel flows, Fluent and CFX have dedicated Open Channel Flow settings. Water level can be directly specified at inlets/outlets, making water surface initialization easy.
What are the key points for setting up Open Channel Flow?
- Inlet: Specify water level (Free Surface Level) and flow rate
- Outlet: Specify water level (Pressure outlet + Open Channel)
- Initialization: Set still water surface with Flat Initialization
- Numerical Scheme: Explicit VOF + Geo-Reconstruct (high accuracy)
Wave Generation and Absorption
For ocean engineering wave simulations, wave-making and wave-absorbing boundary conditions are necessary.
| Method | Overview | Tool Support |
|---|---|---|
| Stokes wave theory | Prescribe wave velocity/water level at inlet | Fluent, STAR-CCM+ |
| 5th-order Stokes wave | Higher-order nonlinear wave | STAR-CCM+ |
| Wave damping zone | Absorb waves with damping | All tools (UDF/Field Function) |
| Numerical wave tank | Wave-maker + damping zone | OpenFOAM (waves2Foam) |
OpenFOAM's waves2Foam is famous, right?
It's a library developed by Jacobsen et al. (DTU, 2012), implementing boundary conditions and relaxation zones for wave generation/absorption. It is widely used for wave force analysis on marine structures.
Parasitic Current Countermeasures
How can parasitic currents be suppressed?
| Countermeasure | Effectiveness | Notes |
|---|---|---|
| Mesh refinement | High | Make cells near the interface sufficiently small |
| Height Function method | Very High | For structured grids, improves curvature accuracy |
| Sharp Surface Force | High | Available in STAR-CCM+ |
| CLSVOF | High | Improves curvature calculation with Level Set |
| Adjusting cAlpha | Medium | Interface compression parameter in OpenFOAM |
The Fusion of Level Set and VOF MethodsโThe Advent of CLSVOF
The Level Set method can accurately calculate geometric interface information (normal, curvature) but has the weakness of incomplete mass conservation. VOF excels at mass conservation but has lower accuracy in describing interface shape. The Coupled Level-Set and VOF (CLSVOF) method proposed by Sussman et al. in the 2000s combines the strengths of both, a hybrid method that calculates curvature using the Level Set function while conserving mass with VOF. In its application to sloshing analysis of aircraft fuel tanks, benchmark results showed that CLSVOF achieved equivalent accuracy with 25% fewer cells than simple VOF.
Related Topics
Experience the theory firsthand with the interactive simulator for this field
All Simulators