核反应堆物理分析第5课.pptVIP

5.3.2 二维扩散方程的数值解法 计算中若不考虑自低能群的向上散射，方程(5-52)可写成 可以看出，当按次序g=1,2,3· · ·,G次序对以上方程求解时，方程右边源项便已知。略去上下标以上方程可以写为： 以二维问题为例讨论以上方程的数值解： 对所取计算平面使用以下 直线分成矩形网格: 共有(N+1)(M+1)个节点， ?xi=xi+1-xi， ?yi=yi+1-yi 将方程(5-64)在(i, j)节点附近 积分，方程第一项为: 二维网格 （i，j）网格示意图 上式导数中应用下面近似公式 并考虑到D有可能不连续得 同理可以得到第二、三、四项积分等于 综合以上各项，方程(5-64)在(i, j)节点附近差分方程为: 对于边界上的点则和前面一样，需要根据边界条件来确定。 对于四周外推边界上，所有网点有： 方程（5-70）便是我们所需要的差分方程组，对于每次源 迭代，方程(5-70)的系数及源项已知。通常采用迭代方法 对该方程组进行求解。这种在源迭代的过程中，对方程 (5-70)求解的过程，为内迭代或通量密度迭代。 求解顺序通常为： 自第一行j＝０开始逐行向上至j=M； 在每一行，从i=0开始，自左向右到i=N。 迭代的方法可以根据具体情况选取。每种迭代方法都有局 限性，有些方法收敛快，而有些方法收敛慢。 从前面的讨论可知，多群扩散方程的数值求解步骤可以 总结为两大步： 源迭代或外迭代 通过源迭代法将多群联立方程组的求解变换为解G个 单群扩散方程，每个单群扩散方程应用差分法通过内迭代方法求解。 中子通量密度迭代或内迭代 由于总群数为G，因而每一次外迭代中都需要对G个方程 作内迭代运算。没有向上散射时，求解的次序可以从高能群 顺序向下运算 临界计算所涉及的问题： 给定堆芯几何形状及材料成分等，求芯部keff 临界搜索 （1）给定芯部成分条件下要求反应堆的芯部尺寸 ---临界尺寸搜索 （2）给定芯部尺寸条件下要求芯部的燃料成分、 控制棒或吸收体的数量和位置--- 临界成分搜索 第 5 章 分 群 扩 散 理 论 单群理论认为堆芯内所有的中子能量相等，而实际上堆芯 内的中子能量分布在很宽的范围。材料的各种截面又和中子 能量密切相关，实际应用中单群理论是不适用的，必须考虑 能量相关的中子慢化与扩散模型。 尽管单群理论给出的结果不够准确，但单群理论简单 明了，在一些情况下能给出解析解，有利于初学者掌握和 理解分群扩散理论概念和方法，并且其解带有普遍意义。 因此, 研究单群理论是非常有意义的. 本章主要内容: 建立多群中子扩散方程、介绍群常数的 计算、以双群模型求解求解反应堆临界问题、介绍多群扩散 方程的数值解法。 5.1 与能量相关的中子扩散方程和分群扩散理论 5.1.1 与能量相关的中子扩散方程 泄漏率： J(r,E) 表示r处能量为E的中子净流密度， 由斐可定律（3-21）， 则单位体积内泄漏率为： 损失率：中子损失率包括吸收损失和散射引起的损失。 产生率： 根据相空间体积元内中子数的变化率等于中子产生率与消失 率之差，可以建立与能量相关的中子扩散方程 在稳态和无源情况下，方程变为： 上式只是对于临界系统才成立。 任意系统稳态中子扩散方程 有效增殖因子keff也称为方程的特征值。 5.1.2 分群扩散理论及多群中子扩散方程 中子能量区域按能量大小 划分为G个能区，最高能量 记为E0 ，最低能量记为EG ， 每一个能量区间称为一个 能群。 能群编号g=1, …，G随中子能量降低而增加。 这样： 能量连续变化的扩散方程变成G个与能量无关的扩散方程。 中子能群的划分 在每一个能群区间对稳态中子扩散方程进行积分，消去方 程中能量，可得G个不含能量变量E的扩散方程，其中 第g群扩散方程为 为了简化以上方程以方便求解，引入一些主要物理量： g群中子通量密度 ?g(r)就是空间r处g群内各种能量中子的总中子通量密度。 g群总截面 g群扩散系数 g群中子产生截面 g群中子裂变谱 群转移截面 ?g’ →g 能量分成G群后，中子扩散方程中的散射源可以写为： 因而， 便表示每秒单位体积内gˊ群中子经受散射碰撞 后，能量落到g群内的中子数。 利用以上定义的群常数，对于g群中子的扩散方程为： 这就是多