中子的输运
一维球几何中子输运方程计算方法研究的开题报告
一维球几何中子输运方程计算方法研究的开题报告
题目:一维球几何中子输运方程计算方法研究
一、研究背景和意义
中子输运是核能反应堆中的重要现象之一,其中中子的漂移和散射
是最基本的现象。
而针对中子输运的计算方法是研究这些现象中最为关
键的一个方面,对于准确模拟中子输运行为和核反应堆的设计分析具有
重要的意义。
目前,一维球几何中子输运方程是求解中子输运运动学的
基本方程,在核反应堆的瞬态和稳态计算中被广泛应用。
因此,对一维
球几何中子输运方程计算方法进行深入的研究是十分必要的。
二、研究内容
本研究将重点研究一维球几何中子输运方程的计算方法,主要包括
以下内容:
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章 中子输运方程-分布的平衡
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. 前言 (3)2. 蒙特卡洛方法概述 (3)2.1 蒙特卡洛方法的基本思想 (4)2.2 蒙特卡洛方法的收敛性、误差 (4)2.2.1 蒙特卡洛方法的收敛性 (4)2.2.2 蒙特卡洛方法的误差 (5)2.3 蒙特卡洛方法的特点 (6)2.4 蒙特卡洛方法的主要应用范围 (7)3. 随机数 (7)3.1 线性乘同余方法 (9)3.2 伪随机数序列的均匀性和独立性 (9)3.2.1 伪随机数的均匀性 (9)3.2.2 伪随机数的独立性 (10)4. 蒙特卡洛方法在粒子输运上的应用 (10)4.1 屏蔽问题模型 (10)4.2 直接模拟方法 (11)4.2.1 状态参数与状态序列 (11)4.2.2 模拟运动过程 (12)4.2.3 记录结果 (15)4.3 蒙特卡洛方法的效率 (16)5. 蒙特卡洛方法应用程序—MCNP (17)5.1 MCNP简述 (17)5.2 MCNP误差的估计 (18)5.3 MCNP效率因素 (19)6. 结论 (19)参考文献 (20)1.前言半个多世纪以来,由于科学技术的发展和电子计算机的发明,蒙特卡洛(Monte Carlo)方法作为一种独立的方法被提出来,并首先在核武器的试验与研制中得到了应用。
蒙特卡洛方法是一种计算方法,但与一般数值计算方法有很大区别。
它是以概率统计理论为基础的一种方法。
由于蒙特卡洛方法能够比较逼真地描述事物的特点及物理实验过程,解决一些数值方法难以解决的问题,因而该方法的应用领域日趋广泛。
蒙特卡洛模拟计算是解决中子在介质中输运较为成熟、有效的方法,对于原子能、辐射防护、剂量学和辐射生物物理学等研究领域实际问题的计算,都可以利用蒙特卡洛方法予以实现。
粒子输运过程可以用玻耳兹曼方程加以描述,然而,以此基础上发展起来的近似数值方法如扩散近似法、离散坐标方法在处理截面与能量相关以及散射各向异性介质、复杂几何条件问题时碰到了较大困难。
精细化瞬态中子输运计算的多级加速理论
精细化瞬态中子输运计算的多级加速理论汇报人:日期:CATALOGUE目录•精细化瞬态中子输运计算方法•多级加速理论•精细化瞬态中子输运计算应用•精细化瞬态中子输运计算展望•多级加速理论在精细化瞬态中子输运计算中的实践与优化01精细化瞬态中子输运计算方法将中子能量分成多个区间,每个区间对应不同的反应率模型。
能量分组时间离散空间离散将时间进行离散化,用有限个时间步长来近似连续的时间演化。
将空间进行网格划分,用有限个网格点来近似连续的空间分布。
03精细化中子输运方程0201用差分算子近似微分算子,将连续的输运方程转化为离散的差分方程。
离散化方程方法有限差分法将空间划分为有限个单元,用分片常数近似连续的反应率函数,将连续的输运方程转化为离散的有限元方程。
有限元法通过随机抽样中子在每个时间步长的运动轨迹,统计中子的数量和能量分布,从而得到中子的输运过程。
蒙特卡洛方法多重网格加速法利用不同尺度的网格来加速求解输运方程,通过在不同尺度的网格上迭代求解输运方程,提高求解速度。
直接求解法直接求解输运方程,得到每个时间步长的中子分布。
并行化加速法将计算任务分配到多个处理器上并行计算,通过同时处理多个计算任务来提高求解速度。
输运方程的加速求解方法02多级加速理论算法原理多级加速算法基于分治策略,将大规模问题分解为若干个小规模子问题,并逐级加速求解。
应用场景精细化瞬态中子输运计算领域存在计算量大、计算时间长的问题,多级加速算法可针对性地优化和加速。
理论背景多级加速算法是一种针对大规模科学计算问题的优化方法,旨在提高计算效率。
将原始大规模瞬态中子输运问题分解为多个较小规模的子问题,并逐级求解。
逐级分解针对每个子问题采用适当的加速方法,如并行计算、GPU加速等,以实现计算效率的提升。
加速方法通过逐级分解和加速,降低计算时间和内存占用,提高整体计算效率。
逐级加速效果1加速效果评估23采用标准测试问题对多级加速理论进行验证和评估,包括计算时间、内存占用、收敛速度等方面。
高速运动介质中中子输运的新处理方法
时
,
一 般把 右 边 各 项 中 的 V 及 V , 的函 数 展开 成
。
及
。
的 函数
使 方 程两 边 统 一 出 现 中 子 在 空 间 运动 的 速 度 变 量
,
但是
,
。
,
实际 上
,
由于 各 种 相 互 作 用 截 面 和 散射 指 示 函 数 不 能 简 单 地 写 成 V 及 砂 的 解 析 函 数
3
;
J “
一 ” `’
、
l
=
丁蔽了 ) 气
仍
\ 了 “
“ x
p
〔一 丽了尸
-
f
哪
;
L ”
一 `
,
-
一
( 4 )
哪
按照 〔
1 〕中作 法
,
不 难看 出 中子 输 运 方程 (
2 )
将改 换成
本文 1 ,
a l年 3
月 16 日 收 到
。
刀 刀
·
彩
一
aN (
t
)
+ 一
V
口
N (
) = S ( V )y
质 中 的 中子 输 运
在 随 流 体 运动 的 坐 标 系 中
,
对 时 间的 导 数 由
.
~ ; 二
D _
,
刀右
一
一了 , -
d多
十 肠
份
一
’
V
一
吮
口 , 一
二
.
苏 口
十 肠 仁, 一
’
a
d
r z
( 12 )
给出
中子输运方程
dx dφ
φ dx
φ (x) = Ae−x / L
(指在B介质内的通量)
2D φ
dφ dx
=
2D Ae−x /
L
Ae−x / L (− 1 ) = − 2D
L
L
1− 2D
∴
β
=
1+
L 2D
L
(β ∞)
9
②当B厚度为a时(含外推距离)
φ (x) = Ash x + Bch x
L
L
用外推边界条件:0 = Ash a + Bch a
=
−
λs 3
∇φ。.
n = cosαi + cos β j + cosγ k
J = J ⋅n
J = −D∇φ
→ 斐克定律(Fick' s law), D为扩散系数。
用Σtr 代替Σ s,Σtr = Σt − Σ s µ0 ≈ Σ s (1 − µ0 )
并用D
=
λtr 3
.
§2.1.3 斐克定律的适用范围
x
⋅ −
1 L
sh a L
L
L
= − 2D cth a LL
( ) ∴ β
=
1−
2D
L
cth
a L
( ) 1+
2D
L
cth
a L
当a
→
∞时,即cth
a L
→ 1,β
→
β∞
讨论:反射层厚度对反照率的影响。
§2.5 扩散方程的积分形式
1.源的迭加性。扩散方程是线性的,中子的粒子性。 2.源互不干涉。因为中子-中子碰撞极少。
211216365_深度学习方法求解中子输运方程的微分变阶理论
复杂偏微分方程的方法正成为数值计算领域前 沿研究热点 & -=. 不同于传统的数据驱动机 器 学 习 方 法 !基 于 物 理 信 息 的 神 经 网 络 "]M''$等 深 度学习框架 不 -!. 需要预先收集大量数据作 为 机 器学 习 的 训 练 集!可 直 接 正 向 求 解 复 杂 高 阶 多 维 方 程 !并 具 有 空 间 几 何 限 制 少 '边 界 条 件 处 理 灵 活 '计 算 流 程 规 范 性 好 '实 验 数 据 同 化 能 力 强 等特点!已在求解 热 工 水 力 '@6 方 程 等 多 个 领 域取得了重 要 突 破 & -!@A. 在 求 解 反 应 堆 中 子 扩 散方程 领 域 也 取 得 了 系 列 进 展 !这 -"@<). 为 中 子 输运方程的求解开辟了一条新的潜在技术途 径 !越 来 越 受 到 业 界 的 重 视 &
<!,-./#%&0*P#&6)3'4'-'#*+514)%51#)9=-A<)(<*514)%(!5,,58)94)##&4)9B#3#%&.15#)'#& *+ ,-./#%&8)#&92A*+'P%&#%)= "494'%/B#%.'*&51#)9=-A<)(<*514)% *!5*//#9#*+5*($-'#&A.4#).#A4.1-%)?)4@#&34'251#)9=-A<))A!514)%
核反应堆中中子输运特性研究
核反应堆中中子输运特性研究在当今的能源领域,核反应堆作为一种重要的能源来源,其内部的物理过程复杂而关键。
其中,中子输运特性的研究对于核反应堆的设计、运行和安全评估具有极其重要的意义。
要理解核反应堆中中子输运特性,首先需要明白中子在反应堆中的作用。
中子是引发核裂变反应的“导火索”,它们与核燃料相互作用,导致原子核分裂并释放出大量的能量。
而中子的输运过程,就是中子在反应堆内的运动、散射、吸收和产生的综合过程。
中子的速度是一个关键因素。
快中子具有较高的能量,它们在与原子核碰撞时,可能引发裂变反应,也可能仅仅是散射而改变方向和降低能量。
慢中子则相对更容易引发裂变,因为它们与原子核的相互作用概率更高。
在核反应堆中,燃料的布置和结构会显著影响中子的输运。
例如,在压水堆中,燃料通常以棒状排列,周围充满了冷却剂。
这种布置方式既有利于热量的传递,也对中子的输运产生了特定的影响。
中子在穿过燃料棒时,可能被吸收引发裂变,也可能从燃料棒散射到冷却剂中,然后再回到燃料棒继续参与反应。
材料的性质对于中子输运也至关重要。
反应堆中的结构材料,如压力容器、控制棒等,会吸收一定数量的中子。
控制棒通常由能够强烈吸收中子的材料制成,通过插入或抽出反应堆堆芯,来控制中子的数量,从而调节反应堆的功率水平。
中子的散射过程是其输运特性中的一个重要环节。
当中子与原子核碰撞时,它们会改变方向和能量。
散射分为弹性散射和非弹性散射。
弹性散射中,中子的总能量保持不变,只是方向和动能发生变化;非弹性散射则会导致中子损失一部分能量。
中子的吸收现象同样不容忽视。
除了在燃料中引发裂变反应被吸收外,中子还可能被其他材料吸收而无法继续参与链式反应。
这就需要在反应堆设计中,精确计算中子的吸收概率,以确保有足够的中子维持链式反应的进行。
为了研究中子输运特性,科学家们采用了多种方法和技术。
理论模型的建立是其中的重要手段之一。
通过数学方程来描述中子的输运过程,可以对反应堆的性能进行预测和分析。
中子动力学
中子动力学中子动力学是一门研究中子在物质中的输运、相互作用和反应过程的学科。
这门学科涉及到很多重要的理论和实验技术,如中子产生、中子散射和中子吸收等。
中子在物质中的输运是中子动力学的一个核心问题。
中子散射、中子重复散射和中子吸收是中子在物质中输运的主要过程。
中子散射是指中子在物质中与原子核碰撞,改变其运动方向和能量的过程。
中子重复散射则是指中子在物质中进行多次散射。
而中子吸收是指中子与物质中的原子核发生核反应,使中子消失掉的过程。
中子在物质中的输运不仅决定着中子的传输性质,也对材料性质产生了很大的影响。
例如,中子的慢化过程可以增加中子与物质中原子核反应的截面,从而提高反应截面和中子吸收截面。
这对于核反应堆和其他中子辐射源的设计和运行具有重要意义。
在中子动力学研究中,还有一个非常重要的问题就是中子的相互作用。
中子与物质中的原子核、电子和其他粒子之间的相互作用是中子在物质中输运和反应的基础。
中子与物质相互作用强度和机制的了解,对于解决中子输运和反应问题非常重要。
中子的相互作用还涉及到材料中的缺陷、杂质、结构等问题,这些都可以通过中子散射和中子吸收实验来研究。
中子动力学的研究还与核能、能源、材料、医学、环保等领域密切相关。
例如,核反应堆的设计和优化需要中子动力学的研究,新能源材料的研究也需要中子散射和中子吸收技术的支持。
医学领域的放射治疗等技术也是基于中子相互作用和中子传输性质而发展起来的。
总之,中子动力学的研究对于推进科学技术进步,促进人类社会的发展具有积极的意义。
通过加强对中子相关技术和理论研究的支持和发展,我们可以更好地掌握中子的特性和行为规律,为解决实际问题提供更可靠的理论和技术基础。
中子输运方程
2
z
dV ྲᒦᔇ
r θ
dΩ
dA
y
ϕ x
ᅄ3.12!āᅎࡴᲝయࢾേࡼာፀᅄ
(点击图片可放大显示)
公式推导:
x − y平面上的dA
Σ sφ (r)dV
⋅
dA cosθ 4πr 2
=
S(r,t)
+
D∇2φ (r,t)
−
Σ aφ (r, t )
→ 扩散方程
S(r) + D∇2φ (r) − Σaφ (r)=0 → 稳态
边界条件: ①φ正值,有限 ②J ,φ连续 ③非凹外边界,J − = 0
线性外推距离(linear extrapolation distance):在单群中子输运理论中,渐进中子通量 密度在边界上的切线延伸到在介质外达到零的一点到介质边界的距离。
围立体角d Ω内的中子期望数]
2.中子角通量密度 ( Angular neutron density) ⇒ 中子通量密度 (标量) φ (r, E, Ω,t) ≡ vn(r, E, Ω,t)
3.中子流角密度( Angular neutron current): j(r, E, Ω,t) = vΩn(r, E, Ω,t) = Ωφ(r, E, Ω,t) 而Ω为单位矢量,故 j = φ
1.课本 72 页, 例 1:无限介质,点源,S 个中子/sec,各向同性。
边界条件:除源r条=件0外:lri→,m0 4φπ(rr
)为有限值。 2 J (r) = S
得解:φ (r)
=
Se −r / L 4πDr
中子输运方程
1建立单个中子的真实运动历史 2对大量中子历史进行跟踪,得到充足随机试验值 3用统计方法估计某个随机变量的数值特征
中子输运方程的积分形式
裂变源或外源
散射项
以上方程建立在特征线上,称为特征线线方程。
则发射密度方程可写成
(8-17)
构造概率分布密度f(xm),随机变量h(xm)的函数 则数学期望为:
计算期望
该模型是发射密度的无偏估计
中子发射密度的三种估计方法:
直接俘获 吸收估计 隐俘获
2
碰撞估计
3
吸收估计方法:
直接俘获估计缺陷:
隐俘获:
隐俘获特点:
举例说明:
碰撞估计:
碰撞估计特点:
中子输运方程源迭代递推
中子输运方程源迭代递推中子输运方程是核反应堆中描述中子输运行为的一种数学模型。
该方程通常采用著名的Boltzmann输运方程来描述中子的输运过程。
求解中子输运方程的一种有效数值方法是使用源迭代(Source Iteration)或者是Krylov子空间迭代方法,其中最著名的是功率迭代(Power Iteration)。
下面是源迭代方法的一般步骤:
1. 初始化:给定初始中子源,通常使用上一个迭代步骤的核反应率来计算。
2. 传输:使用Boltzmann输运方程进行中子传输的数值求解。
这一步通常使用离散的空间和角度网格,将方程离散化为代数方程组。
3. 修正源:根据传输的结果修正中子源。
这个修正通常是通过计算核反应率和散射源来实现的。
4. 迭代:重复步骤2和步骤3,直到解收敛为止。
通常,可以使用一些收敛准则,例如相邻两步之间的核反应率差异小于某个阈值。
5. 计算功率:当迭代收敛后,可以通过整合核反应率和中子通量来计算反应堆的功率。
这只是源迭代方法的一种简单描述,实际应用中会涉及到更多的数值技术和物理参数的处理。
在实际中子输运方程的求解中,使用了一系列的数值方法和近似技术,以便更好地适应不同的反应堆物理学问题。
二维柱几何中子输运方程的并行区域分解方法
二维柱几何中子输运方程的并行区域分解方
法
二维柱几何中子输运方程的并行区域分解方法是指将二维柱几何
中的空间分解成多个独立的(不相连)区域,以实现对物理数值解决
问题的有效分担。
这种分解方法主要应用于解决大型复杂的中子传输
问题,它具有很好的数值性能和易用性,能够很好地支持向量机和其
它并行硬件。
首先,用户需要根据不同的中子输运问题定义几何模型,设定有
限元尺寸和分边界条件,考虑空间的并行化,然后用一种合适的算法
将此几何模型分解成若干不相交的片,这些片充当空间上的并行区域,然后在每个区域内求解中子输运方程的离散方程,最后在分边界处进
行边界交换,完成最终的计算结果。
在实际应用中,二维柱几何中子输运方程的并行区域分解方法有
利于以更具体的空间结构进行代码优化,进而提高计算效率。
它可以
使用并行计算来改善计算性能,同时利用现有计算资源,为实际工程
中的中子热流量计算带来必要的精度和可靠性。
通过上述介绍,可以看出,二维柱几何中子输运方程的并行区域
分解方法具有很好的可扩展性、可靠性和灵活性,可以极大地提高计
算效率和精度,能够很好地支持复杂多变的物理研究。
它不仅可以应
用于二维柱几何中子输运方程,还可以用于其它类型的几何模型,以
解决更加复杂和深刻的数值计算问题。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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种核反应的出射粒子都是带电粒子,通常称为中子引起的带电粒子核反应。总的 来说,中重核的带电粒子反应截面都是比较小的。它们在反应堆临界计算中通常被忽略掉。
一、 一般分析
中子输运理论的基础假设是把中子看做一个点粒子,因此可以用其所在位置、具有的动能、 飞行方向进行精确描述。该假设的合理性在于中子的约化波长比介质的宏观尺寸和中子的平均自 由程要小得多。从德布罗意方程出发可以得到证明。 中子的输运理论是描述大量中子与介质相互作用的规律,它是中子与核发生各种碰撞过程的 综合描述。这是中子输运理论的基础。输运理论通常也称作迁移理论。求解某系统的中子输运方 程,可以得到该介质内中子随能量、空间、角度变化的规律。要进行这样的研究或计算,必须已 知单个中子与多个核素发生单个碰撞过程的大量信息。换句话说,大量的微观核反应数据是核反 应堆理论计算的数据基础。当核反应堆的体积较大时,从任何发生裂变反应产生中子的地方到系 统外边界的距离比介质的中子平均自由程大得多时,用扩散理论也能得到较好的结果,例如核电 站等大型功率堆得物理计算多采用扩散理论。