反应堆的自然循环能力是其固有安全性的一项重要特征参数，利用RELAP5/Mod 3.4程序计算了JRR-3M池式研究堆在无应急冷却系统和有应急冷却系统条件下失去场外电源的事故工况，分析了冷却剂流动方向反转过程中的瞬态自然循环能力，并得到了最大自然循环载热能力。计算结果表明：应急冷却系统的投入明显地降低了燃料和冷却剂的温度，提高了反应堆的安全性；当衰变功率降低至590 kW时关闭辅助泵，利用自然循环能力可使反应堆达到安全状态。通过此研究堆的自然循环能力计算，验证了计算模型、计算方法的可行性，可进一步应用于此类型的研究堆。

Geant4是一款基于C++面向对象技术的蒙特卡罗开发程序包,可以模拟各种已知粒子与物质之间的相互作用。然而该程序包没有提供临界源功能,无法直接用于反应堆物理计算。因此,利用Geant4提供的基础物理模型和粒子跟踪控制等功能,用两种不同方法实现了临界源的设置,实现了基于Geant4的反应堆静态计算程序G4-RSM和反应堆动态计算程序G4-RDM。两个程序均可用于反应堆临界计算,与MCNP计算结果相对误差在5%以内。G4-RDM程序除可用于临界计算外,还可用于模拟堆内事故工况下的中子学瞬态变化。

NECP-SARAX是西安交通大学NECP团队开发的用于快中子反应堆的中子学程序系统。为准确处理快中子反应堆中中等质量核素散射共振以及非弹性散射导致的复杂的中子慢化效应,SARAX程序最初采用连续能量的蒙特卡罗方法计算中子能谱从而获得堆芯计算使用的有效多群截面。由于蒙特卡罗程序计算效率低,且在低能量段统计偏差较大,提出采用基于点截面的超细群方法计算中子慢化能谱,避免了蒙特卡罗方法产生参数时存在的缺陷。堆芯计算采用多群中子输运,通过优化简化几何建模,改进了程序的实用性。采用多种微扰方法计算堆芯各种反应性系数,提出了基于中子输运微扰理论的虚拟密度方法以计算堆内组件变形导致的反应性变化。在进行堆芯瞬态计算时,采用了点堆和改进准静态两种方法,可用于一般快堆和快谱ADS的典型事故分析。OECD发布的一系列快堆基准题测试表明，SARAX程序在快堆计算中具有良好的精度,达到了与国外著名快堆程序相当的水平。有效增殖因子与连续能量的蒙卡计算结果相比偏差在300 pcm以内。同时,由于引入了虚拟密度理论和三维时空动力学模型,程序功能更加完善,可以更好地满足快堆工程设计的需求。

采用目前工业上成熟的非线性迭代计算策略,基于两群粗网有限差分方法和多群UNM 节块方法,开发了针对压水堆工况的三维瞬态扩散计算程序。UNM 方法采用解析基函数作为基函数,通过方程变换解决了解析节块法在临界节块计算不稳定的问题,提高了计算精度。热工计算采用单通道模型和燃料棒一维导热模型,相比目前堆芯在用的经验关系式方法,该模型可以更加准确地计算燃料棒温度分布。采用基于横向积分方程的三节块方法,可以有效减轻控制棒尖齿效应对瞬态计算的影响。为测试程序性能,采用NEACRP等基准算例对程序进行了校验。数值结果表明,开发的程序计算结果正确,适用于压水堆堆芯瞬态过程的模拟。