Part 7:  Multiphase Modeling


1.0 Introduction

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.

2.0 Volume of Fluid Method

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. 

3.0 Calculating Volume Fraction

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:

  1. Macro scale methods
  2. Sub-grid scale methods

3.1 Macro 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.

Volume Fraction of Gas Reactor

Figure 3-1: Volume Fraction of Gas Reactor [1]


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.

3.2 Sub-Grid Scale Methods

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 minuscule 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:

  1. Lift coefficient
  2. Drag coefficient
  3. 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.

4.0 Conclusion

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.

5.0 References


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, <>, Accessed: 2019, Jan, 01., January 2014.


V. R. Raj, “Quadratic Profile Used in QUICK Scheme,” Wikimedia Commons, 12 Nov 2012. [Online]. Available: [Accessed 01 Jan 2019].


Max Pixel, “Cumulus Storm Turbulence Thunderstorm Cloud Roller,” Max Pixel, 01 Jan 2019. [Online]. Available: [Accessed 01 Jan 2019].


S. Wasserman, “Choosing the Right Turbulence Model for Your CFD Simulation,”, 22 Nov 2016. [Online]. Available: [Accessed 01 Jan 2019].


Q. Wang, C. Yan and T. Hui, “Mechanism Design for Aircraft Morphing Wing,” Research Gate, <>, Accessed: 2019, Jan, 01, April 2012.





Want to Learn More?

We can customize each project to meet your goals.

Find the perfect solution for your needs.

Related Articles

Ship Model Scaling

The physics of model scaling make it impossible to perfectly convert between model scale and ship scale. Learn how to get around the scaling correlation problem and convert from model scale to ship scale measurements.

read more

Next Level Hydrodynamics

Time for an upgrade! Far beyond the single purpose testing facility, the modern towing tank expanded into a Swiss Army knife of experiments. Armed with new gadgets and advanced capabilities. Time to see what new tricks the tank holds for us!

read more

The Ship Towing Tank

It may look like a swimming pool, but towing tanks exist for a different purpose. Dragging models down the tank propelled the science of ship design forward across the years. These pools deliver critical measurements for ship hydrodynamics. Discover why we pay such a high price for a fundamental tool of ship design.

read more

Practical CFD Modeling: Mesh Deformation

Mesh deformation is incredibly frustrating, complicated, unstable . . . and unavoidable if you want to incorporate body motions into CFD. Modeling body motion demands mesh deformation, changing the mesh on the fly, while using it to solve transport equations. As you might expect, that brings a host of new challenges. This reviews several new strategies that the CFD engineers needs to consider.

read more

Practical CFD Modeling: Time Variation

When we add the time domain, simulations change from modeling steady scenarios to unsteady, where boundary conditions change over time. Beyond the physics, modeling unsteady flow requires a few changes to the CFD solver. Inner iterations, timestep, Courant Number, and data management all enter into the strategy for the CFD engineer. Today we discuss each of these.

read more