专利名称：一种获得核燃料组件共振参数的方法
技术领域：
本发明属于反应堆物理设计技术领域，具体涉及一种获得核燃料组件共振参数的方法。
背景技术：
核燃料组件共振参数的获取伴随着反应堆物理设计中能量多群处理方法的出现而产生。获取核燃料组件共振参数是反应堆组件设计、堆芯设计以及燃料管理的基础，其精度直接影响后续计算精度。共振计算的主要目的是使用多群数据库的多群截面和共振积分数据，计算实际核燃料组件的共振参数。各类基于确定论的核燃料组件共振参数获取方法大多受到其理论模型的假设、近似的限制，或仅仅能够处理简单几何核燃料组件；或只能近似处理有限的多种规则几何核燃料组件，并未真正实现有效获取多种任意几何、任意数量共振区、任意共振材料组成核燃料组件的共振参数；或由于可获得任意几何，任意数量共振区和材料组成共振参数但效率低下，同样距实际工程应用中存在差距。用于获取核燃料共振参数的子群方法具有较强的几何处理能力、较高的计算效率，但是无法解决多核素共振的问题，只有采用多核素共振干涉表修正的方式，这种方式随着几何复杂程度的提高会逐渐失效，逐渐不满足目前工程实际情况。小波方法具有较强的几何处理能力，理论上可以用来获取任意几何、任意数量共振区核燃料组件的共振参数，但是计算效率低下，离工程实际使用还有一段距离。因此，针对新的需要与挑战，研究一种有效获取任意几何、任意数量共振区、任意共振材料组成核燃料组件的共振参数的方法成为了当务之急。
发明内容
为解决上述问题，本发明提供一种获得核燃料组件共振参数的方法，能够有效获得任意共振材料组成，任意几何，任意数量共振区核燃料组件的共振参数；高效率的处理核燃料共振参数获得过程中的不均匀布置，栅元、组件的不规则几何形状，弹性散射共振问题。本发明的原理是通过根据各个共振核素随能量变化的共振峰的稀疏度采用不同的核燃料共振参数获取方式。在共振能区，随着能量的增大共振峰的峰值将逐渐减小，而共振的频率将逐渐增大，也就是共振峰将逐渐变密，由稀疏的孤立共振峰向密集的共振峰过渡，根据这种变化特点:在共振区高能段共振峰密集的地方，针对截面随能量的变化曲线，不对能量进行划分，而是采用横向划分的方式直接对截面进行划分；在共振区低能段共振峰稀疏共振干涉现象强烈的地方，针对能谱随能量的变化曲线，采用小波函数展开来精确拟合能谱的方式。两种划分方式引入两种不同的获取共振参数的方式，两种获取方式通过两种划分区域之间的散射进行耦合。两种获取方式分别采用现有的任意几何多群中子输运方法作为求解器，得到核燃料的空间相关的能谱和共振自屏截面。最终提高获取任意共振材料组成，任意几何，任意数量共振区核燃料组件共振参数的精度和效率。
为了达到上述目的，本发明采用如下技术方案:一种获得核燃料组件共振参数的方法，包括如下步骤:步骤一:使用NJOY程序，制作常用核素多群数据库，针对子群方法的特点添加了吸收截面插值表和总截面插值表，同时，基于ENDF数据得到核燃料中共振核素的连续能量形式的共振峰随能量的变化曲线；步骤二:在核燃料每个共振核素的变化曲线图中，将孤立峰和密集峰区分，得到两个能量区间，孤立峰所在区间是共振干涉明显的共振能区低能段，密集峰所在区间是共振干涉不明显的共振能区高能段；步骤三:在共振干涉现象不明显的共振能区高能段，通过子群、子群参数的定义，将能群平均共振截面通式转换为子群形式，然后使用一系列不同背景截面下能群平均共振截面求解子群总截面、子群分截面以及子群概率；步骤四:将核燃料每个共振核素的子群参数包括子群截面、子群概率和相应的多群数据耦合:共振能量段使用子群参数，快群和热群使用多群数据，制作子群数据-多群数据率禹合系统；步骤五:针对每个共振能群，使用子群总截面、子群分截面和子群概率简化与多群形式一致的子群二维任意几何特征线中子输运方程；考虑到共振计算的目标是得到有效共振自屏截面，在子群输运方程的求解中，不需要再整个能量段上进行，只需要针对共振能区各个能群独立的进行固定源的子群输运计算，得到子群中子能谱，通过子群中子能谱和子群参数归并得到核燃料的共振参数；考虑到某些核燃料共振核素较多，或者要获得核燃料精细燃耗期间的共振截面，针对共振核素的共振峰的范围和大小，将共振核素分为3-4类，计算过程中的迭代在类的基础上进行，提高迭代的效率；步骤六:使用共振能区高能段到共振能区低能段的群间散射，通过群间散射将共振能区高能段的计算和共振能区低能段的计算耦合；在共振能区高能段使用子群方法得到该能量段的共振参数后，计算共振能区高能段到共振能群低能段的散射源和裂变源；步骤七:利用得到的散射源与裂变源计算源项，在共振干涉现象强烈的共振能区低能区，利用小波展开方法，对连续能量中子能谱的能量变量进行小波展开，得到共振能区的连续能量能谱小波尺度函数展开系数方程组；步骤八:使用二维任意几何特征线方法计算与多群中子输运方程形式一样的小波尺度函数展开系数方程，得到共振区低能段多群中子能谱；为了提高计算的效率，对计算的散射源项进行简化，在迭代中将共振能区高能段到共振能区低能段的群间散射以及共振能区低能段的群内散射在迭代计算中作为常量，提高小波方法的计算效率。最后，使用共振区低能段多群中子能谱归并连续总截面和分截面，得到共振能区低能段的核燃料共振参数。步骤一所述的制作常用核素多群数据库的过程如下:一、使用NJOY程序加工多群数据；二、将多群数据格式转换为NECL多群数据库；三、合成燃耗中所需伪核素；四、制作批处理完成核素的添加、删除、替换功能。步骤四所述的耦合过程如下一、在NECL多群数据库的基础上，将该数据库按照共振核素与非共振核素分为共振核素多群数据库和非共振核素多群数据库；二、对共振核素多群数据库中的核素进行共振能群高能段的子群参数计算，得到子群截面和子群概率参数；三、将共振核素的共振能群高能段的子群参数与该共振核素快群热群的多群截面数据放置在一起，与非共振核素的多群截面数据一起，形成子群数据-多群数据耦合系统。本发明和现有技术相比，具有如下优点:1、避免了首次碰撞概率的求解，可以获得二维任意几何核燃料组件的共振参数；2、相对其他传统的获取方式，本发明可以得到与核燃料组件中空间相关的能谱和共振参数，提高了获取几何任意，共振材料任意，存在温度分布梯度的核燃料组件共振参数的精度；3、针对不同核燃料组件将共振能区划分为两区，在共振干涉现象不明显共振能区高能段，存在大量密集的共振峰，采用横向划分的子群方法，避免对大量密集共振峰的拟合或插值求解，高效率获取这一能区核燃料共振参数，同时避免了共振干涉现象对子群方法获取共振参数的精度影响；在共振干涉现象强烈的共振能区低能段，存在着波峰值很大的孤立峰，采用小波展开函数拟合的方式来获取核燃料共振参数，对于多核素共振干涉既不需要迭代计算也不需要进行共振干涉修正，同时小波展开函数的使用范围仅仅是共振能区低能段的稀疏孤立峰，避免了对高能区密集共振峰的拟合必须采用的高阶数，有效地保证了获取核燃料整个共振区共振参数的效率和精度。
图1为本发明方法流程图。图2为制作多群数据库的流程图。图3为共振能区高能段与共振能区低能段划分示意图。图4为子群方法截面划分不意图。图5为子群-多群耦合数据库系统示意图。图6为子群方法和小波方法散射源耦合示意图。
具体实施例方式下面结合附图和具体实施方式
对本发明结构进行详细说明。本发明的思路是针对不同的核燃料将共振能区分为高能段和低能段，在共振能区高能段采用子群方法，在共振能区低能段采用小波方法，通过高能段到低能段的群间散射耦合，然后求解二维任意几何中子输运方程M0C，最终得到核燃料组件的共振参数。如图1所示，本发明发明的步骤如下:步骤一、制作多群数据库，并且对共振核素进行分区，得到共振能区高能段与共振能区低能段；步骤二、共振核素的共振能区高能段数据参与子群参数计算，得到子群参数；步骤三、对核燃料中存在的共振核素分类，得到3-4类；步骤四、在这些核燃料中存在的共振核素的共振能区高能段进行MOC子群中子输运计算，得到子群通量，归并得到多群截面；
步骤五、在分类的基础上重复进行步骤4，直到迭代收敛；步骤六、进行共振能区高能段的MOC多群输运计算，得到该能段的共振核素多群截面，计算共振能区高能段到共振能区低能段的散射源与裂变源；步骤七、得到散射源与裂变源后，在共振能区低能段求解小波系数方程，得到小波系数矩阵；步骤八、基于小波系数矩阵，得到共振能区低能段的多群截面与多群通量；步骤九、在整个能区进行MOC多群输运计算，最后输出结果。下面进行详细说明:步骤一:使用NJOY程序，制作常用核素多群数据库一NECL，针对子群方法的特点添加了吸收截面插值表和总截面插值表；同时，基于ENDF数据得到核燃料中共振核素的连续能量形式的共振峰随能量的变化曲线。制作常用核素多群数据库一NECL的具体流程图如图2所示，过程如下:一、使用NJOY程序 加工多群数据；二、将多群数据格式转换为NECL多群数据库；三、合成燃耗中所需伪核素；四、制作批处理完成核素的添加、删除、替换功能。其中自主开发的辅助维护程序功能包括将NJOY的输出转化为适合子群方法和小波方法的耦合数据库格式；批处理进行数据删除，增加，修改功能；燃耗伪核素的自动合成功能。步骤二:在核燃料每个共振核素的变化曲线图中，将孤立峰和密集峰区分，得到两个能量区间，孤立峰所在区间是共振干涉明显的共振能区低能段，密集峰所在区间是共振干涉不明显的共振能区高能段。举U238核素说明，如图3所示，在共振区域内，在分界线左边的区域主要存在的是孤立峰，属于共振干涉明显的共振能区低能段。在分界线右边的区域主要存在的是密集峰，属于共振干涉不明显的共振能区高能段。步骤三:在共振能区高能段通过子群方法求解获得子群总截面、子群分截面和子群概率在共振干涉现象不明显共振能区高能段，存在大量密集的共振峰，采用横向划分的子群方法，避免对大量密集共振峰的拟合或插值求解，如图4所示。通过将均匀问题、非均匀问题能群平均共振截面形式一致性的通式展开为子群形式，并利用能群平均共振截面与背景截面的变化关系，求解子群参数多群总截面，背景截面，子群总截面间的关系通过如下公式(I)表示
权利要求
1.一种获得核燃料组件共振参数的方法，其特征在于:包括如下步骤: 步骤一:使用NJOY程序，制作常用核素多群数据库，针对子群方法的特点添加了吸收截面插值表和总截面插值表，同时，基于ENDF数据得到核燃料中共振核素的连续能量形式的共振峰随能量的变化曲线； 步骤二:在核燃料每个共振核素的变化曲线图中，将孤立峰和密集峰区分，得到两个能量区间，孤立峰所在区间是共振干涉明显的共振能区低能段，密集峰所在区间是共振干涉不明显的共振能区高能段； 步骤三:在共振干涉现象不明显的共振能区高能段，通过子群、子群参数的定义，将能群平均共振截面通式转换为子群形式，然后使用一系列不同背景截面下能群平均共振截面求解子群总截面、子群分截面以及子群概率； 步骤四:将核燃料每个共振核素的子群参数包括子群截面、子群概率和相应的多群数据耦合；共振能量段使用子群参数，快群和热群使用多群数据，制作子群数据-多群数据耦合系统； 步骤五:针对每个共振能群，使用子群总截面、子群分截面和子群概率简化与多群形式一致的子群二维任意几何特征线中子输运方程；针对共振能区各个能群独立的进行固定源的子群输运计算，得到子群中子能谱，通过子群中子能谱和子群参数归并得到核燃料的共振参数；针对共振核素的共振峰的范围和大小，将共振核素分为3-4类，计算过程中的迭代在类的基础上进行； 步骤六:使用共振能区高能段到共振能区低能段的群间散射，通过群间散射将共振能区高能段的计算和共振能区低能段的计算耦合；在共振能区高能段使用子群方法得到该能量段的共振参数后，计算共振能区高能段到共振能群低能段的散射源和裂变源； 步骤七:利用得到的散射源 与裂变源计算源项，在共振干涉现象强烈的共振能区低能区，利用小波展开方法，对连续能量中子能谱的能量变量进行小波展开，得到共振能区的连续能量能谱小波尺度函数展开系数方程组； 步骤八:使用二维任意几何特征线方法计算与多群中子输运方程形式一样的小波尺度函数展开系数方程，得到共振区低能段多群中子能谱；对计算的散射源项进行简化，在迭代中将共振能区高能段到共振能区低能段的群间散射以及共振能区低能段的群内散射在迭代计算中作为常量，最后，使用共振区低能段多群中子能谱归并连续总截面和分截面，得到共振能区低能段的核燃料共振参数。
2.根据权利要求1所述的一种获得核燃料组件共振参数的方法，其特征在于:步骤一所述的制作常用核素多群数据库的过程如下: 一、使用NJOY程序加工多群数据； 二、将多群数据格式转换为NECL多群数据库； 三、合成燃耗中所需伪核素； 四、制作批处理完成核素的添加、删除、替换功能。
3.根据权利要求1所述的一种获得核燃料组件共振参数的方法，其特征在于:步骤四所述的耦合过程如下 一、在NECL多群数据库的基础上，将该数据库按照共振核素与非共振核素分为共振核素多群数据库和非共振核素多群数据库；二、对共振核素多群数据库中的核素进行共振能群高能段的子群参数计算，得到子群截面和子群概率参数； 三、将共振核素的共 振能群高能段的子群参数与该共振核素快群热群的多群截面数据放置在一起，与非共振核素的多群截面数据一起，形成子群数据-多群数据耦合系统。
全文摘要
一种获得核燃料组件共振参数的方法，步骤如下制作常用核素多群数据库和共振峰随能量的变化曲线；将孤立峰和密集峰区分，孤立峰所在区间是共振干涉明显的共振能区低能段，密集峰所在区间是共振干涉不明显的共振能区高能段；在共振能区高能段采用子群方法求解子群总截面、子群分截面以及子群概率和相应的多群数据耦合；通过高能段到低能段的群间散射将其计算耦合，得到高能段的共振参数后，计算高能段到低能段的散射源和裂变源；在低能区用小波展开方法得到多群中子能谱；使用共振区低能段多群中子能谱归并连续总截面和分截面，得到低能段的核燃料共振参数；本发明能有效获得任意共振材料组成，任意几何，任意数量共振区核燃料组件的共振参数。
文档编号G06F17/50GK103177154SQ201310045509
公开日2013年6月26日 申请日期2013年2月5日 优先权日2013年2月5日
发明者吴宏春, 贺雷, 曹良志 申请人:西安交通大学