CN107145657A - 针对反应堆中子扩散方程的非均匀几何变分节块方法 - Google Patents
针对反应堆中子扩散方程的非均匀几何变分节块方法 Download PDFInfo
- Publication number
- CN107145657A CN107145657A CN201710287329.4A CN201710287329A CN107145657A CN 107145657 A CN107145657 A CN 107145657A CN 201710287329 A CN201710287329 A CN 201710287329A CN 107145657 A CN107145657 A CN 107145657A
- Authority
- CN
- China
- Prior art keywords
- locking nub
- formula
- neutron
- geometry
- vector
- 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
Links
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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
一种针对反应堆中子扩散方程的非均匀几何变分节块方法,包括如下步骤:1、建立包含节块内部的中子守恒关系以及节块表面的流连续性条件的泛函;2、利用x‑y方向的有限元形状函数和分片常量多项式对节块进行离散,实现非均匀节块的几何和材料的描述功能;3、推导得到响应矩阵方程;4、利用红‑黑扫描的方式对响应矩阵方程进行求解,最终得到整个非均匀求解区域的中子通量密度分布;本发明能够精细描述核反应堆的非均匀栅元结构;该方法将基于变分节块法,采用等参有限元来处理节块内部的精细几何结构,实现非均匀几何的中子扩散方程求解。
Description
技术领域
本发明针对核反应堆堆芯中子学计算领域,提出了一种针对反应堆中子扩散方程的非均匀几何变分节块方法。
背景技术
核反应堆中子学计算研究以核反应堆堆芯为应用对象,其堆芯由许多不同种类的组件构成。根据堆型的不同,组件内部的几何结构和材料布置复杂多变。因此,实际的反应堆中子学问题是一个三维非均匀几何的中子学问题。对核反应堆进行快速、精确的中子学计算,是反应堆设计和校核的基本保障。目前在堆芯物理设计过程中,主要采用均匀几何的扩散方程求解方法。
变分节块法最早由美国西北大学的E.E.Lewis教授提出,是中子学计算方法的杰出代表,具备扎实的工程应用背景。其主要应用作有美国阿贡国家实验室开发的VARIANT程序和NODAL程序,法国原子能委员会的ERANOS程序,以及美国爱达荷国家实验室的INSTANT程序等等。
变分节块法以二阶偶宇称形式的中子扩散方程为出发点,方程呈现椭圆方程的形式,有利于Garlerkin方法的应用,且更适合有限元方法的空间离散。变分节块法的计算思想是:首先通过变分方法在均匀几何求解区域建立包含二阶中子输运方程和自然边界条件的泛函;然后采用标准正交多项式进行Ritz离散,同时利用球谐函数实现角度展开,并构造响应矩阵;最后分别在三维堆芯的各个节块内分别求解响应矩阵方程;节块之间以流和其高阶矩耦合,最终得到问题区域的中子通量密度分布。变分节块法消除了横向积分,离散对象直接针对三维中子通量密度分布。因此,变分节块法不需要精细功率重构技术,只需将最终求得的中子通量密度矩代入展开式就可以得到中子通量密度分布。然而,传统的变分节块法程序仅具备均匀节块的处理能力,不足以精细描述节块内部的非均匀几何,无法避免节块均匀化带来的误差。
随着科学技术的发展和计算机水平的提高,人们已经越来越开始重视减少近似和假设,追求更高精度的中子学计算方法,消除均匀化过程是中子学发展的必然趋势。因此,针对反应堆中子扩散方程的非均匀几何变分节块方法对于中子学计算具有十分重要的意义。
发明内容
为了克服上述现有技术存在的问题,本发明的目的是提出一种针对反应堆中子扩散方程的非均匀几何变分节块方法,它能够精细描述核反应堆的非均匀栅元结构;该方法将基于变分节块法,采用等参有限元来处理节块内部的精细几何结构,实现非均匀几何的中子扩散方程求解。
为了实现上述目的,本发明采取了以下技术方案予以实施:
一种针对反应堆中子扩散方程的非均匀几何变分节块方法,步骤如下:
步骤1:首先根据公式(1)中二阶偶宇称扩散方程建立包含公式(3)中子通量密度φ和中子流密度j的泛函,泛函中包含节块内部的中子守恒关系以及节块表面的流连续性条件:
针对某一特定能群,在扩散近似下,二阶偶宇称扩散方程为:
式中:
φ—节块内部中子通量密度;
—对于x,y,z三个方向的偏导数算子;
Ω—方位角向量;
Σt—中子宏观总截面;
Σa—中子宏观吸收截面;
q—中子源项;
根据变分原理,在由若干节块组成的整个非均匀求解区域上,对应扩散方程的泛函写作各个节块内部及其表面上泛函的叠加贡献:
式中:
F[φ,j]—整个非均匀几何求解区域内的泛函;
Fv[φ,j]—单个节块内部的泛函;
v—节块的编号;
而扩散近似下的各节块泛函
其中Γ是外部边界;
步骤2:有限元形状函数g(x,y)能够用来描述曲边几何结构,因此利用x‐y方向的有限元形状函数g(x,y),分片常量多项式h(x,y),z方向的正交多项式f z(z)和f′z(z)对节块内部中子通量密度φ、节块表面中子流j分别展开,实现非均匀节块的几何和材料描述功能:
式中:
T—转置符号;
f z(z)—节块内部轴向标准正交多项式向量;
g(x,y)—x‐y方向有限元形状函数向量;
—克罗内克积即张量积;
φ—节块内部中子通量密度的展开矩向量;
f′z(z)—节块x‐y表面的轴向标准正交多项式向量;
f γ(γ′)—节块x‐y表面的径向标准正交多项式向量;
j±γ(γ′,z)—节块x‐y表面中子流密度展开矩向量,其中展开矩代表了展开系数的值;j±γ(γ′,z)是关于径向方向上的自变量γ′=x,y和轴向方向的自变量z的函数:当γ=x时γ′=y,j±x(y,z)代表节块左侧和右侧的表面中子流密度展开矩向量;当γ=y时γ′=x,j±y(x,z)代表节块下侧和上侧的表面中子流密度展开矩向量;
j±z(x,y)—节块z表面中子流密度展开矩向量,它是关于径向方向上的自变量x,y的函数;
Δz—节块z方向上的高度;
h(x,y)—节块z表面分片常量,对应于节块内部第e个面积为Ae的有限元网格,它满足其中δe为克罗内克常数,它在第e个有限元处为1,其他位置为0;有限元形状函数g(x,y)选取为等参有限元,能够精细描述压水堆栅元燃料棒的曲边几何形状;
步骤3:将公式(4)至公式(6)中的离散表达式代入各个节块内部的泛函公式(3),得到表征中子通量密度展开系数φ和中子流密度展开系数j之间关系的响应矩阵方程公式(8)、公式(10):
将公式(4)至公式(6)代入公式(3),得节块内部泛函的离散形式:
根据变分原理,对公式(7)取φ的一阶变分为0,得扩散近似形式的矩阵方程:
其中:
其中:
-1—矩阵的求逆;
φ—节块内部中子通量密度展开矩向量;
j—节块表面净中子流密度展开矩向量;
q—中子源项展开矩向量;
—响应矩阵,与节块内部的材料布置、几何形状有关;
对公式(7)取j的一阶变分为0,得节块表面偶宇称中子通量密度展开矩的连续性条件:
联立公式(8)和公式(9):
式中:
j—净中子流密度展开矩向量;
u—中子流源项展开矩向量;
其中:
为了将响应矩阵表达成通用形式,利用变量替换关系式
将公式(10)写为响应矩阵形式:
式中:
j +—出射中子流展开矩向量;
j -—入射中子流展开矩向量;
—响应矩阵,与节块内部的材料布置、几何形状有关,且有
式中:
—单位矩阵;
—响应矩阵,与节块内部的材料布置、几何形状有关;
步骤4:对公式(14)、公式(8)所代表的响应矩阵方程,利用红-黑迭代的方法进行迭代求解,最终得到整个非均匀几何求解区域的中子通量密度分布φ(r)和中子流密度分布j±γ(γ′,z)、j±z(x,y),从而完成针对反应堆中子扩散方程的非均匀几何变分节块方法。
与现有技术相比,本发明具有如下突出优点:
1.传统的节块方法绝大多数仅具备均匀节块的描述能力,本发明通过采用等参有限元对各节块内部的中子通量密度分布进行空间离散,可以有效地描述各个节块的内部的非均匀曲边几何结构,实现非均匀扩散计算的能力。
2.本发明通过在节块的轴向表面采用分片常量进行空间离散,能够实际描述非均匀节块表面的轴向泄漏分布,消除传统节块方法中轴向表面的均匀化过程,更加切合物理实际,同时提高计算精度。
附图说明
图1为等参有限元描述下的压水堆栅元非均匀几何。
具体实施方式
下面结合具体实施方式对本发明作进一步详细说明。该方法采用标准的源迭代的方法进行外迭代。针对群内迭代,具体计算流程包括以下几步:
步骤1:首先根据公式(1)中二阶偶宇称扩散方程建立包含公式(3)中子通量密度φ和中子流密度j的泛函,泛函中包含节块内部的中子守恒关系以及节块表面的流连续性条件:
针对某一特定能群,在扩散近似下,二阶偶宇称扩散方程为:
式中:
φ—节块内部中子通量密度;
—对于x,y,z三个方向的偏导数算子;
Ω—方位角向量;
Σt—中子宏观总截面;
Σa—中子宏观吸收截面;
q—中子源项;
根据变分原理,在由若干节块组成的整个非均匀求解区域上,对应扩散方程的泛函可以写作各个节块内部及其表面上泛函的叠加贡献:
式中:
F[φ,j]—整个非均匀几何求解区域内的泛函,其中φ代表非均匀几何求解区域的中子角通量密度,j代表非均匀几何求解区域的表面中子流密度;
Fv[φ,j]—单个节块内部的泛函,其中φ代表节块内部的中子角通量密度,j代表节块表面的中子流密度;
v—节块的编号;
而扩散近似下的各节块泛函
其中Γ是外部边界;
步骤2:有限元形状函数g(x,y)能够用来描述如图1所示的曲边几何结构,因此利用x‐y方向的有限元形状函数g(x,y),分片常量多项式h(x,y),z方向的正交多项式f z(z)和f′z(z)对节块内部中子通量密度φ、节块表面中子流j分别展开,实现非均匀节块的几何和材料描述功能:
式中:
T—转置符号;
f z(z)—节块内部轴向标准正交多项式向量;
g(x,y)—x‐y方向有限元形状函数向量;
—克罗内克积即张量积;
φ—节块内部中子通量密度的展开矩向量;
f′z(z)—节块x‐y表面的轴向标准正交多项式向量;
f γ(γ′)—节块x‐y表面的径向标准正交多项式向量;
j±γ(γ′,z)—节块x‐y表面的中子流密度展开矩向量,其中展开矩代表了展开系数的值;它是关于径向方向上的自变量γ′=x,y和轴向方向的自变量z的函数:当γ=x时γ′=y,j±x(y,z)代表节块左侧和右侧的表面中子流密度展开矩向量;当γ=y时γ′=x,j±y(x,z)代表节块下侧和上侧的表面中子流密度展开矩向量;
j±z(x,y)—节块z表面的中子流密度展开矩向量,它是关于径向方向上的自变量x,y的函数,脚标±z分别对应了节块的顶侧和低侧表面;
Δz—节块z方向上的高度;
h(x,y)—节块z表面分片常量,对应于节块内部第e个面积为Ae的有限元网格,它满足其中δe为克罗内克常数,它在第e个有限元处为1,其他位置为0;有限元形状函数g(x,y)选取为等参有限元,能够精细描述压水堆栅元燃料棒的曲边几何形状;
步骤3:将公式(4)至公式(6)中的离散表达式代入各个节块内部的泛函公式(3),得到表征中子通量密度展开系数φ和中子流密度展开系数j之间关系的响应矩阵方程公式(8)、公式(10):
将公式(4)至公式(6)代入公式(3),得节块内部泛函的离散形式:
根据变分原理,对公式(7)取φ的一阶变分为0,得扩散近似形式的矩阵方程:
其中:
其中:
-1—矩阵的求逆;
φ—节块内部中子通量密度展开矩向量;
j—节块表面净中子流密度展开矩向量;
q—中子源项展开矩向量;
—响应矩阵,与节块内部的材料布置、几何形状有关;
对公式(7)取j的一阶变分为0,得节块表面偶宇称中子通量密度展开矩的连续性条件:
联立公式(8)和公式(9):
式中:
j—净中子流密度展开矩向量;
u—中子流源项展开矩向量;
其中:
为了将响应矩阵表达成通用形式,利用变量替换关系式
将公式(10)写为响应矩阵形式:
式中:
j +—出射中子流展开矩向量;
j -—入射中子流展开矩向量;
—响应矩阵,与节块内部的材料布置、几何形状有关,且有
式中:
—单位矩阵;
—响应矩阵,与节块内部的材料布置、几何形状有关;
为使表达简便,省略展开基函数中的x,y,z自变量,公式(8)至公式(16)中的各系数矩阵可写作:
其中:
g|±x—g(x,y)函数在节块+x,-x表面上的取值;
g|±y—g(x,y)函数在节块+y,-y表面上的取值;
f z|±z—f z(z)函数在节块+z,-z表面上的取值;
i=x,y,z—关于x,y,z的一阶偏导数符号;
—单位矩阵,其阶数与fz对应;
其他符号的代表意义与前文相同。
步骤4:对于求解区域内各不同种类的节块,根据公式(11),公式(15)至公式(28)分别计算各节块的响应矩阵;
步骤5:针对特定能群,利用红‐黑扫描的方式对响应矩阵方程公式(14)进行迭代求解,得到出、入射偏中子流密度矩j +、j -;在求得出、入射偏中子流密度矩j +、j -后,代入公式(8),即可解得节块内部中子通量密度矩φ,从而可由公式(4)得知节块内部的中子通量密度分布φ(r);
步骤7:进行下一能群的计算,能够求解整个非均匀求解区域内各群的中子通量密度分布和有效增值因子,完成针对反应堆中子扩散方程的非均匀几何变分节块方法。
Claims (1)
1.一种针对反应堆中子扩散方程的非均匀几何变分节块方法,其特征在于:步骤如下:
步骤1:首先根据公式(1)中二阶偶宇称扩散方程建立包含公式(3)中子通量密度φ和中子流密度j的泛函,泛函中包含节块内部的中子守恒关系以及节块表面的流连续性条件:
针对某一特定能群,在扩散近似下,二阶偶宇称扩散方程为:
式中:
φ—节块内部中子通量密度;
—对于x,y,z三个方向的偏导数算子;
Ω—方位角向量;
Σt—中子宏观总截面;
Σa—中子宏观吸收截面;
q—中子源项;
根据变分原理,在由若干节块组成的整个非均匀求解区域上,对应扩散方程的泛函写作各个节块内部及其表面上泛函的叠加贡献:
式中:
F[φ,j]—整个非均匀几何求解区域内的泛函;
Fv[φ,j]—单个节块内部的泛函;
v—节块的编号;
而扩散近似下的各节块泛函
其中Γ是外部边界;
步骤2:有限元形状函数g(x,y)能够用来描述曲边几何结构,因此利用x-y方向的有限元形状函数g(x,y),分片常量多项式h(x,y),z方向的正交多项式f z(z)和f′z(z)对节块内部中子通量密度φ、节块表面中子流j分别展开,实现非均匀节块的几何和材料描述功能:
式中:
T—转置符号;
f z(z)—节块内部轴向标准正交多项式向量;
g(x,y)—x-y方向有限元形状函数向量;
—克罗内克积即张量积;
φ—节块内部中子通量密度的展开矩向量;
f′z(z)—节块x-y表面的轴向标准正交多项式向量;
f γ(γ′)—节块x-y表面的径向标准正交多项式向量;
j±γ(γ′,z)—节块x-y表面中子流密度展开矩向量,其中展开矩代表了展开系数的值;j±γ(γ′,z)是关于径向方向上的自变量γ′=x,y和轴向方向的自变量z的函数:当γ=x时γ′=y,j±x(y,z)代表节块左侧和右侧的表面中子流密度展开矩向量;当γ=y时γ′=x,j±y(x,z)代表节块下侧和上侧的表面中子流密度展开矩向量;
j±z(x,y)—节块z表面中子流密度展开矩向量,它是关于径向方向上的自变量x,y的函数;
Δz—节块z方向上的高度;
h(x,y)—节块z表面分片常量,对应于节块内部第e个面积为Ae的有限元网格,它满足其中δe为克罗内克常数,它在第e个有限元处为1,其他位置为0;有限元形状函数g(x,y)选取为等参有限元,能够精细描述压水堆栅元燃料棒的曲边几何形状;
步骤3:将公式(4)至公式(6)中的离散表达式代入各个节块内部的泛函公式(3),得到表征中子通量密度展开系数φ和中子流密度展开系数j之间关系的响应矩阵方程公式(8)、公式(10):
将公式(4)至公式(6)代入公式(3),得节块内部泛函的离散形式:
根据变分原理,对公式(7)取φ的一阶变分为0,得扩散近似形式的矩阵方程:
其中:
其中:
-1—矩阵的求逆;
φ—节块内部中子通量密度展开矩向量;
j—节块表面净中子流密度展开矩向量;
q—中子源项展开矩向量;
—响应矩阵,与节块内部的材料布置、几何形状有关;
对公式(7)取j的一阶变分为0,得节块表面偶宇称中子通量密度展开矩的连续性条件:
联立公式(8)和公式(9):
式中:
j—净中子流密度展开矩向量;
u—中子流源项展开矩向量;
其中:
为了将响应矩阵表达成通用形式,利用变量替换关系式
将公式(10)写为响应矩阵形式:
式中:
j +—出射中子流展开矩向量;
j -—入射中子流展开矩向量;
—响应矩阵,与节块内部的材料布置、几何形状有关,且有
式中:
—单位矩阵;
—响应矩阵,与节块内部的材料布置、几何形状有关;
步骤4:对公式(14)、公式(8)所代表的响应矩阵方程,利用红-黑迭代的方法进行迭代求解,最终得到整个非均匀几何求解区域的中子通量密度分布φ(r)和中子流密度分布j±γ(γ′,z)、j±z(x,y),从而完成针对反应堆中子扩散方程的非均匀几何变分节块方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710287329.4A CN107145657B (zh) | 2017-04-27 | 2017-04-27 | 针对反应堆中子扩散方程的非均匀几何变分节块方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710287329.4A CN107145657B (zh) | 2017-04-27 | 2017-04-27 | 针对反应堆中子扩散方程的非均匀几何变分节块方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107145657A true CN107145657A (zh) | 2017-09-08 |
CN107145657B CN107145657B (zh) | 2020-02-14 |
Family
ID=59775377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710287329.4A Active CN107145657B (zh) | 2017-04-27 | 2017-04-27 | 针对反应堆中子扩散方程的非均匀几何变分节块方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107145657B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108694299A (zh) * | 2018-07-17 | 2018-10-23 | 中国工程物理研究院核物理与化学研究所 | 基于icem-cfd的二维有限元中子学稳态计算方法 |
CN110287450A (zh) * | 2019-06-25 | 2019-09-27 | 上海交通大学 | 基于积分变分节块法的响应矩阵并行处理实现方法 |
CN112364294A (zh) * | 2020-10-27 | 2021-02-12 | 中国核动力研究设计院 | 堆芯数值计算矩阵的自动转换方法、装置、设备及介质 |
CN112613156A (zh) * | 2020-11-19 | 2021-04-06 | 中国核动力研究设计院 | 一种精细化燃料棒性能分析方法 |
CN112989595A (zh) * | 2021-03-09 | 2021-06-18 | 西安交通大学 | 一种压水堆堆芯瞬态精细功率重构方法 |
CN113254860A (zh) * | 2021-06-04 | 2021-08-13 | 南京航空航天大学 | 一种堆芯栅元中子通量的计算方法 |
CN113673116A (zh) * | 2021-09-01 | 2021-11-19 | 上海交通大学 | 针对均匀几何变分节块法的三维准输运加速方法 |
CN114444326A (zh) * | 2022-02-25 | 2022-05-06 | 上海交通大学 | 针对三维瞬态多群中子扩散方程的非结构网格变分节块方法 |
CN116933553A (zh) * | 2023-08-02 | 2023-10-24 | 上海交通大学 | 数值反应堆中子学的非结构网格体积修正方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5878100A (en) * | 1995-06-30 | 1999-03-02 | Abb Atom Ab | Fuel assembly for a boiling water reactor |
CN105426657A (zh) * | 2015-10-30 | 2016-03-23 | 西安交通大学 | 一种消除反应堆堆芯计算中控制棒尖齿效应的方法 |
CN105808503A (zh) * | 2016-03-07 | 2016-07-27 | 西安交通大学 | 针对反应堆逐棒计算中解析求解栅元不连续因子的方法 |
CN106202867A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种计算快中子堆组件轴向肿胀效应的方法 |
-
2017
- 2017-04-27 CN CN201710287329.4A patent/CN107145657B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5878100A (en) * | 1995-06-30 | 1999-03-02 | Abb Atom Ab | Fuel assembly for a boiling water reactor |
CN105426657A (zh) * | 2015-10-30 | 2016-03-23 | 西安交通大学 | 一种消除反应堆堆芯计算中控制棒尖齿效应的方法 |
CN105808503A (zh) * | 2016-03-07 | 2016-07-27 | 西安交通大学 | 针对反应堆逐棒计算中解析求解栅元不连续因子的方法 |
CN106202867A (zh) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | 一种计算快中子堆组件轴向肿胀效应的方法 |
Non-Patent Citations (2)
Title |
---|
曹良志等: ""先进的压水堆燃料管理计算方法研究及软件研制"", 《西安交通大学学报》 * |
王永平等: ""变分节块法求解中子扩散方程"", 《现代应用物理》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108694299B (zh) * | 2018-07-17 | 2022-03-18 | 中国工程物理研究院核物理与化学研究所 | 基于icem-cfd的二维有限元中子学稳态计算方法 |
CN108694299A (zh) * | 2018-07-17 | 2018-10-23 | 中国工程物理研究院核物理与化学研究所 | 基于icem-cfd的二维有限元中子学稳态计算方法 |
CN110287450A (zh) * | 2019-06-25 | 2019-09-27 | 上海交通大学 | 基于积分变分节块法的响应矩阵并行处理实现方法 |
CN110287450B (zh) * | 2019-06-25 | 2023-03-24 | 上海交通大学 | 基于积分变分节块法的响应矩阵并行处理实现方法 |
CN112364294A (zh) * | 2020-10-27 | 2021-02-12 | 中国核动力研究设计院 | 堆芯数值计算矩阵的自动转换方法、装置、设备及介质 |
CN112613156A (zh) * | 2020-11-19 | 2021-04-06 | 中国核动力研究设计院 | 一种精细化燃料棒性能分析方法 |
CN112989595B (zh) * | 2021-03-09 | 2022-12-09 | 西安交通大学 | 一种压水堆堆芯瞬态精细功率重构方法 |
CN112989595A (zh) * | 2021-03-09 | 2021-06-18 | 西安交通大学 | 一种压水堆堆芯瞬态精细功率重构方法 |
CN113254860A (zh) * | 2021-06-04 | 2021-08-13 | 南京航空航天大学 | 一种堆芯栅元中子通量的计算方法 |
CN113673116A (zh) * | 2021-09-01 | 2021-11-19 | 上海交通大学 | 针对均匀几何变分节块法的三维准输运加速方法 |
CN114444326A (zh) * | 2022-02-25 | 2022-05-06 | 上海交通大学 | 针对三维瞬态多群中子扩散方程的非结构网格变分节块方法 |
CN116933553A (zh) * | 2023-08-02 | 2023-10-24 | 上海交通大学 | 数值反应堆中子学的非结构网格体积修正方法 |
CN116933553B (zh) * | 2023-08-02 | 2024-02-13 | 上海交通大学 | 数值反应堆中子学的非结构网格体积修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107145657B (zh) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107145657A (zh) | 针对反应堆中子扩散方程的非均匀几何变分节块方法 | |
CN103150424B (zh) | 一种获取反应堆堆芯三维中子通量密度精细分布的方法 | |
Yang et al. | A parameter-free generalized moment limiter for high-order methods on unstructured grids | |
Merzari et al. | Cardinal: A lower-length-scale multiphysics simulator for pebble-bed reactors | |
CN105426657B (zh) | 一种消除反应堆堆芯计算中控制棒尖齿效应的方法 | |
Zhao et al. | Improved leakage splitting method for the 2D/1D transport calculation | |
CN103294899A (zh) | 一种计算小型实验反应堆堆芯中子通量分布的方法 | |
CN104198969B (zh) | 一种梯度线圈设计方法 | |
CN106202867B (zh) | 一种计算快中子堆组件轴向肿胀效应的方法 | |
Li et al. | Coupling of RMC and CFX for analysis of pebble bed-advanced high temperature reactor core | |
Graham et al. | Subplane collision probabilities method applied to control rod cusping in 2D/1D | |
CN107066751A (zh) | 针对非均匀几何变分节块方法的平源加速方法 | |
Jamalipour et al. | PWR power distribution flattening using quantum particle swarm intelligence | |
CN107145472A (zh) | 针对非均匀几何变分节块方法的广义矩阵分离加速方法 | |
CN114444326A (zh) | 针对三维瞬态多群中子扩散方程的非结构网格变分节块方法 | |
CN107563038A (zh) | 一种新型的接触热阻有限元求解方法 | |
Sun et al. | A variational nodal formulation for multi-dimensional unstructured neutron diffusion problems | |
CN107103189A (zh) | 一种快速进行反应堆临界曲率搜索的方法 | |
CN107423542A (zh) | 一种适用于逐棒计算的非均匀泄漏修正方法 | |
Jain et al. | DIAMOND–A Method of Characteristics neutron transport code using unstructured meshing | |
Saxena | An adaptive material mask overlay method: Modifications and investigations on binary, well connected robust compliant continua | |
Gibson et al. | On the stability of the Discrete Generalized Multigroup method | |
Marino | Starting Point, Keys and Milestones of a Computer Code for the Simulation of the Behaviour of a Nuclear Fuel Rod | |
Lee et al. | Development of the CAPP code based on the finite element method for the analysis of VHTR cores | |
CN114547952A (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 |