中子扩散理论精选PPT.ppt

关于中子扩散理论第1页，讲稿共62张，创作于星期一2反应堆物理的反应堆物理的核心问题核心问题之一：确定堆内中子通量密度按空间和能量的之一：确定堆内中子通量密度按空间和能量的分布分布 第二章通过求解中子慢化方程，解决了中子通量密度按能量的分布，第二章通过求解中子慢化方程，解决了中子通量密度按能量的分布，(E)(E)E E，即中子能谱，即中子能谱 本章，将研究中子通量密度按空间的分布，即本章，将研究中子通量密度按空间的分布，即(r)(r)r r第2页，讲稿共62张，创作于星期一3Contents引言引言（输运过程、输运理论及扩散现象）（输运过程、输运理论及扩散现象）单能中子扩散方程单能中子扩散方程非增殖介质内中子扩散方程的解非增殖介质内中子扩散方程的解扩散长度、慢化长度和徙动长度扩散长度、慢化长度和徙动长度第3页，讲稿共62张，创作于星期一4输运过程及输运理论输运过程及输运理论中子状态的描述中子状态的描述反应堆物理与屏蔽计算的基本方法反应堆物理与屏蔽计算的基本方法 一、引言一、引言第4页，讲稿共62张，创作于星期一51、输运过程（输运过程（TransportTransport）以及输运理论）以及输运理论对于自然界的微观粒子(中子、光子、电子、离子和分子等)，在介质中会发生无规则碰撞，使得单个粒子的运动形式是杂乱无章的，即某一时刻、在介质中某一位置、具有某种能量与某一运动方向 的粒子，在稍晚些时候，将运动到介质中的另一位置、以另一能量和另一运动方向出现，这一现象称之为粒子输运粒子输运过程过程。对于大量粒子而言，其运动形式呈现出一定的统计规律，使得人类对粒子行为的研究成为可能。那么，用于描述粒子在介质中迁移的统计规律的数学理论，就统称为粒子输粒子输运理论运理论。第5页，讲稿共62张，创作于星期一6u发展简史：发展简史：最早的粒子输运理论是以Clausius(1857)、Maxwell(1860)、Boltzmann(1868)提出的“分子运动论”作为基础发展起来的；1872年，Boltzmann建立了著名的Boltzmann方程方程（又称输运方程），当时是用来描述气体从非平衡态到平衡态的过渡过程；1910年Hilbert论述了Boltzmann方程解的存在性与唯一性，奠定了粒子输运理论的数学基础；1939年发现中子后，随着核反应堆和核武器的出现，中子输运理论中子输运理论得到极快发展；1943年Wick、Marshak、Mark等人提出并发展了球谐函数方法，使得高精度地解析求解Boltzmann中子输运方程成为可能；1946年Von Neumann 和 Ulam等开发了第一个用概率论方法（Monte Carlo方法）计算中子链式反应的程序；1955年Carlson等人提出了离散纵标法(即早期SN方法)；在上述方法的基础上，产生了大批应用程序软件。第6页，讲稿共62张，创作于星期一7中子状态中子状态:位置矢量位置矢量 r(x,y,z)、能量、能量E(或运动速或运动速度度v)、运动方向、运动方向(,)、时间时间t ：单位矢量，模等于：单位矢量，模等于1，方向表示中子的运动，方向表示中子的运动方向，通过极角方向，通过极角 和方和方位角位角 来表示来表示2、中子状态的描述、中子状态的描述第7页，讲稿共62张，创作于星期一8中子角密度中子角密度：在：在r处单位体积内和能量为处单位体积内和能量为E的单位能量间隔内，运动的单位能量间隔内，运动方向为方向为 的单位立体角内的中子数目。的单位立体角内的中子数目。中子角通量密度中子角通量密度：沿沿 方向在单位时间内穿过垂直于这个方向方向在单位时间内穿过垂直于这个方向的单位面积上的中子数目。的单位面积上的中子数目。对中子角密度和中子角通量对所有立体角方向积分，可得前面所定义对中子角密度和中子角通量对所有立体角方向积分，可得前面所定义的中子密度和中子通量密度的中子密度和中子通量密度第8页，讲稿共62张，创作于星期一9中子输运理论的基本问题之一，就是采用中子角密度（或中子角通量密度）来描述中子输运过程。为了得到中子角通量密度，需要建立描述中子输运过程的精确方程，即“玻尔兹曼方程玻尔兹曼方程”。它是一个含有位置r(x,y,z)、能量E(或运动速度v)、运动方向(,)、时间t七个自变量的偏微分积分方程，求解过程非常复杂，只有在极个别的简单情况下，才能求出解析解。第9页，讲稿共62张，创作于星期一确定性方法确定性方法(Deterministic method)数学模型用数学物理方程表示，然后采用数值方法求解数学模型用数学物理方程表示，然后采用数值方法求解优点：计算快速优点：计算快速缺点：模型简化不够精确，大型多维问题需大量计算时间及存储空间等缺点：模型简化不够精确，大型多维问题需大量计算时间及存储空间等典型方法：离散纵标法（典型方法：离散纵标法（SN）非确定性方法非确定性方法(蒙特卡罗方法，蒙特卡罗方法，Monte Carlo method)：基于统计理论，通过计算机的随机模拟来跟踪中子在介质中的运动基于统计理论，通过计算机的随机模拟来跟踪中子在介质中的运动优点：计算精确，可以模拟三维复杂几何模型优点：计算精确，可以模拟三维复杂几何模型缺点：计算耗时，特别是对于深穿透问题（缺点：计算耗时，特别是对于深穿透问题（Deep-penetration）混合方法混合方法研究热点研究热点3、反应堆物理与屏蔽计算基本方法反应堆物理与屏蔽计算基本方法10第10页，讲稿共62张，创作于星期一11中子扩散理论中子扩散理论求出介质内求出介质内中子角通量密度中子角通量密度的分布的分布,才算对介质内中子的分布有才算对介质内中子的分布有了全面了解了全面了解.要做到这一点，需要研究中子输运理论，求解要做到这一点，需要研究中子输运理论，求解中子输运方程中子输运方程。这是。这是一个非常复杂和困难的任务一个非常复杂和困难的任务.在本课程中，我们研究输运理论的简在本课程中，我们研究输运理论的简化形式化形式中子扩散理论中子扩散理论。玻尔兹曼输玻尔兹曼输运方程运方程中子扩散中子扩散方程方程单群中子单群中子扩散方程扩散方程假设中子通量密度角分假设中子通量密度角分布各向同性布各向同性假设中子具有单假设中子具有单一能量一能量本节内容本节内容第11页，讲稿共62张，创作于星期一121.菲克定律菲克定律2.菲克定律的推导菲克定律的推导3.菲克定律和扩散方程的使用范围菲克定律和扩散方程的使用范围4.单能中子扩散方程的建立单能中子扩散方程的建立5.扩散方程的边界条件扩散方程的边界条件 二、二、单能中子扩散方程单能中子扩散方程第12页，讲稿共62张，创作于星期一13分子扩散现象分子扩散现象香水分子的扩散（无风状态）墨滴在静水中的扩散血液中的养分透过细胞膜向细胞内扩散分子扩散是由于分子间的无规则碰撞产生的，使得分子从密度大的地方向密度小的地方扩散，并且分子扩散的速率与分子密度的梯度成正比，也就是服从“菲克定律”。同样，中子的扩散现象也服从“菲克定律”。只不过由于中子密度（1016m-3）比介质原子核密度（1028m-3）要小得多，因而中子的扩散主要是中子与介质原子核碰撞的结果，中子之间的碰撞可以忽略。1、菲克定律、菲克定律(Ficks law)第13页，讲稿共62张，创作于星期一14中子从通量高的地方流向通量低的地方，中子从通量高的地方流向通量低的地方，通量差别越大，中子通量差别越大，中子“流量流量”越大越大菲克定律菲克定律：上式中的上式中的 被称为被称为中子流密度中子流密度（简称中子流、（简称中子流、或流。或流。Current）.中子流密度是一个向量中子流密度是一个向量,其方向是通量场的负梯度方向其方向是通量场的负梯度方向.其数值等于其数值等于垂直于梯度方向垂直于梯度方向的单位面积上每秒穿过的的单位面积上每秒穿过的净中子数目。净中子数目。单位：中子单位：中子/cm/cm2 2.S.S第14页，讲稿共62张，创作于星期一15第15页，讲稿共62张，创作于星期一16l中子流密度是向量，可以写成三个分量之和中子流密度是向量，可以写成三个分量之和l其中三个分量分别称为该方向的分中子流密度，每个其中三个分量分别称为该方向的分中子流密度，每个分量可写成两个分量只差分量可写成两个分量只差JZ+是沿是沿z轴正方向每秒穿过轴正方向每秒穿过x-y平面上单位面积的中子数平面上单位面积的中子数JZ-是沿是沿z轴负方向每秒穿过轴负方向每秒穿过x-y平面上单位面积的中子数平面上单位面积的中子数第16页，讲稿共62张，创作于星期一17l如果某平面与中子流密度如果某平面与中子流密度 方向不垂直，方向不垂直，那么每秒通过该平面上单位面积的净中子那么每秒通过该平面上单位面积的净中子数是数是第17页，讲稿共62张，创作于星期一18第18页，讲稿共62张，创作于星期一19中子流密度与中子通量密度的差别中子流密度与中子通量密度的差别:中子流密度用于描述中子的定向运动，是矢量中子流密度用于描述中子的定向运动，是矢量中子通量密度用于计算中子通量密度用于计算核反应率核反应率，是标量，是标量 两者的量纲相同两者的量纲相同 当所有中子运动方向相同时，中子通量与中子当所有中子运动方向相同时，中子通量与中子流数量（大小）相等。流数量（大小）相等。第19页，讲稿共62张，创作于星期一20场论知识场论知识数量场数量场的的梯度梯度向量场向量场 的的散度散度 哈密顿算符哈密顿算符 拉普拉斯算符拉普拉斯算符 第20页，讲稿共62张，创作于星期一21考虑稳态情况，同时假设：考虑稳态情况，同时假设：(1)介质是介质是无限无限的、的、均匀均匀的；的；(2)在在实验室坐标系实验室坐标系中中散射是各向同性散射是各向同性的；的；(3)介质的吸收截面很小，即介质的吸收截面很小，即 a s(弱吸收介质弱吸收介质)；(4)是随空间位置是随空间位置缓慢变化缓慢变化的函数。的函数。2、菲克定律的推导、菲克定律的推导第21页，讲稿共62张，创作于星期一22以所研究的以所研究的点作为坐标点作为坐标原点原点第22页，讲稿共62张，创作于星期一23考虑上半空间发生的散射使多少中子从考虑上半空间发生的散射使多少中子从上到下穿过上到下穿过 dA首先考虑体积元首先考虑体积元 dV 中的散射中子有多少中的散射中子有多少由于散射中子由于散射中子各向同性各向同性地飞向四面八方地飞向四面八方,飞向飞向 dA的只占一部分的只占一部分.这一份额等于这一份额等于dA的面积与以的面积与以 r 为半径的球面积之比为半径的球面积之比,再乘以再乘以cos 此外并非所有飞向此外并非所有飞向dA的中子都能够到达的的中子都能够到达的dA,沿途的碰撞会使得部分中子沿途的碰撞会使得部分中子“偏离航向偏离航向”，此，此外还有一部分中子会被介质吸收掉外还有一部分中子会被介质吸收掉第23页，讲稿共62张，创作于星期一24第24页，讲稿共62张，创作于星期一25 把上半空间所有地方的散射中子的贡献统把上半空间所有地方的散射中子的贡献统统考虑进来统考虑进来,即对上半空间积分即对上半空间积分,就得到从上就得到从上而下穿过而下穿过dA的总中子数目。的总中子数目。这个数目就是这个数目就是沿负沿负z方向方向的的分中子流密度分中子流密度 乘以乘以dA第25页，讲稿共62张，创作于星期一26由于已经假设中子通量密度是随空间位置缓慢变化的中子通量密度是随空间位置缓慢变化的，将(r)在原点处按泰勒级数展开，取1阶项，代入积分可得第26页，讲稿共62张，创作于星期一27 第27页，讲稿共62张，创作于星期一28推导过程推导过程由于cos和sin从0到2的积分为0，所以含有x，y项的积分将等于0，则可以得到第28页，讲稿共62张，创作于星期一29分步积分洛必达法则第29页，讲稿共62张，创作于星期一30第30页，讲稿共62张，创作于星期一31 令令 ，便完成了，便完成了菲克定律菲克定律之推导，得到之推导，得到斐克定律斐克定律：中子流密度 正比于负的中子通量密度的梯度，其比例常数为扩散系数D第31页，讲稿共62张，创作于星期一32u介质是无限的、均匀的、散射各向同性；介质是无限的、均匀的、散射各向同性；有限介质内，在距离表面几个自由程之外的内部区域，斐克有限介质内，在距离表面几个自由程之外的内部区域，斐克定律是近似成立的；定律是近似成立的；在距真空边界两三个自由程以内的区域，不适用。在距真空边界两三个自由程以内的区域，不适用。u介质的吸收截面很小，即介质的吸收截面很小，即 a 0时，(3-52)解为(3-53)由边界条件(1)可得(3-54)因此根据边界条件(2)可以求出(3-55)中子通量密度为中子通量密度为(3-56)第52页，讲稿共62张，创作于星期一53(3-57)考虑系统对称性系统对称性，用|x|代替上式中的x，可得对所有x均适用的表达式对于无限厚度平面源无限厚度平面源，a，有(3-61)第53页，讲稿共62张，创作于星期一54当a/L=3（介质厚度为中子扩散长度3倍时）时，除在边界附近外，中子通量密度的分布与无限介质内的分布相差不多。薄板泄露较大，边界处中子通量密度下降很快；厚板（大于三个扩散长度），大部分中子在到达边界以前被散射回来，泄露很小=反应堆没有必要采用过厚的反射层反应堆没有必要采用过厚的反射层第54页，讲稿共62张，创作于星期一553.包含两种不同介质的情况包含两种不同介质的情况P76第55页，讲稿共62张，创作于星期一561.1.扩散长度扩散长度(diffusion length)(diffusion length)(3-75)大多数元素的散射截面与能量无关，当热中子能谱按麦克斯韦谱分布时，大多数元素的散射截面与能量无关，当热中子能谱按麦克斯韦谱分布时，热中子吸收截面等于热中子吸收截面等于(3-76)将上式代入将上式代入(3-75)式可得式可得(3-77)四、扩散长度、慢化长度和徙动长度四、扩散长度、慢化长度和徙动长度第56页，讲稿共62张，创作于星期一57考虑无限介质内有一热中子点源的情况，在 的球壳内，第57页，讲稿共62张，创作于星期一58扩散长度的物理意义：扩散长度的物理意义：热中子扩散长度的平方热中子扩散长度的平方,等于无限大介质中的点等于无限大介质中的点源放出的热中子从产生地点到被吸收地点的直源放出的热中子从产生地点到被吸收地点的直线距离的均方值的六分之一线距离的均方值的六分之一.第58页，讲稿共62张，创作于星期一59几种慢化剂在几种慢化剂在293K时热中子扩散参数时热中子扩散参数室温下热中子在石墨中室温下热中子在石墨中(s=0.385cm-1),求解：求解：(1)散射平均自由程（散射平均自由程（2）吸收前的碰撞次数（）吸收前的碰撞次数（3）产生到被吸收的直线距离（）产生到被吸收的直线距离（4）中）中子吸收前走过的路程（吸收平均自由程）子吸收前走过的路程（吸收平均自由程）第59页，讲稿共62张，创作于星期一602.慢化长度慢化长度 源中子能量为E0，热中子能量为Eth。将E0到Eth的中子称为快群中子快群中子。将Eth以下的中子称为热群中子热群中子，同时定义一个移出（减速）截面移出（减速）截面 1使快群转移到热群的中子转移率快群转移到热群的中子转移率源中子慢化到热中子的平均碰撞次数源中子慢化到热中子的平均碰撞次数:因而移出（减速）截面移出（减速）截面 为(3-82)第60页，讲稿共62张，创作于星期一61H2OD2OCBe轻水堆沸水堆高温气冷堆th/(10-4m2)27.5123352904050300写出无限介质点源情况下快群中子的扩散方程无限介质点源情况下快群中子的扩散方程或(3-83)(3-84)将L1定义为慢化长度慢化长度，将L12 称为热中子年龄热中子年龄，用 th表示。引进中子年龄中子年龄定义(3-85)(3-86)随中子能量的降低而增大，当能量为随中子能量的降低而增大，当能量为Eth时，即为热中子年龄。时，即为热中子年龄。并不具有时间的量纲，而具并不具有时间的量纲，而具有长度平方的量纲有长度平方的量纲。与勒类似，是表征能量的另一种变量形式。可以证明，慢化长度平方 L12 或 热中子年龄 th 和 扩散长度的平方 L2 具有类似的物理意义：慢慢化长度的平方或热中子年龄等于无限介质内中子自源点产生出发到在介质内慢化到年龄化长度的平方或热中子年龄等于无限介质内中子自源点产生出发到在介质内慢化到年龄 th(Eth)时所穿行直时所穿行直线距离的均方值的六分之一线距离的均方值的六分之一。第61页，讲稿共62张，创作于星期一感谢大家观看第62页，讲稿共62张，创作于星期一