Research Article  Open Access
C. Baudry, M. Guingo, A. Douce, J. Laviéville, S. Mimouni, M. Boucker, "Numerical Study of the SteadyState Subchannel TestCase with NEPTUNE_CFD for the OECD/NRC NUPEC PSBT Benchmark", Science and Technology of Nuclear Installations, vol. 2012, Article ID 524598, 11 pages, 2012. https://doi.org/10.1155/2012/524598
Numerical Study of the SteadyState Subchannel TestCase with NEPTUNE_CFD for the OECD/NRC NUPEC PSBT Benchmark
Abstract
The multifield computational fluid dynamics (CFD) code NEPTUNE_CFD is applied to carry out a numerical study of the steadystate subchannel testcase of the OECD/NRC NUPEC PWR subchannel and bundle tests (PSBTs) international benchmark, focusing on the simulation of a subset of five selected experimental runs of the centered subchannel configuration. First, using a standard choice for the physical models and a constant, predetermined bubble diameter, the calculated void fraction is compared to experimental data. Besides, the mesh sensitivity of the calculated void fraction is investigated by performing simulations of three grid levels, and the propagation of the experimental uncertainties on the input parameters of the simulations is also studied. Last, calculation results with devoted models for the bubblesize distribution are analyzed. Their impact is visible on the subcooled run, giving void fraction closer to experiments than those obtained with a fixed bubblesize. Voidfraction distribution with bubblesize models is also shown to come closer to experiment for another run with a higher equilibrium quality.
1. Introduction
The OECD/NRC PWR subchannel and bundle Tests (PSBT) benchmark was an international project endorsed by the OECD/NEA (Organisation for Economic Cooperation and Development/Nuclear Energy Agency) and supported by US NRC (United States National Regulatory Commission) and METI (Japanese Ministry of Economy, Trade and Industry), in which a large experimental database of voidfraction measurements performed at NUPEC (Nuclear Power Engineering Corporation) under PWR (Pressurized Water Reactor) thermalhydraulic conditions in different geometric configurations (different types of isolated subchannels or rod bundle) has been made available to the participants for numerical simulation. One of the purposes of this benchmark is to provide experimental data that can be used for the validation of numerical models of voidfraction distribution over a wide range of operating conditions, and for the development of novel approaches. PSBT was not a blind benchmark, in the sense of the experimental value of the variables of interest were provided to the participants. The benchmark was organized in two phases; the first focusing on a voiddistribution benchmark while the topic of the second was the departure from nucleate boiling (DNB) phenomenon. More details on the benchmark organization can be found at [1]. Each phase encompassed several exercises: the current paper focuses on the simulation of the steadystate subchannel exercise of phase I (i.e., Exercise I1). For each simulated run of this exercise, three main results were asked from the participants. (i)A numerical calculation of the crosssection void fraction at a given elevation, to be compared to experimental data. (ii)A 2D view of the voidfraction distribution in this section. Ideally, this result could also be compared to an experimental view obtained by chromotomography. However, this comparison was proven difficult due in particular to the fairly low resolution of the experimental view. Therefore, in the following, only one of these experimental views will be used for comparison with numerical simulations. (iii)The evolution of the crosssection averaged void fraction with respect to the elevation. No experimental data was available for this result. However, it seemed interesting to analyze how the different codes predicted, for instance, the transition to fullydeveloped subcooled boiling for some of the runs.
This paper aims at presenting and analyzing simulations carried out with the multifield computational fluid dynamics (CFD) code NEPTUNE_CFD [2] for Exercise I1. Following the NEA/CSNI best practice guidelines [3], a mesh sensitivity analysis is performed on three grid refinements, and the influence of the experimental uncertainties on the numerical results is assessed. As a first step, the standard set of parameters of the code is applied, with in particular the use of a secondorder (Reynolds stress) turbulence model, as well as a constant, predetermined bubble diameter. In the last part of the paper, the impact of more detailed descriptions of the bubblesize distribution is presented.
2. Experimental Configuration of the SteadyState Subchannel Exercise
The test facility of the steadystate subchannel exercise represents one of the subchannel types found in a PWR assembly. The effective heated length is 1.555 m, while the measurement of the void fraction takes place at 1.4 m from the bottom of the heated section. The external diameter of the rod is 9.5 mm; the rod pitch and the rod gap measure, respectively, 12.6 mm and 3.1 mm. The heating power is uniformly distributed. Additional details about the experimental apparatus may be found in [4]. The void fraction is measured by using the chromotomography (CT) technique, which also gives the local distribution of the timeaveraged void fraction at the measuring section [5]. The initial set of experimental runs of Exercise I1 counted approximately 40 runs, taking place in different subchannel types (center, center with a guide tube, side and corner). However, after discussion between participants and PSBT organizers, five runs have been selected to be more precisely analyzed: (i)the run whose operating conditions are the closest to PWR normal conditions (reference 1.2211), (ii)two runs with a higher inlet temperature and a relatively low wall heat flux (references 1.2223 and 1.2237), (iii)two runs differing only by the value of the inlet temperature (references 1.4325 and 1.4326), with a lower heat flux and a lower pressure than run 1.2211. Run 1.4325 was also chosen due to its particular voidfraction distribution at the measuring section, mainly located in the nonheated corners of the domain.
The value of the controlling parameters of the runs and the crosssection averaged voidfraction measurements are summarized in Table 1. Furthermore, the estimated uncertainties on the flow parameters and on the voidfraction measurement provided by the experimentalists are given in Table 2. It can be observed that run 1.2211 is the only selected run in subcooled conditions in the whole subchannel (with a negative thermal equilibrium quality at the measuring section).


3. The NEPTUNE_CFD Code
3.1. Main Features of the Code
NEPTUNE_CFD is a 3D, multifield CFD code developed in the framework of the NEPTUNE project, financially supported by CEA (Commissariat à l’Énergie Atomique et aux Énergies Alternatives), EDF, IRSN (Institut de Radioprotection et de Sûreté Nucléaire) and AREVANP. It is mainly devoted to the study and the simulation of nuclear reactor applications involving multiphase flows, such as twophase pressurized thermal shock (PTS) and departure from nucleate boiling (DNB) applications. The need to simulate PTS applications have led in particular to the development and the numerical implementation of specific methods of interface detection for free surface flows [6]. The code follows the classical multifield onepressure formulation [7]. For each phase , the basic set of resolved conservation equations (mass, momentum and total enthalpy denoted ) is the following: Additionally, the following relation on volume fractions holds: In these equations, is the volume fraction, is the density, is the velocity, and are, respectively, the laminar and the turbulent stress tensor; is the interfacial mass transfer, represents the momentum transfer from all the other phases to phase , represents the heat transfer between phases, is the pressure, is the temperature, and is the conductivity. Additional jump relations are supplied for the interfacial transfers.
The spatial discretization is a full unstructured finitevolume approach with a collocated arrangement of allvariables. The numerical algorithm used is a semiimplicit, pressurebased method where the system of equations is solved in two major fractional steps: first, a prediction of the velocities based on the momentum equations; then, the coupling between phase fraction, pressure, and energy through mass and energy equations and a simplified form of momentum equations [8].
3.2. Physical Modeling for the Current Study
The turbulence of the liquid phase is modeled using a secondorder, RANS (ReynoldsAveraged NavierStokes) model [9] including bubbleinduced turbulence effects, whereas a turbulent dispersion model is applied on the gas phase [10]. The interfacial transfer of momentum is considered as the sum of different contributions, namely, the drag force (modeled by using the correlation developed by Ishii and Zuber [11]), the added mass force (by using the expression of Zuber [12]) and the formulation of the lift force proposed by Tomiyama et al. [13]. As a general feature of the NEPTUNE_CFD code, no flowregime map is used, as the flow is regarded as bubbly.
At the heating wall, the heat transfer model is an extension of the approach of Kurul and Podowski [14], which is often referred to as the RPI (rensselaer polytechnic institute) model and which consists in splitting the heat flux into three terms: one heating the liquid phase in contact with the wall, one responsible for the bubble generation, and the last one arising from the arrival of liquid water at the wall, caused by bubble departure (the socalled “quenching” flux). When the void fraction in the boundary cells is sufficiently high, a fourth flux is introduced to take into account the convective heat transfer transmitted to the vapor. In the bulk, the recondensation model is based on the RanzMarshall correlation.
For the first series of calculations, a constant, predetermined diameter for the bubbles has been taken equal to 0.3 mm, based on previous experimental observations carried out in PWR conditions [15]. It can be thus expected for this modeling choice to behave more satisfactorily for runs with a relatively low void fraction (i.e., where coalescence and fragmentation phenomena can be neglected). As a second step, the influence of a more refined description of the bubblesize distribution has been studied by performing simulations using the interfacial area models of Yao and Morel [16] and Ruyer and coworkers [17], which have been validated in PWR conditions in a vertical duct geometry.
4. Computational Strategy
4.1. Time Convergence
In order to reach the steady state, a transient algorithm is used in NEPTUNE_CFD. The time step chosen is timedependent, with a maximal CFL number set to 1 for stability reasons and convergence of the iterative numerical algorithm.
4.2. Characteristics of the Grids
Following the NEA/CSNI (Committee on the Safety of Nuclear Installations) best practice guidelines, three grids with different refinements have been used for this study. Their main features are recapitulated in Table 3; crosssection views of the grids are proposed in Figures 1, 2, and 3. The normalized wallnormal distance of the boundarycells centers ranges from approximately 75 (run 1.4326, grid level 3) to approximately 320 (run 1.2211, grid level 1), supporting the use of a highReynolds formulation of a secondorder RANS turbulence model. It can be seen that the refinement factor from grid level 1 to grid level 3 is approximately equal to 3 following two perpendiculars belonging to a crosssection, and 1.3 in the axial direction.

4.3. Boundary and Initial Conditions
The profiles of the inlet velocities and turbulent variables (Reynolds stresses and mean dissipation) have been calculated beforehand through devoted calculations, so as to get developed flow profiles. The procedure followed for these calculations amounts to calculating a case without thermal power in an infinite channel of the same cross section; and then selecting the profiles of the variables where they do not vary anymore. The initial conditions for velocities, temperature and turbulent variables are taken equal to the inlet conditions.
5. Results Obtained with a Constant, Predetermined Bubble Diameter
5.1. Void Fraction at the Measuring Section
The calculations have been carried out on 36 to 144 cores on the EDF R&D Ivanhoe cluster (Westmere 2.93 GHz Infiniband). The time to reach the steady state represents approximately 2 to 4 seconds of “physical” time. First, calculations were run under nominal conditions on the three grid levels to investigate the grid sensitivity of the results. The values of the crosssection averaged void fraction at the 1.40meter section are summarized in Table 4. It can be observed that simulations performed on grid level 1 give voidfraction values come very close to those obtained on grid level 3 (the discrepancy is of order of ±0.01).

Concerning the comparison with the experimental values for the crosssection averaged void fraction: (i)for run 1.2211 (subcooled), the simulations overestimate the void fraction by 0.03; thus, reaching the upper bound of the 1uncertainty on the voidfraction measurement. (ii)for runs 1.2223 and 1.2237 (high inlet temperature and low wall heat flux), the simulations underestimate the experimental void fraction by 0.06 to 0.08; (iii)for runs 1.4325 and 1.4326 (low pressure), the void fraction is overestimated by 0.1. On average, the deviation between numerical results and experimental data is then of the order of 0.06 voidfraction units; and consequently more visible on the subcooled run (1.2211).
5.2. Axial Evolution of the Flow
The axial evolution of the crosssection averaged void fraction for the five runs are presented from Figures 5, 6, 7, 8, and 9, as well as the axial evolution of the wall temperature at a given location of a heating wall (see Figure 4).
It can be seen that for almost every run (except for the subcooled run 1.2211), vapor is generated right from the start of the heated section. Furthermore, concerning runs 1.2223, 1.2237, and 1.4325, a change of the slope in the evolution of the mean void fraction is visible, respectively, for an elevation approximately equal to 0.9 m, 0.6 m, and 0.9 m. This change occurs at the same elevation as an increase of the wall temperature, indicating the beginning of the fullydeveloped subcooled boiling (FDB) [18]. Runs 1.4325 and 1.4326 exhibit approximately the same behavior. For instance, for run 1.4325, different regimes are observed: (i)from the bottom of the subchannel to approximately 0.7 m, the mean void fraction is rising slowly, corresponding to the partial subcooled boiling regime; (ii)from 0.7 m up, the void fraction rises more rapidly, corresponding to the FDB regime.(iii)finally, towards the end of the domain, the wall temperature jumps from 605 K to several thousands Kelvin (not represented): in this region, the calculated nearwall void fraction exceeds 0.8, which is the threshold value used in the code to switch continuously from the nucleate boiling to a pseudofilm boiling. In this model, the totality of the imposed heat flux is transferred to the vapor (rising its temperature above saturation), therefore, bringing about a drop of the heat transfer coefficient. Let us note that this model is not yet considered as validated. Furthermore, Figure 10 represents the evolution of the average liquid temperature and of the saturation temperature for run 1.4325. At the end of the domain, the two temperatures are almost equal, corresponding to the saturated boiling regime.
5.3. Propagation of the Experimental Uncertainties
As the estimated uncertainties on the controlling parameters of the flow are provided (see Table 2), it is interesting to analyze the effects of these uncertainties on the simulations. For this purpose, a series of calculations has been performed in which one of the controlling parameters is taken equal to its nominal values plus or minus the value of the relative uncertainties, the other parameters remaining fixed to the nominal values. This procedure makes it possible to study the effect of one uncertainty independently from the others. In order to keep a reasonable timetoresult, the coarsest grid has been used (grid level 1). The results obtained for the void fraction are recapitulated in Table 5.

These series of calculations show that for a relatively small variation of any parameter, the result on the averaged void fraction is modified by approximately 0.01, whatever its absolute value. Consequently, the relative impact on the voidfraction value of a single parameter variation is more visible on the subcooled case, tha is, the run that shows the lowest voidfraction value. Furthermore, it can be seen than within the bounds provided for the experimental uncertainties, no input parameter seems to have a foremost importance with respect to the others on the calculated void fraction.
As a second step, two additional series of calculations have been carried out: for each series, the set of parameters expected to give the highest or the lowest value of the averaged void fraction has been chosen. The results obtained are reported in Table 6. These simulations show that taking the two “opposite” sets of parameters give numerical results that may span on a range quite wide: for instance, for run 1.4325, the voidfraction value varies on a 0.08wide interval, representing approximately 15% of the nominal value. The impact is more visible for run 1.2211, where a calculated 0.03wide variation represents nearly 50% of the nominal value.

6. Effect of More Refined Descriptions of BubbleSize Distribution
The calculations presented so far have been obtained by using a constant and uniform predetermined bubble diameter equal to 0.3 mm. In order to further investigate the five runs, simulations have been run by using dedicated models for the bubblesize distribution, namely, the model proposed by Yao and Morel [16], further referred to as the YaoMorel model, and the model developed by Ruyer and coworkers [17, 19], further referred to as the RuyerSeiler model. The YaoMorel model follows the socalled “singlesize” approach for bubbly flows [20]. This approach considers that the bubbles have locally the same size, often represented by the Sauter mean diameter, which is directly connected to the local void fraction and to the interfacial area, the latter being obtained by the resolution of an additional transport equation. This equation takes into account bubblesize variation caused by gas compressibility, and contains terms standing for bubble coalescence and fragmentation that need to be modeled.
The RuyerSeiler model can be regarded as more refined since it follows the “momentdensity” approach, which consists of assuming a certain form for the bubble diameter distribution function, and then solving equations on the moments defining this distribution. The RuyerSeiler model assumes a quadratic form for the bubblediameter distribution whose graph extends from the origin and is consequently defined by only one parameter. Another possibility of postulated form for the bubblediameter distribution is a lognormal law, as proposed for instance by Kamp and coworkers [21, 22].
In this section, we specifically focus on two runs: run 1.2211, whose conditions are the closest to PWR normal conditions, and run 1.4325, which exhibits an interesting void fraction distribution according to the experimental 2D view obtained by chromotomography.
6.1. Mesh Sensitivity Analysis
First, the two interfacialarea models have been applied on the three grid levels and for the five runs, to analyze the sensitivity of the results to the grid levels. The averaged void fractions obtained are summarized in Table 7, where it is shown that the discrepancy is of the order of 0.01 voidfraction units from grid level 1 to grid level 3, as it was the case for the computations with a fixed bubble diameter. Figures 11, 12 represent the longitudinal distribution of the averaged void fraction for the three grid levels for run 1.2211, and show that the mesh influence remains in the 0.01limit in the whole domain.

6.2. Analysis of Run 1.2211 with InterfacialArea Models
As shown in Figure 13, the behavior of the two bubblesize models is quite similar concerning the averaged void fraction: From the whole heated length, they predict a void fraction lower than the one obtained with a fixed bubble size. Besides, at the measuring section, the void fraction predicted by the models is closer the experimental value.
This underprediction with respect to the previous calculations can be connected to the size of the bubbles calculated by the models at the measuring section. The calculated average, minimum, and maximum size of the bubbles at the measuring section are summarized in Table 8. It can be observed that the mean bubble size calculated by the devoted models is smaller than the initial choice of 0.3 mm used for the first calculation. This smaller size brings about an increase of the coreflow recondensation phenomenon (due to the higher interfacial area of smaller bubbles), and consequently a decrease of the void fraction.

Crosssection views of the bubblediameter distributions at the measuring section are proposed in Figures 14, 15, showing quite a similar behavior; with the biggest bubbles (0.2 mmdiameter for the RuyerSeiler model, 0.15 mm for the YaoMorel model) found above the heated walls, whereas the smallest are located in the center and the corners of the domain. To provide a detailed representation, the color scaling has been chosen as the most adapted for each case and is consequently not the same on the two figures.
6.3. Analysis of Run 1.4325 with InterfacialArea Models
As can be seen in Table 7, the averaged void fraction for run 1.4325 obtained by using interfacialarea models is very similar to the one obtained with a fixed bubblesize diameter (of the order of 0.45). However, the void fraction distribution predicted by the three calculations present notable discrepancies, as presented by Figures 16, 17, and 18 (with a color scaling adapted to the case).
With a fixed bubble diameter set to 0.3 mm, the void fraction is mostly concentrated in a thin layer adjacent to the heated wall, reaching locally peak values of approximately 0.85 and thus triggering the activation of the model of pseudofilm boiling, as described in Section 5.2. With dedicated bubblesize models, on the opposite, consider the following.(i)The peak value of the void fraction predicted by applying the bubblesize models is quite lower (0.65 for the YaoMorel model; 0.55 for the RuyerSeiler model). Lower voidfraction peak values calculated by the bubblesize models cause the flow to remain in the nucleate boiling remain for the whole elevation, as presented in Figure 19. (ii)Besides, for these calculations (and especially with the RuyerSeiler model), high voidfraction values are found near the nonheated corners of the domain, phenomenon that can be seen on the experimental view. The effect of the lift force, which depends on the bubble size, may account for this phenomenon, as the average bubblesize predicted by the RuyerSeiler model at the measuring section is approximately equal to 1 mm, thus, very different from the calculations with a fixed bubble size. It should be noted that, even if the mean bubble diameter is different, the averaged void fraction is only slightly impacted since recondensation is very small in this case.
7. Conclusion
In this paper, we reported simulations with the multifield CFD code NEPTUNE_CFD in the framework of the first exercise of the OECD/NRC PSBT benchmark, which were carried out on a centered, isolated subchannel geometry and were compared to experimental data. A standard set of physical models was used, as well as a fixed bubble diameter. In nominal conditions, the discrepancy between calculated and experimental mean void fraction at the measurement section is on average of the order of 0.06 voidfraction units. Following the NEA/CSNI best practice guidelines, the impact of the grid on the axial evolution of the mean void fraction has been investigated by using three grid levels, and shown to be of the order of 0.01 voidfraction units. Furthermore, a study of the numerical propagation of the experimental uncertainties has been carried out, providing an envelope of calculated voidfraction varying from 0.05 to 0.1 around the mean value. This sensitivity analysis could be deepened by sampling all input parameters together, in order to study the coupling between these parameters.
As a second step, models allowing to simulate the dispersion in size of the bubbles have been applied. The use of these models impacted only slightly the axial evolution of the mean void fraction (of the order of 0.01 to 0.02 voidfraction units compared to the fixeddiameter calculation), and was more clearly seen on the subcooled run, with reduced difference between calculated and experimental averaged void fraction. Furthermore, for another run with a higher equilibrium quality, voidfraction distribution with bubblesize models is also shown to come closer to experiment. To further investigate the behavior of the different bubblesize models, as a future work, it would be relevant to study the difference between the breakup and coalescence terms in the two models. On a longer term, to improve the quality of the simulations of saturated (high void fraction) cases, novel methods are being developed in the NEPTUNE project which consider the simulation of large bubbles with interface locating technics, while keeping a statistical treatment for the smaller ones [23].
Acknowledgments
This work has been achieved in the framework of the NEPTUNE project, financially supported by CEA (Commissariat à l’Énergie Atomique et aux Énergies Alternatives), EDF, IRSN (Institut de Radioprotection et de Sûreté Nucléaire), and AREVANP.
References
 PSBT page on the OECD website. July 2012, http://www.oecdnea.org/science/wprs/egrsltb/PSBT.
 A. Guelfi, D. Bestion, M. Boucker et al., “NEPTUNE: a new software platform for advanced nuclear thermal hydraulics,” Nuclear Science and Engineering, vol. 156, no. 3, pp. 281–324, 2007. View at: Google Scholar
 Nuclear Energy Agency—Committee on the Safety of Nuclear Installations, “Best practice guidelines for the use of CFD in nuclear reactor safety applications,” Technical Report NEA/CSNI, 2007. View at: Google Scholar
 K. Hori, Y. Akiyama, K. Miyazaki, T. Kurosu, and S. Sugiyama, “Void Fraction in a single channel simulating one subchannel of a PWR fuel assembly,” in TwoPhase Flow Modelling and Experimentation, G. P. Celata and R. K. Shah, Eds., vol. 1, ETS edition, 1995. View at: Google Scholar
 K. Hori, K. Miyazaki, T. Kurosu, S. Sugiyama, J. Matsumoto, and Y. Akiyama, “In bundle void fraction measurement of PWR fuel assembly,” in Proceedings of the 2nd ASME/JSME Nuclear Engineering Joint Conference, F. Peterson, Ed., vol. 1, March 1993. View at: Google Scholar
 P. Coste, J. Pouvreau, C. Morel, J. Laviéville, M. Boucker, and A. Martin, “Modeling turbulence and friction around a large interface in a threedimension twovelocity eulerian code,” in Proceedings of the 12th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH '12), Pittsburgh, Pa, USA, October 2007. View at: Google Scholar
 M. Ishii, ThermoFluid Dynamic, Theory of TwoPhase, Eyrolles, 1975.
 N. Méchitoua, M. Boucker, J. Laviéville, J. M. Hérard, S. Pigny, and G. Serre, “An unstructured finite volume solver for two phase water/vapour flows based on an elliptic oriented fractional step method,” in Proceedings of the 10th International Topical Meeting Nuclear Reactor Thermal Hydraulics (NURETH '10), Seoul, Republic of Korea, 2003. View at: Google Scholar
 S. Mimouni, F. Archambeau, M. Boucker, J. Lavieville, and C. Morel, “A second order turbulence model based on a Reynolds stress approach for twophase boiling flow. Part 1: Application to the ASUannular channel case,” Nuclear Engineering and Design, vol. 240, no. 9, pp. 2233–2243, 2010. View at: Publisher Site  Google Scholar
 M. Lance and M. L. de Bertodano, “Phase distribution phenomena and wall effects in bubbly twophase flows,” Multiphase Science and Technology, vol. 8, pp. 69–123, 1994. View at: Google Scholar
 M. Ishii and N. Zuber, “Drag coefficient and relative velocity in bubbly, droplet or particulate flows,” AIChE Journal, vol. 25, no. 5, pp. 843–855, 1979. View at: Google Scholar
 N. Zuber, “On the dispersed twophase flow in the laminar flow regime,” Chemical Engineering Science, vol. 19, no. 11, pp. 897–917, 1964. View at: Google Scholar
 A. Tomiyama, H. Tamai, I. Zun, and S. Hosokawa, “Transverse migration of single bubbles in simple shear flows,” Chemical Engineering Science, vol. 57, no. 11, pp. 1849–1858, 2002. View at: Publisher Site  Google Scholar
 N. Kurul and M. Podowski, “Multidimensional effects in forced convection subcooled boiling,” in Proceedings of the 9th International Heat Transfer Conference, vol. 1, p. 21, Jerusalem, Israel, 1990. View at: Google Scholar
 S. Mimouni, F. Archambeau, M. Boucker, J. Lavieville, and C. Morel, “A second order turbulence model based on a Reynolds stress approach for twophase boiling flow and application to fuel assembly analysis,” Nuclear Engineering and Design, vol. 240, no. 9, pp. 2225–2232, 2010. View at: Publisher Site  Google Scholar
 W. Yao and C. Morel, “Volumetric interfacial area prediction in upward bubbly twophase flow,” International Journal of Heat and Mass Transfer, vol. 47, no. 2, pp. 307–328, 2004. View at: Publisher Site  Google Scholar
 P. Ruyer, N. Seiler, M. Weiss, and F. P. Weiss, “A bubble size distribution model for the simulation of bubbly flows,” in Proceedings of the 6th International Conference on Multiphase Flows, Leipzig, Germany, 2007. View at: Google Scholar
 J. G. Collier and J. R. Thome, Convective Boiling and Condensation, Oxford Science, 3rd edition, 1994.
 N. Seiler and P. Ruyer, “Advanced model for polydispersion in size in boiling flows,” in 190ème Session Du Comité Scientifique Et Technique De La Société HydroTechnique De France : Modélisation Des Écoulements Diphasiques bouilLants, 2008. View at: Google Scholar
 T. Hibiki and M. Ishii, “Onegroup interfacial area transport of bubbly flows in vertical round tubes,” International Journal of Heat and Mass Transfer, vol. 43, no. 15, pp. 2711–2726, 2000. View at: Publisher Site  Google Scholar
 A. M. Kamp, Ecoulements turbulents à bulles dans une conduite en micropesanteur [Ph.D. thesis], Institut National Polytechnique de Toulouse, 1996.
 A. M. Kamp, A. K. Chesters, C. Colin, and J. Fabre, “Bubble coalescence in turbulent flows: a mechanistic model for turbulenceinduced coalescene applied to microgravity bubbly pipe flow,” International Journal of Multiphase Flow, vol. 27, no. 8, pp. 1363–1396, 2001. View at: Publisher Site  Google Scholar
 R. Denèfle, S. Mimouni, J. P. Caltagirone, and S. Vincent, “Multifield hybrid approach from bubble to slug flow transition modelling,” in Advances in Fluid Mechanics IX, M. Rahman and C. A. Brebbia, Eds., WIT, 2012. View at: Google Scholar
Copyright
Copyright © 2012 C. Baudry et al. 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.