CN115455669A - 高温气冷堆球床均匀化群常数的确定方法及系统 - Google Patents

高温气冷堆球床均匀化群常数的确定方法及系统 Download PDF

Info

Publication number
CN115455669A
CN115455669A CN202211030876.1A CN202211030876A CN115455669A CN 115455669 A CN115455669 A CN 115455669A CN 202211030876 A CN202211030876 A CN 202211030876A CN 115455669 A CN115455669 A CN 115455669A
Authority
CN
China
Prior art keywords
temperature gas
balls
batch
cooled reactor
neutron flux
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.)
Pending
Application number
CN202211030876.1A
Other languages
English (en)
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.)
Huaneng Nuclear Energy Technology Research Institute Co Ltd
Original Assignee
Huaneng Nuclear Energy Technology Research Institute Co Ltd
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 Huaneng Nuclear Energy Technology Research Institute Co Ltd filed Critical Huaneng Nuclear Energy Technology Research Institute Co Ltd
Priority to CN202211030876.1A priority Critical patent/CN115455669A/zh
Publication of CN115455669A publication Critical patent/CN115455669A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本申请提出了高温气冷堆球床均匀化群常数的确定方法及系统，所述方法包括：获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；根据选取的所述多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。本申请提出的技术方案，在不降低计算精度的情况下，提高了高温气冷堆球床均匀化群常数的计算速度，节约了计算资源。

Description

高温气冷堆球床均匀化群常数的确定方法及系统
技术领域
本申请涉及匀化群常数计算领域，尤其涉及高温气冷堆球床均匀化群常数的确定方法及系统。
背景技术
球床式高温气冷堆运行过程中，燃料球和石墨球在球床堆芯中自上而下地缓慢流动，小部分流出堆芯的球将被替换成新的燃料球装填到堆芯，而大部分流出堆芯的球将返回堆芯循环使用，因此堆芯由不同批次(不同燃料富集度、燃耗水平和温度分布)的球(燃料球和石墨球)组成。
球床式高温气冷堆的堆芯扩散计算所需要的混合球床均匀化群常数，需要利用球床区域内各批次球的多群中子通量分布和均匀化群常数归并计算获得，以往都需要对球床的每个区域计算各批次球的碰撞概率，然后求解碰撞概率方程组来得到各批次球的多群中子通量分布，需要花费大量的计算时间和计算资源。
发明内容
本申请提供的高温气冷堆球床均匀化群常数的确定方法及系统，以至少解决在高温气冷堆球床均匀化群常数时需花费大量的计算时间和计算资源的技术问题。
本申请第一方面实施例提出一种高温气冷堆球床均匀化群常数的确定方法，所述方法包括：
获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；
利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；
基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；
根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。
优选的，所述获取高温气冷堆球床堆芯中不同批次的球，包括：
根据高温气冷堆球床堆芯中燃料球的燃料富集度、燃耗水平、温度分布将所述燃料球划分为不同批次的球，同时将高温气冷堆球床堆芯中的石墨球划分为与燃料球不同的一个批次。
进一步的，各组合模型中各批次球的多群中子通量比的计算式如下：
Figure BDA0003817316910000021
式中，
Figure BDA0003817316910000022
为第k种组合模型的第i批次球第g能群的中子通量比，
Figure BDA0003817316910000023
为第k种组合模型的第i批次球第g能群的中子通量密度，i∈[1～I],I为第k种组合模型中批次的总数。
优选的，所述根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数，包括：
获取各组合模型对应的插值权重系数；
根据所述各组合模型对应的插值权重系数及所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域内各批次球的多群中子通量；
根据所述当前高温气冷堆球床区域内各批次球的多群中子通量当确定前高温气冷堆球床区域的均匀化群常数。
进一步的，所述当前高温气冷堆球床区域内各批次球的多群中子通量的计算式如下：
Figure BDA0003817316910000024
式中，φi,g为当前高温气冷堆球床区域内第i批次球第g能群的中子通量比，
Figure BDA0003817316910000025
为第k种组合模型的第i批次球第g能群的中子通量比，ak为第k个组合模型对应的插值权重系数，k∈[1～K]，K为组合模型的个数。
进一步的，所述当确定前高温气冷堆球床区域的均匀化群常数的计算式如下：
Figure BDA0003817316910000026
式中，
Figure BDA0003817316910000027
为当确定前高温气冷堆球床区域第g能群的均匀化宏观反应截面，其中x表示反应类型，包括总、裂变、吸收等反应，fi为第i批次球占区域的体积份额，Σi,x,g为第i批次球第g能群的宏观反应截面，
Figure BDA0003817316910000031
为第i批次球第g'能群中子散射至第g能群的宏观散射截面，φi,g′为当前高温气冷堆球床区域内第i批次球第g'能群的中子通量比，Σi,s,g′→g为第i批次球第g'能群的宏观散射截面。
本申请第二方面实施例提出一种高温气冷堆球床均匀化群常数的确定系统，包括：
获取模块，用于获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；
第一确定模块，用于利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；
选取模块，用于基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；
第二确定模块，用于根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。
优选的，所述获取高温气冷堆球床堆芯中不同批次的球，包括：
根据高温气冷堆球床堆芯中燃料球的燃料富集度、燃耗水平、温度分布将所述燃料球划分为不同批次的球，同时将高温气冷堆球床堆芯中的石墨球划分为与燃料球不同的一个批次。
优选的，所述各组合模型中各批次球的多群中子通量比的计算式如下：
Figure BDA0003817316910000032
式中，
Figure BDA0003817316910000033
为第k种组合模型的第i批次球第g能群的中子通量比，
Figure BDA0003817316910000034
为第k种组合模型的第i批次球第g能群的中子通量密度，i∈[1～I],I为第k种组合模型中批次的总数。
优选的，所述第二确定模块，包括：
获取单元，用于获取各组合模型对应的插值权重系数；
第一确定单元，用于根据所述各组合模型对应的插值权重系数及所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域内各批次球的多群中子通量；
第二确定单元，用于根据所述当前高温气冷堆球床区域内各批次球的多群中子通量当确定前高温气冷堆球床区域的均匀化群常数。
本申请的实施例提供的技术方案至少带来以下有益效果：
本申请提出了高温气冷堆球床均匀化群常数的确定方法及系统，所述方法包括获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。本申请提出的技术方案，在不降低计算精度的情况下，提高了高温气冷堆球床均匀化群常数的计算速度，节约了计算资源。
本申请附加的方面以及优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本申请的实践了解到。
附图说明
本申请上述的和/或附加的方面以及优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
图1为根据本申请一个实施例提供的一种高温气冷堆球床均匀化群常数的确定方法的流程图；
图2为根据本申请一个实施例提供的一种高温气冷堆球床均匀化群常数的确定系统的结构图；
图3为根据本申请一个实施例提供的第二确定模块的结构图。
具体实施方式
下面详细描述本申请的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本申请，而不能理解为对本申请的限制。
本申请提出的高温气冷堆球床均匀化群常数的确定方法及系统，所述方法包括获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。本申请提出的技术方案，在不降低计算精度的情况下，提高了高温气冷堆球床均匀化群常数的计算速度，节约了计算资源。
下面参考附图描述本申请实施例的高温气冷堆球床均匀化群常数的确定方法及系统。
实施例一
图1为根据本申请一个实施例提供的一种高温气冷堆球床均匀化群常数的确定方法的流程图，如图1所示，所述方法包括：
步骤1：获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；
其中，所述获取高温气冷堆球床堆芯中不同批次的球，包括：
根据高温气冷堆球床堆芯中燃料球的燃料富集度、燃耗水平、温度分布将所述燃料球划分为不同批次的球，同时将高温气冷堆球床堆芯中的石墨球划分为与燃料球不同的一个批次。
步骤2：利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；
需要说明的是，所述各组合模型中各批次球的多群中子通量比的计算式如下：
Figure BDA0003817316910000051
式中，
Figure BDA0003817316910000052
为第k种组合模型的第i批次球第g能群的中子通量比，
Figure BDA0003817316910000053
为第k种组合模型的第i批次球第g能群的中子通量密度，i∈[1～I],I为第k种组合模型中批次的总数。
示例的，将高温气冷堆球床堆芯中不同燃料富集度、不同燃耗水平、不同温度的燃料球，划分成不同批次的球，另外再把球床中的石墨球当作一种特殊批次的球；
以排列组合的方式，搭建一系列包含不同数量比例、不同批次球(燃料球和石墨球)的物理模型，采用蒙特卡罗物理计算程序对不同排列组合球的组合模型进行中子输运计算，获得每个组合中各批次球的多群中子通量密度，进而获得每个组合中各批次球的多群中子通量比，其中，所述计算程序可以为：NECP-MCX；
一个组合的数据包含该组合中各批次球的多群中子通量比，每个组合的数据均占据表格的一个位置，从而制成高温气冷堆球型元件组合的多群中子通量比例表。
步骤3：基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；
步骤4：根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。
在本公开实施例中，所述步骤4具体包括：
步骤4-1：获取各组合模型对应的插值权重系数；
步骤4-2：根据所述各组合模型对应的插值权重系数及所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域内各批次球的多群中子通量；
进一步的，所述当前高温气冷堆球床区域内各批次球的多群中子通量的计算式如下：
Figure BDA0003817316910000061
式中，φi,g为当前高温气冷堆球床区域内第i批次球第g能群的中子通量比，
Figure BDA0003817316910000062
为第k种组合模型的第i批次球第g能群的中子通量比，ak为第k个组合模型对应的插值权重系数，k∈[1～K]，K为组合模型的个数。
需要说明的是，
Figure BDA0003817316910000063
公式的成立需满足
Figure BDA0003817316910000064
式中，
Figure BDA0003817316910000065
ζi为球床堆芯区域内第i批次球所占的数量比，ni为球床堆芯区域内第i批次球的数量。
示例的，根据实际工况下球床堆芯区域各批次球的数量信息、温度分布和燃耗水平，从高温气冷堆球型元件组合的多群中子通量比例表选取若干个具有相同批次球的组合数据，并选取合适的插值权重系数进行插值计算，获得实际工况下混合球床某个区域内各批次球的多群中子通量比。
步骤4-3：根据所述当前高温气冷堆球床区域内各批次球的多群中子通量当确定前高温气冷堆球床区域的均匀化群常数。
进一步的，所述当确定前高温气冷堆球床区域的均匀化群常数的计算式如下：
Figure BDA0003817316910000066
式中，
Figure BDA0003817316910000067
为当确定前高温气冷堆球床区域第g能群的均匀化宏观反应截面，其中x表示反应类型，包括总、裂变、吸收等反应，fi为第i批次球占区域的体积份额，Σi,x,g为第i批次球第g能群的宏观反应截面，
Figure BDA0003817316910000071
为第i批次球第g'能群中子散射至第g能群的宏观散射截面，φi,g′为当前高温气冷堆球床区域内第i批次球第g'能群的中子通量比，Σi,s,g′→g为第i批次球第g'能群的宏观散射截面。
综上所述，本实施例提出的一种高温气冷堆球床均匀化群常数的确定方法，所述球床均匀化群常数只需要各批次球的多群相对中子通量密度，即多群中子通量比，因此预制高温气冷堆球型元件组合的多群中子通量比例表，后续每一轮燃料循环计算直接插值计算即可获得实际工况下各批次球的多群中子通量比，省略了复杂的碰撞概率计算流程和碰撞概率方程组求解流程，进而在不降低计算精度的情况下大幅减少了计算时间和计算资源。
实施例二
图2为根据本申请一个实施例提供的一种高温气冷堆球床均匀化群常数的确定系统，如图2所示，包括：
获取模块100，用于获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；
第一确定模块200，用于利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；
选取模块300，用于基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；
第二确定模块400，用于根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。
在本公开实施例中，所述获取高温气冷堆球床堆芯中不同批次的球，包括：
根据高温气冷堆球床堆芯中燃料球的燃料富集度、燃耗水平、温度分布将所述燃料球划分为不同批次的球，同时将高温气冷堆球床堆芯中的石墨球划分为与燃料球不同的一个批次。
在本公开实施例中，所述各组合模型中各批次球的多群中子通量比的计算式如下：
Figure BDA0003817316910000072
式中，
Figure BDA0003817316910000073
为第k种组合模型的第i批次球第g能群的中子通量比，
Figure BDA0003817316910000074
为第k种组合模型的第i批次球第g能群的中子通量密度，i∈[1～I],I为第k种组合模型中批次的总数。
在本公开实施例中，如图3所示，所述第二确定模块400，包括：
获取单元401，用于获取各组合模型对应的插值权重系数；
第一确定单元402，用于根据所述各组合模型对应的插值权重系数及所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域内各批次球的多群中子通量；
第二确定单元403，用于根据所述当前高温气冷堆球床区域内各批次球的多群中子通量当确定前高温气冷堆球床区域的均匀化群常数。
需要说明的是，所述当前高温气冷堆球床区域内各批次球的多群中子通量的计算式如下：
Figure BDA0003817316910000081
式中，φi,g为当前高温气冷堆球床区域内第i批次球第g能群的中子通量比，
Figure BDA0003817316910000082
为第k种组合模型的第i批次球第g能群的中子通量比，ak为第k个组合模型对应的插值权重系数，k∈[1～K]，K为组合模型的个数。
所述当确定前高温气冷堆球床区域的均匀化群常数的计算式如下：
Figure BDA0003817316910000083
式中，
Figure BDA0003817316910000084
为当确定前高温气冷堆球床区域第g能群的均匀化宏观反应截面，其中x表示反应类型，包括总、裂变、吸收等反应，fi为第i批次球占区域的体积份额，Σi,x,g为第i批次球第g能群的宏观反应截面，
Figure BDA0003817316910000085
为第i批次球第g'能群中子散射至第g能群的宏观散射截面，φi,g′为当前高温气冷堆球床区域内第i批次球第g'能群的中子通量比，Σi,s,g′→g为第i批次球第g'能群的宏观散射截面。
综上所述，本实施例提出的一种高温气冷堆球床均匀化群常数的确定系统，进而在不降低计算精度的情况下大幅减少了计算时间和计算资源。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本申请的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本申请的实施例所属技术领域的技术人员所理解。
尽管上面已经示出和描述了本申请的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本申请的限制，本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (10)

1.一种高温气冷堆球床均匀化群常数的确定方法，其特征在于，所述方法包括：
获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；
利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；
基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；
根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。
2.如权利要求1所述的方法，其特征在于，所述获取高温气冷堆球床堆芯中不同批次的球，包括：
根据高温气冷堆球床堆芯中燃料球的燃料富集度、燃耗水平、温度分布将所述燃料球划分为不同批次的球，同时将高温气冷堆球床堆芯中的石墨球划分为与燃料球不同的一个批次。
3.如权利要求1所述的方法，其特征在于，所述各组合模型中各批次球的多群中子通量比的计算式如下：
Figure FDA0003817316900000011
式中，
Figure FDA0003817316900000012
为第k种组合模型的第i批次球第g能群的中子通量比，
Figure FDA0003817316900000013
为第k种组合模型的第i批次球第g能群的中子通量密度，i∈[1～I],I为第k种组合模型中批次的总数。
4.如权利要求1所述的方法，其特征在于，所述根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数，包括：
获取各组合模型对应的插值权重系数；
根据所述各组合模型对应的插值权重系数及所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域内各批次球的多群中子通量；
根据所述当前高温气冷堆球床区域内各批次球的多群中子通量当确定前高温气冷堆球床区域的均匀化群常数。
5.如权利要求4所述的方法，其特征在于，所述当前高温气冷堆球床区域内各批次球的多群中子通量的计算式如下：
Figure FDA0003817316900000021
式中，φi,g为当前高温气冷堆球床区域内第i批次球第g能群的中子通量比，
Figure FDA0003817316900000022
为第k种组合模型的第i批次球第g能群的中子通量比，ak为第k个组合模型对应的插值权重系数，k∈[1～K]，K为组合模型的个数。
6.如权利要求5所述的方法，其特征在于，所述当确定前高温气冷堆球床区域的均匀化群常数的计算式如下：
Figure FDA0003817316900000023
式中，
Figure FDA0003817316900000024
为当确定前高温气冷堆球床区域第g能群的均匀化宏观反应截面，其中x表示反应类型，包括总、裂变、吸收等反应，fi为第i批次球占区域的体积份额，Σi,x,g为第i批次球第g能群的宏观反应截面，
Figure FDA0003817316900000025
为第i批次球第g'能群中子散射至第g能群的宏观散射截面，φi,g′为当前高温气冷堆球床区域内第i批次球第g'能群的中子通量比，
Figure FDA0003817316900000026
为第i批次球第g’能群的宏观散射截面。
7.一种高温气冷堆球床均匀化群常数的确定系统，其特征在于，包括：
获取模块，用于获取高温气冷堆球床堆芯中不同批次的球，然后将所述球进行排列组合得到不同批次组合的各组合模型；
第一确定模块，用于利用蒙特卡罗法确定各组合模型中各批次球的多群中子通量比，并将所述各组合模型中各批次球的多群中子通量比构成多群中子通量比例表；
选取模块，用于基于当前高温气冷堆球床堆芯中各批次球的数量信息、温度分布、燃耗水平在所述多群中子通量比例表中，选取当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比；
第二确定模块，用于根据所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域的均匀化群常数。
8.如权利要求7所述的系统，其特征在于，所述获取高温气冷堆球床堆芯中不同批次的球，包括：
根据高温气冷堆球床堆芯中燃料球的燃料富集度、燃耗水平、温度分布将所述燃料球划分为不同批次的球，同时将高温气冷堆球床堆芯中的石墨球划分为与燃料球不同的一个批次。
9.如权利要求7所述的系统，其特征在于，所述各组合模型中各批次球的多群中子通量比的计算式如下：
Figure FDA0003817316900000031
式中，
Figure FDA0003817316900000032
为第k种组合模型的第i批次球第g能群的中子通量比，
Figure FDA0003817316900000033
为第k种组合模型的第i批次球第g能群的中子通量密度，i∈[1～I],I为第k种组合模型中批次的总数。
10.如权利要求7所述的系统，其特征在于，所述第二确定模块，包括：
获取单元，用于获取各组合模型对应的插值权重系数；
第一确定单元，用于根据所述各组合模型对应的插值权重系数及所述当前高温气冷堆球床堆芯中各组合模型中各批次球的多群中子通量比确定当前高温气冷堆球床区域内各批次球的多群中子通量；
第二确定单元，用于根据所述当前高温气冷堆球床区域内各批次球的多群中子通量当确定前高温气冷堆球床区域的均匀化群常数。
CN202211030876.1A 2022-08-26 2022-08-26 高温气冷堆球床均匀化群常数的确定方法及系统 Pending CN115455669A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211030876.1A CN115455669A (zh) 2022-08-26 2022-08-26 高温气冷堆球床均匀化群常数的确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211030876.1A CN115455669A (zh) 2022-08-26 2022-08-26 高温气冷堆球床均匀化群常数的确定方法及系统

Publications (1)

Publication Number Publication Date
CN115455669A true CN115455669A (zh) 2022-12-09

Family

ID=84300385

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211030876.1A Pending CN115455669A (zh) 2022-08-26 2022-08-26 高温气冷堆球床均匀化群常数的确定方法及系统

Country Status (1)

Country Link
CN (1) CN115455669A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115983049A (zh) * 2023-03-20 2023-04-18 西安交通大学 一种应用于球床式高温气冷堆的不连续因子计算方法
CN117133490A (zh) * 2023-07-25 2023-11-28 华能核能技术研究院有限公司 一种缩短球床式高温气冷初装堆芯建立过程的方法及系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115983049A (zh) * 2023-03-20 2023-04-18 西安交通大学 一种应用于球床式高温气冷堆的不连续因子计算方法
CN115983049B (zh) * 2023-03-20 2023-10-31 西安交通大学 一种应用于球床式高温气冷堆的不连续因子计算方法
CN117133490A (zh) * 2023-07-25 2023-11-28 华能核能技术研究院有限公司 一种缩短球床式高温气冷初装堆芯建立过程的方法及系统
CN117133490B (zh) * 2023-07-25 2024-05-07 华能核能技术研究院有限公司 一种缩短球床式高温气冷初装堆芯建立过程的方法及系统

Similar Documents

Publication Publication Date Title
CN115455669A (zh) 高温气冷堆球床均匀化群常数的确定方法及系统
Rhodes et al. CASMO-5 development and applications
CN110580935B (zh) 一种获取全堆有效共振自屏截面的方法
Bostelmann et al. SCALE Analysis of a Fluoride Salt-Cooled High-Temperature Reactor in Support of Severe Accident Analysis
Guo et al. Optimization of reactivity control in a small modular sodium-cooled fast reactor
CN114203270A (zh) 一种适用于弥散颗粒燃料的子群参数计算方法及系统
Guo et al. Challenges and progress of uncertainty analysis for the pebble-bed high-temperature gas-cooled reactor
Yuan et al. Research on the on-the-fly homogenization method based on RMC code for criticality calculations
Wang et al. Evaluation of pebble-bed homogenized cross sections in HTGR fuel cycle simulations
Tran Neutronic characteristics of an OTTO refueling PBMR
CN115510703A (zh) 球床反应堆倒料的模拟计算方法
Ortiz et al. Nuclear fuel lattice optimization using neural networks and a fuzzy logic system
Nakano et al. Plutonium isotopic composition of high burnup spent fuel discharged from light water reactors
Reiss et al. Full-core SSCWR calculations applying a fast computational method
Wang et al. Effects of homogeneous geometry models in simulating the fuel balls in HTR-10
Rohde et al. Development of the coupled 3D neutron kinetics/thermal-hydraulics code DYN3D-HTR for the simulation of transients in block-type HTGR
CN116432381A (zh) 混合燃料球床平衡堆芯计算方法
Constantin et al. Void reactivity and pin power calculation for a CANDU cell using the SEU-43 fuel bundle
Oppe et al. Modeling of continuous reload HTR systems by the PANTHERMIX code system
CN117133490B (zh) 一种缩短球床式高温气冷初装堆芯建立过程的方法及系统
Reitsma et al. An overview of the FZJ-tools for HTR core design and reactor dynamics, the past, present and future
NAICKER et al. Uncertainty Propagation for the HTTR in Terms of the IAEA CRP on HTGR Uncertainty Analysis in Modelling
CN117409875A (zh) 一种球床式高温气冷堆反射层中子泄漏修正方法及系统
Hudson et al. A Pebble Bed Reactor cross section methodology
US20230197302A1 (en) Core performance calculation apparatus

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