1、第三章第三章 中子扩散理论中子扩散理论 中子在介质中的输运过程中中子在介质中的输运过程中的运动状态由位置矢量的运动状态由位置矢量r(x,y,z),能量能量 E, 和运动方向和运动方向表示。表示。通过极角通过极角和方位角和方位角来来表示表示中子角密度函数中子角密度函数n(r,E, )定义：定义：在在r处单位体积内和能量为处单位体积内和能量为E的的单位能量间隔内，运动方向为单位能量间隔内，运动方向为 的单位立体角内的中子数目。的单位立体角内的中子数目。 中子角通量密度中子角通量密度定义为定义为:对中子角密度和中子角通量密度积分便可得到与运动方向无对中子角密度和中子角通量密度积分便可得到与运动方向无
2、关的标量中子密度和标量中子通量密度关的标量中子密度和标量中子通量密度这些量是反应堆物理经常需要计算的量这些量是反应堆物理经常需要计算的量。ddrddrrdSdsinsin222)(),(),(EvErnErdErnErn4),(),(dErEr4),(),(方向方向 的表示的表示要求解反应堆内中子密度和中子通量密度的分布一般要求解反应堆内中子密度和中子通量密度的分布一般采用两种方法：采用两种方法：u确定论方法确定论方法-根据边界条件和初始条件解数学物理方程根据边界条件和初始条件解数学物理方程 得出所求问题的精确解或近似解。得出所求问题的精确解或近似解。 适用于问题的几何结构不太复杂的情况。适用
3、于问题的几何结构不太复杂的情况。u非确定论方法非确定论方法又称为又称为Monte Carlo方法，是基于统计方法，是基于统计 概率理论的方法，适用于问题的几何结构概率理论的方法，适用于问题的几何结构 比较复杂的情况。比较复杂的情况。本章是用确定论方法研究中子的输运过程建立描述中子在本章是用确定论方法研究中子的输运过程建立描述中子在介质输运过程的介质输运过程的中子扩散方程中子扩散方程。中子扩散方程是研究中子。中子扩散方程是研究中子在介质内运动的基本方程，它是研究反应堆理论的重要工在介质内运动的基本方程，它是研究反应堆理论的重要工具和基础。具和基础。3.1 单能中子扩散方程单能中子扩散方程 中子的
4、扩散和气体分子的扩散很相似，中子的扩散和气体分子的扩散很相似，它们都从浓度高的区域向浓度底的区域它们都从浓度高的区域向浓度底的区域扩散，扩散的速率与粒子的密度的梯度扩散，扩散的速率与粒子的密度的梯度成正比，既都服从成正比，既都服从“斐克扩散定律斐克扩散定律”。由于在热堆中子密度（由于在热堆中子密度（1016/m3）比介质）比介质的原子核密度（的原子核密度（ 1028/m3 ）小很多，因）小很多，因此它与气体分子的扩散又有不同，主要此它与气体分子的扩散又有不同，主要区别在于：区别在于：分子扩散是由于分子间的分子扩散是由于分子间的碰撞引起，而中子的扩散主要是由中子碰撞引起，而中子的扩散主要是由中子
5、与原子核之间碰撞的结果与原子核之间碰撞的结果，中子之间的，中子之间的相互碰撞可以忽略不计。相互碰撞可以忽略不计。中子与介质原子核中子与介质原子核 的散射碰撞的散射碰撞3.1 .1 斐克定律斐克定律 下面我们通过中子扩散过程来推导下面我们通过中子扩散过程来推导稳态情况下中子扩散方程，并假设：稳态情况下中子扩散方程，并假设：u介质是介质是 无限的、均匀的无限的、均匀的u在实验室坐标系中散射是各向同性在实验室坐标系中散射是各向同性u介质的吸收截面很小即介质的吸收截面很小即asu中子通量密度是随时间位置缓慢中子通量密度是随时间位置缓慢 变化的函数变化的函数 设在设在r处的体积元处的体积元 内中子通量密
6、度为内中子通量密度为?(r),每秒发生散射的中子数目为每秒发生散射的中子数目为 ,每秒自体积元每秒自体积元内散射出来沿着内散射出来沿着方向未经碰撞到达方向未经碰撞到达dA上的中子数是上的中子数是cosdldAdV dVrs)(dAdlerlsscos)(41从从-到到0积分积分式中式中?(r)不是不是r的函数的函数, 是一个未知函数是一个未知函数,所以上述积分无法所以上述积分无法计算计算, 我们可以将我们可以将?(r)(r)按按r r的函数展开的函数展开这里这里 沿沿方向的方向倒数方向的方向倒数,可以表示如下可以表示如下:x, y, z为在在x, y, z轴的投影轴的投影,完成以上积分可得沿完
7、成以上积分可得沿方向每秒穿过方向每秒穿过dA上的中子数为：上的中子数为：dAdlerdAlsscos)(40 dldlrr)()(dldzyxdldzzdldyydldxxdldzyxcossinxsinsinycoszrrsdldrdA1)(cos4对对 的半空间积分的半空间积分,就可以得到每秒沿就可以得到每秒沿z轴正方向自下轴正方向自下而上穿过而上穿过dA的中子数的中子数 。完成积分可得：完成积分可得：对对 的半空间积分的半空间积分,就可以得到每秒沿就可以得到每秒沿z轴负方向自上轴负方向自上而下穿过而下穿过dA的中子数的中子数 。0zeddzyxrdAddldrdAdAJzyxsrrsez
8、zsincos)(1)(41)(cos420200)(zJzrrrJsz)(614)()(0zezrrrJsz)(614)()(zJ单位时间内沿着单位时间内沿着z方向穿过方向穿过dA平面单位面积的平面单位面积的净中子数净中子数Jz为为 叫做叫做z方向的中子流密度或净中子流密度，若方向的中子流密度或净中子流密度，若dA的取向的取向与与x轴垂直，沿着轴垂直，沿着x方向穿过方向穿过dA平面单位面积平面单位面积净中子数净中子数Jx为为同样，沿着同样，沿着y方向穿过方向穿过dA平面单位面积平面单位面积净中子数净中子数Jy为为如果所讨论的面元并不垂直于任何坐标轴，那么单位时间内如果所讨论的面元并不垂直于任
9、何坐标轴，那么单位时间内穿过穿过dA平面单位面积平面单位面积净中子数净中子数J为三个分量之和为三个分量之和zrrJrJrJszzz)(3)()()()(rJzxrrJsx)(3)(yrrJsy)(3)(cos)(cos)(cos)(3zryrxrJs可以把上式写成矢量形式即可以把上式写成矢量形式即式中式中 斐克定律斐克定律矢量矢量J称为中子流密度，称为中子流密度，Jx ，Jy， Jz 是它在是它在 x,y,z 轴上的投影，轴上的投影，它表示空间任何一个点上中子宏观净流动的方向和梯度。它表示空间任何一个点上中子宏观净流动的方向和梯度。强调：强调：J即不同于中子束强度即不同于中子束强度 I，也不同
10、于中子通量密度，也不同于中子通量密度?(r,)。 它是由许多具有不同方向的微分中子束它是由许多具有不同方向的微分中子束矢量合成的量，矢量合成的量， 表示该处中子的净流动情况情况。它与中子通量密度表示该处中子的净流动情况情况。它与中子通量密度 ?(r,(r,) )的关系为的关系为斐克定律斐克定律表示：中子流密度表示：中子流密度J正比于负的中子通量密度梯度，正比于负的中子通量密度梯度， 其比例常数叫作其比例常数叫作扩散系数扩散系数，并用，并用D表示。斐克定律可表示。斐克定律可 写成写成nJJnkjincoscoscosgradkJjJiJJszyx3drrJ4),()(DgradJ3sD推导过程中
11、使用了推导过程中使用了在实验室坐标系中中子的散射是各向在实验室坐标系中中子的散射是各向同性同性的假设，实际计算中应对散射的各向异性进行修正，的假设，实际计算中应对散射的各向异性进行修正，必须用输运的平均自由程必须用输运的平均自由程 tr代替散射平均自由程代替散射平均自由程s，扩，扩散系数散系数D可写为可写为 为平均散射角余弦。为平均散射角余弦。斐可定律表明斐可定律表明： 任一处净中子流动的方向与中子通量密任一处净中子流动的方向与中子通量密度分布的梯度的方向相反。度分布的梯度的方向相反。grad? 的方向指向的方向指向?的增加方的增加方向，所以向，所以 J的方向指向的方向指向?减少最快的方向。减
12、少最快的方向。3trD01strA3203.1 .2 单能中子扩散方程的建立单能中子扩散方程的建立 核反应堆理论所基于的一个基本原理就是核反应堆理论所基于的一个基本原理就是”中子数守中子数守恒恒”,即在即在一定的体积内一定的体积内,中子数对时间的变化率应等于该体积中中子数对时间的变化率应等于该体积中子的产生率减去该体积内中子的吸收率和泄露率子的产生率减去该体积内中子的吸收率和泄露率. 中子数的中子数的守恒方程可以表达为守恒方程可以表达为中子的扩散方程就是基于这一平衡原理建立的。中子的扩散方程就是基于这一平衡原理建立的。u泄露率泄露率 利用高斯散度公式利用高斯散度公式)()()(),(ALSdV
13、trndtdV吸收率泄露率产生率SdSntrJ),(泄露率VVSdVtrJdivdVtrJdSntrJ),(),(),(泄露率u产生率产生率 设中子源分布函数用设中子源分布函数用S(r,t)表示，在体积表示，在体积V内中子产生率内中子产生率u吸收率吸收率 在体积在体积V内中子吸收率内中子吸收率中子数的守恒方程可以表达为中子数的守恒方程可以表达为去掉等式两边的积分可得去掉等式两边的积分可得方程叫做方程叫做连续方程，连续方程，在反应堆理论计算中具有非常主要的在反应堆理论计算中具有非常主要的地位。无论斐可定律是否适用，该方程都是普遍成立。地位。无论斐可定律是否适用，该方程都是普遍成立。VdVtrS)
14、,(产生率VadVtr ),(吸收率VVaVVdVtrdivJdVtrdVtrSdVtrndtd),(),(),(),(),(),(),(),(trdivJtrtrSttrna利用利用可得可得在斐可定律成立的基础上，连续方程可以写为：在斐可定律成立的基础上，连续方程可以写为：这是这是单能的中子扩散方程单能的中子扩散方程，如中子通量密度不随时间变化，如中子通量密度不随时间变化，上式就变为：上式就变为：称为稳态单能的中子扩散方程，称为稳态单能的中子扩散方程，这个方程是以斐可定律为基础这个方程是以斐可定律为基础得到，它的应用受到斐可定律适用范围的限制，仅适用于单能得到，它的应用受到斐可定律适用范围的
15、限制，仅适用于单能中子情况。中子情况。 是拉普拉斯算符，在不同坐标系的表示式为：是拉普拉斯算符，在不同坐标系的表示式为：zDzyDyxDxDgraddivdivJ)(2222222DzyxD泄露率),(),(),(),(12trtrDtrSttrva0)()()(2rSrrDa22222222 zyx直角坐标系2222222211 zrrrr柱坐标系22222222222sin1cot112 rrrrrr球坐标系3.1 .3 中子扩散方程的边界条件中子扩散方程的边界条件 必须用边界条件来确定扩散方程的解中的任意积分常数必须用边界条件来确定扩散方程的解中的任意积分常数,边界条件的数目应恰好使方程
16、由唯一的解。解扩散方程常用边界条件的数目应恰好使方程由唯一的解。解扩散方程常用的边界条件有：的边界条件有：u扩散方程适用范围，中子通量密度必须是正的、有限实数扩散方程适用范围，中子通量密度必须是正的、有限实数u在两种不同扩散性质的介质交界面上，垂直于分界面的中在两种不同扩散性质的介质交界面上，垂直于分界面的中子流密度和中子通量密度相等。子流密度和中子通量密度相等。两式相加减得扩散方程的边界条件：两式相加减得扩散方程的边界条件：BxAxJJBxAxJJBBAAdxdDdxdDBA在两种介质的分界面上的中子扩散在两种介质的分界面上的中子扩散u介质与真空交界外表面上从介质与真空交界外表面上从 真空返
17、回介质的中子流等于零真空返回介质的中子流等于零即即或或假设从交界面处将中子通量密度假设从交界面处将中子通量密度的分布曲线按它在交界面处的斜率向真空作直线外推，则在的分布曲线按它在交界面处的斜率向真空作直线外推，则在离开交界面距离离开交界面距离d处的位置上中子通量密度为零，处的位置上中子通量密度为零， 我们有我们有d称为称为直线外推距离直线外推距离00 xxJ064000 xtrxxdxdJtrxdxd2300ddxdx00trd32 应用输运理论和扩散理论的应用输运理论和扩散理论的外推距离求得的扩散方程的解外推距离求得的扩散方程的解以上以上d值是不准确的，因为值是不准确的，因为d值是根据扩散定
18、律推导而来，而值是根据扩散定律推导而来，而扩散定律不适用于真空交界处。更精确的中子输运理论所得扩散定律不适用于真空交界处。更精确的中子输运理论所得到的平面到的平面d值为值为 。在自由外表面的边界条件可以用。在自由外表面的边界条件可以用更简单的形式表示：更简单的形式表示：在自由表面外推距离在自由表面外推距离d处，中子通量密度为零。处，中子通量密度为零。trd7104. 03.1 .4 斐克定律和扩散理论的适用范围斐克定律和扩散理论的适用范围 在推导斐克定律时，我们做了一些假设，所以斐克定律在推导斐克定律时，我们做了一些假设，所以斐克定律 的应用范围是有限制的。的应用范围是有限制的。u假定了扩散介
19、质是无限的假定了扩散介质是无限的 在有限的介质内，在距离其表面几个自由程以外的全部在有限的介质内，在距离其表面几个自由程以外的全部 区域斐克定律时成立的，而在距真空边界两三个自由程区域斐克定律时成立的，而在距真空边界两三个自由程以内区域，它是不适用的。以内区域，它是不适用的。u推导中中子通量密度展成泰勒级数并只取到了一级项推导中中子通量密度展成泰勒级数并只取到了一级项 这要求在所讨论点的几个平均自由程内，中子通量密度这要求在所讨论点的几个平均自由程内，中子通量密度 必须缓慢变化或它的梯度变化不大。必须缓慢变化或它的梯度变化不大。 在控制棒附近或两在控制棒附近或两 种扩散性质明显不同的介质交界面
20、附近的几个平均自由程种扩散性质明显不同的介质交界面附近的几个平均自由程内，斐克定律不适用。此外，斐克定律只适用于内，斐克定律不适用。此外，斐克定律只适用于as弱弱吸收介质。吸收介质。 u推导中并没有考虑中子源的贡献，中子流密度的贡献只推导中并没有考虑中子源的贡献，中子流密度的贡献只 是来自中子与介质核的散射碰撞是来自中子与介质核的散射碰撞 在强中子源两三个平均自由程的区域内，斐克定律不在强中子源两三个平均自由程的区域内，斐克定律不 适用。适用。3.2 非增殖介质内中子扩散方程的解非增殖介质内中子扩散方程的解稳态单能的中子扩散方程稳态单能的中子扩散方程无源情况下，即除中子源所在的位置以外的无源区
21、域无源情况下，即除中子源所在的位置以外的无源区域,扩散扩散方程有以下形式方程有以下形式: 或或L称为中子的扩散长度称为中子的扩散长度,它表征中子在介质中扩散特征的一它表征中子在介质中扩散特征的一个重要的量。以上方程称为波动方程或亥姆霍兹方程，个重要的量。以上方程称为波动方程或亥姆霍兹方程， 加加上适当的边界条件就可以得出以上数理方程的解。下面列上适当的边界条件就可以得出以上数理方程的解。下面列出一些常见的简单几何形状下波动方程的普遍解。出一些常见的简单几何形状下波动方程的普遍解。0)()()(2rSrrDa0)()(22LrraDkL2210)()(22rr解的形式一维平板球 或一维圆柱22B
22、2B222ddx sincosABxCBxAeCesincoshBxBxABxCBx或222d2 dddrrr sincosBrBrACrreeBrBrACrrsinhcoshBrBrACrr222d1 dddrr r 00()()AI BrCY Br00AI () CK ()BrBr注：分别是第一类和第二类零阶贝塞尔函数；注：分别是第一类和第二类零阶贝塞尔函数； 分别是第一类和第二类零阶修正贝塞尔函数。分别是第一类和第二类零阶修正贝塞尔函数。 （见附表（见附表8）。）。在一些几何形状情况下波动方程在一些几何形状情况下波动方程 的解的解下面我们讨论几种特殊情况下扩散方程的解，它可以帮助下面我们
23、讨论几种特殊情况下扩散方程的解，它可以帮助我们掌握和熟悉扩散方程的求解和如何使用边界条件我们掌握和熟悉扩散方程的求解和如何使用边界条件：u无限介质内点源的情况无限介质内点源的情况 在介质中有一个每秒各向同性放射出在介质中有一个每秒各向同性放射出S个中的点源，采用个中的点源，采用 球坐标，原点选择在点源上。球对称的扩散方程为：球坐标，原点选择在点源上。球对称的扩散方程为：这个方程在这个方程在r=0处不成立，其边界条件为：处不成立，其边界条件为：（1）除）除r=0处以外，中子通量在各处均为有限值；处以外，中子通量在各处均为有限值；（2）中子源条件：）中子源条件： 引入新变量引入新变量 ，代入扩散方
24、程可将扩散方程化为：，代入扩散方程可将扩散方程化为： )0(, 0)()(2222rLrdrrdrdrdSrJrr)(4lim20ru 0222Ludrud方程的解为：方程的解为：所以：所以：C=0, 所以所以 由由根据中子源边界条件：根据中子源边界条件：得到得到最后得到无限介质内的中子通量密度为：最后得到无限介质内的中子通量密度为：LrLrCeAeu/reCreArLrLr/)(reArLr /)(LrerrLDAdrrdDrJ/2)11()()(SeLrDArJrLrrr/020)1(4lim)(4limDSA404)(/rDrSerLru无限平面源位于有限厚度无限平面源位于有限厚度 介质
25、内的情况介质内的情况设源为强度为设源为强度为S的平面中子源，的平面中子源，扩散方程为扩散方程为边界条件为：边界条件为：（1）（2）中子源条件）中子源条件当当x为正值时，扩散方程的解为：为正值时，扩散方程的解为：由边界条件（由边界条件（1）可得：）可得：00)()(222xLxdxxd0)2/(,)2/(aax时当2/)(lim0SxJxLxLxCeAex/)(LaAeC/平面源位于有限厚介质的情况平面源位于有限厚介质的情况通量密度可以表达为：通量密度可以表达为：根据边界条件（根据边界条件（2）可以得到：）可以得到：中子通量密度的解为：中子通量密度的解为：由于对称性，用由于对称性，用|x|代替代
26、替x可得到对所有可得到对所有x适用的中子通量密适用的中子通量密度的解度的解用用 乘分子和分母，并利用双曲函数性质乘分子和分母，并利用双曲函数性质可得：可得：)(/ )(/LxaLxeeAx1/)1(2LaeDSLALaLxaLxeeeDSLx/ )(/12)(LaLxaLxeeeDSLx/|)|(/ |12)()2/cosh(2/|)|2sinh(2)(LaLxaDSLxLae2)(21sinhuueeu)(21coshuueeu通过实际的边界向外泄露的中子流密度等于通过实际的边界向外泄露的中子流密度等于对于无限介质平面源情况，对于无限介质平面源情况，a，有，有我们我们可以把扩散长度看作中子通
27、量密度的衰减长度可以把扩散长度看作中子通量密度的衰减长度，由图中，由图中可以看出当介质厚度为扩散长度的三倍时，除在边界附近，可以看出当介质厚度为扩散长度的三倍时，除在边界附近，中子通量密度的分布与无限介中子通量密度的分布与无限介质内的分布相差不多。质内的分布相差不多。对于单能的情况，反射层厚度对于单能的情况，反射层厚度大于三个扩散长度时，其效果大于三个扩散长度时，其效果就大致和无限厚度相当。因此，就大致和无限厚度相当。因此，没有必要使用过厚的反射层。没有必要使用过厚的反射层。)2/cosh(2)/cosh(2/LaLdSdxdDJdaLxeDSLx/ |2)(不同厚度介质内的中子通量密度分布不
28、同厚度介质内的中子通量密度分布u包含两种不同介质的情况包含两种不同介质的情况 在不同介质的交界面上扩散在不同介质的交界面上扩散方程必须满足交界面边界条件方程必须满足交界面边界条件边界条件：边界条件：x为正值时，扩散方程的解是：为正值时，扩散方程的解是： 和和2|, 0)(1)(0,2|, 0)(1)(222222121212axxLdxxdxaxxLdxxd2/)(lim)2( |102SxJxxx中子源条件：）趋近于零。（时，）当（2/222/1121)4()2/()2/()3(axaxdxdDdxdDaa)/sinh()/cosh(11111LxCLxA22/2/22LxLxeCeA双区介
29、质内中子通量密度分布双区介质内中子通量密度分布由边界条件（由边界条件（1）可得）可得C2=0，边界条件（边界条件（2）可得：）可得：由边界条件（由边界条件（3）和（和（4）可得：）可得：图中虚线部分代表的是没有介质图中虚线部分代表的是没有介质2时，中子通量密度的分布时，中子通量密度的分布。1112DSLC)2/sinh()2/cosh()2/sinh()2/cosh(2121112112121111LaLDLaLDLaLDLaLDDSLA11211122212)2/sinh()2/cosh()2exp(2LaLDLaLDLaLSLA双区介质内中子通量密度分布双区介质内中子通量密度分布3.3*
30、反照率反照率 介质介质 A 介质介质 B J+当平板介质外再围上一层扩散介质后，当平板介质外再围上一层扩散介质后，中子通量密度分布的下降将比于真空交中子通量密度分布的下降将比于真空交 J-界时减缓许多。这就是堆芯使用反射界时减缓许多。这就是堆芯使用反射层的原因。层的原因。反射层的效率可以通过反反射层的效率可以通过反射系数或反照率表示射系数或反照率表示：根据扩散定律，反照率可写为：根据扩散定律，反照率可写为：通常反照率采用反射介质的性质来表示。反照率不仅取决于通常反照率采用反射介质的性质来表示。反照率不仅取决于反射介质的材料特征，而且还取决于系统尺寸和几何形状。反射介质的材料特征，而且还取决于系
31、统尺寸和几何形状。 JJdxdDdxdDdxdDdxdDJJ21212424对于对于无限平板反射层无限平板反射层，反照率等于反照率等于对于对于有限厚度的反射层有限厚度的反射层a时，时，。反照率的重要应用在于反照率的重要应用在于用来作为于反射层介质相邻的分界面用来作为于反射层介质相邻的分界面上的边界条件，以代替反射层介质上的边界条件，以代替反射层介质。如果能精确知道堆芯水。如果能精确知道堆芯水反射层的反照率，在作芯部计算时可以在芯部于反射层上应反射层的反照率，在作芯部计算时可以在芯部于反射层上应用下列边界条件以代替反射层：用下列边界条件以代替反射层：这样，就不必对反射层部分进行这样，就不必对反射
32、层部分进行计算计算，从而节省大量计算时间。，从而节省大量计算时间。LxCeLDLD2121)coth(21)coth(21LaLDLaLDJJ3.4 扩散长度、慢化长度和徙动长度扩散长度、慢化长度和徙动长度u扩散长度扩散长度 大多数元素散射截面与能量无关，而吸收截面服从大多数元素散射截面与能量无关，而吸收截面服从1/v 律，当热中子能谱按麦克斯韦分布时，热中子吸收截面等于律，当热中子能谱按麦克斯韦分布时，热中子吸收截面等于a,0 是能量为是能量为 En =0.0253 eV 的中子吸收截面，的中子吸收截面，Tn 为中为中子温度，子温度，ga 是非是非 1/v 修正因子，修正因子，代入上式代入上
33、式 32traaDLanaagT29320,2)(1293202nanTgTLL 为了阐明扩散长度的物理意义，我们计算热中子从产生地点为了阐明扩散长度的物理意义，我们计算热中子从产生地点到被吸收地点穿行距离的均方值到被吸收地点穿行距离的均方值 对于无限介质中的点源对于无限介质中的点源，在球，在球壳内每秒被吸收的中子数是壳内每秒被吸收的中子数是所以均方值所以均方值 （空间二次距）可以（空间二次距）可以表示成表示成将点源的中子通量密度代入可得：将点源的中子通量密度代入可得： 或或 对于对于平面源的情况有平面源的情况有2rdrrra)(422r020222)(4)(4(drrrdrrrrraa200
34、326/LerdrerrLrLr2261rL2221xL点源空间二次矩的计算点源空间二次矩的计算 从计算可以看出，扩散长度从计算可以看出，扩散长度L的大小直接影响堆内热中子的大小直接影响堆内热中子的泄露。的泄露。 L愈大，则热中子自产生地点到被吸收地点所移愈大，则热中子自产生地点到被吸收地点所移动的直线平均距离也愈大，因而热中子泄露到反应堆外的动的直线平均距离也愈大，因而热中子泄露到反应堆外的概率也就愈大概率也就愈大。u慢化长度慢化长度 我们还希望能计算出中子我们还希望能计算出中子在介质中从产生地（快中子）在介质中从产生地（快中子）到慢化成热中子时所穿行的到慢化成热中子时所穿行的直线距离。这与
35、堆内中子的慢直线距离。这与堆内中子的慢化过程中的泄露有关，同样考虑一个点源的情况，源中子化过程中的泄露有关，同样考虑一个点源的情况，源中子能量为能量为E0 ，我们把，我们把E0 到到 Eth 的中子称为快群中子。的中子称为快群中子。0),()(1EEthdEErx慢化长度计算慢化长度计算而把而把Eth 以下的中子称为热群中子，同时定义一个移出截面以下的中子称为热群中子，同时定义一个移出截面1 使使 设设s 为快中子的宏观散射截面，为快中子的宏观散射截面， s ?1 便是每秒单位体积便是每秒单位体积内快中子发生的碰撞数内快中子发生的碰撞数，因此一个源中子从初始能量，因此一个源中子从初始能量E0
36、降降低到低到 Eth 平均所需要的碰撞次数为平均所需要的碰撞次数为因此因此快群中子转移到热群中子转移率快群中子转移到热群中子转移率为为 得到得到以下的中子数到处每秒单位体积内减速在thErr )(11/ )/ln(0thsEEN thsEE0111ln1thsEE01ln我们可以求出无限介质点源情况下我们可以求出无限介质点源情况下快群中子的扩散方程快群中子的扩散方程该方程类似于方程（该方程类似于方程（3-48），）， L1 称为称为慢化长度慢化长度，它具有，它具有长度的量刚。反应堆物理中长度的量刚。反应堆物理中L21 称为中子年龄，用称为中子年龄，用th 表示。表示。 即为即为慢化长度慢化长度
37、。中子的年龄。中子的年龄 (E)定义为定义为当当E=Eth ,(E)便等于热中子年龄便等于热中子年龄th ， th 是随着中子能量是随着中子能量降低或中子慢化时间的增大而增大的函数，它有年龄的意义。降低或中子慢化时间的增大而增大的函数，它有年龄的意义。0121211121112ln3100EEDLLDthtrathEdEEEDdEdEEEDEsEEs)()()()()(0但它并不具有时间的量刚，而具有长度平方的量刚。因为能但它并不具有时间的量刚，而具有长度平方的量刚。因为能量愈低则中子离开源点的距离愈大。量愈低则中子离开源点的距离愈大。 可以证明慢化长度平方或热中子年龄和扩散长度的平方可以证明
38、慢化长度平方或热中子年龄和扩散长度的平方具有相似的物理意义，具有相似的物理意义，即慢化长度的平方即慢化长度的平方L21或热中子年龄或热中子年龄th 等于在无限长度介质内中子自源点产生发出在介质中慢化到等于在无限长度介质内中子自源点产生发出在介质中慢化到年龄年龄th（Eth）时所穿行的直线距离的均方值的六分之一）时所穿行的直线距离的均方值的六分之一，它具有长度平方的量纲。它具有长度平方的量纲。u徙动长度徙动长度 反应堆计算中经常用到的量反应堆计算中经常用到的量徙动面积徙动面积定义为：定义为：M称为徙动长度称为徙动长度。从热中子年龄和扩散长度的意义由上式可得：。从热中子年龄和扩散长度的意义由上式可
39、得：thLM22)(61222dsrrMdsMrrr两边取均方值两边取均方值由于由于rs 和和 rd 的方向彼此不相的方向彼此不相关，因而两者夹角的余弦的关，因而两者夹角的余弦的平均值等于零，所以平均值等于零，所以徙动面积徙动面积M2 是中子由裂变产生的快中子变成热中子并在介质是中子由裂变产生的快中子变成热中子并在介质中扩散最终被吸收所穿行直线距离均方值的六分之一。中扩散最终被吸收所穿行直线距离均方值的六分之一。 徙动长度徙动长度M是影响芯部中子泄露程度的重要参数，是影响芯部中子泄露程度的重要参数， M值值愈大，则中子愈容易泄露愈大，则中子愈容易泄露。cos2222sddsMrrrrr222261)(61MdsrrrM徙动长度的计算徙动长度的计算
