CARR自给能探测器实验设计与验证研究.pdf

517Nuclear Scienceand EngineeringJun.20232023年6 月No.3Vol.43工程第3期核科学与第43卷CARR自给能探测器实验设计与验证研究乔雅馨，陶?杰，吴小飞（中国原子能科学研究院，北京10 2 413）摘要：辐照靶件的释热计算对堆内实验设计影响重大。本文依据在CARR堆上开展的自给能探测器实验，分析了不同核评价数据库在辐射俘获反应释热计算上的差异来源，并对较早版本数据库的计算结果进行了一定的修正。实验结果表明，以ENDF/B-V.0 等库为代表的数据库计算结果更为合理，并建议研究者在进行后续堆内实验设计时，关注KERMA因子的能量平衡情况，优先采用数据更为完备的核评价数据库进行计算。关键词：自给能探测器；释热计算；核评价数据库；CARR中图分类号：TL32文章标志码：A文章编号：0 2 58-0 9 18（2 0 2 3）0 3-0 517-0 5Study on Self-Powered Neutron Detector ExperimentDesign and Verification in CARRQIAO Yaxin,TAO Jie,WU Xiaofei*(China Institute of Atomic Energy,Beijing 102413,China)Abstract:Heat releasing calculation for irradiated targets have a tremendous influence onin-core experiment design.Based on the self-powered neutron detector experiment in CARR,the difference of heat releasing calculation results on radioactive capture reaction with severalevaluated nuclear data libraries has been studied,and a modification to the calculation withearlier database edition has been implemented.The experimental results show that the heatreleasing calculation results with ENDF/B-V II.0 are more reasonable,and a evaluated nucleardata library edition with more data is always suggested to be used for future in-core experimentdesign.Users should check if the KERMA factors and energy balance are adequate or not toget reliable results.Key words:Self-powered neutron detector;Heat releasing calculation;Evaluated nuclear datalibrary;CARR收稿日期：2 0 2 2-0 2-2 8作者简介：乔雅馨（19 9 0 一），女，吉林乾安人，工程师，博士研究生，现从事反应堆物理方面研究基金项目：国家重点研发计划资助项目（2 0 0 2 YFB1902600）通讯作者：吴小飞，E-mail：w u x i a o f e i c i a e.a c.c n518自给能探测器能够连续测绘和监督活性区内部中子通量的分布和变化，还可以累积反应堆燃料组件的燃耗深度数据，适合在高温高湿和强腐蚀的堆芯环境中长期工作，近年来被广泛用作大型压水堆的堆内探测器，是堆芯仪表系统设备国产化任务中的重要一环。某新型压水堆设计采用（10 3Rh）自给能探测器，在正常运行工况下，探测器环境为中子注量率71013n/cms水平，温度350 水平，这些国产化探测器用于工程实践之前，需要通过实验获取关键参数，以证明其在反应堆实际测量条件下的核性能指标满足设计要求。1设计需求中国先进研究堆（CARR）是一座轻水冷却慢化、重水反射的反中子阱型高性能研究堆，CARR重水箱内设置有不同尺寸的2 1根垂直孔道，目前已开展多种同位素辐照及燃料、材料辐照考验应用。堆内辐照实验开展前，首先应根据实验需求完成靶件的物理热工分析，确定堆内辐照位置及配套结构设计，通过反复迭代优化，在得到符合技术需求的设计方案的同时，保证反应堆的运行安全与辐射安全。CARR堆垂直孔道内介质为轻水，为了模拟受试探测器在压水堆内的高温高通量环境，首先通过计算选择满足通量需求的实验位置，之后在选定孔道内架设干孔道，以铝支架作为辅助热源，通过调整铝材质量来满足温度需求。2计算方法辐照靶件的释热，包括中子碰撞热和光子热，以及生成核素的衰变热等 2 ，其中，中子和沉积是靶件释热的主要来源，通常采用蒙特卡洛程序MCNP计算。MCNP程序使用ACE格式的数据库，通常由NJOY程序对ENDF格式的评价核数据库处理而成。由于CARR堆在设计阶段大量采用的是ENDF/B-VI.8数据库，因此，在实验设计计算时，首先选用ENDF/B-VI.8库进行分析，继而选用其他版本数据库进行验证。对受试探测器及实验场所的几何建模如图1和图2 所示。根据技术需求，CARR堆可提供的垂直孔道为一根70mm的同位素辐照湿孔道，在该孔道内架设一根底部封闭、顶部通大气的干孔道，作为辅助热源的铝支架设计为双层，受试探测器固定在铝支架上，探测器的外径约为2mm，包括了发射体、绝缘体和收集体。堆内孔道外壁轻水一千孔道外壁空气双层铝支架受试探测器热电偶图1实验孔道内径向布置图Fig.1The radial schematic of experiment channel受试探测器/热电偶电缆堆内孔道外壁轻水千孔道外壁双层铝支架空气图2 实实验孔道内轴向布置图Fig.2The axial schematic of experiment channel3计算结果分析与验证3.1六种数据库的计算结果比较分别计算了实验场所处受试探测器和铝支架的释热率。除受试探测器发射体材料（10 3Rh）和铝支架材料（2 7 A1）外，计算中用到的其他材料均选用ENDF/B-VI.8数据库，10 3Rh和2 7 A1则分别采用ENDF/B-VI.8、END F/B-VII.0 3、CEND L-3.2 4 、JEFF-3.35、JEND L-4.0 6 和BROND-3.17等六种国际上主要的评价核数据库。六种数据库计算的铝支架释热率基本一519V栅元体积；径迹长度；粒子权重；其中：W-致，而受试探测器中10 3Rh的释热率分为明显的两类，具体如表1所示。经过比较，与两类结果对应的数据库共同的区别在于，相比 ENDF/B-VI.8、CENDL-3.2 和 JENDL-4.0 库，ENDF/B-VII.0、JEFF-3.3和BROND-3.1库中的Rh-103材料增加了光子产生数据，包括文件12(photon production yielddata）和文件 14（photon angular distributions)。表1不同数据库计算13Rh释热率Table110Rh heat releasing calculation withdifferent evaluated nuclear data libraries核评价库名计算释热率/WENDF/B-VI.8351.38CENDL-3.2348.96JENDL-4.0336.12ENDF/B-VII.07.48JEFF-3.38.08BROND-3.17.36由于两类数据库Rh-103释热率的计算差异巨大，而实验可接受的温度范围无法覆盖如此大的差异，如果计算结果比实际释热值偏大，按照计算结果设计的方案就会导致实验失败；而如果计算结果比实际释热值偏小，按照计算结果设计的方案则可能导致实验超温，危及反应堆的运行安全。3.2差异来源分析为了分析数据库带来的计算差异，下面分别计算受试探测器处的中子注量率、光子注量率以及10 3Rh的中子释热率和光子释热率。计算结果如表2 所示（未按反应堆功率归一）。表2#不同数据库计算中子/光子注量率及13Rh释热Table 2 Flux&13Rh heat releasing calculationwith different evaluated nuclear data libraries核评价库名ENDF/B-VI.8ENDF/B-VII.0中子注量率3.33 10-53.35 10-5光子注量率2.23 1055.94 10-5中子释热率1.64 10-46.49 10-10光子释热率1.44 1063.63 106注：计算结果为MCNP自动归一至“单个初始粒子”，不具备实际物理意义。从计算结果可以看出，不同数据库计算的中子注量率基本一致，ENDF/B-VI.8库的光子注量率要略小于ENDF/B-VI.0库的光子注量率，与此对应的 ENDF/B-VI.8库计算 10 3Rh 光子释热率也略小于ENDF/B-VII.0库的光子释热率。但是，由于 ENDF/B-VI.8 库的中子释热率比ENDF/B-VII.0库的中子释热率大了近6 个量级，导致在最终的计算结果（表1）中，ENDF/B-VI.8库的总释热率比ENDF/B-VII.0库大了近2 个量级。在 MCNP 中，注量率和释热率是通过式(1）式（2）计算的,JaejJaifayJaow(,o,E,t)=w工(1)VH,=Pa J E di dy J do,(E)H(E)(r,0,E,t)=WTo,(E)H(E)Pam(2)Q:(E)-微观总截面；H(E)-MCNP定义的释热数（heatingnumber);P原子密度;m栅元质量。从计算公式以及计算结果对比可以看出，造成中子释热差异的来源可能有两项，一个是反应总截面(E)，另一个是释热数 H(E)。对于不同版本的数据库，由于有微观实验数据支持，反应总截面一般差异不大。图3绘制了不同版本数据库反应总截面随能量变化的关系图。B68C32B80F3310001001010-410-2110210410610510-310-110103105107能量/eV图3不同版本数据库反应总截面Fig.3The total reaction cross-section for differentevaluated nuclear data libraries520对于中子，其释热数H(E)的计算公式为，H(E)=E-Z,P,(E)Ei,out(E)-Q,+Ei,(E)(3)其中:p;(E)一对应入射能量E的反应i发生的概率；Eiou(E)一反应i的平均出射中子能量；9,一反应i的Q值;Ei,(E)一反应i的平均出射能量。通过比较得知，MCNP中定义的中子释热数 H(E)正是 KERMA（K i n e t i c En e r g y R e l e a s e dinMaterials）因子k,(E)的变体 8 ，其通常用于计算混合物中的释热率，在辐照损伤领域有大量应用。KERMA因子可以视为另一种微观反应截面，单位为能量截面（eV.barns），有两种计算办法，一种是“直接法”，k,(E)=Z,Em(E)0,(E)(4)其中：i一混合物中的材料种类；j一反应种类;E一被第1种次级粒子带走的总动能。由于求和需要对包括反冲核在内的所有带电产物进行计算，而早期评价数据库中并不包含如此详细的谱信息，因此，可以通过“能量平衡方法”来计算KERMA因子。能量平衡方法的基本思想是，从可获得的能量中减去分配给中子和光子的能量即可得到由带电粒子带走的能量，计算公式如下：k,(E)=(E+Qj-Ejn-Ejir)o,(E)(5)其中：E一一次级中子的总能量；E，一包含了光子产额的次级光子的能量。NJOY在HEATR模块采用“能量平衡方法”计算KERMA因子，给出MCNP计算所需的中子释热数H(E)，并储存在连续能量截面库ACE文件的ESZ部分 9 。图4给出由NJOY计算的不同数据库的中子释热数H(E)，可以看出，其差值正是6 个量级，与中子释热率的计算差异一致。造成这种差异的来源是，ENDF/B-VI.8等数据库中由于缺少光子产生数据，即NJOY在计算KERMA因子时，将公式（5）中的Ej视为零，导致k,计算结果偏大，在物理上，这相当于假设所有的光子能量不经过输运均沉积在本地，使得本应由光子带走的能量计入到了中B68C32B80F3310-110-210-31010一1010-710-1010-810610-410-211010-1110910-710510-310-1能量/MeV图4日由NJOY计算的不同数据库的中子释热数Fig.4 Heating number calculated by NJOY withdifferent evaluated nucleardata libraries子释热中去。事实上，NJOY手册中也提到，对于辐射俘获（n，）反应，由于评价的俘获谱中一般不包括出射光子对反冲组分的影响，因此，很多KERMA因子会计算为负值或者极大值。对于负数的KERMA因子，NJOY的开发者曾给出解释，认为虽然计算结果为负，但整体上保持了能量平衡 8 。在之前的工作中，胡也等人介绍了不同版本ENDF/B核评价数据库中Ag、In 等核素中子释热率截面差异，建议基于ENDF/B-VI的控制棒释热率计算值更加具有工程参考意义 10 ,KonnoC等人检查了基于ENDF/B-VI.1计算的KERMA因子，发现超过2 0 0 种核素存在计算值不合理的情况，这些核素的KERMA因子为负值或者极大值 11。赵秋娟等人研究了中子KERMA因子的计算，认为由于评价数据及其表达不合理，能量不平衡问题普遍存在于各评价库中，并提出可以用计算的运动学上限来替代出现能量平衡问题的KERMA因子 12 。因此，本文采用NJOY计算ENDF/B-VI.8库的KERMA因子及其运动学上下限，并用运动学上限值来代替KERMA因子，更新ACE文件中的中子释热数H(E)，重新计算受试探测器的释热率。计算结果如表3所示（未按反应堆功率归一），可以看出，中子释热计算结果得到了明显的改善。但是，由于中子释热在总释热率中的占比很小，修正KERMA因子仅可以使521得计算结果回到更为合理的量级，而由于光子释热仍未得到修正，因此在设计实验方案时需要留出更多的安全裕量。光子释热计算结果的改善依赖于微观评价数据的修订。3.3实验结果验证确定因评价数据库带来的计算差异原因后，依据ENDF/B-VII.0库进行了实验方案设计，并完成了受试探测器在 CARR 堆上的核性能测试。在设计方案中，受试探测器所处位置中子注量率环境为7 10 13n/cms，内层支架的计算温度为342 36 1（轴向温度分布范围），外层支架的计算温度为2 8 52 9 2（轴向温度分布范围）。受试探测器实验期间，最高中子注量率达到了8 10 13n/cms水平，最高温度达到了40 0 水平，达到了设计阶段需求的全部技术指标，并验证了ENDF/B-VI.0库对于10 3Rh释热计算结果的合理性。表3修正释热数后计算中子/光子注量率及10 3Rh释热Table 3Flux&1Rh heat releasing calculation withdifferent evaluated nuclear data libraries aftermodifying heating number核评价库名ENDF/B-VI.8中子注量率3.280 2810-5光子注量率2.2174010-5中子释热率6.164 95 10-9光子释热率1.382 60 10-6注：计算结果为MCNP自动归一至“单个初始粒子”，不具备实际物理意义。4总结与结论本文介绍了采用不同核评价数据库在CARR堆进行自给能探测器辐照的实验设计与验证，分析了来源于核评价数据库的计算差异产生的原因，并对计算结果进行了一定的修正。实验验证结果表明，在CARR堆垂直孔道内架设干孔道，采用合适质量的铝支架作为辅助热源，可以模拟大型压水堆高温高通量环境，开展国产自给能探测器的核性能考验实验。CARR堆已开展多种同位素辐照生产，燃料、材料辐照考验应用，这类应用涉及的核素种类多，范围广，且多数不属于重要裂变产物，在早期版本的评价库中有大量数据缺失，核反应释热计算不准确可能是一个共性问题。因此，建议研究者在进行后续堆内实验设计时关注KERMA因子的能量平衡情况，优先采用数据更为完备的核评价数据库进行计算，也建议评价者对当前应用较多的核素，如V、R h、Lu 等的数据进行进一步修订。参考文献：1王文滋，孙景海，杨友琏堆芯自给能中子通量探测器J.核技术，19 8 0，1：1-14.2沈峰。堆内辐照过程中辐照靶件的核发热和传热研究J.原子能科学技术，2 0 0 3，37：345-348.3BROWN D.A.,et al.,ENDF/B-VIl.0:The 8th MajorRelease of the Nuclear Reaction Data Library with CIELO-project Cross Sections,New Standards and ThermalScattering Data J.Nuclear Data Sheets,2018,148:1.4 Ge Zhigang,et al.,CENDL-3.2:The new version ofChinese general purpose evaluated nuclear data library C.EPJ.Web.Conf,2020:239,09001.5Plompen A.J.M.,et al.,The joint evaluated fission andfusion nuclear data library,JEFF-3.3J,Eur.Phys.J.A,2020,56:181.6SShibata K.,Iwamoto O.,T.Nakagawa et al.,JENDL-4.O:A New Library for Nuclear Science and Engineering JJ.J.Nucl.Sci.Technol,201l,48:1.7IgnatyukA.V.,Fursov B.I.,The latest BROND-3development C.Proc.Int.Conf.Nuclear Data forScienceandTechnology,2007:759.8 MacFarlane R.E.,et al.,The NJOY Nuclear Data ProcessingSystem,Version 2016RJ.LA-UR-17-20093,Los AlamosNational Laboratory,2016.9 Conlin Jeremy Lloyd,Romano Paul.A Compact ENDF(ACE)Format Specification RJ.LA-UR-19-29016,LosAlamos National Laboratory,2019.10胡也，陈义学，杨寿海，等基于MonteCarlo方法的压水堆相关组件内热源计算与分析 J.核科学与工程，2013,33:168-174.11J Konno C.,et al.,Remarks on KERMA Factors in ACE filesJ.NuclearData Sheets，2 0 14,118:450.12 赵秋娟，吴海成，葛智刚连续能量点截面库CENACE-1.0的能量平衡检查J原子能科学技术，2 0 14，48：9 6 8-9 7 3.