CN108647402A - 一种获取六角形几何压水堆堆芯反射层参数的方法 - Google Patents

一种获取六角形几何压水堆堆芯反射层参数的方法 Download PDF

Info

Publication number
CN108647402A
CN108647402A CN201810353519.6A CN201810353519A CN108647402A CN 108647402 A CN108647402 A CN 108647402A CN 201810353519 A CN201810353519 A CN 201810353519A CN 108647402 A CN108647402 A CN 108647402A
Authority
CN
China
Prior art keywords
locking nub
reactor core
fuel assembly
close
fuel
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
Application number
CN201810353519.6A
Other languages
English (en)
Other versions
CN108647402B (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.)
Xian Jiaotong University
Original Assignee
Xian 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 Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201810353519.6A priority Critical patent/CN108647402B/zh
Publication of CN108647402A publication Critical patent/CN108647402A/zh
Application granted granted Critical
Publication of CN108647402B publication Critical patent/CN108647402B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

一种获取六角形几何压水堆堆芯反射层参数的方法，根据反射层组件相对燃料组件的位置，将反射层节块分为中心反射层节块与非中心反射层节块；中心反射层节块再根据与燃料组件交界面的个数分为单面紧贴燃料、双面紧贴燃料与三面紧贴燃料的三类中心反射层节块；以各个中心反射层节块为中心，参考实际堆芯布置进行六角形阵列排布，建立反射层组件二维模型；进行非均匀输运计算，获得各个反射层节块的中子通量密度分布；进行均匀化计算，获得均匀化少群截面及不连续因子，本发明针对六角形组件式压水堆的几何及材料特性，对堆芯不同位置处的反射层组件进行建模并计算出相应的少群均匀化常数，为下游堆芯计算提供了条件。

Description

一种获取六角形几何压水堆堆芯反射层参数的方法
技术领域
本发明涉及核反应堆堆芯核设计及安全领域，具体涉及一种获取六角形几何压水堆堆芯反射层参数的方法。
背景技术
随着核电行业的不断发展，各种新式堆型蓬勃发展，但商业应用最广泛的还是压水堆。相比于快中子堆，压水堆作为一种热谱堆芯，具有中子自由程较短，中子能谱较软，各向异性较不明显等特点。因此，传统的压水堆燃料管理计算程序多采用“两步法”对压水堆运行的过程进行模拟仿真，继而保证堆芯运行的安全。“两步法”，是指先对采用了反射边界条件的单个或多个组件进行非均匀输运计算，获得该问题的精细中子能谱，再藉由该能谱进行空间归并能群压缩获得组件大小节块的等效均匀化常数；通过第一步获得的不同节块不同状态点建立的少群参数库，对整个堆芯进行低价扩散计算。
“两步法”中，较为关键的一个环节，即是反射层参数的获取。在较为合适的计算条件下获得较为精确的反射层参数，对于堆芯组件功率分布，有效增殖因子等各类参数的计算有着极为重要的影响。然而，随着针对不同用途的核设计的不断提出，压水堆组件不再局限于传统的矩形设计，六角形组件式堆芯也有着愈来愈广泛的应用。
传统压水堆“两步法”计算中，多采用一维反射层模型。在传统压水堆堆芯内，大多数反射层组件只有一个面与燃料组件交界，一维模型已能保证足够精度；然而，在六角形组件式堆芯中，只有极少数反射层组件仅有一面挨着燃料组件，很多有两个面与燃料组件交界，甚至六个边界面中的三个都紧挨着不同的燃料组件。因此，沿用传统一维反射层模型对六角形组件式压水堆堆芯来说，是不适用的。一种考虑了堆芯几何材料布置的二维反射层组件计算方法能够有效的完成对六角形几何压水堆的“两步法”模拟与计算。
发明内容
为了克服上述现有技术存在的问题，本发明的目的在于提供一种获取六角形几何压水堆堆芯反射层参数的方法，针对六角形组件式压水堆的几何及材料特性，对堆芯不同位置处的反射层组件进行建模并计算出相应的少群均匀化常数，为下游堆芯计算提供了条件。
为了实现上述目的，本发明采用了以下技术方案予以实施：
一种获取六角形几何压水堆堆芯反射层参数的方法，包括如下步骤：
步骤1：读取压水堆堆芯几何尺寸，根据围板、吊篮结构材料的大小确认反射层模型建模的区域，并根据燃料组件对边距划分六角形节块；
步骤2：根据步骤1中划分的六角形节块，对处于堆芯不同位置的反射层组件按相对燃料组件的距离划分为紧贴燃料组件与非紧贴燃料组件的两类反射层节块，其中，由于非紧贴燃料组件的反射层节块相对于紧贴燃料组件的反射层节块的通量水平要低，建立二维模型时只要将其包含在另一个紧贴燃料组件的反射层节块的计算模型中即可；对于紧贴燃料组件的反射层节块(以下简称为中心反射层节块，非紧贴燃料组件的反射层节块称为非中心反射层节块)则根据紧贴排布的燃料组件个数分为单面紧贴燃料的中心反射层节块、两面紧贴燃料的中心反射层节块和三面紧贴燃料的中心反射层节块；
步骤3：对于步骤2划分出的反射层节块，根据实际堆芯中燃料组件的对边距的大小，确认反射层模型中燃料组件的个数与排数；沿堆芯径向辐射方向要包含一至四层燃料组件(实际堆芯内装填的燃料组件对边距越小，反射层模型中应包含的燃料组件越多)；至此，应包含在一个六角形几何压水堆反射层模型内的六角形节块包含以下三类：中心反射层节块一个；非中心反射层节块若干；燃料组件节块若干；
步骤4：将步骤3中的三类六角形节块以中心反射层六角形节块为中心，参照实际堆芯布置进行六角形阵列排布，堆芯径向辐射方向必须包含步骤3中确认的燃料组件排数；堆芯周向则根据六角形组件压水堆燃料管理计算程序的几何处理能力，权衡计算效率与精度(模型内的节块数越多，越接近实际堆芯布置，越精确)进行反射层模型建模；
步骤5：根据压水堆堆芯的几何结构与材料布置，精确描述各节块内的非均匀几何；至此，六角形几何压水堆的反射层组件二维模型建立完毕；基于此模型进行非均匀输运计算，获得所需的各个反射层节块的中子通量密度随空间、能量的分布；
步骤6：由步骤5中获得的中子通量密度随空间、能量的分布，根据等效均匀化常数计算的三大守恒条件(各能群的反应率守恒、节块的各个界面上的界面净流守恒和反应堆的特征值守恒)获得均匀化少群截面，并根据不连续因子的定义求得相应反射层组件的不连续因子。
与现有技术相比，本发明有如下突出优点：
考虑处于不同堆芯位置的反射层组件能谱受燃料组件的影响，并根据压水堆堆芯特性考虑了模型中燃料组件的数目及厚度，使最终获得的反射层组件能谱相较于简单的一维模型更为接近真实堆芯的实际能谱。
附图说明
图1是本发明方法流程图。
图2是十二分之一堆芯示意图。
图3是六角形几何压水堆堆芯节块划分示意图。
图4是二维反射层模型示意图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细说明。
具体步骤如图1所示。本发明是一种获取六角形几何压水堆堆芯反射层参数的方法，主要包括六角形节块的划分与分类、燃料组件节块数目的选取以及二维反射层组件模型的建立，具体步骤如下：
1：读取压水堆堆芯几何尺寸，根据围板、吊篮等结构材料的大小确认反射层建模的区域(以田湾核电引进的俄罗斯VVER-1000为例，堆芯布置见图2)，并根据燃料组件对边距划分六角形节块；
2：根据1中划分的六角形节块，对处于堆芯不同位置的反射层组件按相对燃料组件的距离划分为紧贴燃料组件与非紧贴燃料组件的两类反射层节块，见图3。其中，由于非紧贴燃料组件的反射层节块相对于紧贴燃料组件的反射层节块的通量水平要低，建立二维模型时只要将其包含在另一个紧贴燃料组件的反射层节块的计算模型中即可；对于紧贴燃料组件的反射层节块(以下简称为中心节块，非紧贴燃料组件的反射层节块称为非中心节块)则根据紧贴排布的燃料组件个数分为单面紧贴燃料的中心节块、两面紧贴燃料的中心节块和三面紧贴燃料的中心节块。对于不同类别的中心节块，需要建立不同的几何模型(紧贴中心节块的燃料组件节块必须包含在建立的二维模型中)；
3：对于2划分出的反射层节块，根据实际堆芯中燃料组件的对边距的大小，确认反射层模型中燃料组件的个数与排数。沿堆芯径向辐射方向的燃料组件节块的排数不少于下式中确定的N：
式中，P为燃料组件对边距，cm；N最终取值范围在1～4之间。
至此，应包含在一个六角形几何压水堆反射层模型内的六角形节块应包含以下三类：中心节块一个；非中心节块若干；燃料组件节块若干；
4：将3中的三类六角形节块以中心反射层六角形节块为中心，参照实际堆芯布置进行六角形阵列排布，堆芯径向辐射方向必须包含3中确认的燃料组件排数；堆芯周向则根据六角形组件压水堆燃料管理计算程序的几何处理能力，权衡计算效率与精度(模型内的节块数越多，越接近实际堆芯布置，越精确)进行建模。以图3中编号为1的单面中心节块和编号为2的非中心节块为例，最终的二维模型见图4；
5：根据压水堆堆芯的几何结构与材料布置，精确描述4中所建模型中所包含的各节块内的非均匀几何。至此，六角形几何压水堆的反射层组件二维模型建立完毕。基于此模型进行非均匀输运计算，获得所需的各个反射层节块的中子通量密度随空间、能量的分布；
6：由5中获得的中子通量密度随空间、能量的分布，根据等效均匀化常数计算的三大守恒条件(各能群的反应率守恒、节块的各个界面上的界面净流守恒和反应堆的特征值守恒)中的各能群反应率守恒原则，并以组件内非均匀中子标通量密度的输运解近似替代反应堆堆芯中的中子标通量密度：
获得均匀化少群截面；
并根据不连续因子的定义：
求得该节块各个边界面的不连续因子。
在以本发明建立的二维反射层组件模型中，燃料组件的数目的确认参考了传统M310堆芯组件边长与实际要计算六角形几何堆芯反射层组件对边距的比值，确保了燃料组件侧外部边界条件不会太大地影响模型内部的中子能谱，保证了获取的反射层参数具有足够的精度。
本发明在获取六角形几何堆芯反射层参数方面，具有可靠的精度及较高的效率，适应工程实际中的计算要求。

Claims (1)

1.一种获取六角形几何压水堆堆芯反射层参数的方法，其特征在于：包括如下步骤：
步骤1：读取压水堆堆芯几何尺寸，根据围板、吊篮结构材料的大小确认反射层模型建模的区域，并根据燃料组件对边距划分六角形节块；
步骤2：根据步骤1中划分的六角形节块，对处于堆芯不同位置的反射层组件按相对燃料组件的距离划分为紧贴燃料组件与非紧贴燃料组件的两类反射层节块，其中，由于非紧贴燃料组件的反射层节块相对于紧贴燃料组件的反射层节块的通量水平要低，建立二维模型时只要将其包含在另一个紧贴燃料组件的反射层节块的计算模型中即可；以下简称为中心反射层节块，非紧贴燃料组件的反射层节块称为非中心反射层节块，对于紧贴燃料组件的反射层节块，则根据紧贴排布的燃料组件个数分为单面紧贴燃料的中心反射层节块、两面紧贴燃料的中心反射层节块和三面紧贴燃料的中心反射层节块；
步骤3：对于步骤2划分出的反射层节块，根据实际堆芯中燃料组件的对边距的大小，确认反射层模型中燃料组件的个数与排数，沿堆芯径向辐射方向的燃料组件节块的排数不少于下式中确定的N：
N最终取值范围在1～4之间；至此，应包含在一个六角形几何压水堆反射层模型内的六角形节块包含以下三类：中心反射层节块一个；非中心反射层节块若干；燃料组件节块若干；
步骤4：将步骤3中的三类六角形节块以中心反射层六角形节块为中心，参照实际堆芯布置进行六角形阵列排布，堆芯径向辐射方向必须包含步骤3中确认的燃料组件排数；堆芯周向则根据六角形组件压水堆燃料管理计算程序的几何处理能力，权衡计算效率与精度进行反射层模型建模；
步骤5：根据压水堆堆芯的几何结构与材料布置，精确描述各节块内的非均匀几何；至此，六角形几何压水堆的反射层组件二维模型建立完毕；基于此模型进行非均匀输运计算，获得所需的各个反射层节块的中子通量密度随空间、能量的分布；
步骤6：由步骤5中获得的中子通量密度随空间、能量的分布，根据等效均匀化常数计算的三大守恒条件即各能群的反应率守恒、节块的各个界面上的界面净流守恒和反应堆的特征值守恒获得均匀化少群截面，并根据不连续因子的定义求得相应反射层组件的不连续因子。
CN201810353519.6A 2018-04-19 2018-04-19 一种获取六角形几何压水堆堆芯反射层参数的方法 Active CN108647402B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810353519.6A CN108647402B (zh) 2018-04-19 2018-04-19 一种获取六角形几何压水堆堆芯反射层参数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810353519.6A CN108647402B (zh) 2018-04-19 2018-04-19 一种获取六角形几何压水堆堆芯反射层参数的方法

Publications (2)

Publication Number Publication Date
CN108647402A true CN108647402A (zh) 2018-10-12
CN108647402B CN108647402B (zh) 2019-04-09

Family

ID=63746835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810353519.6A Active CN108647402B (zh) 2018-04-19 2018-04-19 一种获取六角形几何压水堆堆芯反射层参数的方法

Country Status (1)

Country Link
CN (1) CN108647402B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110598324A (zh) * 2019-09-12 2019-12-20 西安交通大学 一种核反应堆弥散型板型燃料元件堆芯流固耦合计算方法
CN111523234A (zh) * 2020-04-23 2020-08-11 西安交通大学 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法
CN114707189A (zh) * 2022-06-02 2022-07-05 西安交通大学 一种等效模拟压水堆堆芯内燃料组件弯曲的方法
CN115295187A (zh) * 2022-07-08 2022-11-04 中国核动力研究设计院 一种六角形组件堆芯物理计算坐标转换方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4038027A1 (de) * 1990-11-29 1992-06-11 Forschungszentrum Juelich Gmbh Deckenreflektor eines mittelgrossen hochtemperaturreaktors
CN101719384A (zh) * 2009-12-03 2010-06-02 中国核动力研究设计院 具有快中子转换区的核反应堆堆芯
CN103294899A (zh) * 2013-05-10 2013-09-11 西安交通大学 一种计算小型实验反应堆堆芯中子通量分布的方法
CN103871529A (zh) * 2014-03-26 2014-06-18 清华大学 一种球床式高温气冷堆的底反射层结构
CN105590658A (zh) * 2015-12-29 2016-05-18 中国科学院合肥物质科学研究院 一种利用中子管驱动的次临界微型反应堆
CN106981319A (zh) * 2017-01-10 2017-07-25 上海核工程研究设计院 一种采用无铍熔盐氧化铍慢化的熔盐堆堆芯
CN107066745A (zh) * 2017-04-20 2017-08-18 西安交通大学 获取快中子堆堆芯瞬态过程三维中子通量密度分布的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4038027A1 (de) * 1990-11-29 1992-06-11 Forschungszentrum Juelich Gmbh Deckenreflektor eines mittelgrossen hochtemperaturreaktors
CN101719384A (zh) * 2009-12-03 2010-06-02 中国核动力研究设计院 具有快中子转换区的核反应堆堆芯
CN103294899A (zh) * 2013-05-10 2013-09-11 西安交通大学 一种计算小型实验反应堆堆芯中子通量分布的方法
CN103871529A (zh) * 2014-03-26 2014-06-18 清华大学 一种球床式高温气冷堆的底反射层结构
CN105590658A (zh) * 2015-12-29 2016-05-18 中国科学院合肥物质科学研究院 一种利用中子管驱动的次临界微型反应堆
CN106981319A (zh) * 2017-01-10 2017-07-25 上海核工程研究设计院 一种采用无铍熔盐氧化铍慢化的熔盐堆堆芯
CN107066745A (zh) * 2017-04-20 2017-08-18 西安交通大学 获取快中子堆堆芯瞬态过程三维中子通量密度分布的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HESSAM MALMIR等: ""Comparison between triangular and hexagonal modeling of a hexagonal-structured reactor core using box method"", 《ANNALS OF NUCLEAR ENERGY》 *
SIEGFRIED MITTAG等: ""DISCONTINUITY FACTORS FOR NON-MULTIPLYING MATERIAL IN TWO-DIMENSIONAL HEXAGONAL REACTOR GEOMETRY"", 《ANNALS OF NUCLEAR ENERGY》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110598324A (zh) * 2019-09-12 2019-12-20 西安交通大学 一种核反应堆弥散型板型燃料元件堆芯流固耦合计算方法
CN110598324B (zh) * 2019-09-12 2020-08-25 西安交通大学 一种核反应堆弥散型板型燃料元件堆芯流固耦合计算方法
CN111523234A (zh) * 2020-04-23 2020-08-11 西安交通大学 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法
CN111523234B (zh) * 2020-04-23 2021-10-19 西安交通大学 一种基于轴向展开的模拟压水堆堆芯三维中子通量的方法
CN114707189A (zh) * 2022-06-02 2022-07-05 西安交通大学 一种等效模拟压水堆堆芯内燃料组件弯曲的方法
CN114707189B (zh) * 2022-06-02 2022-08-19 西安交通大学 一种等效模拟压水堆堆芯内燃料组件弯曲的方法
CN115295187A (zh) * 2022-07-08 2022-11-04 中国核动力研究设计院 一种六角形组件堆芯物理计算坐标转换方法
CN115295187B (zh) * 2022-07-08 2024-03-12 中国核动力研究设计院 一种六角形组件堆芯物理计算坐标转换方法

Also Published As

Publication number Publication date
CN108647402B (zh) 2019-04-09

Similar Documents

Publication Publication Date Title
CN108647402B (zh) 一种获取六角形几何压水堆堆芯反射层参数的方法
Ryu et al. Solution of the BEAVRS benchmark using the nTRACER direct whole core calculation code
CN107689256B (zh) 一种核反应堆次临界度测量方法
CN109830317B (zh) 一种基于跟踪计算的堆芯功率本征正交分解在线重构方法
CN109859867B (zh) 一种基于本征正交分解的堆芯三维中子通量实时重构方法
CN105404723A (zh) 一种精确计算燃料组件棒功率分布的方法
CN107038294B (zh) 针对轻水堆的基于等效一维棒模型的共振自屏计算方法
CN106126925A (zh) 一种改进反应堆堆芯三维中子通量密度精细分布的方法
Theler Unstructured grids and the multigroup neutron diffusion equation
CN105512387B (zh) 一种燃料组件精细功率重构方法
CN108447574A (zh) 一种核电站堆外探测器中子计数率计算方法及其系统
Yu et al. The development and validation of nuclear cross section processing code for reactor-RXSP
CN107122564A (zh) 一种压水堆pin‑by‑pin计算少群常数参数化方法
Ray et al. A collision probability and MOC-based lattice and burnup code for analysis of LWR fuel assemblies
Bahadir Improved PWR radial reflector modeling with simulate5
Palau et al. recent progess in the vetv of the french apollo3 code 3d full core analysis of the uh1. 2 integral experiment using idtmethod of characteristics
Aleshin et al. Calculations of 3D full-scale VVER fuel assembly and core models using MCU and BIPR-7A codes
Gorodkov Using core symmetry in Monte Carlo dominance ratio calculations
Bays et al. Upgrading limiting peak-power analysis techniques with modern validation and uncertainty quantification for the advanced test reactor
Singh et al. Development of a Lattice Neutron Transport Code Ν-Trac for 1-D Slab, 2-D Square and 2-D Hexagonal Geometries
Hunt A review of criticality safety models used in evaluating arrays of fissile materials
Merzari et al. Cardinal: A Lower Length-Scale Simulator for Fluoride Cooled High Temperature Reactors
Rätz et al. Experimental validation of control rod related perturbations of moderator regions in an SCWR-like fuel lattice
Ferrer et al. Particle Conservation and the Linear Source Method of Characteristics for Isotropic and Anisotropic Sources
Choi et al. Preliminary Simulation Results of Peach Bottom Unit 2 Cycles 1 and 2 with MPACT

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