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

Practical CFD Modeling: Volume of Fluid Modeling

Computational Fluid Dynamics (CFD) can model multiple fluids with the volume of fluid method. (VOF) The volume of fluid method opens new horizons for advanced modeling, which requires additional planning from the CFD engineer. Dive into the boundary conditions, meshing strategy, stability concerns, and more. Discover the world of VOF modeling.

read more

Practical CFD Modeling: Turbulence

Turbulence demands modeling just like any other equation in computational fluid dynamics (CFD). As the CFD engineer, you need to describe boundary conditions for your turbulence equations. This article describes how to define boundary conditions for turbulence and provides typical values for normal simulations.

read more

How to Cut the Bunker Bill: Preparing for IMO 2020

Ready or not, here it comes! No matter who you talk to, IMO 2020 promises to be a time of uncertain fuel prices. If fuel prices will go up, your fuel consumption needs to go down. Here are eight practical ways to reduce your fuel consumption.

read more

Guts of CFD: Near Wall Effects

Turbulence does tricky things near walls. Boundary layers and laminar sublayers compact interesting flow patterns into a very small space. Small it may be, but experience proved we cannot ignore it. The boundary layer forms on the body, which is our object of interest, arguably the most critical region. Turbulence is most critical near the wall, and we need to consider near wall effects.

read more

Guts of CFD: Turbulence

How we address turbulence is the defining feature of modern computational fluid dynamics (CFD). No modern computer has the power to directly compute the full details of turbulence (as of 2019). Instead, we make approximations and develop empirical models. What type of approximation, and which models should you select?

read more