No two particles may occupy the same point in space. No two fluids may occupy the same cell in a mesh. . . unless we use multiphase modeling. In computational fluid dynamics (CFD), we often need to model scenarios that involve more than one fluid. The volume of fluid method (VOF) allows the CFD solver to combine and mix multiple fluids together on a single mesh. This article introduces the basic concepts of the VOF method and its implications for CFD modeling.

VOF assumes that every cell in the domain contains a fraction of each fluid. (Fractions of 0 or 1 regularly occur.) The solver creates a composite fluid for each cell, based on the fraction of each fluid within that cell. VOF developed its own nomenclature to keep the terminology straight.

Individual material fluids are called phases.

The entire composite of all materials is referred to as the generic fluid.

What determines the portion of each phase in every cell? The VOF method introduces a new transport equation called the volume fraction (VF). If the model requires more than two fluid phases, there are multiple volume fractions. Multiphase modeling always creates one less transport equation than the total number of fluid phases modeled.

- 2 fluid phases >> 1 VF equation
- 3 fluid phases >> 2 VF equations
- 4 fluid phases >> 3 VF equations
- …

These transport equations capture the change in each of the fluid phases and couple that change with the momentum and continuity equations.

The volume fraction tracks the movement of the phases, while the remaining equations capture the forces on the fluid. But there are two options on the exact method of combining volume fraction with momentum.

Normally, the continuity equation always combines the fluid phases and works with a single composite fluid. One continuity equation for the entire composite fluid. But the modifications to the momentum equation change depending on the situation. VOF generally follows two methods regarding the momentum equation:

- Macro scale methods
- Sub-grid scale methods

Macro scale methods occur in cases where the interface between the fluids stretches across multiple cells. The interface can be defined by tracking the volume fraction across cells in the mesh. Common examples are free surface flow in ships, large mixing vessels in industrial flow (Figure 3‑1), river flow in dams, channel flow in pipes. For all these cases, the interface was large and stretched across multiple cells. Macro scale methods apply to the majority of VOF simulations.

For the macro scale methods, the momentum equations for all the phases get combined. The entire simulation employs only a single set of momentum equations for the composite fluid. This works best for cases where fluid mixing is not a concern. Situations where we are not concerned with the differences in momentum between the two different phases.

The sub-grid method works on the tiny end. It models scenarios where the interface between fluid phases is smaller than a single cell. Or at least too small to be clearly defined in 1-2 cells. Imagine tiny gas bubbles dissolving into a liquid, or one liquid phase diffusing into another. Those would be appropriate for sub-grid methods. You see these methods most often in industrial processes, involving mixing of fluid phases. Particularly gases and liquids where the industrial process intentionally creates miniscule gas bubbles.

In sub-grid methods, the solver creates a separate set of momentum equations for each fluid phase. This requires more computing time, since you just increased the number of momentum equations to solve. As part of the VOF method, additional terms add into these momentum equations to capture the momentum transfer between the different fluid phases.

For that momentum transfer, this method assumes each fluid phase is a droplet too small to resolve in the grid. As the CFD engineer, you need to specify the details about the fluid interface so the solver can calculate the moment transfer:

- Lift coefficient
- Drag coefficient
- Interface area / bubble diameter

The exact extent of the details depends on your specific solver. The CFD engineer needs to have an educated guess at these specifics, which is why this method typically only applies to industrial mixing processes. You need an application where the details can be tightly controlled.

Volume of fluid modeling (VOF) expands the capabilities of CFD to allow limitless combinations of different fluid phases. The volume fraction allows us to combine different phases into a single composite fluid. With VOF, fluids mix and dance across our domain. The world of VOF encompasses everything from droplets of diesel spraying in cylinder all the way up to tsunami waves crashing against the city of Tokyo. The advanced capabilities of VOF inspire awe.

[1] | S. Tao, F. Yuqing, L. Graeme and J. Kaixi, “CFD simulation of bubble recirculation regimes in an internal loop airlift reactor,” in 27th International Mineral Processing Congress, <https://www.researchgate.net/publication/289001212_CFD_simulation_of_bubble_recirculation_regimes_in_an_internal_loop_airlift_reactor>, Accessed: 2019, Jan, 01., January 2014. |

[2] | V. R. Raj, “Quadratic Profile Used in QUICK Scheme,” Wikimedia Commons, 12 Nov 2012. . Available: https://commons.wikimedia.org/wiki/File:Quadratic_profile.jpg. . |

[3] | Max Pixel, “Cumulus Storm Turbulence Thunderstorm Cloud Roller,” Max Pixel, 01 Jan 2019. . Available: https://www.maxpixel.net/Cumulus-Storm-Turbulence-Thunderstorm-Cloud-Roller-567678. . |

[4] | S. Wasserman, “Choosing the Right Turbulence Model for Your CFD Simulation,” Engineering.com, 22 Nov 2016. . Available: https://www.engineering.com/DesignSoftware/DesignSoftwareArticles/ArticleID/13743/Choosing-the-Right-Turbulence-Model-for-Your-CFD-Simulation.aspx. . |

[5] | Q. Wang, C. Yan and T. Hui, “Mechanism Design for Aircraft Morphing Wing,” Research Gate, <https://www.researchgate.net/figure/y-plus-value-for-the-CFD-model-10-degree-of-extension-angle-of-attack-6_fig4_268478784>, Accessed: 2019, Jan, 01, April 2012. |