针对裂变产额和半衰期的燃耗计算灵敏度和不确定度分析方法

基于广义微扰理论推导了裂变产额和半衰期的燃耗灵敏度系数理论模型，该模型考虑了原子核密度和中子通量的相互影响，并开发了燃耗计算中有效增殖因数和原子核密度等响应参数对核数据的灵敏度和不确定度分析程序。基于评价核数据中裂变产物独立产额的标准差数据，产生了针对压缩燃耗数据库的裂变产额协方差矩阵，以提高不确定度的计算精度。基于ENDF/B-Ⅶ.1数据库量化了UAM基准题TMI-1栅元无限增殖因数及重要裂变产物和重核的原子核密度由裂变产额和半衰期引入的不确定度。数值结果表明，对于栅元无限增殖因数，裂变产额和半衰期引入的不确定度很小；对于部分裂变产物的原子核密度，裂变产额和半衰期会引入较大的不确定度。     

中子扩散三维非均匀变分节块法及平源加速方法研究

作为三维全堆芯非均匀输运计算的基础,基于扩散近似提出了三维非均匀变分节块法,能够采用节块内部的等参有限元和节块轴向表面的分片常量精细描述燃料棒的非均匀几何结构,消除均匀化过程。同时,为了降低计算代价,相应地提出了平源加速方法(FS)。数值结果表明,该方法具有可靠的精度,且FS能够在不影响精度的前提下有效地降低计算内存和计算时间。     

基于抽样方法的特征值不确定度分析

核数据是反应堆物理计算的基础数据,研究其不确定度对反应堆物理计算引入的不确定度,对提高反应堆的安全性和经济性具有重要意义。本文基于抽样理论研究了反应堆物理计算不确定度分析的方法,研发了不确定度分析程序UNICORN。基于ENDF/B-Ⅶ.1评价数据库,使用NJOY程序开发了多群协方差数据库。采用UNICORN程序和多群协方差数据库对三哩岛燃料棒和基准题RB31的k∞进行了不确定度分析,得到核数据库中各分反应道截面的不确定度对k∞造成的不确定度。结果表明:238 U(n,γ)截面对三哩岛燃料棒k∞造成的不确定度最大,相对不确定度达0.4%左右;协方差数据库的不同来源会对不确定度分析结果造成一定影响。     

基于经典微扰理论的特征值灵敏度和不确定度分析

核数据不确定度作为组件/栅元计算不确定度的重要来源,备受重视和研究。本文采用经典微扰理论,推导输运计算中keff对于核数据的灵敏度系数和不确定度的计算方法。基于ENDF/B-Ⅶ.1制作多群协方差数据库,并根据所采用的组件输运求解程序的截面模型对分反应道协方差矩阵进行归并。开发灵敏度和不确定度分析程序COLEUS,对传统压水堆燃料栅元进行计算分析。数值结果表明,栅元计算的keff对235 U每次裂变中子产额的扰动最为敏感,238 U俘获截面对keff不确定度的贡献最大。目前的核数据的不确定度会给keff带来0.4%~0.5%的不确定度。     

基于AutoCAD二次开发实现中子输运方程特征线法求解

在先进反应堆的组件设计计算中,特征线方法(MOC)是沿生成的特征线求解中子输运方程,理论上不受几何形状的限制,但需对组件进行几何描述和射线追踪等预处理,现有的MOC程序在几何预处理上实际还存在很多限制.为彻底消除特征线方法在几何方面的限制,借助AutoCAD二次开发功能来实现MOC方法的几何预处理.在此基础上开发了MOC程序AutoMOC,对各种问题的计算表明,程序不仅在几何处理上具有很高的灵活性,同时,其计算结果与MCNP等现有程序计算结果符合良好.     

氮化硅-铼嬗变模拟靶材的制造和性能研究

采用无压烧结的方法(在1550℃保温3小时)添加5wt.％MgO,5wt.％Y2Oa和不同含量的Re制造出了Si3N4-Re嬗变模拟靶材.研究了Re含量的不同对嬗变靶材致密度和抗弯强度的影响,并通过物相分析,显微组织和断口形貌观察对相关机理进行了分析.     

子通道程序SUBSC的开发及稳态验证

采用面向对象模块化编程技术开发了面向大规模热工水力计算的自主化子通道程序SUBSC,利用SUBSC和COBRA程序分别计算了典型压水堆1/4组件,结果表明,两者计算结果吻合很好。为进一步验证SUBSC程序,计算了PSBT稳态5×5棒束基准题,结果表明,在各种工况下SUBSC程序计算得到的通道平均含汽率与实验测量值吻合很好,最大相对偏差仅为0.7%,证明了程序具有较高的计算精度。为提高SUBSC程序的计算效率,引入不完全LU分解预处理的再启动GMRES算法求解质量守恒方程,对多组件的计算结果表明,SUBSC程序具备大规模热工水力计算能力。     

CASMO-4E/SIMULATE-3程序系统计算BEAVRS基准题

BEAVRS基准题是麻省理工学院计算反应堆物理小组2013年公布的压水堆三维全堆芯高保真计算基准题。本文使用传统两步法程序系统CASMO-4E/SIMULATE-3对其进行建模与跟踪计算。结果表明:在热态零功率(HZP)工况下,径向探测器反应率最大相对误差为-18.8%,与蒙特卡罗程序MC21程序的-16.1%相当;在燃耗深度为38.7 MW·d/tHM、堆芯功率为20.3%FP时,径向探测器反应率最大相对误差为-12.9%;在堆芯第1循环内,堆芯各燃耗点处临界硼浓度与测量值误差在40ppm以内,满足工程精度要求。     

基于预估修正改进准静态方法的中子时空动力学研究

研究了一种基于预估修正思想的改进准静态方法,用于求解中子时空动力学方程。通过隐式差分法预估下一时刻中子注量率以产生形状函数,避免了传统的改进准静态方法中形状函数的迭代计算,并依靠幅函数修正中子注量率,得到最终中子注量率的分布。计算表明,在较大时间步下,该方法计算速度较传统的改进准静态方法有明显的提高,且能保持较高的计算精度。     

强各向异性散射中子输运方程加速方法

采用修正的PN方程加速中子输运方程的SN数值求解过程.并将该加速方法应用到国际通用的ANISN程序中,对其进行了数值检验.数值计算结果表明,本文加速方法可以使迭代次数大量减少,特别是对于强各向异性散射问题具有明显的加速效果.