Linear stability analysis of RELAP5 two-fluid model in nuclear reactor safety results

• Published in: Annals of nuclear energy 2020-12, Vol.149, p.107720, Article 107720
• Main Authors: Saraswat, S.P., Munshi, P., Allison, C.
• Format: Article
• Language: English
• Subjects: Dispersion analysis; Linear stability of two fluid model; Nuclear reactor safety analysis; Numerical stability; RELAP5
• ISSN: 0306-4549; 1873-2100
• DOI: 10.1016/j.anucene.2020.107720

•The linear stability of the RELAP5 TFM in the nuclear reactor accident conditions.•The numerical stability of the RELAP5 TFM in nuclear reactor accident conditions.•The effect of the addition of diffusion terms on the linear stability.•The effect of the addition of a bubble collision force on the linear stability.•Sensitivity analysis of RELAP5 TFM. System thermal-hydraulic code RELAP5 is based on a two-fluid, non-equilibrium, and non-homogeneous hydrodynamic model for simulation of transient two-phase behavior. The code model includes six governing equations to incorporate the mass, energy, and momentum of the two fluids. In this paper, linear stability analysis is performed to check the ill-posedness of the RELAP5 specific two-fluid model (TFM) for all normal and accident conditions of a standard pressurized water reactor (PWR). The analysis gives information about the soundness of the model and identifies the range of parameters where the solutions obtained from the model will be numerically convergent. The linear two-phase fluid dynamic stability (by dispersion analysis) of the RELAP5 one-dimensional two-fluid model and numerical stability of the difference equation formulation (by Von Neumann method) is presented. The present analysis shows that the two-fluid model becomes ill-posed for some fluid conditions, where the results are less accurate, so sensitivity analysis plays an important role. The ill-posed nature implies that results thus obtained (by finite difference method) have to be interpreted carefully because of the sensitive nature of reactor safety analysis. It is also identified that the variation in various parameters (like slip ratio, system pressure, void fraction, and phasic velocities) can affect the error growth rates. It has been demonstrated that the basic system of one-dimensional two-phase flow equations, that possesses complex characteristics, exhibits unbounded instabilities in the short-wavelength limit and constitutes an improperly posed initial value problem. The semi-implicit numerical method, which is unconditionally stable for hyperbolic systems, becomes unstable for non-hyperbolic systems. For some of the fluid conditions, even after the introduction of artificial viscosity terms (in the difference equation formulation) that damp the high-frequency spatial components of the solution, are not sufficient for regularization of the two-fluid model. Thus, there is a need for the addition of newer terms, e.g. bubble collision