反应堆物理计算中的边界流响应矩阵方法研究

发布时间：2020-07-08 07:17

【摘要】：在反应堆物理计算中，传统的等效均匀化计算一般采用体积-通量权重法，在反射边界条件下，该方法可以确保等效前后各群反应率及平均通量守恒。但是，在边界条件产生明显变化时，该方法会有较大的计算误差。基于响应矩阵方法的基本思路，本文提出了一种新的栅元等效均匀化计算方法，通过利用栅元边界入射流与出射流之间的响应关系进行等效。与在传统方法中保证等效前后平均通量守恒不同的是，新的等效均匀化方法的目标是保证入射流与出射流以及各种不同反应率之间的响应关系守恒。与体积-通量权重法相比，新的等效均匀化方法几乎不增加计算量，而在带反射层的平板问题和二维C5G7MOX基准例题上分别进行的测试计算结果表明，在特征值和栅元功率分布的计算精度方面，新的栅元等效均匀化方法具有明显的改善作用。 基于响应矩阵方法的原理，针对反应堆堆芯中存在大量重复结构这一特点，本文进一步研究了基于边界流响应矩阵的组件等效和全堆迭代计算方法。其基本思想是对组件进行等效，但并不是等效为一个均匀材料，而是利用组件表面出射流与入射流之间的响应关系进行等效。在特征值问题计算过程中，本文采用将组件各群裂变截面都除以堆芯keff的方法，使组件边界出射流与入射流的耦合关系为线性关系，从而将特征值问题转化为固定源问题。这样只需在全堆迭代前计算得到各组件的响应矩阵，再通过各组件的位置关系、利用边界流的传递与响应在全堆芯范围内进行迭代，即可求得堆芯keff以及各个组件的边界流，进而得到堆芯各处的通量分布。本文基于上述方法编写了计算程序，并在二维C5G7MOX基准例题等计算模型上进行了验证计算，取得了较高的计算效率和精度。 与其他输运计算方法相比，响应矩阵方法在计算速度方面可以提高至少两个量级。本文进一步研究了在不同网格数的情况下响应矩阵方法的计算时间，验证了该方法在计算速度上的优势。而组件边界流展开的近似处理使得该方法在栅元功率分布计算精度上有待提高，为此本文提出了形状因子修正方法以及特征响应矩阵方法。形状因子修正方法可以部分改善栅元功率分布精度，而在反射边界条件下计算响应矩阵的特征响应矩阵方法则具有进一步的研究潜力。
【学位授予单位】：清华大学
【学位级别】：硕士
【学位授予年份】：2014
【分类号】：TL329
【图文】：

图 1.1 传统的反应堆堆芯物理计算流程示意图元均匀化过程中，首先对组件内各种类型的栅元进行等效均化时，根据适当的几何模型，对能群采用多群近似（如 69 的输运计算方法，求出空间通量分布，并群的同时进行空

栅元几何等效示意图

二维C5G7-MOX基准例题1/4堆芯布置图

【参考文献】

相关期刊论文 前5条

1 成红武，谢仲生，尹邦华，章宗耀;组件均匀化常数计算[J];核动力工程;1995年03期

2 李文淮;刘志宏;;组件均匀化中的边界条件误差修正方法研究[J];核科学与工程;2011年03期

3 汤春桃;;线性源近似的中子输运方程特征线解法[J];原子能科学技术;2011年12期

4 赵晶;彭良辉;刘召远;刘志宏;张剑;;SN2D程序在C5G7基准例题上的验证与分析[J];原子能科学技术;2013年S1期

5 刘志宏;刘召远;彭良辉;赵晶;张剑;;采用边界流响应矩阵的组件等效计算方法研究[J];原子能科学技术;2013年S1期

本文编号：2746258