蒙特卡罗共轭输运法计算堆外探测器三维空间响应函数

堆外中子探测器的读数包含有堆芯功率分布的信息,探测器读数与堆芯功率分布之间的联系可由探测器空间响应函数建立.准确的堆外探测器三维空间响应函数是实现由堆外探测器读数重构出堆芯功率分布的关键之一.本文利用蒙特卡罗程序(MCNP),采用共轭输运的计算方法,以压水堆为例,计算分析堆外探测器的三维空间响应函数,并与以往二维离散纵标法的计算结果进行比较.结果表明采用三维计算是必要的,而蒙特卡罗共轭输运法可以方便地实现三维复杂几何下的空间响应函数计算.     

CAP1000核电厂堆外探测器响应函数计算方法研究

堆外探测器响应函数代表了堆芯活性区各组件对堆外探测器读数的贡献,反映了堆芯功率分布与探测器读数的关系。本文利用二维离散纵标法(S_N)程序DORT,研究其共轭输运方法,建立CAP1000反应堆模型,分析其堆外探测器径向和轴向响应函数及其特性,并与采用DORT程序正向输运计算的结果进行比较。研究表明,共轭输运方法可以极大简化计算量,且计算结果与正向输运方法结果符合较好。     

CAP1000核电厂堆外探测器响应函数计算方法研究

堆外探测器响应函数代表了堆芯活性区各组件对堆外探测器读数的贡献,反映了堆芯功率分布与探测器读数的关系。本文利用二维离散纵标法(SN)程序 DORT,研究其共轭输运方法,建立 CAP1000反应堆模型,分析其堆外探测器径向和轴向响应函数及其特性,并与采用 DORT 程序正向输运计算的结果进行比较。研究表明,共轭输运方法可以极大简化计算量,且计算结果与正向输运方法结果符合较好。     

堆外探测器空间响应函数计算方法研究

通过堆外探测器空间响应函数与反应堆功率分布的内积可得到探测器电流。因此,获得堆外探测器的空间响应函数对堆外探测器电流信号的刻度和修正具有重要意义。本文利用二维输运计算程序DORT和基于ENDF/B-Ⅶ.1制作的多群数据库MATXS-47,采用求解共轭中子输运方程的方法对压水堆核电厂Indian Point 2的堆外探测器径向空间响应函数进行了计算,计算结果与文献值吻合,表明本文所采用的响应函数计算方法是正确的。     

去射线效应堆外探测器响应函数计算研究

堆外探测器响应函数表征了堆芯活性区各位置处的裂变中子对堆外探测器响应的贡献,通过共轭S_N输运计算可快速得到堆外探测器的响应函数。然而,堆外探测器远离堆芯且相对于堆芯体积很小,S_N方法的计算结果会受到射线效应的影响。为解决堆外探测器响应函数计算中的射线效应问题,研究了共轭首次碰撞源射线效应消除方法。此外,为克服共轭首次碰撞源方法在三维堆芯计算中面临的计算量大、内存需求高等问题,研究了共轭首次碰撞源的并行化计算方法和动态内存管理方法。基于韩国Kori-1压水堆的计算结果表明:共轭首次碰撞源S_N方法和多群蒙特卡罗方法具有相当的计算精度,但计算效率高1个量级。     

高温气冷堆堆外探测器空间响应函数的计算和特性分析

建立高温气冷堆核电厂示范工程( HTR-PM)反应堆堆外探测器空间响应函数的计算模型.基于共轭输运理论,分别使用蒙特卡罗程序( MCNP)和三维离散纵标(SN)程序TORT计算高温气冷堆堆外探测器空间响应函数.对堆外探测器空间响应函数主要特性的分析及对2种计算结果的比较表明:SN程序TORT的计算结果和MCNP的计算结果一致;基于共轭中子输运理论建立堆外探测器读数和堆芯功率分布的映射关系(探测器空间响应函数)可行;MCNP的计算效率较低,得到三维分布的空间响应函数存在较大的统计误差;堆外探测器读数主要由正对探测器的堆芯局部区域的高能中子产生.     

三维平几何各向异性散射问题的节块离散纵标输运计算

本文给出了三维（ｘ、ｙ、ｚ）几何各向异性散射问题的节块离散纵标输运计算模型，并研制了相应软件ＮＯＴＲＡＮ／３Ｄ，得到了满意的计算结果。     

角度自适应离散纵标输运计算方法

精确的屏蔽计算方法是核装置辐射屏蔽设计的重要基础,离散纵标法(SN)是主要的屏蔽计算方法之一。本文基于价值理论的目标导向与角度自适应相结合的方法,有效地减弱了角度的离散误差。求解输运共轭方程获得目标函数的重要性分布,采用局部角度离散误差与目标函数的重要性加权,产生后验误差估计,为角度自适应过程提供判断依据。角通量密度的映射采用多项式权重法和球谐函数拟合法。数值结果表明,对于具有直孔道或曲折孔道的屏蔽问题,在相同精度下离散角度数减少了1～2个数量级,极大地减少了计算量。角度自适应方法以较少的离散方向获得了准确的计算结果,有效地减弱了角度离散误差对屏蔽计算精度的影响。     

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

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

二维中子输运问题的离散纵标短特征线方法研究

离散纵标法是求解中子输运方程的主要数值方法之一,空间变量离散及误差控制对保证输运计算精度至关重要。传统有限差分离散方法对于特定模型会产生非物理振荡问题,粗网精度不足使得低阶差分方法的应用具有局限性。本文研究了二维常数和线性短特征线方法,短特征线空间离散基于中子输运的特征线解,根据输运方程的空间矩守恒构造网格角通量密度完成输运方程求解。选取固定源和临界问题进行测试验证并分析了网格敏感性。数值结果表明,线性短特征线离散对网格敏感性较低,较常数短特征线和低阶差分方法具有更高的计算精度及效率。     

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

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

CR-39测量中子能谱响应函数的计算

计算了在一定的化学和电化学蚀刻条件下,由聚乙稀、铅(金)和CR-39组成的叠层探测器的响应函数(1～20 MeV)。计算结果表明,叠层探测器的响应函数分布近似为高斯分布,由铅组成的探测器的响应函数的半高宽比由金组成的探测器的响应函数的半高宽窄,由前者组成的探测器测量中子能谱分辨率高于后者。     

三维离散纵标方法在堆内构件释热率计算中的初步应用

以美国H.B.Robinson-2#机组反应堆压力容器(RPV)基准实验的参数为输入数据,采用三维离散纵标方法程序(TORT)计算压力辐照监督管处中子能谱及典型核素的活度值。计算得到的辐照监督管处中子能谱与基准实验结果趋势一致、吻合较好;典型核素活度的计算值与测量值之比(C/M)为1.04±0.04。用TORT对福建宁德核电站堆内构件释热率分布进行初步计算,并与蒙特卡罗方法(MCNP)的计算结果相比较,两种方法的结果表现出良好的一致性。最后对TORT程序应用于堆内释热率计算进行讨论。     

探测器时空响应函数对信噪比的影响

理论推导了时空分辨探测系统的信噪比表达式,研究了时空响应函数对探测像素信号信噪比的影响。计算结果表明时空响应函数对信噪比的影响依赖于辐射源的性质和探测器的增益特性。     

An Improvement of the Transverse Leakage Treatment for the Nodal SN Transport Calculation Method in Hexagonal-Z Geometry

To overcome the divergent behavior of the NSHEX code, a nodal S_N code for hexagonal geometry, for some transport calculations, an improvement has been made in the calculation of the axial leakage. The axial leakage, previously calculated by using the quadratic transverse leakage approximation (QLA), is calculated by a new method of analytically treating the spatial distribution within a node, based on the axial homogeneity of the ordinary core. The verification tests were performed for the KNK-II model geometry of the NEACRP 3-D Neutron Transport Benchmarks and the large assembly-size KNK-II model. It is found that k_ett values obtained by introducing the new method agree with the reference Monte Carlo calculation results within 0.1% Δk/k for the KNK-II model, although the QLA method did not converge for two cases. Furthemore the new method succeeded in calculations for the large assembly-size model, in which the QLA method failed for all cases. Thus the new method has been found accurate and convergence achieved for the cases in which the QLA method failed.