核反应堆物理分析第2章课件

1、第二章中子慢化和慢化能谱 反应堆内裂变中子的平均能量为2 MeV。由于中子散射碰撞而降低速度的过程成为慢化过程。热堆内，弹性散射对慢化过程起主要作用。在慢化过程，热堆内中子密度按能量具有稳定的分布，称之为中子慢化能谱 2.1中子的弹性散射过程 2.1.1 弹性散射时能量的变化我们可以求得在质心的速度VCM设靶核静止，则在碰撞前质心系中中子和靶核的速度为：可以看出在质心系内，中子与靶核的总动量为零：碰撞前后动量和动能守恒：由此我们可以解得：由此可以看出在质心系里，碰撞前后中子和靶核的速度大小不变，只是运动的方向发生了变化。在实验室系里，碰撞后中子的速度为：由余弦定律可得：将已知各量代入可得：在实

2、验室系碰撞前后中子能量之比为：让可得：有以上结果可以看出：（1） 碰撞前后中子能量没有损失。（2） 一次碰撞中中子的最大能量损失为（3）中子在一次碰撞中损失的最大能量与靶核的质量有关。 A=1，则 =0，Emin=0，即中子与氢碰撞后能量全部 损失掉。 A=235，则 =0.983，Emin=0.02E，即中子与235U碰撞 后能量最大损失约为碰撞前中子能量的2%。所以应该 选择轻核元素作为慢化剂。为了获得实验室系和质心散射角之间的关系，由图2-2：利用（2-13）式代替cosc ， 可得到实验室系中散射角和碰撞前后中子能量的关系2.1.2 散射后中子能量的分布中子的能量变化与散射角度之间有对

3、应关系，根据碰撞后中子散射角分布的概率可以求得碰撞后中子能量的分布概率。 我们有关系式： (2-18)实验和理论计算（量子力学）表明，对一般的轻元素，当能量E 小于几个MeV时，在质心系内中子的散射是各向同性的，即碰撞后中子在任一立体角内出现的概率相等。在这种情况下，一个中子被散射到立体角dc内的概率：因而积分可得 (2-19)由（2-13）式微分可得 (2-20)由（2-19）、（2-20）和（2-18）可得散射后能量是均匀分布的，与碰撞后能量无关。并且2.1.3 平均对数能降 对数能降 u 定义为： 或 E=E0e-u其中 E0为选定的参考能量，一般选 E0=2 MeV 或 E0=10 M

4、eV。随中子能量的减小，中子的对数能降在增大，其变化与能量相反。一次碰撞后对数能降的增加量为：由（2-14）式可知，一次碰撞最大的对数能降为在研究中子的慢化过程时，有一个常用的量，就是每次碰撞中子能量的自然对数的平均变化值，叫做平均对数能降在质心系内各向同性的情况下：积分后可得：当 A 10 如用Nc 表示中子从能量E1 慢化到能量E2平均碰撞次数，则使中子能量由2 MeV慢化到0.0253 eV时分别所需要的与H核、石墨核以及235U核的平均碰撞次数为:2.1.4 平均散射角余弦在质心系中中子每次碰撞平均散射角余弦为：这是预期结果，因为在质心系中中子散射是各向同性。在实验室系中中子每次碰撞平

5、均散射角余弦为：由于中子在实验室系和质心系中有对应关系，因此由（2-16）和（2-19）可得因而，尽管在质心系是各向同性的，但在实验室系确是各向异性的，而且在实验室系中子散射后沿它原来运动方向的概率较大。平均散射角余弦的大小表示了各向异性的程度。在实验室系平均散射角余弦随着靶核质量数的减小而增大，靶核的质量越小，中子散射后各向异性（向前运动）的概率就越大。2.1.5 慢化剂的选择反应堆中要求慢化剂具有较大宏观散射截面s和平均对数能降 。通常把乘积s叫做慢化剂的慢化能力。我们还要求慢化剂有较小的吸收截面，定义s / a叫做慢化比。慢化剂慢化能力 s /m-1慢化比 s / aH2OD2OBe石墨

6、1.5310-21.7710-31.610-36.310-47021001501702.1.6 中子的平均寿命在无限介质中，裂变中子慢化到热中子所需要的平均时间称为慢化时间。中子在时间dt内与原子核发生的碰撞数为在dt时间里对数能降的增量等于n,即 或于是，中子由裂变能慢化到热能的慢化时间为用平均值s来替代s(E)，对上式积分可得ts的估计值ts 一般在10-4 到10-6秒量级。介质中的热中子在自产生至被俘获以前所经历的平均时间，称为扩散时间，热中子的平均寿命。平均寿命为：对于吸收截面满足1/v规律的介质，有a(E)v=a0v0 式中a0是当v0 =2200m/s 时的热中子宏观截面。上式表

7、明对于1/v介质热中子的平均寿命与中子能量无关。 Table 2-2 几种慢化剂的慢化和扩散时间 快中子自裂变产生到慢化成为热中子，直到最后被俘获的平均时间，称为中子的平均寿命。慢化剂慢化时间/s扩散时间/sH2OD2OBeBeO石墨6.310-65.110-55.810-57.510-51.410-41.410-40.1373.8910-36.7110-31.6710-22.2 无限均匀介质内中子的慢化能谱慢化密度：它的定义为在r处每秒每单位体积内慢化到能量E以下的中子数， 用q（r，E）表示。 r处能量为E的中子每秒发生散射的次数为： s(r, E)(r, E)，而散射函数f(EE) 表示

8、能量为E的中子散射后能量变为E的概率，因而在r处每秒每单位体积内能量为E的中子慢化到能量E以下的中子数为 q（r，E）应等于E E 的所有能量中子慢化到E以下的中子数目的总和，也就是对E 的积分，即将散射函数用（2-21）表示，并对上式积分可得： 慢化密度q（r，E）给出了r处中子被慢化并通过给定定能量E的慢化率。下面讨论无限介质内慢化方程。对于无限大介质，中子通量密度只与能量有关，与坐标r无关。散射到能量E附近dE能量间隔内的中子由两部分组成： 中子源产生直接进入该能量间隔 的中子。用S（E）dE 表示， S（E）为中子源强分布函数。 中子与介质原子核的散射结果。 根据中子平衡的稳定条件，单

9、位时间和体积内，散射到能量微元dE内的中子数和源中子数之和应该等于从这个能量微元散射出去和被吸收的中子数总数。 稳态无限介质内的中子慢化方程为： 方程的解便是中子慢化能谱。无吸收单核素无限介质情况 最简单的情况： 只含有一种核素的无吸收介质。 中子源S(E0) 为均匀分布。我们只讨论慢化区（1 eV 0.1 MeV）内的弹性散射慢化问题，此区不包括由于裂变反应直接产生 中子源，这时散射在质心系内是各向同性。慢化方程可写为：EE0时可以证明它的渐进解的形式为为了确定常数C，把（2-44）代入（2-41）这里 是平均对数能降。因而在渐进情况下，慢化能谱为：对无吸收纯氢介质上式便是慢化方程的严格解。

10、对无吸收情况，单能源，q(E)=S0 ，上式变为无吸收混合物无限介质情况 对于混合物s(E)=Nisi ，其中Ni为混合物中i元素的核子数。这时，慢化方程（2-42）对于混合物介质变为这里i表示混合物的组分，如果所有元素截面都等于常数，则上式的渐进解为（E E0 ）i 为中子与i种元素碰撞的平均对数能降。定义混合物的平均对数能降为 （2-50）中子通量密度的慢化能谱分布可表示为： （2-51）它和单核素的慢化能谱分布（2-46）一样，只是用混合物平均对数能降 来代替。（2-50）可以用来计算混合物的平均对数能降。 （2-51）式就是无限无吸收介质内中子慢化能谱分布。它表明无限无吸收介质内在慢化

11、区内服从1/E分布或称之为费米谱分布，我们常把它作为反应堆内中子能谱的近似。无限介质弱吸收情况 对于无限介质弱吸收情况，即as ， 设在能量间隔E 到 （E-dE）范围内，慢化密度由于中子被吸收减小了dq，它等于在dE内被吸收的中子数，因此在弱吸收情况，即as 我们近似认为(E)基本上和无吸收情况下相同，这样上式变为将上式从E 到 E0积分，同时利用q(E0)=S0， 有根据逃脱共振俘获概率p（E）的定义，有2.3 均匀介质中的共振吸收当中子能量慢化到100 keV以下中能区，反应堆内的很多重要的材料如U， Pu， Th 等多表现出强烈的共振吸收特征，具有很高并且很密集的共振峰。在慢化过程中必

12、然有一部分中子被共振吸收。共振吸收对反应堆内的链式反应过程有非常重要的影响。上一部分讨论的是中子在介质具有弱吸收下的慢化过程，这一部分我们将讨论中子在有强共振吸收下的慢化过程。2.3.1 均匀介质中的有效共振积分及逃脱共振俘获概率先以 一个理想的简单情况下的共振吸收为例说明共振吸收的基本物理特征。设有一无限大均匀介质，介质中均匀分布着中子源，每秒单位体积内放出S0个能量为E0的快中子。 这时中子通量密度与空间无关只是能量的函数。同时假设强共振峰不但可分辨而且峰与峰之间的间距足够大。这种情况下，到达共振峰i前的慢化中子通量密度可以用（2-46）或（2-51）式来表示。为了方便我们进一步假设中子的

13、源强为 S0=s ，到达共振峰i前的慢化中子通量密度可以简化为：共振峰在这种情况下（源强为 S0=s ）的吸收反应率为NA 为单位体积内共振吸收剂的核子数，i为共振峰宽度。称作共振峰i有效共振积分。也既吸收反应率为因而逃脱共振俘获概率等于两边取对数并利用 可得到裂变中子慢化到热中子的逃脱共振俘获概率p等于所有共振峰的pi的乘积其中I为整个共振区的有效共振积分 根据实验结果可以得到不同浓度情况下238U 和 232Th的有效共振积分的经验公式：对于238U 在室温有对于232Th 在室温有2.4 热中子能谱和热中子平均截面2.4 .1 热中子能谱 在压水堆中通常将Ec=0.625 eV定义为分界

14、能或缝合能, Ec能量以下的中子称为热中子。所谓热中子是指中子与所在的介质的原子或分子处于热平衡状态的中子。处于热平衡状态的热中子，它们的能量分布也服从麦克斯韦-波耳兹曼分布，即实际上，热中子的能谱分布与介质原子核的麦克斯韦并不完全相同。因为：在反应堆中，所有的热中子都是从高能慢化而来，然后与介质达到热平衡，这样子较高能区的中子数就较多。由于介质也要吸收中子，因此必然有一部分中子还没有慢化成热中子以前就被介质吸收了，其结果又造成了能量较低部分的中子份额减少，高能中子的份额较大。 这一现象称为热中子能谱的“硬化”。 精确计算热中子能谱是比较复杂问题，因为在处理能量低于1电子伏的中子与慢化剂核的散

15、射时已不能把慢化剂核看成静止的，自由的，必须考虑到慢化剂核的热运动等因素。 在实际计算中，可以近似认为热中子能谱仍然具有麦克斯韦的分布的形式只是热中子最概然能量En=kTn/2比介质原子核的最概然Em=kTm/2要高。这相当于把介质的麦克斯韦分布谱向右移动，使Tm 增大到Tn。 Tn称为中子温度，中子温度的数值一般要比介质温度高。 中子温度高于介质温度的差值将随着介质慢化能力的减少和吸收截面的增加而增大，Tn 与Tm的关系可近似地 用以下公式表达：s 为栅元或介质的慢化能力；a(kTm) 为中子能量等于kTm的栅元或介质的宏观吸收截面；Tm 为介质温度。假定栅元或介质内各元素核的吸收截面服从1/v率，则：对于一些弱吸收的纯慢化剂，中子温度可以用以下近似公式计算： 当 A25时，0 25时，其中热中子反应堆内中子能谱分布高能区（能量大于0.1 MeV）,中子能谱近似地可以用裂变中子谱来描述。在慢化区，中子能量密度的能谱近似按照1/E规律变化。