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

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

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

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

离散纵标屏蔽计算空间自适应方法研究

离散纵标法是目前国际上通用的核装置辐射屏蔽计算的主要方法之一。空间变量离散误差对离散纵标屏蔽计算的精度至关重要。实际屏蔽计算中存在较强的非均匀性,传统的空间离散方法不能高效地在全局达到较高精度。本文采用线性间断有限元方法进行空间离散,由基于两网格的误差估计和基于残差的误差估计驱动空间网格自适应细化;基于树状结构的六面体网格,对传统的输运扫描进行改进,在粗-细网映射中保持零阶、一阶空间矩守恒。数值结果表明：线性间断有限元方法具有较好的射线传递特性和空间收敛性;自适应算法能较好地对通量密度间断、通量密度梯度大、光学厚度大的区域进行网格细化。在达到相同计算精度时,自适应方法所需的网格数较均匀细分的方法减少约1个数量级,能有效提高屏蔽计算效率。     

三维离散纵标方法在中国实验快堆堆芯屏蔽计算中的应用

离散纵标(又称SN)方法是反应堆屏蔽计算中常用的方法,随着计算能力的发展和离散纵标计算方法的不断完善,使得离散纵标方法在反应堆的屏蔽计算中得到了广泛的应用。本文以中国实验快堆(CEFR)堆芯为研究对象,使用三维离散纵标方法对区域功率份额、组件功率、DPA、寿期内堆芯围板积分快中子注量及寿期内小栅板联箱积分快中子注量进行了计算研究,并与二维离散纵标法和俄罗斯设计报告结果进行对比。研究结果表明:三维离散纵标方法能够减少二维程序几何等效过程中导致的误差,计算结果可靠,可应用于大型快堆堆芯的屏蔽设计中。     

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

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

离散纵标与蒙特卡罗耦合方法在反应堆屏蔽计算中的应用

建立基于离散纵标方法计算程序DOT以及蒙特卡罗方法计算程序MCNP的耦合方法。通过为MCNP程序SSR命令提供二进制源描述文件,实现在不重新编译MCNP程序条件下的离散纵标与蒙特卡罗耦合。在600 MW核电厂反应堆堆坑小室屏蔽计算中应用耦合方法的计算结果验证该方法的正确性。在耦合方法中研究探索源偏倚策略的实现方式,计算结果表明,源偏倚技术应用与离散纵标与蒙特卡罗耦合方法计算中,能有效降低统计方差,提高计算效率。     

三维离散纵标方法在反应堆精细屏蔽计算中的应用研究

随着计算能力的发展,三维离散纵标方法(SN)已逐渐为反应堆辐射屏蔽计算的主流方法之一。本文就三维SN方法应用于反应堆的精细化屏蔽建模与分析的关键问题进行了研究,主要包括精细化几何建模、堆芯固定源模型的创建和数据库截面参数敏感性分析等内容。在此基础上,本文以典型的压水反应堆为对象,构建了精细的三维SN计算分析模型,以压力容器快中子注量率为算例,完整实现了反应堆的精细化三维SN建模与分析,并将三维SN结果和蒙特卡罗方法的计算结果进行了比较分析。对比结果表明,精细化三维SN方法具有较高的计算精度,验证了三维SN方法在反应堆精细屏蔽计算问题中的有效性和正确性。     

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

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

基于离散纵标法的三维耦合燃耗与活化计算方法的发展

燃耗与活化分析在反应堆的燃耗管理与辐射屏蔽设计分析中起关键性的作用.基于一维、二维的输运程序的燃耗与活化分析方法难以解决复杂几何和强烈各向异性散射问题.本文通过耦合三维离散纵标(SN)方法粒子输运程序以及指数欧拉法活化计算程序,发展了快速精确的三维耦合燃耗与活化计算方法.该方法考虑了共振自屏效应动态修正,并采用重要核素...     

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

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

三维离散纵标方法在RPV快中子注量率计算中的初步应用

离散纵标(又称SN)方法是反应堆压力容器快中子注量率计算中最常用的方法之一.计算机的飞速发展和三维离散纵标方法的不断完善,使得三维离散纵标法应用于工程设计成为可能.本文以某压水堆为研究对象,采用基于三维全堆芯中子学程序SCIENCE的计算结果,使用三维离散纵标方法程序计算了压力容器内表面快中子注量率三维分布,并分别与采...     

三维离散纵标-蒙特卡罗耦合方法在快中子注量率计算中的应用

采用三维离散纵标-蒙特卡罗耦合方法(SN-MC耦合方法)对压水堆压力容器的快中子注量率进行计算分析,并与离散纵标方法(SN方法)、蒙特卡罗方法(MC方法)的计算结果以及测量值进行比较。结果表明,SN-MC耦合方法的计算结果和测量值相比偏差在20%以内,满足美国核管会(NRC)法规导则RG1.19对快中子注量率计算方法偏差的要求;当SN-MC耦合方法和MC方法达到相同计数误差水平时,其计算速度比MC方法快2~10倍。     

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

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

一阶中子输运方程的离散纵标有限元方法研究

利用离散纵标(SN)方法离散SN一阶多群中子输运方程,建立了基于剖分单元的有限元离散与非结构网格扫描方法相结合的求解模型。针对给定的三角形单元应用Galerkin变分,获得线性方程组,通过引入定解条件,求解线性代数方程组,获得该三角形单元所有节点上的角通量,然后对其他三角形单元进行扫描,从而解出所有节点处中子通量密度。根据上述理论模型,编制了相应的计算程序FEGT,对一系列例题进行校验的数值结果表明,该程序的计算精度满足要求。     

基于蒙特卡罗-离散纵标双向耦合方法的堆坑粒子注量率计算

为有效解决大型复杂核设施屏蔽计算问题,研究了三维蒙特卡罗(MC)-离散纵标（S_N）双向耦合方法,通过自主开发接口程序实现MC粒子概率分布与S_N角通量密度之间的相互转换,实现MC-S_N双向耦合计算。将基于MC-S_N双向耦合方法的程序用于某反应堆堆坑底部粒子注量率计算。利用MC程序建立堆芯及堆坑处的精细模型进行计算,三维S_N程序用于堆芯下表面与压力容器底面之间区域的计算。通过MC-S_N-MC两步耦合计算,给出堆坑通道及小室内的中子和光子注量率。三维MC-S_N双向耦合方法计算结果与单一MCNP程序结果吻合较好,初步验证了该方法是解决大型复杂核装置屏蔽问题的有效工具。     

离散纵标共轭输运法计算堆外探测器空间响应函数的研究

通过对共轭输运法的研究,采用二维离散纵标程序计算某三代核电厂反应堆满功率运行状态下的堆外探测器响应函数,计算结果与原有数据相比,径向响应因子的相对偏差在3%以内,轴向响应因子的相对偏差在1%以内。表明所采用的计算模型和方法正确。将共轭输运法计算堆外探测器响应的方法在岭澳二期核电项目上进行验证,并与岭澳二期核电工程中采用正向输运法计算的响应因子进行对比,结果表明2种方法的计算偏差不超过5%,验证了共轭输运法与正向输运法计算结果的等价性,但共轭输运法计算效率更高。     

离散纵标输运计算方法在压水堆核电厂41Ar活化源项分析中的应用

压水堆核电厂功率运行期间,反应堆压力容器外的环形空腔空气中所含的⁴⁰Ar被中子活化,形成具有放射性的⁴¹Ar。文章采用二维离散纵标输运计算程序DORT分析了反应堆堆腔区域的中子注量率分布情况,采用NJOY评价核数据处理程序,根据DORT分析得到的通量作为权重通量,利用基础评价核数据库ENDF/B-Ⅶ.0制作⁴⁰Ar中子俘获反应的微观截面,在此基础上,分析了百万千瓦级压水堆核电厂每台机组反应堆堆腔空气中⁴⁰Ar中子活化生成⁴¹Ar的生成率以及电厂⁴¹Ar的环境排放源项。文章给出的⁴¹Ar源项分析方法可作为压水堆核电厂设计中确定⁴¹Ar源项的最佳估算值的参考。     

瑞利阻尼对大质量直接积分时程计算方法的影响

使用直接积分时程计算方法对反应堆冷却剂系统进行大质量点激励时,瑞利阻尼会对大质量点额外附加瑞利阻尼力,使计算结果偏小。本文从理论上讨论该瑞利阻尼力的形成原因,并给出解决方法。数值计算结果表明,此方法可在实际工程计算中应用。     

基于蒙特卡罗-离散纵标双向耦合方法的快中子注量基准分析

蒙特卡罗（MC）-离散纵标（S_N）双向耦合方法是解决大型复杂核装置屏蔽问题的有效方法。本文针对三维MC-S_N双向耦合方法在大型压水堆核电站屏蔽计算中的应用,进行了基准验证分析。基于美国核管会（NRC）发布的NUREG/CR-6115压水堆基准模型,采用自主开发的三维MC-S_N双向耦合屏蔽计算分析方法,利用MCNP4C精确计算堆芯到热屏蔽精细模型以及位于压力容器内部计算区域的精确模型,三维S N 程序TORT用于进行热屏蔽到第2下降区外表面间的计算。通过自主研发的接口程序实现MC粒子概率分布与S_N角通量密度间的相互转换,实现MC和S_N 双向耦合计算。三维MC-S_N双向耦合方法计算结果与基准报告提供的MCNP、DORT结果符合良好,初步验证了该方法解决大型复杂核装置屏蔽问题的可行性。