删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

高能质子在散裂靶中的能量沉积计算与实验验证

本站小编 Free考研考试/2021-12-29

摘要:高能质子在散裂靶中的能量沉积是散裂靶中子学研究的重要内容之一, 准确掌握高能质子在散裂靶中引起的能量沉积分布与瞬态变化是开展散裂靶热工流体设计的重要前提. 本文采用MCNPX, PHITS与FLUKA三种蒙特卡罗模拟程序, 计算并比较了高能质子入射重金属铅靶、钨靶的能量沉积分布及不同粒子对总能量沉积的占比贡献; 针对高能质子入射金属钨靶的能量沉积实验数据空白, 采用热释光探测器阵列测量了250 MeV质子束入射厚钨靶的能量沉积分布, 实验结果表明蒙特卡罗模拟程序在散裂靶中能量沉积的计算结果具有较高的可靠性.
关键词: 能量沉积/
高能质子/
散裂靶/
蒙特卡罗程序

English Abstract


--> --> -->
近年来, 我国先后启动了中国散裂中子源(China Spallation Neutron Source, CSNS)项目[1], 以及未来先进核裂变能-加速器驱动次临界系统 (accelerator driven sub-critical system, ADS)项目[2]. CSNS项目采用1.6 GeV的高能质子束流轰击钽包钨靶[3], 在靶上发生散裂反应并产生大量中子, 为中子散射及相关领域提供科学研究与实验平台; 在ADS项目中, 高能质子轰击重金属靶, 为ADS次临界反应堆提供持续中子外源. 在散裂靶中, 高能质子主要通过与靶原子的激发电离过程, 以及与靶核发生散裂反应过程损失能量. 掌握高能质子在散裂靶中能量沉积的产生机制, 分析热量沉积的空间分布与瞬态变化, 是开展散裂靶热工流体设计的重要前提.
目前, 高能质子束在散裂靶中的能量沉积数据主要通过MCNPX[4], PHITS[5], FLUKA[6]等蒙特卡罗粒子输运程序计算获得, 模拟计算的准确度依赖于程序中物理模型对各个能量沉积途径的处理及采用的计算方法. MCNPX程序中的核子级联模型包括BERTINI, ISABEL, CEM2K和INCL4, 各个模型的主要区别在于对级联相互作用过程与对核密度等的近似处理; PHITS采用JAM模型作为强子级联模型, 在低能端采用JENDL核数据; FLUKA采用PEANUT模型与ENDF核数据进行相关计算.
受实验条件如束流能量、束流时间等的限制, 高能质子在散裂靶中的能量沉积测量数据非常有限. 国际上采用量热法进行靶内能量沉积测量[7-12], 即通过测量靶内各处由于能量沉积导致的温度变化而获得能量沉积分布, 但是只有束流较强时才能使靶体温度上升到可准确测量的范围. 文献[13]表明在较低束流强度条件下, 在散裂靶内布置热释光探测器(thermoluminescence detector, TLD), 通过测量TLD的剂量分布可以分析散裂靶内的能量沉积. 本文采用MCNPX, PHITS与FLUKA三种蒙特卡罗程序, 首先计算了不同能量高能质子入射重金属铅靶、钨靶的总能量沉积及其在入射深度的分布, 并通过与已有的测量数据进行比较, 评估不同级联模型的计算结果; 其后计算了不同能量的质子在散裂靶中产生的不同类型粒子, 并获得各种粒子对总能量沉积的占比贡献; 最后利用兰州重离子加速器[14]提供的250 MeV质子束流, 采用TLD阵列测量250 MeV质子束入射钨靶的能量沉积分布, 测量结果验证了蒙特卡罗程序计算质子入射散裂靶能量沉积的可靠性.
2
2.1.计算方法
-->高能质子在散裂靶中的能量沉积计算分别采用MCNPX, PHITS和FLUKA三种粒子输运程序开展. 在采用MCNPX程序进行计算时, 分别采用了BERTINI模型、ISABEL模型、CEM2K模型与INCL4模型作为级联模型. 在采用PHITS程序和FLUKA程序进行计算时, 则分别采用了JAM级联模型与PEANUT模型. 计算所选用的散裂靶模型为直径为20 cm, 长度为50 cm的圆柱体. 铅的密度为11.3 g/cm3, 同位素204Pb, 206Pb, 207Pb和208Pb的质量丰度分别为1.4%, 24.1%, 22.1%和52.4%; 钨的密度为19.0 g/cm3, 同位素180W, 182W, 183W, 184W和186W的质量丰度分别为0.12%, 26.5%, 14.31%, 30.64%和28.43%. 为了对比于国际上已有的实验数据, 验证模拟计算方法的可靠性, 质子能量分别为800, 1000, 1200 MeV, 从圆柱体靶端面垂直入射, 质子在入射横截面上的抽样概率服从二维高斯分布, 半高宽均为2.4 cm.
2
2.2.质子入射铅靶的能量沉积计算结果
-->图1(a)给出了MCNPX, PHITS与FLUKA三种粒子输运程序对入射能量为800, 1000, 1200 MeV的质子入射铅靶的总能量沉积. 同时, 图1(a)还给出了Belyakov-Bodin等[7]采用热电偶法在铅靶中测量到的能量沉积实验值. 可以看出, 三种粒子输运程序的计算结果相对于实验值整体偏大; 与MCNPX相比, PHITS与FLUKA的计算结果更大; 即使采用MCNPX程序进行计算, 不同的级联模型所对应的模拟结果也各不相同: 采用ISABEL模型的计算值较大, INCL4与BERTINI模型次之, CEM2K模型的计算结果最小. 对1000 MeV的质子入射能量而言, PHITS的模拟值比测量值高约16.1%, FLUKA的模拟值比测量值高15.6%, MCNPX-CEM2K的模拟值比测量值高约9%, 采用CEM2K级联模型获得的能量沉积值最接近实验测量值, 计算结果的统计误差可忽略不计.
图 1 (a)不同模拟程序对铅靶总能量沉积计算的对比; (b) 铅靶中能量沉积线性密度的轴向分布
Figure1. (a) Comparison of total energy deposition in lead target calculated by different Monte Carlo codes; (b) axial distribution of linear density of energy deposition in lead target.

散裂靶单位长度中的能量沉积值被称为能量沉积线性密度. 图1(b)给出了采用MCNPX-CEM2K模型计算出的能量沉积线性密度在铅靶中的分布. 通过与文献[7]中的测量数据进行比较, 可以看出, 高能质子在铅靶长度方向的能量沉积线性密度计算结果与测量结果在整体上能较好符合; 随着质子入射深度的增加, 能量沉积线性密度计算值的偏差略有增大.
另外, 我们采用MCNPX程序中的CEM2K级联模型计算得到了不同能量质子入射铅靶时, 不同粒子对总能量沉积值的占比贡献, 如表1所示. 计算数据表明, 质子作为入射粒子, 通过与靶核外电子的电离、激发过程及与靶核的散裂反应过程损失能量, 对总能量沉积的贡献最大; 铅是原子序数较高的材料, 对光子的阻止本领很强, 这导致了质子与散裂靶相互作用产生的光子在靶中的能量沉积较大; 随着入射质子能量的增大, 散裂反应产生的次级带电粒子增多, 预平衡与蒸发过程的剩余核会达到更高的自旋和激发态, 这导致了质子对总能量沉积的占比逐渐减小, 轻带电粒子与光子的占比贡献逐步增大.
粒子800 MeV1000 MeV1200 MeV
沉积能
量/MeV
对总能量沉积值
的占比/%
沉积能
量/MeV
对总能量沉积值
的占比/%
沉积能
量/MeV
对总能量沉积值
的占比/%
全部粒子497.9100572.0100648.7100
质子413.983.13435.676.15455.070.15
光子47.59.5474.813.08104.016.03
π00.20.040.40.060.50.08
带电π介子(± π)6.61.3313.72.3921.93.37
13.72.7521.23.7129.04.47
5.01.018.51.4812.31.90
氦–32.80.575.10.908.01.23
α粒子6.11.2310.11.7714.62.26
中子2.00.402.70.473.30.51
初级质子
电离作用
268.852.73247.948.63233.145.72


表1CEM2K级联模型计算质子入射铅靶产生的不同粒子对总能量沉积的占比贡献
Table1.The calculated contribution of different particles to the total energy deposition in lead target by CEM2K-Cascade-Mode

为了对比不同级联模型的次级粒子计算结果, 我们还采用MCNPX中的BERTINI, ISABEL, INCL4级联模型分别对1000 MeV质子进行了模拟, 计算结果如表2. 从表2可以看出, 4种级联模型对初级质子电离作用的能量沉积计算结果几乎相等, 这是因为初级质子的电离作用仅受到入射粒子和散裂靶的影响, 与级联过程无关; 相对其他两种级联模型, CEM2K对次级质子沉积能量计算值最小, 对光子、氘、氚、氦-3的能量沉积计算值最大, 这是因为CEM2K融合了DUBNA级联模型、Exciton激子模型、GEM蒸发模型、RAL裂变模型、Fermi-breakup模型及光核反应模型等, 在计算剩余核、裂变产物、轻核产生及碰撞-反弹过程时更为准确[15].
粒子BERTINIISABELCEM2 KINCL4
沉积能
量/MeV
对总能量沉积值
的占比/%
沉积能
量/MeV
对总能量沉积值
的占比/%
沉积能
量/MeV
对总能量沉积值
的占比/%
沉积能
量/MeV
对总能量沉积值
的占比/%
全部粒子580.3100603.0100572.0100.00594.5100
质子474.081.67493.081.76435.676.15500.684.20
光子63.710.9876.212.6474.813.0860.710.21
π00.30.050.30.050.30.060.30.04
带电π介
子(± π)
14.92.5613.92.3113.72.3915.52.60
6.41.103.70.6221.23.714.00.67
3.00.521.90.328.51.481.80.30
氦–30.40.060.10.025.10.900.20.03
α粒子14.82.5611.81.9610.11.779.11.53
中子2.90.491.90.322.70.472.50.42
初级质子电
离作用
245.148.06243.647.79247.948.63245.048.06


表2BERTINI, ISABEL, CEM2K与INCL4级联模型计算1000 MeV质子入射铅靶产生的不同粒子对总能量沉积值的占比贡献
Table2.The calculated contribution of different particles to the total energy deposition in lead target by 1000 MeV protons with BERTINI, ISABEL, CEM2K, and INCL4 cascade mode.

2
2.3.质子入射钨靶的能量沉积计算结果
-->同样, 对不同能量的质子入射钨靶产生的能量沉积进行计算, 图2(a)为采用不同计算模型获得的总能量沉积随入射质子能量的变化. 与铅靶的情况类似, PHITS与FLUKA的计算结果相对于MCNPX各个级联模型偏大. 图2(b)为采用MCNPX-CEM2K模型计算的能量沉积线性密度分布. 在现阶段, CSNS项目采用了钽包钨作为散裂靶, ADS项目则采用钨镍合金靶, 然而, 由于国内外尚无高能质子入射钨靶产生的能量沉积实验数据, CSNS与ADS的热工流体设计只能根据计算数据开展. 对于中高能质子入射钨靶的能量沉积计算亟需实验数据验证和可靠程度评估.
图 2 (a)不同模拟程序对钨靶总能量沉积计算的对比; (b)钨靶中能量沉积线性密度分布
Figure2. (a) Comparison of total energy deposition in tungsten target calculated by different Monte Carlo code; (b) axial distribution of linear density of energy deposition in tungsten target.

2
3.1.实 验
-->目前, 国内仅有兰州重离子加速器国家实验室的高能质子实验平台可提供能量为250 MeV的质子束流, 流强约为107 proton/s. 由于束流强度偏低, 我们采用TLD[13,16]的测量方法, 通过测量钨靶中TLD的剂量分布, 从而获得能量沉积数据, 并与模拟计算的结果进行对比.
质子在钨靶中的能量沉积的验证实验示意图如图3所示. 整个钨靶厚度为4 cm, 由厚度分别为1.0, 1.0, 0.5, 0.5, 0.5, 0.5 cm的几个靶块组成. 实验使用的TLD是由防化研究院生产的GR100系列, 其主要成分是LiF(Mg, Ti). 在钨靶的质子入射端面放置了一层TLD, 用于分析质子束流强度与分布. 在钨靶靶块之间布置了若干TLD, 用于分析钨靶不同深度处的能量沉积分布. 为了固定TLD的位置, 加工制作了1 mm厚的硬铝支架, 2个TLD的间距是3 mm. 置于钨靶内部的TLD大多数采用“十”字分布. 为了降低加速器输出质子脉冲稳定性和重复性影响, 实验仅接受了1个质子脉冲束团的辐照, 这也能确保所有TLD都能工作在剂量线性响应区域.
图 3 质子在钨靶中的能量沉积测量示意图
Figure3. Schematic of the energy deposition measurement in a tungsten target incident by protons.

经过辐照的TLD采用RGD6型热释光仪[17]进行数据读出. RGD6型热释光仪通过测量TLD的积分光产额并与标准60Co伽玛源剂量的标定光产额进行对比, 其输出量的单位是μGy. 第一层TLD的剂量读出结果如图4所示. 可以看出, 在钨靶几何中心水平面偏左下方区域的TLD剂量读数值较大, 这个区域对应着加速器输出质子束流的直接照射位置. 根据第一层TLD测得的质子直射区域, 结合图3的实验布置, 可以看出仅有 (0, –1) 位置的TLD在质子入射深度的各层对应位置有测量点. 我们采用MCNPX对单个能量为250 MeV的质子直射TLD时的总能量沉积开展了模拟计算. 具有不同线性能量转移值(linear energy transfer, LET)的带电粒子对TLD照射相同剂量后, 相对于等同剂量的60Co伽玛源射线照射, TLD的积分发光量会发生变化, 这被称为TLD的相对发光效率[18]. 第一层TLD被能量为250 MeV的质子直接入射, 对应相对发光效率值[19]约为1.02. TLD测量值的单位是μGy, 计算值的单位是MeV/g, 通过单位转换与发光效率修正, 对比得到在本次实验中单个质子脉冲辐照钨靶时, 到达(0, –1)位置TLD有效面积上的质子数目约为4.74 × 106个.
图 4 第一层TLD的剂量读出值
Figure4. TLD dose readouts at the first layer.

2
3.2.实验结果分析
-->为了获得不同深度位置TLD的相对发光效率值, 利用SRIM[20]程序计算了质子在穿透不同厚度的钨靶后的剩余能量, 并转换为等效水LET值如图5所示. 根据文献[19], 当入射粒子的等效水LET值在0.4—1.2 keV/μm范围时, LiF(Mg, Ti)相对发光效率约为1.02.
图 5 不同深度的钨靶中质子的平均能量与水等效LET值
Figure5. Average proton energy in the tungsten target and equivalent LET in water.

由于质子是重金属钨靶中能量沉积的主要贡献者, 结合TLD有效面积接受的入射质子数目与相对发光效率, 可以将TLD剂量测量值μGy进行单位转换, 得到单个质子在TLD中引起的能量沉积值, 即MeV/(g·proton).
我们选用MCNPX-CEM2K级联模型对 (0, –1) 位置不同质子入射深度的TLD进行了模拟计算. 首先, 根据图3中测量实验布置的几何条件建立计算模型; 然后, 利用程序中的能量沉积记录卡对分布在钨靶内部不同质子入射深度的TLD中总能量沉积进行了记录. TLD能量沉积模拟计算值与实验测量值, 以及其二者的比较如图6所示.
图 6 钨靶中TLD的能量沉积测量值与计算值
Figure6. Energy deposition comparison between measurement and calculation of TLD in tungsten target.

可以看出, 总的说来, MCNPX-CEM2K模型能够较好地模拟钨靶中热量沉积的结果. 通过细致地对比可以发现, 在钨靶较浅位置, 实验测量值比模拟计算值略小; 随着钨靶深度的增加, 测量值与模拟值的偏差逐步增大; 在质子射程末端位置, 二者相差最大. 出现较大差别的原因有以下几点: 首先, TLD的有效面积较小, 入射质子束细微的非垂直入射会对较深钨靶位置的测量值带来较大影响; 第二, 随着钨靶深度的增加, 较高LET的次级带电粒子对总能量沉积的占比贡献增大, 质子对总能量沉积的贡献份额减小, 这就需要对TLD的相对发光效率进行更细致的修正. 实验测量数据的误差为10%, 主要来源于质子束斑的不规则、同批次TLD出厂性能的差异、剂量读出设备的系统误差等.
本文分别采用MCNPX, PHITS和FLUKA三种蒙特卡罗粒子输运程序计算了高能质子在金属散裂靶中的能量沉积分布, 通过与已有实验测量结果的比较, 表明这三种蒙特卡罗程序能够对铅靶的能量沉积进行较准确的模拟计算, 相对于实验测量值, 模拟计算的结果高约10%—15%, 即蒙特卡罗程序的模拟计算结果能够比较保守地估计铅靶中的热量沉积. 对3个不同模拟程序及同一模拟程序的不同级联模型的计算结果进行比较, 表明MCNPX程序中的CEM2K级联模型的计算结果最接近实验测量值. 由于迄今国内外尚无高能质子在钨靶中的能量沉积测量数据, 我们利用热释光探测器LiF(Mg, Ti), 测量了能量为250 MeV的质子束在钨靶中的能量沉积分布, 结果表明除了在质子射程末端有较大的差异外, 能量沉积的实验测量结果与MCNPX的计算值有较好的吻合. 通过与实验数据的比较, 验证了蒙特卡罗模拟程序对能量沉积计算具有较高的可靠性.
相关话题/计算 测量 实验 程序 数据

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 外加横向电场作用下石墨烯纳米带电子结构的密度泛函紧束缚计算
    摘要:采用基于密度泛函理论的紧束缚方法计算研究了外加横向电场对边缘未加氢/加氢钝化的扶手椅型石墨烯纳米带的电子结构及电子布居数的影响.计算结果表明,石墨烯纳米带的能隙变化受其宽带影响.当施加沿其宽度方向的横向外加电场时,纳米带的能带结构及态密度都会产生较大的变化.对于具有半导体性的边缘未加氢纳米带, ...
    本站小编 Free考研考试 2021-12-29
  • 自由分子区内纳米颗粒的热泳力计算
    摘要:基于非平衡态分子动力学模拟方法,研究了自由分子区内纳米颗粒的热泳特性.理论研究表明,纳米颗粒与周围气体分子之间的非刚体碰撞效应会明显地改变其热泳特性,经典的Waldmann热泳理论并不适用,但尚未有定量的直接验证.模拟计算结果表明:对于纳米颗粒而言,当气-固相互作用势能较弱或气体温度较高时,气 ...
    本站小编 Free考研考试 2021-12-29
  • 船载系泊状态下基于原子重力仪的绝对重力测量
    摘要:重力场是反映地球质量分布及变化的重要参数,动态重力测量在地质调查、地球物理、资源勘探等领域有着重要应用.目前动态重力测量均基于相对测量原理,动态相对重力仪存在零点漂移问题,影响其测量性能.动态绝对重力仪可以为相对重力仪提供同步同址校准,解决其长漂问题,因此备受关注.本文基于原子重力仪和惯性稳定 ...
    本站小编 Free考研考试 2021-12-29
  • 量子真空计量标准中的非极性稀薄气体折射率测量研究
    摘要:为进一步提高真空量值的复现性和准确性,最新研究采用量子技术实现对真空量值的测量与表征.该方法利用Fabry-Perot谐振腔实现腔内气体折射率的精密测量,并反演出气体密度,进而获得对应的真空量值,其中气体折射率的测量是影响真空量值准确性的关键.本文基于第一性原理,利用从头计算理论计算了在已知压 ...
    本站小编 Free考研考试 2021-12-29
  • SiC表面圆环槽边缘效应实验研究
    摘要:基于固体边缘效应,对碳化硅(SiC)表面激光加工圆环形沟槽的润湿特性进行实验研究,通过分析去离子水在圆环槽上的润湿性能及其在边缘处的铺展行为,获得了环槽深度与环槽宽度对液滴在边缘处最大表观接触角的影响规律.结果表明,SiC圆环槽阻碍液滴铺展,光滑基体表面上接触角为70°,激光加工圆环槽深度为2 ...
    本站小编 Free考研考试 2021-12-29
  • 一种用于Z箍缩实验的软X射线成像系统
    摘要:基于塑料闪烁体转换和光学条纹相机的方法建立了一套用于Z箍缩实验中的软X射线条纹图像诊断系统,解决了以往实验中使用的X射线条纹相机易被电磁环境干扰以及相机电极部件易被实验产生的高速粒子损伤的问题.诊断系统的光谱响应范围主要集中在0.2—10keV,系统的空间分辨率经过理论评估小于120μm,通过 ...
    本站小编 Free考研考试 2021-12-29
  • 六方氮化硼单层中一种(C<sub>N</sub>)<sub>3</sub>V<sub>B</sub>缺陷的第一性原理计算
    摘要:二维六方氮化硼(hBN)的点缺陷最近被发现可以实现室温下的单光子发射,而成为近年的研究热点.尽管其具有重要的基础和应用研究意义,hBN中发光缺陷的原子结构起源仍然存在争议.本文采用基于密度泛函理论的第一性原理计算,研究hBN单层中一种B空位附近3个N原子被C替代的缺陷(CN)3VB.在hBN的 ...
    本站小编 Free考研考试 2021-12-29
  • 基于旋转不变技术信号参数估计的激光扫频干涉测量方法
    摘要:激光扫频干涉测量技术具有无测距盲区、非接触、单次测量多目标的能力.通过傅里叶变换可提取目标拍频频率,进而解算距离.然而受激光器调频带宽限制,通过傅里叶变换得到的目标分辨率受限于固有分辨率.为解决该问题,本文提出采用基于旋转不变技术的信号参数估计(ESPRIT)算法对测量信号进行频谱分析.实验通 ...
    本站小编 Free考研考试 2021-12-29
  • In掺杂<i>h</i>-LuFeO<sub>3</sub>光吸收及极化性能的第一性原理计算
    摘要:h-LuFeO3是一种窄带隙铁电半导体材料,已被证明在铁电光伏领域有较好的应用前景.然而,较低的极化强度使光生电子-空穴对复合率大,限制了h-LuFeO3基铁电光伏电池效率的提高.为改善h-LuFeO3的极化强度,提高光吸收性质,本文利用第一性原理计算方法研究了In原子在h-LuFeO3不同位 ...
    本站小编 Free考研考试 2021-12-29
  • 吉瓦级强流相对论多注电子束二极管的优化设计与实验研究
    摘要:多注相对论速调管放大器可在较高的工作频段实现GW级功率微波产生,在很多领域得到了发展和应用.多注相对论速调管中强流相对论多注电子束相互之间存在空间电磁场的作用,使得多注电子束从二极管引入多注漂移管,以及在多注漂移管中的传输运动受到影响,导致电子束会轰击到管壁上,早期实验中多注电子束的传输通过率 ...
    本站小编 Free考研考试 2021-12-29