Numerical study of air ingress transition to natural circulation in a high temperature helium loop

The generation-IV high temperature gas cooled reactors (HTGRs) are designed with many passive safety features, one of which is the ability to passively remove heat under a loss of coolant accident (LOCA). However, several common reactor designs do not prevent against a large break in the coolant system and may therefore experience a depressurized LOCA. This would lead to air entering into the reactor system via several potential modes of ingress: diffusion, gravity currents, and natural circulation. At the onset of a LOCA, the initial rate of air ingress is expected to be very slow because it is governed by molecular diffusion. However, after several hours, natural circulation would commence, thus, bringing the air into the reactor system at a much higher rate. As a consequence, air ingress would cause the high temperature graphite matrix to oxidize, leading to its thermal degradation and decreased passive heat (decay) removal capability. Therefore, it is essential to understand the transition of air ingress from molecular diffusion to natural circulation in an HTGR system. This paper presents results from a computational fluid dynamics (CFD) model to study the air ingress transition behavior. These results are validated against an h-shaped high temperature helium loop experiment (McIlroy et al., 2010). Details are provided to quantitatively predict the transition time from molecular diffusion to natural circulation.