DOI QR코드

DOI QR Code

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)
  • Received : 2010.10.29
  • Published : 2010.12.31

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

References

  1. IAEA. Renaissance is in motion. 2010 [cited; Available from:http://www.world-nuclear-news.org/nn_laea_renaissance_is_in_motion_0609102.html.
  2. 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.
  3. 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.
  4. Smith, B., Lawrence, Assessment of CFD codes used in nuclear reactor safety simulations. Nuclear Engineering and Technology, 2010. 42(4): p. 339 - 364. https://doi.org/10.5516/NET.2010.42.4.339
  5. 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.
  6. 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. https://doi.org/10.1016/j.nucengdes.2007.02.046
  7. 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.
  8. 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. https://doi.org/10.1016/j.nucengdes.2008.12.007
  9. 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.
  10. 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.
  11. Guelfi, A., et al., NEPTUNE: A new software platform for advanced nuclear thermal hydraulics. Nuclear Science and Engineering, 2007. 156: p. 281 - 324.
  12. 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.
  13. Koncar, B. and I. Tiselj, Influence of near-wall modelling on boiling flow simulation. Nuclear Engineering and Design, 2010. 240(2): p. 275-283. https://doi.org/10.1016/j.nucengdes.2008.07.019
  14. 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.
  15. 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. https://doi.org/10.1615/MultScienTechn.v18.i1.10
  16. 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. https://doi.org/10.1016/0021-9991(81)90145-5
  17. 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. https://doi.org/10.1006/jcph.1994.1155
  18. Peng, D., et al., A PDE-basedfast local level set method. Journal of Computational Physics, 1999. 155(2): p. 410-438. https://doi.org/10.1006/jcph.1999.6345
  19. Enright, D., et al., A hybrid particle level set method for improved interface capturing. Journal of Computational Physics, 2002. 183(1): p. 83-116. https://doi.org/10.1006/jcph.2002.7166
  20. 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.
  21. 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. https://doi.org/10.1006/jcph.2000.6537
  22. 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. https://doi.org/10.1016/S0021-9991(03)00087-1
  23. 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. https://doi.org/10.1006/jcph.2001.6888
  24. 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. https://doi.org/10.1016/0021-9991(87)90187-2
  25. 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. https://doi.org/10.1016/0010-4655(91)90071-R
  26. Gibbs, J.W., On the equilibrium of heterogeneous substances. Connecticut Academy Transactions, 1877. 3: p. 343.
  27. 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. https://doi.org/10.1063/1.1744102
  28. Ginzburg, V.L. and L.D. Landau, On the theory of superconductivity. Zhurnal Eksperimentalnoye i Teoreticheskoye Fiziki, 1950. 20: p. 1064.
  29. 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. https://doi.org/10.1115/1.2822535
  30. Chorin, A.J., Numerical Solution of Navier-Stokes Equations. Mathematics of Computation, 1968. 22(104): p. 745 - 762. https://doi.org/10.1090/S0025-5718-1968-0242392-2
  31. Versteeg, H.K. and W. Malalasekera, An introduction to computational fluid dynamics The finite volume method. 2nd ed. 2007: Prenticel Hall.
  32. 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. https://doi.org/10.1080/10407788608913463
  33. Kuhn, S., Transport mechanisms in mixed convective flow over complex surfaces. 2008, Swiss Federal Institute of Technology (ETH): Zurich.
  34. 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. https://doi.org/10.1007/s001620050083
  35. Henn, D.S. and R.I. Sykes, Large-eddy simulation of flow over wavy surfaces. 1999, Cambridge Journals Online. p. 75-112.
  36. 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. https://doi.org/10.1006/jcph.2001.6746
  37. 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. https://doi.org/10.1016/S0010-4655(99)00473-7
  38. Yabe, T., et al., An exactly conservative semi-lagrangian scheme (CIP-CSL) in one dimension. Monthly Weather Review, 2001. 129(2): p. 332-344. https://doi.org/10.1175/1520-0493(2001)129<0332:AECSLS>2.0.CO;2
  39. 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. https://doi.org/10.1006/jcph.2000.6625
  40. 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. https://doi.org/10.1016/j.jcp.2009.04.003
  41. 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. https://doi.org/10.1016/S0377-0427(02)00535-6
  42. 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. https://doi.org/10.1016/0010-4655(95)00124-7
  43. 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. https://doi.org/10.1016/0021-9991(92)90240-Y
  44. 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.
  45. 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.
  46. 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. https://doi.org/10.1103/PhysRevE.57.4323
  47. Chen, L.-Q., Phase-field models for microstructure evolution. Annual Review of Materials Research, 2002. 32(1): p. 113-140. https://doi.org/10.1146/annurev.matsci.32.112001.132041
  48. Stephan, P. and J. Hammer, A new model for nucleate boiling heat transfer. Warme- und Stoffubertragung, 1994. 30: p. 119–125.