强激光与粒子束, 2017, 29 (3): 036024, 网络出版: 2017-03-29  

JFNK在高温堆扩散计算中的应用

Application of Jacobian-free Newton-Krylov method for high temperature reactor neutron diffusion equation calculation
作者单位
清华大学 核能与新能源技术研究院, 先进核能技术协同创新中心, 先进反应堆工程与安全教育部重点实验室, 北京 100084
摘要
研究了JFNK框架下高温堆中子扩散问题的求解方法。研究结果表明,JFNK方法在求解与源迭代相同形式中子扩散方程时,相对残差下降趋势为逐渐加快并趋于稳定,有利于更高求解精度的实现。使用通量归一化附加方程可以获得更好的JFNK非线性迭代特性,但在算例中其部分牛顿修正方程求解时间偏多,总计算时间高于显式有效增殖系数附加方程法,需要研究更高效的JFNK预处理方法对线性求解环节进行改善。
Abstract
This paper studies the application of solving high temperature reactor (HTR) neutron diffusion equation with Jacobian-free Newton-Krylov (JFNK) method. Results show that when solving neutron diffusion equation, the relative residual norm of JFNK method decreases slowly at the beginning. Then the rate of convergence become faster and finally reaches a relatively stable value. This feature is conducive to a high-accuracy solution. In the test of two kinds of additional equations, the neutron diffusion equation with flux normalization condition has a better nonlinear convergence behavior. However, due to the longer computational time in solving linear equations, its total computational time is more than the one with k expression. More efficient preconditioning methods should be studied to improve linear equations.
参考文献

[1] Zhang Zuoyi, Wu Zongxin, Wang Dazhong, et al. Current status and technical description of Chinese 2× 250 MWth HTR-PM demonstration plant[J]. Nuclear Engineering and Design, 2009, 239(7): 1212-1219.

[2] Zhang Zuoyi, Wu Zongxin, Sun Yuliang, et al. Design of Chinese modular high-temperature gas-cooled reactor HTR-PM[C]//Proceedings of the 2nd International Topical Meeting on High Temperature Reactor Technology. 2004.

[3] 吴宏春,曹良志,孙幸光,等.中子输运方程数值计算的加速方法研究[C]//第十二届反应堆数值计算与粒子输运学术会议论文集. 2008.(Wu Hongchun, Cao Liangzhi, Sun Xingguang, et al. Research of acceleration methods in neutron transport equation numerical calculation. CORPHY-12 Proceedings. 2008)

[4] Gill D F, Azmy Y. Newton’s method for solving k-eigenvalue problems in neutron diffusion theory[J]. Nuclear Science and Engineering, 2011, 167(2): 141-153.

[5] 吴文斌,李庆,孙伟,等.采用JFNK方法求解三维中子扩散方程[J].核动力工程, 2015(s2): 17-19.(Wu Wenbin, Li Qing, Sun Wei, et al. Solving 3-D Neutron diffusion equation with JFNK algorithm. Nuclear Power Engineering, 2015(s2): 17-19)

[6] Gill D F, Azmy Y Y, Warsa J S, et al. Newton’s method for the computation of k-eigenvalues in SN transport applications[J]. Nuclear Science and Engineering, 2011, 168(1): 37-58.

[7] Fichtl E D, Warsa J S, Densmore J D. The Newton-Krylov method applied to negative-flux fixup in SN transport calculations[J]. Nuclear Science and Engineering, 2010, 165(3): 331-341.

[8] Knoll D A, Keyes D E. Jacobian-free Newton-Krylov methods: A survey of approaches and applications[J]. Journal of Computational Physics, 2004, 193(2): 357-397.

[9] Eisenstat S C, Walker H F. Choosing the forcing terms in an inexact Newton method[J]. SIAM Journal on Scientific Computing, 1996, 17(1): 16-32.

[10] Hageman L A, Young D M. Applied iterative methods[M]. North Chelmsford, Massachusetts: Courier Corporation, 2004.

[11] Gill D F. Newton-Krylov methods for the solution of the k-eigenvalue problem in multigroup neutronics calculations[M]. The Pennsylvania State University, 2009.

卢佳楠, 郭炯, 李富. JFNK在高温堆扩散计算中的应用[J]. 强激光与粒子束, 2017, 29(3): 036024. Lu Jia’nan, Guo Jiong, Li Fu. Application of Jacobian-free Newton-Krylov method for high temperature reactor neutron diffusion equation calculation[J]. High Power Laser and Particle Beams, 2017, 29(3): 036024.

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!