压水堆高保真换料循环计算功能开发与验证应用_王习宁.pdf
《压水堆高保真换料循环计算功能开发与验证应用_王习宁.pdf》由会员分享,可在线阅读,更多相关《压水堆高保真换料循环计算功能开发与验证应用_王习宁.pdf(7页珍藏版)》请在咨信网上搜索。
1、文章编号:0258-0926(2023)02-0030-07;DOI:10.13832/j.jnpe.2023.02.0030压水堆高保真换料循环计算功能开发与验证应用王习宁,刘宙宇*,周欣宇,温兴坚,曹璐,张思凡,许晓北,易思宇,李帅铮,李帆,苏鑫西安交通大学核科学与技术系,西安,710049摘要:在自主开发的数值反应堆物理计算程序 NECP-X 基础上开发了压水堆的换料循环计算功能,并针对某 M310 机组首循环、第 2 循环和第 3 循环的启动物理实验,以及针对前 2 个循环的燃耗进行了精细建模计算。计算值与实测值的比较结果表明:首循环、第 2 循环和第 3 循环启动物理实验的临界硼浓度
2、、控制棒价值、温度系数计算结果误差均较小,符合验收准则;不同燃耗深度下的临界硼浓度、堆芯功率分布与实测值的比较结果显示,稳定燃耗点处最大硼浓度偏差为39ppm(1ppm=106),最大的组件功率误差小于4.5%,随着燃耗的加深,堆芯功率的分布逐渐展平,误差逐渐减小。计算结果表明 NECP-X 程序已经具备商用压水堆启动物理实验和多燃料循环的计算能力。关键词:数值反应堆;NECP-X;多循环换料;大型压水堆验证中图分类号:TL326文献标志码:ADevelopment and Validated Application of Calculation Functionof High Fidelit
3、y Refueling Cycle for Pressurized Water ReactorWang Xining,Liu Zhouyu*,Zhou Xinyu,Wen Xingjian,Cao Lu,Zhang Sifan,Xu Xiaobei,Yi Siyu,Li Shuaizheng,Li Fan,Su XinDepartment of Nuclear Science and Technology,Xian Jiaotong University,Xian,710049,ChinaAbstract:The refueling cycle calculation function for
4、 pressurized water reactor (PWR)isdeveloped on the basis of the self-developed numerical nuclear reactor physics calculation codeNECP-X.Startup physics experiments are conducted for the first,second and third cycles of anM310 reactor,and fine modeling calculation is carried out for the first two cyc
5、les.By comparing thecalculated values with the measured values,it shows that the errors of calculation results of criticalboron concentration,control rod worth and temperature coefficient in the startup physicsexperiments for the first,second and third cycles are relatively small,which meet the acce
6、ptancecriteria.The results of comparison of the critical boron concentration and core power distributionwith the measured values at different burnup levels show that the maximum boron concentrationdeviation at the stable burnup point is 39ppm(1ppm=106),and the maximum assembly powererror is less tha
7、n 4.5%.With the increase of burnup level,the core power distribution flattens outand the error decreases gradually.The calculation results show that NECP-X already has thecalculation function for the startup physics experiments and multi-fuel cycle of commercial PWRs.Key words:Numerical reactor,NECP
8、-X,Multi-cycle refueling,Large PWR validation 收稿日期:2022-03-25;修回日期:2022-12-26基金项目:国家自然科学基金(11735011)作者简介:王习宁(1996),男,硕士研究生,现从事核反应堆物理研究工作,E-mail:*通讯作者:刘宙宇,E-mail: 第 44 卷第 2 期核 动 力 工 程Vol.44 No.22 0 2 3 年 4 月Nuclear Power EngineeringApr.2023 0 引言当前国内外可进行反应堆高保真计算的程序主 要 分 为:以 MCNP1、JMCT2、RMC3、SuperMC4为代
9、表的蒙特卡洛程序和以DeCART5、MPACT6、NECP-X7 为代表的确定论程序。这些程序大多实现了对启动物理实验和首循环的模拟计算8-9,但对第 2 循环和后续循环的计算相对较少。为了更好地将高保真数值模拟计算应用于实际的堆芯设计和校验计算,需要开发并验证高保真数值模拟程序多循环堆芯验证计算功能。本文基于我国自主开发的数值反应堆物理计算程序 NECP-X 开发了压水堆的换料和多燃料循环计算功能,并通过某 M310 机组的首循环、第2 循环、第 3 循环模拟计算,给出了 M310 机组3 个循环的启动物理实验计算结果和首循环、第2 循环燃耗模拟计算结果,并与实测值进行了比较。1 NECP-
10、X 程序介绍NECP-X 程序是由西安交通大学核工程计算物理实验室研发的数值反应堆物理计算程序,是基于 69 群多群数据库和连续能量数据库,共振处理采用全局-局部耦合共振自屏计算方法10-11,通过 2 维/1 维耦合输运计算方法进行输运计算12。同时,采用粗网有限差分(CMFD)方法加速13,通过采用空间区域分解、角度区域分解与特征线并行的多重并行策略,提高计算效率,已实现了基于栅元的物理热工耦合和非均匀全堆瞬态计算14。NECP-X 程序基于高精度压缩燃耗数据库,利用切比雪夫有理近似方法求解点燃耗方程,并结合经典的预估校正方法和子步方法,实现了中子输运和点燃耗的耦合计算15。2 压水堆高保
11、真堆芯多循环计算功能的设计和开发 2.1 高保真计算换料与两步法计算换料的区别 2.1.1 高保真建模对换料计算的影响两步法基于均匀的节块进行堆芯计算,每个燃料组件分为 4 个均匀的节块,因此基于节块实现换料,映射关系比较简单。不同循环同一位置的节块大小一致,直接对节块燃耗、温度等进行赋值即可。高保真一步法计算换料的映射关系比较复杂,对于高保真而言,并没有像传统两步法一样的均匀化过程,而是详实描述真实模型的各部分,如燃料、包壳、慢化剂、吸收体和格架等精细结构,不同类型的栅元(如燃料棒和导向管)拥有不同的材料区数目和材料编号。由于换料之后,组件的位置会发生变化,同一栅元位置不同循环的组件网格和材
12、料数目不同,导致换料的映射关系变得复杂。2.1.2 高保真循环计算换料文件设计与功能设计两步法计算效率和内存要求较低,一般直接串行计算的效率已经能够满足要求。但是对于高保真而言,需要通过大规模区域分解并行才能有效解决计算时间和内存的问题,但区域分解并行使得不同中央处理器(CPU)核心拥有不同的映射关系,而彼此之间这些信息并不“共享”,由此给换料计算带来了问题:不同循环区域分解方式可能不同,使用的 CPU 核心数不同,几何网格映射关系随之改变。为了适应更灵活的并行方案,换料计算时不同循环的网格映射关系变得非常复杂。2.2 高保真循环计算换料文件设计与功能设计高保真计算时,一般以模块(即 1/4
13、组件,且按照 4 个模块在整个组件中的方位,分别为西北-WN、东北-EN、西南-WS、东南-ES)为基本单元。针对高保真多循环计算,设计了其换料功能,一个燃料循环计算结束后,每个 CPU 核心将其计算区域中每个模块内所有材料区的核子密度输出至换料文件。换料计算时,程序建立换料之后的堆芯几何模型,并进行并行区域分解,每个计算核心负责全堆芯的部分模块。然后,每个 CPU 核心读取换料文件,获取对应模块每个材料区的核子密度。其中换料文件的设计、换料后几何模型建立和换料文件的读取是整个换料过程的核心环节。2.2.1 换料文件生成高保真计算的换料文件数据量很大而且众多的材料信息分别存储在不同的 CPU
14、核心中。以某 M310 机组为例,堆芯燃料组件 157 盒,全堆芯的燃耗区(各层的燃料)总数达 2479680,总材料区总数为 14678776,假设每一个燃料区的核素数目为 200 个,单个燃耗点的换料文件大小达 9.7 GB。高保真计算时每个 CPU 核心只包含了各自计算区域的信息,并行的 HDF5 格式输出将会显著降低输出效率。王习宁等:压水堆高保真换料循环计算功能开发与验证应用31 以 M310 机组为例,采用 2652 个 CPU 核心在线生成换料的 HDF5 文件,花费时间数小时。如果每个 CPU 单独输出换料文件,将对后续计算重新建立映射关系带来极大的困难。为了解决上述问题,本文
15、提出了并行在线输出、离线合并换料文件的方法,其具体操作过程为:在每个燃耗点,每个 CPU 核心输出该 CPU核心包含区域的材料文件,包括核素号、核子密度和温度等信息,以及该材料所处模块的燃料组件的全局位置坐标信息;离线将每个 CPU 核心输出的换料文件合成为一个文件,供换料时读取。同样针对 M310 机组问题,每个 CPU 核心输出核子密度文件时间为 11.3 s,单核离线合并文件需要 40 min,但由于仅在停堆的燃耗点进行合成,并且不需要占用大量的并行资源,总体效率和资源利用率很高。2.2.2 换料功能设计在高保真计算的中子学建模中,每个组件的材料区和几何网格划分是严格对应的,但是在换料之
16、后,不同循环中相同位置的组件网格划分可能不一致,因此需要建立新的网格映射关系,通过保持同一个燃料组件在不同燃料循环的网格不变性以实现换料过程中各个燃耗区和计算网格的映射,本文设计并建立了换料计算时网格的划分方式,首先读入所有几何信息,根据当前循环计算参数设置新料的网格,对于旧料,从换料文件中获取该燃料组件第一次网格划分时的参数进行网格划分,从而保证该循环计算时的网格与旧料的原始网格信息一致。程序计算时,在区域分解完成后开始读入换料文件。各个 CPU 核心均直接从换料文件获取各个 CPU 核心需要的材料信息,避免出现单个计算核心读入全部换料文件导致内存瓶颈问题。3 数值结果与分析 3.1 问题描
17、述与建模本文将高保真多循环计算功能应用于某M310 机组的首、第 2、第 3 循环计算分析中。M310 机组首循环中,有 1.8%、2.4%、3.1%3 种富集度,第 2 和第 3 循环中加入富集度 3.7%的新燃料组件。轴向活性区分为 52 层,上下反射层分别为 9 层和 7 层,同时对反射层进行显式建模,具体包括水隙、围板、吊篮等堆芯外围结构,轴向包括端塞、阻力塞、管座等结构,采用全堆芯建模计算。M310 机组的功率历史燃耗点较多,真实模拟每个燃耗点的计算成本非常高,本文计算时简化了功率历史,但保证计算时每个待分析比较的燃耗点的功率水平、控制棒棒位和实际测量时的情况完全一致,以确保计算结果
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 压水堆 高保真 循环 计算 功能 开发 验证 应用 王习宁
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。