CN116933553B - 数值反应堆中子学的非结构网格体积修正方法 - Google Patents

数值反应堆中子学的非结构网格体积修正方法 Download PDF

Info

Publication number
CN116933553B
CN116933553B CN202310960863.2A CN202310960863A CN116933553B CN 116933553 B CN116933553 B CN 116933553B CN 202310960863 A CN202310960863 A CN 202310960863A CN 116933553 B CN116933553 B CN 116933553B
Authority
CN
China
Prior art keywords
volume
response matrix
neutron
correction
volume correction
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
Application number
CN202310960863.2A
Other languages
English (en)
Other versions
CN116933553A (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.)
Shanghai Jiaotong University
Original Assignee
Shanghai Jiaotong 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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN202310960863.2A priority Critical patent/CN116933553B/zh
Publication of CN116933553A publication Critical patent/CN116933553A/zh
Application granted granted Critical
Publication of CN116933553B publication Critical patent/CN116933553B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

一种数值反应堆中子学的非结构网格体积修正方法,根据变分原理构造二阶偶对称中子输运方程的泛函形式,经变量展开和离散处理得到未经修正的系数矩阵,通过对系数矩阵进行体积修正后,利用变分原理得到体积修正后的节块响应矩阵方程,再采用裂变源迭代的方式对体积修正后的节块响应矩阵方程进行求解被。本发明能够在不引入额外计算代价的条件下,使计算结果对于问题的体积快速收敛,从而有效地提升计算精度和计算效率,减少了计算资源的使用。

Description

数值反应堆中子学的非结构网格体积修正方法
技术领域
本发明涉及的是一种反应堆数值模拟领域的技术,具体是一种数值反应堆中子学的非结构网格体积修正方法。
背景技术
在数值反应堆计算中,为模拟小型反应堆等具有强烈非均匀性的先进反应堆型,需要采用具有较强几何描述能力的先进中子学算法,如非结构网格变分节块法(UVNM)。然而,在非均匀的问题中往往具有曲边类型的几何,直接采用非结构网格对曲边几何进行描述往往无法做到非结构网格与真实几何体积的完全守恒。体积不守恒现象在非均匀中子输运问题中会导致较大误差,而单纯通过加密网格逼近实际几何的真实体积则会大幅提高问题求解所需的计算时间与计算资源,限制了UVNM在反应堆数值模拟中的实用性。
发明内容
本发明针对现有非结构网格变分节块法的网格体积在曲边几何中与真实体积不守恒以及无法在较少的网格数目下实现对复杂问题的精确模拟的不足,提出一种数值反应堆中子学的非结构网格体积修正方法,可以在不引入额外计算代价的条件下,使计算结果对于问题的体积快速收敛,从而有效地提升计算精度和计算效率,减少了计算资源的使用。
本发明是通过以下技术方案实现的:
本发明涉及一种数值反应堆中子学的非结构网格体积修正方法,根据变分原理构造二阶偶对称中子输运方程的泛函形式,经变量展开和离散处理得到未经修正的系数矩阵,通过对系数矩阵进行体积修正后,利用变分原理得到体积修正后的节块响应矩阵方程,再采用裂变源迭代的方式对体积修正后的节块响应矩阵方程进行求解,具体包括以下步骤:
步骤一:首先根据变分原理构造二阶偶对称中子输运方程的泛函形式:
其中:为中子的位置;/>为中子的运动方向;ψ为节块内部的偶阶中子角通量;χγ为定义在节块表面的奇阶中子角通量;φ(r)为中子标通量;Σx(r)为各类反应截面;S(r)为包含群间散射源和裂变源的源项。
步骤二:应用Ritz离散将二阶偶对称中子输运方程中的变量分别展开后对泛函进行离散,得到未经修正的系数矩阵,具体为:ψ(r,Ω)=fT(r)ψ(Ω),φ(r)=fT(r)φ,S(r)=fT(r)s,,其中:向量ψ(Ω),ωγ(Ω),φ以及s分别为ψ(r,Ω),Ω·nγχγ(r,Ω),φ(r)以及S(r)的展开矩;向量f(r)和hγ(r)分别为在节块内部以及表面采用Gram-Schmidt正交化得到的正交归一的基函数,未经修正的系数矩阵包括:i,j=x,y,z,Eγ(r)=∫γfhTdΓ,IV=∫VffTdV,其中:Ωi代表角度沿坐标轴方向的分量,fT和hT分别为向量fT(r)和/>的简写。
步骤三:对系数矩阵进行体积修正,具体包括:
3.1)根据实际的曲边几何边界,在所给出的非结构网格中选定需要修正的区域,并计算当前区域的体积修正因子其中:该待修正区域中非结构网格的总体积之和为VM,实际几何体积为VA
3.2)采用修正因子对与体积相关的系数矩阵进行修正,得到体积修正后的系数矩阵其中:AC和/>代表修正后的系数矩阵,Σt代表总反应截面,V代表某一特定非结构网格区域。
3.3)重复步骤3.1和3.2直至遍历问题域中所有待修正区域。
步骤四:利用变分原理令泛函的一阶变分为0,得到满足非结构网格体积与实际曲边几何的守恒关系的体积修正后的节块响应矩阵方程ψn=Hnq-Cnωn其中:下标n为某一特定的角度方向,n=1,...,N;/>以及/>分别为出射和入射方向的中子偏流矩;响应矩阵/>
步骤五:采用裂变源迭代的方式对体积修正后的节块响应矩阵方程进行求解,具体为:采用群内迭代收敛群内中子流密度,采用多群迭代收敛多群偶阶中子通量,采用裂变源迭代收敛裂变源,最终求解得到全局的偶阶中子通量密度。
本发明涉及一种实现上述方法的系统,包括:参数初始化单元、响矩阵构造单元、响应矩阵修正单元以及响应矩阵方程求解单元,其中:参数初始化单元根据用户提供的材料以及几何信息对参数进行初始化并将材料截面信息和非结构网格信息传递给响应矩阵构造单元;响应矩阵构造单元根据Ritz离散得到未经修正的系数矩阵A以及Eγ(r),并将未经修正的系数矩阵传递给响应矩阵修正单元;响应矩阵修正单元遍历所有修正区域得到体积修正因子,并计算修正后的系数矩阵AC再根据满足非结构网格体积与实际曲边几何的守恒关系得到体积修正方法的响应矩阵Hn、Cn、Bn、Rn;响应矩阵方程求解单元采用裂变源迭代的方式对体积修正方法的响应矩阵进行求解,最终得到全局的偶阶中子通量密度。
技术效果
本发明根据实际的几何边界对问题的非结构网格选定区域进行修正;遍历所有待修正区域计算区域的体积修正因子,通过所提出的方法对系数响应矩阵进行体积修正,实现非结构网格与实际曲边几何的体积守恒。与现有技术相比,本发明在不引入额外计算量的条件下通过对非结构网格体积进行修正来获得更高的计算精度,从而减少在反应堆中子学性能以及安全性分析等应用中的计算资源消耗,并提高计算效率。
附图说明
图1为本发明的方法流程图;
图2为实施例圆形区域体积修正前后网格形状示意图;
图中:a.初始的未修正网格,b.体积修正后的网格;
图3为单栅元模型不同网格密度示意图;
图中:a)64网格,b)192网格,c)448网格,d)696网格;
图4为单栅元模型应用体积修正算法前后计算结果随网格密度的变化趋势示意图。
具体实施方式
本实施例涉及一种数值反应堆中子学的非结构网格体积修正方法,应用于UO2单栅元模型中对其可行性以及性能进行测试,如图2所示,为四种不同网格密度下的UO2单栅元模型,边界条件为全反射。
本实施例计算中UVNM所采用的节块内部空间展开阶数为2阶;节块表面的空间展开阶数为1阶;SN的方位角方向展开阶数为30阶;极角方向展开为18阶;有效增殖因数keff、中子通量以及裂变源的收敛限设置分别为1.0E-06,1.0E-05,1.0E-05,所有算例均在10个MPI进程下并行执行。
单栅元模型共有两个区域:燃料区域和慢化剂区域,经过计算得到燃料区域和慢化剂区域的实际几何体积分别为0.916cm2以及0.672cm2;根据实际体积与所输入非结构网格的总体积计算修正因子,根据修正后的系数矩阵构造出响应矩阵方程;最后通过裂变源迭代求解得到keff与中子通量分布。
如图3所示,为采用与不采用体积修正方法在不同网格密度下的计算结果对比,其中的参考解1.32557为多群蒙特卡罗的计算结果。计算结果表明,在低网格密度下不采用修正的keff计算结果与参考解的误差为1,125pcm,而采用体积修正方法的keff计算结果与参考解的误差仅为-72pcm。在不采用体积修正方法时,网格数目需要加密到696才能达到与采用体积修正方法相当的精度,而此时696网格数所需的计算时间约为64网格数计算时间的12倍,即体积修正方法在该问题中实现了约12倍的加速比。
与现有技术相比,本方法通过对非结构网格区域的响应矩阵进行体积修正,从而在不引入额外计算代价的情况下显著提高计算精度,节约了数值反应堆的中子学计算代价。在UO2单栅元模型的64网格问题中的误差从1,125pcm减少至-72pcm,在相同精度下实现了约12倍的加速比。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。

Claims (3)

1.一种数值反应堆中子学的非结构网格体积修正方法,其特征在于,根据变分原理构造二阶偶对称中子输运方程的泛函形式,经变量展开和离散处理得到未经修正的系数矩阵,通过对系数矩阵进行体积修正后,利用变分原理得到体积修正后的节块响应矩阵方程,再采用裂变源迭代的方式对体积修正后的节块响应矩阵方程进行求解;
所述的二阶偶对称中子输运方程的泛函形式为:
其中:为中子的位置;/>为中子的运动方向;ψ为节块内部的偶阶中子角通量;χγ为定义在节块表面的奇阶中子角通量;φ(r)为中子标通量;Σx(r)为各类反应截面;S(r)为包含群间散射源和裂变源的源项;
所述的未经修正的系数矩阵,具体为:ψ(r,Ω)=fT(r)ψ(Ω),φ(r)=fT(r)φ,S(r)=fT(r)s,其中:向量ψ(Ω),ωγ(Ω),φ以及s分别为ψ(r,Ω),Ω·nγχγ(r,Ω),φ(r)以及S(r)的展开矩;向量f(r)和hγ(r)分别为在节块内部以及表面采用Gram-Schmidt正交化得到的正交归一的基函数,未经修正的系数矩阵包括:Eγ(r)=∫γfhTdΓ,IV=vVffTdV,其中:Ωi代表角度沿坐标轴方向的分量,fT和hT分别为向量fT(r)和/>的简写;
所述的体积修正,具体包括:
1)根据实际的曲边几何边界,在所给出的非结构网格中选定需要修正的区域,并计算当前区域的体积修正因子其中:该待修正区域中非结构网格的总体积之和为VM,实际几何体积为VA
2)采用修正因子对与体积相关的系数矩阵进行修正,得到体积修正后的系数矩阵
其中:AC代表修正后的系数矩阵,Σt代表总反应截面,V代表某一特定非结构网格区域;
3)重复步骤1和步骤2直至遍历问题域中所有待修正区域;
所述的体积修正后的节块响应矩阵方程ψn=Hnq-Cnωn其中:下标n为某一特定的角度方向,n=1,...,N;/>以及/>分别为出射和入射方向的中子偏流矩;响应矩阵/>Rn=(Gn+I)-1(Gn-I),
2.根据权利要求1所述的数值反应堆中子学的非结构网格体积修正方法,其特征是,所述的采用裂变源迭代的方式对体积修正后的节块响应矩阵方程进行求解,具体为:采用群内迭代收敛群内中子流密度,采用多群迭代收敛多群偶阶中子通量,采用裂变源迭代收敛裂变源,最终求解得到全局的偶阶中子通量密度。
3.一种实现权利要求1或2所述数值反应堆中子学的非结构网格体积修正方法的系统,其特征在于,包括:参数初始化单元、响矩阵构造单元、响应矩阵修正单元以及响应矩阵方程求解单元,其中:参数初始化单元根据用户提供的材料以及几何信息对参数进行初始化并将材料截面信息和非结构网格信息传递给响应矩阵构造单元;响应矩阵构造单元根据应用Ritz离散得到未经修正的系数矩阵A以及Eγ(r),并将未经修正的系数矩阵传递给响应矩阵修正单元;响应矩阵修正单元遍历所有修正区域得到体积修正因子,并计算修正后的系数矩阵AC再根据满足非结构网格体积与实际曲边几何的守恒关系得到体积修正方法的响应矩阵Hn、Cn、Bn、Rn;响应矩阵方程求解单元采用裂变源迭代的方式对体积修正方法的响应矩阵进行求解,最终得到全局的偶阶中子通量密度。
CN202310960863.2A 2023-08-02 2023-08-02 数值反应堆中子学的非结构网格体积修正方法 Active CN116933553B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310960863.2A CN116933553B (zh) 2023-08-02 2023-08-02 数值反应堆中子学的非结构网格体积修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310960863.2A CN116933553B (zh) 2023-08-02 2023-08-02 数值反应堆中子学的非结构网格体积修正方法

Publications (2)

Publication Number Publication Date
CN116933553A CN116933553A (zh) 2023-10-24
CN116933553B true CN116933553B (zh) 2024-02-13

Family

ID=88380545

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310960863.2A Active CN116933553B (zh) 2023-08-02 2023-08-02 数值反应堆中子学的非结构网格体积修正方法

Country Status (1)

Country Link
CN (1) CN116933553B (zh)

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145657A (zh) * 2017-04-27 2017-09-08 西安交通大学 针对反应堆中子扩散方程的非均匀几何变分节块方法
CN107145472A (zh) * 2017-04-27 2017-09-08 西安交通大学 针对非均匀几何变分节块方法的广义矩阵分离加速方法
CN107423542A (zh) * 2017-04-13 2017-12-01 西安交通大学 一种适用于逐棒计算的非均匀泄漏修正方法
CN107742050A (zh) * 2017-11-14 2018-02-27 中国空气动力研究与发展中心计算空气动力研究所 一种面向混合网格高精度阻力预测的熵修正方法
CN110275733A (zh) * 2019-06-27 2019-09-24 上海交通大学 基于有限体积法求解声子玻尔兹曼方程的gpu并行加速方法
CN110287450A (zh) * 2019-06-25 2019-09-27 上海交通大学 基于积分变分节块法的响应矩阵并行处理实现方法
CN110909511A (zh) * 2019-11-22 2020-03-24 电子科技大学 一种无曲面体积分的无粘低速绕流数值模拟方法
CN110929443A (zh) * 2019-12-02 2020-03-27 中国水利水电科学研究院 一种基于高精度地形概化的二维洪水模拟方法
CN110990926A (zh) * 2019-12-02 2020-04-10 中国水利水电科学研究院 一种基于面积修正率的城市地表建筑水动力学仿真方法
CN111783343A (zh) * 2020-07-08 2020-10-16 中国人民解放军国防科技大学 非结构网格方向模板选择方法与流场求解方法
CN112632821A (zh) * 2020-12-21 2021-04-09 中国核动力研究设计院 基于区域分解获取堆芯中子标通量的方法
CN113673116A (zh) * 2021-09-01 2021-11-19 上海交通大学 针对均匀几何变分节块法的三维准输运加速方法
CN114387421A (zh) * 2022-01-17 2022-04-22 哈尔滨工程大学 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法
CN114444413A (zh) * 2022-01-21 2022-05-06 西安交通大学 一种板状燃料堆芯亚通道级三维热工水力分析方法
CN114444326A (zh) * 2022-02-25 2022-05-06 上海交通大学 针对三维瞬态多群中子扩散方程的非结构网格变分节块方法
CN114997035A (zh) * 2022-06-20 2022-09-02 西北核技术研究所 一种基于体积和非计数区修正的全局减方差方法
CN115691708A (zh) * 2022-09-06 2023-02-03 上海交通大学 压水堆crud轴向线性截面变化的中子学模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220222033A1 (en) * 2021-01-11 2022-07-14 Rovi Guides, Inc. Customized volume control in small screen media players

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107423542A (zh) * 2017-04-13 2017-12-01 西安交通大学 一种适用于逐棒计算的非均匀泄漏修正方法
CN107145657A (zh) * 2017-04-27 2017-09-08 西安交通大学 针对反应堆中子扩散方程的非均匀几何变分节块方法
CN107145472A (zh) * 2017-04-27 2017-09-08 西安交通大学 针对非均匀几何变分节块方法的广义矩阵分离加速方法
CN107742050A (zh) * 2017-11-14 2018-02-27 中国空气动力研究与发展中心计算空气动力研究所 一种面向混合网格高精度阻力预测的熵修正方法
CN110287450A (zh) * 2019-06-25 2019-09-27 上海交通大学 基于积分变分节块法的响应矩阵并行处理实现方法
CN110275733A (zh) * 2019-06-27 2019-09-24 上海交通大学 基于有限体积法求解声子玻尔兹曼方程的gpu并行加速方法
CN110909511A (zh) * 2019-11-22 2020-03-24 电子科技大学 一种无曲面体积分的无粘低速绕流数值模拟方法
CN110990926A (zh) * 2019-12-02 2020-04-10 中国水利水电科学研究院 一种基于面积修正率的城市地表建筑水动力学仿真方法
CN110929443A (zh) * 2019-12-02 2020-03-27 中国水利水电科学研究院 一种基于高精度地形概化的二维洪水模拟方法
CN111783343A (zh) * 2020-07-08 2020-10-16 中国人民解放军国防科技大学 非结构网格方向模板选择方法与流场求解方法
CN112632821A (zh) * 2020-12-21 2021-04-09 中国核动力研究设计院 基于区域分解获取堆芯中子标通量的方法
CN113673116A (zh) * 2021-09-01 2021-11-19 上海交通大学 针对均匀几何变分节块法的三维准输运加速方法
CN114387421A (zh) * 2022-01-17 2022-04-22 哈尔滨工程大学 基于cmfd加速输运计算的球床反应堆堆芯参数获取方法
CN114444413A (zh) * 2022-01-21 2022-05-06 西安交通大学 一种板状燃料堆芯亚通道级三维热工水力分析方法
CN114444326A (zh) * 2022-02-25 2022-05-06 上海交通大学 针对三维瞬态多群中子扩散方程的非结构网格变分节块方法
CN114997035A (zh) * 2022-06-20 2022-09-02 西北核技术研究所 一种基于体积和非计数区修正的全局减方差方法
CN115691708A (zh) * 2022-09-06 2023-02-03 上海交通大学 压水堆crud轴向线性截面变化的中子学模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于非均匀变分节块法的pin-by-pin计算加速算法研究;张滕飞 等;原子能科学技术;第53卷(第7期);1160-1164 *
非结构变形网格和离散几何守恒律;刘君 等;航空学报;第37卷(第8期);2395-2407 *

Also Published As

Publication number Publication date
CN116933553A (zh) 2023-10-24

Similar Documents

Publication Publication Date Title
Cai et al. Adaptive radial-basis-function-based multifidelity metamodeling for expensive black-box problems
Loseille et al. Achievement of global second order mesh convergence for discontinuous flows with adapted unstructured meshes
Zhao et al. Reliability‐based topology optimization using stochastic response surface method with sparse grid design
Li et al. Orthogonal bases for polynomial regression with derivative information in uncertainty quantification
Uryū et al. Binary neutron stars: equilibrium models beyond spatial conformal flatness
Kuzmin et al. Algebraic flux correction II. Compressible Euler equations
CN114444326B (zh) 针对三维瞬态多群中子扩散方程的非结构网格变分节块方法
CN116933553B (zh) 数值反应堆中子学的非结构网格体积修正方法
Deluzet et al. A two field iterated asymptotic-preserving method for highly anisotropic elliptic equations
Duan et al. An efficient ADER discontinuous Galerkin scheme for directly solving Hamilton-Jacobi equation
Kaye et al. A posteriori error estimator for adaptive local basis functions to solve Kohn-Sham density functional theory
McDonald et al. High-order CENO finite-volume schemes for multi-block unstructured mesh
Munguía et al. Goal-Oriented Mesh Adaptation for Flows in Thermochemical Nonequilibrium
Du et al. Uniform convergence of a nonlinear energy-based multilevel quantization scheme
Frambati et al. Practical unstructured splines: Algorithms, multi-patch spline spaces, and some applications to numerical analysis
Wang et al. Delaunay graph based inverse distance weighting for fast dynamic meshing
Masiello et al. High Order Method of Characteristics for 2-D Unstructured Meshes
Takagi et al. High-Order Low-Dissipation Shock-Resolving TENO-THINC Schemes for Hyperbolic Conservation Laws
Frank et al. Adaptive finite element simulation of the time-dependent simplified PN equations
Brezillon et al. Development and application of a flexible and efficient environment for aerodynamic shape optimisation
He et al. Aerodynamic data fusion with a multi-fidelity surrogate modeling method
Cheng et al. Discontinuity-resolving shock-capturing schemes on unstructured grids
Huhtanen et al. Exclusion and inclusion regions for the eigenvalues of a normal matrix
Ferrer et al. A robust arbitrarily high order transport method of the characteristic type for unstructured tetrahedral grids
Davies et al. Numerical experiments using deflation with the HISQ action

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