Improvement of RELAP5 code for density wave instability analysis in parallel narrow rectangular channels. (December 2018)

• Record Type: Journal Article
• Authors: Lian, Qiang; Tian, Wenxi; Qiu, Suizheng; Su, Guanghui
• Is Part Of: Annals of nuclear energy. Volume 122(2018)
• Page Start: 241
• Page End: 255
• Publication Date: 2018-12
• Subjects: Density wave instability -- RELAP5 modification -- Narrow rectangular channel -- Stability boundary; Nuclear energy -- Periodicals; Nuclear engineering -- Periodicals; 621.4805
• DOI: 10.1016/j.anucene.2018.08.049
• Languages: English
• ISSNs: 0306-4549

Abstract

Highlights: RELAP5 is modified to predict DWO in parallel narrow rectangular channels. Two-phase friction model and nucleate boiling model are the most influential models. The modified results are in better agreement with experimental data. The modifications are effective on prediction of threshold power and period. Abstract: As one of the most common types of oscillatory instability, density wave oscillation (DWO) is of practical importance in the general analysis of two-phase boiling systems. Due to the confinement effect of narrow space, the two-phase phenomena in thin rectangular channels are probably different from those in conventional pipes with large diameters. As the widely accepted system analysis code, RELAP5 has the capability to simulate the DWO phenomena in narrow channels, but the accuracy is limited under the condition of different pressures, flow rates, or subcoolings. In this paper, modifications for related constitutive models are proposed based on single-variable sensitivity analysis. Two-phase friction model and nucleate boiling heat transfer model are selected as the most effective models to improve the predictability of RELAP5. The influence of both models mentioned above is analyzed and it reveals that they are related to two-phase pressure drop and the vapor void fraction distribution along the channel, respectively. The modified calculation is assessed against the experimental results under the wide range of pressure 1–10 MPa, and these comparisons are performed in the map of stability boundary with two dimensionless parameters, phase change number Npch and subcooling number Nsub. Then the heating power and the period of flow oscillation under threshold condition are assessed in detail. The results show that a better agreement can be obtained between the modified calculations and the measured experimental data.