中子的输运

合集下载

一维球几何中子输运方程计算方法研究的开题报告

一维球几何中子输运方程计算方法研究的开题报告

一维球几何中子输运方程计算方法研究的开题报告
题目:一维球几何中子输运方程计算方法研究
一、研究背景和意义
中子输运是核能反应堆中的重要现象之一,其中中子的漂移和散射
是最基本的现象。

而针对中子输运的计算方法是研究这些现象中最为关
键的一个方面,对于准确模拟中子输运行为和核反应堆的设计分析具有
重要的意义。

目前,一维球几何中子输运方程是求解中子输运运动学的
基本方程,在核反应堆的瞬态和稳态计算中被广泛应用。

因此,对一维
球几何中子输运方程计算方法进行深入的研究是十分必要的。

二、研究内容
本研究将重点研究一维球几何中子输运方程的计算方法,主要包括
以下内容:
1. 一维球几何中子输运方程简介及基本假设。

2. 一维球几何中子输运方程的数学推导及求解方法。

3. 常用的一维球几何反应堆中子输运计算程序及其特点。

4. 一维球几何中子输运方程数值计算方法的优化及改进。

三、研究方法
本研究将采用文献调研和实验方法共同开展,首先通过文献查阅及
分析来系统学习已有的理论及计算方法;接着通过搭建计算系统进行实
验模拟,从数值计算的角度分析目前计算方法中存在的不足和不稳定性,并提出针对性的优化和改进方法。

四、研究成果预期
本研究旨在探索一维球几何中子输运方程计算方法的新思路和新方法,并通过优化和改进原有方法,提高计算效率和计算准确度,为模拟中子输运行为和核反应堆设计等领域提供更加准确和可靠的计算手段。

中子输运方程源迭代加速方法研究

中子输运方程源迭代加速方法研究

学校代码:10200研究生学号:2017101955分类号:O24密级:无中子输运方程源迭代加速方法研究Study on source iteration acceleration for neutron transportequation作者:李悦指导教师:安恒斌研究员一级学科:数学二级学科:计算数学研究方向:并行计算学位类型:学术硕士东北师范大学学位评定委员会2020年5月东北师范大学硕士、博士研究生学位论文检测终稿确认单根据《东北师范大学研究生学位论文检测办法》规定,我校硕士、博士研究生通过论文答辩后,须对其学位论文进行学术不端行为检测,检测结果将作为是否授予学位的重要依据。

提交中国知网硕士、博士论文数据库的学位论文也将以此版本为准。

申请学位的研究生及指导教师郑重承诺:提交检测的学位论文为最终版本,学位论文一经检测,不以任何理由更换检测论文。

学生签名:指导教师签名:年月日年月日独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。

据我所知,除了文中特别加以标注和致谢的地方外,论文中不包含其他人已经发表或撰写过的研究成果,也不包含为获得东北师范大学或其他教育机构的学位或证书而使用过的材料。

与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。

学位论文作者签名:日期:学位论文版权使用授权书本学位论文作者完全了解东北师范大学有关保留、使用学位论文的规定,即:东北师范大学有权保留并向国家有关部门或机构送交学位论文的复印件和磁盘,允许论文被查阅和借阅。

本人授权东北师范大学可以将学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或其它复制手段保存、汇编学位论文。

(保密的学位论文在解密后适用本授权书)学位论文作者签名:指导教师签名:日期:日期:摘要在核反应堆数值模拟等各类应用中,需要求解中子输运方程,且中子输运方程在其求解过程中占据着主要的计算量.中子输运方程描述中子群体在介质内的输运过程.该方程是关于空间,时间,能量和运动方向的偏微分—积分方程.由于方程的复杂性,无法得到与之对应的解析解,因此一般情况下只能对中子输运方程进行数值求解.中子输运方程离散后是一个大型非线性(线性)方程组,其数值求解过程十分复杂,且所需计算量很大,求解时间较长.对于离散后的中子输运方程,需要采用多个层次的迭代方法进行求解.其中源迭代方法是求解中子输运方程的关键算法.由于源迭代过程占据整个求解过程的主导计算量,且源迭代方法在某些情况下收敛非常慢,甚至不收敛,所以针对源迭代的加速研究尤为重要.目前为止,已有很多关于源迭代加速方法的研究.其中非线性扩散综合加速方法(NDA)可以用来加速各向异性散射情况下的中子输运方程的求解.但随着散射的各向异性增强,NDA方法的加速效果显著下降.Anderson加速方法(AA)近几年被广泛地应用于提高不动点迭代收敛速度.由于源迭代过程可以看作一类不动点,所以可以自然地应用AA方法提高源迭代的收敛速度.本文比较了NDA方法和AA方法在加速源迭代收敛速度方面的效果.分析了两类方法的适用性.进一步考虑了两类加速方法的结合,获得了AA-NDA方法以及AA-NDA(2)方法.AA方法与NDA方法的结合改善了NDA方法随着各向异性增强加速效果下降的不足.通过一维情形下的数值算例,对比了AA方法, NDA方法以及AA-NDA,AA-NDA(2)方法对于求解定源中子输运方程的不同加速效果,验证了AA-NDA方法及AA-NDA(2)方法的有效性.另外,在三维非结构网格情形下,比较了AA加速源迭代方法和无AA加速源迭代方法求解中子输运方程的不同数值结果,以及不同核处理器情形下各个方法并行计算的数值结果.关键词:NDA加速;Anderson加速;中子输运方程;非结构网格;源迭代方法AbstractIn various applications such as nuclear reactor numerical simulation,it is neces-sary to solve the neutron transport equation.The solution of the neutron transport equation occupies the main calculation in many important applications of numeri-cal simulation.The neutron transport equations are used to describe the transport process of the neutron population in the medium.This equation is a kind of par-tial differential integral equations of multidimensional independent variables about space,time,energy and direction.Because of the complexity of the equation,it can only be solved numerically in general.After the neutron transport equation is dis-cretized,a large scale nonlinear equations are obtained.The source iteration is the most popular method for solving the discretized neutron transport u-ally,the source iteration converges slowly.By now there has been much research on the acceleration for source iteration.One of widely used acceleration method is the nonlinear diffusion acceleration method(NDA).NDA method can be used to accelerate the solution of fixed-source problem,a typical anisotropic scattering model.However,the effect of NDA method will deteriorate when the scattering source is less isotropic.Anderson acceleration method(AA)has been widely used to improve the convergence rate of fixed point iteration in recent years.Because the source iteration process can be regarded as a kind of fixed point iteration,AA method can be used naturally to improve the convergence rate of source iteration. In this paper,the comparison study is conducted about the effectiveness of NDA method and AA method,and analysis about the applicability of the methods is given.Furthermore,by the combination of the methods,two kinds of new methods, AA-NDA,AA-NDA(2)is obtained.The combination method of AA and NDA can be used to remedy the degraded performance of NDA when scattering is anisotrop-ic.Through numerical examples in one-dimensional case,the different acceleration effect of AA method,NDA method,AA-NDA,AA-NDA(2)method are compared for solving the neutron transport equation of fixed source.In addition,in the case of three-dimensional unstructured grid parallel computing,The results of AA accel-erated source iteration method and non-AA accelerated source iteration method are compared.Key words:NDA acceleration;Anderson acceleration;neutron transport e-quation;unstructured grid;source iteration method目录中文摘要 (I)英文摘要 (II)目录 (III)第一章绪论 (1)1.1选题意义 (1)1.2中子输运方程 (2)1.3论文的主要工作 (4)1.4论文的组织结构 (5)第二章中子输运方程及数值计算方法 (7)2.1散射源的球谐函数展开 (7)2.1.1球谐函数与Legendre多项式 (7)2.1.2一维情形下的散射概率函数的球谐函数展开 (9)2.1.3三维情形下的散射概率函数的球谐函数展开 (10)2.2中子输运方程的频率(能量)离散与多群近似 (12)2.3中子输运方程的方向离散与离散纵坐标近似 (13)2.3.1泄露算子 (13)2.3.2SN离散方法 (14)2.4中子输运方程的空间离散 (18)2.4.1迎风格式 (18)2.4.2菱形差分 (20)第三章中子输运方程的迭代求解 (21)3.1源迭代方法 (21)3.1.1不动点迭代 (21)3.1.2源迭代方法及理论基础 (22)3.2几类源迭代加速方法简介 (26)第四章源迭代非线性综合加速方法与Anderson加速方法及改进 (29)4.1非线性扩散综合加速方法 (29)4.2Anderson加速方法 (36)4.2.1Anderson加速不动点迭代 (36)4.2.2Anderson加速源迭代 (38)4.3AA-NDA方法 (40)4.4AA-NDA(2)方法 (43)第五章一维情形下的数值结果 (47)第六章三维情形下中子输运方程的求解 (67)6.1三维均匀矩形网格定态多群中子输运方程 (67)6.2三维非结构网格上输运方程及离散 (68)第七章三维情形下的数值结果 (71)展望 (77)参考文献 (79)致谢 (83)第一章绪论S1.1选题意义输运理论的发展已有一百多年的历史.1872年玻尔兹曼(L.Boltzmann)导出了分子分布函数随时间和空间变化的微分—积分方程,奠定了分子运动理论的基础.事实上,对于中子,光子,电子等一些粒子都可以导出类似的粒子守恒方程,称做输运方程.中子输运方程就是其中一类具体的方程[11,12,21,22,23].中子输运理论是研究中子在介质内运动的过程和规律的基础理论.随着核能利用的蓬勃发展,中子输运在核科学技术领域中已经成为一个独立的基础理论科学,并在核反应堆物理,屏蔽,核技术的工程应用和军事技术等领域中获得广泛的应用.近年来,随着计算机和计算方法的发展,数值方法已逐步成为中子输运理论的主要研究手段[28−30].因此对中子输运过程的数值模拟是最关键的问题之一.中子输运方程刻画中子群体在介质内的输运过程.中子输运方程的基本物理量是中子角通量,也称为中子密度.中子角通量是关于空间、时间、能量和运动方向相关的函数.因此,中子输运方程是与空间,时间,能量和方向相关的微分—积分方程.由于该方程具有很强的非线性,因此一般情况下只能进行数值求解.目前离散求解中子输运方程时,能量一般采用多群近似,空间一般采用差分法,运动方向采用离散纵坐标法(SN)[24−26,33].中子输运方程离散后是一个大型的非线性(线性)方程组,其数值求解过程非常复杂且所需计算量很大,耗费时间较长.对于离散后的中子输运方程,需要采用多个层次的迭代方法进行求解.其中源迭代是求解中子输运方程的关键算法[27].由于源迭代过程占据整个求解过程的主导计算量,且源迭代方法在某些情况下收敛非常慢,甚至失效.所以对源迭代的加速研究尤为重要.自1979年起出现了诸多关于源迭代的加速方法.如非线性扩散加速方法(Nonlinear Diffusion Acceleration NDA)[1,6,35−37],Anderson加速方法[2],角度多重网格(ANMG)加速方法[4],扩散综合加速方法[8],Chebyshev加速方法[9,10,31,32,34],粗网格再平衡法[11,12]等.近期有关源迭代加速方法研究,相关工作主要为矩估计加速(高低阶加速)方法.其中非线性扩散加速(NDA)已经成为求解中子输运方程(定源问题,k特征值问题)比较代表性的方法.NDA方法的主要思想就是利用低阶方程加速(更细网格上)高阶方程.NDA方法对于各向异性较弱问题的中子输运方程其加速效果很好,但对于各向异性比较强的情况NDA方法的加速效果明显弱化,甚至不收敛.近几年Anderson加速不动点迭代方法受到了许多关注.针对计算非线性的不动点问题,如果G是一个在解的邻域范围内的压缩映射,不动点问题进行逐次迭代后将会收敛.但收敛可能非常缓慢.Anderson加速是提高不动点迭代收敛的一类有效的方法.Anderson加速算法的主要思想是通过一组向量线性组合前k+1个迭代步中的近似解来表示最新的迭代步的值.其中,这个组合向量通过极小化每一步相应的残差得到.由于中子输运方程的源迭代可以看做是个不动点迭代,所以可以应用Anderson加速到中子输运方程的源迭代过程中.通过对NDA进行理论分析可以发现,当NDA加速各向异性较强的中子输运方程时,加速效果不好是因为在整个源迭代加速过程中仅仅0阶矩和1阶矩被加速,而其它高阶矩都没有被加速.在理论上高阶矩被扩展的低阶系统加速是可能的,但实际操作是非常复杂的,尤其在多维空间上更难进行.Anderson加速可以看作对每一个矩都进行了加速,所以可以考虑两者组合的加速方式.即在NDA的基础上进行Anderson加速得到新的算法AA-NDA.AA-NDA改进了NDA当各向异性增强时加速弱化的缺点,可预见的加速效果更加有效.同时由于Anderson加速相当于是在迭代外面的独立的加速过程,所以把NDA与AA结合操作起来不是很困难.近几年,在非结构网格的几何区域分解情况下,中子输运方程的并行计算成为了一个新的热点[13,14,15].由于SN方法在求解中子输运方程时在各个方向是独立的,所以在几何区域分解基础上,可采用并行的求解方法,在此基础上设计各个方向同时扫描的并行算法.源迭代加速算法的研究可以更好地提高中子输运方程的数值求解过程中的收敛速度,为核技术的工程应用,军事医学等领域的实际应用提供理论基础.S1.2中子输运方程中子输运方程的实质是中子在介质内迁移的守恒关系表达式.中子在介质内的运动轨迹和速度是个随机的过程,它的迁移是受周围原子核散射的结果.为了简化问题,我们做以下的假设[11]:(1)把中子看做是一个点粒子,即中子可以被速度和位置完全描述.(2)忽略中子与中子之间的碰撞.由于在一般的核技术问题以及核反应堆问题中,中子的密度比介质内的原子核密度小的很多,所以可以忽略中子与中子之间的碰撞.即把中子在介质内的输运过程主要看作是中子与介质原子核碰撞的结果.(3)因为中子不带电荷所以不受电磁的影响,所以认为中子在介质内两次碰撞之间穿行的路线是直线.(4)中子与介质内核的碰撞、发射可以认为是瞬时的.中子输运方程就是在上述假设下,描述中子群体在介质内的输运过程.求解中子输运方程可得到介质内中子密度分布函数.在输运理论中,中子是一个点粒子.所以中子在介质中的运动状态可以由速度、位置来确定.因此在任意时刻t,中子的运动状态由位置向量⃗r(x,y,z),能量E,和方向⃗Ω(θ,ϕ)六个变量来表示.其中对于不同的坐标系情况下,⃗r和⃗Ω的表示不同.对于单个中子,它是按照杂乱无章的折线轨迹在介质内随机运动的.每个中子运动不停直到中子从反应堆表面逸出或被吸收为止,整个过程是随机的.实际上我们研究的不是个别中子的轨迹和所在的位置,而是在空间不同点处中子密度的宏观期望分布.当中子密度足够大的时候,可以用一种处理大量中子行径的宏观理论推导出中子输运方程.研究中子输运过程的基本原则,就是中子数目守恒/中子平衡.即中子密度随时间的变化率等于它的产生率减去泄露率和移出率,也就是ðnðt=产生率(Q)−泄漏率(L)−移出率(R)中子输运方程的建立主要涉及以下几项:(1)泄漏项:L=⃗Ω·∇ϕ(⃗r,⃗Ω,E,t)(1.1) (2)移出项:R=σtϕ(⃗r,⃗Ω,E,t)(1.2) (3)产生项:Q=∫︁∞0∫︁4πσs(⃗r,⃗Ω′→⃗Ω,E′→E)ϕ(⃗r,⃗Ω′,E′,t)dE′d⃗Ω′+q(⃗r,E)(1.3)由守恒关系可以得到中子输运方程的表达式为:1 v ðϕðt+⃗Ω·∇ϕ+σt(E)ϕ=q(⃗r,E,⃗Ω)+∫︁∞∫︁4πσs(⃗r,E′→E,⃗Ω′→⃗Ω)ϕ(⃗r,E′,⃗Ω′,t)dE′d⃗Ω′(1.4)其中σt为输运总截面,σs(⃗r,E′→E,⃗Ω′→⃗Ω)为中子从能量E′方向⃗Ω′散射到能量E方向⃗Ω的散射截面,q为外源项.E为能量变量,⃗Ω(θ,ϕ)为方向变量,⃗r(x,y,z)为位置变量,t为时间变量.ϕ=ϕ(⃗r,⃗Ω,E,t)为角通量密度.当系统处于平衡状态(稳态)时,ðnðt=0.所以稳态下的中子输运方程为:⃗Ω·∇ϕ+σt ϕ=q(⃗r,E,⃗Ω)+∫︁∞∫︁4πσs(⃗r,E′→E,⃗Ω→⃗Ω′)ϕ(⃗r,E′,⃗Ω′)dE′d⃗Ω′(1.5)由此可见,稳态情形下中子输运方程是含有六个自变量的微分—积分方程.在数学上,这样的方程解析求解是十分困难的,只能数值求解.但由于输运过程的复杂性,即便是应用数值求解仍然是十分复杂的.中子输运方程是中子密度分布函数应该满足的方程.所以对于具有相同的参数σs,σt的反应堆,不论初始状态、以及边界条件如何,都应该满足同一个方程.但从数学上讲,由于中子输运方程求出的解中包含任意的积分常数,所以要给定适当的定解条件.即给出适当的边界条件以及初始条件来保证能给出问题的唯一解.由于中子输运方程的散射截面参数σs通常是在已知的角通量的情况下求得的,因此可以认为中子输运方程是非线性的.由于中子输运方程大规模、非线性、多变量的特点,只能对其进行数值求解.求解中子输运方程通常采用的方法是源迭代方法.由于方程的特点,源迭代收敛的速度非常的慢,甚至在一些情况下是不收敛的.国内外一直有很多关于中子输运方程源迭代加速方法的研究.但在很多情形下,例如在各向异性情形下,很多源迭代加速方法效果不够理想.所以对于源迭代的加速方法的研究以及不同方法的适当组合尤为重要.S1.3论文的主要工作本文是对中子输运方程的源迭代加速方法的研究.论文的主要工作是以Anderson 加速方法(AA)为基础的,结合非线性扩散综合加速方法(NDA)对源迭代过程进行加速研究.通过数值实验对比各个方法以及新方法的不同加速效果.针对三维非结构网格情形下,开展中子输运方程的AA加速研究.具体内容为:1.分析并推导中子输运方程的球谐函数展开形式及不同变量的各种离散方法.在此方程的基础上设计不同源迭代加速方法的数值实验.2.把Anderson加速方法,非线性扩散综合加速方法(NDA)应用到中子输运方程中.并给出相关的算法描述形式.3.以Anderson加速方法为基础与非线性扩散综合加速方法结合形成新的加速方法,包括AA-NDA方法,AA-NDA(2)方法.将这两种加速方法应用到中子输运方程的源迭代过程中,并给出相关的算法表达式.4.通过数值实验分析对比AA方法,NDA方法,AA-NDA方法,AA-NDA(2)方法对于源迭代方法的加速效果.5.推导三维非结构网格情形下的中子输运方程的离散过程,用源迭代方法和Anderson 加速方法求解三维非结构网格下的中子输运方程,并给出相关的算法描述形式.6.设计三维非结构网格情况下源迭代方法和Anderson加速方法的数值实验,在三维非结构网格情况下,验证Anderson加速方法对于源迭代过程的加速效果;设计并行情况下的数值实验,测试不同处理器核数情况下,AA方法和源迭代方法求解中子输运方程的数值结果.S1.4论文的组织结构本文分为七章.第一章绪论部分介绍了中子输运方程的背景,方程的求解方法及各类加速方法.第二章介绍了中子输运方程散射源的球谐函数展开及方程中各个变量的离散方法,其中能量离散采用多群近似方法;方向离散采用SN离散方法;空间离散采用菱形差分(迎风格式)方法.第三章介绍了中子输运方程的数值求解方法—源迭代方法,以及各类源迭代加速方法的研究现状.第四章介绍了Anderson加速方法(AA),非线性综合加速方法(NDA),以及在Anderson加速方法的基础上,结合NDA方法形成的AA-NDA方法和AA-NDA(2)方法.详细地描述了几类方法的算法形式,分析了几类方法的原理及特点.第五章介绍了一维情形下各类方法的数值实验结果,对比了AA方法,NDA方法,AA-NDA方法,AA-NDA(2)方法对源迭代的不同加速效果.第六章介绍了三维情形下中子输运方程的数值计算,重点介绍三维均匀矩形网格和三维非结构网格上输运方程的离散.第七章介绍了三维非结构网格情形下AA方法对源迭代过程的加速效果,且设计并行情形下的数值实验,测试不同处理器核数情况下,源迭代方法及AA加速的源迭代方法求解中子输运方程的数值结果.第二章中子输运方程及数值计算方法S2.1散射源的球谐函数展开中子输运方程是描述中子在系统内经过产生、迁移和消失达到的守恒关系.稳态情况下的中子输运方程为:⃗Ω·∇ϕ+σt (E)ϕ=q(⃗r,E,⃗Ω)+∫︁∞∫︁4πσs(⃗r,E′→E,⃗Ω′→⃗Ω)ϕ(⃗r,E′,⃗Ω′)dE′d⃗Ω′(2.1)其中:E为能量变量,⃗Ω(θ,ϕ)为方向变量,⃗r(x,y,z)为位置变量,t为时间变量,ϕ=ϕ(⃗r,⃗Ω,E)为角通量密度,q为外源项,σt(E)为中子总截面,σs(⃗r,E′→E,⃗Ω→⃗Ω′)为中子从能量E′,方向⃗Ω′散射到能量E,方向⃗Ω的散射截面.散射概率函数为方程右端的第二项,代表散射中子通量,如式(2.2)所示:∫︁∞0∫︁4πσs(⃗r,E′→E,⃗Ω′→⃗Ω)ϕ(⃗r,E′,⃗Ω′,t)dE′d⃗Ω′(2.2)接下来的章节将介绍对散射概率函数的具体展开过程.因为展开过程中会涉及到一些多项式展开的基础知识,所以下面我们将对这些所需的基础知识进行简单介绍. S2.1.1球谐函数与Legendre多项式对于一维情况下的散射概率函数,我们将对其中的散射截面σs以及中子角通量函数ϕ进行Legendre多项式展开.所以下面我们将介绍Legendre多项式的定义以及其正交性、递推性[16].定义2.1.1在区间[−1,1]上,取权函数为ρ(x)≡1.由{1,x,x2,...,x n,···}正交化得到的多项式为勒让德(Legendre)多项式,记为P n(x),其表达形式为:P0(x)=1,P n(x)=12n n!d ndx n(x2−1)n,n=1,2,...P n(x)的前几项为:P0(x)=1;P1(x)=x;P2(x)=(3x2−1)/2;P3(x)=(5x3−3x)/2;...Legendre多项式具有以下重要的两个性质.这两个性质在中子输运方程的变形推导中将会用到.性质2.1.1(正交性)任意两个Legendre多项式是正交的,即∫︁1−1P n(x)P m(x)dx=⎧⎪⎨⎪⎩0,m=n;22n+1,m=n.性质2.1.2(递推关系)Legendre多项式具有三项递推关系.(n+1)P n+1(x)=(2n+1)xP n(x)−nP n−1(x),n=1,2,...通过递推关系式,由P0(x)=1,P1(x)=x,可以得到P2(x)=(3x2−1)/2,P3(x)=(5x3−3x)/2,P4(x)=(35x4−30x3+3)/8,...接下来介绍伴随Legendre多项式P mn的表达形式及递推关系式,这在球谐函数的介绍中将会用到.定义2.1.2在球坐标系下求解拉普拉斯方程就可以得到伴随勒让德多项式,具体形式如下所示:P m n (x)=12n n!(1−x2)m/2d n+mdx n+m(x2−1)n m=0,1,...,n.伴随勒让德多项式的前三项为:P1 1(x)=(1−x2)1/2或者P11(cosθ)=sinθP1 2(x)=3(1−x2)1/2x或者P12(cosθ)=3sinθcosθP2 2(x)=3(1−x2)1/2或者P22(cosθ)=3sin2θ伴随勒让德多项式具有下述的递推关系:性质2.1.3(递推关系)与勒让德多项式类似,伴随勒让德多项式具有如下递推关系:xP mn (x)=12n+1[︀(n−m+1)P mn+1(x)+(n+m)P n−1(x)]︀(x2−1)dP mn(x)dx=nxP mn(x)−(n+m)P mn−1(x)或者sinθP mn =12n+1[︀(n−m+1)P mn+1+(n+m)P mn−1]︀sinθP mn=12n+1(P m+1n+1−P m+1n−1)对于三维情况下的散射概率函数,与一维情况类似,我们将对其中的散射截面σs以及中子角通量函数ϕ进行球谐函数多项式展开.所以下面我们介绍球谐函数及其相关的的定义以及正交性.定义2.1.3形如以下形式的函数称为球谐函数:Y n,m(⃗Ω)=Y n,m(θ,ϕ)=⎧⎨⎩P|m|n(cosθ)sin|m|ϕ,m=−1,...,−n;P mn(cosθ)cos mϕ,m=0,1,2,...,n,n=0,1,...其中P mn为伴随勒让德多项式,向量⃗Ω由θ,ϕ两个变量表示.球谐函数具有以下正交性的特点.这个性质在中子输运方程的变形推导中将会用到.性质2.1.4球谐函数Y n,m(⃗Ω)的正交性:∫︁ΩY n,m(⃗Ω)Y n′,m′(⃗Ω)d⃗Ω=⎧⎪⎨⎪⎩0,n′=n或m′=m;2π2n+1(n+|m|)!(n−|m|)!(1+δm),n′=n且m′=m.其中δm为克罗内克尔符号δm 0=⎧⎨⎩0,m=0;1,m=其他值.S2.1.2一维情形下的散射概率函数的Legendre多项式展开中子输运方程中散射项的计算比较复杂.散射项使得各个方向,不同频率的中子的计算耦合在一起.为了方便算法讨论与设计,可以将散射项用球谐函数(Legendre 多项式)进行展开.下面讨论散射项的球谐函数(Legendre多项式)展开方式.为了介绍方便,本小节先介绍空间一维情形下,散射项的Legendre多项式的展开.在一维情况下,式(2.1)中的方向变量⃗Ω是由μ=cosθ来唯一表示.假设散射概率函数只跟散射前后的运动方向之间的夹角μ0=⃗Ω·⃗Ω′有关,则一维平面稳态中子输运方程可表示为:μ·ðϕðx+σt(E)ϕ=q(x,E,μ)+∫︁∞∫︁4πσs(x,E′→E,μ0)ϕ(x,E′,μ′)dE′dμ′(2.3)式中散射概率函数为:∫︁∞0∫︁4πσs(x,E′→E,μ0)ϕ(x,E′,μ′)dE′dμ′(2.4)对于一维问题,利用勒让德多项式展开散射概率函数中的散射截面σs,有:σs(x,E′→E,μ0)=∞∑︁l=02l+12σsl(x,E′→E)P l(μ0)(2.5)其中P l(μ0)为勒让德多项式,σsl(x,E′→E)为展开系数.展开系数由勒让德多项式的正交性求得.为计算具体的表达式,在(2.5)式两边分别乘以勒让德多项式P l(μ0),由正交性可得展开系数σsl的表达式,如下:σsl(x,E→E′)=∫︁1−1σs(x,E′→E)P l(μ0)dμ0再利用勒让德多项式展开角通量函数ϕ(x,E,μ),可得:ϕ(x,E,μ)=∞∑︁l=02l+12φl(x,E)P l(μ)(2.6)其中φl(x,E)为ϕ(x,E,μ)的第l阶矩通量.为计算具体的表达式,在(2.6)式的两端分别乘以勒让德多项式P l(μ),由勒让德多项式的正交性,得到矩通量φl的表达式,如下:φl(x,E)=∫︁1−1ϕ(x,E,μ)P l(μ)dμ将(2.5)式,(2.6)式代入散射函数(2.4)式中,得到Legendre多项式展开后的散射函数,如下:∫︁∞0∫︁4πσs(x,E′→E,μ0)ϕ(x,E′,μ′)dE′dμ′=∫︁∞0∫︁4π∞∑︁l=02l+12σsl(x,E′→E)P l(μ0)×∞∑︁l=02l+12φl(x,E)P l(μ)dE′dμ′=∫︁∞∑︁l=02l+12σsl(x,E→E′)φl(x,E)P l(μ)dE′(2.7)将(2.7)式代入到(2.3)式中,得到散射源Legendre多项式展开后的一维中子输运方程:μ·ðϕðx+σt(E)ϕ=q(x,E,μ)+∫︁∞∞∑︁l=02n+12σsl(x,E→E′)φl(x,E)P l(μ)dE′(2.8)S2.1.3三维情形下的散射概率函数的球谐函数展开三维情形下,散射概率函数的球谐函数展开方式与一维类似.仍假设散射概率函数只跟散射前后的运动方向之间的夹角μ0=⃗Ω·⃗Ω′有关,即散射概率函数可以表示为:∫︁∞0∫︁4πσs(⃗r,E′→E,μ0)ϕ(⃗r,E′,⃗Ω′)dE′d⃗Ω′(2.9)其中的散射截面可以表示为:σs(⃗r,E′→E,⃗Ω→⃗Ω′)=σs(⃗r,E′→E,μ0)利用勒让德多项式展开上式:σs(⃗r,E′→E,μ0)=∞∑︁l=02l+12σsl(⃗r,E′→E)P l(μ0)(2.10)其中P l(μ0)为勒让德多项式,σsl(⃗r,E′→E)为展开系数.在(2.10)式两边分别乘以P l(μ0),由Legendre多项式的正交性可得展开系数σsl的表达式,如下:σsl(⃗r,E→E′)=∫︁1−1σs(⃗r,E′→E,μ0)P l(μ0)dμ0进一步利用球谐函数展开角通量函数ϕ(⃗r,E,⃗Ω),有ϕ(⃗r,E,⃗Ω)=ϕ(⃗r,E,θ,ϕ)=∞∑︁n=0(2n+1)4π∞∑︁m=−nφn,m Y n,m(θ,ϕ)(2.11)其中Y n,m为球谐函数,φn,m为展开系数,由下式决定:φn,m=∫︁⃗Ωϕ(⃗r,E,⃗Ω)Y(⃗Ω)d⃗Ω图2.1:能量分群示意图将(2.10)式,(2.11)式分别带入(2.9)式中,利用勒让德多项式加法公式及正交性可以得到:∫︁∞0∫︁4πσs(⃗r,E′→E,⃗Ω→⃗Ω′)ϕ(⃗r,E′,⃗Ω′)dE′d⃗Ω′=∞∑︁n=0(2n+1)4π∞∑︁m=−nY n,m(θ,ϕ)∫︁dE′σsl(⃗r,E′→E)φn,m(2.12)将式(2.12)代入到(2.1)式中,得到散射源球谐函数展开后的三维中子输运方程:⃗Ω·∇ϕ+σt (E)ϕ=q(⃗r,E,⃗Ω)+∞∑︁n=0(2n+1)4π∞∑︁m=−nY n,m(θ,ϕ)∫︁dE′σsl(⃗r,E′→E)φn,m(2.13)S2.2中子输运方程的频率(能量)离散与多群近似数值求解中子输运方程之前,要对各个维度的变量进行离散[22].由于中子输运方程的多变量(能量,方向,空间)的特点,所以中子输运方程的变量离散十分重要.本小节我们研究能量变量的离散过程.对于能量变量,我们一般采用多群近似进行离散.多群近似是对连续的能量范围[E min,E max]进行若干的划分,其中每一个区间称为一个群,如图2.1所示.为了方便介绍能量变量的离散过程,减少数学上的麻烦,便于了解方法的步骤及形式,我们以一维情形下的中子输运方程为例介绍多群近似的过程.上一小节我们得到了散射源球谐函数展开后的一维中子输运方程:μ·ðϕðx+σt(E)ϕ=q(x,E,μ)+∫︁∞∞∑︁l=02n+12σsl(x,E→E′)φl(x,E)P l(μ)dE′(2.14)。

中子输运方程和扩散方程区别

中子输运方程和扩散方程区别

中子输运方程和扩散方程区别摘要:一、引言二、中子输运方程和扩散方程的定义及基本原理三、中子输运方程和扩散方程的区别1.适用范围2.物理意义3.数学形式4.求解方法四、实际应用案例五、结论正文:一、引言中子输运方程和扩散方程都是在核物理和核工程领域中具有重要意义的方程。

它们在描述中子在物质中的行为方面具有密切的联系,但又有明显的区别。

本文将详细阐述这两者之间的区别,并介绍各自的适用范围、物理意义、数学形式和求解方法。

二、中子输运方程和扩散方程的定义及基本原理1.中子输运方程:中子输运方程是描述中子在物质中传播和散射过程的偏微分方程。

它反映了中子在物质中的空间分布和能量变化。

中子输运方程基于neutronBoltzmann 方程推导而来,适用于中子在物质中的各种输运过程。

2.扩散方程:扩散方程是描述中子在物质中由于碰撞引起的能量和方向变化的过程。

它主要关注中子在物质中的扩散行为,反映了中子在物质中的传输特性。

扩散方程基于Fick 定律推导而来,适用于中子在物质中的扩散过程。

三、中子输运方程和扩散方程的区别1.适用范围:中子输运方程适用于描述中子在物质中的各种输运过程,包括扩散、散射等。

扩散方程则主要关注中子在物质中的扩散行为。

2.物理意义:中子输运方程反映了中子在物质中的空间分布和能量变化,强调了中子的宏观输运特性。

扩散方程则关注中子在物质中的扩散过程,体现了中子在物质中的微观行为。

3.数学形式:中子输运方程是一偏微分方程,描述中子在物质中的空间分布和能量变化。

扩散方程则为一组微分方程,描述中子在物质中的扩散过程。

4.求解方法:中子输运方程的求解方法主要有数值方法、解析方法等。

扩散方程的求解方法主要有稳态法、非稳态法等。

四、实际应用案例:1.中子输运方程:在核反应堆设计中,中子输运方程用于预测中子在反应堆中的分布情况,以确保反应堆的安全和高效运行。

此外,在核燃料棒设计和核辐射防护方面也有广泛应用。

2.扩散方程:在核燃料元件设计中,扩散方程用于预测中子在燃料元件中的扩散行为,以优化燃料元件的结构和材料。

第9章 中子输运方程-分布的平衡

第9章 中子输运方程-分布的平衡
下面给出输运方程的简要推导,通过推导我们会知道它只是一种平衡关系。当然你也可 以进行更加精细的推导,包括量子力学公式,但是我们认为这对于理解输运方程的基本物理 概念是不必要的。中子“数密度”定义如下:
n(r, v, t)d 3rd 3v ≡ 在时间 t , r 处体元 d 3 r ,速度在 v 处 d 3 v 内, 中子的期望值
和 φ (r) = ∫φ(r, E)dE
(9.17)
这是研究中子扩散的起点。为了简单起见,今后将不再标出上划线,但是(9.16)式的原型应 该牢记。我们在第 10 章会再次研究这个方程的不同解法,包括对于边界条件的处理。注意 到(9.16)反映的仍然是一个平衡关系——从一个体元流入和流出的差值(扩散)与同一个体 元中产生和吸收的差值所形成的平衡。
Ⅱ. 中子输运方程——介绍
中子输运方程是在研究中子在介质中相互作用和迁移时最基本的方程。对于结构——速 度相空间中的中子,这个方程及方程的解是一个随时间变化的分布函数。得到这个分布函数 就可以解决在反应堆理论中几乎所有感兴趣的问题。但是,通常并不需要知道分布函数本身, 知道在相空间中的坐标如速度方向、能量或位置的积分就已经足够了。
(9.20)
F (E' Ω' → E, Ω) = F (Ω' → Ω)δ (E − E0 )
(9.21)
f (E) = Σ , 常数
(9.23)
这里δ ( x) 是狄拉克 δ 函数——它的值在 x ≠ 0 处皆为 0,在 x = 0 的值是无穷大(当被积分 区域包括这个点时, δ 函数的积分等于 1)。以下是 δ ( x) 函数的一些性质:
(2) 散射
∫∫∫ ΣS (E' )φ(r, E', Ω', t)d 3rdE'dΩ'∆tF (E'Ω' → EΩ)dEdΩ

中子输运方程和扩散方程区别

中子输运方程和扩散方程区别

中子输运方程和扩散方程区别摘要:1.中子输运方程和扩散方程的定义与含义2.中子输运方程和扩散方程的物理背景与应用领域3.中子输运方程和扩散方程的数学表达式及求解方法4.中子输运方程和扩散方程的区别与联系5.泄漏迭代法在求解中子扩散方程中的应用正文:一、中子输运方程和扩散方程的定义与含义中子输运方程和扩散方程都是物理学中描述粒子传输过程的方程。

中子输运方程主要应用于中子在物质中的输运过程,而扩散方程则广泛应用于粒子在各种介质中的扩散现象。

二、中子输运方程和扩散方程的物理背景与应用领域中子输运方程主要用于研究中子在核反应堆中的传输过程,对于核反应堆的设计、仿真和安全验证具有重要意义。

扩散方程则广泛应用于粒子在气体、液体和固体等介质中的扩散现象,如气体分子的扩散、污染物在环境中的扩散等。

三、中子输运方程和扩散方程的数学表达式及求解方法中子输运方程的数学表达式通常是基于积分形式的,描述了中子在物质中的输运过程。

求解方法主要有常微分方程求解法、有限元法等。

而扩散方程的数学表达式则是基于偏微分方程的,描述了粒子在介质中的扩散现象。

求解方法包括经典数值解法、有限差分法等。

四、中子输运方程和扩散方程的区别与联系中子输运方程和扩散方程在物理背景、应用领域和数学表达式上都有所区别,但它们都是描述粒子传输过程的方程,具有一定的联系。

在实际应用中,可以根据问题的具体特点选择合适的方程进行求解。

五、泄漏迭代法在求解中子扩散方程中的应用泄漏迭代法是一种求解中子扩散方程的有效方法,通过迭代计算可以逐步逼近中子扩散方程的解。

该方法在核反应堆物理计算等领域具有广泛的应用,对于提高计算精度和效率具有重要意义。

总结:中子输运方程和扩散方程是描述粒子传输过程的两种重要方程,它们在物理背景、应用领域和数学表达式上有所区别,但也具有一定的联系。

在实际应用中,可以根据问题的具体特点选择合适的方程进行求解。

线性源近似的中子输运方程特征线解法

线性源近似的中子输运方程特征线解法

Li a o r e App o i a i n S he e f r M e h d o a a t rs is ne r S u c r x m to c m o t o f Ch r c e i tc
TANG Chun t o —a
( angh uce Sh aiN l arEngi e r n Re e r h and D e i I tt t n e i g s a c sgn nsiu e,Shan ghai2 23 00 3, Chi a) n
is o r s ond n mo fc ton o ne tve o r e it i to we e t c r e p ig diia i f r ga i s u c d s rbu i n r pr po e . The o sd
OE CD/ NEA G7 M OX D e c ma k a d as l— e i e C5 一 2 b n h r n efd f d BW R i i o e p o l m r n m n - r r b e we e c
Absr c : M e ho of ha a t r s is ( OC) f s l i n ut o ta po t qu to ta t t d c r c e itc M or o v ng e r n r ns r e a in
b s d o u t uc ur d a e n ns r t e m e h ha le d c m e ne o t e und m e t l e h ds o s s a r a y be o o f h f a n a m t o f r
i a o h r b sc a s m p i n f r M 0C. A i e r s u c ( ) c a a t rs is s h me a d s n t e a i s u to o l a o r e LS n h r c e itc c e n

中子输运方程和扩散方程区别

中子输运方程和扩散方程区别

中子输运方程和扩散方程区别1. 物理意义中子输运方程和扩散方程都是描述粒子(中子)在介质中传播的方程,但它们有着不同的物理意义。

中子输运方程描述的是中子在介质中由于碰撞和扩散作用而产生的输运过程。

它涉及到中子的速度分布、通量分布和中子密度随时间、空间的变化。

中子输运方程是概率密度函数的时间演化和空间演化的耦合,描述了中子在空间和时间上的分布变化。

扩散方程则描述的是粒子的扩散过程,即粒子从高浓度区域向低浓度区域的传播。

它涉及到粒子的浓度分布、通量与扩散系数之间的关系以及扩散过程中的各项热力学参数。

扩散方程是浓度梯度驱动的方程,描述了粒子分布的空间变化。

2. 数学形式中子输运方程和扩散方程在数学形式上也有所不同。

中子输运方程的一般形式为:∂ρ/∂t + div(ρvΦ) + ∇·(ρvε NBC) = ρsterdam蹋U,,式中ρ为中子密度,v为中子速度,Φ为通量,N为宏观因子,C为弹性截面,ε为源项。

这个方程包括了中子的时间演化、空间扩散和产生-吸收等过程。

扩散方程的数学形式为:∂c/∂t = D ∇²c + f(c)其中,c为粒子的浓度,D为扩散系数,f(c)为反应项,描述了粒子浓度的变化。

这个方程仅描述了粒子的扩散过程,没有考虑到粒子的产生-吸收等过程。

3. 边界条件中子输运方程和扩散方程在边界条件上也有所不同。

中子输运方程的边界条件通常需要考虑中子的入射、反射和泄漏等情况,具体形式可以为:-div(ρvΦ) = ρ_s - ρ_r ,边界上中子的入射通量和泄漏通量等于中子的反射通量和中子在边界上的产生量之和。

扩散方程的边界条件通常需要考虑粒子在边界上的流入流出情况,具体形式可以为:c(x=0,t)=c_0 ,边界上的粒子浓度等于初始浓度。

-D ∂c/∂x |_{x=L} = Q ,边界上的粒子通量等于粒子产生量减去粒子吸收量。

4. 解法中子输运方程和扩散方程的解法也有所不同。

快中子与慢中子在材料中的输运特性

快中子与慢中子在材料中的输运特性

快中子与慢中子在材料中的输运特性快中子与慢中子在材料中的输运特性是研究核材料与辐射防护领域的一项重要课题。

中子具有无电荷、质量较大、穿透力强等特点,因此在核材料的相互作用以及辐射的传输过程中起到至关重要的作用。

了解中子在材料中的输运特性,对于核材料的设计、辐射防护以及核能利用具有重大意义。

快中子是指能量较高的中子,其能量通常在几百千电子伏特到几十兆电子伏特之间。

快中子在材料中的输运特性主要受到散射、吸收和漫反射等因素的影响。

在材料中,快中子与原子核的碰撞主要通过库仑散射和核子弹性散射完成。

库仑散射是指快中子与原子核之间的库仑力相互作用,当中子与原子核发生反弹碰撞时,快中子的能量会减少。

核子弹性散射是指快中子与原子核碰撞后,两者都保持动能守恒,但方向发生改变。

在材料中,快中子还会与材料中的杂质原子发生碰撞,从而引起能量的散失和方向的改变。

慢中子是指能量较低的中子,其能量通常在几个电子伏特以下。

慢中子在材料中的输运特性主要受到吸收、散射和俘获等因素的影响。

在材料中,慢中子与原子核的碰撞主要通过弹性散射、非弹性散射和共振吸收等过程完成。

弹性散射是指慢中子与原子核碰撞后,两者都保持能量和动量守恒,但方向发生改变。

非弹性散射是指慢中子与原子核碰撞后,两者的能量和动量都会发生改变。

共振吸收是指慢中子与具有特定能级的原子核共振发生吸收,这种共振吸收对于核材料的设计和利用具有重要意义。

快中子和慢中子的输运特性不同主要是由于它们能量的差异导致的。

对于快中子而言,其能量较高,因此在材料中的输运过程中会有较多的库仑散射和核子弹性散射,同时还会受到杂质原子的散射的影响。

相比之下,慢中子能量较低,因此在材料中的输运过程中主要是弹性散射、非弹性散射和共振吸收等过程起主导作用。

这种不同的输运特性决定了中子在材料中的传输方式和传输效果。

对于核材料的设计,需要根据其所需的中子束流特性,选择适合的中子源和材料,以实现所需的核反应。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

当A>10时,上式近似简化为
可见:当质量数大于10时,平均对数损耗是由靶核的质量数A决定的,与中子的能量近乎无关。
2.3 散射过程
3)慢化能力与慢化剂的选择
• 从中子慢化的角度,慢化剂应为轻元素,具有大的平均对数能降ξ; • 同时,慢化剂应具有较大的散射截面Σs; • 小的吸收截面Σa; 通常把乘积ξΣs叫作慢化剂的慢化能力。 ξΣs/Σa,叫作慢化比,它是表明慢化剂优劣的一个重要参数。
二、 中子与介质相互作用的物理过程
2.1 裂变过程
当一个中子与高原子序数的铀或钚等同位素发生碰撞时可能引发裂变反应,使之分裂成质量 相近的两个裂变碎片。由于裂变碎片的中子-质子比偏高而不稳定,通过放出一个到几个中子使其 回到稳定态;或者将一个中子转变为一个质子,同时放出一个β粒子使其达到稳定。在235U的裂变 过程中,每吸收一个中子,平均放出2.5个中子(通常称为 值),同时释放出大约200MeV的裂 变能。 裂变过程释放出的中子分为两种:一种是在极短的时间之内(约10-15s)放出的,通常称为瞬 发中子,占全部裂变中子的99%以上;另一种是缓发中子,它是在瞬发中子停止发射后,伴随裂 变产物的β衰变陆续发射的,要持续几分钟。通常将缓发中子按不同的衰变期、能量、份额分为6 组,其总份额称为β。对235U的裂变,β=0.0075,对239Pu则β要小得多,Am,Cm等锕系核素的β值 更小。
2.2 辐射俘获过程
中子与原子核发生碰撞时,若入射中子动能比靶核的每一核子平均作用能小时,就把入射中 子看作与整个核发生作用。这时核反应分为两个阶段,首先入射中子被靶核吸收形成复合核,处 于激发态的复合核具有的激发能等于被俘获中子的动能加上中子的结合能;接着该复合核以发射 粒子的方式退激发返回到基态。这叫做俘获反应。 低能中子与靶核可能发生如下的俘获反应: 放出γ射线的辐射俘获反应(n,γ)、放出α粒子的(n,α)反应、放出质子的(n,p)反应、 裂变即(n,f)反应。只有少数轻核能与中子发生(n,α)和(n,p)反应,只有较重的裂变核 才能与中子发生(n,f)反应。 因此,对反应堆中低能的中子来说,除了(n,f)之外,与(n,γ)相比,发生(n,α)和 (n,p)的概率要小得多。
一、 一般分析
2. 发展简史
• 最早的粒子输运理论是以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方法); • 在上述方法的基础上,产生了大批应用程序软件。
2.3 散射过程
时,入射中子的能量大于靶核的第一激发态和基态之间相差的能量, 就可能使靶核受到激发,发生非弹性散射反应。当具有某动能的入射中子与散射核发生非弹性散 射时,该中子先是被散射核吸收形成复合核,然后释放出一个动能较低的中子。入射中子损失的 那部分动能转变成散射核的激发能,使其处于激发态。随后散射核放出γ射线,回到基态。发生这 种非弹性散射的先决条件是中子的能量必须大于核的基态到第一激发态的能级间隔。中重核这个 能级间隔约为0.1MeV,轻核的要大得多。因此只有较高能量的中子才能引起轻核的非弹性散射反 应。
好的慢化剂不但要具有较大的慢化能力,还要有较大的慢化比
2.3 散射过程
从表中可以看出: • 水的慢化能力最大→轻水堆具有较小的堆芯体积,水吸收截面较大→水堆用富集铀; • 重水具有良好的慢化性能,可用天然铀,但价格昂贵(CANDU); • 石墨慢化性能较好,但慢化能力较小→石墨堆体积较庞大(HTGR); • 另外,慢化剂的选择还应从其它角度考虑如,辐照性能、价格等。
2.3 散射过程
一个中子与原子核发生碰撞后仍然释放出一个中子,只不过出射中子的动能低于入射中子的 动能,中子损失的这部分能量被靶核得到了。这就是散射过程。散射过程又可分为弹性散射和非 弹性散射两种。它们在反应堆物理计算中扮演重要角色。 1.弹性散射 如果中子与靶核发生碰撞时,入射中子的能量小于靶核的第一激发态和基态之间相差的能量, 就不可能使靶核受到激发,只能发生弹性散射反应。发生弹性散射时质心系内中子的动能不变。 转换为实验室系,则中子因碰撞靶核而损失了动能,这部分动能并没有使靶核受到激发,而成为 静止靶核的反冲能。弹性散射过程可以看作是“弹性球”式的碰撞。中子将其能量转移给散射核的 大小决定于中子与靶核的质量比和中子散射的角度。当散射角相同时,散射核的质量与中子的质 量几乎相等,若快中子与其发生对心碰撞,可将其全部能量传递给氢核而成为热中子。
2.3 散射过程
非弹性散射的几个特点:
• • • • • 只有快中子才能产生非弹性散射。 主要发生在中子发射后的10-8~10-7s时间间隔内。 重核比较容易激发非弹性散射,发射的γ射线能量也较低。 中子能量高且靶核质量大,非弹性散射截面就大; 同位素中子源发射的中子能量低 ,发生非弹性散射的概率小,可以忽略;加速器中子源发射的 中子能量为14MeV,较易发生非弹性散射。
2.5 中子的平均寿命
1. 扩散时间 快中子慢化成热中子后,将在介质内扩散一段时间。定义无限介质内热中子在自产生至被俘获以 前所经过的平均时间,称为扩散时间/热中子平均寿命,用 t 表示。在无限均匀的介质中,热中子 寿命在数值上等于在该介质中热中子的平均扩散自由程与热中子平均速度的比值。
平均扩散自由程是热中子从产生到被吸收为止的自由行程平均值,在原子核物理学中有:
2.3 散射过程
1)碰撞前后中子的动能比值E2/E1 每次弹性碰撞后,快中子的能量损失与靶核的质量数A、入射中子的能量E1以及散射角θ有关。 由动量守恒和能量守恒定律,可得
对上式讨论: 当θ=0o时,中子动能没损失,E2最大,E2=E1。 当θ=180o时,即发生正碰撞时,中子的动能损失最大,E2最小,有
2.5 中子的平均寿命
3. 中子的平均寿命 在反应堆动力学计算中需要用到快中子自裂变产生到慢化成为热中子,直至最后被俘获的平 均时间,称为中子的平均寿命,用l 表示。显然
l t ts
• • • • • 对于热中子反应堆,中子平均寿命主要由热中子的平均寿期,即扩散时间决定; 对于压水堆,中子的平均寿期l≈10-4s; 对于快中子反应堆,中子的平均寿期l≈10-7s; 上面的计算基于无限介质,没有考虑泄露的影响; 对于具有实际大小的反应堆系统,计算时应考虑泄露的影响,对中子平均寿命进行修正,所得 寿命将比不计泄露时的短。
一、 一般分析
1. 输运过程与输运理论
由于中子与原子核的无规则碰撞,中子在介质内的运动是一种杂乱无章的运动,即初始在堆内 某一位置具有某种能量及某一运动方向的中子,稍晚些时候,将运动到堆内另一位置以另一能 量和另一运动方向出现。这一现象称为中子在介质内的输运过程(Transport)。描述这一过程 的精确方程为玻尔兹曼输运方程(Boltzmanm equation)。 输运理论:围观粒子(中子、光子、电子、离子和分子等)在介质中的迁移统计规律的数学理 论:不是研究个别粒子的运动,而是研究大量粒子运动所表现的非平衡统计运动规律。
于是有:
2.5 中子的平均寿命
2. 慢化时间 慢化时间ts:在无限介质内,裂变中子由裂变能E0慢化到热能Eth所需要的平均时间。
s ts 2
1 1 E0 Eth
2.5 中子的平均寿命
下图是几种热中子反应堆慢化剂的慢化时间及扩散时间
慢化时间:10-4~10-6s; 扩散时间:10-2~ 10-4s
2.4 其它核反应过程
在核反应堆中除发生上述主要的核反应之外,还可以发生其它核反应,例如(n,2n),(n, 3n),(n,p),(n,d),(n,t),(n,3He),(n,α)等。其共同特点是要求入射中子 能量比较高,而且截面数值较小。前两种反应是增殖中子的。一般来说,只有入射中子能量在几 个MeV以上才能发生(n,2n);发生(n,3n)则要求入射中子具有10MeV左右甚至更高的能量, 而且(n,3n)的截面要比(n,2n)的小得多。所以在裂变对的物理研究中,(n,3n)的作用 要比(n,2n)的小,常把它忽略。但对于聚变堆来说,系统中的高能中子显著增多,发生这两 种核反应的可能性大得多。发生(n,2n)反应的最重要的核素当属9Be,因为它的(n,2n)反应 截面较大,常用做聚变堆的反射层材料。 上述其它5种核反应的出射粒子都是带电粒子,通常称为中子引起的带电粒子核反应。总的 来说,中重核的带电粒子反应截面都是比较小的。它们在反应堆临界计算中通常被忽略掉。
一、 一般分析
中子输运理论的基础假设是把中子看做一个点粒子,因此可以用其所在位置、具有的动能、 飞行方向进行精确描述。该假设的合理性在于中子的约化波长比介质的宏观尺寸和中子的平均自 由程要小得多。从德布罗意方程出发可以得到证明。 中子的输运理论是描述大量中子与介质相互作用的规律,它是中子与核发生各种碰撞过程的 综合描述。这是中子输运理论的基础。输运理论通常也称作迁移理论。求解某系统的中子输运方 程,可以得到该介质内中子随能量、空间、角度变化的规律。要进行这样的研究或计算,必须已 知单个中子与多个核素发生单个碰撞过程的大量信息。换句话说,大量的微观核反应数据是核反 应堆理论计算的数据基础。当核反应堆的体积较大时,从任何发生裂变反应产生中子的地方到系 统外边界的距离比介质的中子平均自由程大得多时,用扩散理论也能得到较好的结果,例如核电 站等大型功率堆得物理计算多采用扩散理论。
相关文档
最新文档