界面流法计算反应堆六角形燃料组件中子通量密度分布

利用界面流法计算两维六角形轻水堆燃料组件中子通量密度分布。子区内中子源在空间上采用二次分布近似，还考虑了六角形组件周边水隙对组件内中子通量密度的影响。根据提出的模型，编制了TPHEX－E程序，并对一些轻水堆六角形组件问题作了计算，计算结果与蒙特卡罗方法计算结果进行了比较，符合良好。本程序可用于六角形轻水堆燃料组件计算。     

六角形轻水堆组件中子通量密度分布的计算

介绍利用穿透概率法求解二维六角形轻水堆燃料组件中子通量密度分布。子区内中子源及通量密度在空间上采用二次分布 ,子区表面通量密度在空间上采用平通量密度分布 ,在方向上采用简化 6P1近似。根据提出的模型 ,编制了TPHEX D程序 ,并对一些轻水堆六角形组件问题作了计算 ,计算结果与MC结果进行了比较 ,符合良好。本程序可用于六角形轻水堆燃料组件计算。     

穿透概率法求解二维六角形轻水堆燃料组件中子通量密度分布

研究利用穿透概率法求解二维六角形轻水堆燃料组件内中子通量密度分布。子区内中子源采用线性分布,子区表面通量密度在方向上采用简化6P1近似。提出了六角形组件周边水隙的处理方法。根据提出的模型,编制了TPHEX－C程序,并对六角形组件进行了计算,结果与蒙特卡罗方法计算的结果符合良好。     

三维六角形组件压水堆堆芯燃料管理计算及程序系统研究

介绍所研制的WWER型压水堆堆芯燃料管理计算程序系统TPFAP-H／CSIM-H，六角形组件均匀化计算程序TPFAP-H是在压水堆正方形组件程序TPFAP的基础上，采用穿透概率法与响应矩阵方法相结合计算六角形组件内中子能谱分布，并考虑六角形栅元特点改造开发而成的CSIM-H是以先进六角形节块扩散程序为基础．参照SIMULATE程序功能而研制的物理-热工水力耦合的三维六角形节块PWR堆芯燃料管理程序两者通过接口程序LINK连接起来，可以考虑燃耗，功率、慢化剂密度变化．控制棒、氙等参数的多种反馈效应对IAEA的WWER-1000型Kalinin核电厂基准问题的校算的结果表明，临界硼浓度、功率和燃耗分布等结果与国际各研究机构的结果吻合良好，偏差均在工程要求之内。     

燃料组件少群参数计算程序KYLIN-1的基准验证

利用国际原子能机构(IAEA)压水堆棒状燃料组件基准问题和板型燃料组件基准问题对燃料组件少群参数计算程序KYLIN-1进行了验证分析。结果表明:KYLIN-1程序计算得到的燃料组件无限增殖系数(k∞)和重要核素核密度结果与国际上其他机构的计算结果符合良好;棒状燃料组件相对功率分布计算结果与参考程序符合较好。     

快堆燃料组件内精细功率分布的计算

堆芯热通道因子是堆芯热工设计及安全分析的一项重要参数,确定热通道因子需用中子学计算给出较准确的燃料组件内元件棒功率分布。在三维六角形几何节块扩散理论基础上,使用多项式重构的方法计算节块内中子通量密度分布和功率密度分布。针对快堆六角形燃料组件的特点,用小六角形积分的方法计算组件内元件棒功率,得到组件内各元件棒功率分布。在NAS程序基础上,编制了元件棒功率分布计算模块NAS PIN。通过与蒙特卡罗程序的校验可发现,二者计算结果符合较好,计算精度可满足工程设计的需要。     

SONG—输运计算模块程序开发

多功能栅格计算程序SONG采用特征线方法(MOC)及粗网有限差分(CMFD)加速方法进行中子输运计算,具备在数据库能群结构下全组件精细几何计算能力。同时具有与MOC相适应的几何预处理模块,采用基于组件的模块化射线追踪,可处理方形、六角形组件及棒状、板状燃料元件。通过模块化的流程与数据结构设计,开发形成了几何无关的MOC输运求解模块,同时形成了可扩展的组件几何预处理模块。不同形状组件的几何处理模块与输运求解模块具有统一的数据接口。通过相关问题的计算表明,SONG程序具备多几何组件处理能力,同时输运计算结果具有较好的精度、效率及稳定性。     

棒状氢化锆慢化钍基熔盐堆燃料组件稳态核热耦合程序开发

基于蒙特卡罗粒子输运程序MCNP与自主开发的子通道热工水力学程序SubTH,开发了棒状氢化锆慢化钍基熔盐堆燃料组件稳态核热耦合程序MCNP-SubTH,解决核热耦合程序因网格类型不同难以耦合的问题,程序具有普适性。MCNP-SubTH通过外耦合的方式进行MCNP和SubTH之间的数据交换,将MCNP计算得到的功率场加载到SubTH的求解文件中,然后将SubTH计算得到的密度和温度场更新到MCNP的输入卡中,实现程序迭代计算。分模块验证了MCNP-SubTH的准确性,并用MCNP-SubTH对棒状氢化锆慢化钍基熔盐堆燃料组件进行了稳态核热耦合计算,验证了核热耦合方法的有效性。     

压水堆六角形燃料组件均匀化 计算软件包TPFAP-HEX

介绍了所研制的具有工程实用价值的压水堆六角形燃料组件均匀化计算软件包。该组件中子空间能谱的计算采用穿透概率法与响应矩阵法相结合的方法，在六角形几何内求解中子积分输运方程。在此方法中，栅元内中子源采用空间线性或二次近似，栅元表面中子通量密度角分布采用简化6P     

超临界六角形双排燃料组件性能分析

利用物理-热工水力耦合计算程序系统(MCATHAS)分析2种六角形双排超临界燃料组件,充分考虑了超临界水冷堆(SCWR)中冷却剂、慢化剂轴向温度、密度的剧烈变化和功率分布的相互影响。计算结果表明,双排六角形组件具有均匀慢化和充分慢化性能,文中提出的D6-1型组件在仅采用一种燃料成分、不添加可燃毒物的情形下,其径向功率峰值因子低于1.10。另外,研究表明,由于组件间隙具有较大热周和较小流通面积,需要在实际工程应用中增加隔热涂层以降低组件外盒壁的导热率。     

基于二维特征线方法的模块化六角形几何预处理方法

本文采用模块化特征线跟踪的方法,借鉴矩形几何的处理方式,针对六角形几何开发了几何预处理算法,得到的几何信息为特征线计算提供前端输入。实际编程发现该方法耗时在ms级别,效率很高。同时扩展了对六角形阵列和组件的几何描述,针对该几何预处理方法,研究了对应的边界条件描述方法,为特征线计算程序的编写奠定了基础。     

界面流方法二维六角形组件计算

讨论了用界面流方法计算二维六角形组件中子通量分布。从积分输运方程出发,导出了一种简便的数学模型,在子区内采用平源通量近似,并假设中子发射和散射为各向同性。在子区表面上,中子通量的空间分布为常数,中子角通量分布通过伴随勒让特多项式展开表示,采用DP_1近似。推导出界面流方程组,给出了泄漏、穿透几率矩阵及其矩阵元素的表达式及计算方法。根据提出的数学模型,编制了TPHEX程序,对二维六角形组件进行了计算,本程序可用于水堆六角形燃料组件计算。     

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

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

燃料管理工作总结

2004是核能开发“十五”计划的第四个年度,现役核电厂燃料管理技术研究在本年度完成的主要工作如下。1)六角形组件燃料管理程序包,学习程序TPHEX、HBRP、ADCB、LINK、SIXTUS及SIMSIX。发现参数计算程序TPHEX在输入含有两类富集度燃料组件中燃料棒位置时,程序出错,经调试和研究,     

快堆组件内精细功率分布的改进方法与程序

在应用中子扩散方程六角形节块法程序时，往往需要计算组件内的精细功率分布。本文简要叙述了一种利用节块法程序的结果信息，再构造六角形组件内精细功率分布的计算方法，根据此方法编制了计算机程序ＨＥＸＰＩＮ。与细网有限差分程序的结果相比较表明：ＨＥＸＰＩＮ在堆芯功率峰值处的最大计算误差小于１％。     

基于变分节块法六角形Quasi-diffusion程序开发及验证

传统扩散理论在中子各向异性强的堆芯计算中具有较低的精度，Quasi-diffusion方程相比于传统扩散方程引入更少的近似，通过艾丁顿因子描述传统扩散理论不能反映的中子流各向异性特点。国内外对六角形几何三维Quasi-diffusion方程研究有所不足。针对艾丁顿因子计算的非线性特点，本文基于“两步法”的思想，将艾丁顿因子看作是一个特殊的少群参数，采用中子能谱适应性更好的蒙特卡罗程序SERPENT计算，并对传统扩散变分节块法进行拓展，开发了六角形组件堆芯计算程序VNMQD。采用3D VVER1000基准题、RBWR单组件问题、3D BN600简化模型对程序进行了验证，结果表明：VNMQD程序开发正确，对于非均匀性较强的堆芯，VNMQD比传统扩散方法计算精度更高、计算效率接近，实现了计算精度和计算效率的平衡。     

钴调节棒更换后RFSP-IST程序通量计算不确定性分析

在中核核电运行管理有限公司秦山第三核电厂（简称秦山三核）调节棒组件变更设计的物理分析中,用于堆芯计算程序RFSP-IST的钴调节棒增量截面由DRAGON产生,它的方法模型与秦山三核安全分析报告RFSAR（2007版）所采用的超栅元计算程序MULTICELL不完全相同,因此有必要对调节棒组件变更前后RFSP-IST程序通量计算不确定性进行分析。基于秦山三核1、2号机组的相关历史运行数据,采用95/95单边上限不确定性分析方法,对调节棒组件变更前后RFSP-IST程序通量计算不确定性进行分析。数值计算结果表明,调节棒组件变更设计及超栅元增量截面计算程序变更未对RFSP IST程序通量计算不确定性产生影响。     

基于NECP-Bamboo程序的商用压水堆乏燃料组件核素成分分析

乏燃料组件核素成分的精确计算是乏燃料临界安全分析等工作的输入条件，放射性源项计算是乏燃料组件核素成分分析的典型应用。国内现有程序由于存在数据库中核素种类不全、辐照过程无法完全模拟等弊端，限制了乏燃料后处理安全分析的可靠性和经济性。本文基于完全自主化的压水堆堆芯分析软件NECP-Bamboo,研发了商用压水堆乏燃料组件核素成分计算程序Bamboo-SFuel,利用辐照后实验(PIE)实测数据对核素成分进行了定量验证与分析，通过与Scale程序包计算结果进行对比验证了程序源项计算的精度，还探究了不同燃耗数据库对核素成分和源项计算结果的影响。数值结果表明，Bamboo-SFuel能精确分析不同辐照条件下商用压水堆乏燃料组件的核素成分和放射性源项，使用NECP-Bamboo程序中不同核素数目的燃耗数据库对重要核素成分计算结果影响不大，但对总的放射性源项计算结果影响较大；基于内置的包含1 547种核素的燃耗数据库，该程序可同时给出可靠的乏燃料临界安全分析和辐射安全分析关注的重要核素成分。     

Bamboo程序在方形组件压水堆中的适用性验证研究

本文对国产自主化核设计程序Bamboo程序在方形组件压水堆中的适用性开展了研究。主要内容包括：利用Bamboo程序对三种国内典型的方形组件压水堆进行建模,并将堆芯临界硼浓度、堆芯功率分布、堆芯燃耗分布、控制棒价值和反应性系数等参数的计算结果,同堆芯实测结果或SCIENCE程序计算结果进行对比验证。结果表明,Bamboo程序在典型方形组件压水堆堆芯计算中具有良好的精度,能够满足核电厂的堆芯核设计需求。研究结果为Bamboo程序进一步的工程应用奠定了基础。     

中子输运计算界面流方法的微扰计算

以六角形几何中子积分输运计算界面流算法及其相对应的数学共扼方程计算为基础,利用微扰原理计算了当反应堆六角形组件中栅元核参数发生微量变化时系统反应性的变化。计算结果表明,本文所开发的基于六角形几何中子积分输运算法的微扰计算方法是正确的。