第5章分群理论2014.ppt

* 采用因式分解法得 同理，若消去?2,c,可得 (5-26) (5-27) 式中 (5-28) (5-29) 对于临界反应堆，k’?大于1，?2、?2为正的实数。由(5-26)及(5-27)看出， ?1,c及?2,c均满足波动方程 其中B2等于?2、-?2。用-B2?代替?2?，则连立求解(5-18)、(5-19)可得双群理论的增殖因子 (5-30) (5-31) * 如果假设快群不发生裂变，则?f1,c=0，则上式可简化为 其中 (5-32) (5-33) 方程(5-26)及(5-27)的解是下列两个波动方程组解的线性组合 ?1,c, ?2,c的一般解可写为 (5-34) (5-35) (5-36) (5-37) * 可以证明：4个待定常数A,C, A’, C’中只有两个是独立的。 都是方程的允许解，因而 ，根据方程(5-19)可得 令 s1=A’/A，由上式可得 (5-38) 同理可得 s1、s2叫做耦合系数，其值由系统材料属性决定，因此：待定常数A,C, A’, C’中只有两个是独立的。 (5-39) * 芯部中子通量密度的普遍解为 (5-40) (5-41) 例如，对于侧面带有反射层的圆柱形反应堆，中子通量密度为 其中 * 2. 反射层方程的解 (5-22) (5-23) 方程(5-22)为齐次方程，解为 (5-42) 方程(5-23)为非齐次方程，解可写为 (5-43) 方程(5-43)代入方程(5-23)得 整理得 对于侧面带反射层的圆柱形反应堆，反射层中子通量密度为 * 5.2.3 双群临界方程及中子通量密度分布 双群临界方程 双群临界方程的四个边界条件 方括号表示函数取芯部与反射层交界面上的值。 根据方程(5-40),(5-41,(5-42),(5-43),根据边界条件得 其中 (5-45) (5-46) (5-47) (5-48) * 为使方程组有非零解，根据克莱姆法则，系数判别式须为零，可得双群临界方程 (5-51) (4-46) (4-47) (4-62) 中子通量密度分布 * * 5.3 多群扩散方程的数值解法 解析方法仅适用于简化模型的求解（两群、两区、一侧有反射层） 数值方法是实际工程中唯一有效的求解扩散方程的方法 多群扩散方程是一个齐次特征值问题，不能直接离散求解，通常采用“外-内迭代方法”求解。 外迭代（源迭代）：通过迭代求特征值的过程 内迭代：对外迭代过程中出现的扩散方程进行数值求解的过程 * 5.3.1 源迭代法 无外源情况下，多群扩散方程可写为 (5-52) (5-53) 源迭代法过程： 假定一个初始裂变源分布及初始特征值，即 初始裂变源分布及初始特征值代入(5-52)，从高能到低能，逐群求解中子扩散方程，得中子通量密度分布?g(1)。 将?g(1)代入方程(5-53)，得第二代裂变中子源分布Q(1)(r)。 根据Q(1)(r)求得第二代特征值keff(1)，得第二代迭代源项。 将第二代源项代入(5-52)求解，得下一代中子通量密度分布。 … 直到收敛。 * 对第n次迭代 (5-54) (5-55) 根据keff的物理意义得 (5-56) 已经证明 (5-57) 收敛准则： 特征值收敛准则 中子通量密度分布收敛准则 (5-58) (5-59) * 内迭代 前面我们已经说过，解多群扩散方程是用所谓源迭代法来做的，即预先设定一个源项和一个Keff。在计算每一群通量分布时，源项是已知的。计算每一群通量分布时是用的也是迭代法，这个迭代过程称为内迭代。 源迭代也称为为外迭代 * 5.3.2 二维扩散方程的数值解法 计算中不考虑自低能群的向上散射(up-scattering)，方程(5-52)可写成 (5-60) (5-61) (5-62) 略去上下标 (5-63) 以二维（x,y）问题为例，讨论其解法。 共有(N+1)(M+1)个节点，?xi=xi+1-xi， ?yi=yi+1-yi 对所取计算平面使用以下直线族 (5-64) * 二维网格 (i, j)网格示意图 对方程(5-64)在(i, j)节点附近积分，方程第一项为 (5-65) 应用以下近似公式 得 (5-66) * 方程第二、三、四项分别为 (5-67) (5-68) (5-69) 综合(5-66)、(5-67)、(5-68)和(5-69)得(i, j)节点差分方程 (5-70) * 边界条件：外推边界处 对于每次源迭代，方程(5-70)的系数及源项已知。通常采用迭代方法对该方程组进行求解。这种在源迭代的过程中，对方程(5-70)求解的过程，为内迭代或通量密度迭代。 求解顺序一般为： 自第一行j＝０开始逐行向上至j=M； 在每一行，从i=0开始，自左向