第3章核反应堆物理基础.pptVIP

第三章 单速中子扩散 §3.1 单速中子扩散方程 斐克定律 斐克定律的推导 * * §3.2 非增殖介质内中子扩散方程的解 §3.3 扩散长度 §3.1 单速中子扩散方程 §3.4 与能量相关的中子扩散方程 基本概念 中子角密度: 中子角通量密度: 稳态情况 中子角密度: 中子角通量密度: 一、中子流密度: 单位时间,沿着一面元ds法线方向n穿过该面元单位面积的中子数，用Jn表示。 显然,一束经过准直的中子,其运动方向与法线方向一致,中子流密度在数值上与中子角通量密度相等。即有 如果之间有一夹角 如果中子运动在4 立体角方向，那么运动在 方向， 立体角的中子对中子流密度的贡献为 4 方向的贡献为 可看着 分别为两半球立体角积分 中子扩散:  中子扩散服从斐克扩散定律 中子经过逐次散射碰撞从高中子密度区向低中子密度区迁移的过程，称为中子扩散。  （1）介质无限、均匀 斐克定律： 描述中子流密度与中子通量密度之间的关系 斐克定律的推导 条件假设： （2）在实验室坐标系中散射各向同性 （3）介质的吸收截面很小，即 （4）中子通量密度是随空间位置缓慢变化的稳态函数 坐标系如图所示，在xy平面上有一面积dA，在r 处有一体积元dV。dV内每秒发生散射的中子数为 这些散射中子向dA面积运动的数目为 而 因此，自dV向着dA运动的中子数是 因此，每秒自dV内散射出来未经碰撞而到达dA上的中子数是 这些中子 不是全部都能达到dA，有些在未到达之前就被介质吸收或从这个方向散射出去，实际能到达dA的几率为 因为 每秒自xy平面上方整个半空间区域散射穿过dA面积上的中子数为 在球坐标系中 其中， 在球坐标系中 为沿负z方向的中子流密度 表示每秒自上而下（负z方向）穿过dA单位面积上的中子数，穿过dA面积上的中子数为 积分号下的 是一个待求未知函数， 由于它是随空间缓慢变化的函数，可以在原点按泰列级数展开 对 积分从0到2π，故含x,y项的积分将等于零 对 积分从0到2π，故含x,y项的积分将等于零，消去等式两边的dA，得到 经计算，可以得出 用同样的方法，可以得到沿正z方向的中子流密度 这样，单位时间沿z方向的净中子数为 叫做沿z方向的中子流密度或净中子流密度。 同样,可以求出x方向和y方向上的中子流密度分别为 如果所讨论的面积元并不垂直于任一坐标轴的方向，它的法线n与x、y、z 轴分别成α、β、γ角度，那么单位时间穿过这一面积元单位面积的净中子数J便等于 （3-9） （3-14） 可以把（3-9）式写成所讨论面积元的单位法线矢量n与矢量J的乘积，即 （3-14）式称为斐克定律。它表示；中子流密度J正比于负的通量密度梯度，其比例常数D叫做扩散系数. 令 二、单速扩散方程的建立 中子密度对时间的变化率应等于中子的产生率减去中子的吸收率和泄漏率,即 产生率（S）-泄漏率（L）-吸收率（A） （3-14） 首先计算泄漏率 计算中子泄漏示意图 点 处有一个小体积元 上下两表面面积 单位时间内由下表面进入 dV的中子数 由上表面流出 dV的中子数 通过平行于xy平面的上下两个平面，单位时间从dv中泄漏出去的中子数为 用同样方法可以求出通过平行于yz平面和xz平面两个表面，从dv中泄漏出去的中了数分别是 中子从dV内泄漏的总数应等于以上三项之和，这样，单位时间，单位体积泄漏出去的中子数 （3-15） 若护散系数D与空间位置无关，那么便得到 每秒每单位体积内被吸收的中子数为 设每秒每单位体积内产生的中子数为 （3-16） 稳态单速中子扩散方程 泄漏率 吸收率 得单速中子扩散方程 产生率 根据单位时间单位体积内的中子数守恒 产生率（S）-泄漏率（L）-吸收率（A） 是拉普拉斯算符，在反应堆计算常用的几种坐标系中 直角坐标系 柱坐标系 球坐标系 三、扩散方程的边界条件 由扩散方程可以得到普遍解，只有加上问题具体的物理或几何 特性才能给出唯一解，边界条件就是问题具体特性的主要内容 经常用到的几种边界条件： 垂直于界面的中子流密度相等 （1）在扩散方程适用的范围内，中子通量密度必须是正的和有限的实数 （2）在两种不同介质的界面上，中子通量密度相等 （3）与真空交界的介质界面上，由真空返回的中子流等于零 外推边界距离d §3.2 非增殖介质内中子扩散方程的解 稳态情况下的扩散方程为 如果s=0，即对于除中子源所在位置以外的无源区域，扩散方程具有如下的齐次形式 其中 和 L2具有长度平方的量纲，通常我们称L为中子的扩散长度，它是表征中子在介质中扩散特性的一个重要的量 方程 称为波动方程或亥姆霍兹方程，可以用数理方程一些标准方法求得它的普遍