核反应堆物理课设解决方案.docVIP

课程设计(综合实验)报告 名 称： 核反应堆物理分析 题 目：利用双群理论求解堆芯参数 院 系： 1111111111111 班 级： 111111111111111 学 号： 111111111111 学生姓名： 11111111111 指导教师： 111111111 设计周数： 成 绩： 课程设计 1.0 0.0 1.0 0.0 0.0075 0.1651 0.0 0.0 0.0121 0.121 0.0004 0.020 0.0241 - 0.0493 - 1.26 0.354 1.130 0.166 利用数值求解下面问题，计算其快群和热群中子通量密度空间分布及有效增殖因子。 提示： 采用内外迭代方法求解； 给出源迭代过程； 给出内迭代过程并给出堆芯双群方程、反射层双群方程以及边界条件，详见课本第五章； 建立数学模型，由于问题具有对称性，因此，只计算1/4堆芯即可。 可用赛德尔迭代法求解；（详见课本第五章）； 给出计算流程 编写程序，给出结果，并用OREGIN画图把结果表现出来； 结果的分析讨论； 写出报告。 2.1差分方程的建立 为简化计算，进行以下假设： （1）二维XY几何模型是一个整体，即堆芯和反射层交界处的中子通量密度和中子流连续； （2）认为该问题中反射层厚度已包含外推边界； （3）认为该问题中堆芯和反射层对称分布。 在带反射层的多区双群扩散方程的数值解法中，通常是通过外-内迭代过程来求解。外迭代是一个通过迭代求特征值的迭代过程，所以又成为源迭代。内迭代则是对源迭代中出现的扩散方程进行具体的数值求解。 计算中若不考虑中子自低能群的向上散射，多群扩散方程可以写为： （1） 略取上下标，该方程可以写成以下标准形式： （2） 本设计讨二维（x，y）情况下的差分方程。对所取平面首先用 x = ，，…，，…， y = ，，…，，…， 直线族把平面分成许多矩形网格（如图一）。交点（，）称为节点或网点，如图一中共有（N+1）(M+1)个节点， 称为网距。 图1 网点的划分 图2 （i，j）节点示意图 按照数值解法的相关要求和定义，讨论如图二所示的（i，j）节点，根据二维几何情况下扩散方程导出（i，j）节点的差分方程为 （3） 式中系数为 （4） 对于边界上的点需要根据边界条件来确定，本设计由于反射层已经包括了外推距离，所以外围上的点的中子通量密度均设为0。 上述（3）和（4）式便是求解单群中子扩散的差分方程组，其系数，，…，等均可以事先计算求得，对每次外迭代，其右端源项也是已知的。因此，它是一个含有的代数方程组，对于二维问题，其系数矩阵是一个五队角的矩阵，可以很方便的应用通常的线性方程组的求解方法进行求解。 2.2源迭代过程 根据连续性的假设条件，堆芯的双群方程与放射层得双群方程可以归结为下面两个方程： 快群中子扩散方程： （5） 热群中子扩散方程： （6） 设 （7） 式（5）是一个关于的齐次方程组。式（5）除了中子通量密度外，和也是未知参数，并且这两个参数又与中子通量密度有关。因此我们假定一个初始的裂变源分布，认为在整个堆芯都等于1或某个给定常数，并猜测一个初始的值，同时把作为初始迭代源项代入到方程（5）的右端中。这样，方程（1）的右端便是一个已知项，方程（1）便成为一个普通的非齐次方程组，它可以用通常的数值方法求解。假定由它求出的中子通量密度分布为，把代入方程（6），由方程（6）数值求解求出，将求得的、带入（7）式便可求得第二次迭代裂变中子源分布，并由它可以求得有效增殖因数的新的估计值和第二代迭代源项来，……。依此类推，逐次的迭代下去，对于第n次迭代计算有