CN114840988B - 压水堆堆芯自动建模方法 - Google Patents
压水堆堆芯自动建模方法 Download PDFInfo
- Publication number
- CN114840988B CN114840988B CN202210427670.6A CN202210427670A CN114840988B CN 114840988 B CN114840988 B CN 114840988B CN 202210427670 A CN202210427670 A CN 202210427670A CN 114840988 B CN114840988 B CN 114840988B
- Authority
- CN
- China
- Prior art keywords
- calculation
- reactor core
- fuel
- component
- modeling
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 17
- 239000000446 fuel Substances 0.000 claims abstract description 51
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 230000001052 transient effect Effects 0.000 claims description 14
- 239000002574 poison Substances 0.000 claims description 9
- 231100000614 poison Toxicity 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 8
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 claims description 7
- 229910052796 boron Inorganic materials 0.000 claims description 7
- 239000002826 coolant Substances 0.000 claims description 6
- 239000013589 supplement Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 5
- 230000001502 supplementing effect Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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/12—Geometric CAD characterised by design entry means specially adapted for CAD, e.g. graphical user interfaces [GUI] specially adapted for CAD
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
- G21C17/001—Mechanical simulators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/20—Configuration CAD, e.g. designing by assembling or positioning modules selected from libraries of predesigned modules
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Human Computer Interaction (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Architecture (AREA)
- Plasma & Fusion (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种压水堆堆芯自动建模方法,参照核电厂购买燃料组件时使用的类型。定义组件“型号”,建立或补充完善组件型号库;将按需“制造”的组件摆放或配插到堆芯中之后,只需要提出堆芯运行要求,即可由建模程序自动搜索并生成内核计算软件所需的输入文件、自动并发调用内核计算程序、自动完成计算结果的收集整理,实现了前台建模输入与后台计算方法实现完全脱耦,并与工程实际保持高度一致,也实现了全自动高可靠性处理和全范围高精细建模。本发明能简化建模计算流程,满足现有压水堆两步法计算程序建模计算的需求。解决了压水堆堆芯建模过程中的人员要求高、操作复杂易出错、工作效率低等问题。
Description
技术领域
本发明涉及核反应堆物理数值仿真计算领域,具体涉及一种压水堆堆芯自动建模方法。
背景技术
压水堆堆芯分析软件在商用压水堆核电站的堆芯设计、运行分析、安全审评等多个过程中起到了关键的作用。考虑到计算资源和计算代价,现有的压水堆堆芯分析软件多采用两步法技术,具体过程包括:先利用栅格计算软件,针对整个堆芯内所有不同的燃料组件横剖面(称之为燃料栅格,一般有10-100种不等)进行二维燃料栅格的建模计算,再对紧临活性区大约一个燃料组件宽度内的非活性区(称之为反射层,包括径向反射层和轴向反射层)进行反射层栅格建模计算(一般有1-40种不等),获得每一种栅格的均匀化少群常数(包括宏观截面、重要核素的微观截面、扩散系数、不连续因子、精细功率分布形状因子等)与栅格状态参数(包括燃耗深度、硼浓度、燃料有效温度、慢化剂温度、控制棒插入/抽出的状态等)之间的离散对应关系;再利用参数化技术,按照栅格类型逐一处理,获得均匀化少群常数与栅格状态参数之间的连续函数关系,贮存起来制成少群常数库;然后利用堆芯稳态过程计算软件,在三维全堆芯范围内进行中子输运-热量传输-水力流动-核素燃耗-运行控制等多物理多尺度耦合模拟与分析,获得稳态过程的仿真结果,并为极短时间的瞬态过程仿真提供数据基础;利用堆芯瞬态过程计算软件,在区分瞬发中子与缓发中子的条件下,完成瞬态过程的多物理多尺度耦合模拟与分析,获得瞬态过程的仿真结果。
上述传统建模方式要求用户不仅要输入堆芯实体参数,还需要了解两步法,才能安排好两个程序的建模计算任务,整个过程复杂繁琐、手动操作特别多、效率低下、容易出错却极难查错。同时,在建模过程中,尚有大量建模分析经验,仍主要通过手把手教学的方式在传承,建模计算的学习成本很高。
因此,针对以上提到的问题,需要发明一种效率更高,能够充分利用已有经验、不再要求软件用户必须掌握方法理论的自动建模方法。
发明内容
为了克服上述现有技术存在的问题,本发明的目的在于提供一种压水堆堆芯自动建模方法,提高压水堆两步法计算程序建模的效率,实现专家建模经验的信息化,降低建模过程的出错率,保证建模计算的可靠性和精度。
为了实现上述目的,本发明采取了以下技术方案:
一种压水堆堆芯自动建模方法,包括如下步骤:
步骤1:用户建立或补充完善组件型号库:参照核电厂购买燃料组件时使用的类型,定义组件“型号”,并按照实际需求,把允许定义的组件型号分为5种,分别为燃料组件、探测器组件、分离式可燃毒物组件、控制棒组件和阻力塞组件,其中,探测器组件、分离式可燃毒物组件、控制棒组件、阻力塞组件称为相关组件;将每种组件型号的几何结构尺寸和材料构成参数保存在“型号文件”中,所有组件型号文件构成组件型号数据库;
步骤2:用户按照建立核反应堆堆芯方案的需求,依据步骤1给出的组件型号文件作为模板,定义实际摆放到堆芯中的组件,并逐一命名;
步骤3:用户将步骤2定义的实际摆放到堆芯中的组件,按照命名,填写到堆芯布置中进行摆放;将燃料组件与相关组件摆放在同一个位置,意味着实际上是将同一位置的燃料组件和相关组件配插在一起;
步骤4:用户设置堆芯的运行参数,包括控制棒分组方式、控制棒移动、化学可溶硼浓度、堆芯压力、运行时间、功率水平、冷却剂入口温度和冷却剂入口流量;
步骤5:根据步骤3中实际摆放到堆芯中的组件的布置,完成组件实体插配;再针对堆芯内所有位置上的组件找出所有不同类型的组件横剖面(称之为燃料栅格),通过旋转和折叠对称性判断尽可能减少燃料栅格数量;然后,将径向反射层区域的横向剖面和轴向反射层区域的纵向剖面均划分成与燃料栅格等大的若干个区域,将每一个区域都作为一个反射层栅格,并与堆芯内使用次数最多的燃料栅格直接拼接在一起,构成反射层超栅格;最后,对所有燃料栅格和反射层超栅格,依据其几何参数和材料参数,生成二维栅格计算程序所需的输入文件;
步骤6:调用二维栅格计算程序,通过并发计算,完成所有的燃料栅格和反射层超栅格计算,并根据栅格计算中的工况点,选取燃耗深度Bu,硼浓度CB,燃料温度Tf,慢化剂温度Tm作为状态参数,以Σ=f1(Bu,CB,Tm)+f2(Bu,Tf)的形式进行少群参数参数化计算;将计算结果存为少群常数库;其中Σ为少群常数,f1和f2为待定系数的多项式;
步骤7:根据步骤3中组件实体的布置,生成堆芯稳态和瞬态计算软件所需的输入文件:根据堆芯中的组件布置和运行方式,结合步骤5选出的燃料栅格和反射层超栅格类型,得出三维堆芯各空间位置与少群常数库的对应关系,根据上述对应关系,生成堆芯稳态和瞬态计算软件所需的输入文件;
步骤8:调用堆芯稳态和瞬态计算软件进行计算,并完成计算结果的收集和整理,将计算结果绘制为图表,方便用户查看。
与现有技术相比,本发明有如下突出优点:
1.前台建模输入与后台计算方法实现完全脱耦,前台与工程实际保持高度一致。通过自动建模的方式实现了专家经验的信息化,用户建模时只需按照肉眼所见的实体完成对压水堆堆芯布置及运行需求的数据描述,即可开展堆芯计算分析,而无需了解具体计算方法的理论。
2.全自动高可靠性处理。用户只需要完成堆芯描述,之后的建模处理、内核程序大量输入文件的填写、所有内核程序的运行调用、中间数据文件的转移和拼接、最终计算结果的整理等所有操作,均由建模程序自动完成,提高了建模计算效率,降低了建模过程的出错率,保证了建模的可靠性。
3.全范围高精细建模。通过建模规则的程序化,实现了对全堆几何与材料的高精度建模处理。包括燃料组件中燃料棒与水隙的精确刻画,以及其他复杂结构精细建模,如格架与慢化剂材料的自动化打混处理,围板、吊篮、中子屏蔽体等结构的自动化精细化建模等。减低建模复杂度的同时保证了建模的精度,能更好的满足工程实际需求。
总的来说,本发明是一种高效可靠、可满足现有压水堆两步法计算程序建模计算的方法。
具体实施方式
下面结合具体实施方式对本发明作进一步详细说明:
一种压水堆堆芯自动建模方法,包括如下步骤:
步骤1:用户建立或补充完善组件型号库:参照核电厂购买燃料组件时使用的类型,定义组件“型号”,并按照实际需求,把允许定义的组件型号分为5种,分别为燃料组件、探测器组件、分离式可燃毒物组件、控制棒组件和阻力塞组件。其中,探测器组件、分离式可燃毒物组件、控制棒组件、阻力塞组件称为相关组件。将每种组件型号的几何结构尺寸和材料构成参数保存在“型号文件”中,所有组件型号文件构成组件型号数据库。
步骤2:用户按照建立核反应堆堆芯方案的需求,依据步骤1给出的组件型号文件作为模板,定义实际摆放到堆芯中的组件,并逐一命名;
步骤3:用户将步骤2定义的实际摆放到堆芯中的组件,按照命名,填写到堆芯布置中进行摆放。可以将燃料组件与相关组件摆放在同一个位置,意味着实际上是将同一位置的燃料组件和相关组件配插在一起,比如燃料组件中配插分离式可燃毒物组件,或者燃料组件中配插阻力塞组件。
步骤4:用户设置堆芯的运行参数,包括控制棒分组方式、控制棒移动、化学可溶硼浓度、堆芯压力、运行时间、功率水平、冷却剂入口温度、冷却剂入口流量。
步骤5:根据步骤3中实际摆放到堆芯中的组件的布置,完成组件实体插配,如将分离式可燃毒物组件插配到燃料组件实体的导向管中,得到带毒物的燃料组件;再针对堆芯内所有位置上的组件找出所有不同类型的组件横剖面(称之为燃料栅格),通过旋转和折叠对称性判断尽可能减少燃料栅格数量;然后,将径向反射层区域的横向剖面和轴向反射层区域的纵向剖面均划分成与燃料栅格等大的若干个区域,将每一个区域都作为一个反射层栅格,并与堆芯内使用次数最多的燃料栅格直接拼接在一起,构成反射层超栅格;最后,对所有燃料栅格和反射层超栅格,依据其几何参数和材料参数,生成二维栅格计算程序所需的输入文件。
步骤6:调用二维栅格计算程序,通过并发计算,完成所有的燃料栅格和反射层超栅格计算,并根据栅格计算中的工况点,选取燃耗深度Bu,硼浓度CB,燃料温度Tf,慢化剂温度Tm作为状态参数,以Σ=f1(Bu,CB,Tm)+f2(Bu,Tf)的形式进行少群参数参数化计算。将计算结果存为少群常数库。其中Σ为少群常数,f1和f2待定系数的多项式。
步骤7:根据步骤3中组件实体的布置,生成堆芯稳态和瞬态计算软件所需的输入文件:根据堆芯中的组件布置和运行方式,结合步骤5选出的燃料栅格和反射层超栅格类型,得出三维堆芯各空间位置与少群常数库的对应关系,根据上述对应关系,生成堆芯稳态和瞬态计算软件所需的输入文件。
步骤8:调用堆芯稳态和瞬态计算软件进行计算,并完成计算结果的收集和整理,将计算结果绘制为图表,方便用户查看。
Claims (1)
1.一种压水堆堆芯自动建模方法,其特征在于:包括如下步骤:
步骤1:用户建立或补充完善组件型号库:参照核电厂购买燃料组件时使用的类型,定义组件“型号”,并按照实际需求,把允许定义的组件型号分为5种,分别为燃料组件、探测器组件、分离式可燃毒物组件、控制棒组件和阻力塞组件,其中,探测器组件、分离式可燃毒物组件、控制棒组件、阻力塞组件称为相关组件;将每种组件型号的几何结构尺寸和材料构成参数保存在“型号文件”中,所有组件型号文件构成组件型号数据库;
步骤2:用户按照建立核反应堆堆芯方案的需求,依据步骤1给出的组件型号文件作为模板,定义实际摆放到堆芯中的组件,并逐一命名;
步骤3:用户将步骤2定义的实际摆放到堆芯中的组件,按照命名,填写到堆芯布置中进行摆放;将燃料组件与相关组件摆放在同一个位置,意味着实际上是将同一位置的燃料组件和相关组件配插在一起;
步骤4:用户设置堆芯的运行参数,包括控制棒分组方式、控制棒移动、化学可溶硼浓度、堆芯压力、运行时间、功率水平、冷却剂入口温度和冷却剂入口流量;
步骤5:根据步骤3中实际摆放到堆芯中的组件的布置,完成组件实体插配;再针对堆芯内所有位置上的组件找出所有不同类型的组件横剖面称之为燃料栅格,通过旋转和折叠对称性判断尽可能减少燃料栅格数量;然后,将径向反射层区域的横向剖面和轴向反射层区域的纵向剖面均划分成与燃料栅格等大的若干个区域,将每一个区域都作为一个反射层栅格,并与堆芯内使用次数最多的燃料栅格直接拼接在一起,构成反射层超栅格;最后,对所有燃料栅格和反射层超栅格,依据其几何参数和材料参数,生成二维栅格计算程序所需的输入文件;
步骤6:调用二维栅格计算程序,通过并发计算,完成所有的燃料栅格和反射层超栅格计算,并根据栅格计算中的工况点,选取燃耗深度Bu,硼浓度CB,燃料温度Tf,慢化剂温度Tm作为状态参数,以Σ=f1(Bu,CB,Tm)+f2(Bu,Tf)的形式进行少群参数参数化计算;将计算结果存为少群常数库;其中Σ为少群常数,f1和f2待定系数的多项式;
步骤7:根据步骤3中组件实体的布置,生成堆芯稳态和瞬态计算软件所需的输入文件:根据堆芯中的组件布置和运行方式,结合步骤5选出的燃料栅格和反射层超栅格类型,得出三维堆芯各空间位置与少群常数库的对应关系,根据上述对应关系,生成堆芯稳态和瞬态计算软件所需的输入文件;
步骤8:调用堆芯稳态和瞬态计算软件进行计算,并完成计算结果的收集和整理,将计算结果绘制为图表,方便用户查看。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210427670.6A CN114840988B (zh) | 2022-04-22 | 2022-04-22 | 压水堆堆芯自动建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210427670.6A CN114840988B (zh) | 2022-04-22 | 2022-04-22 | 压水堆堆芯自动建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114840988A CN114840988A (zh) | 2022-08-02 |
CN114840988B true CN114840988B (zh) | 2024-02-23 |
Family
ID=82566390
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210427670.6A Active CN114840988B (zh) | 2022-04-22 | 2022-04-22 | 压水堆堆芯自动建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114840988B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009134498A2 (en) * | 2008-02-11 | 2009-11-05 | Westinghouse Electric Company Llc | A methodology for modeling the fuel rod power distribution within a nuclear reactor core |
CN113094947A (zh) * | 2021-03-29 | 2021-07-09 | 西安交通大学 | 一种核反应堆堆芯核热耦合分析方法 |
CN113326648A (zh) * | 2021-06-17 | 2021-08-31 | 中国核动力研究设计院 | 考虑环境效应的栅元均匀化常数计算方法、系统及终端 |
-
2022
- 2022-04-22 CN CN202210427670.6A patent/CN114840988B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009134498A2 (en) * | 2008-02-11 | 2009-11-05 | Westinghouse Electric Company Llc | A methodology for modeling the fuel rod power distribution within a nuclear reactor core |
CN113094947A (zh) * | 2021-03-29 | 2021-07-09 | 西安交通大学 | 一种核反应堆堆芯核热耦合分析方法 |
CN113326648A (zh) * | 2021-06-17 | 2021-08-31 | 中国核动力研究设计院 | 考虑环境效应的栅元均匀化常数计算方法、系统及终端 |
Non-Patent Citations (2)
Title |
---|
压水堆核电站功率控制系统对象研究及仿真分析;林桦;王云伟;张往锁;;自动化仪表;20170516(第05期);全文 * |
基于NECP-X程序的三维复杂几何小型压水堆全堆芯一步法计算;曹璐等;核动力工程;20181230;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114840988A (zh) | 2022-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2008157937A (ja) | 認可基準を満たす原子炉炉心のコンピュータ設計方法およびシステム | |
CN109409022B (zh) | 一种用于核电堆芯物理仿真可视化建模调试测试方法 | |
MX2007013337A (es) | Metodo y sistema para diseñar un nucleo de reactor nuclear para operaciones de potencia nominal ascendente. | |
Vojackova et al. | Safety analyses of reactor VVER 1000 | |
Zhang et al. | Development and verification of the coupled thermal–hydraulic code–TRACE/SCF based on the ICoCo interface and the SALOME platform | |
CN109359433B (zh) | 一种堆芯换料设计的自动化方法 | |
Rais et al. | Neutronics modeling of the CROCUS reactor with SERPENT and PARCS codes | |
CN114840988B (zh) | 压水堆堆芯自动建模方法 | |
CN107122564B (zh) | 一种压水堆pin-by-pin计算少群常数参数化方法 | |
Martínez-Quiroga et al. | NIRK3D and 3DKIN: General Description and current status of the new 3D kinetics capabilities of RELAP5/SCDAPSIM/MOD4. 0 | |
Lizorkin et al. | Development of codes and KASKAD complex | |
Bao | Development of the model for the multi-physics analysis of Molten Salt Reactor Experiment using GeN-Foam code | |
Kim et al. | Preliminary coupling of MATRA code for multi-physics analysis | |
CN113326648B (zh) | 考虑环境效应的栅元均匀化常数计算方法、系统及终端 | |
Litskevich et al. | Verification of the current coupling collision probability method with orthogonal flux expansion for the assembly calculations | |
DeHart et al. | Preliminary results for the OECD/NEA time dependent benchmark using Rattlesnake, Rattlesnake-IQS and TDKENO | |
Hegyi et al. | Coupling the ATHLET 3.0 and the KIKO3DMG multigroup 3D kinetic code developed for the fast spectrum Gen-IV reactors | |
Allelein et al. | First results for fluid dynamics, neutronics and fission product behavior in HTR applying the HTR code package (HCP) prototype | |
Zhang et al. | R&D Roadmap for Enhanced Resilient Plant Systems, Metrics, Scenarios, Risk Analysis, and Modeling and Simulation | |
CN113901029B (zh) | 一种支持堆芯多专业计算软件的统一建模方法及系统 | |
Armishaw et al. | The ANSWERS Code MONK–A new approach to scoring, tracking, modelling and visualisation | |
Stauff et al. | Updated status of the ART neutronic fast reactor tools integration to the NEAMS Workbench | |
Xhonneux et al. | The enhanced fuel management models of the HTR code package (HCP) module SHUFLE | |
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 | |
Bousquet | Development of advanced methods for safety assessment of sodium cooled fast reactors |
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 |