Test-Case Number 1: Rise of A Spherical Cap Bubble in A Stagnant Liquid (PN)

Test-case number 1: Rise of a spherical cap bubble in a
stagnant liquid (PN)

10 April 2003
Olivier Lebaigue, DER/SSTH/LMDL, CEA/Grenoble, F-38054 Grenoble cedex 9, France

Phone: +33 (0)4 38 78 36 70, Fax: +33 (0)4 38 78 50 36, E-Mail: olivier.lebaigue@cea.fr
Christophe Duquennoy, EDF-SEPTEN, 12-14 av. Dutrievoz,

69628 Villeurbanne cedex, France, E-Mail: christophe.duquennoy@edf.fr
Stéphane Vincent, TREFLE - UMR CNRS 8508, ENSCPB

Université Bordeaux 1, 33607 Pessac cedex, France
Phone: +33 (0)5 40 00 27 07, Fax: +33 (0)5 40 00 66 68, E-Mail: vincent@enscpb.fr
1 Practical significance and interest of the test-case

This test-case could usually be considered as a very preliminary one for a new numerical
method. An extensive tester may want to reproduce most parts of the Clift, Grace and
Weber map (Clift et al. , 1978). However, this selected case deserves special attention
for the result not only consists in a final shape of the bubble (that is nevertheless a real
criteria of comparison) but also in a precise build-up of the bubble velocity, starting from
rest, exhibiting an overshoot before reaching its final asymptotic value.
To get the proper results, mainly the correct terminal velocity, and to reproduce the
overshoot, a numerical method has to take accurately into account buoyancy, viscous
stresses and surface tension effects. In particular, this test allows validating the numerical
model that takes care of jump conditions at the interface (see e.g. (Scardovelli & Zaleski,
1999)). However, the test is less severe than the ”Free rise of a liquid inclusion in a stagnant
liquid”, a test-case proposed by Lemonnier and Hervieu, presented in this volume.
2 Definitions and model description

The situation of the test-case is relative to a fluid inclusion rising in another fluid. The
inclusion and the surrounding fluid are initially at rest. Gravity induced buoyancy is the
only force inducing the motion. The test-case consists in the computation of the transient
build-up of the velocity of the rising inclusion that finally reaches a constant value.
The physical model is reduced to Navier-Stokes equation in both phases, a constant

surface tension at the interface. No phase-change takes place at the interface. As the
solution does not depend on a possible compressibility of one or both of the phases, the
test-case can be conducted in both cases (compressible / incompressible), depending on
specific features of the numerical method to be evaluated.
Reference calculations are available in non-dimensional units; however, a typical set of

dimensional physical parameter is suggested. The length scale of the problem is the initial
diameter de of the inclusion. The velocity scale for speed of displacement of the center of
mass U is,
p
Uc = gde , (1)
2/4 Test-case number 1 by O. Lebaigue, C. Duquennoy and S. Vincent
where g is the acceleration of the gravity. The time scale is therefore

p
tc = de /g. (2)
Reduced parameters are τ = t/tc and u = U/Uc . According to these definitions, the
non-dimensional reference calculation consist in the reduced time evolution of the speed
of displacement of the center of mass:
u = U/Uc = f (τ ) = f (t/tc ). (3)
The computation can be conducted either for an axisymmetrical domain or in a true

three-dimensional domain. As the limited extend of the domain has an impact on the
terminal velocity of the inclusion (see e.g. Harmathy (1960)), the size of the computational
domain has to be increased as long as an effect on the results is noted. As a rough first
estimate, we suggest that that the computational domain has a minimal extend equal to
ten diameters in all directions. According to the work of Harmathy (1960), the shape of
the bubble is not affected by the domain extension whereas the terminal rising velocity
modification can be estimated through the semi-empirical relation
Ucconf ined
2
de
≈1− , (4)
Uc∞ D
where D is a characteristic dimension of the domain in a plan perpendicular to the gravity

direction.
The physical parameters, namely ρL and ρV , respectively the density of the surround-
ing fluid and the fluid of the inclusion, µL and µV the dynamic viscosities and σ the surface
tension, are chosen to get proper values of the non-dimensional quantities for which ref-
erence computations are available: the Morton number M o, the Bond number Bo and
the ratio of densities ρL /ρV and viscosities µL /µV . The Morton number and the Bond
number are defined as usually by,
g µ4L
Mo = , (5)
ρL σ 3
and
ρL g d2e
Bo = . (6)
σ
3 Summary of the requested calculations

• Compute the displacement of an inclusion with the following non-dimensional prop-
erties M o = 0.056, Bo = 40, ρL /ρV = 100 and µL /µV = 100.
• As an example, we suggest the following physical properties: ρL = 1000 kg.m−3 ,

ρV = 10 kg.m−3 , µL = 0.273556 P a.s−1 , µV = 0.00273556 P a.s−1 , σ = 0.1 N.m,
g = 10 kg.s−2 , de = 0.02 m.
• Extract the position of the center of mass of the inclusion and then deduce its speed of
displacement. The first point of comparison is the value of the reduced asymptotic
velocity. This value can be obtained even with a peculiar point of the interface,
such as the apex. Of course, in this later case, the temporal evolution around the
overshoot (Figure 1) cannot be capture.
Test-case number 1 by O. Lebaigue, C. Duquennoy and S. Vincent 3/4
Figure 1: Reduced time evolution of the speed of displacement for two mesh sizes. After a figure of
Blanco-Alvarez (1995).
Figure 2: Recirculation zone and steam lines at reduced time τ = 2.86. After a figure of Benkenida
(1999).
4/4 Test-case number 1 by O. Lebaigue, C. Duquennoy and S. Vincent
• In addition to the main result, additional features consist in comparisons of the

non-dimensional values of the over-shoot in the build-up of velocity (Figure 1).
• Further comparisons are the shape of lines of current, the equilibrium shape of the
inclusion and the size of the recirculation zone (Figure 2). This late characteristic
requires that the inclusion has risen a length of more than ten diameters (Hnat &
Buckmaster, 1976).
References
Benkenida, A. 1999. Développement et validation d’une méthode de simulation
d’écoulements diphasiques sans reconstruction d’interfaces. Application la dynamique
des bulles de Taylor. Ph.D. thesis, Institut National Polytechnique de Toulouse, France.
Blanco-Alvarez, A. 1995. Quelques aspects de l’écoulement d’un fluide visqueux autour

d’une bulle déformable : une analyse par simulation directe. Ph.D. thesis, Institut
National Polytechnique de Toulouse, France.
Clift, R., Grace, J.R., & Weber, M.E. 1978. Bubbles, drops and particles. Academic Press.
Harmathy, T.Z. 1960. Velocity of large drops and bubbles in media of infinite or restricted
extent. American Institute of Chemical Engineers J., 6, 281–288.
Hnat, J.G., & Buckmaster, J.D. 1976. Spherical cap bubbles and skirt formation. Physics
of Fluids, 19, 182–194.
Scardovelli, R., & Zaleski, S. 1999. Direct Numerical Simulation of free-surface and inter-
facial flow. Annual Review of Fluid Mechanics, 31, 567–603.

Test-Case Number 1: Rise of A Spherical Cap Bubble in A Stagnant Liquid (PN)

Uploaded by

Copyright:

Available Formats

Test-Case Number 1: Rise of A Spherical Cap Bubble in A Stagnant Liquid (PN)

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Test-Case Number 1: Rise of A Spherical Cap Bubble in A Stagnant Liquid (PN)

Uploaded by

Copyright:

Available Formats

Test-case number 1: Rise of a spherical cap bubble in a

stagnant liquid (PN)

Olivier Lebaigue, DER/SSTH/LMDL, CEA/Grenoble, F-38054 Grenoble cedex 9, France

Christophe Duquennoy, EDF-SEPTEN, 12-14 av. Dutrievoz,

Stéphane Vincent, TREFLE - UMR CNRS 8508, ENSCPB

1 Practical significance and interest of the test-case

2 Definitions and model description

The physical model is reduced to Navier-Stokes equation in both phases, a constant

Reference calculations are available in non-dimensional units; however, a typical set of

where g is the acceleration of the gravity. The time scale is therefore

u = U/Uc = f (τ ) = f (t/tc ). (3)

The computation can be conducted either for an axisymmetrical domain or in a true

where D is a characteristic dimension of the domain in a plan perpendicular to the gravity

3 Summary of the requested calculations

• As an example, we suggest the following physical properties: ρL = 1000 kg.m−3 ,

• In addition to the main result, additional features consist in comparisons of the

Blanco-Alvarez, A. 1995. Quelques aspects de l’écoulement d’un fluide visqueux autour

You might also like