CN111523234B - 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法 - Google Patents
一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法 Download PDFInfo
- Publication number
- CN111523234B CN111523234B CN202010328810.5A CN202010328810A CN111523234B CN 111523234 B CN111523234 B CN 111523234B CN 202010328810 A CN202010328810 A CN 202010328810A CN 111523234 B CN111523234 B CN 111523234B
- Authority
- CN
- China
- Prior art keywords
- neutron
- layer
- reactor core
- pressurized water
- axial
- 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
- 230000004907 flux Effects 0.000 title claims abstract description 111
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000008859 change Effects 0.000 claims abstract description 4
- 239000000463 material Substances 0.000 claims description 10
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 238000003780 insertion Methods 0.000 claims description 4
- 230000037431 insertion Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 6
- 230000010354 integration Effects 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000005293 physical law 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法,首先将所要模拟的三维压水堆堆芯沿轴向划分若干层,针对每一层将三维的中子角通量展开,轴向采用线性模型近似模拟中子通量轴向变化;根据轴向线性展开后的三维中子通量,径向基于特征线方法建立二维中子输运模型;根据边界条件获得0阶和1阶的中子角通量的关系,以此简化特征线方法建立的二维中子输运模型;最后,轴向从下到上沿每一层求解简化后的二维中子输运模型,获得每一层的三维压水堆堆芯的中子通量分布。本发明相比现有技术,避免了二维中子输运模型和一维中子输运模型中横向积分产生的中子泄漏项,求解稳定性好,可用于数值反应堆的输运模块计算,提高数值模拟的稳定性。
Description
技术领域
本发明涉及核反应堆堆芯设计和安全领域,具体涉及一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法。
背景技术
随着核电行业的不断发展,为了应对核电厂延寿及安全分析,对压水堆堆芯的数值模拟的精度要求和稳定性要求越来越高。为了满足压水堆数值模拟的精度和稳定性要求,数值反应堆这类高保真方法越来越多的应用在压水堆的数值模拟当中。
高保真方法采用全堆芯直接求解,计算的网格数目较多,其输运模块的计算三维中子输运模型的负担较大。目前常采用横向积分的方式,将三维中子输运模型转化为二维中子输运模型和一维中子输运模型分别计算,并通过泄漏项对两个模型进行耦合。
由于泄漏项的存在,使得二维一维数值模拟过程中的中子源为负,分别计算一维中子输运模型和二维中子输运模型的过程中中子通量会出现不符合物理规律的负值,导致数值模拟计算发散。
发明内容
为了克服上述现有技术中存在的问题,本发明提供一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法,该方法避免采用传统的横向积分的方式耦合计算一维中子输运模型和二维中子输运模型,与传统方法相比,避免数值模拟过程中因泄漏项导致的负中子源和负中子通量的问题,增加计算过程的稳定性。
为了实现上述目的,本发明采用了以下技术方案予以实施:
一种针对压水堆堆芯的三维中子通量数值模拟方法,包括如下步骤:
步骤1:读取所要模拟的压水堆堆芯的几何信息、材料信息和边界条件;
步骤2:根据步骤1获取的压水堆堆芯的几何信息和材料信息,在控制棒的插入深度位置以及在轴向材料中子吸收截面发生变化的位置,将压水堆堆芯几何沿轴向切分为若干层,针对每一层将三维中子角通量展开,轴向每一层采用线性模型近似模拟中子角通量轴向变化,如公式(1)所示;线性模型的基函数如公式(2)所示;压水堆堆芯轴向的每一层的上下表面的中子角通量如公式(3)所示;
其中,
p—压水堆堆芯轴向的层数编号;
g—中子输运模型的能群编号;
m—中子输运模型的角度方向编号;
n—中子角通量密度沿轴向展开的阶数编号;
x—中子角通量密度的径向的网格横坐标;
y—中子角通量密度的径向的网格纵坐标;
z—中子角通量密度的轴向的网格竖坐标;
b0(z)—第0阶的轴向展开的基函数;
b1(z)—第1阶的轴向展开的基函数;
bn(z)—第n阶的轴向展开的基函数;
z+—压水堆堆芯轴向的每一层的上表面的位置坐标值;
z-—压水堆堆芯轴向的每一层的下表面的位置坐标值;
Δz—压水堆堆芯轴向的每一层的高度值;
步骤3:三维中子角通量沿轴向线性展开后,对于每一层,将三维中子角通量输运模型以轴向基函数为权重在轴向的每一层内进行积分,并采用特征线方法建立两个二维的中子通量输运模型,如公式(4.1)和公式(4.2)所示;
其中,
ξm—沿m方向的中子角通量与竖轴的夹角余弦;
Δzp—压水堆堆芯轴向的每一层的高度;
s—压水堆堆芯每一层内径向的特征线的编号;
步骤4:根据轴向每一层在反应堆堆芯的位置,建立边界面上第0阶和第1阶的二维中子角通量的关系;对于真空边界,使真空边界表面的入射中子角通量为零;对于全反射边界,全反射边界表面的入射中子角通量由对应的反射角的出射中子角通量反射得到;对于连续边界,入射中子角通量与相邻层的这个表面的出射中子角通量相同;如公式(5.1)、公式(5.2)和公式(5.3)所示;
真空边界条件
全反射边界条件
连续边界条件
其中,
m′—m方向的在上下表面的反射角方向;
步骤5:根据步骤4中第0阶和第1阶的二维中子通量的关系,获得不同边界条件下的0阶二维中子输运模型,如公式(6.1)、公式(6.2)和公式(6.3)所示;
真空边界条件
全反射边界条件
连续边界条件
步骤6:轴向从下往上沿每一层求解基于特征线方法建立的0阶二维的中子输运模型,获得压水堆堆芯的三维中子通量分布。
与现有技术相比,本发明有如下突出优点:
本发明避免采用传统的横向积分的方式耦合计算一维中子输运模型和二维中子输运模型,与传统方法相比,沿轴向采用线性模型近似展开三维中子通量,基于线性展开后的三维中子通量,沿径向采用特征线方法建立多个二维的中子输运模型,避免数值模拟过程中因泄漏项导致的负中子源和负中子通量的问题,增加计算过程的稳定性。
附图说明
图1是本发明方法流程图。
图2是压水堆堆芯轴向线性模型近似的计算对象示意图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明。
具体步骤如图1所示,本发明将所要模拟的三维压水堆堆芯沿轴向划分为若干层,在每一层基于特征线方法建立两个二维中子输运模型;根据每一层的边界条件,简化二维中子输运模型;从下往上沿每一层求解二维中子输运模型,获得压水堆的三维中子通量。具体步骤如下:
步骤1:读取所要模拟的压水堆堆芯的几何信息、材料信息和边界条件;
步骤2:根据步骤1获取的压水堆堆芯的几何信息和材料信息,在控制棒的插入深度位置以及在轴向材料中子吸收截面发生变化的位置,将压水堆堆芯几何沿轴向切分为若干层,针对每一层将三维中子角通量展开,轴向每一层采用线性模型近似模拟中子角通量轴向变化,如公式(1)所示;线性模型的基函数如公式(2)所示;压水堆堆芯轴向的每一层的上下表面的中子角通量如公式(3)所示;
其中,
p—压水堆堆芯轴向的层数编号;
g—中子输运模型的能群编号;
m—中子输运模型的角度方向编号;
n—中子角通量密度沿轴向展开的阶数编号;
x—中子角通量密度的径向的网格横坐标;
y—中子角通量密度的径向的网格纵坐标;
z—中子角通量密度的轴向的网格竖坐标;
b0(z)—第0阶的轴向展开的基函数;
b1(z)—第1阶的轴向展开的基函数;
bn(z)—第n阶的轴向展开的基函数;
z+—压水堆堆芯轴向的每一层的上表面的位置坐标值;
z-—压水堆堆芯轴向的每一层的下表面的位置坐标值;
z—压水堆堆芯轴向的每一层的中间位置坐标值;
Δz—压水堆堆芯轴向的每一层的高度值;
步骤3:三维中子角通量沿轴向线性展开后,对于每一层,将三维中子角通量输运模型以轴向基函数为权重在轴向的每一层内进行积分,并采用特征线方法建立两个二维的中子通量输运模型,如公式(4.1)和公式(4.2)所示;
其中,
ξm—沿m方向的中子角通量与竖轴的夹角余弦;
Δzp—压水堆堆芯轴向的每一层的高度;
s—压水堆堆芯每一层内径向的特征线的编号;
步骤4:根据轴向每一层在反应堆堆芯的位置,建立边界面上第0阶和第1阶的二维中子角通量的关系;对于真空边界,使真空边界表面的入射中子角通量为零;对于全反射边界,全反射边界表面的入射中子角通量由对应的反射角的出射中子角通量反射得到;对于连续边界,入射中子角通量与相邻层的这个表面的出射中子角通量相同;如公式(5.1)、公式(5.2)和公式(5.3)所示;
真空边界条件
全反射边界条件
连续边界条件
其中,
m′—m方向的在上下表面的反射角方向;
步骤5:根据步骤4中第0阶和第1阶的二维中子通量的关系,获得不同边界条件下的0阶二维中子输运模型,如公式(6.1)、公式(6.2)和公式(6.3)所示;
真空边界条件
全反射边界条件
连续边界条件
步骤6:轴向从下往上沿每一层求解基于特征线方法建立的0阶二维的中子输运模型,获得压水堆堆芯的三维中子通量分布。
图2展示的是压水堆堆芯。首先在控制棒的插入深度位置以及在轴向材料中子吸收截面发生变化的位置,将压水堆堆芯几何沿轴向切分若干层;对于每一层的轴向采用线性模型近似展开三维中子通量,并沿径向基于特征线方法建立两个二维中子输运模型;其中,Ωm表示特征线的方向。
Claims (1)
1.一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法,其特征在于:包括如下步骤:
步骤1:读取所要模拟的压水堆堆芯的几何信息、材料信息和边界条件;
步骤2:根据步骤1获取的压水堆堆芯的几何信息和材料信息,在控制棒的插入深度位置以及在轴向材料中子吸收截面发生变化的位置,将压水堆堆芯几何沿轴向切分为若干层,针对每一层将三维中子角通量展开,轴向每一层采用线性模型近似模拟中子角通量轴向变化,如公式(1)所示;线性模型的基函数如公式(2)所示;压水堆堆芯轴向的每一层的上下表面的中子角通量如公式(3)所示;
其中,
p—压水堆堆芯轴向的层数编号;
g—中子输运模型的能群编号;
m—中子输运模型的角度方向编号;
n—中子角通量密度沿轴向展开的阶数编号;
x—中子角通量密度的径向的网格横坐标;
y—中子角通量密度的径向的网格纵坐标;
z—中子角通量密度的轴向的网格竖坐标;
b0(z)—第0阶的轴向展开的基函数;
b1(z)—第1阶的轴向展开的基函数;
bn(z)—第n阶的轴向展开的基函数;
z+—压水堆堆芯轴向的每一层的上表面的位置坐标值;
z-—压水堆堆芯轴向的每一层的下表面的位置坐标值;
Δz—压水堆堆芯轴向的每一层的高度值;
步骤3:三维中子角通量沿轴向线性展开后,对于每一层,将三维中子角通量输运模型以轴向基函数为权重在轴向的每一层内进行积分,并采用特征线方法建立两个二维的中子通量输运模型,如公式(4.1)和公式(4.2)所示;
ξm>0
ξm<0
其中,
ξm—沿m方向的中子角通量与竖轴的夹角余弦;
Δzp—压水堆堆芯轴向的每一层的高度;
s—压水堆堆芯每一层内径向的特征线的编号;
步骤4:根据轴向每一层在反应堆堆芯的位置,建立边界面上第0阶和第1阶的二维中子角通量的关系;对于真空边界,使真空边界表面的入射中子角通量为零;对于全反射边界,全反射边界表面的入射中子角通量由对应的反射角的出射中子角通量反射得到;对于连续边界,入射中子角通量与相邻层的这个表面的出射中子角通量相同;如公式(5.1)、公式(5.2)和公式(5.3)所示;
真空边界条件
全反射边界条件
连续边界条件
其中,
m′—m方向的在上下表面的反射角方向;
步骤5:根据步骤4中第0阶和第1阶的二维中子通量的关系,获得不同边界条件下的0阶二维中子输运模型,如公式(6.1)、公式(6.2)和公式(6.3)所示;
真空边界条件
全反射边界条件
连续边界条件
步骤6:轴向从下往上沿每一层求解基于特征线方法建立的0阶二维的中子输运模型,获得压水堆堆芯的三维中子通量分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010328810.5A CN111523234B (zh) | 2020-04-23 | 2020-04-23 | 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010328810.5A CN111523234B (zh) | 2020-04-23 | 2020-04-23 | 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111523234A CN111523234A (zh) | 2020-08-11 |
CN111523234B true CN111523234B (zh) | 2021-10-19 |
Family
ID=71903609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010328810.5A Active CN111523234B (zh) | 2020-04-23 | 2020-04-23 | 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111523234B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113672849A (zh) * | 2021-08-26 | 2021-11-19 | 中国核动力研究设计院 | 一种基于轴向通量展开的一步法输运计算方法及系统 |
CN114528719B (zh) * | 2022-04-22 | 2022-07-08 | 西安交通大学 | 一种基于二维堆芯的用于压水堆的在线能群压缩方法 |
CN117593474B (zh) * | 2024-01-18 | 2024-04-09 | 西安交通大学 | 一种压水堆堆芯三维功率分布重构方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150424A (zh) * | 2013-02-05 | 2013-06-12 | 西安交通大学 | 一种获取反应堆堆芯三维中子通量密度精细分布的方法 |
CN106126925A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种改进反应堆堆芯三维中子通量密度精细分布的方法 |
CN106126932A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种压水堆节块法控制棒尖齿效应的处理方法 |
CN106773666A (zh) * | 2016-11-11 | 2017-05-31 | 中国电力科学研究院 | 一种针对压水堆一回路系统的模型参数获取方法 |
CN108647402A (zh) * | 2018-04-19 | 2018-10-12 | 西安交通大学 | 一种获取六角形几何压水堆堆芯反射层参数的方法 |
CN110717275A (zh) * | 2019-10-12 | 2020-01-21 | 西安交通大学 | 一种针对压水堆堆芯的三维中子通量数值模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101614772B1 (ko) * | 2014-12-19 | 2016-04-22 | 한전원자력연료 주식회사 | 신경망 회로를 이용한 원자로 노심의 축방향 출력분포 합성방법 및 그 방법이 적용된 노심보호계통 |
-
2020
- 2020-04-23 CN CN202010328810.5A patent/CN111523234B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150424A (zh) * | 2013-02-05 | 2013-06-12 | 西安交通大学 | 一种获取反应堆堆芯三维中子通量密度精细分布的方法 |
CN106126925A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种改进反应堆堆芯三维中子通量密度精细分布的方法 |
CN106126932A (zh) * | 2016-06-24 | 2016-11-16 | 西安交通大学 | 一种压水堆节块法控制棒尖齿效应的处理方法 |
CN106773666A (zh) * | 2016-11-11 | 2017-05-31 | 中国电力科学研究院 | 一种针对压水堆一回路系统的模型参数获取方法 |
CN108647402A (zh) * | 2018-04-19 | 2018-10-12 | 西安交通大学 | 一种获取六角形几何压水堆堆芯反射层参数的方法 |
CN110717275A (zh) * | 2019-10-12 | 2020-01-21 | 西安交通大学 | 一种针对压水堆堆芯的三维中子通量数值模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111523234A (zh) | 2020-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111523234B (zh) | 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法 | |
CN111523233B (zh) | 一种针对三维压水堆堆芯的中子输运计算方法 | |
CN109830317B (zh) | 一种基于跟踪计算的堆芯功率本征正交分解在线重构方法 | |
Nagy | The prism method for terrain corrections using digital computers | |
KR20060015760A (ko) | 핵 반응로의 코어를 위한 컴퓨터 모델링 방법 | |
CN103294899A (zh) | 一种计算小型实验反应堆堆芯中子通量分布的方法 | |
CN106126925A (zh) | 一种改进反应堆堆芯三维中子通量密度精细分布的方法 | |
CN111209690B (zh) | Fcm燃料元件中triso燃料颗粒随机分布的建模方法 | |
CN107038293B (zh) | 针对板状燃料的基于等效一维板模型的共振自屏计算方法 | |
CN110717275B (zh) | 一种针对压水堆堆芯的三维中子通量数值模拟方法 | |
CN108090292A (zh) | 一种宽弦风扇叶片二维有限元建模方法 | |
CN113139306A (zh) | 一种复合材料螺旋桨空化噪声的数值预报方法 | |
CN107886573B (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
CN105512387A (zh) | 一种燃料组件精细功率重构方法 | |
CN114528719B (zh) | 一种基于二维堆芯的用于压水堆的在线能群压缩方法 | |
CN113672849A (zh) | 一种基于轴向通量展开的一步法输运计算方法及系统 | |
CN114757123B (zh) | 一种用于板形核燃料堆芯的跨维度流固耦合分析方法 | |
CN108038262A (zh) | 一种考虑sssi效应的楼层反应谱简化计算方法 | |
CN114547952B (zh) | 一种基于深度学习-界面流耦合的压水堆中子输运方法 | |
CN114282403B (zh) | 一种耦合生境适宜模型的高效高精度栖息地模拟方法 | |
CN110705184A (zh) | 一种反应堆堆芯精细化数值求解的虚拟体积力动量源法 | |
Baiwei et al. | Ship resistance performance optimization design based on CAD/CFD | |
CN115563894A (zh) | 一种复杂分层及拼接数字岩心的渗透率预测方法 | |
CN113139307A (zh) | 一种复合材料螺旋桨的振动噪声数值预报方法 | |
Frepoli | Numerical simulation of the inlet plenum flow field for a LWR reactor pressure vessel |
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 |