CFD analysis and experimental validation of steady state mixed convection sodium flow

Abstract Formation and destruction of thermal stratification can occur under certain flow conditions in the upper plenum of sodium cooled fast breeder reactors (SFR). The flow patterns in the hot sodium pool of SFR upper plenum are very complex and include zones of free and wall-bounded jets, recirculation and stagnation areas. The interaction of sodium flow and thermal stratification has been analyzed experimentally at CEA in the SUPERCAVNA facility. The facility consists of a rectangular cavity with heated side walls. The flow is driven by a wall-bounded cold jet at the bottom of the cavity. Experimental data of the temperature distribution inside the cavity are available for steady-state and transient flow conditions. In the present study, the steady state experiments are analyzed with the CEA CFD reference code TrioCFD and the commercial code FLUENT employing Reynolds Averaged Navier-Stokes (RANS) equations. The SUPERCAVNA modeling is based on preliminary separate effects’ studies of wall-bounded jets and free convection benchmark simulations for low Prandtl number fluids. It was found that a two-dimensional treatment is sufficient to reproduce correctly the measured thermal stratification for steady-state SUPERCAVNA experiments. However, it is necessary to take into account conjugated heat transfer between cavity sodium flow and side-wall heating channel. Using temperature-dependent physical properties was also found to be an important factor in simulating the experiments correctly. Applying the Boussinesq approximation to study the impact of buoyancy on the vertical flow momentum was found to be justified. Turbulence modelling is necessary for a successful simulation of the experiments. For practical reasons, the analysis was restricted to turbulence models of the k-e family. High Re k-e models, in either standard, realizable or RNG formulations, do not lead to significant differences in the calculated temperature fields. Non-linear eddy viscosity modelling might improve the quality of the simulation.

[1]  T. Shih,et al.  A New K-epsilon Eddy Viscosity Model for High Reynolds Number Turbulent Flows: Model Development and Validation , 1994 .

[2]  F. Leroux,et al.  Interaction of recirculation and stable stratification in a rectangular cavity filled with sodium , 1988 .

[3]  Ashraf Balabel,et al.  On the performance of linear and nonlinear k–ε turbulence models in various jet flow applications , 2011 .

[4]  Masud Behnia,et al.  Fine Mesh Solutions Using Stream Function-Vorticity Formulation , 1990 .

[5]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[6]  G. Grötzbach,et al.  Challenges in low-Prandtl number heat transfer simulation and modelling , 2013 .

[7]  R. Stieglitz,et al.  Flow measurement techniques in heavy liquid metals , 2010 .

[8]  H. Ninokata,et al.  A turbulence model study for simulating flow inside tight lattice rod bundles , 2005 .

[9]  Martin Wosnik,et al.  A similarity theory for the turbulent plane wall jet without external stream , 2000, Journal of Fluid Mechanics.

[10]  John Derek Jackson,et al.  Turbulent mixed convection heat transfer to liquid sodium , 1983 .

[11]  D. Tenchine,et al.  Some thermal hydraulic challenges in sodium cooled fast reactors , 2010 .

[12]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[13]  Wolfgang Rodi,et al.  The Turbulent Wall Jet Measurements and Modeling , 1983 .

[14]  S. Orszag,et al.  Development of turbulence models for shear flows by a double expansion technique , 1992 .

[15]  M. Ye,et al.  Investigation on the applicability of turbulent-Prandtl-number models in bare rod bundles for heavy liquid metals , 2017 .

[16]  Emilio Baglietto,et al.  status and perspective of turbulence hat transfer modelling for the industrial application of liquid metal flows , 2015 .

[17]  Johan Persson,et al.  An experimental study of a two-dimensional plane turbulent wall jet , 1998 .

[18]  W. Kays,et al.  Turbulent Prandtl number : where are we ? , 1994 .

[19]  Philippe Emonot.  Méthodes de volumes éléments finis : applications aux équations de Navier Stokes et résultats de convergence , 1992 .

[20]  Jean-Luc Guermond,et al.  International Journal for Numerical Methods in Fluids on Stability and Convergence of Projection Methods Based on Pressure Poisson Equation , 2022 .

[21]  T. Wetzel,et al.  Numerical investigation of turbulent aided mixed convection of liquid metal flow through a concentric annulus , 2017 .