弥散颗粒毒物的多尺度耦合燃耗算法

弥散颗粒毒物具有特殊的空间结构和较强的空间自屏效应,在燃耗过程中容易出现微观分层现象。在数值模拟中,直接精细化求解将带来巨大的计算量和网格密度,具有一定的挑战。为此,本文提出了一种全新的基于多尺度耦合的燃耗计算方法:通过微观精细球层模型和宏观均匀栅元模型的耦合,将弥散颗粒介质的精细求解问题简化成对一个简单常规介质的快速求解,解决了弥散颗粒毒物在全局范围内精细燃耗求解过程中所面临的计算量大和网格密度高等问题,准确地表征弥散颗粒毒物的燃耗特征,为弥散颗粒介质的求解提供了一个新思路。经过初步验证,该算法在有效增殖因数、中子通量密度和核素核数密度等中子物理参数中有较好的表现。

弥散颗粒燃料及可燃毒物双重非均匀特性初步研究

弥散颗粒燃料及可燃毒物由于其固有安全性及自屏效应而被广泛关注,但其双重非均匀性为中子学计算带来挑战。为了研究弥散颗粒系统的双重非均匀性大小,评价体积均匀化方法的适用性,本文针对弥散不同类型、不同相体积、不同颗粒尺寸的颗粒以及不同富集度燃料基体的栅元系统进行了分析,评价栅元系统的颗粒模型与体积均匀化模型在零燃耗下的反应性偏差。分析结果显示,对于弥散燃料颗粒,体积均匀化方法的计算偏差随弥散颗粒尺寸的增加、燃料富集度的增加、以及弥散颗粒相体积的减小而增大;对于弥散可燃毒物颗粒,体积均匀化方法的计算偏差随弥散颗粒的颗粒尺寸的增加、基体燃料富集度的减小、弥散颗粒相体积的增加、以及弥散颗粒吸收截面的增大而增大。同时本文给出了弥散颗粒的双重非均匀性大小的大致顺序,针对双重非均匀性最小和最大的两种毒物颗粒也进行了详细分析,给出了是否需要考虑其双重非均匀性的大致判定条件,为弥散颗粒系统的数值计算提供指导。

双重非均匀系统反应性等效物理转换方法

本文提出用环形反应性等效物理转换(RRPT)方法和杂交反应性等效物理转换(HRPT)方法分别处理含有单颗粒类型的双重非均匀系统和同时含有弥散燃料及弥散可燃毒物的双颗粒类型的双重非均匀系统,通过典型算例,计算了系统的k inf及反应性,并以蒙特卡罗程序颗粒模型(GM)的基准解为参考值,给出了本文方法、反应性等效物理转换(RPT)方法、改进反应性等效物理转换(IRPT)方法及体积均匀化(VH)方法的计算结果与参考值之间的计算偏差。结果表明,对于仅含有弥散燃料和吸收截面相对较小的弥散可燃毒物的双重非均匀系统,RPT方法和RRPT方法的计算精度相当,计算精度高于VH方法;对于含有吸收截面相对较大的弥散可燃毒物的双重非均匀系统,RRPT方法和HRPT方法的计算精度均高于RPT方法、IRPT方法及VH方法。

双重非均匀系统环形RPT计算方法研究

弥散燃料与弥散可燃毒物由于具有双重非均匀性,采用传统体积均匀化方法(VHM)会带来较大的计算偏差。反应性等效物理转换(RPT)方法被应用于含弥散燃料的双重非均匀系统,具有方法简单且计算精度较高的特点。本文首先对传统RPT方法和改进RPT(IRPT)方法进行了分析和验证,结果表明,这2种方法对于含有弥散可燃毒物的双重非均匀系统燃耗过程中依然存在相对较大的计算偏差;然后提出环形RPT(RRPT)方法和2步环形RPT(TRRPT)方法分别用于处理含单一颗粒类型和含2种颗粒类型的双重非均匀系统,通过含不同类型可燃毒物的算例验证并与蒙卡颗粒模型基准解对比可知,本文提出的RRPT方法和TRRPT方法可用于处理含弥散燃料和弥散可燃毒物的双重非均匀系统,相比传统方法具有更高计算精度和更广适用范围。