基于中子电报方程的均匀裸堆瞬态中子输运问题解析解研究

斐克定律作为扩散理论的基础,其本质是在中子输运方程的P₁近似方程(电报方程)的基础上进一步忽略了中子流密度关于时间的导数项,因此在瞬态过程中无法准确描述实际的中子动力学过程。本文基于单能中子电报动力学方程,针对一维无限大平板裸堆中子动力学问题,利用分离变量法推导了其解析解,并与扩散中子动力学方程的解析解进行了对比分析。研究发现：在瞬态变化过程中,电报方程的解与扩散方程的解相比,空间项仍保持余弦函数的形式,但是时间项的变化则更为复杂,一是时间项的级数组合形式受到问题的几何和材料的影响,二是高阶谐波表现出随时间振荡变化的规律。该研究结果可以为后续基于中子电报方程的数值理论研究提供参照和依据。     

点堆随机动力学方程伊藤解与中子数概率分布解析解的对比分析

为了研究反应堆弱中子源启动过程中的中子数密度和缓发中子先驱核随机涨落现象,我们推导和建立了点堆随机动力学方程组,在传统的点堆动力学方程组中引入了伊藤随机项。为了验证方程组的伊藤解方法和计算精度,我们在简化物理条件和方程形式下,对定态系统的中子数密度分别用随机动力学方程伊藤解和中子数概率分布函数解析解进行了对比分析。结果表明,伊藤解是一种有效、具有较高计算精度的方法,计算精度满足sigma的标准,置信水平在95%以上。     

新概念熔盐堆物理计算方法研究及程序设计

考虑新概念熔盐堆燃料盐的流动特性,从基本的粒子守恒方程出发,推导了熔盐堆的中子动力学模型,并采用数值方法对3种工况下熔盐堆的临界问题进行计算,考察流动对有效增殖系数、快中子分布、热中子分布及缓发中子先驱核分布的影响。结果表明：质量流量对有效增殖系数的影响很小,对热中子分布的影响比对快中子分布的影响大,而质量流量越大,缓发中子先驱核移出堆芯的比率也越大。     

反应堆时空动力学方程近似解

导出了反应堆时空动力学方程的近似解析解.此解的中子密度空间分布曲率是随堆功率变化的,因此它比点堆模型方程的解析解精确.可用于反应堆现场运行所需的监督性快速计算.对反应堆的安全运行有重要的指导意义.     

熔盐堆物理热工耦合程序开发及验证分析

熔盐堆(Molten Salt Reactor,MSR)是第四代反应堆6种堆型中唯一的液态燃料反应堆,与固态燃料-液体冷却剂反应堆相比,原理上有较大不同。在熔盐堆中,流动的熔盐既是燃料又是冷却剂与慢化剂,中子物理学与热工水力学相互耦合;由于熔盐的流动性,缓发中子先驱核会随燃料流至堆芯外衰变,造成缓发中子的丢失,导致堆芯反应性降低。正是由于熔盐堆的这些新特性,造成熔盐堆内缓发中子先驱核、温度等参数变化与固态燃料反应堆有所不同,需要研究熔盐堆在各种工况下的相关物理参数变化。本文主要工作是考虑缓发中子先驱核的流动性对熔盐堆的影响,研究适用于熔盐堆的二维圆柱几何时空中子动力学程序及与之耦合的热工水力学程序;利用该程序对熔盐堆中子物理学和热工水力学进行耦合计算,验证熔盐堆相关实验数据;并且计算了熔盐堆无保护启停泵及堆芯入口温度过冷过热工况,用于分析熔盐堆的安全特性。计算结果表明,程序能够对熔盐反应堆实验(Molten Salt Reactor Experiment,MSRE)的相关实验数据进行较好的模拟计算,并且验证了熔盐堆的固有安全性。     

改进准静态本征函数法求解反应堆时空中子动力学方程

本文提出一种改进准静态本征函数法并用于求解时空中子动力学方程。在改进准静态近似条件下,利用因子分解法将中子通量密度分解为幅函数和形状函数,进而采用本征函数法求解形状函数轴向和径向的近似解析解。在求解幅函数过程中,将幅函数化简为点堆模型下的形式进行求解,最终导出反应堆三维时空中子动力学方程的解析解。与其他解析法和数值解相比,改进准静态本征函数法的适用范围更广,计算速度更快。     

液体燃料熔盐堆三维稳态分析程序开发

液体燃料熔盐堆的物理热工特性与固体燃料反应堆有很大的不同,在分析计算中必须考虑燃料流动特性的影响,一般分析固体反应堆的程序均不能直接用于分析液体燃料熔盐堆。根据熔盐堆的流动特性,建立了液体燃料熔盐堆的三维中子动力学模型和流动传热模型,开发了针对液体燃料熔盐堆的三维稳态核热耦合程序,并以此分析了稳态情况下MOSART堆的物理热工特性。结果表明,堆芯流速对快中子和热中子影响较小,对堆芯温度和缓发中子分布影响较大。     

熔盐堆中燃料流动对缓发中子的影响分析

熔盐堆具有良好的中子经济性、固有安全性、可在线后处理、防核不扩散等特点,是六种第四代先进反应堆堆型中唯一的液体燃料反应堆。然而,熔盐堆中采用流动的熔盐作为液体燃料,从而缓发中子先驱核会随着燃料的流动流出堆芯并在堆芯外发生衰变,这不同于固体燃料反应堆。文中针对了一座实际运行过的熔盐实验堆(Molten Salt Reactor Experiment,MSRE),基于中子动力学模型,采用圆柱体均匀堆的近似处理方法推导了液体燃料反应堆的缓发中子先驱核浓度数学模型,研究了恒定流速下的反应性损失及不同燃料熔盐流速对缓发中子分布的影响。结果表明缓发中子在越靠近堆芯中心区域的位置就越多,同时熔盐流速的变化对衰变周期越短的缓发中子先驱核组数的影响比较小。通过本研究,可以了解熔盐堆中缓发中子随着燃料流动的变化情况,为熔盐堆安全分析提供参考依据。     

液态燃料熔盐堆三维动力学程序开发及验证

液态燃料熔盐堆的燃料熔盐在一回路中循环流动,一回路高温熔盐既是燃料,又是冷却剂,大部分核裂变能直接释放在燃料熔盐之中。随着燃料熔盐流动,一部分缓发中子先驱核（Delayed Neutron Precursors,DNP）在堆芯外一回路中衰变引起反应性损失。液态燃料熔盐堆中子物理与热工流体紧密耦合,传统固态燃料反应堆堆芯核热耦合程序不再适用于液态燃料熔盐堆。针对液态燃料熔盐堆特点,建立了包含带对流项的DNP输运方程和带热内热源热工流体方程的液态燃料熔盐堆动力学模型,并基于节块展开法,开发了堆芯三维动力学程序ThorCORE3D。使用美国橡树岭国家实验室建造运行的熔盐实验堆（Molten Salt Reactor Experiment,MSRE）稳态和瞬态实验基准题,对ThorCORE3D程序进行了初步验证。结果表明：ThorCORE3D程序计算值与MSRE实验值吻合良好,适用于液态燃料熔盐堆稳态设计与瞬态分析。     

石墨慢化通道式熔盐堆有效缓发中子份额计算方法研究

有效缓发中子份额(βeff)是研究反应堆动力学特性的关键参数。在液态燃料熔盐堆(MSR)中,燃料流动引起缓发中子先驱核(DNP)在堆内的再分布,并使部分DNP在堆外回路衰变,从而导致βeff的计算方法与固态燃料反应堆不同。为评估石墨慢化通道式熔盐堆内燃料流动引起的反应性损失,研究缓发中子随燃料的流动行为,同时为堆设计和安全分析提供依据,分别基于解析方法和数值方法推导了计算βeff的数学模型,计算了熔盐实验堆(MSRE)在额定工况下的DNP损失份额和堆内DNP浓度分布,并分析了燃料在堆外流动时间和入口流量对βeff的影响。结果表明：两种方法均可对DNP行为提供合理描述;固定燃料在堆外流动时间,βeff随入口流量的增加而减小;固定入口流量,βeff随燃料在堆外流动时间的增加而减小,80 s后趋于稳定。