Modeling of flashing-induced flow instabilities for a natural circulation driven novel modular reactor

An analytical study based on frequency domain analysis is presented on the flashing-induced flow instability in a natural circulation test facility, which was designed to investigate the flow instability for a BWR-type novel modular reactor (NMR). To address the flashing phenomena at low pressure conditions, such as initial startup transients or accidents, the liquid enthalpy change in the P-T diagram due to reduced hydrostatic head in the riser or chimney was treated as an axially uniform heat flux. Based on the drift flux model, the system transfer function was obtained through small perturbations about the steady state in the frequency domain. The D-partition method was used to determine the neutral stability boundary in the dimensionless stability plane, which was constituted of the subcooling number and phase change number. From the frequency domain analysis, the flashing stability boundary and the density wave oscillations boundary could be predicted. Some parametric studies had been performed on the system pressure and the inlet flow resistance coefficients in the stability analysis. The results showed that the flashing stability boundary was more sensitive to the system pressure than the density wave oscillations. In addition, the theoretical stability boundaries were benchmarked against the experimental stability boundaries from quasi-steady state tests. Although the general stability boundary agreed well with the experiments, certain discrepancies still existed due to the assumptions of thermal equilibrium in current study. In the future, the thermal non-equilibrium conditions including subcooled boiling will be taken into account in the flashing induced stability analysis.