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

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

应用不连续因子修正的六角形解析节块方法

在六角形解析节块方法中引入不连续因子,对该方法进行改进研究,并研制了分析程序HANDF和HANDF-C。应用该程序对UO2和MOX燃料组件组成的基准题进行了计算,结果表明,改进后的方法能够有效地提高反应堆堆芯(特别是非均匀性较强的堆芯)功率分布和有效增殖因子的计算精度。     

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

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

超临界水冷堆环形燃料组件核热耦合分析

在超临界水冷堆预概念设计中,组件设计是十分重要的,将影响堆芯性能。超临界水冷堆中水密度变化剧烈的特性要求必须进行核热耦合分析。从中子学及热工性能角度,使用三维核热耦合程序对环形燃料组件进行了优化设计。应用中子学计算程序FENNEL-N对环形燃料组件进行三维扩散计算,可得到组件内单棒功率分布,应用热工计算程序SUBSC对组件进行子通道分析。在计算过程中,分析了燃料棒间距及燃料棒与组件壁盒之间的间隙对组件性能的影响。计算结果显示,增大棒间距和棒壁间隙能提高组件k_inf,但会增大组件内功率峰因子;子通道受热不均匀性对组件热工性能影响较大,通过加入定位格架的方式能展平冷却剂出口温度,降低最大包壳温度。对环形燃料组件的安全分析表明,从中子学角度该组件是安全的。     

基于燃料棒层面核热耦合的超临界水堆堆芯研究

针对超临界水堆堆芯内流体物性分布非均匀性显著、核热反馈强烈的特点,建立了适用于超临界水堆运行环境的、基于燃料棒层面的精细化堆芯中子学/热工水力耦合方法,开发了子通道程序NCEDSCWR、节块扩散计算程序MRAPS、多功能程序COUPLE,结合西屋公司组件能谱计算程序PARAGON,构建了堆芯中子学/热工耦合分析程序系统SCAP。以具有121盒燃料组件的超临界水堆堆芯进行模拟分析,研究了堆芯三维功率分布和流体物性分布的特点以及反应性参数与重要同位素密度等随燃耗的变化规律。结果表明,本文提出的精细化核热耦合方法和开发的程序系统可以应用于超临界水堆堆芯的研究与分析,相关研究结果对超临界水堆堆芯设计具有一定的指导意义。     

堆芯瞬态耦合模拟软件CTSS的开发与验证

瞬态堆芯耦合模拟软件CTSS V1.0是以节块法堆芯中子学计算软件NACK V1.0、热工水力子通道软件CORTH V2.0、燃料元件性能分析软件FUPAC V1.1为模块的耦合软件,用于模拟典型压水堆堆芯性能,计算瞬态运行物理、热工、燃料等专业参数。堆芯三维时空中子动力学软件NACK V1.0采用粗网节块法进行堆芯扩散计算,为子通道模块和燃料性能分析模块提供堆芯精细功率。CORTH V2.0用于计算反应堆堆芯冷却剂的温度和密度。FUPAC V1.1用于模拟燃料棒在堆内的热力学行为以及计算燃料棒有效温度。NEACRP-L-335压水堆基准问题验证计算结果表明,CTSS V1.0的计算结果与国际基准程序PARCS总体符合较好。     

六角形组件内精细功率分布的计算方法与程序

六角形组件内精细功率分布的计算方法与程序王侃（清华大学）谢仲生，陈学俊（西安交通大学）关键词六角形组件，精细功率分布，扩散方程，节块法１引言由节块法程序仅能得到关于每个组件的平均信息。而从热工计算或安全分析的角度考虑，要求知道组件内各个栅元的通量或功...     

压水堆堆芯燃料装载格式的自动优化-MSOFFL程序

利用优化方法寻找压水堆最优的燃料装载格式,优化的目标函数为堆芯的功率峰因子,用物理原则指导下的直接搜索法寻找优化解。通过堆芯内燃料组件互换来搜索功率峰因子最小的燃料装载格式。搜索分两段进行:第一阶段是布置堆芯内部的新组件,以得到具有较小功率峰因子的合理装载格式;第二阶段是在上述合理装载格式基础上,布置堆芯内有燃耗的组件,以进一步降低堆芯的功率峰因子,得到优化装载格式。在优化过程中,可分别采用或随意组合采用3/2群粗网扩散理论和两群节块格林函数扩散理论计算每次换料后的反应堆两维功率分布及其它反应堆状态量,不仅保证了结果的准确性,而且可以节省CPU时间。除上述优化计算外,还能作燃耗、调临界硼浓度等计算。计算程序是采用FORTRAN-Ⅳ算法语言编制。     

基于格林函数节块法的物理与热工-水力耦合方法研究

开发了三维物理与热工-水力耦合的PWR堆芯瞬态分析程序NGFMN-K/COBRA-Ⅳ/COBRA-Ⅳ(NCC)。少群时空中子动力学计算采用格林函数节块法程序NGFMN-K,隐式耦合子通道程序COBRA-Ⅳ实现瞬态计算。采用P10H8B功率重构方法给出热组件栅元功率分布,耦合另一个COBRA-Ⅳ程序模块,进行热组件子通道分析得到安全参数。对NEACRP-L-335 C1弹棒基准问题的计算表明,NCC程序的计算结果与参考结果符合很好,说明程序计算正确,可用于评估事故结果。     

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

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

Development of a Method for Evaluation of Pin-wise Power Distribution in Fuel Assemblies of Fast Reactors

In the design of fast reactor core with higher burnup and higher linear power, prediction accuracy of burnup history of fuel pin should be upgraded so as to assure fuel integrity without extra design margin under increased neutron fluence and burnup. A method is studied to predict fuel pin-wise power and its burnup history in fast reactors accurately based on an analytic solution of diffusion theory equation on hexagonal geometry with boundary condition from core calculation by finite-differenced diffusion calculation code. The present method is applied to a fast reactor core model, and its accuracy in predicting fuel pin power is tested. The result is compared with the reference solution by the finite difference calculation with very fine mesh. It is found that the present method predicts the power peaking factors in fuel assemblies accurately. The fuel pin-wise nuclide depletion calculation is also done using neutron fluxes for each fuel pin. The result shows that the fuel pin-wise depletion calculation is very important in predicting the burnup history of the fuel assembly in detail.     

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

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

反应堆堆芯先进中子学模拟软件SCAP-N研发

堆芯中子学计算是反应堆设计分析的基础,为提高堆芯中子学计算的模拟分辨率与计算精度,开发了反应堆堆芯先进中子学模拟软件(SCAP-N)。该程序首先根据轴向特征对堆芯进行分层,并逐层进行二维堆芯非均匀输运计算,再采用超级均匀化方法(SPH)获得栅元等效均匀化截面,最后进行三维堆芯逐棒(pin-by-pin)输运计算,获得堆芯有效增殖因子与精细棒功率分布。为提高程序计算效率,采用分布式/共享式(MPI/OPENMP)混合并行方式对程序进行了并行化开发。利用虚拟反应堆(VERA)系列基准例题及美国先进非能动压水堆(AP1000)启动物理试验实测数据对程序进行了测试验证。结果表明,相比于商用核设计程序系统,SCAP-N程序采用的逐棒输运技术能够提高堆芯中子学的计算精度。与同类型高精度中子学程序相比,SCAP-N具有更高的计算效率,可进一步提高核电厂的经济性及运行灵活性。     

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

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

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

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

基于任意三角形网格的中子扩散变分节块法

传统的基于矩形和六角形几何的堆芯计算程序已不适用于具有复杂几何的新型反应堆堆芯计算,本文开展了基于任意三角形网格的多群中子扩散变分节块方法研究。首先,采用ANSYS软件对计算区域进行三角形网格剖分,并利用坐标变换将任意三角形变换为正三角形;其次,采用Galerkin变分技术建立包含节块中子平衡方程的泛函,将三角形节块内变量利用正三角形内正交基函数进行展开;最后,利用变分原理,获得中子通量密度与节块边界上分中子流的响应关系,并基于传统的源迭代法对其进行求解。基于上述理论模型开发了程序TriVNM,并采用不同几何基准题进行了验证。结果表明,TriVNM计算的堆芯keff和归一化功率分布与参考解吻合较好,该计算方法适用于复杂几何堆芯扩散计算。     

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

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

船用反应堆堆芯精细功率重构计算

针对现有船用反应堆安全分析仿真软件不能计算堆芯精细功率分布这一缺陷,开展了堆芯径向和轴向精细功率分布重构计算和分析。采用三次样条插值法对轴向精细功率进行重构计算,采用双线性和双三次插值法对径向精细功率进行重构计算,并与采用细网差分的专业物理程序的计算结果进行比较。结果表明,本工作精细功率重构计算简单、可靠,有较高的精度,对船用反应堆安全运行分析和监督管理具有重要的参考价值。