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

• Published in: Nuclear engineering and design Vol. 326; pp. 333 - 343
• Main Authors: Bieder, U., Ziskind, G., Rashkovan, A.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.01.2018; Elsevier BV; Elsevier
• Subjects: Boussinesq approximation; Breeder reactors; Computational fluid dynamics; Computer simulation; Convection; Eddy viscosity; Equilibrium flow; Fast breeder; Fluid flow; Formulations; Free convection; Heat transfer; Low Prandtl number fluid; Mathematical analysis; Mathematical models; Mixed convection; Modelling; Navier-Stokes equations; Nuclear Experiment; Nuclear Theory; Physical properties; Physics; Prandtl number; Reynolds averaged Navier-Stokes method; Sodium; Sodium cooled reactors; Stagnation; Steady state; Studies; Temperature dependence; Temperature distribution; Temperature effects; Thermal stratification; Turbulence; Turbulence models; Unsteady flow; Vertical flow; Viscosity; Wall bounded jet; Low Prandtl number fluid; Wall bounded jet; Sodium; Thermal stratification; Mixed convection
• ISSN: 0029-5493; 1872-759X
• DOI: 10.1016/j.nucengdes.2017.11.028

•Mixed convection liquid metal flow is analyzed by CFD.•Conjugated heat transfer modelling between cavity and heating channel.•Modeling strategy is determined by detailed pre-calculations.•Comparison of calculations to experiments.•Successful simulations with RANS type (k-ε family) turbulence models. 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-ε family. High Re k-ε 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.