改进准静态本征函数法求解反应堆时空中子动力学方程

一第50卷第7期

原子能科学技术Vol.50,No.7一2016年7月Atomic Ener gy Science and Technolo gy Jul.2016

改进准静态本征函数法求解反应堆

时空中子动力学方程

李明芮1,黎浩峰1,2,陈文振1,邢一晋1,

3(1.海军工程大学核能科学与工程系,湖北武汉一430033;2.海军司令部核安全部,北京一100841;3.海军装备部驻重庆地区军事代表局,重庆一400000)摘要:本文提出一种改进准静态本征函数法并用于求解时空中子动力学方程三在改进准静态近似条件下,利用因子分解法将中子通量密度分解为幅函数和形状函数,进而采用本征函数法求解形状函数轴向和径向的近似解析解三在求解幅函数过程中,将幅函数化简为点堆模型下的形式进行求解,最终导出反应堆三维时空中子动力学方程的解析解三与其他解析法和数值解相比,改进准静态本征函数法的适用范围更广,计算速度更快三

关键词:改进准静态本征函数法;时空中子动力学;因子分解法

中图分类号:TL326一一一文献标志码:A一一一文章编号:1000-6931(2016)07-1245-06

收稿日期:2015-07-29;修回日期:2015-12-21基金项目:国家自然科学基金资助项目(11301540)作者简介:李明芮(1989 ),男,山东枣庄人,博士研究生,从事核反应堆安全分析研究doi :10.7538/y zk.2016.50.07.1245Im p roved Q uasi -static Ei g enfunction Method for Solvin g S p atial -time Neutron Kinetics E q uation

LI Min g -rui 1,LI Hao -fen g 1,2,CHEN Wen -zhen 1,XING Jin 1,3

(1.De p artment o f Nuclear Ener g y Science and En g ineerin g ,Naval Universit y o f En g ineerin g ,Wuhan 430033,China ;

2.Nuclear Sa f et y De p artment o f Naval Command ,Bei j in g 100841,China ;

3.Naval E q ui p ment Bureau o f Chon gq in g ,Chon gq in g 400000,China )Abstract :一An im p roved q uasi -static method was p ro p osed to solve the s p ace -time neu -tron kinetics e q uation.Under the im p roved q uasi -static a pp roximation condition ,the neutron flux densit y was decom p osed into the ran g e -function and sha p e -function b y the factorization method.Then the a pp roximate anal y tical solutions in the axial and radial directions were obtained usin g the ei g enfunction method to solve the sha p e -function.In the solvin g p rocess ,the ran g e -function was sim p lified as the form of p oint reactor model ,and the anal y tical solution of three -dimensional s p ace -time neutron kinetic e q ua -tion for reactor was https://www.sodocs.net/doc/9917221463.html, p ared with other anal y tical method and numerical solution ,the im p roved q uasi -static ei g enfunction method has wider a pp lication sco p e and faster com p utation s p eed.