Heat transfer coefficients model for SIMMER-III and SIMMER-IV

Brear, D. J.*; Kondo, Satoru; Sogabe, Joji; Tobita, Yoshiharu*; Kamiyama, Kenji

JAEA-Research 2024-009, 134 Pages, 2024/10

JAEA-Research-2024-009.pdf:2.45MB

The SIMMER-III/SIMMER-IV computer codes are being used for liquid-metal fast reactor (LMFR) core disruptive accident (CDA) analysis. The sequence of events predicted in a CDA is often influenced by the heat exchanges between LMFR materials, which are controlled by heat transfer coefficients (HTCs) in the respective materials. The mass transfer processes of melting and freezing, and vaporization and condensation are also controlled by HTCs. The complexities in determining HTCs in a multi-component and multi-phase system are the number of HTCs to be defined at binary contact areas of a fluid with other fluids and structure surfaces, and the modes of heat transfer taking into account different flow topologies representing flow regimes with and without structure. As a result, dozens of HTCs are evaluated in each mesh cell for the heat and mass transfer calculations. This report describes the role of HTCs in SIMMER-III/SIMMER-IV, the heat transfer correlations implemented and the calculation of HTCs in all topologies in multi-component, multi-phase flows. A complete description of the physical basis of HTCs and available experimental correlations is contained in Appendices to this report. The major achievement of the code assessment program conducted in parallel with code development is summarized with respect to HTC modeling to demonstrate that the coding is reliable and that the model is applicable to various multi-phase problems with and without reactor materials.

Effect of gas density and surface tension on liquid film thickness in vertical upward disturbance wave flow

Zhang, H.*; Mori, Shoji*; Hisano, Tsutomu*; Yoshida, Hiroyuki

International Journal of Multiphase Flow, 159, p.104342_1 - 104342_15, 2023/02

https://doi.org/10.1016/j.ijmultiphaseflow.2022.104342

Lattice Boltzmann modeling and simulation of forced-convection boiling on a cylinder

Saito, Shimpei*; De Rosis, A.*; Fei, L.*; Luo, K. H.*; Ebihara, Kenichi; Kaneko, Akiko*; Abe, Yutaka*

Physics of Fluids, 33(2), p.023307_1 - 023307_21, 2021/02

https://doi.org/10.1063/5.0032743

A Boiling phenomenon in a liquid flow field is known as forced-convection boiling. We numerically investigated the boiling system on a cylinder in a flow at a saturated condition. To deal with such a phenomenon, we developed a numerical scheme based on the pseudopotential lattice Boltzmann method. The collision was performed in the space of central moments (CMs) to enhance stability for high Reynolds numbers. Furthermore, additional terms for thermodynamic consistency were derived in a CMs framework. The effectiveness of the model was tested against some boiling processes, including nucleation, growth, and departure of a vapor bubble for high Reynolds numbers. Our model can reproduce all the boiling regimes without the artificial initial vapor phase. We found that the Nukiyama curve appears even though the focused system is the forced-convection system. Also, our simulations support experimental observations of intermittent direct solid-liquid contact even in the film-boiling regime.

Development of unstructured mesh-based numerical method for sodium-water reaction phenomenon in steam generators of sodium-cooled fast reactors

Uchibori, Akihiro; Watanabe, Akira*; Takata, Takashi; Ohshima, Hiroyuki

Journal of Nuclear Science and Technology, 54(10), p.1036 - 1045, 2017/10

https://doi.org/10.1080/00223131.2017.1344155

To evaluate a sodium-water reaction phenomenon in a steam generator of sodium-cooled fast reactors, a computational fluid dynamics code SERAPHIM, in which a compressible multicomponent multiphase flow with sodium-water chemical reaction is computed, has been developed. The original SERAPHIM code is based on the difference method. In this study, unstructured mesh-based numerical method was developed to advance a numerical accuracy for the complex-shaped domain including multiple heat transfer tubes. Numerical analysis of an underexpanded jet experiment was performed as part of validation of the unstructured mesh-based numerical method. The calculated pressure profile showed good agreement with the experimental data. Applicability of the numerical method for the actual situation was confirmed through the analysis of water vapor discharging into liquid sodium. The effect of use of the unstructured mesh was also investigated by the two analyses using structured and unstructured mesh.

Development of prediction technology of two-phase flow dynamics under earthquake acceleration, 16; Experimental and numerical study of pressure fluctuation effects on bubble motion

Kato, Yuki; Yoshida, Hiroyuki; Yokoyama, Ryotaro*; Kanagawa, Tetsuya*; Kaneko, Akiko*; Monji, Hideaki*; Abe, Yutaka*

Proceedings of 23rd International Conference on Nuclear Engineering (ICONE-23) (DVD-ROM), 8 Pages, 2015/05

Development of prediction technology of two-phase flow dynamics under earthquake acceleration

Yoshida, Hiroyuki; Nagatake, Taku; Takase, Kazuyuki; Kaneko, Akiko*; Monji, Hideaki*; Abe, Yutaka*

Mechanical Engineering Journal (Internet), 1(4), p.TEP0025_1 - TEP0025_11, 2014/08

https://doi.org/10.1299/mej.2014tep0025

Numerical analysis of three-dimensional two-phase flow behavior in a fuel assembly

Takase, Kazuyuki; Yoshida, Hiroyuki; Ose, Yasuo*; Akimoto, Hajime

WIT Transactions on Engineering Sciences, Vol.50, p.183 - 192, 2005/00

no abstracts in English

Numerical simulation of the interfacial growth of the stratified wa vytwo-phase flow in the horizontal rectangular channel

Ebihara, Kenichi; Watanabe, Tadashi

Proceedings of 2004 ASME International Mechanical Engineering Congress and Exposition (IMECE '04) (CD-ROM), 9 Pages, 2004/11

This paper describes the numerical simulation of the interfacial growth of the stratified wavy two-phase flow in the horizontal rectangular channel. The influence of the channel width upon the growth of the interfacial wave was evaluated by the several simulations for the different channel width. The numerical simulation model adopted in this paper is a one-component two-phase fluid model of the lattice Boltzmann method. The wave growth was observed and the dimensionless numbers that characterize the two-phase flow state were measured during the computations. The relation between the wave growth and the dimensionless numbers, which characterize the flow state, was compared with that in the flow pattern map proposed by Taitel and Dukler. It was verified in the case of the wide channel width that the simulated relation was almost in agreement with that in the flow pattern map. It was shown that the narrower the channel width became, the more mass flow rate of the rare phase the interfacial growth needed and the obtained relation deviated from that in the flow pattern map.

Direct numerical simulation of turbulent heat transfer in plane impinging jet; Effects of impingement distance on heat transfer in confined space

Hattori, Hirofumi*; Sato, Hiroshi; Nagano, Yasutaka*

Nihon Kikai Gakkai Rombunshu, B, 70(696), p.1919 - 1926, 2004/08

https://doi.org/10.1299/kikaib.70.1919

no abstracts in English

Numerical study of zonal flow dynamics and electron transport in electron temperature gradient driven turbulence

Li, J.; Kishimoto, Yasuaki

Physics of Plasmas, 11(4), p.1493 - 1510, 2004/04

https://doi.org/10.1063/1.1669397

The electron temperature gradient (ETG) driven turbulence in tokamak core plasmas is numerically investigated based on three-dimensional gyrofluid model with adiabatic ion response. Attentions are focused on the zonal flow dynamics in ETG fluctuations and the resultant electron heat transport. A high electron energy confinement mode is found in the weak magnetic shear regime, which is closely relevant with self-organization behavior of turbulence through the enhanced zonal flow dynamics rather than the weak shear stabilization of ETG fluctuations. It is demonstrated that the weak shear is favorable for the enhancement of zonal flows in ETG turbulence.