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

堆外中子探测器的读数包含有堆芯功率分布的信息,探测器读数与堆芯功率分布之间的联系可由探测器空间响应函数建立.准确的堆外探测器三维空间响应函数是实现由堆外探测器读数重构出堆芯功率分布的关键之一.本文利用蒙特卡罗程序(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的计算效率较低,得到三维分布的空间响应函数存在较大的统计误差;堆外探测器读数主要由正对探测器的堆芯局部区域的高能中子产生.     

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

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

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

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

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

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

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

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

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

以子区内中子源为常源近似条件下的中子积分输运方程为前向方程,推导了中子积分输运方程界面流算法在六角形几何情况下的数学共扼方程;介绍了该数学共扼方程求解的内、外迭代策略,并对前向方程计算程序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程序应用于堆内释热率计算进行讨论。     

CEFR探测器三维空间响应函数敏感性研究

本文利用离散纵标程序研究了网格划分、能群数目、角度离散数和散射截面展开阶数等计算条件对堆芯区域共轭中子注量率分布形状的影响。基于求解共轭输运方程计算了中国实验快堆钠池内探测器三维空间响应函数,分析了控制棒位置、燃耗累积、燃料组件装载等因素对探测器空间响应函数的影响。结果表明：共轭中子注量率分布形状对网格宽度的敏感性低,可利用粗网格条件下共轭中子注量率分布形状求解探测器空间响应函数;控制棒对空间响应函数的影响与控制棒在堆芯中所处位置有关,组件内空间响应函数受控制棒影响程度与组件和控制棒的相对位置有关;探测器空间响应函数受新装载组件影响较明显,但对新装载组件位置的敏感性低。本文结果为大型快堆探测器空间响应函数计算提供了参考,为快堆动态刻棒提供了技术支持。     

三维结构网格SN输运求解的指数短特征线空间离散方法

空间离散格式选取和空间网格划分对于离散纵标法求解粒子输运问题的计算精度至关重要。短特征线离散格式根据粒子出射、入射贡献关系分割计算网格,可有效消除因空间离散失真导致的非物理振荡问题。指数短特征线基于指数分布假设适合求解深穿透、强衰减的屏蔽问题,已有研究中变量代换处理过程存在几何运算冗余、非线性求解成本高等问题,影响计算效率和实际应用前景。本文基于三维结构网格提出新的指数短特征线离散格式,改变分割变换方案避免了单元投影映射过程,应用并调整指数矩算法、优化指数因子求解等处理改善计算效率。通过4个简单几何固定源问题定性定量比较分析了若干离散格式的计算精度与效率,指数短特征线方法粗网格精度远优于其他方法,出色的粗网格精度可补偿其较高的单网格计算成本,恰当的网格划分可发挥离散方法的计算效率。     

蒙特卡罗共轭输运法计算反应堆压力容器快中子注量率

为了保证压力容器（RPV）在核电厂寿期内的安全,通过理论方法准确评估其受到的快中子积分注量率非常重要。本文提出了一种应用共轭输运理论解决深穿透问题的计算方法,并将该方法的计算结果与基准题HBR-2给出的实测值及确定论方法的结果进行了比较。结果表明：本文计算结果与基准题给出的实测数据吻合良好,大多反应率计算相对误差小于10%,最大相对误差不超过35%;70%以上的计算结果准确性优于确定论方法,表明本文提出的解决蒙特卡罗深穿透问题的方法是有效且准确的。