第4章-均匀反应堆临界理论2014

1、1 均匀反应堆的临界理论均匀反应堆的临界理论 主讲：张竞宇主讲：张竞宇 2 Contents 前言前言 均匀裸堆的单群扩散理论均匀裸堆的单群扩散理论 有反射层反应堆的单群扩散理论有反射层反应堆的单群扩散理论 功率分布展平概念功率分布展平概念 3 在上一章中我们讨论了中子在在上一章中我们讨论了中子在非增殖介质非增殖介质内扩散内扩散 的规律和中子扩散方程的解法。现在我们进一步的规律和中子扩散方程的解法。现在我们进一步 将其用于讨论由核燃料和慢化剂等组成的将其用于讨论由核燃料和慢化剂等组成的有限均有限均 匀增殖介质匀增殖介质（反应堆系统）内的中子扩散问题。（反应堆系统）内的中子扩散问题。 中心问题是

2、讨论反应堆的临界。中心问题是讨论反应堆的临界。 4 在反应堆临界理论中，主要研究两方面的问题：在反应堆临界理论中，主要研究两方面的问题： l 各种形状的反应堆达到临界状态的条件（临界条件）各种形状的反应堆达到临界状态的条件（临界条件）: e.g., 临界时系统的体积大小（临界体积）和燃料成分（富临界时系统的体积大小（临界体积）和燃料成分（富 集度）及其装载量（临界质量）集度）及其装载量（临界质量）。 l 临界状态下系统内中子通量密度（或功率）的空间分布。临界状态下系统内中子通量密度（或功率）的空间分布。 5 确定轮方法计算实确定轮方法计算实 际的反应堆系统际的反应堆系统 l物理过程与中子能量的

3、复杂物理过程与中子能量的复杂 依赖关系依赖关系 “ “分群理论分群理论” l几何的复杂性几何的复杂性 “ “几何近几何近 似似”处理处理 l材料的复杂性材料的复杂性 “ “均匀化均匀化” 处理处理 均匀反应堆：均匀反应堆：是指这样一种堆，其堆芯的各种材料（燃料、是指这样一种堆，其堆芯的各种材料（燃料、 慢化剂、结构材料等等）是均匀地混合在一起的。因此整慢化剂、结构材料等等）是均匀地混合在一起的。因此整 个堆芯的材料特性是一致的，核截面等数据都是一样的。个堆芯的材料特性是一致的，核截面等数据都是一样的。 6 均匀堆与非均匀堆均匀堆与非均匀堆 l世界上数以千计的反应堆中，只有一个名叫世界上数以千计

4、的反应堆中，只有一个名叫“水锅炉水锅炉” 的实验堆是均匀堆。其堆芯是硫酸铀酰的水溶液。的实验堆是均匀堆。其堆芯是硫酸铀酰的水溶液。 l其他的都是非均匀堆，堆芯中的燃料和慢化剂是分开其他的都是非均匀堆，堆芯中的燃料和慢化剂是分开 的，不混在一起。的，不混在一起。 既然如此，为何还要研究均匀反应堆？既然如此，为何还要研究均匀反应堆？ Why? 7 研究思路：从容易的着手，逐步精确化研究思路：从容易的着手，逐步精确化 1.均匀堆比较容易描述，建立的物理数学模型均匀堆比较容易描述，建立的物理数学模型 比较简单。但是，从中引出的基本概念有普遍比较简单。但是，从中引出的基本概念有普遍 应用价值。应用价值。

5、 2.工程设计中，对实际的非均匀反应堆进行分析工程设计中，对实际的非均匀反应堆进行分析 时，也要先进行时，也要先进行 “均匀化均匀化”，化为均匀堆。，化为均匀堆。 8 1.均匀裸堆的单群中子扩散方程的建立均匀裸堆的单群中子扩散方程的建立 2.均匀裸堆的单群扩散方程的解均匀裸堆的单群扩散方程的解 3.热中子反应堆的临界条件热中子反应堆的临界条件 4.几种几何形状裸堆的几何曲率和中子通量密度几种几何形状裸堆的几何曲率和中子通量密度 分布分布 5.反应堆曲率和临界计算任务反应堆曲率和临界计算任务 6.单群理论的修正单群理论的修正 均匀裸堆的单群理论均匀裸堆的单群理论 9 裸堆裸堆：无反射层的反应堆：

6、无反射层的反应堆 单群单群：全部中子都在一个能群里。或者可以进一：全部中子都在一个能群里。或者可以进一 步假设堆内里所有中子都是热中子。步假设堆内里所有中子都是热中子。 1 1. .单群理论的建立单群理论的建立 10 ),(),(),(),( ),(1 0 2 trStrktrtrD t tr aa 对于由燃料与慢化剂组成的均匀增殖介质反应堆系统，单位时对于由燃料与慢化剂组成的均匀增殖介质反应堆系统，单位时 间、单位体积内的间、单位体积内的裂变中子源强裂变中子源强为：为： ),(),(trtrS fF 根据无限介质增殖因子定义根据无限介质增殖因子定义 ),(),(trktrS aF 代入代入单

7、群中子扩散方程单群中子扩散方程可得可得 l D及及 a是对中子能谱平均后的数值；是对中子能谱平均后的数值； l 在反应堆运行初期，须考虑外源中子，大多数情况下忽略外中子，认为在反应堆运行初期，须考虑外源中子，大多数情况下忽略外中子，认为 裂变中子是反应堆内中子的唯一来源裂变中子是反应堆内中子的唯一来源 ? ? ? ? 11 0), 2 (), 2 (t a t a 2 2 1)( )( 1 )( )( L k dt tdT tTDx x )()(),(tTxtx )()0 ,( 0 xx a/2a/2 0 x 无限平板反应堆 ),(),(),( ),(1 2 txktxtxD t tx aa

8、(4-3) 无外源无限平板反应堆单群扩散方程无外源无限平板反应堆单群扩散方程 初始条件为初始条件为 (4-4) 边界条件为边界条件为 (4-5) ),( 1 ),( ),(1 2 2 tx L k tx t tx D (4-6) 由式(4-3)得 利用分离变量法求解分离变量法求解，方程具有如下形式的解： (4-7) 将(4-7)式代入(4-6)式 (4-8) 2 2. .均匀裸堆的单群扩散方程的解均匀裸堆的单群扩散方程的解 2 a D L 12 , 5 , 3 , 1n a n Bn , 3 , 2 , 1 ) 12( n a n Bn BxCBxAxsincos)( 0 2 cos Ba A

9、 2 2 )( )( B x x 上式两端必须等于某一常数，设为-B2，有 0)()( 22 xBx 或 (4-9) 波动方程(4-9)式的通解为 由于初始通量密度分布0(x)关于x=0平面对称，因此只能选择满足对称条件的解，即 BxAxcos)( 由边界条件(4-5)式可导出(x)满足如下的边界条件：(a/2)=0 因此要求或(4-10) x a n AxBAx nnnn ) 12( coscos)( 对于特征值特征值Bn，波动方程的解 n(x) 称为此问题的特征函数称为此问题的特征函数。 13 2 2222 / (1)(1) a n nn DL l DL BDL B nn ltk n Ce

10、T /) 1( n nn n l k dt tdT tT 1)( )( 1 2 2 1)( )( 1 n n n B L k dt tdT tTD nn ltk n n ex a n Atx /)1( 1 ) 12( cos),( 对应于每一个Bn值和n(x)，都有一个Tn(t)与之对应 该式可转换为 式中 22 1 n n BL k k (4-12) (4-13) (4-14) 方程(4-12)解为 其中C为待定常数。对于一维平板反应堆，其中子通量密度的完全解中子通量密度的完全解就是对n=1到 n=所有项的总和，即 (4-15) 14 1. 对于一定几何形状和体积的反应堆芯部，若对于一定几何

11、形状和体积的反应堆芯部，若B12对应的对应的k11，则，则(k1-1)为正值，中子通量密度为正值，中子通量密度 (x,t)将随时间不断增加，系将随时间不断增加，系 统处于超临界状态。统处于超临界状态。 3. 若调整堆芯尺寸或改变材料成分，使若调整堆芯尺寸或改变材料成分，使k1 =1，则其余，则其余( kn-1)都将为负值。都将为负值。 中子通量密度中子通量密度 (x,t)第一项将与时间无关，而其它各项将随时间而衰减。第一项将与时间无关，而其它各项将随时间而衰减。 当时间足够长时，当时间足够长时，n1各项将衰减到零，系统处于稳态，中子通量密度各项将衰减到零，系统处于稳态，中子通量密度 按基波形式

12、按基波形式(B=B1)分布，系统处于临界状态分布，系统处于临界状态。 , 3 , 2 , 1 ) 12( n a n Bn 22 1 n n BL k k nn ltk n n ex a n Atx /)1( 1 ) 12( cos),( 3 3. .热中子反应堆的临界条件热中子反应堆的临界条件 三种情况：三种情况： 15 0)()( 22 rBr g 1 1 2 1 2 1 BL k k 重要结论：重要结论： (1) 裸堆单群近似的临界条件为： (4-17) B12为波动方程的最小特征值，记为Bg2，称为特征曲率，称为特征曲率；k1为有效增殖因 子。 (2) 反应堆处于临界状态时，中子通量密

13、度按最小特征值中子通量密度按最小特征值Bg2对应的基波函对应的基波函 数分布数分布，也就是说，此时反应堆的中子通量密度空间分布满足波动方程 (4-18) 16 x a Ax cos)( 1 )(1 22 1 a L k k 无限平板反应堆的临界条件无限平板反应堆的临界条件为 (4-19) 若系统材料组成给定若系统材料组成给定，则只有一个唯一的尺寸a0能使k1=1，即为临界大小； 当aa0时，则k11，为超临界；当aa0时，k10，超临界； 0，次临界 | |表示反应堆偏离临界状态的程度。 反应性单位：反应性单位：PCM，1PCM=10-5； ；元， 元，1元反应性元反应性=1 eff （反应堆

14、动力学）（反应堆动力学） (4-49) 33 单群是一种非常近似的方法。对于热中子反应堆，直接以上临界条件有较大误差单群是一种非常近似的方法。对于热中子反应堆，直接以上临界条件有较大误差。 修正：修正：用用M2=L2+ 来替换上式中的来替换上式中的L2 1 1 22 1 g BL k k 2 2 1 L k Bm 1 1 22 1 g BM k k 2 2 1 M k Bm 热中子反应堆的 修正单群理论修正单群理论。 例题：例题：P100 6 6. .单群理论的修正单群理论的修正 34 反应堆的最佳形状反应堆的最佳形状 这里所谓最佳最佳，是指用同样材料，做成的 反应堆临界体积最小。 组成反应堆

15、的材料确定后，材料曲率就定组成反应堆的材料确定后，材料曲率就定 了，堆的临界几何曲率了，堆的临界几何曲率Bg也定了。但是，也定了。但是， 对同一个几何曲率对同一个几何曲率Bg值，可以有不同的形值，可以有不同的形 状和尺寸，对应不同的体积状和尺寸，对应不同的体积。 35 222 2222 gxyz B B + B + B abc g 22 B() R 22222 2.405 ()() grz BBB RH 36 反应堆的最佳形状（续） 几何形状最佳尺寸最小临界体 积 长方体abc161/ 圆柱D1.08H148/ 球R130/ 3 B 3 B 3 B 37 1. 反射层的作用以及如何选择反射层反

16、射层的作用以及如何选择反射层 2. 一侧带有反射层的反应堆一侧带有反射层的反应堆 带有反射层的球形堆带有反射层的球形堆 一侧带有反射层的柱形堆一侧带有反射层的柱形堆 3. 反射层节省反射层节省 三、有反射层反应堆的单群扩散理论三、有反射层反应堆的单群扩散理论 38 减少芯部中子的泄漏，从而减小芯部的临界体积，节省一减少芯部中子的泄漏，从而减小芯部的临界体积，节省一 部分核燃料；部分核燃料； 提高反应堆的平均输出功率，这是由于反射层的原因，其提高反应堆的平均输出功率，这是由于反射层的原因，其 芯部中子通量密度分布比裸堆的中子通量密度分布更加平芯部中子通量密度分布比裸堆的中子通量密度分布更加平 坦

17、的缘故。坦的缘故。 如何选择反射层？如何选择反射层？ 反射层材料散射截面要大反射层材料散射截面要大 反射层材料吸收截面要小反射层材料吸收截面要小 良好的慢化能力良好的慢化能力 热中子堆热中子堆常用的反射层材料有：常用的反射层材料有：H2O, D2O, 石墨，铍等。石墨，铍等。 快堆快堆常用反射层材料：常用反射层材料：重核重核 1 1. .反射层的作用以及如何选择反射层反射层的作用以及如何选择反射层 39 1 裸堆裸堆 2 有反射层反应堆有反射层反应堆 40 2 2 1 c c L k k B 0)()()( 2 rkrrD caccaccc 芯部稳态单群扩散方程芯部稳态单群扩散方程 (角标角标

18、 c) (4-52) u 该方程只有对于临界系统临界系统才成立。 u 对于任意给定材料成分及几何形状与尺寸的反应堆系统，它不一定处于稳态 引入一个特征参数引入一个特征参数k来进行调整使其达到临界。 0)()()( 2 r k k rrD caccaccc 或者写为 0)()( 22 rBr ccc Lc 为芯部的扩散长度。 (4-53) (4-54) 多区问题多区问题: 分别写出每一区的扩散方程分别写出每一区的扩散方程 证明证明: K为芯部的有效增殖因子为芯部的有效增殖因子 2 2. .一侧带有反射层的反应堆一侧带有反射层的反应堆 41 2 2 1 r r L k 0)()( 22 rkr r

19、rr 反射层稳态单群扩散方程反射层稳态单群扩散方程 （角标为（角标为 r） (4-56) 式中 (4-57) Lr为反射层的扩散长度。 边界条件边界条件为： (1) 在芯部和反射层的交界面上 rc rrcc DD (4-58) (2) 反射层的外推边界上中子通量密度为零外推边界上中子通量密度为零 42 1. 带有反射层的球形堆带有反射层的球形堆（芯部半径为（芯部半径为R，反射层厚度为，反射层厚度为T） r rk A r rk Cr rr r )cosh()sinh( )( r rTRk Cr r r )(sinh )( r rB Ar c c )sin( )(4-59) 芯部方程式解芯部方程式

20、解： 反射层方程式解：反射层方程式解： (4-60) )(tanh TRkCA r )sinh( 1)cosh()sin()cos( )sinh()sin( 22 Tk RR Tkk CD R RB R RBB AD R Tk C R RB A r rr r ccc c rc 根据：反射层外推边界反射层外推边界r=R+T处中子通量密度为零，有处中子通量密度为零，有 上式代入(4-60)可得 (4-61) 方程(4-59)、(4-61)有两个常数，由芯部与反射层交界面处边界条件芯部与反射层交界面处边界条件确定 sinh x=e(x)-e(-x)/2 cosh x=ex+e(-x)/2 tanh

21、x=sinh x / cosh x 43 1cot()1coth cccr rr RT DB RB RD LL 带反射层球形反应堆单群临界方程反射层球形反应堆单群临界方程 给出了反应堆临界时几何尺寸（给出了反应堆临界时几何尺寸（R，T）和）和 材料特性（材料特性（ Dc，Lr，Dr）之间的关系）之间的关系 当堆芯材料、反射层材料和尺寸已经确定时，用 此公式计算堆芯临界尺寸的值。 当堆芯尺寸、反射层材料和尺寸已经确定时，用 它来计算达到临界所需的堆芯材料成分 材料定，曲率知，算尺寸 尺寸知，算曲率，配材料。 44 2. 侧面带有反射层的圆柱形反应堆侧面带有反射层的圆柱形反应堆 (半径半径R，高度

22、，高度H，侧面反射层厚度，侧面反射层厚度T) 0) 2 ,() 2 ,( H r H r rc 0),(zTR r ; rrccrc DD 0),(),( 22 zrBzr ccc 芯部 反射层0),(),( 22 zrkzr rrr (4-63) (4-64) 边界条件边界条件为 (1) 在z=H/2处 (2) 在r=R+T处 (3) 在r=R处 (4-65) (4-66) (4-67) 45 RR 0 芯部周围有反射层芯部周围有反射层 部分泄露出芯部的中子在反部分泄露出芯部的中子在反 射层内被散射而返回芯部，减少了中子损失，提高射层内被散射而返回芯部，减少了中子损失，提高 了中子的不泄露概

23、率。了中子的不泄露概率。 因此在芯部材料性质相同情况下，临界体积就要因此在芯部材料性质相同情况下，临界体积就要 比裸堆的临界体积小。比裸堆的临界体积小。 反射层节省反射层节省 ：反应堆加上反射层所引起的临界尺寸：反应堆加上反射层所引起的临界尺寸 的减少。的减少。 (4-81) 2 2. .反射层节省反射层节省 46 反射层节省 47 22 22 222 405. 2 2 405. 2 effeffzr zrc HRHR BBB 对于圆柱形反应堆圆柱形反应堆，反射层节省通常分别用径向和轴向的反射层径向和轴向的反射层节省来表示 ) 22 (, 0 0 HH RR zr (4-82) 可以把有可以把

24、有反射层反应堆的几何曲率用芯部外形尺寸增大反射层反应堆的几何曲率用芯部外形尺寸增大 或或 2 来表示来表示。 带反射层球形堆反射层球形堆 2 2 R Bg 圆柱形反应堆圆柱形反应堆 R Reff eff、 、Heff称为称为等效半径等效半径、等效高度等效高度 48 )tanh( r r L T L rr rccc L T L R DRBRBDcoth1)cot(1 r rcc L T LBBtanh)tan( 带反射层球形反应堆反射层节省反射层球形反应堆反射层节省 (1) 反射层厚度很小，即反射层厚度很小，即 T Lr = Lr。 过大的增加反射层厚度是没有太大意义的！过大的增加反射层厚度是没

25、有太大意义的！ r crc L T RBLBcoth)cot( 临界方程临界方程 设Dc=Dr 设反射层节省反射层节省 ，将R=R0- =/Bc- 代入上式得 r rc c L T LB B tanharctan 1 反射层节省反射层节省 因为因为 很小，可得很小，可得 tanhx = e(x)-e(-x) / ex+e(-x) 49 讨论题讨论题 一个裸堆，加上反射层以后，几何曲率与 材料曲率还相等吗？ 50 讨论题答案 不相等了。不相等了。 裸堆如果原先是临界的，加上反射层以后，裸堆如果原先是临界的，加上反射层以后， 就超临界了。就超临界了。 如欲仍然维持临界，只能如欲仍然维持临界，只能 改变堆芯成分，降低材料曲率改变堆芯成分，降低材料曲率 减小堆芯尺寸，提高几何曲率减小堆芯尺寸，提高几何曲率 51 V H dVr V K )( 1 max 反应堆内的中子通量密度空间分布是不均匀的，而功率密度和中子通量反应堆