数字化反应堆高保真中子学程序SHARK研发

SHARK程序是由中国核动力研究设计院新近研发的基于全堆芯确定论非均匀输运理论体系的数字化反应堆软件。该软件从多群数据库的截面与共振数据出发,采用改进子群方法刻画有效共振截面的复杂非均匀效应,采用二维/一维或准三维特征线方法开展堆芯层面非均匀输运计算。目前该程序的定态微观问题计算能力已建立完毕。数值结果显示,SHARK程序对于商用压水堆相关基准问题具有良好的计算精度和效率。     

一步法输运计算程序KuaFu开发与验证

为了拓展一步法输运计算方法在结构复杂先进反应堆中的应用,基于构建实体几何理论及二维/一维耦合方法,采用C++、Python混合编程开发了一步法输运计算程序KuaFu,并应用粗网有限差分方法(CMFD)、大规模并行技术对二维/一维耦合方法进行加速。通过C5G7基准题对几何建模的可视化功能、并行功能及计算精度进行评估,获得计算结果与蒙特卡罗程序(MCNP)的相对误差。计算结果表明,程序具有较好的可视化功能和用户友好性;KuaFu程序与MCNP参考解符合较好,计算精度良好。     

基于DRAGON程序的二维/一维耦合堆芯输运计算程序开发

DRAGON程序能够进行三维特征线法(MOC)计算,但需大量的计算时间与内存,虽在其中引入大量近似并应用多种加速算法,仍无法满足工程应用的需求。基于DRAGON程序中现有的二维MOC计算模块,开发二维MOC与一维节块展开法(NEM)耦合的堆芯输运计算程序DRAGON_HEU,在三维粗网有限差分(CMFD)加速算法的全局框架下通过轴向和径向泄漏项将二维全堆非均匀栅元MOC计算与轴向均匀栅元扩散计算耦合,实现三维堆芯Pin-by-Pin计算。运用C5G7-3D三维扩展基准题验证DRAGON_HEU,计算结果表明:DRAGON_HEU能够节约大量计算时间且具有很高的计算精度。     

六角形布置反应堆高保真物理计算方法研究与验证

俄罗斯商用压水堆VVER和大多数实验堆均采用了六角形紧凑型栅格布置,为了实现VVER和六角形实验堆的高保真数值模拟分析,本文基于数值反应堆物理计算程序（NECP-X）开展了六角形堆芯高保真计算方法研究和程序开发。首先,将全局-局部耦合共振自屏计算方法拓展至六角形堆芯,实现六角形堆芯燃料棒的全堆芯高精度共振计算;其次,基于2D/1D耦合输运计算方法研究了六角形堆芯的高保真计算方法;最后,为了提高全堆芯计算的计算效率,研究了基于区域分解松耦合的非结构网格的粗网有限差分（CMFD）加速方法,可以实现以矩形、六角形和其他多边形栅元为基础的pin-by-pin CMFD 加速。为了验证六角形堆芯高保真计算方法的精度和效率,计算了六角形C5G7基准问题,并分析了六角形输运计算方法的计算精度和CMFD方法的加速效果;将NECP-X程序应用于西安脉冲堆的2D全堆芯计算,与蒙特卡罗程序的结果对比表明NECP-X程序计算得到的特征值和功率分布均具有较高精度。因此,本文建立的六角形堆芯高保真计算方法可以应用于六角形堆芯的分析计算。      

基于数字化技术的反应堆测量系统修正方法研究

在分析影响测量系统准确度的各种影响因素基础上,基于数字化技术的计算优势在反应堆系统参数测量中提出"仪表修正"和"系统修正"2种修正方法。2种修正方法用数字化技术实现后,提高了瞬态和事故工况下测量的准确度,扩展了仪表的应用范围,可进一步提高反应堆的安全性和经济性。     

数字化核测量系统反应堆倍周期计算的敏感性分析

反应堆倍周期是核反应堆工程中的一个重要参数。在反应堆启动和功率提升过程中,操纵员可通过反应堆倍周期来了解反应堆的运行状态,并据此控制反应性。数字化核测量系统通过对与反应堆功率成正比的电压信号进行采样和处理,计算得到反应堆倍周期。在实际的应用中,电压信号往往包含测量噪声,对计算结果带来较大的不确定性。针对数字化核测量系统的倍周期计算问题,对其敏感性进行了分析,并给出相应的算例。     

NECP-X中基于改进泄漏项分割技术的二维/一维耦合输运方法

二维/一维耦合输运方法较好地平衡了效率与精度,因此被广泛应用于一步法全堆芯输运计算。二维/一维耦合输运方法中,由于泄漏项在方程右端,导致二维特征线法(MOC)计算时方程右端总源项在迭代过程中可能成为负值,造成迭代发散。本文针对二维/一维耦合输运计算中的负源项问题,提出了一种改进的泄漏项分割方法。新的泄漏项分割方法可在不造成计算精度损失和仅增加有限内存的条件下,显著提高二维/一维耦合输运方法的稳定性。通过强泄漏算例、C5G7基准题、VERA-3A基准题等进行测试,表明该方法对提高二维/一维耦合输运方法稳定性具有显著的效果。     

基于CAD技术的特征线中子输运计算程序开发

特征线方法(Method of Characteristics,MOC)能否应用于复杂几何关键在于能否将特征线方法与有效的几何处理方法结合起来。本文在菱形差分特征线理论基础上,基于FDS团队自主研发的核与辐射输运计算自动建模软件MCAM的几何处理引擎,研发了基于CAD技术的特征线中子输运计算程序,并利用相关基准例题对程序进行了数值验证,其结果与参考值吻合良好,表明本文方法和程序的可行性、正确性与可靠性。     

KYCORE程序中子输运计算的改进

KYCORE程序是在中国核动力研究设计院开发的二维组件计算程序KYLIN-2基础上开发的三维堆芯数值计算程序,其中中子输运部分采用径向特征线方法(MOC)与轴向离散坐标法(SN)直接角通量耦合的方法实现高精度计算,并通过粗网有限差分方法(CMFD)加速实现快速收敛。KYCORE程序因为计算流程的简化,导致可能出现不收敛,因此在计算方法和网格划分上做了改进,提高了计算的稳定性和包容性。通过与C5G7扩展基准题和蒙特卡罗程序的计算对比,数值验证了KYCORE输运部分计算的稳定性与准确性。     

基于组件模块化特征线方法的中子输运计算研究

栅元模块化特征线方法(MOC)在处理压水堆组件水隙等问题上存在几何处理上的困难。为了克服这些问题,采用Fortran语言开发了基于组件模块化特征线方法(AMRT)的中子输运计算程序NECP-Medlar,并采用两重粗网有限差分方法进行加速。2D C5G7基准题和典型压水堆组件问题的数值计算结果表明,该程序具有良好的计算精度和较高的计算效率,并且能够通过直接计算组件之间的水隙,较精确地描述组件中子通量密度的分布。     

基于TXS平台的数字化反应堆保护系统实现特性

数字化反应堆保护系统(Reactor Protection System,RPS)作为重大改进项在CPR1000压水堆机组岭澳核电站首次引入,设计方案经历了多次变更和优化,最终伴随着岭澳核电站成功商运后完成了方案固化。本文介绍了基于阿海珐公司提供的安全级数字化仪控平台(TELEPERM XS,TXS)在岭澳核电站反应堆保护系统最终实现方案,进而重点分析反应堆保护系统方案实现特性。上述论述对后续核电站反应堆保护系统方案实现及其国产化研制提供了借鉴和参考意义。     

数字化反应堆技术在设计阶段的应用研究

为了解决核反应堆设计效率低、资源集成度不高、创新能力受限等问题,提出了数字化反应堆技术研发的实施构想。通过对数字化反应堆技术在核反应堆研发设计阶段的应用研究情况,发现搭建灵活可扩展的数字化基础平台框架、建立基于系统工程的三维协同设计体系以及基于知识工程的大数据管理是数字化反应堆技术应用在反应堆工程设计阶段的重要研究内容。      

弥散颗粒型燃料特征线方法输运计算研究

弥散颗粒型燃料的中子输运问题因其特有的随机性和双重非均匀性难以直接使用现有输运方法进行求解。Sanchez-Pomraning方法借助更新方程,对特征线方法进行改进,使其能应用于弥散颗粒型燃料的输运计算中。本文对二维圆柱形弥散颗粒燃料输运问题进行了计算,数值结果表明：程序在不同颗粒填充率、不同颗粒尺寸、燃料颗粒与毒物颗粒共存的问题下均能保证较好的计算精度,反应性特征值绝对偏差大多低于100 pcm,仅在QUADRISO毒物颗粒填充时绝对偏差达到163 pcm。本文方法能满足弥散颗粒型燃料的输运求解要求,为新型燃料的设计研究工作提供了可靠的结果。     

基于离散纵标输运计算方法的三维燃耗程序发展研究

为了精确描述和分析具有强烈各向异性中子注量率空间分布的反应堆燃耗过程,本文实现了三维SN 输运计算与燃耗计算的耦合,发展了相应的三维输运燃耗耦合计算程序.该程序系统采用接口程序自动耦合三维SN输运计算程序和同位素燃耗计算程序的方法实现对三维中子学计算模型的精细燃耗计算,获得燃料同位素成分、燃耗反应性、中子注量率空间分布等参数随燃耗时间的变化量.采用IAEA 基准校核例题对程序系统进行了校核,计算结果初步证明了所开发的三维燃耗程序系统的正确性.     

基于CSG和OpenMP的复杂几何输运计算程序开发及验证

为了适应材料几何布置越来越复杂的小型研究堆计算需求,基于构建实体几何理论和矩阵特征线方法,开发了具有复杂几何输运计算能力的2维特征线程序MOCAGE,并采用OpenMP并行编程模型对几何前处理中的特征线追踪进行并行化设计。通过不规则几何问题以及3种不同控制棒布置形式的HTTR基准题对程序的特征线追踪能力与计算精度进行评估,给出了计算结果与MCNP5多群计算参考值的相对误差。结果表明:所开发的程序能够正确实现对复杂对象的几何建模并进行特征线追踪,计算结果与参考值符合较好,精度满足程序验证要求,采用OpenMP并行编程能显著减少几何预处理时间。     

一种求解反应堆中子扩散方程和输运方程的方法

结合反应堆中子物理方程的特点和非结构网格技术,提出了一种在任意网格下求解反应堆中子物理方程的方法。以反应堆中子扩散方程和一阶离散纵标的输运方程为例,从空间离散、方程离散和边界条件实施等方面介绍了该方法的实施过程。利用该方法编制的程序计算了BN-600基准题、CFR1000概念堆和BWR栅元3种情况的有效增殖因数,并与其他程序的计算结果进行了对比,初步验证了该方法的合理性。     

SONG—多功能栅格物理程序概述

面向下一代反应堆堆芯分析的多功能栅格物理程序(SONG)采用特征线方法(MOC)求解任意几何区域内的中子输运方程,可灵活模拟由棒状元件或板状元件组成的矩形栅格或六角栅格的两维燃料组件模型或两维堆芯模型。程序采用指数矩阵方法求解拓展的裂变产物链和全锕系嬗变链,以满足超长寿期、超深燃耗的钍铀燃料循环或铀钚燃料循环计算的功能需求。相关验证计算结果表明SONG程序满足预期要求。     

基于特征线方法的三维中子输运程序(Ⅱ)——数值验证

介绍了开发的三维中子输运特征线方法程序TCM,使用数值实例对TCM的准确性进行了验证,同时比较了不同插值阶数、离散角度数目和网格大小对计算结果的影响。结果表明,三维特征线方法程序TCM能够求解任意几何形状及反射边界条件的中子输运问题,计算结果具有较高的精度。     

基于特征线方法的三维中子输运程序(Ⅰ)——边界条件的插值处理

特征线方法在处理三维中子输运问题的反射边界条件时存在空间问题和角度问题,本文提出了采用平面插值和球面插值方法来分别处理这两个问题,插值方法不仅可以有效解决空间问题和角度问题,而且三维特征线方法在采用插值方法处理边界条件后对离散求积组和射线布置没有要求。