第3章-核反应堆临界理论汇编.ppt

第3章 核反应堆临界理论 研究均匀裸堆的临界问题。主要研究下面两个问题： （1）各种形状的反应堆达到临界状态的条件（临界条件），临界时系统的体积大小和燃料成分及其装载量； （2）临界状态下系统内中子通量密度（或功率密度）的分布。 研究方法 首先，研究均匀裸堆。 其次，研究带反射层的反应堆。 最后，研究非均匀堆。 方法是先单群，后多群研究。 什么是均匀裸堆？ 均匀：燃料、慢化剂、结构材料等堆芯内一切材料均匀混合： 裸堆：没有反射层； 中子源：有增殖介质。 在非均匀堆的研究上，从理论上给出了均匀堆和非均匀堆的中子通量密度分布的差别，并对非均匀堆与均匀堆的四因子公式差别作了简要分析。 3.1 均匀裸堆的单群理论 3.1.1 单群扩散方程 中子源S主要是增殖介质的裂变源 称为曲率，它等于 ，后者表征了中子通量密度空间分布的弯曲程度 裸堆单群临界方程 裸堆单群临界方程 材料曲率 3.1.2 平板裸堆 图3-1 无限宽有限厚的平板均匀裸堆 3.1.2 平板裸堆 均匀平板的波动方程为 边界条件 1外边界 2 中心处 解 由于B2＞0，所以方程 的通解可写为 (x)＝AcosBx+CsinBx (3-8) 式中，A和C为待定系数。由边界条件(2)，则有C＝0。因而式(3—8)可写成下列形式 (x)＝AcosBx 由边界条件(1) 得 因为A≠0，故有 n=1,3,5 因而 这些常数Bn，称为该方程的本征值，对应的函数 n(x)＝Acos(Bnx)称为本征值Bn对应的本征函数。 解 对于临界的反应堆，随着时间变大，除去第一个模态(n=1)外的所有模态(n1)都衰减了。渐进的（或持久的）中子通量密度取决于n=1的模态。于是， 几何曲率 临界系统中的中子通量密度分布为 系数A 系数A由功率条件决定。一个系统只要是临界的，则中子通量密度分布中的A可为任意常数。这个事实表明，从理论上讲，堆功率是可任意提高的。实际上由于热工传热条件，燃料制造工艺等工程问题的限制，堆功率总是限止在某一数值以下。因而系数A也由功率条件限止为某一数值。 无限长平板堆单位面积所对应的体积所发出的功率为 积分后得 3.1.3 有限高圆柱形均匀裸堆 设一有限高圆柱形 均匀裸堆，高为H， 半径为R，如图3-2 所示，采用圆柱形 坐标，坐标原点位 于轴线的半高度上。 在圆柱形坐标系中 拉普拉氏算符的表达式为 中子通量密度分布是对称的，故与 无关。因而有限高圆柱均匀裸堆的波动方程可以写为 (1)不计外推长度时，反应堆外边界上，中子通量密度为零 (2) 中子通量密度分布对称 二阶偏微分方程 采用分离变量法 代入方程(3-16)得到 上式等号左边只与r有关，等号右边只与z有关。因为r和z是两个彼此无关的独立变量，要使式(3-22)保持相等关系，只有两边都是常数才行。该常数记为 ,则得到 方程（3-23） 方程 解 由边界条件(2)，式(3-20)得A2＝0，则 方程（3-24） 方程 及 令 化为零阶贝塞尔方程 通解 C1，C2是两个待定常数。J0(x)，Y0(x)分别是零阶第一类及第二类贝塞尔函数，它们随x的变化见图3-3。 图3-3 零阶贝塞尔函数 零阶贝塞尔方程解 当x＝0时，J0(0)＝1，而Y0(0)→-∞。 r＝0处堆内中子通量密度有限，所以C2＝0。则式(3-32)可写成 由边界条件式(3-17) 正如图3-3所示，在x＝x1，x2…等处，函数J0(x)为零，即J0(xn)＝0。故 必须取如下的本征值 曲率 同样，当反应堆达到临界时，只有n＝1这个本征值有物理意义。 因为 所以 临界时，r向的中子通量密度分布为 反应堆的中子通量密度分布为 常数A 取环状体积元dv＝2 drdz，在(r,z)处该体积元所发出来的功率为 从零到R对r积分，从-H/2到H/2对z积分可得到功率 代入 的表达式，上式可化为 常数A（续） 分成r积分和z积分，得 3.1.4 单群临界方程与临界条件 当反应堆临界时，材料曲率应该与几何曲率相等，即 材料曲率与几何曲率的关系与堆状态有关 ： 次临界 ＞ keff＜1 临界 ＝ keff＝1 超临界 ＜ keff＞1 临界方程