基于指数变换的瞬态扩散方程求解方法研究
Study on the Transient Neutron Diffusion Equation Solving Method Based on Exponential Transform

作者: 庞博 , 赵强 , 郝琛 , 宇炎 :哈尔滨工程大学核安全与仿真技术国防重点学科实验室,黑龙江 哈尔滨;

关键词: 中子时空动力学指数变换时间积分法共轭平方法Neuron Kinetics Model Exponential Transform Time Integration Method Conjugate Gradient Squared Method

摘要: 本文重点研究了基于指数变换的瞬态扩散方程求解方法,时间上采用指数变换减弱刚性问题的影响,并采用时间积分法提高计算精度。几何上采取细网有限差分以提高计算精度,并且采用平方共轭梯度法进行最后的方程组求解。本文基于上述方法开发了三维多群中子时空动力学计算程序。

Abstract: Transient diffusion equation solution method based on exponential transform is studied mainly. Time integral aspect exponential transform weakens the influence of the stiff problem, and the time integration method is used to improve calculation accuracy. Fine mesh finite difference method is taken to improve the geometry calculation accuracy, and conjugate gradient squared method is used for solving final equations. Based on the above methods, multi-group neutron kinetics model calculation program is developed.

文章引用: 庞博 , 赵强 , 郝琛 , 宇炎 (2016) 基于指数变换的瞬态扩散方程求解方法研究。 核科学与技术, 4, 41-49. doi: 10.12677/NST.2016.42006

参考文献

[1] 胡大璞. 求解中子时空动力学方程的动态反应性法[J]. 核动力工程, 1994, 15(1): 45-53.

[2] 焦惠先, 李惠云. 小型压水堆弹棒事故反应性分析[J]. 核科学与工程, 1987, 7(3-4): 219.

[3] 赵文博. 瞬态节块格林函数方法及其与热工–水力耦合研究[D]: [博士学位论文]. 北京: 清华大学, 2012.

[4] 黄祖恰. 核反应堆动力学基础[M]. 北京: 北京大学出版社, 2007: 420-421.

[5] Reed, W.H. and Hansen, K.F. (1970) Alternating Direction Methods for the Reactor Kinetics Equations. Nuclear Science and Engineering, 41, 431-442.

[6] 谢仲生. 核反应堆物理数值计算[M]. 北京: 原子能出版社, 1997: 160.

[7] Zimin, V.G. and Ninokata, H. (1998) Nodal Neutron Kinetics Model Based on Nonlinear Iteration Procedure for LWR Analysis. Annals of Nuclear Energy, 25, 507-528.
http://dx.doi.org/10.1016/S0306-4549(97)00078-9

分享
Top