Scholarly article on topic 'Modeling of Multisize Bubbly Flow and Application to the Simulation of Boiling Flows with the Neptune_CFD Code'

Modeling of Multisize Bubbly Flow and Application to the Simulation of Boiling Flows with the Neptune_CFD Code Academic research paper on "Chemical engineering"

Share paper
OECD Field of science

Academic research paper on topic "Modeling of Multisize Bubbly Flow and Application to the Simulation of Boiling Flows with the Neptune_CFD Code"

Hindawi Publishing Corporation

Science and Technology of Nuclear Installations

Volume 2009, Article ID 953527, 8 pages


Research Article

Modeling of Multisize Bubbly Flow and Application to the Simulation of Boiling Flows with the Neptune_CFD Code

Christophe Morel1 and Jerome M. Lavieville2

1 Commissariat a l'Energie Atomique, DEN/DER/SSTH/LMDL, 17 Rue des Martyrs, 38054 Grenoble Cedex 9, France

2 R&D Division, Electricité de France, 6 Quai Watier, 78400 Chatou, France

Correspondence should be addressed to Christophe Morel, Received 19 December 2007; Accepted 9 June 2008 Recommended by Fabio Moretti

This paper describes the modeling of boiling multisize bubbly flows and its application to the simulation of the DEBORA experiment. We follow the method proposed originally by Kamp, assuming a given mathematical expression for the bubble diameter pdf. The original model is completed by the addition of some new terms for vapor compressibility and phase change. The liquid-to-interface heat transfer term, which essentially determines the bubbles condensation rate in the DEBORA experiment, is also modeled with care. First numerical results realized with the Neptune.CFD code are presented and discussed.

Copyright © 2009 C. Morel and J. M. Lavieville. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

This paper describes the modeling of boiling multisize bubbly flows and its application to the simulation of the DEBORA experiment (see [1]). First of all, the two-fluid model we use for the calculation of boiling bubbly flows is presented (see [2, 3]). In its first version, this two-fluid model assumed that the bubbles can be characterized by a single bubble diameter which is the classical Sauter mean diameter. Of course, this diameter is not a constant, but is determined as a function of the bubbles interfacial area concentration and the void fraction. As a consequence, the mean bubble diameter evolves according to several physical phenomena: bubble coalescence and breakup, but also phase change phenomena, including bubble nucleation and vapor compressibility inside the bubbles. All these phenomena are taken into account as different source terms in the interfacial area concentration balance equation.

In real flows, a bubble diameter spectrum is often observed, rather than a single bubble diameter. In order to take this diameter spectrum into account, Kamp et al. (see [4-6]) introduced a log-normal law in order to mathematically describe the bubble probability density function (pdf). As a consequence, the bubble diameter pdf is a function of

only two parameters which can be analytically expressed as functions of two particular moments of the bubble diameter distribution function. The interfacial area concentration is one of these two moments. Therefore, only one additional moment transport equation is necessary in order to be able to reconstruct the bubble diameter distribution function, hence extending our previous models to the multisize effects appearing in boiling bubbly flows. The original work of Kamp was devoted to the study of bubble coalescence in adiabatic flows under microgravity. Here, we extend his model ability to boiling bubbly flows by taking into account the vapor compressibility as well as phase change terms in the different balance equations.

The outline of the paper is the following. The classical two-fluid model for boiling bubbly flows is summarized in Section 2. In Section 3, the original model of Kamp is briefly recalled. This model is extended in Section 4 in order to take into account the vapor compressibility as well as the phase change effects, including wall nucleation, into the different balance equations. Section 5 is devoted to the modeling of the liquid-to-interface heat transfer term including the multisize effect of the bubbles. Numerical simulations of one DEBORA experimental test realized with the Neptune_CFD code are illustrated in Section 6, before concluding and

making some perspectives for future work in Section 7. This activity was realized in the context of the European NURESIM project.

2. Two-Fluid Model for Boiling Bubbly Flow

The two-fluid model we use for our boiling bubbly flow calculations is constituted of the following six balance equations (e.g., [3]):

(i) two mass balance equations:

dakpk dt

+ y-(akpkVk) = r* k = L, G,

where t is the time, at, pk, Vk denote the volumetric fraction of phase k, its averaged density and velocity and Tk is the interfacial mass transfer per unit volume and unit time; the phase index k takes the values L for the liquid phase and G for the gas phase;

(ii) two momentum balance equations:

+ v.(akPkVkVk)

= -ak Vp + Mi + akpkg + V- [ak(lk + fk )] k = L, G,

where p is the pressure, g is the gravity acceleration, Mk

is the interfacial momentum transfer per unit volume and

unit time, and T_k and tJ denote the molecular and turbulent

stress tensors, the latter being also called the Reynolds stress


(iii) two total enthalpy balance equations:

n V2 akPk[hk + —

+ V^ akPk(hk + V-) Vk)

= ak~dt + akpkg -Vk

+ hihk, + + + q'Wk - v-a(q + qJ)] k = L, G,

where hk is the phase-averaged enthalpy for phase k and hki is the interfacial-averaged enthalpy. We have assumed that the two phases are governed by the same averaged pressure field p and we make no distinction between the pressures in the two phases or between the bulk pressure and the interface pressure for simplicity. The three terms Tk, Mk, and q^a-i denote the interfacial transfer terms of mass, momentum, and heat, the quantity ai being the interfacial area concentration. The terms q'Wk denote the wall-to-fluid heat transfer per unit volume and unit time for each phase. The heterogeneous (wall) nucleation of bubbles is included as a part of the term Tk (the bubble heterogeneous nucleation corresponds to an interfacial transfer of mass between phases, even if it is located onto the heated wall surface). The two terms qk and qJ denote the molecular and turbulent heat fluxes inside phase k. Except for the liquid-to-interface heat transfer term q'I^iai (see Section 5), the modeling of

the different averaged transfer terms (interfacial, turbulent, and wall transfer terms) is largely discussed in our previous papers (see [2, 3, 7]). Here, we only present what has been changed through the introduction of the bubble multisize modeling.

3. Kamp's Model for the Coalescence of Multisize Bubbles

Kamp et al. (see [4-6]) propose to model the bubble size distribution function by a log-normal law. The bubbles are assumed to remain spherical and are characterized by their diameter d. The y's moment of the bubble distribution function and the log-normal law for the bubble diameter pdf are given by

Sy=n P(d)dY8d, P(d) =



where P(d) denotes the bubble diameter pdf and the function f = nP is the bubble diameter distribution function, n being the total number of bubbles per unit volume of the two-phase mixture. It should be noted that n and P(d) are local instantaneous quantities (defined in the sense of the ensemble average). In fact, n and P(d) should be noted as n(x, t) and P(d;x, t) but we have omitted their dependence in the position x and time t for simplicity. The diameter pdf P depends on (x, t) through the dependence of the two parameters doo(x, t) and a(x, t). It should be noted that d, x, and t are independent variables. The first four remarkable moments of the bubble diameter distribution function are related to more usual quantities by the following relations:

, ^ Si dio = —

, a¡=nS2,

a=—, (5)

where di0 is a mean bubble diameter, ai is the interfacial area concentration, and a is the void fraction (vapor volumetric fraction). An infinity of mean diameters can be constructed by using the following definition:


The diameter d\0 is one of them, the Sauter mean diameter d32 constructed with the second and third moments is another one. The two parameters of the diameter pdf are analytically expressed by the following functions of the particular moments Si and S2 and of the void fraction a (see [5]):

doo ^ —TT e nS2

-5a 2/2

Then, with these two parameters being determined, any moment defined by the first relation (4) can be expressed analytically by the following relation:

Sv = -dY-Vw-9».

With the void fraction field a being part of the solution of the equations of the two-fluid model (1)-(3), it is sufficient then to use two additional balance equations for the two particular moments S1 and S2 to completely close the geometrical quantities (provided that these two equations are themselves closed). The departure point is the so-called Liouville equation for the bubble diameter distribution function which can be derived in a very general manner as

dnp(d) + V-[nP(d)v(d)] = -d [GnP(d)]+ f(d, x, t),

where v(d) is the velocity of a bubble having a diameter d, the quantity G(d; x, t) is the Lagrangian derivative of the bubble diameter measured along the bubble path. It represents the continuous bubble size variation due to vapor compressibility and phase change. The last term f(d, x, t) represents all phenomena which are responsible for discontinuous bubble diameter changes, namely the bubble coalescence and breakup as well as bubble nucleation and collapse. Multiplying (9) by dY and integrating the resulting equation with respect to d from zero to infinity, the following balance equation for the y's moment is obtained:

— + V-(Syvy ) = yGy-iSy-i + fy,

which is used for y = 1,2. In this equation, vy, Gy, and fy are defined by the following relations:

JT v(d)dYP(d)8d STdYP(d)8d '

J T G(d)dYP(d)Sd J TdYP(d)Sd '


In the original model of Kamp, the term proportional to Gy-1 is neglected and the coalescence is the only effect taken into account in the modeling of fy. At the end, the modeled equations for S1 and S2 read (see [5])

+ V.(SiVg)

= Vt(S)"OT<21'3 - 2)S1/3f (1,Poo),

dS2 dt

+ V-(S2Vg)

^1/3Щ(22/3 - 2)S5/3 f (2,5,Poo),

where the two velocities v1 and v2 have been assumed to

be equal to the averaged gas velocity vG. In the right-

hand side of (12), eL is the dissipation rate of liquid

turbulence, Ct is the ratio of the dispersed phase turbulent

fluctuations by the continuous phase turbulent fluctuations, and f (1, a,Poo) and f (2, a,Poo) are numerical functions of the width parameter a and the coalescence probability (based on the molecular chaos assumption) P00 of two bubbles having a diameter d00 (see the original work of Kamp for more information).

4. Extension of Kamp's Model to Vapor Compressibility and Phase Change

In the original work of Kamp, the term G in (9) was neglected. Here we calculate the corresponding terms Gy-1 appearing in (10) when G is due to the vapor compressibility (Section 4.1) and to phase change (Section 4.2).

4.1. Vapor Compressibility. Let D(x, t) be the bubble diameter in physical space (the notation d stands for the bubble diameter in phase space). In the absence of phase change, the bubble mass is conserved along its trajectory. As a consequence, the bubble diameter variation and the vapor density variation measured along the bubble path are related through

DgD^D Dt dt

+ v-VD

id+и, <13»

where p denotes the vapor density and v the bubble velocity. The quantity G appearing in (9) is the conditional expectation of the Lagrangian derivative given by (13) conditioned by the equality D(x, t) = d:

DgD Dt

| D = dj

3Ш + == d) (14)

3 \ p\dt d 1 ( dpG „

-3 -^Ы + VG'VPG

If we assume that the vapor density p and the bubble velocity v do not depend on the bubble diameter D, the conditional average appearing in the second expression of (14) can be replaced by the unconditional one. The last equality in (14) simply assumes that we neglect nonlinear effect when we average the term under the brackets, therefore keeping only first order effects. This is clearly an approximation where the microscopic quantities p and v characterizing the bubble are replaced by the macroscopic (averaged) ones pG and vG.

Substituting the expression (14) into (10) and (11) gives the following expressions for the effect of the vapor compressibility on the two particular moments considered:

?1 + v. dt

(S1 Vg)

= G0S1


= ( dPG 3pG V dt


5 + ^^g)


= 2GiSi

dpG dt

4.2. Phase Change. In the previous paragraph, we made the assumption of no phase change in order to derive the compressibility terms in the equations for S1,2. Here we assume a constant vapor density p = pG in order to derive the phase change terms in the same equations. At the end, the two effects are taken into account simultaneously by assuming that the corresponding terms can be linearly superimposed into the S1,2 balance equations. For a constant vapor density, the bubble diameter variation in physical space is simply given by

DgD^SD 2m

—— = + v.vD = —, Dt dt - pG

where m denotes the bubble mass gain per unit surface per unit time due to vaporization or condensation through the bubble surface. Substituting (16) into (14) gives

g = /DD | D = d) = (2m | D = d) = 2 m | D = d). \ Dt / \pG I Pg

The closure issue consists to propose an expression for the conditional expectation (m | D = d). The simplest choice is to assume that (m | D = d) is simply given by the ratio (re -rN )/ai, where re is the vapor source per unit volume per unit time and rN is the part of this source due to newly nucleated bubbles, therefore re - rN is the part of the phase change through the surfaces of the already existing bubbles. Making this simple choice, we obtain the following expressions for the effect of the phase change on the two particular moments considered:

Í + V'(S2Vg)

n c 2 re - rn c

= G0S0 =--So,

phase change Pg ai

= 2GxSI = - su

phase change

where S0 is the volumetric number of bubbles n obtained by making y = 0 into the relation (8). A more sophisticated model taking into account the dependence of (m | D = d) on the diameter d will be presented in the next section.

4.3. A Few Words on Nucleation. In subcooled boiling flows, as this is the case for the DEBORA experiment, vapor bubbles are essentially nucleated onto the heated wall surface (heterogeneous nucleation) rather than in the liquid bulk, because the liquid is subcooled in the major part of the flow. Therefore, we can neglect a possible homogeneous nucleation in comparison to the heterogeneous nucleation. The newly nucleated bubbles are supposed to be generated

with a unique size: the so-called detachment diameter dd. In our model, this detachment diameter is calculated with the help of Unal's correlation (see [8]), used together with the model of Kurul and Podowski [9] for the active nucleation site density. These models are used to estimate the vapor produced in the meshes adjacent to the wall due to nucleation rN. Dividing rN by the mass of a nucleated bubble pGnd3dI6, the volumetric number of bubbles source term is obtained. The corresponding source terms in the S1 and S2 equations are simply obtained by multiplying this volumetric number of bubbles source term by the diameter dd or by its square d2d, according to the order of the moment considered.

5. Heat and Mass Transfers through the Surface of the Existing Bubbles

For a relatively low speed flow, the local instantaneous interfacial balance of total energy (see [10, 11]) can be simplified:

(qG - qr)-nG - m(hGi - hLI) = 0.

In this relation, q ■nG and q ■nG are the heat flux densities on the two sides of the interface, nG being the unit vector normal to the interface directed towards the exterior of the bubble. The relation (19) is valid at each point of the bubble interface. The quantities hGI and hLI are the vapor and liquid enthalpies at the bubble interface, they are generally assumed to be equal to the saturation enthalpies. The relation (19) gives the link between the rate of phase change m and the two heat flux densities between the two phases and the interface. In a bubbly flow with relatively small bubbles, the temperature of the vapor inside the bubbles is generally close to the saturation temperature. In subcooled liquid, the rate of bubble condensation is essentially determined by the heat exchange on the liquid side, rather than the one on the vapor side. Neglecting qG^nG in comparison to qL^nG in the previous equation and introducing a Nusselt number NuL to characterize the liquid-to-interface heat transfer, (19) becomes


hGi - hLi (hGi - hLi)D

(0L - 0sat), (20)

where XL is the liquid thermal conductivity and 0L and 0sat represent the liquid and saturation temperatures, respectively. The Nusselt number is given by the Ranz-Marschall law [12]

Nul = 2 + 0.55Re1/2PrL/3,

where Re and PrL are the bubble Reynolds number and the liquid Prandtl number, respectively. Neglecting the fluctuations of the enthalpies hGI and hLI, the averaging of the first relation (20) gives

q'Liai = Tg^gi - hLi) with

rG = J(m I D = d)nd2nP(d)8d, (22)

q'LIaI = J(qLn I D = d)nd2nP(d)8d,

where rG and q'ilal are the average interfacial mass and heat transfer terms appearing in (1) and (3). The term rG is related to the vapor phase because m is the bubble mass gain and q'ilal is the liquid-to-interface heat transfer per unit volume per unit time, the heat transfer density qLl being expressed per unit area per unit time, ai being the interfacial area concentration. Equation (22) is very useful because it directly relates the mass transfer rate to the liquid-to-interface heat transfer rate (under the assumption of a nearly saturated vapor inside the bubbles). The averaging of the second relation (20) gives

qLIai =

/Nu/X/ \ D

(di - 0sat) I D = d)nd2nP(d)Sd.

The relations (21) and (23) are the basic tools to calculate the average interfacial liquid-to-interface heat transfer. In what follows, we give three different approximated expressions for this term. The first one is a single-size approximation of the liquid-to-interface heat flux, it is given by

qïiai = I ndNulXl(dl - 6sat)nP(d)Sd Nui(ds2)-

'-X/(©/ - ©Sat).

In this expression, it is implicitly assumed that all the bubbles have the same diameter which is given by the Sauter mean diameter d32 = 6a/al (single-size approximation). The relation (21) is directly calculated as a function of the Sauter mean diameter and the average relative velocity of the bubbles into the Reynolds number. The second (more precise) approximation consists in using the log-normal law for the bubble diameter pdf (4). The result is

qiiai =

2nXiS1 ■

■0.55nXj.-Ц^ Ivr |1/2nd002

00 exH 85

x (©L - ©sat).

The dependence of the bubble Reynolds number Re on the bubble diameter d has been taken into account in the derivation of (25), but not the variation of the relative velocity \vr|, also appearing in the Reynolds number, as a function of the bubble size. When using (25), the norm of the relative velocity is assumed to be given by the averaged fields. In order to alleviate this approximation, a third model has been constructed by assuming that the norm of the relative velocity can be accurately given by the terminal velocity

4 Apgd 3piCD'

where Ap is the absolute value of the density difference \pL -pG\ and CD is the drag coefficient, which can be given for spherical bubbles by (e.g., [13])

Cd = 24(1 + 0.1Re3/4). Re

Assuming that the bubble Reynolds number is sufficiently high in order to approximate the preceding relation by CD = 2.4/Re1/4, we obtain the following relation for the square root of the terminal velocity:

APg 1.8pcv}/4


Reevaluating the liquid-to-interface heat transfer from (4), (21), (23), and (28) gives a third model:

qiiai =

2nXLS1 + 0.55nXL l/2

x nd13/7 exp ( 169 52


(©/ - ©sat).

Now, having the model given by (20)-(21) for the bubble mass gain m, we can evaluate more precisely the corresponding terms in the equations for the two moments Si and S2, namely the terms G0 and G1 with G given from (17). The results are the following:

lt + V-(S1Vg)

phase change

Pg hgi - h/i d00 xi pr/

n /а2 exp

Pr//3 ( Apg

_ Pg hgi - h/i vi '

x (©/ - ©sat),


n ( а2


lt + V-(S2V-g) 8 X/

phase change

-n(©L - ©sat)

2.2 X/ Pr//3/ Apg

Pg hGi - h/i v//2 V 1.8p/v/

X (©/ - ©sat).

2/7 ( 18 \ nd607 exp ^ 4"8 $2j

To conclude the modeling sections, we can summarize the complete modeling of boiling bubbly flows with multisize bubbles in its present state. The basic equations are the mass, momentum and total enthalpy balance equations written for the two phases, they are given by (1)-(3). A multisize model is used in conjunction with the two-fluid model to evaluate the evolution of the bubble diameter spectrum under several physical effects: coalescence, vapor compressibility, and phase-change. It should be noted that a bubble breakup term is removed in the present state of the model. A log-normal law is used to describe the bubble diameter pdf. This law involves two parameters d00 and a which are analytically determined from two particular moments of the bubble diameter pdf (7). Therefore, a set of two additional transport equations are written for these two moments. These moments balance equations have the

general form (10) with their right-hand side including a coalescence model (12), a vapor compressibility model (15), and two possible phase change models (18) or (30). The model also includes wall nucleation source terms which are not presented here (see [2]). At the end, we have proposed several expressions for the liquid-to-interface heat transfer term appearing in the liquid enthalpy balance (3). This term is responsible for the bubble condensation in the core of the flow. The first one (24) is the classical expression used in usual (single size) bubbly flow models. The two others (25) and (29) take into account the bubble diameter spectrum as modeled by the log-normal law. The last one (29) also takes into account the dependence of the bubble relative velocity on the bubble diameter. All these models should be tested. In the following section, we give the results of a comparison of Neptune_CFD simulations with an experimental test of the DEBORA database.

6. Numerical Simulations of a Boiling Bubbly Flow

The DEBORA experiment, carried out at the French Commissariat a l'Energie Atomique (see [1]), has been chosen to evaluate our model. In this experiment, the R-12 has been adopted as the working fluid to simulate the PWR (pressurized water reactor) conditions under low pressure. Some liquid R-12 flows upwardly inside a vertical pipe having an internal diameter equal to 19.2 mm. A part of the tube wall of 3.5 m long is electrically heated. Numerous vapor bubbles are nucleated onto the heated wall surface, and condense in the core of the flow, where the liquid is subcooled. At the end of the heated section, the radial profiles of the void fraction, bubble diameter, and vertical bubble velocity are measured by means of an optical probe and the liquid temperature is measured by means of a thermocouple.

Here we evaluate the different models presented above on a particular DEBORA experimental test. The controlling parameters of this test are the following.

(i) Test pressure: 14.59 bar.

(ii) Inlet mass flux density (liquid only): 2027 kg/m2s.

(iii) Liquid inlet temperature: 28.52°C.

(iv) Wall to fluid heat flux density: 76200 W/m2.

// '' ■// /


// '' . ■ // ' □

V/ D □

0.002 0.004 0.006 0.008 0.01

□ Experiment

- Eqs. (18) & (24)

...... Eqs. (18) & (25)

Eqs. (18) & (29)

---Eqs. (30) & (29), hki = hk

- Eqs. (30) & (29), hki = h^t

Figure 1: comparison of the different models on the void fraction profile.

0.004 0.006 r (m)

□ Experiment

- Eqs. (18) & (24)

...... Eqs. (18) & (25)

Eqs. (18) & (29) --- Eqs. (30) & (29), hki = hk ■■ Eqs. (30) & (29), hki = hksA

Figure 2: comparison of the different models on the vapor velocity profile.

All the models have been implemented in the multiphase Neptune_CFD code. With the flow being axi-symmetric, we use a two-dimensional axi-symmetric calculation grid. This grid is characterized by 80 axial meshes and 10 radial meshes. We assume that this relatively coarse grid is sufficient for a first test of the different models. The models tested are the three expressions of the liquid-to-interface heat flux (24), (25), and (29) and the multisize model for bubbly flows with the two different expressions (18) and (30) for the phase change terms. The compared quantities are the void fraction (Figure 1), the vapor mean velocity (Figure 2), the liquid temperature (Figure 3), the interfacial area concentration (Figure 4), the Sauter mean diameter d32 (Figure 5), and the mean diameter d20 defined by (6) (Figure 6). All the radial

profiles (experimental and calculated) have been taken at the end of the heated length.

Each figure compares the results of five calculations together with the experimental values.

According to these comparisons, the model giving better results on the void fraction, the vapor mean velocity, the liquid temperature, and the interfacial area concentration seems to be the one using (29) and (30) with the interfacial averaged enthalpies being equal to the phase averaged ones. It can also be seen that the void fraction is strongly overestimated for all the models tested, but the interfacial area concentration is overestimated too, therefore giving approximately the right values for the bubble Sauter mean diameter and the mean diameter d20. The overestimation of

'3 a £

0.005 r (m)

о Experiment

- Eqs. (18) & (24)

...... Eqs. (18) & (25)

Eqs. (18) & (29)

---Eqs. (30) & (29), hki = hk

- Eqs. (30) & (29), hki = hksl

Figure 3: comparison of the different models on the liquid temperature profile.

Й 1000

0.004 0.006 r (m)

0.004 0.006 r (m)

□ Experiment Eqs. (18) & (29)

- Eqs. (18) & (24) ---Eqs. (30) & (29), hki = hk

-■ Eqs. (18) & (25) ■-■- Eqs. (30) & (29), hki = hksat

Figure 5: comparison of the different models on the Sauter mean diameter profile.

0.5 - .

0.004 0.006 r (m)

on 3000

Eqs. (18) & (29)

---Eqs. (30) & (29), hki = hk

■■ Eqs. (30) & (29), hki = hM

Figure 4: comparison of the different models on the interfacial area profile.

□ Experiment

- Eqs. (18) & (24)

...... Eqs. (18) & (25)

□ Experiment

- Eqs. (18) & (24)

...... Eqs. (18) & (25)

Eqs. (18) & (29)

---Eqs. (30) & (29), hki = hk

- Eqs. (30) & (29), hki = hsA

Figure 6: comparison of the different models on the mean diameter d20 profile.

the void fraction could be attributed to the inadequacy of the Ranz and Marschall law (21) to describe the condensation of bubbles in a subcooled liquid. This issue will be studied in future work.

7. Conclusions

A multisize model for boiling bubbly flows has been presented in detail. This model is not completely achieved (e.g., a bubble breakup term is missing in the equations). Nevertheless, we made first calculations to evaluate the capabilities of this model in its present state. Five different calculations have been done and compared to a single experimental test of the DEBORA experiment. It has been shown that the results are sensitive to the expressions used for

the liquid-to-interface heat transfer (24), (25), or (29) and to the phase change terms in the geometrical equations (18) or (30). The model using the particular expressions (29) and (30), with the interfacial averaged enthalpies being given by the phase averaged ones, gives better results. Unfortunately, the void fraction and the interfacial area density are strongly overestimated. Further investigations are needed to clarify these questions.

Several issues can be raised for future developments:

(i) a bubble breakup model should be added for completeness of the model;

(ii) a bubble collapse model is missing too;

(iii) the physical model given by Ranz and Marschall [12] (21) could not be adapted for condensing bubbles

in subcooled liquid. The validity of this model should be evaluated first (perhaps by using DNS calculations).


Neptune_CFD is a three-dimensional two-fluid code developed more especially for nuclear reactor applications within the framework of the Neptune project, financially supported by CEA (Commissariat à l'Energie Atomique), EDF (Electricité de France), IRSN (Institut de Radioprotection et de Surete Nucleaire), and AREVA-NP. This activity was realized in the context of the NURESIM project financed by the European Union.


[1] J. Garnier, E. Manon, and G. Cubizolles, "Local measurements on flow boiling of refrigerant 12 in a vertical tube," Multiphase Science and Technology, vol. 13, no. 1-2, pp. 1-111, 2001.

[2] W. Yao and C. Morel, "Volumetric interfacial area prediction in upward bubbly two-phase flow," International Journal of Heat and Mass Transfer, vol. 47, no. 2, pp. 307-328, 2004.

[3] C. Morel, S. Mimouni, J. M. Lavieville, and M. Boucker, "R113 boiling bubbly flow in an annular geometry simulated with the NEPTUNE code," in Proceedings of the 11th International Topical Meeting on Nuclear Reactor Thermalhydraulics (NURETH-11), Avignon, France, October 2005, paper: 248.

[4] A. M. Kamp, Ecoulements turbulents a bulles dans une conduite en micropesanteur, Ph.D. thesis, Institut National Polytechnique de Toulouse, Toulouse, France, 1996.

[5] A. M. Kamp, A. K. Chesters, C. Colin, and J. Fabre, "Bubble coalescence in turbulent flows: a mechanistic model for turbulence-induced coalescene applied to microgravity bubbly pipe flow," International Journal of Multiphase Flow, vol. 27, no. 8, pp. 1363-1396, 2001.

[6] C. Colin, X. Riou, and J. Fabre, "Turbulence and shear-induced coalescence in gas-liquid pipe flows," in Proceedings of the 5th International Conference on Multiphase Flow (ICMF '04), Yokohama, Japan, May-June 2004, paper no. 425.

[7] C. Morel, J. Pouvreau, J. M. Lavieville, and M. Boucker, "Numerical simulations of a bubbly flow in a sudden expansion with the NEPTUNE code," in Proceedings of the 3rd International Symposium on Two-Phase Flow Modeling and Experimentation, Pisa, Italy, September 2004.

[8] H. C. Unal, "Maximum bubble diameter, maximum bubble-growth time and bubble-growth rate during the subcooled nucleate flow boiling of water up to 17.7MN/m2" International Journal of Heat and Mass Transfer, vol. 19, no. 6, pp. 643-649, 1976.

[9] N. Kurul and M. Z. Podowski, "Multidimensional effects in forced convection subcooled boiling," in Proceedings of the 9th International Heat Transfer Conference, vol. 1, pp. 21-26, Jerusalem, Israel, August 1990, paper no.-04.

[10] J. M. Delhaye, Thermohydraulics of Two-Phase Systems for Industrial Design and Nuclear Engineering, McGraw-Hill, New York, NY, USA, 1981.

[11] M. Ishii and T. Hibiki, Thermo-Fluid Dynamics of Two-Phase Flow, Springer, New York, NY, USA, 2006.

[12] W. E. Ranz and W. R. Marschall, "Evaporation from drops," Chemical Engineering Progress, vol. 48, pp. 173-180, 1952.

[13] M. Ishii, "Two-fluid model for two-phase flow," Multiphase Science and Technology, vol. 5, no. 1-4, pp. 1-63, 1990.

Copyright of Science & Technology of Nuclear Installations is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.