Direct numerical simulation of internally heated natural convection in a hemispherical geometry

Research output: Contribution to journalArticlepeer-review

Standard Standard

Direct numerical simulation of internally heated natural convection in a hemispherical geometry. / Bian, Boshen; Dovizio, Daniele; Villanueva, Walter.
In: International Journal of Heat and Mass Transfer, Vol. 220, 124997, 01.03.2024.

Research output: Contribution to journalArticlepeer-review

HarvardHarvard

APA

CBE

MLA

VancouverVancouver

Bian B, Dovizio D, Villanueva W. Direct numerical simulation of internally heated natural convection in a hemispherical geometry. International Journal of Heat and Mass Transfer. 2024 Mar 1;220:124997. Epub 2023 Nov 30. doi: 10.1016/j.ijheatmasstransfer.2023.124997

Author

Bian, Boshen ; Dovizio, Daniele ; Villanueva, Walter. / Direct numerical simulation of internally heated natural convection in a hemispherical geometry. In: International Journal of Heat and Mass Transfer. 2024 ; Vol. 220.

RIS

TY - JOUR

T1 - Direct numerical simulation of internally heated natural convection in a hemispherical geometry

AU - Bian, Boshen

AU - Dovizio, Daniele

AU - Villanueva, Walter

PY - 2024/3/1

Y1 - 2024/3/1

N2 - Internally heated (IH) natural convection can be found in nature, industrial processes, or during a severe accident in a light water reactor. In this accident scenario, the nuclear reactor core and some internal structures can melt down and relocate to the lower head of the reactor pressure vessel (RPV) and interact with the remaining coolant. Subsequent re-heating and re-melting under decay and oxidation heat creates a transition from a debris bed to a molten pool. The molten pool, which can involve more than hundred tons of dangerously superheated oxidic and metallic liquids, imposes thermo-mechanical loads on the vessel wall that can lead to a thermal and/or structural failure of the vessel and subsequent release of radioactive materials to the reactor pit, and can possibly make its way to the environment. This study uses Direct Numerical Simulation (DNS) to investigate homogeneous IH molten pool convection in a hemispherical domain using Nek5000, an open-source spectral element code. With a Rayleigh number of 1.6 × 10^11, the highest reached through DNS in this confined hemispherical geometry, and a Prandtl number of 0.5, which corresponds to a prototypic corium, the study provides detailed information on the thermo-fluid behavior. The results show a turbulent flow with three distinct regions, consistent with the general flow observations from the BALI experiments. The study also presents detailed information on turbulence, such as turbulent kinetic energy (TKE), turbulent heat flux (THF), and temperature variance. Additionally, the study provides 3D heat flux distributions along the boundaries. The heat fluxes along the top boundary fluctuate due to the turbulent eddies in the vicinity, while along the curved boundary the heat fluxes increase nonlinearly from the bottom to the top.

AB - Internally heated (IH) natural convection can be found in nature, industrial processes, or during a severe accident in a light water reactor. In this accident scenario, the nuclear reactor core and some internal structures can melt down and relocate to the lower head of the reactor pressure vessel (RPV) and interact with the remaining coolant. Subsequent re-heating and re-melting under decay and oxidation heat creates a transition from a debris bed to a molten pool. The molten pool, which can involve more than hundred tons of dangerously superheated oxidic and metallic liquids, imposes thermo-mechanical loads on the vessel wall that can lead to a thermal and/or structural failure of the vessel and subsequent release of radioactive materials to the reactor pit, and can possibly make its way to the environment. This study uses Direct Numerical Simulation (DNS) to investigate homogeneous IH molten pool convection in a hemispherical domain using Nek5000, an open-source spectral element code. With a Rayleigh number of 1.6 × 10^11, the highest reached through DNS in this confined hemispherical geometry, and a Prandtl number of 0.5, which corresponds to a prototypic corium, the study provides detailed information on the thermo-fluid behavior. The results show a turbulent flow with three distinct regions, consistent with the general flow observations from the BALI experiments. The study also presents detailed information on turbulence, such as turbulent kinetic energy (TKE), turbulent heat flux (THF), and temperature variance. Additionally, the study provides 3D heat flux distributions along the boundaries. The heat fluxes along the top boundary fluctuate due to the turbulent eddies in the vicinity, while along the curved boundary the heat fluxes increase nonlinearly from the bottom to the top.

U2 - 10.1016/j.ijheatmasstransfer.2023.124997

DO - 10.1016/j.ijheatmasstransfer.2023.124997

M3 - Article

VL - 220

JO - International Journal of Heat and Mass Transfer

JF - International Journal of Heat and Mass Transfer

SN - 1879-2189

M1 - 124997

ER -