Browse > Article
http://dx.doi.org/10.5516/NET.2010.42.6.620

MULTI-SCALE MODELING AND ANALYSIS OF CONVECTIVE BOILING: TOWARDS THE PREDICTION OF CHF IN ROD BUNDLES  

Niceno, B. (Laboratory for Thermal Hydraulics, Nuclear Energy and Safety Department Paul Scherrer Institut)
Sato, Y. (Laboratory for Thermal Hydraulics, Nuclear Energy and Safety Department Paul Scherrer Institut)
Badillo, A. (Laboratory for Thermal Hydraulics, Nuclear Energy and Safety Department Paul Scherrer Institut)
Andreani, M. (Laboratory for Thermal Hydraulics, Nuclear Energy and Safety Department Paul Scherrer Institut)
Publication Information
Nuclear Engineering and Technology / v.42, no.6, 2010 , pp. 620-635 More about this Journal
Abstract
In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.
Keywords
Multi-Scale Modeling; Convective Boiling; Critical Heat Flux; Immersed Boundary Method; CIP-CSL2; Phase Field;
Citations & Related Records
Times Cited By KSCI : 1  (Citation Analysis)
Times Cited By Web Of Science : 1  (Related Records In Web of Science)
Times Cited By SCOPUS : 3
연도 인용수 순위
1 Yabe, T. and T. Aoki, A universal solver for hyperbolic equations by cubic-polynomial interpolation I. Onedimensional solver. Computer Physics Communications, 1991. 66(2-3): p. 219-232.   DOI
2 Gibbs, J.W., On the equilibrium of heterogeneous substances. Connecticut Academy Transactions, 1877. 3: p. 343.
3 Cahn, J.W. and J.E. Hilliard, Free energy of a nonuniform system. I. interfacial free energy. Journal of Chemical Physics, 1958. 28(2): p. 258 - 267.   DOI
4 Ginzburg, V.L. and L.D. Landau, On the theory of superconductivity. Zhurnal Eksperimentalnoye i Teoreticheskoye Fiziki, 1950. 20: p. 1064.
5 Bestion, D., et al., Review of available data for validation of NURESIM two-phase CFD software applied to CHF investigations. Science and Technology of Nuclear Installations, Article ID 214512, 2009: p. 14.
6 Guelfi, A., et al., NEPTUNE: A new software platform for advanced nuclear thermal hydraulics. Nuclear Science and Engineering, 2007. 156: p. 281 - 324.
7 Jamet, D., et al. Towards a multi-scale approach of twophase flow modeling in the context of DNB modeling. in OECD/NEA Workshop on Experiments and CFD code Application to Nuclear Reactor Safety (XCFD4NRS). 2008. Grenoble, France.
8 Koncar, B. and I. Tiselj, Influence of near-wall modelling on boiling flow simulation. Nuclear Engineering and Design, 2010. 240(2): p. 275-283.   DOI
9 Ahn, H.S., et al., Pool boiling CHF enhancement by micro/nanoscale modification of zircaloy-4 surface. Nuclear Engineering and Design. In Press, Corrected Proof.
10 Theofanous, T.G. and T.N. Dinh, High heat flux boiling and burnout as microphysical phenomena: mounting evidence and opportunities. Multiphase Science and Technology, 2006. 18(1): p. 1 - 26.   DOI
11 Hirt, C.W. and B.D. Nichols, Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of Computational Physics, 1981. 39(1): p. 201-225.   DOI   ScienceOn
12 Sussman, M., P. Smereka, and S. Osher, A Level set approach for computing solutions to incompressible twophase flow. Journal of Computational Physics, 1994. 114(1): p. 146-159.   DOI
13 Peng, D., et al., A PDE-basedfast local level set method. Journal of Computational Physics, 1999. 155(2): p. 410-438.   DOI
14 Enright, D., et al., A hybrid particle level set method for improved interface capturing. Journal of Computational Physics, 2002. 183(1): p. 83-116.   DOI
15 IAEA. Renaissance is in motion. 2010 [cited; Available from:http://www.world-nuclear-news.org/nn_laea_renaissance_is_in_motion_0609102.html.
16 Stosic, Z.V. The nuclear renaissance and AREVA’s reactor designs for the 21st century: EPR and SWR-1000. in International Conference Nuclear Energy for New Europe. 2007. PortoroZ, Slovenia.
17 Guelfi, A., et al., Advanced boiling two-phase flow modeling: A major challenge for nuclear industry”, in AIRH-SHF-SFT Colloquium on two-phase boiling flow. 2008: Grenoble, France.
18 Smith, B., Lawrence, Assessment of CFD codes used in nuclear reactor safety simulations. Nuclear Engineering and Technology, 2010. 42(4): p. 339 - 364.   DOI
19 Frank, T., C. Lifante, and A.D. Burns, Recent advances in modeling and simulation of boiling processes in ANSYS CFD, in ANSYS Conf. and 27th CADFEM Users Meeting. 2009: Leipzig, Germany.
20 Ustinenko, V., et al., Validation of CFD-BWR, a new twophase computational fluid dynamics model for boiling water reactor analysis. Nuclear Engineering and Design, 2008. 238(3): p. 660-670.   DOI   ScienceOn
21 Wehle, F., et al., Progress in AREVA NP's LWR thermal hydraulics methodology, in Workshop on Modeling and Measurements of Two-Phase Flows and Heat Transfer in Nuclear Fuel Assemblies. 2006: Stockholm, Sweeden.
22 Le Corre, J.-M., S.-C. Yao, and C.H. Amon, A mechanistic model of critical heat flux under subcooled flow boiling conditions for application to one- and three-dimensional computer codes. Nuclear Engineering and Design, 2010. 240(2): p. 235-244.   DOI
23 Podowski, M.Z. and R.M. Podowski, Mechanistic multidimensional modeling of forced convection boiling heat transfer. Science and Technology of Nuclear Installations, Article ID 387020, 2009: p. 10.
24 Toda, K., Y. Ogata, and T. Yabe, Multi-dimensional conservative semi-Lagrangian method of characteristics CIP for the shallow water equations. Journal of Computational Physics, 2009. 228(13): p. 4917-4944.   DOI
25 Yabe, T., et al., The next generation CIP as a conservative semi-Lagrangian solver for solid, liquid and gas. Journal of Computational and Applied Mathematics, 2002. 149(1): p. 267-277.   DOI
26 Xiao, F., T. Yabe, and T. Ito, Constructing oscillation preventing scheme for advection equation by rational function. Computer Physics Communications, 1996. 93(1): p. 1-12.   DOI
27 Brackbill, J.U., D.B. Kothe, and C. Zemach, A continuum method for modeling surface tension. Journal of Computational Physics, 1992. 100(2): p. 335-354.   DOI   ScienceOn
28 Brereton, G. and D. Korotney, Coaxial and oblique coalescence of two rising bubbles. In Dynamics of Bubbles and Vortics Near a Free Surface. Vol. 119. 1991, New York: ASME.
29 Shi, B., S. Sinha, and V. Dhir, K, Molecular dynamics simulation of the density and surface tension of water by particle-particle particle-mesh method. Journal of Chemical Physics, 2006. 124(204715): p. 7.
30 Karma, A. and W.-J. Rappel, Quantitative phase-field modeling of dendritic growth in two and three dimensions. Physical Review E, 1998. 57(4): p. 4323.   DOI
31 Chen, L.-Q., Phase-field models for microstructure evolution. Annual Review of Materials Research, 2002. 32(1): p. 113-140.   DOI
32 Stephan, P. and J. Hammer, A new model for nucleate boiling heat transfer. Warme- und Stoffubertragung, 1994. 30: p. 119–125.
33 Lay, J.H. and V.K. Dhir, Shape of a vapor stem during nucleate boiling of saturated liquids. Journal of Heat Transfer 1995. 117(2): p. 394 - 401.   DOI   ScienceOn
34 Chorin, A.J., Numerical Solution of Navier-Stokes Equations. Mathematics of Computation, 1968. 22(104): p. 745 - 762.   DOI   ScienceOn
35 Versteeg, H.K. and W. Malalasekera, An introduction to computational fluid dynamics The finite volume method. 2nd ed. 2007: Prenticel Hall.
36 Hutchinson, B.R. and G.D. Raithby, A multigrid method based on the additive correction strategy. Numerical Heat Transfer, Part B, 1986. 9: p. 27.   DOI
37 Kuhn, S., Transport mechanisms in mixed convective flow over complex surfaces. 2008, Swiss Federal Institute of Technology (ETH): Zurich.
38 Cherukat, P., et al., Direct Numerical Simulation of a Fully Developed Turbulent Flow over a Wavy Wall. Theoretical and Computational Fluid Dynamics, 1998. 11(2): p. 109-134.   DOI
39 Henn, D.S. and R.I. Sykes, Large-eddy simulation of flow over wavy surfaces. 1999, Cambridge Journals Online. p. 75-112.
40 Xiao, F. and T. Yabe, Completely conservative and oscillationless semi-lagrangian schemes for advection Transportation. Journal of Computational Physics, 2001. 170(2): p. 498-522.   DOI
41 Tanaka, R., T. Nakamura, and T. Yabe, Constructing exactly conservative scheme in a non-conservative form. Computer Physics Communications, 2000. 126(3): p. 232-243.   DOI
42 Yabe, T., et al., An exactly conservative semi-lagrangian scheme (CIP-CSL) in one dimension. Monthly Weather Review, 2001. 129(2): p. 332-344.   DOI
43 Yabe, T., F. Xiao, and T. Utsumi, The Constrained Interpolation Profile Method for Multiphase Analysis. Journal of Computational Physics, 2001. 169(2): p. 556-593.   DOI
44 Youngs, D.L., Time-dependent multi-material flow with large fluid distortion. In Numerical Methods for Fluid Dynamics. Numerical Methods for Fluid Dynamics, ed. K.W. Morton and M.J. Baines. 1982, New York: Academic Press. 273-285.
45 Sussman, M. and E.G. Puckett, A coupled level set and volume-of-fluid method for computing 3D and axisymmetric incompressible two-phase flows. Journal of Computational Physics, 2000. 162(2): p. 301-337.   DOI
46 Sussman, M., A second order coupled level set and volumeof-fluid method for computing growth and collapse of vapor bubbles. Journal of Computational Physics, 2003. 187(1): p. 110-136.   DOI
47 Nakamura, T., et al., Exactly conservative semi-Lagrangian scheme for multi-dimensional hyperbolic equations with directional splitting technique. Journal of Computational Physics, 2001. 174(1): p. 171-207.   DOI
48 Takewaki, H. and T. Yabe, The cubic-interpolated pseudo particle (CIP) method: application to nonlinear and multidimensional hyperbolic equations. Journal of Computational Physics, 1987. 70(2): p. 355-372.   DOI