一种稳定性增强及高精度数值方法在RELAP5中的实现与评价

李伟; 苏光辉; 秋穗正; 田文喜; 陈平; 李垣明

doi:10.11889/j.0253-3219.2016.hjs.39.110601

您当前的位置：

首页 >

文章列表页 >

一种稳定性增强及高精度数值方法在RELAP5中的实现与评价

核能科学与工程 | 更新时间：2021-05-24

- 一种稳定性增强及高精度数值方法在RELAP5中的实现与评价
- Implementation and assessment of a stability-enhancing and high-resolution numerical scheme in RELAP5
- 核技术 2016年39卷第11期文章编号：110601
- 作者机构：
  
  1.中国核动力研究设计院核反应堆系统设计技术重点实验室成都 610213
  2.西安交通大学核科学与技术学院西安 710049
- 作者简介：
  
  [ "李伟,男,1986年出生,2015年于西安交通大学核科学与技术专业获博士学位,主要从事核燃料设计、反应堆热工水力及安全分析", "LI Wei, male, born in 1986, graduated from Xi’an Jiaotong University with a doctoral degree in 2015, major in nuclear science and technology, mainly focusing on nuclear fuel research and design, reactor thermal-hydraulic and safety analysis" ]
- 基金信息：
- DOI：10.11889/j.0253-3219.2016.hjs.39.110601
  中图分类号： TL33
- 收稿日期：2016-05-30，
  
  修回日期：2016-09-01，
  
  纸质出版日期：2016-11
- 稿件说明：
移动端阅览
李伟, 苏光辉, 秋穗正, 等. 一种稳定性增强及高精度数值方法在RELAP5中的实现与评价[J]. 核技术, 2016,39(11):110601

Wei LI, Guanghui SU, Suizheng QIU, et al. Implementation and assessment of a stability-enhancing and high-resolution numerical scheme in RELAP5[J]. Nuclear techniques, 2016, 39(11): 110601
李伟, 苏光辉, 秋穗正, 等. 一种稳定性增强及高精度数值方法在RELAP5中的实现与评价[J]. 核技术, 2016,39(11):110601 DOI： 10.11889/j.0253-3219.2016.hjs.39.110601.

Wei LI, Guanghui SU, Suizheng QIU, et al. Implementation and assessment of a stability-enhancing and high-resolution numerical scheme in RELAP5[J]. Nuclear techniques, 2016, 39(11): 110601 DOI： 10.11889/j.0253-3219.2016.hjs.39.110601.

摘要

在计算稳定性的改进方面，修正RELAP5程序中的虚拟质量力（Virtual mass force）形式，同时添加了新的界面压力项（Interface pressure）；在计算精度的改进方面，采用具有总变差减小（Total variation diminish，TVD）特点的高精度通量限制器（Flux limiter）方法取代RELAP5程序原本的一阶迎风方法来离散质量和能量守恒方程中的对流项。在模拟水平管道内空泡份额微扰随时间发展的数值实验中，相比改进前的RELAP5，改进后的RELAP5计算得到的微扰幅度并未增长；在模拟液相沉降的数值实验中，改进前的RELAP5程序计算得到了不真实的空泡份额分布，而改进后的RELAP5在不同的网格数量下能够得到收敛的稳定解。对Ransom水龙头数值实验和Marviken CFT 15大破口喷放实验的计算表明，改进后的具有二阶TVD格式的RELAP5程序能够得到更接近实验数据的计算结果。

Abstract

Background

RELAP5 has been widely used in the field of nuclear thermal-hydraulics and safety analysis. However

improvement issues in terms of numerical stability and accuracy exist due to its basic single-pressure two-fluid model.

Purpose

This study aims to enhance the numerical stability and numerical accuracy of RELAP5 by modifying the mathematical structure of RELAP5 basic model.

Methods

Firstly

the default mathematical expression of the virtual mass force was modified in the RELAP5 code. Secondly

an interface pressure term was incorporated in the code

and a total variation diminish (TVD)-type flux limiter method was implemented in the RELAP5 code in place of the default 1st-order upwind scheme for the convective terms in mass and energy conversation equations.

Results

In a numerical test where the evolvements of void fraction perturbation were concerned

the modified RELAP5 predicted no growth of the perturbation amplitude in contrast to a rapidly developing divergence for the default RELAP5. In addition

for the phase segmentation problem

convergence was accomplished on finer mesh by the modified RELAP5

while non-physical distribution of the void fraction was rendered by the default RELAP5 especially for finer mesh. When applied to the Ransom water faucet numerical test and Marviken CFT 15 experiments

the improved RELAP5 with 2nd-order of accuracy was more reliable in that the predictions were closer to the experimental results.

Conclusion

Compared with the original RELAP5

the modified or improved RELAP5 showed evident numerical stability and accuracy during the assessment in this work.

关键词

Keywords

references

L Ge , H T Wang , G L Zhang , . Thermal-hydraulic design and transient analysis of passive cooling system for CPR1000 spent fuel storage pool . Nuclear Science and Techniques , 2016 . 27 ( 1 ): 8 DOI: 10.1007/s41365-016-0017-6 http://doi.org/10.1007/s41365-016-0017-6 .

R W Lyczkowski , D Gidaspow , C W Solbrig , . Characteristics and stability analyses of transient one-dimensional two-phase flow equations and their finite difference approximations . Nuclear Science and Engineering , 1978 . 66 ( 3 ): 378 - 396 . DOI: 10.13182/NSE78-4 http://doi.org/10.13182/NSE78-4 .

P D Lax . On the nation of hyperbolicity . Communications on Pure and Applied Mathematics , 1980 . 33 ( 3 ): 395 - 397 . DOI: 10.1002/cpa.3160330309 http://doi.org/10.1002/cpa.3160330309 .

D Wang , J H Mahaffy , J Staudenmeier , . Implementation and assessment of high-resolution numerical methods in TRACE . Nuclear Engineering and Design , 2013 . 263 ( 0 ): 327 - 341 . DOI: 10.1016/j.nucengdes.2013.05.015 http://doi.org/10.1016/j.nucengdes.2013.05.015 .

张小英 , 丁斐 , 陈佳跃 . 反应堆一维两流体模型二阶精度数值解法研究 . 核动力工程 , 2013 . 34 ( 4 ): 27 - 32 .

Xiaoying ZHANG , Fei DING , Jiayue CHEN . Study on numerical solution of the two dimensional two fluid model of the reactor . Nuclear Power Engineering , 2013 . 34 ( 4 ): 27 - 32 . DOI: 10.3969/j.issn.0258-0926.2013.04.007DOI:10.3969/j.issn.0258-0926.2013.04.007 http://doi.org/10.3969/j.issn.0258-0926.2013.04.007DOI:10.3969/j.issn.0258-0926.2013.04.007 .

D Drew , L Cheng , R Lahey . The analysis of virtual mass effects in two-phase flow . International Journal of Multiphase Flow , 1979 . 5 ( 4 ): 233 - 242 . DOI: 10.1016/0301-9322(79)90023-5 http://doi.org/10.1016/0301-9322(79)90023-5 .

I Tiselj , S Petelin . Modelling of two-phase flow with second-order accurate scheme . Journal of Computational Physics , 1997 . 136 ( 2 ): 503 - 521 . DOI: 10.1006/jcph.1997.5778 http://doi.org/10.1006/jcph.1997.5778 .

L Hogon , U Lee , K Kim , . Application of hyperbolic two-fluids equations to reactor safety code . Nuclear Engineering and Technology , 2003 . 35 ( 1 ): 45 - 54 . http://cn.bing.com/academic/profile?id=1577674790&encoded=0&v=paper_preview&mkt=zh-cn .

RELAP5 Code Development Team. RELAP5/MOD3 code manual Vol I:code structure, system models, and solution methods[R]. Idaho Falls:Idaho National Laboratory, 2001

H K Versteeg , W Malalasekera . An introduction to computational fluid dynamics:the finite volume method , : Essex, UK Pearson Education , 2007 .

N P Waterson , H Deconinck . Design principles for bounded higher-order convection schemes-a unified approach . Journal of Computational Physics , 2007 . 224 ( 1 ): 182 - 207 . DOI: 10.1016/j.jcp.2007.01.021 http://doi.org/10.1016/j.jcp.2007.01.021 .

Ransom V. Numerical benchmark tests[A]. Hewitt G, Delhay J, Zuber N. Multiphase science and technology[M]. Washington DC:Hemisphere Publishing Corporation, 1987

Kim K, Kim H J. Assessment of RELAP5/MOD 2 critical flow model using Marviken test data 15 and 24[R]. Washington, DC:Office of Nuclear Regulatory Research, U.S. Nuclear Regulatory Commission, 1992

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

二次开发的RELAP5在倾斜条件下对PRS运行分析的适用性研究

小型堆二次侧非能动余热排出系统特性计算分析

基于EPICS架构的熔盐堆实时建模与仿真平台开发及验证

基于SiPM的塑料闪烁光纤剂量计的研制

用于产氚率测量的碳酸锂探测片液闪样品制备方法研究