CN114547952B - 一种基于深度学习-界面流耦合的压水堆中子输运方法 - Google Patents
一种基于深度学习-界面流耦合的压水堆中子输运方法 Download PDFInfo
- Publication number
- CN114547952B CN114547952B CN202210438438.2A CN202210438438A CN114547952B CN 114547952 B CN114547952 B CN 114547952B CN 202210438438 A CN202210438438 A CN 202210438438A CN 114547952 B CN114547952 B CN 114547952B
- Authority
- CN
- China
- Prior art keywords
- neutron
- emergent
- triangular
- pressurized water
- water reactor
- 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.)
- Active
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/25—Design optimisation, verification or simulation using particle-based methods
-
- 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/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- 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
Abstract
本发明公开了一种基于深度学习‑界面流耦合的压水堆中子输运方法,通过蒙特卡罗方法模拟局部简单几何问题的中子与材料俘获、裂变或散射后的中子迁移过程,利用蒙特卡罗方法的模拟结果和深度学习方法学习不同材料和几何单元的中子入射‑出射分布规律,基于界面流守恒建立简单几何单元之间的连接关系和中子输运计算过程。该方法颠覆了确定论方法中必须先共振计算获得多群截面、再进行复杂几何多群输运计算的传统过程,形成一种综合了蒙特卡罗方法、深度学习和确定论方法的新型中子输运计算方法。
Description
技术领域
本发明涉及核反应堆堆芯设计和安全领域,具体涉及一种基于深度学习-界面流耦合的压水堆中子输运方法。
背景技术
在目前的核能发展中,针对复杂几何压水堆的中子输运计算,蒙特卡罗方法具有很高的计算精度和几何适应性,然而大量的粒子抽样计算导致计算效率很低。确定论方法相对蒙特卡罗方法具有很高的计算效率,但现阶段的共振计算方法和输运计算方法均无法满足复杂几何压水堆的高精度高效率计算要求。深度学习方法具有很强的复杂规律的学习能力,可以利用深度学习方法学习蒙特卡罗方法中大量中子与压水堆燃料、慢化剂、结构材料碰撞的迁移规律,实现深度学习方法在复杂几何压水堆的高精度高效率输运计算。
发明内容
提出一种基于深度学习-界面流耦合的压水堆中子输运方法,通过蒙特卡罗方法模拟三角形网格基本单元内中子与材料的俘获、裂变或散射后的中子迁移过程,利用蒙特卡罗方法的模拟结果训练深度学习模型不同材料的中子入射-出射分布规律,基于界面流守恒建立压水堆剖分后三角形网格基本单元之间的连接关系和中子输运计算扫描过程,形成一种综合了蒙特卡罗方法、深度学习和界面流方法的新型中子输运计算方法。
为了实现上述目的,本发明采用了以下技术方案予以实施:
一种基于深度学习-界面流耦合的压水堆中子输运方法,包括如下步骤:
步骤1:读取所要模拟的压水堆堆芯的几何信息和材料信息;
步骤2:针对步骤1读取的几何信息和材料信息,将压水堆燃料栅元剖分为三角形网格,每个三角形网格基本单元内是均匀材料,并利用特征参数表征三角形网格基本单元:
(a)几何特征参数:三角形的三个边的边长;
(b)材料特征参数:包括燃料特征参数、慢化剂特征参数和结构材料特征参数,燃料特征参数为温度和富集度;慢化剂特征参数为温度和硼浓度;结构材料特征参数为温度;
步骤4:基于步骤3获得的样本数据,针对慢化剂和结构材料的三角形网格基本单元,建立并训练以入射中子能量、入射中子角度、三角形网格基本单元几何特征参数和材料特征参数为输入,以三个边的出射中子能量分布和出射中子角度分布为输出的慢化剂-结构材料全连接深度学习网络模型;基于步骤3获得的样本数据,针对燃料的三角形网格基本单元,建立以压水堆特征值、入射中子能量、入射中子角度、三角形网格基本单元几何特征参数和材料特征参数为输入,以三角形网格基本单元内总的裂变中子数、三个边的出射中子能量分布、出射中子角度分布为输出的燃料全连接深度学习网络模型;
步骤5:根据三角形网格基本单元的界面流连续以及压水堆燃料栅元剖分为三角形网格基本单元时自然形成的几何相邻关系,从压水堆边界任意一个三角形网格基本单元出发,按照几何相邻关系依次扫描,直至获得压水堆内所有三角形网格基本单元的扫描顺序;设定压水堆初始特征值、总的初始裂变源项、所有三角形网格基本单元的三条边的初始入射中子能量分布和入射中子角度分布;
步骤6:根据三角形网格基本单元的入射中子能量和入射中子角度,利用步骤4获得的慢化剂-结构材料全连接深度学习网络模型,获得三个边的出射中子能量分布和出射中子角度分布;扫描所有慢化剂和结构材料三角形网格基本单元的每个边的入射中子能量和入射中子角度,对获得的出射中子能量分布和出射中子角度分布分别加和,获得各个边的总出射中子能量分布和总出射中子角度分布;
步骤7:根据压水堆特征值、燃料三角形网格基本单元的入射中子能量和入射中子角度,利用步骤4获得的燃料全连接深度学习网络模型,获得该入射中子对应的总的裂变中子数、三个边的出射中子能量分布和出射中子角度分布;扫描燃料区域的所有三角形网格基本单元的每个边的入射中子能量和入射中子角度,对获得的裂变中子数加和、出射中子能量分布和出射中子角度分布分别加和,获得压水堆总的裂变中子数、各个边的总出射中子能量分布和总出射中子角度分布;
步骤8:根据公式(1)更新压水堆新的特征值和裂变中子数;根据三角形网格基本单元的界面流连续和连接关系,利用三角形基本单元的出射中子分布更新三角形边的入射中子值;
步骤9:判断特征值和总裂变中子数是否收敛,若相邻两次计算的特征值的相对偏差大于10-5或相邻两次计算的总裂变中子数的相对偏差大于10-5,则未收敛,重复步骤6至步骤8,重新进行计算,直到收敛,至此完成了整个压水堆的输运计算。
与现有技术相比,本发明有如下突出优点:
本发明不是针对整个压水堆进行深度学习网络模型的训练,而是利用深度学习方法进行三角形网格基本单元内部中子碰撞和迁移复杂规律的学习,克服了针对整个压水堆计算难以产生大量样本的问题,并保证获得的深度学习网络模型具有更高的几何普适性。同时,本发明通过三角形网格基本单元的深度学习网络模型与界面流方法的结合,实现了针对复杂反应堆的高效率高精度计算。
附图说明
图1是本发明方法流程图。
图2是将一个单栅元剖分为三角形网格的示意图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明。
本发明一种基于深度学习-界面流耦合的压水堆中子输运方法,具体步骤如图1所示,本发明通过蒙特卡罗方法模拟局部简单几何问题的中子与介质俘获、裂变或散射后的中子迁移过程,并在压水堆中具体到每一个燃料栅元,如图2所示,将燃料栅元剖分为三角形网格单元。利用蒙特卡罗方法的模拟结果和深度学习方法学习不同材料和三角形网格单元的中子入射-出射迁移规律,基于界面流守恒建立简单几何单元之间的连接关系和中子输运计算过程。具体步骤如下:
步骤1:读取所要模拟的压水堆堆芯的几何信息和材料信息;
步骤2:针对步骤1读取的几何信息和材料信息,将压水堆燃料栅元剖分为三角形网格,每个三角形网格基本单元内是均匀材料,并利用特征参数表征三角形网格基本单元:
(a)几何特征参数:三角形的三个边的边长;
(b)材料特征参数:包括燃料特征参数、慢化剂特征参数和结构材料特征参数,燃料特征参数为温度和富集度;慢化剂特征参数为温度和硼浓度;结构材料特征参数为温度;
步骤4:基于步骤3获得的样本数据,针对慢化剂和结构材料的三角形网格基本单元,建立并训练以入射中子能量、入射中子角度、三角形网格基本单元几何特征参数和材料特征参数为输入,以三个边的出射中子能量分布和出射中子角度分布为输出的慢化剂-结构材料全连接深度学习网络模型;基于步骤3获得的样本数据,针对燃料的三角形网格基本单元,建立以压水堆特征值、入射中子能量、入射中子角度、三角形网格基本单元几何特征参数和材料特征参数为输入,以三角形网格基本单元内总的裂变中子数、三个边的出射中子能量分布、出射中子角度分布为输出的燃料全连接深度学习网络模型;
步骤5:根据三角形网格基本单元的界面流连续以及压水堆燃料栅元剖分为三角形网格基本单元时自然形成的几何相邻关系,从压水堆边界任意一个三角形网格基本单元出发,按照几何相邻关系依次扫描,直至获得压水堆内所有三角形网格基本单元的扫描顺序;设定压水堆初始特征值、总的初始裂变源项、所有三角形网格基本单元的三条边的初始入射中子能量分布和入射中子角度分布;
步骤6:根据三角形网格基本单元的入射中子能量和入射中子角度,利用步骤4获得的慢化剂-结构材料全连接深度学习网络模型,获得三个边的出射中子能量分布和出射中子角度分布;扫描所有慢化剂和结构材料三角形网格基本单元的每个边的入射中子能量和入射中子角度,对获得的出射中子能量分布和出射中子角度分布分别加和,获得各个边的总出射中子能量分布和总出射中子角度分布;
步骤7:根据压水堆特征值、燃料三角形网格基本单元的入射中子能量和入射中子角度,利用步骤4获得的燃料全连接深度学习网络模型,获得该入射中子对应的总的裂变中子数、三个边的出射中子能量分布和出射中子角度分布;扫描燃料区域的所有三角形网格基本单元的每个边的入射中子能量和入射中子角度,对获得的裂变中子数加和、出射中子能量分布和出射中子角度分布分别加和,获得压水堆总的裂变中子数、各个边的总出射中子能量分布和总出射中子角度分布;
步骤8:根据公式(1)更新压水堆新的特征值和总裂变中子数;根据三角形网格基本单元的界面流连续和连接关系,利用三角形基本单元的出射中子分布更新三角形边的入射中子值;
步骤9:判断特征值和总裂变中子数是否收敛,若相邻两次计算的特征值的相对偏差大于10-5或相邻两次计算的总裂变中子数的相对偏差大于10-5,则未收敛,重复步骤6至步骤8,重新进行扫描计算,直到收敛,至此完成了整个压水堆的输运计算。
Claims (1)
1.一种基于深度学习-界面流耦合的压水堆中子输运方法,其特征在于:包括如下步骤:
步骤1:读取所要模拟的压水堆堆芯的几何信息和材料信息;
步骤2:针对步骤1读取的几何信息和材料信息,将压水堆燃料栅元剖分为三角形网格,每个三角形网格基本单元内是均匀材料,并利用特征参数表征三角形网格基本单元:
(a)几何特征参数:三角形的三个边的边长;
(b)材料特征参数:包括燃料特征参数、慢化剂特征参数和结构材料特征参数,燃料特征参数为温度和富集度;慢化剂特征参数为温度和硼浓度;结构材料特征参数为温度;
步骤4:基于步骤3获得的样本数据,针对慢化剂和结构材料的三角形网格基本单元,建立并训练以入射中子能量、入射中子角度、三角形网格基本单元几何特征参数和材料特征参数为输入,以三个边的出射中子能量分布和出射中子角度分布为输出的慢化剂-结构材料全连接深度学习网络模型;基于步骤3获得的样本数据,针对燃料的三角形网格基本单元,建立以压水堆特征值、入射中子能量、入射中子角度、三角形网格基本单元几何特征参数和材料特征参数为输入,以三角形网格基本单元内总的裂变中子数、三个边的出射中子能量分布、出射中子角度分布为输出的燃料全连接深度学习网络模型;
步骤5:根据三角形网格基本单元的界面流连续以及压水堆燃料栅元剖分为三角形网格基本单元时自然形成的几何相邻关系,从压水堆边界任意一个三角形网格基本单元出发,按照几何相邻关系依次扫描,直至获得压水堆内所有三角形网格基本单元的扫描顺序;设定压水堆初始特征值、总的初始裂变源项、所有三角形网格基本单元的三条边的初始入射中子能量分布和入射中子角度分布;
步骤6:根据三角形网格基本单元的入射中子能量和入射中子角度,利用步骤4获得的慢化剂-结构材料全连接深度学习网络模型,获得三个边的出射中子能量分布和出射中子角度分布;扫描所有慢化剂和结构材料三角形网格基本单元的每个边的入射中子能量和入射中子角度,对获得的出射中子能量分布和出射中子角度分布分别加和,获得各个边的总出射中子能量分布和总出射中子角度分布;
步骤7:根据压水堆特征值、燃料三角形网格基本单元的入射中子能量和入射中子角度,利用步骤4获得的燃料全连接深度学习网络模型,获得该入射中子对应的总的裂变中子数、三个边的出射中子能量分布和出射中子角度分布;扫描燃料区域的所有三角形网格基本单元的每个边的入射中子能量和入射中子角度,对获得的裂变中子数加和、出射中子能量分布和出射中子角度分布分别加和,获得压水堆总的裂变中子数、各个边的总出射中子能量分布和总出射中子角度分布;
步骤8:根据公式(1)更新压水堆新的特征值和裂变中子数;根据三角形网格基本单元的界面流连续和连接关系,利用三角形基本单元的出射中子分布更新三角形边的入射中子值;
步骤9:判断特征值和总裂变中子数是否收敛,若相邻两次计算的特征值的相对偏差大于10-5或相邻两次计算的总裂变中子数的相对偏差大于10-5,则未收敛,重复步骤6至步骤8,重新进行计算,直到收敛,至此完成整个压水堆的输运计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210438438.2A CN114547952B (zh) | 2022-04-25 | 2022-04-25 | 一种基于深度学习-界面流耦合的压水堆中子输运方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210438438.2A CN114547952B (zh) | 2022-04-25 | 2022-04-25 | 一种基于深度学习-界面流耦合的压水堆中子输运方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114547952A CN114547952A (zh) | 2022-05-27 |
CN114547952B true CN114547952B (zh) | 2022-07-08 |
Family
ID=81667191
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210438438.2A Active CN114547952B (zh) | 2022-04-25 | 2022-04-25 | 一种基于深度学习-界面流耦合的压水堆中子输运方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114547952B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115935770B (zh) * | 2022-12-16 | 2023-07-25 | 西安交通大学 | 针对核反应堆屏蔽设计的基于神经网络的中子输运方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5267276A (en) * | 1991-11-12 | 1993-11-30 | The University Of Chicago | Neutron transport analysis for nuclear reactor design |
WO2015160398A2 (en) * | 2014-01-16 | 2015-10-22 | University Of Florida Research Foundation, Inc. | System and method for analysis of fissionable materials by multispectral active neutron interrogation analysis |
CN106354946B (zh) * | 2016-08-30 | 2017-07-04 | 中国科学院合肥物质科学研究院 | 一种基于过渡区域的蒙特卡罗与确定论耦合粒子输运方法 |
CN108733903B (zh) * | 2018-05-08 | 2019-06-11 | 西安交通大学 | 一种按中子能量进行个性化处理的中子输运数值模拟方法 |
CN109992864B (zh) * | 2019-03-22 | 2020-02-04 | 成都理工大学 | 非常规双重介质储层体积压裂数值模拟及参数优化方法 |
CN111523233B (zh) * | 2020-04-23 | 2021-12-28 | 西安交通大学 | 一种针对三维压水堆堆芯的中子输运计算方法 |
-
2022
- 2022-04-25 CN CN202210438438.2A patent/CN114547952B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114547952A (zh) | 2022-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107273582B (zh) | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 | |
CN111414722B (zh) | 一种核反应堆堆芯物理与热工耦合的模拟方法 | |
KR100763723B1 (ko) | 원자로 노심을 위한 컴퓨터 모델링 방법 | |
CN111523233B (zh) | 一种针对三维压水堆堆芯的中子输运计算方法 | |
Chen et al. | ARES: a parallel discrete ordinates transport code for radiation shielding applications and reactor physics analysis | |
CN107038293B (zh) | 针对板状燃料的基于等效一维板模型的共振自屏计算方法 | |
CN114547952B (zh) | 一种基于深度学习-界面流耦合的压水堆中子输运方法 | |
CN114913936B (zh) | 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法 | |
Shaukat et al. | Multiobjective core reloading pattern optimization of PARR-1 using modified genetic algorithm coupled with Monte Carlo methods | |
Clarno et al. | Capturing the effects of unlike neighbors in single-assembly calculations | |
CN106202867B (zh) | 一种计算快中子堆组件轴向肿胀效应的方法 | |
Ortensi et al. | Deterministic modeling of the high temperature test reactor | |
Tang et al. | Development and verification of an MOC code employing assembly modular ray tracing and efficient acceleration techniques | |
Robert et al. | Evaluation of Serpent capabilities for hyperfidelity depletion of pebble bed cores | |
Jain et al. | DIAMOND–A Method of Characteristics neutron transport code using unstructured meshing | |
Collins et al. | Coupled Simulations of Fast-Spectrum Molten Salt Reactors | |
Forget et al. | A heterogeneous coarse mesh solution for the 2-D NEA C5G7 mox benchmark problem | |
George et al. | Development of a Full-core Reactivity Equivalence for FeCrAl Enhanced Accident Tolerant Fuel in BWRs | |
Oktavian et al. | Triangle-Based Quasi-Diffusion Nodal Method for Fast Reactor Analysis | |
Song et al. | Multi-objective optimization of method of characteristics parameters based on genetic algorithm | |
Singh et al. | Development of a Lattice Neutron Transport Code Ν-Trac for 1-D Slab, 2-D Square and 2-D Hexagonal Geometries | |
Li et al. | Analysis and improvement of global-local self-shielding calculation scheme for AIC control rods | |
Rooijen et al. | Analysis of a VVER-1000 in-core fuel management benchmark with DRAGON and DONJON | |
Nezondet et al. | Investigations on the iso-geometric analysis method applied to neutron transport calculations in two-and three-dimensional spatial domains | |
Connolly et al. | Extending MPACT to 2D Hexagonal Geometry |
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 |