Multi-Phase Flow: Rising Bubble

Overview

The aim of this test case is to validate the following parameters of laminar transient multi-phase flow with a pre-referenced case [1]:

Bubble vertical velocity
Bubble center of mass
Bubble profile

The mesh was created locally using the blockMesh tool and then imported on to the SimScale platform.

Import validation project into workspace

Geometry

The domain of interest is a 2-dimensional rectangular space $2 m$

2 m

high and $1 m$

1 m

wide as shown in Fig.1.

Analysis type and Domain

A uniformly-spaced hexahedral mesh was generated using the blockMesh tool (see Fig.2.). In order to keep the flow two-dimensional, the mesh was designed to have only one layer in the $y$

y

direction.

Tool Type : OPENFOAM®

Analysis Type : Multi-phase (Laminar)

Mesh and Element types :

Table 1: Mesh Metrics¶
Mesh type	Cells in x	Cells in y	Cells in z	Number of nodes	Type
snappyHexMesh	120	240	1	28800	2D hex

Simulation Setup

Fluids:

The domain is filled with fluid $f_{1}$

f_{1}

(see Table 2). A bubble of fluid $f_{2}$

f_{2}

with diameter $0.5 m$

0.5 m

is placed at $t = 0$

t = 0

in the domain at a height of $0.5 m$

0.5 m

(with bubble center as reference) on the central vertical line. Two cases are then simulatedwith different fluid properties as per [1]. Gravity was assumed to be $0.98 m s^{- 2}$

0.98 m s^{- 2}

in the negative $y$

y

direction.

Table 2: Fluid Properties for two cases¶
Case	$r h o_{1}$ $r h o_{1}$	$r h o_{2}$ $r h o_{2}$	$m u_{1}$ $m u_{1}$	$m u_{2}$ $m u_{2}$	$s i g m a$ $s i g m a$
Case 1	$1000$ $1000$	$100$ $100$	$10$ $10$	$1$ $1$	$24.5$ $24.5$
Case 2	$1000$ $1000$	$1$ $1$	$10$ $10$	$0.1$ $0.1$	$1.96$ $1.96$

Boundary Conditions:

Table 3: Boundary Conditions for multi-phase flow simulation¶
Parameter	Top and Bottom (A,B)	Left and Right (C,D)	Front and Back (E,F)
Velocity	$0.0 m s^{- 1}$ $0.0 m s^{- 1}$	Slip	Empty
Pressure	Fixed Flux Pressure $(0)$ $(0)$	Fixed Flux Pressure $(0)$ $(0)$	Empty
Phase Fraction	Zero Gradient	Zero Gradient	Empty

Results

A comparison of the bubble profiles for the two cases with the referenced case solutions is given in Fig.3. below.

Fig.3. Visualization of Bubble Profiles and comparison with reference for case 1 (left) and case 2 (right)

Shown in Fig.4. are graphs for the development of bubble center of mass and bubble velocity over time compared with referenced results.

Fig.4. Development of bubble center of mass and bubble velocity over time compared with referenced results.

References

[1]	(1, 2) FeatFlow Benchmark Suite: Rising Bubble

Disclaimer

This offering is not approved or endorsed by OpenCFD Limited, producer and distributor of the OpenFOAM software and owner of the OPENFOAM® and OpenCFD® trade marks. OPENFOAM® is a registered trade mark of OpenCFD Limited, producer and distributor of the OpenFOAM software.

Last updated: November 25th, 2019

Did you find this article helpful?

How can we do better?

We appreciate and value your feedback.

Success

Documentation