一种基于反距离加权插值的核热强耦合方法技术

本发明专利技术公开了一种基于反距离加权插值的核热强耦合方法，步骤如下：1.根据核燃料组件几何结构构建网格模型；2.设置对应的边界条件；3.进行计算流体力学稳态计算，在达到当前残差收敛标准后，使用反距离加权插值完成网格模型间数据传递，开始进行中子物理计算；4.通过对应的温度均方根插值算法更新核素的温度数据库；5.中子物理计算完毕后，提取计算结果中的有效增殖因子k

【技术实现步骤摘要】
一种基于反距离加权插值的核热强耦合方法


[0001]本专利技术涉及核反应堆堆芯设计与安全分析领域，具体涉及一种基于反距离加权插值的核热强耦合方法。

技术介绍

[0002]反应堆多物理场耦合是对反应堆进行综合分析的一项关键技术，在传统方式下，将中子物理计算与热工流体计算分割开来，得到的结果通常是偏保守的，无法精确反映反应堆或堆内组件的真实物理情况与现象。因此将中子学计算程序与热工流体计算程序耦合起来进行核热耦合的稳态分析将为反应堆设计提供更加精细的物理参数与相关设计准则。
[0003]主流的蒙特卡罗方法输运计算程序主要有MCNP、OpenMC、Serpent以及SCALE；主流的计算流体力学(Computational Fluid Dynamics,CFD)软件主要有：ANSYS Fluent、ANSYS CFX以及STAR
‑
CCM+。目前核热耦合程序大多采用MCNP与Fluent、MCNP与STAR
‑
CCM+、MCNP与OpenFOAM进行耦合，其中二者的网格映射是一大难点，一般使用的方法为网格一一映射：即将MCNP与Fluent使用同一套网格，二者的网格数量、网格栅元的中心坐标与各个面的坐标均完全相同；但该方法缺点是在目标结构非常复杂时，会造成MCNP生成网格的困难，并且MCNP在输入卡上有最大栅元数目与最大曲面数目的限制，所以在流体力学计算需要较大数目网格量时，为了保证MCNP的成功运行，使用网格一一映射方法不得不将流体力学计算的网格数量进行下调，这样可能会造成流体力学计算网格质量较差，误差较大的情况。
[0004]还有一种较为常用的方法为体积权重方法，该方法较上面的方法有较大改进，中子学计算与流体力学计算不在使用同一套网格，而是在MCNP上使用粗网格，在流体力学计算上采用细网格，这样可以在一定程度上克服上述中受MCNP网格最大数量限制的因素，但该方法在求解体积与权重因子时会造成计算量大，求解困难的情况；并且该方法容易受到几何模型与网格划分方式的影响，若几何模型较为复杂或者MCNP与流体力学计算网格中有任意一方采用非结构化网格，则在精细耦合求解时会在求解体积上进行大量的运算，由此而造成一定的计算困难。
[0005]反距离加权插值算法也可称为距离倒数乘算法。其计算结果不依赖于相关的数学条件。目标节点的值是其与插值节点之间的距离和插值节点的值进行加权平均后计算得到的。当插值节点与目标节点相距较远时，证明该插值节点对于目标节点的影响较小，故该权值较小；反之，当插值节点与目标节点相距较近时，则证明该插值节点对于目标节点的影响较大，故该权值较大。
[0006]在使用反距离加权插值的三维核热耦合方法时，关键点在于网格的中心点坐标，所以反距离加权插值可以使用粗细两种网格进行计算。其中ANSYS Fluent计算采用精细的细网格，MCNP采用较为粗糙的粗网格进行计算。这样可以有效提高MCNP的运行速度并且提高ANSYS Fluent计算结果的精度，并且突破MCNP的最大栅元数目的限制。在反距离加权插值算法进行的核热耦合的过程中，几何体的相对形状则是次要的，只需要保证粗细网格的原点与相对位置是相同的即可，该方法在处理复杂几何时会非常有效。
[0007]综上所述，反距离加权插值算法的核热耦合具有以下几点优势：
[0008]1.突破MCNP最大栅元数与最大曲面数限制；
[0009]2.可以划分更加精细的ANSYS Fluent网格模型来进行更加精确的流体力学计算；
[0010]3.提高MCNP运行速度；
[0011]4.反距离加权插值算法与几何体的形状无关，更加高效地处理复杂几何；
[0012]中国先进实验快堆的板型燃料组件的流道为闭式流道，组件内冷却剂流道间几乎不存在冷却剂质量和动量的交换，组件间的水隙也较小，板型燃料组件中冷却剂流速很高，传热性能也随之较大；其与相同体积下的棒束型燃料组件相比，中子通量密度更高，工作温度较低，燃耗较深，为了进行更加精细的求解板型燃料元件与组件中的中子物理与热工水力情况，开发了基于反距离加权插值的精细三维核热强耦合方法。

技术实现思路

[0013]为了克服上述现有技术存在的问题，本专利技术的目的在于提供一种基于反距离加权插值的核热强耦合方法，应用本专利技术的方法所得到的数据可以为反应堆堆内燃料元件、组件、堆芯等相关设计研究工作提供参考。
[0014]为了达到上述目的，本专利技术采用如下技术方案：
[0015]一种基于反距离加权插值的核热强耦合方法，包括如下步骤：
[0016]步骤一：使用前处理建模软件、计算流体力学软件以及蒙特卡洛方法输运计算程序进行计算对象的模型构造与网格划分；使用前处理建模软件构造核燃料组件几何模型，并将该几何模型导入计算流体力学软件中，再使用计算流体力学软件划分核燃料组件的精细结构化网格并进行流体力学计算；使用蒙特卡罗方法输运计算程序对核燃料组件进行建模，模型中包含燃料元件、冷却剂与包壳；通过计算流体力学软件中的用户自定义函数对应用于蒙特卡罗方法输运计算程序中的核燃料组件模型完成粗网格自动划分；其中，通过计算流体力学软件中的用户自定义函数对应用于蒙特卡罗方法输运计算程序中核燃料组件模型完成粗网格自动划分的具体步骤如下：
[0017]步骤1：向计算流体力学软件的用户自定义函数中指定核燃料组件几何模型的尺寸以及材料；
[0018]步骤2：指定在使用计算流体力学软件的用户自定义函数对核燃料组件几何模型进行网格划分时坐标原点位置；
[0019]步骤3：指定X方向、Y方向以及Z方向的网格划分数目，通过计算流体力学软件的用户自定义函数获取核燃料组件的网格中心点坐标、曲面位置坐标；
[0020]步骤4：依据步骤3中得到的网格中心点坐标以及曲面位置坐标，生成用于中子物理计算的蒙特卡洛输运计算程序的曲面输入卡；
[0021]步骤5：依据步骤4中生成的曲面输入卡，指定相对应的布尔运算规则，生成用于中子物理计算的蒙特卡洛方法输运计算程序的栅元输入卡；
[0022]步骤6：依据步骤5中生成的栅元输入卡，通过指定不同区域材料标号生成用于中子物理计算的蒙特卡洛方法输运计算程序的材料输入卡；
[0023]步骤7：指定本征值计算中模拟粒子数、迭代次数以及裂变源项分布生成本征值计算输入卡；
[0024]步骤8：指定裂变能量统计计数，生成裂变能量计数输入卡；
[0025]步骤9：对步骤4中生成的曲面输入卡、步骤5中生成的栅元输入卡、步骤6中生成的材料输入卡、步骤7中生成的本征值计算输入卡以及步骤8中生成的裂变能量计数输入卡进行合并，得到最终的的蒙特卡洛方法输运计算程序的输入卡；
[0026]步骤二：分别设定计算流体力学软件、蒙特卡洛方法输运计算程序中对于核燃料组件的边界条件；其中计算流体力学软件设置组件为速度入口边界条件以及压力为零的出口边界条件，组件四周为对称边界条件；蒙特卡洛方法输运计算程序设置组件四周为全反射边界条件，上本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于反距离加权插值的核热强耦合方法，其特征在于：使用有限容积法求解三维流体力学方程并使用蒙特卡罗方法求解三维中子输运方程；通过构建流体力学细网格模型以及蒙特卡罗方法输运计算粗网格模型，考虑冷却剂温度、燃料温度对中子物理的反馈作用，使用反距离加权插值算法完成粗、细网格模型之间的数据映射，经耦合迭代后实现精细的核热耦合计算；该方法能够为核反应堆设计、安全分析与运行提供高精度的核热耦合计算结果；该方法包含如下步骤：步骤一：使用前处理建模软件、计算流体力学软件以及蒙特卡洛方法输运计算程序进行计算对象的模型构造与网格划分；使用前处理建模软件构造核燃料组件几何模型，并将该几何模型导入计算流体力学软件中，再使用计算流体力学软件划分核燃料组件的精细结构化网格并进行流体力学计算；使用蒙特卡罗方法输运计算程序对核燃料组件进行建模，模型中包含燃料元件、冷却剂与包壳；通过计算流体力学软件中的用户自定义函数对应用于蒙特卡罗方法输运计算程序中的核燃料组件模型完成粗网格自动划分；其中，通过计算流体力学软件中的用户自定义函数对应用于蒙特卡罗方法输运计算程序中的核燃料组件模型完成粗网格自动划分的具体步骤如下：步骤1：向计算流体力学软件的用户自定义函数中指定核燃料组件几何模型的尺寸以及材料；步骤2：指定在使用计算流体力学软件的用户自定义函数对核燃料组件几何模型进行网格划分时坐标原点位置；步骤3：指定X方向、Y方向以及Z方向的网格划分数目，通过计算流体力学软件的用户自定义函数获取核燃料组件的网格中心点坐标、曲面位置坐标；步骤4：依据步骤3中得到的网格中心点坐标以及曲面位置坐标，生成用于中子物理计算的蒙特卡洛输运计算程序的曲面输入卡；步骤5：依据步骤4中生成的曲面输入卡，指定相对应的布尔运算规则，生成用于中子物理计算的蒙特卡洛方法输运计算程序的栅元输入卡；步骤6：依据步骤5中生成的栅元输入卡，通过指定不同区域材料标号生成用于中子物理计算的蒙特卡洛方法输运计算程序的材料输入卡；步骤7：指定本征值计算中模拟粒子数、迭代次数以及裂变源项分布生成本征值计算输入卡；步骤8：指定裂变能量统计计数，生成裂变能量计数输入卡；步骤9：对步骤4中生成的曲面输入卡、步骤5中生成的栅元输入卡、步骤6中生成的材料输入卡、步骤7中生成的本征值计算输入卡以及步骤8中生成的裂变能量计数输入卡进行合并，得到最终的的蒙特卡洛方法输运计算程序的输入卡；步骤二：分别设定计算流体力学软件、蒙特卡洛方法输运计算程序中对于核燃料组件的边界条件；步骤三：使用计算流体力学软件进行流体力学计算，待连续性残差降至1E
‑
3时即认为流体力学计算已经收敛，然后通过用户自定义函数提取流体力学计算结果中目标几何体各个网格的温度与密度，再利用反距离加权插值算法进行计算，将流体力学计算模型中的精细网格上的数据映射进蒙特卡洛方法输运计算程序的粗糙网格中；
反距离加权插值算法：如果d(x,x
i
)≥err，则：否则：u(x)＝u
i
其中：式中：i为插值节点的索引；d(x,x...

【专利技术属性】
技术研发人员：张大林，周星光，王心安，王成龙，田文喜，苏光辉，秋穗正，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：