CN114387421A - 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法 - Google Patents

基于cmfd加速输运计算的球床反应堆堆芯参数获取方法 Download PDF

Info

Publication number
CN114387421A
CN114387421A CN202210048455.5A CN202210048455A CN114387421A CN 114387421 A CN114387421 A CN 114387421A CN 202210048455 A CN202210048455 A CN 202210048455A CN 114387421 A CN114387421 A CN 114387421A
Authority
CN
China
Prior art keywords
grid
mesh
representing
central
fine
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210048455.5A
Other languages
English (en)
Other versions
CN114387421B (zh
Inventor
郝琛
文宇晨
岳鑫
康乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202210048455.5A priority Critical patent/CN114387421B/zh
Publication of CN114387421A publication Critical patent/CN114387421A/zh
Application granted granted Critical
Publication of CN114387421B publication Critical patent/CN114387421B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205—Re-meshing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明提供一种基于CMFD加速输运计算的球床反应堆堆芯参数获取方法。首先将球床区域和燃料球进行网格划分，分别得到CMFD计算的粗网和输运计算的细网，再通过距离关系使得粗网与细网能够严格对应并得到粗网的截面等信息。由于圆柱几何的特殊性，求得圆柱粗网下中心网格与邻居网格间的宽度和面积并从中子扩散方程开始推导，最终得到圆柱几何下CMFD线性系统。由于圆柱自身特殊的几何特性，用中子流密度代替中子流计算。并对球床网格进行二次划分的特殊处理，进而减少了求粗网均匀化参数的计算误差。为了保证高阶输运计算与CMFD方法的严格等价，引入了节块不连续因子和扩散系数修正因子，解决了传统CMFD可能计算发散的问题。

Description

基于CMFD加速输运计算的球床反应堆堆芯参数获取方法
技术领域
本发明涉及中子输运计算领域，具体涉及一种基于CMFD加速输运计算的球床反应堆堆芯参数获取方法。
背景技术
目前核反应堆建模与计算软件中OpenMOC、NECP-X、CBZ等软件都有常规轻水反应堆在高保真输运计算下的数值加速功能，其中最常用的数值加速方法是粗网有限差分方法CMFD。在传统的CMFD方法中，发展出基于广义等价理论的CMFD加速方法gCMFD，基于偏流的CMFD加速方法pCMFD等。
传统CMFD技术利用粗网格加速细网格，在中子输运计算迭代收敛之前，利用中子输运计算归并出粗网格的等效均匀化参数，再利用均匀化参数定义出来的一个等价的有限差分方程，通过求解该有限差分方程得到粗网格的交界面中子流和网格体积平均通量，从而修正中子输运计算的细网格参数，为精细化中子输运方程提供良好的初值，减少其迭代次数，进而加速中子输运计算的迭代收敛。
在传统CMFD基础上发展了基于偏流的pCMFD，其使用部分中子流而不是净中子流，从全局计算得到部分中子流更新每个局部问题边界的入射角通量。pCMFD在加速高阶计算的同时，可以作为耦合局部计算的全局框架。因为在传统CMFD和pCMFD中，仅仅只有一个节块不连续因子(NDF)来保证中子流的守恒，因此在gCMFD中增加一个扩散系数修正因子(MDF)来调整扩散系数，更好的保证了中子流的守恒，有利于计算的收敛。
MOC具有良好的几何适应性，用MOC方法求解球床堆的中子输运方程是极佳的方法，因此MTTA技术通过将许多条短特征线合并成长特征线进行输运计算，能够加速计算结果的收敛。
随着核电技术发展，新一代堆型-球床式高温气冷堆被提出，但是由于球床堆在球床区域复杂的几何结构，针对该几何结构，采用的现有的数值加速方法为：使用特征线方法(MOC)来求解球床堆中的中子输运方程，并采用并行加速，长特征线加速技术(MTTA)等方法来进行数值加速计算。
对于现有的数值加速方法，通常有以下几点缺陷：
(1)对不同的几何需要不同的网格划分方式；
(2)直接在输运计算过程中进行加速，加速效果差；
(3)仅仅只有一个不连续因子保证通量的正值，可能造成计算误差；
(4)对于球床堆中的复杂模型，难以寻找粗网与细网的对映关系；
(5)在光学厚度过大时，可能会导致结果收敛不佳。
发明内容
针对上述技术问题，本发明实施例提供一种基于CMFD加速输运计算的球床反应堆的堆芯参数获取方法，能够至少解决上述技术问题之一。
本发明采用的技术方案为提供一种基于CMFD加速输运计算的球床反应堆的堆芯参数获取方法，所述方法包括如下步骤：
S100，构建球床反应堆堆芯模型，所述球床反应堆堆芯模型包括球床模型和燃料球模型；
S110，对所述球床模型按照第一设定划分规则进行网格划分，得到球床的粗网格集C＝(C1，C2，…，Cn)；以及对所述燃料球模型按照第二设定划分规则进行划分，得到第一细网格集I1＝(I11，I12，…，I1m)；n为球床的粗网格集中包括的粗网格数量，m为第一细网格集中包括的细网格数量；
S120，对每个球床的粗网格按照第一设定划分规则进行网格划分，得到第二细网格集I2＝(I21，I22，…，I2k)，k为第二细网格集中包括的细网格数量；
S130，对于C中的任一粗网格Ci，i的取值为1到n，执行如下步骤：
遍历I1，如果Ci与Ij为相交或者包含关系，则Ci内被Ij占据的部分Ci∩Ij替换为燃料球的材料属性，j的取值为1到m；
S140，基于第一细网格集和第二细网格集中的细网截面信息进行中子输运计算，得到细网格中的中子通量密度、中子流和反应堆堆芯的有效增值系数；
S150，基于得到的细网格的中子通量密度，获取C中的粗网格均匀化截面信息；
S160，构建圆柱几何下的CMFD线性系统，并基于输运计算的中子流和获取的均匀化截面化信息对圆柱几何下的CMFD线性系统进行求解，得到粗网格的中子通量密度和CMFD方程的有效增值系数；
S170，基于得到的粗网格的中子通量密度和CMFD方程的有效增值系数进行中子输运计算，得到细网格中的中子通量密度和堆芯的有效增值系数；
S180，如果(K当-K上)/K当＜a，执行S190，否则，返回S150；其中，K当为当前得到的反应堆堆芯的有效增值系数，K上为上一次得到的反应堆的有效增值系数，a设定误差阈值；
S190，结束计算程序，输出当前得到的细网格的中子通量密度和反应堆的有效增值系数，作为目标堆芯参数；
所述圆柱几何下的CMFD线性系统为：
Figure BDA0003472757590000031
其中，Neighbors表示中心网格的邻居网格编号neighbors∈{W，E，B，T，S，N}，W，E，B，T，S，N分别表示中心网格的西侧、东侧、底侧、顶侧、南侧和北侧；
Figure BDA0003472757590000032
表示中心网格c向邻居网格n方向的耦合扩散系数；
Figure BDA0003472757590000033
表示中心网格c的平均中子通量；
Figure BDA0003472757590000034
表示邻居网格n向中心网格c方向的耦合扩散系数；
Figure BDA0003472757590000035
表示邻居网格n的平均中子通量；
Ac，n表示中心网格c与邻居网格n的交界面面积；
Vc表示中心网格c的体积；
∑s，g′→g，c表示中心网格c由g’能群到g能群的散射截面；
Figure BDA0003472757590000036
表示中心网格c在第g’能群的平均中子通量；
λ表示特征值；
χg表示第g能群的裂变能谱；
ν表示平均裂变中子数；
∑f，g’，c表示中心网格c第g’能群的裂变截面；
∑t，g，c表示中心网格c第g能群中的总宏观截面；
Figure BDA0003472757590000037
表示中心网格c在第g能群的平均中子通量。
本发明实施例提供的基于CMFD加速输运计算的球床反应堆的堆芯参数获取方法，首先将球床区域和燃料球进行网格划分，分别得到CMFD计算的粗网和输运计算的细网，再通过距离关系使得粗网与细网能够严格对应并得到粗网的截面等信息。由于圆柱几何的特殊性，求得圆柱粗网下中心网格与邻居网格间的宽度和面积并从中子扩散方程开始推导，最终得到圆柱几何下CMFD线性系统。由于圆柱自身特殊的几何特性，用中子流密度代替中子流计算。并对球床网格进行二次划分的特殊处理，进而减少了求粗网均匀化参数的计算误差。为了保证高阶输运计算与CMFD方法的严格等价，引入了节块不连续因子和扩散系数修正因子，解决了传统CMFD可能计算发散的问题。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的基于CMFD加速输运计算的球床反应堆的堆芯参数获取方法的流程示意图；
图2(a)和图2(b)为球床网格划分示意图；
图3(a)和图3(b)为燃料球网格划分示意图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1是本发明实施例提供的基于CMFD加速输运计算的球床反应堆的堆芯参数获取方法。本发明实施例提供的了一种基于CMFD加速输运计算的球床反应堆的堆芯参数获取方法，所述球床反应堆包括圆柱状的球床和设置在所述球床内部的燃料球，如图1所示，所述方法可包括如下步骤：
S100，构建球床反应堆堆芯模型，所述球床反应堆堆芯模型包括球床模型和燃料球模型。
在本发明实施例中，可采用CSG方法构建球床反应堆堆芯模型。CSG方法通过点、线、面的交集、并集和补集等逻辑操作对反应堆模型进行几何建模。具体地：基于输入信息构建构建球床反应堆堆芯模型，得到初始反应堆模型各部分的几何类型，如圆柱，立方体，球体等。输入信息可包括几何信息和材料信息，所述几何信息包括反应堆球床和燃料球的几何尺寸、几何空间布置，所述材料信息包括反应堆球床和燃料球各类材料的能群数和各个宏观截面信息。S110，对所述球床模型按照第一设定划分规则进行网格划分，得到球床的粗网格集C＝(C1，C2，…，Cn)；以及对所述燃料球模型按照第二设定划分规则进行划分，得到第一细网格集I1＝(I11，I12，…，I1m)；n为球床的粗网格集中包括的粗网格数量，m为第一细网格集中包括的细网格数量。
在本发明实施例中，第一设定划分规则为沿圆柱的径向、轴向和周向进行等距划分，即在r,θ,z三个方向上划分；第二设定划分规则为沿球体的半径、方位角和极角进行等径划分。其中，按照第一设定划分规则划分得到的球床粗网格可如图2(a)和图2(b)所示，图2(a)为圆柱几何的划分示意图，图2(b)为得到的球床粗网格示意图，每个粗网格可包括六个面，分别为东面、西面、顶面、底面、南面和北面，每个粗网格在径向上的长度、在轴向上的长度和圆心角均相同。按照第二设定划分规则划分的燃料球的细网格可如图3(a)和图3(b)所示，图3(a)为球形几何划分示意图，图3(b)为得到的燃料球细网格示意图。S120，对每个球床的粗网格按照第一设定划分规则进行网格划分，得到第二细网格集I2＝(I21，I22，…，I2k)，k为第二细网格集中包括的细网格数量。
S130，对于C中的任一粗网格Ci，i的取值为1到n，执行如下步骤：
遍历I1，如果Ci与Ij为相交或者包含关系，则Ci内被Ij占据的部分Ci∩Ij替换燃料球的材料属性，j的取值为1到m。
因为CMFD方程中的粗网所需要的截面信息需要由细网信息进行均匀化得到，所有这一步骤是为了准确描述粗网格和细网格产生相应关系时，能够将面上包含的截面信息进行交换，得到更为准确的粗网格信息。
具体地，Ci与Ij的关系可通过如下步骤确定：
S10，获取Ij的圆心到Ci的各个面的垂直距离djE，djw，djT，djB，djS，djN，其中，djE，djw，djT，djB，djS，djN分别表示Ij的圆心到Ci的东面、西面、顶面、底面、南面和北面的距离。
S20，如果min(djE，djw，djT，djB，djS，djN)＜Rj，则表示Ij与Ci相交；
S30，如果min(djE，djw，djT，djB，djS，djN)＞Rj，则表示Ij与Ci相离；
S20，如果
Figure BDA0003472757590000051
则表示Ij包含于Ci，p∈(E，W，T，B，S，N)。在本发明实施例中，W，E，B，T，S，N分别表示中心网格的西侧、东侧、底侧、顶侧、南侧和北侧，如图2(b)所示。
由于每个燃料球是按照等径划分的，因此，在Ij和Ci的关系为相交时，可基于Ij的圆心到Ci的最短距离确定Ci∩Ij。
S140，基于第一细网格集和第二细网格集中的细网截面信息进行中子输运计算，得到细网格中的中子通量密度、中子流和反应堆堆芯的有效增值系数。
本领域技术人员知晓，可采用现有方法将第一细网格集和第二细网格集中的细网截面信息作为构建的中子输运方程的输入参数进行中子输运计算，得到细网格中的中子通量密度、中子流和反应堆堆芯的有效增值系数。
S150，基于得到的细网格的中子通量密度，获取C中的粗网格均匀化截面信息。
在S150中，可通过如下方程获取粗网格均匀化参数：
Figure BDA0003472757590000061
Figure BDA0003472757590000062
Figure BDA0003472757590000063
Figure BDA0003472757590000064
Figure BDA0003472757590000065
其中：
φc,g表示粗网c内能群g中的粗网格中的中子通量；
φi,g表示细网i内能群g中的细网格中的中子通量；
∑t,c,g表示粗网c内能群g中的粗网总宏观截面；t表示总宏观截面；
∑t,i,g表示细网i内能群g中的细网总宏观截面；
v∑f，c，g表示能群g中粗网的宏观裂变截面；
v∑f，i，g表示能群g中细网的宏观裂变截面；
∑s，c，g′→g表示能群g’到g中粗网散射截面；
∑s，i，g′→g表示能群g’到g中细网散射截面；
χc，g表示能群g中粗网裂变能谱；
χi，g表示能群g中细网裂变能谱；
Δvi表示细网格的体积。
S160，构建圆柱几何下的CMFD线性系统，并基于输运计算的中子流和获取的均匀化截面化信息对圆柱几何下的CMFD线性系统进行求解，得到粗网格的中子通量密度和CMFD方程的有效增值系数。
在本发明实施例中，构建的圆柱几何下的CMFD线性系统可为：
Figure BDA0003472757590000071
其中，Neighbors表示中心网格的邻居网格编号neighbors∈{W，E，B，T，S，N}；
Figure BDA0003472757590000072
表示中心网格c向邻居网格n方向的耦合扩散系数；
Figure BDA0003472757590000073
表示中心网格c的平均中子通量；
Figure BDA0003472757590000074
表示邻居网格n向中心网格c方向的耦合扩散系数；
Figure BDA0003472757590000075
表示邻居网格n的平均中子通量；
Ac，n表示中心网格c与邻居网格n的交界面面积；
Vc表示中心网格c的体积；
∑s，g′→g，c表示中心网格c由g’能群到g能群的散射截面；
Figure BDA0003472757590000076
表示中心网格c在第g’能群的平均中子通量；
λ表示特征值；
χg表示第g能群的裂变能谱；
ν表示平均裂变中子数；
∑f，g’，c表示中心网格c第g’能群的裂变截面；
∑t，g，c表示中心网格c第g能群中的总宏观截面；
Figure BDA0003472757590000077
表示中心网格c在第g能群的平均中子通量。
其中，如图2(b)所示，Ac，n满足如下条件：
Figure BDA0003472757590000081
Figure BDA0003472757590000082
Ac，E＝ΔrcΔzc；
Ac，W＝ΔrcΔzc；
Ac，B＝rc′ΔrcΔθc；
Ac，T＝rc′ΔrcΔθc；
其中：
Δrc表示中心网格c径向的长度；
Δzc表示中心网格c轴向的长度；
Δθc表示中心网格c周向的长度；
rc′表示中心网格c的几何中心在径向坐标轴的值；
Figure BDA0003472757590000083
表示中心网格c在径向坐标轴左边的值；
Figure BDA0003472757590000084
表示中心网格c在径向坐标轴右边的值。
进一步地，在本发明实施例中，圆柱几何下的CMFD线性系统中的耦合耦合扩散系数可通过下式表示：
Figure BDA0003472757590000091
Figure BDA0003472757590000092
Figure BDA0003472757590000093
Figure BDA0003472757590000094
Figure BDA0003472757590000095
Figure BDA0003472757590000096
其中，jc，E、jc，W、jc，T、jc，B、jc，S和jc，N分别表示中心粗网格c向n方向的中子流，n∈{W，E，B，T，S，N}；
fc，n表示中心网格c向位于n方向的邻居网格的节块不连续因子(n∈{W，E，S，N，T，B})；
gc，n表示中心网格c向位于n方向的邻居网格的扩散系数修正因子；
Dc表示中心网格c的扩散系数；
Dn表示邻居网格n的扩散系数；
Δrc表示中心网格c径向的长度；
Δzc表示中心网格c轴向的长度；
ΔzB中心网格c在B方向的邻居网格的轴向长度；
ΔzT表示中心网格c在T方向的邻居网格的轴向长度；
Δθc表示中心网格c周向的长度
Δθn表示在E方向的邻居网格的周向长度；
Δrc表示中心网格c的径向长度；
ΔrE表示在E方向的邻居网格的径向长度；
rc′表示中心网格c的几何中心在径向坐标轴的值；
rW′表示位于中心网格c的W方向上的邻居网格的几何中心在径向坐标轴的值；
rE′表示位于中心网格c的E方向上的邻居网格的几何中心在径向坐标轴的值。
进一步地，在本发明实施例中，fc，n和gc，n满足下述条件：
如果
Figure BDA0003472757590000101
且
Figure BDA0003472757590000102
则
Figure BDA0003472757590000103
gc，n＝fc，n；否则，
Figure BDA0003472757590000104
gc，n＝1；
其中，
Figure BDA0003472757590000105
表示输运计算的细网格中的中心网格c向n方向的中子流；hc，n表示中心粗网格c和位于n方向的邻居网格的中心距离；
Figure BDA0003472757590000106
表示输运计算中的中心细网格c向位于n方向的邻居网格的交界面的中子通量；
Figure BDA0003472757590000107
表示输运计算中的中心细网格c的平均中子通量。
在传统轻水堆的方形粗网网格中，Δo总是恒定的，因此导致粗网的体积恒定。而在圆柱几何下，基于其几何特殊性，在θ方向的网格会呈现出逐渐增大的趋势，因此圆柱几何的中子流难以得到，故本发明实施例采用了中子流密度进行求解耦合扩散系数，以保证高阶中子输运计算的中子流与低阶CMFD线性系统的中子流一致。
进一步地，hc，n满足如下条件：
Figure BDA0003472757590000108
Figure BDA0003472757590000109
hc，E＝rc′Δθc，hE，c＝rc′ΔθE；
hc，W＝rc′Δθc，hW，c＝rc′ΔθW；
hc，B＝Δzc，hB，c＝ΔzB；
hc，T＝Δzc，hT，c＝ΔzT；
其中：
Δzc表示中心网格c轴向的长度；
ΔzB位于中心网格c的B方向上的邻居网格的轴向长度；
ΔzT表示位于中心网格c的T方向上的邻居网格的轴向长度；
Δθc表示中心网格c周向的长度；
ΔθE表示位于中心网格c的E方向上的邻居网格的周向长度；
ΔθW表示位于中心网格c的W方向上的邻居网格的周向长度；
Figure BDA0003472757590000111
表示中心网格c在径向坐标轴左边的值；
Figure BDA0003472757590000112
表示中心网格c在径向坐标轴右边的值；
rc′表示中心网格c的几何中心在径向坐标轴的值；
rS′表示位于中心网格c的S方向上的邻居网格的几何中心在径向坐标轴的值；
rN′表示位于中心网格c的N方向上的邻居网格的几何中心在径向坐标轴的值。
S170，基于得到的粗网格的中子通量密度和CMFD方程的有效增值系数进行中子输运计算，得到细网格中的中子通量密度、中子流和反应堆堆芯的有效增值系数。
将得到的粗网格的中子通量密度和CMFD方程的有效增值系数作为中子输运方程的输入参数，得到更新后的细网格中的中子通量密度、中子流和反应堆堆芯的有效增值系数。
S180，如果(K当-K上)/K当＜a，执行S190，否则，返回S150；其中，K当为当前得到的反应堆堆芯的有效增值系数，K上为上一次得到的反应堆的有效增值系数，a设定误差阈值。
在本发明一示意性实施例中，0＜a＜1，优选，a＝1/10-8。
S190，结束计算程序，输出当前得到的细网格的中子通量密度和反应堆的有效增值系数，作为目标堆芯参数。
本发明实施例中，使用的英文缩写对应的中文名称如下表所示：
Figure BDA0003472757590000113
Figure BDA0003472757590000121
本发明的实施例还提供了一种非瞬时性计算机可读存储介质，该存储介质可设置于电子设备之中以保存用于实现方法实施例中一种方法相关的至少一条指令或至少一段程序，该至少一条指令或该至少一段程序由该处理器加载并执行以实现上述实施例提供的方法。
本发明的实施例还提供了一种电子设备，包括处理器和前述的非瞬时性计算机可读存储介质。
本发明的实施例还提供一种计算机程序产品，其包括程序代码，当所述程序产品在电子设备上运行时，所述程序代码用于使该电子设备执行本说明书上述描述的根据本发明各种示例性实施方式的方法中的步骤。
虽然已经通过示例对本发明的一些特定实施例进行了详细说明，但是本领域的技术人员应该理解，以上示例仅是为了进行说明，而不是为了限制本发明的范围。本领域的技术人员还应理解，可以对实施例进行多种修改而不脱离本发明的范围和精神。本发明开的范围由所附权利要求来限定。

Claims (10)

1.一种基于CMFD加速输运计算的球床反应堆堆芯参数获取方法，其特征在于，所述球床反应堆堆芯包括圆柱状的球床和设置在所述球床内部的燃料球；所述方法包括如下步骤：
S100，构建球床反应堆堆芯模型，所述球床反应堆堆芯模型包括球床模型和燃料球模型；
S110，对所述球床模型按照第一设定划分规则进行网格划分，得到球床的粗网格集C＝(C1，C2，…，Cn)；以及对所述燃料球模型按照第二设定划分规则进行划分，得到第一细网格集I1＝(I11，I12，…，I1m)；n为球床的粗网格集中包括的粗网格数量，m为第一细网格集中包括的细网格数量；
S120，对每个球床的粗网格按照第一设定划分规则进行网格划分，得到第二细网格集I2＝(I21，I22，…，I2k)，k为第二细网格集中包括的细网格数量；
S130，对于C中的任一粗网格Ci，i的取值为1到n，执行如下步骤：
遍历I1，如果Ci与Ij为相交或者包含关系，则Ci内被Ij占据的部分Ci∩Ij替换为燃料球的材料属性，j的取值为1到m；
S140，基于第一细网格集和第二细网格集中的细网截面信息进行中子输运计算，得到细网格中的中子通量密度、中子流和反应堆堆芯的有效增值系数；
S150，基于得到的细网格的中子通量密度，获取C中的粗网格均匀化截面信息；
S160，构建圆柱几何下的CMFD线性系统，并基于输运计算的中子流和获取的均匀化截面化信息对圆柱几何下的CMFD线性系统进行求解，得到粗网格的中子通量密度和CMFD方程的有效增值系数；
S170，基于得到的粗网格的中子通量密度和CMFD方程的有效增值系数进行中子输运计算，得到细网格中的中子通量密度和堆芯的有效增值系数；
S180，如果(K当-K上)/K当＜a，执行S190，否则，返回S150；其中，K当为当前得到的反应堆堆芯的有效增值系数，K上为上一次得到的反应堆的有效增值系数，a设定误差阈值；
S190，结束计算程序，输出当前得到的细网格的中子通量密度和反应堆的有效增值系数，作为目标堆芯参数；
所述圆柱几何下的CMFD线性系统为：
Figure FDA0003472757580000021
其中，Neighbors表示中心网格的邻居网格编号neighbors∈{W，E，B，T，S，N}，W，E，B，T，S，N分别表示中心网格的西侧、东侧、底侧、顶侧、南侧和北侧；
Figure FDA0003472757580000022
表示中心网格c向邻居网格n方向的耦合扩散系数；
Figure FDA0003472757580000023
表示中心网格c的平均中子通量；
Figure FDA0003472757580000024
表示邻居网格n向中心网格c方向的耦合扩散系数；
Figure FDA0003472757580000025
表示邻居网格n的平均中子通量；
Ac,n表示中心网格c与邻居网格n的交界面面积；
Vc表示中心网格c的体积；
∑s,g′→g,c表示中心网格c由g’能群到g能群的散射截面；
Figure FDA0003472757580000026
表示中心网格c在第g’能群的平均中子通量；
λ表示特征值；
xg表示第g能群的裂变能谱；
v表示平均裂变中子数；
∑f,g′,c表示中心网格c第g’能群的裂变截面；
∑t,g,c表示中心网格c第g能群中的总宏观截面；
Figure FDA0003472757580000027
表示中心网格c在第g能群的平均中子通量。
2.根据权利要求1所述的方法，其特征在于，第一设定划分规则为沿圆柱的径向、轴向和周向进行划分；第二设定划分规则为沿球体的半径、方位角和极角进行划分。
3.根据权利要求2所述的方法，其特征在于，Ci与Ij的关系通过如下步骤确定：
S10，获取Ij的圆心到Ci的各个面的垂直距离djE，djw，djT，djB，djS，djN，其中，djE，djw，djT，djB，djS，djN分别表示Ij的圆心到Ci的东面、西面、顶面、底面、南面和北面的距离；
S20，如果min(djE，djw，djT，djB，djS，djN)＜Rj，则表示Ij与Ci相交；
S30，如果min(djE，djw，djT，djB，djS，djN)＞Rj，则表示Ij与Ci相离；
S20，如果
Figure FDA0003472757580000028
则表示Ij包含于Ci，p∈(E，W，T，B，S，N)。
4.根据权利要求1所述的方法，其特征在于，圆柱几何下的CMFD线性系统中的耦合耦合扩散系数通过下式表示：
Figure FDA0003472757580000031
Figure FDA0003472757580000032
Figure FDA0003472757580000033
Figure FDA0003472757580000034
Figure FDA0003472757580000035
Figure FDA0003472757580000036
其中，jc,E、jc,W、jc,T、jc,B、jc,S和jc,N分别表示中心粗网格c向n方向的中子流，n∈{W，E，B，T，S，N}；
fc,n表示中心网格c向位于n方向的邻居网格的节块不连续因子(n∈{W，E，S，N，T，B})；
gc,n表示中心网格c向位于n方向的邻居网格的扩散系数修正因子；
Dc表示中心网格c的扩散系数；
Dn表示邻居网格n的扩散系数；
Δrc表示中心网格c径向的长度；
Δzc表示中心网格c轴向的长度；
ΔzB中心网格c在B方向的邻居网格的轴向长度；
ΔzT表示中心网格c在T方向的邻居网格的轴向长度；
Δθc表示中心网格c周向的长度
Δθn表示在E方向的邻居网格的周向长度；
Δrc表示中心网格c的径向长度；
ΔrE表示在E方向的邻居网格的径向长度；
rc′表示中心网格c的几何中心在径向坐标轴的值；
rW′表示位于中心网格c的W方向上的邻居网格的几何中心在径向坐标轴的值；
rE′表示位于中心网格c的E方向上的邻居网格的几何中心在径向坐标轴的值。
5.根据权利要求4所述的方法，其特征在于，fc,n和gc,n满足下述条件：
如果
Figure FDA0003472757580000041
且
Figure FDA0003472757580000042
则
Figure FDA0003472757580000043
gc,n＝fc,n；否则，
Figure FDA0003472757580000044
gc,n＝1；
其中，
Figure FDA0003472757580000045
Figure FDA0003472757580000046
表示输运计算的细网格中的中心网格c向n方向的中子流；hc,n表示中心粗网格c和位于n方向的邻居网格的中心距离；
Figure FDA0003472757580000047
表示输运计算中的中心细网格c向位于n方向的邻居网格的交界面的中子通量；
Figure FDA0003472757580000048
表示输运计算中的中心细网格c的平均中子通量。
6.根据权利要求1所述的方法，其特征在于，Ac,n满足如下条件：
Figure FDA0003472757580000049
Figure FDA00034727575800000410
Ac,E＝ΔrcΔzc；
Ac,W＝ΔrcΔzc；
Ac,B＝rc′ΔrcΔθc；
Ac,T＝rc′ΔrcΔθc；
其中：
Δrc表示中心网格c径向的长度；
Δzc表示中心网格c轴向的长度；
Δθc表示中心网格c周向的长度；
rc′表示中心网格c的几何中心在径向坐标轴的值；
Figure FDA0003472757580000051
表示中心网格c在径向坐标轴左边的值；
Figure FDA0003472757580000052
表示中心网格c在径向坐标轴右边的值。
7.根据权利要求5所述的方法，其特征在于，hc,n满足如下条件：
Figure FDA0003472757580000053
Figure FDA0003472757580000054
hc,E＝rc′Δθc,hE,c＝rc′ΔθE；
hc,W＝rc′Δθc,hW,c＝rc′ΔθW；
hc,B＝Δzc,hB,c＝ΔzB；
hc,T＝Δzc,hT,c＝ΔzT；
其中：
Δzc表示中心网格c轴向的长度；
ΔzB位于中心网格c的B方向上的邻居网格的轴向长度；
ΔzT表示位于中心网格c的T方向上的邻居网格的轴向长度；
Δθc表示中心网格c周向的长度；
ΔθE表示位于中心网格c的E方向上的邻居网格的周向长度；
ΔθW表示位于中心网格c的W方向上的邻居网格的周向长度；
Figure FDA0003472757580000055
表示中心网格c在径向坐标轴左边的值；
Figure FDA0003472757580000056
表示中心网格c在径向坐标轴右边的值；
rc′表示中心网格c的几何中心在径向坐标轴的值；
rS′表示位于中心网格c的S方向上的邻居网格的几何中心在径向坐标轴的值；
rN′表示位于中心网格c的N方向上的邻居网格的几何中心在径向坐标轴的值。
8.根据权利要求1所述的方法，其特征在于，在S150中，通过如下方程获取粗网格均匀化参数：
Figure FDA0003472757580000061
Figure FDA0003472757580000062
Figure FDA0003472757580000063
Figure FDA0003472757580000064
Figure FDA0003472757580000065
其中：
φc,g表示粗网c内能群g中的粗网格中的中子通量；
φi,g表示细网i内能群g中的细网格中的中子通量；
∑t,c,g表示粗网c内能群g中的粗网总宏观截面；t表示总宏观截面；
∑t,i,g表示细网i内能群g中的细网总宏观截面；
v∑f,c,g表示能群g中粗网的宏观裂变截面；
v∑f,i,g表示能群g中细网的宏观裂变截面；
∑s,c,g′→g表示能群g’到g中粗网散射截面；
∑s,i,g′→g表示能群g’到g中细网散射截面；
χc,g表示能群g中粗网裂变能谱；
χi,g表示能群g中细网裂变能谱；
Δvi表示细网格的体积。
9.一种非瞬时性计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令或至少一段程序，所述至少一条指令或所述至少一段程序由处理器加载并执行以实现如权利要求1～7中任意一项的所述方法根据权利要求1所述的方法。
10.一种电子设备，其特征在于，包括处理器和权利要求8中的非瞬时性计算机可读存储介质。
CN202210048455.5A 2022-01-17 2022-01-17 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法 Active CN114387421B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210048455.5A CN114387421B (zh) 2022-01-17 2022-01-17 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210048455.5A CN114387421B (zh) 2022-01-17 2022-01-17 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法

Publications (2)

Publication Number Publication Date
CN114387421A true CN114387421A (zh) 2022-04-22
CN114387421B CN114387421B (zh) 2024-10-11

Family

ID=81201542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210048455.5A Active CN114387421B (zh) 2022-01-17 2022-01-17 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法

Country Status (1)

Country Link
CN (1) CN114387421B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114970293A (zh) * 2022-08-02 2022-08-30 西安交通大学 基于二维一维耦合傅里叶分析的粗网扩散系数计算方法
CN115048811A (zh) * 2022-07-12 2022-09-13 中国核动力研究设计院 基于松弛因子的三维中子输运方程计算方法及系统
CN115983049A (zh) * 2023-03-20 2023-04-18 西安交通大学 一种应用于球床式高温气冷堆的不连续因子计算方法
CN116933553A (zh) * 2023-08-02 2023-10-24 上海交通大学 数值反应堆中子学的非结构网格体积修正方法
CN115238481B (zh) * 2022-07-08 2024-01-09 中国核动力研究设计院 一种轴向分区控制棒组件截面计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103150424A (zh) * 2013-02-05 2013-06-12 西安交通大学 一种获取反应堆堆芯三维中子通量密度精细分布的方法
KR101439999B1 (ko) * 2013-04-08 2014-09-12 한국원자력연구원 3차원 원자로 노심 수송해석에 기반한 소격격자 노심 해석용 상수 생산장치 및 방법
CN110598304A (zh) * 2019-09-06 2019-12-20 西安交通大学 一种空间核电推进系统球床反应堆物理热工耦合分析方法
CN113094947A (zh) * 2021-03-29 2021-07-09 西安交通大学 一种核反应堆堆芯核热耦合分析方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103150424A (zh) * 2013-02-05 2013-06-12 西安交通大学 一种获取反应堆堆芯三维中子通量密度精细分布的方法
KR101439999B1 (ko) * 2013-04-08 2014-09-12 한국원자력연구원 3차원 원자로 노심 수송해석에 기반한 소격격자 노심 해석용 상수 생산장치 및 방법
CN110598304A (zh) * 2019-09-06 2019-12-20 西安交通大学 一种空间核电推进系统球床反应堆物理热工耦合分析方法
CN113094947A (zh) * 2021-03-29 2021-07-09 西安交通大学 一种核反应堆堆芯核热耦合分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
唐霄;李庆;柴晓明;涂晓兰;王侃;: "KYCORE程序中子输运计算的改进", 核动力工程, no. 2, 30 December 2017 (2017-12-30) *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115238481B (zh) * 2022-07-08 2024-01-09 中国核动力研究设计院 一种轴向分区控制棒组件截面计算方法
CN115048811A (zh) * 2022-07-12 2022-09-13 中国核动力研究设计院 基于松弛因子的三维中子输运方程计算方法及系统
CN115048811B (zh) * 2022-07-12 2023-11-24 中国核动力研究设计院 基于松弛因子的三维中子输运方程计算方法及系统
CN114970293A (zh) * 2022-08-02 2022-08-30 西安交通大学 基于二维一维耦合傅里叶分析的粗网扩散系数计算方法
CN114970293B (zh) * 2022-08-02 2022-10-11 西安交通大学 基于二维一维耦合傅里叶分析的粗网扩散系数计算方法
CN115983049A (zh) * 2023-03-20 2023-04-18 西安交通大学 一种应用于球床式高温气冷堆的不连续因子计算方法
CN115983049B (zh) * 2023-03-20 2023-10-31 西安交通大学 一种应用于球床式高温气冷堆的不连续因子计算方法
CN116933553A (zh) * 2023-08-02 2023-10-24 上海交通大学 数值反应堆中子学的非结构网格体积修正方法
CN116933553B (zh) * 2023-08-02 2024-02-13 上海交通大学 数值反应堆中子学的非结构网格体积修正方法

Also Published As

Publication number Publication date
CN114387421B (zh) 2024-10-11

Similar Documents

Publication Publication Date Title
CN114387421A (zh) 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法
KR100763723B1 (ko) 원자로 노심을 위한 컴퓨터 모델링 방법
Yang et al. MOCUM: A two-dimensional method of characteristics code based on constructive solid geometry and unstructured meshing for general geometries
CN109344459B (zh) 一种蒙特卡罗程序的可视建模与转换方法及系统
Cho et al. Solution of the C5G7MOX benchmark three-dimensional extension problems by the DeCART direct whole core calculation code
Robert et al. Evaluation of Serpent capabilities for hyperfidelity depletion of pebble bed cores
Mazumdar et al. Solution of the neutron transport equation by the Method of Characteristics using a linear representation of the source within a mesh
Herrero et al. 3D whole core fine mesh multigroup diffusion calculations by domain decomposition through alternate dissections
Filiciotto et al. Nonconforming Three-Dimensional Model for PWR Control Rod Movements Without Homogenization and Cusping Effect
Colameco Next generation iterative transport-diffusion methodology (ITDM), for LWR core analysis
Litskevich et al. Verification of the current coupling collision probability method with orthogonal flux expansion for the assembly calculations
CN114239437B (zh) 一种精细化瞬态中子输运计算的时间步加速方法
Guo et al. The MOCP code: method of characteristics code for pebble-bed HTR
CN113326648B (zh) 考虑环境效应的栅元均匀化常数计算方法、系统及终端
DeHart et al. Preliminary results for the OECD/NEA time dependent benchmark using Rattlesnake, Rattlesnake-IQS and TDKENO
Oktavian et al. Triangle-based Quasi-diffusion Nodal Method for Fast Reactor Analysis
CN114970293B (zh) 基于二维一维耦合傅里叶分析的粗网扩散系数计算方法
JPH10239480A (ja) 炉心性能計算方法および装置
Zerkak et al. LWR Multi-physics developments and Applications within the framework of the NURESIM European Project
CN113312792B (zh) 堆芯逐棒计算中栅元等效均匀化常数的计算方法及系统
Lu et al. Development and Validation of a Three‐Dimensional Diffusion Code Based on a High Order Nodal Expansion Method for Hexagonal‐z Geometry
JP3023185B2 (ja) 原子炉の炉心性能計算装置
Bekar et al. KENO-VI Primer: Performing Calculations Using SCALE’s Criticality Safety Analysis Sequence (CSAS6) with Fulcrum
CN113504993B (zh) 基于有效共振截面的控制棒尖齿效应处理方法及系统
CN113673006B (zh) 一种物理热工耦合可视化全堆芯复杂几何建模方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant