双热中子束流治疗孔道低浓化BNCT堆初步设计

以我国已经建成的高浓铀为燃料的BNCT堆为研究对象,将其堆芯低浓化并且添加水平热中子双束流治疗孔道,开展双热中子束流BNCT堆堆芯低浓化初步设计,计算分析该BNCT堆的keff、控制棒价值、顶铍效率、堆芯能谱、堆芯径向通量、轴向通量、辐照管通量等参数,得到双热中子束流治疗孔道低浓化BNCT堆初步设计方案.     

改进Flower型超临界水冷快堆负空泡反应性初步分析设计

超临界水冷快堆结构简化,堆芯功率密度高并且不存在沸腾危机,是一种比较有前途的先进核能系统.但潜在的正空泡反应性系数是需要认真考虑的安全问题.本文建立了改进Flower型超临界水冷快堆堆芯模型,并通过对堆芯分区布置,冷却剂密度分层、MOX燃料合理设计、燃料富集度轴向分层布置,blanket组件盒厚度增加及内部通道采用贫铀冷却的方案,获得了负的空泡反应性系数,初步达到了堆芯的安全设计要求.     

三束流驱动ADS次临界反应堆堆芯初步设计

在 "启明星-1"实验装置的设计基础上,设计了三束流驱动ADS次临界反应堆堆芯的布置方案,计算并比较了三束流驱动和单束流驱动ADS次临界堆内径向与轴向中子通量密度分布.结果表明,三束流驱动ADS的次临界反应堆的布置方案通过调整源点位置能降低堆内轴向功率形状因子,能够有效的展平堆内功率径向分布.     

医院中子照射器堆芯中子学参数计算方法

 建立了利用WIMS+CITATION计算医院中子照射器I型堆(IHNI-I)堆芯中子学参数的模型。栅元群常数计算采用WIMS束棒几何模型,控制棒、顶铍反射层、底铍反射层、侧铍反射层以及堆芯每一环燃料元件作为不同栅元类型;全堆芯计算采用CITATION程序R-z几何模型。计算了堆芯的功率分布、顶铍反应性价值、控制棒价值、温度系数、堆芯燃耗等中子学参数,计算结果与文献数据一致。本文所建立的计算模型可用于IHNI-I堆芯的物理计算。     

双排棒组件超临界水堆堆芯方案设计

结合国际上多种超临界水堆堆芯设计方案的优点,提出了一种新的压力容器式低泄漏堆芯设计方案,其特点是,堆芯中采用了双排棒正方形闭式燃料组件和三区低泄漏换料.双排棒燃料组件由两排燃料棒包围一个慢化剂水棒构成,可以使得慢化均匀;三区低泄漏换料可以大大延长堆芯寿期,降低压力容器快中子注量.通过堆芯三维物理热工耦合计算发现,该方案寿期内的最大包壳温度(MCST)为684℃,堆芯寿期为300个有效满功率天,且功率分布平坦.在此基础上,对所有组件进行了更为保守的子通道热工水力计算,得出MCST为685.3℃,进一步表明所提堆芯设计方案在物理热工方面是可行的.     

轴对称通道内超临界水对流换热数值模拟

为了研究高性能轻水堆堆芯冷却通道内超临界工质流动与换热特性,采用Fluent软件对竖直上升环形通道内超临界水与管壁间的换热过程进行数值模拟研究。通过计算结果与实验结果的对比,评估了不同网格数量、壁面第1层网格参考高度和湍流模型对数值模拟准确性的影响。对5个不同湍流模型进行对比计算。针对3种热流密度进行堆芯流动换热计算。模拟结果表明:随着网格数量的增加和壁面第1层网格参考高度的降低,预测准确度逐渐提高,但当壁面第1层网格参考高度小于0.07后,预测准确度与其无关;与其他湍流模型相比较,标准k-ω模型和剪切压力传输(SST)k-ω模型都能预测传热恶化现象,SST k-ω模型较标准k-ω模型计算结果更为准确;近壁面区域物性的剧烈变化是引起传热强化及传热恶化的主要原因。     

铀氢锆堆群常数库的生成及堆芯的物理和安全参数的计算

应用栅元计算程序ＷＩＭＳ－Ｄ／４及含有氢化锆中氢核素数据的ＷＩＭＳ－Ｎ２核数据库,计算了我国第１座商用铀氢锆堆（西安脉冲堆）的栅元参数,生成了两群群常数随燃耗变化的数据库．用所产生的群常数库和ＴＲＩＧＡＣ程序对西安脉冲堆堆芯重要的物理和安全参数———堆芯通量分布、Ｋｅｆ随燃耗的变化、负温度系数以及功率峰因子等进行了计算,这些参数对堆芯燃料管理和核安全分析具有重要的意义．     

高温堆双区堆芯交混区内球分布实验研究

双区堆芯的设计代表了高温气冷堆技术的发展方向之一。为确定高温气冷堆双区堆芯交混区内燃料球分布情况开展实验研究,解决与高温堆堆芯结构和物理设计密切相关的问题。实验结果表明:交混区内燃料球和石墨球的概率分布具有良好的对称性;堆芯上部的交混区内燃料球分散度较小,中部较大,交混区并没有无限扩散而具有确定的范围,说明球流运动是一种有限的分散性;在加球面上加入任一球其在交混区内的概率密度分布满足Gauss分布规律;根据Gauss分布对称中心可以确定交混区内球流均值流线并从理论上确定交混区大小。     

控制棒位置及局部温度扰动对熔盐堆有效倍增因数影响的研究

熔盐堆是第四代反应堆中唯一一种以液体为燃料的反应堆,因此对于熔盐堆反应性的研究不同于其他反应堆。基于蒙特卡洛方法,利用 Monte Carlo N-Particle(MCNP)软件模拟控制棒在堆芯径向不同位置及轴向不同插入深度对熔盐堆堆芯有效倍增因数的影响。随后将熔盐堆堆芯由上到下划分成八个区域,分别计算熔盐与石墨在八个不同区域发生多普勒效应以及发生膨胀效应对整个堆芯的有效倍增因数的影响。结果表明控制棒位于熔盐堆不同位置对反应堆有效倍增因数影响不同,沿径向21.21 cm处插入深度80 cm时控制棒有效利用价值最高。熔盐在不同区域发生多普勒效应时,顶部和底部对反应堆有效倍增因数影响相对较大。不同区域熔盐发生膨胀效应时,轴向中心处对有效倍增因数的影响相对较大。石墨发生局部温度扰动对有效倍增因数的影响较小。     

控制棒位置及局部温度扰动对熔盐堆有效倍增因数的影响

熔盐堆是第四代反应堆中唯一一种以液体为燃料的反应堆,因此对于熔盐堆反应性的研究不同于其他反应堆。基于蒙特卡洛方法,利用Monte Carlo N-Particle(MCNP)软件模拟控制棒在堆芯径向不同位置及轴向不同插入深度对熔盐堆堆芯有效倍增因数的影响。随后将熔盐堆堆芯由上到下划分成八个区域,分别计算熔盐与石墨在八个不同区域发生多普勒效应,以及发生膨胀效应对整个堆芯的有效倍增因数的影响。结果表明控制棒位于熔盐堆不同位置对反应堆有效倍增因数影响不同,沿径向21.21 cm处插入深度80 cm时控制棒有效利用价值最高。熔盐在不同区域发生多普勒效应时,顶部和底部对反应堆有效倍增因数影响相对较大。不同区域熔盐发生膨胀效应时,轴向中心处对有效倍增因数的影响相对较大。石墨发生局部温度扰动对有效倍增因数的影响较小。     

超临界水氧化水膜反应器的气固两相流分析

基于欧拉-拉格朗日方程的拉格朗日离散相模型,并建立颗粒的受力模型,对超临界水膜反应器内的气固两相流特性进行数值模拟。模拟并分析了颗粒直径、入口速度对超临界水膜反应器内的颗粒运动规律、分布特性以及颗粒接近多孔璧的参数的影响,并且分析了超临界水膜反应器中一定直径的颗粒的流体随体性。对研究和优化超临界水膜反应器的气固两相流特性具有一定的指导意义。     

中子通量测量系统可靠性

核电堆芯中子通量测量系统是对核反应堆进行中子通量实时精准测量的专用工具,其可靠性对核反应堆的运行也有重要影响.针对中子通量测量系统多状态的特点,提出一种将Markov模型和GO法结合的方法.首先通过Markov状态空间法对中子通量测量系统的可靠性进行分析建模,求得系统处于各个运行状态下的概率,然后根据系统原理构建GO图,用GO法实现了多状态下的中子通量测量系统可靠性分析.最后以某核电站中子通量测量系统为算例,验证了该方法,为中子通量测量系统可靠性的分析提供了有效的手段.     

HTGR控制棒区不连续因子受边界条件的影响

为改善高温气冷堆控制棒区物理计算精度,采用不连续因子理论修正扩散方程,使用扩散差分法程序C ITAT ION进行求解。对于1/60和1/30扇区堆芯模型,分析了边界条件和均匀化区对不连续因子计算的影响。结果表明,以离散纵标法程序SN 2D求得的精细网格的输运解为基准,在最大误差2%时,用该方法计算得到的有效增殖因子、反应性价值和活性区的中子注量率分布均与参考解相符。采用不同均匀化模型的不连续因子计算,受边界条件的影响都不大。     

一种基于MCNP的堆芯物理—热工耦合计算方法

核反应堆的温度分布计算在核反应堆的安全设计中十分重要.我们利用MCNP计算的线功率密度带入MATLAB编写的单通道模型热工程序,配合热工计算得到的冷却剂密度分布调整MCNP的模型重新计算,能够将堆芯物理计算和热工计算进行结合,得到了更为准确的单通道温度分布、功率分布及热管因子.     

压水堆燃料元件棒传热的分析

反应堆中堆芯传热问题非常复杂.本文对压水堆燃料元件棒内的传热做了简化,将燃料芯块视为含内热源的稳态传热,芯块与包壳之间的间隙以及包壳看作稳态导热,探讨燃料元件棒内热量的传递,得到各区间的温度分布、热流量分布.     

交混翼对方形棒束通道内超临界流动传热影响的初步研究

以CFD软件FLUENT为工具对超临界水堆的三维带交混翼的方形单通道模型进行了数值模拟.模型考虑了交混翼的仰角对于通道内温度分布及流场的影响,结果表明,交混翼造成流体绕流,使得流道内冷热流体混合,从而使得流道内流体的温度分布均匀,有效改善了燃料棒表面温度分布情况,降低热点温度;交混翼在15°、30°、45°和50°四种仰角情况下,30°仰角情况更有利于温度区域均匀.     

An accurate solution of point kinetics equations of one-group delayed neutrons and an extraneous neutron source for step reactivity insertion

The continuous indication of the neutron density and its rate of change are important for the safe startup and operation of reactors. The best way to achieve this is to obtain analytical solutions of the neutron kinetics equations because none of the developed numerical methods can well satisfy the need for real-time or even super-time computation for the safe startup and operation of reactors in practice. In this paper, an accurate analytical solution of point kinetics equations with one-group delayed neutrons and an extraneous neutron source is proposed to calculate the change in neutron density, where the whole process from the subcritical stage to critical and supercritical stages is considered for step reactivity insertions. The accurate analytical solution can also be used as a benchmark of all numerical methods employed to solve stiff neutron kinetics equations.