Multiphase flow simulation using CLSVOF methods
(May 2020 - June 2020)
As part of a course on advanced CFD, I wrote a sequential program to solve for three-dimensional incompressible two-phase flows using the CLSVOF method, originally proposed by Sussman and Puckett (2000). The code also implements the improvements suggested by Griebel and Klitz (2016) and Son (2011). Results from two test-cases (Griebel and Klitz, 2016) are shown below. It can be seen that the mass is conserved to an excellent degree, along with proper interface-tracking.
Two-dimensional stretching of a circle in a shear velocity field
The figures below track the interface between a circle (phase I) and a fluid (phase II) over time, when a shearing velocity field (periodic with time-period = 8 sec) is applied. The aim is to recover the initial interface (t = 0) at the end of the time-period (i.e., at t = 8).
t = 0, 8 sec
t = 2 sec
t = 4 sec
t = 6 sec
Three-dimensional deformation of a sphere in a shear velocity field
The figures below track the interface between a sphere (phase I) and a fluid (phase II) over time, when a shearing velocity field (periodic with time-period = 3 sec) is applied. The aim is to recover the initial interface (t = 0) at the end of the time-period (i.e., at t = 3).
t = 0 sec
t = 0.48 sec
t = 0.98 sec
t = 1.49 sec
t = 1.98 sec
t = 2.48 sec
t = 3 sec
