LES investigations of turbulent heat transfer structures with exponential power escalation

Abstract

The BORAX-type accident is a key scenario in the safety design of pool-type experimental nuclear reactors. It considers a reactivity insertion large enough to initiate an exponential power escalation with a period as small as tens of milliseconds, necessitating effective heat removal by the coolant system. This study presents a computational analysis of single-phase transient heat transfer under thermal-hydraulic conditions relevant to such reactors, focusing on turbulent channel flow between planar fuel plates with highly subcooled water. Our investigation uses Large Eddy Simulation (LES) performed with TrioCFD, the open-source CFD code developed by the French Atomic Energy Commission (CEA). Several modeling options are tested: incompressible flow with temperature as a passive scalar, incompressible flow with varying viscosity, and quasicompressible flow with temperature-dependent viscosity and mass density. To evaluate the accuracy of our physical models, all simulations were performed under uniform conditions, using the same experimental parameters, mesh, and numerical strategies. The specific examined scenario involves a highly subcooled turbulent channel flow at moderate pressure, experiencing an exponential power increase. We ensure detailed turbulence resolution at the Batchelor scale in the wall normal direction close to the heating wall. This communication presents a comparative study of LES simulations of the single-phase flow under exponential power excursion at the wall. Subsequent analyses focused on the impact of different velocity-temperature couplings on the prediction of wall heat transfer, incorporating comparisons with experimental infrared thermometry measurements. Early in the transient phase, temperature peaks are confined within turbulent streaks, as reported in the existing literature. However, as the exponential power transient unfolds, the viscous layer and streak structures destabilize, leading to a more dispersed distribution of hot spots. Preliminary findings suggest a transition in heat transfer modes from turbulence-driven to conduction during the final stages of the single-phase transient.