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

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

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

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

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

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

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

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

三维六角形节块多群中子扩散程序NDHEX

本文介绍用DIF3D (NOD)求解二、三维六角形几何系统下中子扩散方程的理论模型及数值计算方法。六角形节块内的中子通量密度分布采用高次多项式近似表示,最后导出通量矩方程及偏流的响应矩阵方程。应用粗网再平衡和渐近源外推方法加速收敛。参考此方法编制了计算程序NDHEX,并对一些六角形基准问题进行了计算。结果表明:NDHEX的计算结果与DIF3D(NOD)的计算结果符合很好;与差分程序相比,具有更高的精度与计算效率。它可用于快堆计算。     

轻水反应堆燃料组件中子通量和功率分布的计算程序

本文介绍用穿透几率法计算二维轻水堆燃料组件内中子通量分布的两种计算模型和程序.在子区内及表面上中子通量采用线性空间分布近似,子区表面上角通量分别采用准 DP_1和 QP_1近似。对一些轻水堆组件基准问题作了验证计算。计算结果与 S_N、节块 S_N 以及积分输运理论等方法进行比较,其结果符合良好。这些程序可用于轻水堆燃料组件的计算。     

注量展开节块法求解三维六角形几何中子扩散方程

提出了一种在三维六角形几何节块内数值求解中子扩散方程的节块法该方法把节块内各群中子注量分布用解析基函数近似展开为了改善节块耦合关系．提出了,一种新的节块边界条件：面平均偏流零次矩和一次矩同时保持连续。此外．将响应矩阵技术应用于迭代求解过程,使得该方法具有较高的计算效率基于本文提出的模型,发展了三维六角形组件中子扩散计算程序FEMHEX。通过对二维、三维VVER基准问题校验计算表明,该方法能高效．准确的给出有效增殖系数以及节块功率分布。     

解析基函数展开方法求解二维三角形几何中子扩散方程

提出了一种在二维三角形几何内数值求解中子扩散方程的节块方法.节块内的各群通量分布用解析基函数近似展开,节块之间采用面偏流零次矩和一次矩进行耦合;给出了三角形几何下的节块扫描方案;采用响应矩阵技术进行迭代求解,开发了二维三角形组件中子扩散计算程序ABFEM-T.通过基准问题的校验计算,表明该方法能准确地给出有效增值系数及节块功率分布,可用于复杂的非结构几何区域的中子扩散问题的求解.     

穿透几率法求解X-Y几何多群中子积分输运方程

本文介绍应用穿透几率法求解X-Y几何多群中子输运方程的程序—TPXY。把系统分成若干子区,子区内部中子通量呈线性分布,子区表面通量角分布采用4P_1近似,空间呈线性分布。对一系列组件基准问题作了校验计算,计算结果与S_N,碰撞几率及其它穿透几率方法的结果进行了比较,均符合良好。本程序可用于轻水堆组件计算。     

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

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

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

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

中子输运计算界面流方法的数学共扼方程

以子区内中子源为常源近似条件下的中子积分输运方程为前向方程,推导了中子积分输运方程界面流算法在六角形几何情况下的数学共扼方程;介绍了该数学共扼方程求解的内、外迭代策略,并对前向方程计算程序TPHEX进行了改造,得到了常源近似情况下数学共扼方程计算程序TFHEX_J0。通过算例校验表明TPHEX_J0与TPHEX程序所计算的系统本征值符合良好,TFHEX_J0程序的计算结果是可靠的。     

三区非均匀栅元X-Y几何块方法中子输运计算

基于中子积分输运理论,应用综合界面流和碰撞几率技巧的块方法,导出了处理三区非均匀栅元结构的二维（Ｘ－Ｙ）几何多群中子输运问题的数值模型。即对于由若干栅元组成的按Ｘ－Ｙ几何排列的堆芯结构,对每一类栅元剖分为圆柱形元件（如燃料棒、控制棒、可燃毒物棒等）、包壳和慢化剂三个均匀区,用碰撞几率（ＣＰ或ＰＩＪ）方法计算各区的中子通量分布;对于相邻栅元用ＤＰ１近似的中子流来耦合;因此,块方法具有精度高、速度快、能灵活处理各种几何问题的优点,是目前动力堆组件计算最有前途的方法之一。基于块方法基本理论,发展了三区栅元模型,导出了计算方法,编制了ＦＯＲＴＲＡＮ计算机程序。为验证其精度和适用性,对两个例题进行了计算,并与其它程序的计算结果进行了比较,证明功率分布和本征值均符合较好。     

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.     

Analysis of NEA C5MOX benchmark with computer code “cohint” based on interface current approach

Interface current approach to solution of neutron transport equation has been earlier used for LWR lattice problems. The analysis of the C5MOX benchmark is a opportunity to test its applicability to heterogeneous reactor problem. Computer code COHINT, which incorporates a routine for solution of neutron transport equation in X-Y geometry by interface current approach was used for this analysis by representing the fuel rod as a square one. Region interface angular fluxes are represented by a four-term expansion. It is found that the average pin power error is about 2.28 % (peak pin error 4.1 %) relative to reference calculations. Further improvement is possible by introduction of the capability to represent circular rods with in a square cell in COHINT.     

基于先进节块法的六边形三维中子物理程序的开发与验证

六边形燃料组件在液态金属冷却快堆尤其是钠冷快堆中被广泛应用,针对这类堆型的设计与安全分析需要对堆芯中子通量与中子流进行三维全堆芯耦合计算。经过多年发展,目前已有多种解析节块法、积分节块法、节块展开法等先进节块法能在笛卡尔坐标系下较为精确求解多维中子扩散方程。本文通过径向半解析节块法耦合轴向高阶节块展开法的综合节块方法开发了反应堆三维中子物理计算软件SA HNHEX,并对VVER 440二维、三维基准题进行建模与仿真计算。计算结果与参考值符合较好,初步验证了使用该方法进行反应堆堆芯中子扩散计算的正确性。     

A Method of Calculation of Detailed Power Distribution in Fuel Assemblies

Abstract

A method has been developed that effectively estimates the detailed distribution of power generation in the fuel or blanket assemblies in nuclear reactors. A two-dimensional, one-group diffusion model is applied to a region of homogeneous composition enclosed in a contour devoid of concavity viewed from outside. The diffusion equation is reduced to the form of Helmholtz equation, and a non-homogeneous boundary condition of Dirichlet or Neumann type is given on the contour, using neutron fluxes previously obtained in coarse mesh diffusion criticality calculations covering the whole reactor. This boundary value problem in two-dimensional space is made to yield a solution in the form of a potential due to a single or double layer. The method is applied to a hexagonal cell of a fast reactor. The results of calculation are amply accurate in comparison with the corresponding values from the usual fine-mesh diffusion scheme and with much shorter computing time.     

Analysis on the neutron kinetics for a molten salt reactor

The neutron kinetics of the molten salt reactor is significantly influenced by the fuel salt flow, which leads to the analysis methods for the conventional reactors using solid fuels not being applicable for the molten salt reactors. In this study, a neutron kinetic model considering the fuel salt flow is established based on the neutron diffusion theory, which consists of two-group neutron diffusion equations for the fast and thermal neutron fluxes and six-group balance equations for delayed neutron precursors. The temperature feedback on the neutron kinetics is considered by introducing a heat transfer model in the core, in which the group constants which are dependent on the temperature are calculated by the code DRAGON. The mathematical equations are discretized and numerically calculated by developing a code, in which the fully implicit scheme is adopted for the time-dependent terms, and the power law scheme is for the convection–diffusion terms. The neutron kinetics is conducted during three transient conditions including the rods drop transient, the pump coastdown transient and the inlet temperature drop transient. The relative power changes and the distributions of the temperature, neutron fluxes and delayed neutron precursors under these three different transient conditions are obtained in the study. The results provide some valuable information for the research and design of this new generation reactor.