CN113704996A - 一种基于轴向通量展开准三维输运计算方法及系统 - Google Patents
一种基于轴向通量展开准三维输运计算方法及系统 Download PDFInfo
- Publication number
- CN113704996A CN113704996A CN202110993778.7A CN202110993778A CN113704996A CN 113704996 A CN113704996 A CN 113704996A CN 202110993778 A CN202110993778 A CN 202110993778A CN 113704996 A CN113704996 A CN 113704996A
- Authority
- CN
- China
- Prior art keywords
- flux
- axial
- dimensional
- source region
- flat source
- 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
Links
- 230000004907 flux Effects 0.000 title claims abstract description 146
- 238000004364 calculation method Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 61
- 238000004590 computer program Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 230000004992 fission Effects 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 description 7
- 238000010168 coupling process Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009257 reactivity Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种基于轴向通量展开准三维输运计算方法及系统,包括:建立三维中子输运方程,三维中子输运方程的输入项为先进反应堆数据;在轴向上,采用一阶差分形式进行通量展开,对出现的负通量进行置零操作,根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,根据轴向、径向泄漏项计算得到各平源区的标通量。本发明对于负通量进行置零,并通过置零后统计得到边界流计算泄漏项,通过泄漏项重新计算该平源区标通量。从而避免简单置零造成流与通量不匹配,影响计算收敛性。
Description
技术领域
本发明涉及核反应堆堆芯设计和反应堆物理数值计算领域,具体涉及一种基于轴向通量展开准三维输运计算方法及系统。
背景技术
作为核反应堆系统分析计算的基础,反应堆物理分析计算通过求解中子输运方程,获得堆芯反应性和全堆精细功率分布。为快速开展先进核动力堆芯研发,需要研发先进的高精度反应堆物理设计软件。为模拟复杂结构堆芯,国内外正广泛开展基于精确物理模型和精细几何建模的“一步法”反应堆物理计算方法研究。
一步法输运计算方法是指对于中子输运方程,在空间、角度上少做甚至不做近似,从而得到具有高精度高保真度的三维中子输运方程计算结果。对于一步法输运计算,主要有三维特征线方法、二维/一维耦合方法、轴向通量展开的准三维等三种方法。三维特征线方法,是指采用特征线方法直接对于三维方程进行三维特征线求解,这种方法避免了近似,理论上具有最高的计算精度。然而,对于三维特征线求解,计算量较大,在现有计算条件下应用困难。现有技术尚未实现基于三维特征线方法的全堆芯全燃耗的一步法输运计算。二维/一维耦合方法,是指对于三维中子输运方程沿轴向、径向分别进行积分,采用逐层求解的思想,将三维方程转化为二维方程与一维方程耦合形式,将三维问题分解为二维、一维问题,采用泄漏项进行耦合。该方法针对大型压水堆轴向非均匀性较弱的特点,进行了适当近似,降低了计算量。但是该方法引入泄漏项造成了稳定性差问题,目前主要应用于压水堆计算,难以应用于新型先进反应堆。轴向通量展开的准三维方法,理论上避免了二维/一维耦合方法的稳定性问题,同时相比直接三维特征线方法具有明显的效率优势。
对于轴向通量展开准三维方法,一阶近似的情况下引入了差分关系,造成可能出现负通量情况,负通量会影响右端源项,造成迭代过程的不稳定。
发明内容
本发明所要解决的技术问题是基于轴向通量展开准三维方法,引入了差分关系,造成可能出现负通量情况,负通量影响右端源项,造成迭代过程的不稳定,直接影响到轴向通量展开准三维方法数值计算的稳定性,目的在于提供一种基于轴向通量展开准三维输运计算方法及系统,解决了提高基于轴向通量展开准三维方法数值计算稳定性的问题。
本发明通过下述技术方案实现:
一种基于轴向通量展开准三维输运计算方法,包括:建立三维中子输运方程,所述三维中子输运方程的输入项为先进反应堆数据;在轴向上,采用一阶差分形式进行通量展开,对出现的负通量进行置零操作,根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,根据所述轴向、径向泄漏项计算得到各平源区的标通量。
由于当前对于轴向通量展开准三维方法,一阶近似的情况下引入了差分关系,造成可能出现负通量情况,负通量会影响右端源项,造成迭代过程的不稳定。本发明方法对于负通量进行置零,并通过置零后统计得到边界流计算泄漏项,通过泄漏项重新计算该平源区标通量。从而避免简单置零造成流与通量不匹配,影响计算收敛性。
进一步的,根据所述轴向、径向泄漏项计算得到各平源区的标通量,具体为:将轴向、径向泄漏项引入所述三维中子输运方程中,采用准三维特征线方法进行求解,得到各平源区的标通量。
进一步的,对出现的负通量进行置零操作,具体为:对出射面边界角通量为负的通量进行置零操作。
进一步的,根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,还包括:根据置零后的边界面通量,经过累加积分,得到各个平源区的轴向、径向泄漏项。
进一步的,从所述三维中子输运方程出发,在轴向上,采用一阶差分形式进行通量展开,引入差分关系进行通量展开,得到第L层最终形式的二维中子输运方程形式下式:
其中,m表示角度,g表示能群,x、y、z分别表示空间所在位置的x、y、z坐标,表示总截面,η表示幅角正弦,μ表示极角余弦,表示层L能群g方向m上的平均角通量,和表示L层上下面边界角通量,表示总源项,为裂变源与散射源之和,ΔzL表示L层层高,μ>0表示从下到上进行逐层扫描,μ<0表示从上到下进行扫描。
进一步的,对于出射面角通量直接置零,其次根据置零后的边界面通量计算各个平源区i的轴向、径向泄漏项,最后通过三维中子输运方程的简化式,计算得到各平源区的标通量,如下:
本发明的第二种实现方式,采用上述的任一种基于轴向通量展开准三维输运计算方法,包括:方程建立模块:用于建立三维中子输运方程,所述三维中子输运方程的输入项为先进反应堆数据;负通量置零模块:用于在轴向上,采用一阶差分形式进行通量展开,对出现的负通量进行置零操作;标通量计算模块:用于根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,根据所述轴向、径向泄漏项计算得到各平源区的标通量。
本发明的第三种实现方式,一种设备,其特征在于,所述设备包括:一个或多个处理器;存储器,用于存储一个或多个程序,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器执行上述的一种基于轴向通量展开准三维输运计算方法。
本发明的第四种实现方式,一种存储有计算机程序的计算机可读存储介质,其特征在于,该程序被处理器执行时实现上述的一种基于轴向通量展开准三维输运计算方法。
本发明与现有技术相比,具有如下的优点和有益效果:
本发明提出负通量置零处理方法,进一步提高了轴向通量展开准三维输运方法的计算稳定性,为准三维方法在数字化反应堆一步法输运计算的应用提供重要支撑。
附图说明
为了更清楚地说明本发明示例性实施方式的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。在附图中:
图1为轴向通量展开差分关系示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
本实施例1是一种应用于轴向通量展开准三维输运计算方法,处理准三维方法中引入差分关系进行通量展开造成的边界负通量,通过对负通量进行置零,从而提高该方法的数值计算稳定性。
在准三维输运计算基于特征线方法中将空间区域划分为若干小区域,即平源区,各个平源区近似认为中子通量及中子源项一致;平源区内的中子通量即为平源区标通量;通过泄漏项计算平源区标通量的过程即采用准三维特征线方法求解中子输运方程。简单置零过程仅将边界角通量置为0,并不处理标通量,因此破坏了原始中子输运方程求解过程中流与通量的一致性,本实施例1对于负通量进行置零,并通过置零后统计得到边界流计算泄漏项,通过泄漏项重新计算该平源区标通量,从而避免简单置零造成流与通量不匹配影响计算收敛性。
对于轴向通量展开的准三维输运方法,从三维中子输运方程出发,三维中子输运方程的输入项为先进反应堆的CSG几何模块、截面数据,引入差分关系进行通量展开,得到第L层最终形式的二维中子输运方程形式如式(1)所示准三维输运计算基本公式,采用特征线方法对该式进行求解。实际计算中,在极角为正从下往上扫描时,利用L-边界角通量得到相应源项,通过二维特征线计算得到第L层平均角通量通过差分关系计算得到L+边界角通量如式(2)所示轴向差分关系式,在式(2)计算过程中会出现负通量问题。在差分网格层高较高的情况下,出射面边界角通量可能为负。差分网格即图1所示的层L,层高即Δz,图1右端所示的角通量(从下往上计算时)及(从上往下计算时),也即式(2)中的本实施例1首先对于出射面角通量直接置零,其次根据置零后的边界面通量(即和)计算各个平源区i的轴向、径向泄漏项和最后通过三维中子输运方程,进行求解,如式(3)所示。计算得到各平源区的标通量φg,i,如式(3)所示通量置零方法,根据置零后的边界面通量计算得到泄漏项再根据式(3)所示简化的三维中子输运方程求解通量,通过这种方式保证了通量与流的守恒关系。
本实施例1对于负通量进行置零后,通过累加积分得到轴向、径向泄漏项,通过中子平衡方程重新计算平源区标通量。从而提高轴向通量展开准三维输运方法计算稳定性,为准三维方法在数字化反应堆一步法输运计算的应用提供重要支撑。
本实施例1对于通量展开引入的差分关系,采用置零的方式处理潜在存在的负通量,并通过泄漏项计算标通量,从而保证标通量与净流的匹配。通过这种方式,提高了准三维输运方法的数值稳定性,对于轴向通量展开准三维方法的应用具有重要价值。其中,负通量为引入差分关系后计算得到的通量为负值,与实际物理不符,需额外处理,从而保证整个迭代过程的收敛;标通量即中子通量,与计算过程中产生的角通量对应,角通量与角度相关,标通量与角度无关;净流是中子输运方法中的通用名词,在边界上针对角通量在各个角度进行积分,即可计算得到该边界面的净流,再根据净流计算泄漏项。
实施例2
本实施例2是在实施例1的基础上,一种基于轴向通量展开准三维输运计算系统,系统采用实施例1的基于轴向通量展开准三维输运计算方法,实施例2的系统具体包括:
方程建立模块:用于建立三维中子输运方程,三维中子输运方程的输入项为先进反应堆的CSG几何模块、截面数据。
负通量置零模块:用于在轴向上,采用一阶差分形式进行通量展开,对出现的负通量进行置零操作;
标通量计算模块:用于根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,根据所述轴向、径向泄漏项计算得到各平源区的标通量。
实施例3
本实施例3在实施例1的基础上。本实施例3提供了一种设备,所述设备包括:一个或多个处理器;
存储器,用于存储一个或多个程序,当一个或多个程序被一个或多个处理器执行时,使得一个或多个处理器执行一种基于轴向通量展开准三维输运计算方法。
其中,一种基于轴向通量展开准三维输运计算方法按照实施例1中的方法步骤执行。在此不再一一赘述。
实施例4
本实施例4在实施例1基础上,本实施例4提供了一种存储有计算机程序的计算机可读存储介质,该程序被处理器执行时实现一种基于轴向通量展开准三维输运计算方法。
其中,一种基于轴向通量展开准三维输运计算方法按照实施例1中的方法步骤执行。在此不再一一赘述。
本领域普通技术人员可以理解实现上述事实和方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,涉及的程序或者所述的程序可以存储于一计算机所可读取存储介质中,该程序在执行时,包括如下步骤:此时引出相应的方法步骤,所述的存储介质可以是ROM/RAM、磁碟、光盘等等。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于轴向通量展开准三维输运计算方法,其特征在于,包括:
建立三维中子输运方程,所述三维中子输运方程的输入项为先进反应堆数据;
在轴向上,采用一阶差分形式进行通量展开,对出现的负通量进行置零操作,根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,根据所述轴向、径向泄漏项计算得到各平源区的标通量。
2.如权利要求1所述的方法,其特征在于,根据所述轴向、径向泄漏项计算得到各平源区的标通量,具体为:将轴向、径向泄漏项引入所述三维中子输运方程中,采用准三维特征线方法进行求解,得到各平源区的标通量。
3.如权利要求1所述的方法,其特征在于,对出现的负通量进行置零操作,具体为:对出射面边界角通量为负的通量进行置零操作。
4.如权利要求1所述的方法,其特征在于,根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,还包括:根据置零后的边界面通量,经过累加积分,得到各个平源区的轴向、径向泄漏项。
8.一种基于轴向通量展开准三维输运计算系统,其特征在于,采用权利要求1-7任一所述的方法,包括:
方程建立模块:用于建立三维中子输运方程,所述三维中子输运方程的输入项为先进反应堆数据;
负通量置零模块:用于在轴向上,采用一阶差分形式进行通量展开,对出现的负通量进行置零操作;
标通量计算模块:用于根据置零后的边界面通量计算各个平源区的轴向、径向泄漏项,根据所述轴向、径向泄漏项计算得到各平源区的标通量。
9.一种设备,其特征在于,所述设备包括:一个或多个处理器;存储器,用于存储一个或多个程序,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器执行如权利要求1-7任一所述的方法。
10.一种存储有计算机程序的计算机可读存储介质,其特征在于,该程序被处理器执行时实现如权利要求1-7任一所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110993778.7A CN113704996A (zh) | 2021-08-27 | 2021-08-27 | 一种基于轴向通量展开准三维输运计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110993778.7A CN113704996A (zh) | 2021-08-27 | 2021-08-27 | 一种基于轴向通量展开准三维输运计算方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113704996A true CN113704996A (zh) | 2021-11-26 |
Family
ID=78655784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110993778.7A Pending CN113704996A (zh) | 2021-08-27 | 2021-08-27 | 一种基于轴向通量展开准三维输运计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113704996A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114528719A (zh) * | 2022-04-22 | 2022-05-24 | 西安交通大学 | 一种基于二维堆芯的用于压水堆的在线能群压缩方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017193367A1 (zh) * | 2016-05-13 | 2017-11-16 | 清华大学 | 双频容性耦合等离子体刻蚀鞘层二维数值模拟方法及系统 |
CN112632822A (zh) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | 基于三维泄漏项分割的核反应堆中子通量获取方法及装置 |
-
2021
- 2021-08-27 CN CN202110993778.7A patent/CN113704996A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017193367A1 (zh) * | 2016-05-13 | 2017-11-16 | 清华大学 | 双频容性耦合等离子体刻蚀鞘层二维数值模拟方法及系统 |
CN112632822A (zh) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | 基于三维泄漏项分割的核反应堆中子通量获取方法及装置 |
Non-Patent Citations (2)
Title |
---|
赵晨等: "Analysis and comparison of the 2D/1D and quasi-3D methods with the direct transport code SHARK", 《NUCLEAR ENGINEERING AND TECHNOLOGY》, pages 21 * |
郑勇: "矩形特征线方法加速技术及三维中子输运计算方法研究", 《中国博士学位论文全文数据库 工程科技II辑》, pages 99 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114528719A (zh) * | 2022-04-22 | 2022-05-24 | 西安交通大学 | 一种基于二维堆芯的用于压水堆的在线能群压缩方法 |
CN114528719B (zh) * | 2022-04-22 | 2022-07-08 | 西安交通大学 | 一种基于二维堆芯的用于压水堆的在线能群压缩方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yamamoto et al. | Derivation of optimum polar angle quadrature set for the method of characteristics based on approximation error for the Bickley function | |
CN110580376B (zh) | 一种基于Revit和Dynamo的桩基础模型创建方法 | |
CN107273582A (zh) | 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法 | |
CN110717275B (zh) | 一种针对压水堆堆芯的三维中子通量数值模拟方法 | |
Yang et al. | MOCUM: A two-dimensional method of characteristics code based on constructive solid geometry and unstructured meshing for general geometries | |
CN113672849A (zh) | 一种基于轴向通量展开的一步法输运计算方法及系统 | |
CN112800641B (zh) | 一种基于区域分解并行的广义粗网有限差分加速方法 | |
CN107038293B (zh) | 针对板状燃料的基于等效一维板模型的共振自屏计算方法 | |
CN107038294A (zh) | 针对轻水堆的基于等效一维棒模型的共振自屏计算方法 | |
CN114387421A (zh) | 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法 | |
CN113704996A (zh) | 一种基于轴向通量展开准三维输运计算方法及系统 | |
CN114048271A (zh) | 北斗网格数据模型在数据库中的存储方法和装置 | |
CN112632822B (zh) | 基于三维泄漏项分割的核反应堆中子通量获取方法及装置 | |
Häggblom et al. | Transmission probability method of integral neutron transport calculation for two-dimensional rectangular cells | |
CN105512387B (zh) | 一种燃料组件精细功率重构方法 | |
Sanchez et al. | Diffusion piecewise homogenization via flux discontinuity ratios | |
Yamamoto et al. | Recent developments in the GENESIS code based on the Legendre polynomial expansion of angular flux method | |
CN116340709A (zh) | 一种六角形组件压水堆堆芯稳态中子扩散计算方法 | |
CN114723808B (zh) | 一种基于fpga的光条亚像素中心提取方法 | |
CN108694299A (zh) | 基于icem-cfd的二维有限元中子学稳态计算方法 | |
Connolly et al. | A coarse mesh radiation transport method for 2-D hexagonal geometry | |
CN115048811B (zh) | 基于松弛因子的三维中子输运方程计算方法及系统 | |
Kang et al. | The advanced multilevel predictor-corrector quasi-static method for pin-resolved neutron kinetics simulation | |
CN115048810A (zh) | 基于线性源加速的三维中子输运方程计算方法及系统 | |
CN106251401A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211126 |