chapt4.1-均匀裸堆临界计算综述.pptVIP

* 第四章：反应堆临界理论 4.1均匀裸堆的单群理论 4.1.1均匀裸堆的单群扩散方程 4.1.2热中子反应堆的临界条件 4.1.3各种几何形状裸堆计算 4.1.4单群理论修正-徙动长度 4.2有反射层均匀堆 4.3非均匀反应堆 * 4.1.1均匀裸堆单群扩散方程 均匀——燃料与慢化剂等堆内一切材料均匀混合 裸堆——没有反射层的反应堆。 一群——中子能量相同，（~单速中子） * 4.1.1均匀裸堆单群扩散方程 单群中子扩散方程： 假设所有中子都是裂变瞬发中子，热中子反应堆内单位时间、单位体积内产生的中子数 其中 L2=D/?a 扩散长度 即得堆芯单群中子扩散方程 边界条件： 1）在反应堆外推边界处，中子通量密度为0； 2）在不同介质交界面上，中子通量密度相等，中子流密度相等 * 单群中子扩散方程的解 采用分离变量法： 单群中子扩散方程变为 上式左边只与r有关，而右端只与t有关，因此等式两端必须都等于某一常数（设为-B2），则有 为典型的波动方程，又称亥姆霍兹方程，B2称为方程的本征值（或特征值） 对于给定几何形状系统，在满足边界条件下，可得到一系列从小到大排列的本征值Bn2，及对应的本征函数?n (n=1,2,3…)。 n=1的本征函数称为基波特征函数或一阶谐波，其余为高阶谐波； n=1的本征值B12称为反应堆的几何曲率Bg2(它只与堆的大小和形状有关) * 由于特征函数?n(r)的正交性，对于每一个n值的项都是线性独立的，因而对于每一个Bn2值和?n(r) ，有一个Tn(t)与之对应： 两端同乘L2/(1+L2 Bn2) 与时间及空间相关的中子通量密度： * 例：无限宽有限厚平板均匀裸堆 x a （2）在平板的外推边界面（x=d）上通量为0。 边界条件： （1）中子通量对称分布，在x=0处净中子流密度为0。根据斐克定律有： 外推长度d=2/(3?s)=2D~厘米量级，而一般热中子堆的堆芯尺寸在米的量级，故可认为外推边界就是堆芯几何边界 通解： ?(x)=A’cosBx+CsinBx ?(x)=A’cosBx Bn= (2n-1)?/a n=1,2,3… 无限平板堆中子通量密度： * 4.1.2热中子反应堆的临界条件 通过前一小节的讨论可以知道： 1）特征值Bn2随n增加而增大。最小特征值是n=1时的B12值； 2）当n增加时，kn单调递减，也就是说对应于最小特征值B12之k1是k1,…, kn中的最大值。另外，考虑到Bn2与系统尺才有关，当系统尺寸加大时， Bn2便减小，因而改变系统的尺寸就可以改变Bn2值，从而也就改变了kn值。 Bn= (2n-1)?/a n=1,2,3… 平板堆： * 第二种情况： 着k1 ＞1，则(k1-1)＞0，这时中子通量密度将随时间不断地增长，反应堆将处于超临界状态。 第三种情况：若通过调整反应堆堆芯尺寸，使k1恰好等于1，则其余k2 ，……，kn(n1)的值都将小于1，那么?(r,t)中n=1的项将与时间无关，而n1的各项将随时间而衰减。因而当时间足够长时，n1各项都已衰减到零，系统达到稳态，这时中子通量密度按基波形式分布，反应堆处于临界状态。因而反应堆维持临界的条件是k1=1，显然k1便是前面定义的有效增殖系数。 第一种情况： 对于一定几何形状和体积的反应堆芯部，若B12对应的k1小于1，那么，其余的k2，……， kn都将小于1，这时的(kn-1)都是负值。这时，?(r,t)将随时间t按指数规律衰减，因而系统处于次临界状态。 中子通量密度分布： * 结论： 裸堆单群近似的“临界条件” 又称为单群理论的临界方程： 最小特征值，即几何曲率 当反应堆处于临界状态时，中子通量密度按最小特征值Bg2所对应的基波特征函数分布，也就是说稳态反应堆的中子通量密度空间分布满足波动方程： 这里B12系波动方程的最小特征值Bg2(几何曲率) 热中子泄漏率 热中子不泄漏几率： * 例题 x a 设有如图所示一维石墨慢化反应堆，k?=1.06，L2=0.03米，?tr=0.028米。试求：(1)达到临界时反应堆的厚度H和中子通量密度的分布；(2)设取H=2.5米，试求反应堆的有效增殖系数k。 (1) 反应堆临界时 对应的几何曲率 外推距离 d=0.7104 ?tr=0.7104×0.028?0.02米 临界时反应堆的厚度 H=a-2d=2.222-0.04=2.182米 中子通量密度的分布 （2）若H=2.5米，则反应堆几何曲率 反应堆的不泄漏几率 有效增殖系数 * 4.1.3几种几何形状裸堆计算 1）球形反应堆：半径R(含外推距离) 应用球坐标系统，并把原点取在球心上。根据对称性条件，这时的波动方程为 通解： 边界条件： 1）r?0，??有限值：