Monte Carlo simulation of Cherenkov light generated by underwater Co-60 source
-
摘要: 随着核应用领域的不断拓宽,放射源丢失事故发生的概率也随之增加。机载伽马谱仪可有效搜寻地面放射源,然而对于放射源丢失于水域的情况,由于伽马射线经由水层屏蔽后可探测性降低,故利用放射源在水中产生的切伦科夫辐射对其进行搜寻显得十分重要。采用MCNP与Geant4相结合的方法,以及在Geant4程序中采用接续计算技巧,对Co-60源在水中的切伦科夫光产生以及传输进行了计算,计算表明,切伦科夫光经水中传播后,主要波段在300~600 nm,强度呈由边缘到中心渐强的特征分布,分布范围大致与放射源在水中的深度一致,在水中传输300 m后其光通量约为100 cm-2,可利用光谱特征和强度分布特征对其进行测量。
-
关键词:
- 放射源丢失 /
- 切伦科夫辐射 /
- MCNP耦合Geant4计算 /
- Geant4接续计算 /
- 光学测量
Abstract: With wide application of nuclear techniques, accidents of lost radioactive sources increase. The airborne gamma spectrometer can be used for searching the lost radiation sources on the ground level. However, for radioactive sources lost in water, the use of gamma spectrometer is limited as a result of the shielding of gamma rays by water. So detection of underwater radioactive source based on Cherenkov light generated by the radioactive source is becoming important. With applications of combined simulation of Geant4 and MCNP, and continuation simulation method in Geant4, distributions and transmission of Cherenkov light generated by underwater Co-60 sealed source were simulated. The simulation reveals that wavelength of Cherenkov light is between 300~600 nm through transmission in water. The light intensity becomes stronger from the edge to the center, and the distribution range approximately equals to the depth of the radioactive source in water. The light flux is about 100 Cherenkov photons·cm-2 after 300 m transmission in water. The Cherenkov light can be detected by the characteristics of its wavelength spectrum and intensity distribution. -
为满足物理实验需要,ICF激光驱动器装置的终端光学组件工作在以3ω为主的三波长激光同时辐照环境下,且3ω通量往往都接近元件的损伤阈值,终端光学组件是装置稳定、高效和安全运行的薄弱环节[1-5]。由于同时存在波长效应和波长间的能量耦合效应,多波长同时辐照下的初始损伤和损伤增长的物理机制比单波长要复杂很多[6-9]。研究多波长激光同时辐照下的初始损伤和损伤增长的实验现象,分析多波长激光同时辐照下损伤过程的能量耦合机制,有助于我们理解激光诱导光学元件损伤的物理机制,预测光学元件的使用寿命,对系统进行优化设计,提高系统的负载能力并降低运行成本[10-14]。本文研究了3ω单独辐照、3ω+2ω和3ω+1ω双波长同时辐照下熔石英的初始损伤和损伤增长实验规律,获得了在3ω阈值附近,低能量密度的2ω和1ω对初始损伤和损伤增长的影响规律,利用基于飞秒双脉冲的阴影成像方法获得了冲击波速度,分析了多波长激光同时辐照下光学元件损伤的能量耦合机制。
1. 实验系统
本文所采用损伤测试实验系统及实验方法见文献[15]。为了使研究结果对工程应用具有参考意义,选择各波长能量密度比与工程装置终端组件相近的情况作为实验条件。损伤坑面积采用辐照后离线显微成像的方式获得。采用图1所示的泵浦−探针阴影成像系统进行冲击波速度测量,Nd:YAG纳秒多波长激光作为泵浦脉冲,钛宝石飞秒激光作为探针脉冲。飞秒脉冲经1∶1分光镜分成两束,其中一束经延迟线产生延迟时间Δt。利用飞秒激光器外同步触发信号发生器(DG645),信号发生器控制纳秒激光与飞秒激光的延时ΔT。控制CCD曝光时间,采用单发双脉冲成像方法,单次成像可获得同一冲击波的两个波前,经长度定标可得到Δt内冲击波前的平均速度。利用光电探头和示波器获取ΔT,由于飞秒种子和纳秒激光器输出的时间抖动,ΔT存在±10 ns的抖动。实验样品厚度2 mm,透镜聚焦250 mm。
实验样品基材为Corning 7980,尺寸为100 mm×70 mm×5 mm和100 mm×10 mm×2 mm,表面精抛光,实验前用氢氟酸(HF)和氟化氨(NH4F )溶液在超声波环境下进行表面刻蚀,经去离子水清洗,并经高纯酒精脱水处理后擦净。实验样品在1ω,2ω和3ω单独辐照下的最大零几率损伤阈值分别约为64,45和19 J/cm2,损伤增长阈值分别约为15,10和6 J/cm2[15]。
2. 实验研究
2.1 双波长同时辐照下的初始损伤几率
表1为不同能量密度组合的3ω+2ω,3ω+1ω同时辐照下熔石英的初始损伤几率,虽然从3ω@17 J/cm2+2ω@9 J/cm2,3ω@17 J/cm2+1ω@10 J/cm2能量密度组合已经出现损伤,但随着能量密度逐渐增加,损伤几率没有明显变化,考虑样品差异和损伤的随机性、离散性,可以认为2ω和1ω对损伤几率没有明显贡献。在3ω@24 J/cm2+2ω@9 J/cm2,3ω@24 J/cm2+1ω@20 J/cm2时损伤几率分别增加到10%和14%,此时,可认为2ω和1ω对损伤几率有贡献。
表 1 双波长同时辐照下熔石英的初始损伤几率Table 1. Damage probability of fused silica irradiated by dual wavelength laserwavelengths energy density/(J·cm−2) damage probability/% wavelengths energy density/(J·cm−2) damage probability/% 3ω+2ω 17+3 0 3ω+1ω 17+5 0 17+6 0 17+10 3 17+9 2 17+15 4.5 17+12 5.5 17+20 6 20.5+3 7 20.5+5 4.5 20.5+6 5 20.5+10 4 20.5+9 1.5 20.5+15 3 20.5+12 3.5 20.5+20 8.5 24+3 3 24+5 5 24+6 8 24+10 7 24+9 10 24+15 8 24+12 12 24+20 14 27.5+3 31 - - 2.2 双波长同时辐照下的初始损伤面积
表2为不同能量密度(组合)的3ω单独、3ω+2ω和3ω+1ω双波长同时辐照下的损伤几率和损伤坑平均面积的实验结果。3ω单独辐照下,随着3ω能量密度增加,损伤几率和损伤坑平均面积都快速增加。在双波长同时辐照下,随着2ω和1ω的能量密度增加,虽然损伤几率没有明显变化,但损伤坑平均面积呈单调增加趋势。
表 2 2ω和1ω对损伤几率和损伤坑平均面积的影响Table 2. Influence of 2ω and 1ω wavelengths on damage probability and damage degreewavelengths energy density/(J·cm−2) damage probability/% average area/μm2 3ω 20.5 1.5 2 043 24 4 3 180 30.5 28 4 565 31 33.5 4 737 33 47 6 455 3ω+2ω 20.5+3 4 2 105 20.5+6 5 2 517 20.5+9 1.5 4 029 20.5+12 3.5 5 134 3ω+1ω 20.5+5 4.5 1 903 20.5+10 4 2 480 20.5+15 3 2 507 20.5+20 8.5 3 701 2.3 双波长同时辐照下的损伤增长
表3和表4是3ω+2ω,3ω+1ω双波长同时辐照下熔石英的损伤增长实验规律。可以看出,在3ω@5 J/cm2+2ω@2.4 J/cm2和3ω@6 J/cm2+1ω@2.4 J/cm2的情况下,2ω和1ω对损伤增长开始有明显贡献。
表 3 3ω+2ω同时辐照下熔石英的损伤增长几率Table 3. Damage growth probability of fused silica irradiated by 3ω and 2ω simultaneouslyNo. 3ω energy density/(J·cm−2) 2ω energy density/(J·cm−2) damage growth probability/% 1 4.4 2.4 0 2 5 0 0 3 5 1.2 0 4 5 2.4 5 5 6 0 0 6 6 1.2 11 7 6 2.4 17 8 6.5 0 13 9 6.5 1.2 18 10 6.5 2.4 26 表 4 3ω+1ω同时辐照下熔石英的损伤增长几率Table 4. Damage growth probability of fused silica irradiated by 3ω and 1ω simultaneouslyNo. 3ω energy density/(J·cm−2) 1ω energy density/(J·cm−2) damage growth probability/% 1 4.4 2.4 0 2 5 0 0 3 5 1.2 0 4 5 2.4 0 5 6 0 0 6 6 1.2 0 7 6 2.4 7 8 6 3.6 12 9 6.5 0 10 10 6.5 1.2 14 11 6.5 2.4 20 2.4 双波长同时辐照下的超快成像
图2是脉冲能量均为30 mJ的1ω,3ω单独辐照和同时辐照下,熔石英后表面损伤的典型阴影成像图,Δt=28 ns。其中1,2和3分别为体内一次冲击波前、体内二次冲击波前和空气冲击波前,1′,2′和3′为Δt时间间隔后的冲击波前。选择体内一次冲击波为研究对象,图2(a),(b)和(c)对应的冲击波前在Δt时间间隔内的平均速度分别为4.58×103,6.07×103和6.42×103 m/s。
3. 分析与讨论
双波长同时辐照下,在3ω单独辐照损伤阈值和损伤增长阈值附近,当能量密度比2ω∶3ω≥0.2∶1(损伤增长)和≥0.5∶1(初始损伤)时,2ω会对损伤增长和初始损伤产生明显贡献;当能量密度比1ω∶3ω≥0.4∶1(损伤增长)和≥0.8∶1(初始损伤)时,1ω会对损伤增长和初始损伤产生明显贡献;在初始损伤中,能量密度比2ω∶3ω<0.5∶1的2ω和1ω∶3ω<0.8∶1的1ω虽然对损伤几率影响较小甚至没有影响,但会引起损伤坑平均面积增加。
对于多波长同时辐照下熔石英的初始损伤和损伤增长实验现象一个可能合理的解释:
(1)初始损伤及损伤增长过程均可分为两个阶段。第一阶段,缺陷(损伤坑)吸收激光产生自由电子并形成低温等离子体(其膨胀不足以产生损伤(增长));第二阶段,低温等离子体高效吸收后续激光,内能增加,形成高温等离子体,其剧烈膨胀产生损伤(增长)。
(2)双波长同时辐照下,当3ω能量密度在其单独辐照的阈值附近,能量密度远低于其自身损伤阈值和损伤增长阈值的2ω和1ω被3ω产生的自由电子、低温等离子体和高温等离子体吸收,从而引起损伤几率(程度)增加。即,如果2ω和1ω被3ω产生的自由电子或低温等离子体大量吸收,形成高温等离子体,则主要引起损伤几率增加;如果2ω和1ω主要被3ω产生的高温等离子体吸收,导致等离子体内能进一步增加,则主要引起损伤程度增加。从这个角度来看,多波长同时辐照下,只要其中一个波长对元件产生了损伤或损伤增长现象,其产生的自由电子和等离子体就会对其他波长激光产生吸收,从而对初始损伤和损伤增长造成影响。
根据Taylor等人[16]提出的等离子体点爆炸球面波模型
v∝γ(Eρ)1/5t2/5 (1) 式中:v是冲击波前速度;E是等离子体能量;γ为无量纲常数;ρ是介质密度。根据2.4的实验测试结果,
6.425>6.075+4.585 ,即双波长同时辐照下等离子体积累的能量大于两个波长单独辐照下等离子体积累的能量之和,表明波长间的耦合效应促进了激光能量向材料(等离子体)沉积的效率。当然,这需要忽略缺陷、激光及相互作用过程的差异。这种现象可以解释为,较高能量密度的3ω被缺陷吸收产生自由电子和等离子体,自由电子和等离子体吸收1ω,使得原本对1ω不吸收(或吸收很小)的缺陷产生了吸收;同时,吸收1ω后,电子密度和等离子体温度增加,又会增加其对3ω脉冲后续能量的吸收效率;使得原本单独吸收3ω不会损伤的产生了损伤,原本会产生损伤的增加了损伤程度。4. 结 论
多波长激光同时辐照下熔石英元件的初始损伤和损伤增长过程中,同时存在激光与材料相互作用的波长效应和波长间的能量耦合效应,自由电子和等离子体是波长间能量耦合效应的吸收媒介,会促进激光能量向材料沉积的效率。当3ω能量密度在其损伤阈值和损伤增长阈值附近,能量密度远低于其自身阈值的2ω和1ω会被3ω产生的自由电子和等离子体吸收,虽然可能不会对损伤几率和损伤增长阈值产生影响,但会增加初始损伤程度和损伤增长系数。
-
表 1 Co-60外壳表面伽马、电子源强
Table 1. Source intensity of gamma radiation and electrons
position gamma radiation
/(particles·s-1)electrons
/(particles·s-1)tube surface 1.358 4×1014 1.116 1×1012 top surface 6.822 4×1012 5.741 7×1010 bottom surface 2.126 1×1013 1.748 1×1011 -
[1] 闻良生, 龚频, 黄茜, 等. 小型旋翼机机载辐射环境监测系统的设计与实现[J]. 强激光与粒子束, 2016, 28: 106004. doi: 10.11884/HPLPB201628.160036Wen Liangsheng, Gong Pin, Huang Xi, et al. Design and implementation of minitype rotorcraft airborne radiation monitoring system. High Power Laser and Particle Beams, 2016, 28: 106004 doi: 10.11884/HPLPB201628.160036 [2] 倪卫冲, 刘士凯, 高国林, 等. AGS-863航空伽马能谱勘查系统机载试验[J]. 中国核科学技术进展报告, 2011, 2(1): 335-343. https://cpfd.cnki.com.cn/Article/CPFDTOTAL-EGVD201110001060.htmNi Weichong, Liu Shikai, Gao Guolin, et al. Airborne testing of AGS-863 airborne gamma spectrometry survey system. Progress Report on China Nuclear Science & Technology, 2011, 2(1): 335-343 https://cpfd.cnki.com.cn/Article/CPFDTOTAL-EGVD201110001060.htm [3] 翁渝民. 单光子计数-弱信号检测的有力手段[J]. 物理, 1980, 9(1): 20-24. https://www.cnki.com.cn/Article/CJFDTOTAL-WLZZ198001007.htmWeng Yumin. Single photon counting-efficient technique for weak single measurement. Physics, 1980, 9(1): 20-24 https://www.cnki.com.cn/Article/CJFDTOTAL-WLZZ198001007.htm [4] 舒迪昀, 汤晓兵, 侯笑笑, 等. 基于Cerenkov效应水下放射源搜寻技术的可行性分析研究[J]. 原子能科学技术, 2015, 49(4): 582-588. https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201504002.htmShu Diyun, Tang Xiaobing, Hou Xiaoxiao, et al. Analysis of feasibility for searching underwater radioactive source using Cerenkov effect. Atomic Energy Science and Technology, 2015, 49(4): 582-588 https://www.cnki.com.cn/Article/CJFDTOTAL-YZJS201504002.htm [5] 刘斌, 贾清刚, 张天奎, 等. 水下切伦科夫光光斑的蒙特卡罗模拟[J]. 强激光与粒子束, 2013, 25(1): 196-200. doi: 10.3788/HPLPB20132501.0196Liu Bin, Jia Qinggang, Zhang Tiankui, et al. Monte Carlo simulation of Cherenkov light spot produced by underwater radioactive source. High Power Laser and Particle Beams, 2013, 25(1): 196-200 doi: 10.3788/HPLPB20132501.0196 [6] Agostinelli S, Allison J, Amako K, et al. Geant4—a simulation toolkit[J]. Nuclear Instruments and Methods in Physics Research A, 2007, 506(3): 250-303. https://www.sciencedirect.com/science/article/pii/S0168900203013688 [7] Allison J, Amako K, Apostolaki J, et al. Geant4 developments and applications[J]. IEEE Trans Nuclear Science, 2006, 53(1): 270-278. https://ieeexplore.ieee.org/document/1610988/ [8] Zhang Qingmin, Hu Zhigang, Deng Bangjie, et al. A simple iterative method for compensating response delay of self-powered neutron detector[J]. Nuclear Science and Engineering, 2017, 186(1): 293-302. [9] GB7465-2009. 高活度钴60密封放射源[S]. 中华人民共和国国家标准, 2009.GB7465-2009. High activity cobalt-60 sealed radioactive sources. PRC standard, 2009 [10] Pope R M, Fry E S. Absorption spectrum(380~700 nm) of pure water. Integrating cavity measurement[J]. Appl Opt, 1997, 36(33): 8710-8723. https://pubmed.ncbi.nlm.nih.gov/18264420/ [11] Quickenden T I, Irvin J A. The ultraviolet absorption spectrum of liquid water[J]. J Chem Phys, 1980, 72(8);4416-4428. [12] 曹婷婷, 罗时荣. 天空直射光谱和天空光谱的测量与分析[J]. 物理学报, 2006, 56(9): 5554-5557. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB200709095.htmCao Tingting, Luo Shirong. Measurement and analysis of direct sunlight and skylight spectra. Acta Phsica Sinica, 2006, 56(9): 5554-5557 https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB200709095.htm [13] 徐英莹, 金伟其. 夜晚天空光谱辐射测量研究及光谱去噪分析[J]. 光谱学与光谱分析, 2012, 32(6): 1456-1459. https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201206006.htmXu Yingying, Jin Weiqi. Measurement of night sky spectral radiation and analysis of spectral denoising. Spectroscopy and Spectral Analysis, 2012, 32(6): 1456-1459 https://www.cnki.com.cn/Article/CJFDTOTAL-GUAN201206006.htm 期刊类型引用(4)
1. 丁文豪,魏望和,朱焜,陶聪,冷劲松. W波段半矩形环螺旋线行波管设计. 机电工程技术. 2025(02): 60-63+82 . 百度学术
2. 詹晓非,朱增伟. 太赫兹矩形折叠波导慢波结构铜叠片电铸厚度均匀性研究. 真空电子技术. 2022(04): 67-72 . 百度学术
3. 王自成,唐伯俊,谢文球,田宏,董芳. 0.22THz高效率行波管的互作用计算. 强激光与粒子束. 2016(02): 88-93 . 本站查看
4. 谢辅强,丁桂甫,赵小林,程萍. 140GHz带状注弯折波导SU-8工艺与传输特性(英文). 强激光与粒子束. 2015(08): 148-154 . 本站查看
其他类型引用(6)
-