留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

D-T中子源旋转靶活化计算与分析

高晖 沈姚崧 袁永刚

王可, 段艳涛, 石立华, 等. 一种基于双环差分结构的脉冲磁场传感器[J]. 强激光与粒子束, 2022, 34: 043003. doi: 10.11884/HPLPB202234.210337
引用本文: 高晖, 沈姚崧, 袁永刚. D-T中子源旋转靶活化计算与分析[J]. 强激光与粒子束, 2018, 30: 096002. doi: 10.11884/HPLPB201830.180039
Wang Ke, Duan Yantao, Shi Lihua, et al. A pulsed magnetic field sensor based on dual-loop differential structure[J]. High Power Laser and Particle Beams, 2022, 34: 043003. doi: 10.11884/HPLPB202234.210337
Citation: Gao Hui, Shen Yaosong, Yuan Yonggang. Activation calculation and analysis for rotary target chamber irradiated by D-T neutron[J]. High Power Laser and Particle Beams, 2018, 30: 096002. doi: 10.11884/HPLPB201830.180039

D-T中子源旋转靶活化计算与分析

doi: 10.11884/HPLPB201830.180039
详细信息
    作者简介:

    高晖(1972—), 女,副研究员,主要从事辐射防护研究; qianliyan2005@caep.cn

  • 中图分类号: TL64

Activation calculation and analysis for rotary target chamber irradiated by D-T neutron

  • 摘要: 为实现14 MeV D-T中子源旋转靶活化特性计算与分析,首次采用一种新的技术途径,将自主研发的活化程序BURNDOT与蒙特卡罗自动建模系统MCAM、蒙特卡罗粒子输运程序MCNP相结合,通过中子输运、材料活化、光子输运模拟计算的耦合,考察了材质、栅元、主要活化核素对靶室活化特性的影响。结果表明,约以辐照后68 h为界,材料铜、316不锈钢先后作为旋转靶室活度的主要贡献者,前者其产生的缓发γ剂量因62Cu, 64Cu核素的存在而达到活化剂量最大值,后者因有长半衰期核素55Fe, 57Co, 54Mn等的存在,但其产生的剂量率值低于安全限值10 μSv/h。采用新方法的计算结果与采用欧洲活化程序FISPACT-2007的计算结果符合较好。
  • 脉冲磁场的测量是电磁脉冲研究领域的重点之一[1-6]。目前,脉冲磁场传感器的发展迅速,其中,屏蔽小环磁场传感器因其结构简单、性能可靠,特别是对于中低频信号的测量精度较高而备受关注。但是磁偶极子天线的输出端结构破坏了天线的对称性由于磁偶极子天线不可避免地存在输出端结构,破坏了天线的对称性,因此对于多极性的电磁脉冲环境测量通常很难满足测量要求,极易被径向电场干扰导致磁场波形失真[7]。尤其是对雷电回击、高空核爆的近场磁场测量,现有的脉冲磁场传感器很难避免电场干扰[8]。因此,本文研制了一种基于双环差分结构的脉冲磁场传感器,通过构造差分结构可以去除绝大部分径向电场响应分量,从而获得了较为纯净的磁场响应,实现了对雷电近场环境的准确测量。

    在电气小条件下,磁偶极子天线可由激励源和等效的集总电感表示,因此单环磁场传感器的等效电路如图1所示[9]。其中L为磁偶极子的等效集总电感,ZL为负载阻抗,ω为测得的磁场信号的频率,Bi为入射磁场的磁感应强度。小环天线将被拾取磁场强度部分作为激励源,则磁场响应的大小和方向取决于入射磁场和传感器的灵敏度。此时传感器的灵敏度可以用等效面积Ae表示。激励源具有电压量纲,其拾取量可以表示为BitAe

    图  1  磁偶极子天线等效电路
    Figure  1.  Equivalent circuit of single loop magnetic field sensor

    则负载电压VL可以表示为

    V(jω)=jωAeZLjωL+ZLB(jω)
    (1)

    由式(1)可知,根据磁偶极子天线的工作频率范围,磁场传感器可以分为微分模式,过渡模式和自积分模式[10]。当ZLωL,单环传感器的负载电压VL与微分磁场成正比,可表示为

    VL=AetB
    (2)

    感应近场条件下,交变电场和交变磁场在相位上正交,因而在磁场测量中极易受到电场干扰。在屏蔽小环天线的基础上引入双环差分结构,即将两个相同的屏蔽小环天线平行镜像对称放置,如图2所示。同时将增大环天线的匹配阻抗从而减小环感应电流,以降低环间互扰带来的不利影响。其中,EH分别表示入射的脉冲电场和脉冲磁场分量,用下标1和2对双环进行区分,并将其输出端电压U按照贡献来源不同分为UEUH,分别表示电场和磁场响应的贡献。

    图  2  双环差分原理示意图
    Figure  2.  Schematic diagram of dual-loop differential principle

    根据电磁感应原理,则有

    {U1=UE1+UH1U2=UE2+UH2
    (3)

    两个开路截断的屏蔽小环天线镜像对称放置时,由于轴向磁通量反向而径向电通量相同,则双环天线的磁场响应电压分量应反相,而电场响应分量方向相同。则有

    U1U2=UH1UH2=2UH
    (4)

    联立式(2),(4)可知,差分后的输出电压即为的微分磁场与天线等效面积乘积的2倍,且并不包含电场响应分量。

    基于双环差分测量原理,本文提出一种脉冲磁场传感器如图3所示,图3(a)为工作原理图,图3(b)为实物图。传感器主要由双环天线、光传输系统组成,光传输系统又由光发射机、光纤和光接收机组成。双环天线的输出端分别与光发射机的输入端连接,光发射机的输出端与光接收机的输入端通过光纤连接,光接收机的输出端接入示波器即可实现磁场测量。

    图  3  双环式差分磁场传感器
    Figure  3.  Dual-loop differential magnetic field sensor

    由式(2)可知,环天线的输出电压与磁场微分成正比。为了得到可直接观测的磁场波形,每一个环天线的输出端串联了一个RC积分电路,在光发射机内完成电压信号的积分运算;其次,利用差分电路把两个积分磁场信号差分以去除电场响应分量,可以实现纯净磁场信号的输出。

    双环天线由两个屏蔽小环天线组成,主要负责将分布的磁场参量转换成终端电压。小环天线由一段弯成环形的半钢同轴电缆制成,左半环末端的内导体与右半环屏蔽层短路焊接。为了便于拆装,测试用的小环天线与光发射机之间选择BNC接口连接。

    双环天线应当选择两个尺寸相同的环天线,根据式(2)可知,单环传感器的量程由负载电压VL和环天线的等效面积Ae的范围共同决定。其中,负载电压VL主要由光发射机的输入电压等级决定;等效面积Ae则主要与环天线半径r有关,在进行初步估算时,可以近似认为Aer2

    同时,双环天线之间应保证镜像对称放置条件且保持适宜的环间距L。其中镜像对称是差分计算的基本前提,而环间距是影响环间互扰的主要因素之一。当L过小时,双环天线互扰加剧;当L过大时,两个环天线所处的磁场环境不同,差分后引入误差。因此,应当在互扰程度可接受的前提下取较小的环间距。

    针对选定的天线尺寸,为了匹配合适的环间距,本文利用CST Studio进行仿真计算,其仿真设置如图4所示。首先,利用一个环天线作为标准环,测得标准环在1.2/50 μs双指数标准雷电波下的响应波形,记为标准波形;其次,在相同场环境中于相距标准环距离L处增加一个镜像放置的干扰环后,测得标准环上的响应波形,记为干扰波形。

    图  4  双环天线互扰分析示意图
    Figure  4.  Mutual interference analysis of dual-loop antenna

    为了评估环间距对互扰程度的影响,本文提出一种综合失真率ε,用以表征双环天线的互扰程度大小,其计算公式给出如下

    ε = εu+εt
    (5)

    式中:εu为幅值失真率,用以表征测量波形相较于标准波形在幅值上的失真程度;εt为时间失真率,用以表征测量波形相较于标准波形在上升时间和半宽时间上的失真程度。

    εt=|T1iT1s|/2T1d+|T2iT2s|/2T2d
    (6)
    εu=|UiUs|/Ud
    (7)

    式中:UsUi分别为标准波形和干扰波形的幅值;Ud表示标准波形与干扰波形之间的幅值差的绝对值;T1sT1i分别为标准波形和干扰波形的上升时间;T1d表示标准波形与干扰波形之间的上升时间差的绝对值;T2sT2i分别为标准波形和干扰波形的半宽时间;T2d表示标准波形与干扰波形之间的半宽时间差的绝对值。

    互扰仿真时取环直径为10 cm,共设置五组不同的环间距L,分别为1,3,5,7,9 cm。仿真结果见表1。由表1可以看到,当环间距小于5 cm时,综合失真率ε大于5 %。综合考虑,当r=10 cm时,取L=5 cm。最后,根据确定的环半径和环间距,制成固定模具,用于保持双环天线相对位置固定。

    表  1  仿真结果
    Table  1.  Results of simulation
    L/cmεu/%εt/%ε/%
    110.522.5813.11
    34.632.206.83
    53.271.774.45
    70.781.512.29
    90.631.442.07
    下载: 导出CSV 
    | 显示表格

    脉冲磁场的测量通常伴随着复杂的变化的电磁环境,为了避免传输过程中的干扰,目前常用的是基于光电隔离技术的光传输解决方案[11-13]。本文采用的光传输系统主要由光发射机、单模光纤、光接收机组成。

    光发射机用于将输入的磁场微分电信号转换为输出的差分磁场光信号,在其内部依次完成输入信号的积分、差分、放大、非线性补偿、滤噪、电光转换,最终实现光信号的输出。

    光发射机内置两块8.4 V锂电池,可以实现±8.4 V直流供能。另外,光发射机采用铝制外壳,以保护信号在其内部传递不受干扰,同时严格限制孔缝结构,使其可以在强电磁场中工作。经试验验证,其屏蔽效能超过80 dB。

    单模光纤用于连接光发射机和光接收机,其长度可以根据所需布线长度而调整。通常利用光纤衰减系数来评估光纤对光信号功率的衰减值。对于单模光纤,光纤衰减系数通常取0.18 dB/km[14]

    光接收机通常置于屏蔽良好的室内或者屏蔽柜工作,其输入端通过ST接口与光纤相连,输出端可以通过BNC接口与示波器等采集设备相连。光接收机用于实现磁场光信号到磁场电信号的转换,光信号在其内部经由光电转换、放大、补偿、阻抗转换模块,最终在其输出端复现出光发射机输入级的电信号。

    本文中,我们提出了双环差分测量原理,并设计研制了一种用于近场测量的脉冲磁场传感器,这种新型的磁场传感器,相较于传统的单环传感器具有较强的抗干扰能力,特别是在强电场环境中,能够实现对磁场的准确测量。同时,双环差分结构可以适用于多种环天线,从而匹配不同磁场环境的测量要求,因此在脉冲磁场测量领域有着较为广泛的应用前景。

    依据国军标GJB1389A,在实验室搭建了一套邻近雷击环境模拟装置[15],如图5所示。模拟装置主要由Marx发生器,场照射器,截断装置和电容分压器组成。其中,Marx发生器用于产生1.2/50 μs标准雷电波冲击高压;场照射器用于承载被测物并产生相对均匀的电磁脉冲环境;截断装置用于产生截断脉冲高压,形成陡峭的下降沿;电容分压器用于接地系统的高压测量[16]

    图  5  邻近雷击环境模拟装置
    Figure  5.  Mutual interference analysis of dual-loop antenna

    由于场照射器采用平行板式,因此在正负电极两端存在巨大瞬时电位差,导致在垂直方向存在强大的瞬变电场。在利用传统单环磁场传感器测量水平磁场时,即使采取较好的屏蔽手段也极易出现测量结果被干扰的现象,因此考虑采用双环差分结构去除电场响应分量,从而实现磁场的准确测量。

    对比图6三个波形可以看到:(1) 单环测量波形在0.4~1.5 μs之间存在着强烈的干扰,与双环测量波形在波形前端并不吻合;(2) 单环测量波形前端的干扰在1.5 μs之后消失,与双环测量波形在波形后端基本吻合;(3) CST仿真波形前端与双环天线相似且并不存在振荡分量;(4) 三个波形在波形后端具有很高的一致性,基本符合被测磁场环境。

    图  6  磁场测量结果对比图
    Figure  6.  Comparison of magnetic field measurement results

    综上所述,单环测量波形在0.4~1.5 μs之间的波形畸变显然不是磁场响应而应该是电场干扰,可能是由于截断装置击穿空气间隙时产生的高能快速变化的电场导致。对比单环测量和双环测量结果可以看到,双环测量波形通过差分去除了这部分干扰分量从而得到了较为纯净的磁场响应。因此,试验和仿真结果可以证明,双环传感器相较于传统的单环传感器具有更强的抗干扰能力。

    本文提出了双环差分测量原理,并设计研制了一种用于近场测量的脉冲磁场传感器,这种新型的磁场传感器,相较于传统的单环传感器具有较强的抗干扰能力,特别是在强电场环境中,能够实现对磁场的准确测量。同时,双环差分结构可以适用于多种的环天线,从而匹配不同磁场环境的测量要求,因此在脉冲磁场测量领域有着较为广泛的应用前景。

  • 图  1  旋转靶室几何模型剖视图

    Figure  1.  Sectional view of geometrical modeling of rotary target chamber

    图  2  活化计算流程图

    Figure  2.  Flowchart of activation calculation

    图  3  中子能谱

    Figure  3.  Neutron spectra of cells

    图  4  不同栅元中子注量分布

    Figure  4.  Distribution of neutron flux of cells

    图  5  不同材料活度随时间变化曲线

    Figure  5.  Curves of activity generated from different materials vs time

    图  6  停止辐照初时(T=0)不同栅元活度贡献

    Figure  6.  Proportion of activity of cells at the first time after neutron radiation

    图  7  铜生核素的活度随时间变化曲线

    Figure  7.  Curves of nuclide activity generated from copper vs time

    图  8  316不锈钢生核素的活度随时间变化曲线

    Figure  8.  Curves of nuclide activity generated from 316 stainless steel vs time

    图  9  铝生核素的活度随时间变化曲线

    Figure  9.  Curves of nuclide activity generated from aluminum vs time

    图  10  停止辐照初时(t=0))缓发γ三维剂量场分布

    Figure  10.  Three dimensional delayed γ dose distribution right after neutron radiation

    图  11  辐照后缓发γ剂量率随距离变化曲线

    Figure  11.  Curves of delayed γ dose rate at different distance from point A after neutron radiation over 8 h

    图  12  A处缓发γ剂量率随栅元变化曲线

    Figure  12.  Curves of delayed γ dose rate from different cells at A location after neutron radiation

    表  1  旋转靶室栅元结构参数

    Table  1.   Structure parameters of cells consisting the rotary target chamber

    cell No. material volume/cm3 mass/g function
    1 316 stainless steel 401.730 3 166.20 activation
    2 316 stainless steel 90.714 714.94 activation
    3 pure copper 69.115 619.27 activation
    4 pure copper 90.085 807.16 activation
    5 316 stainless steel 11.819 93.15 activation
    6 pure aluminum 18.148 49.00 activation
    7 pure aluminum 147.800 399.06 activation
    8 pure aluminum 64.280 173.56 activation
    9 316 stainless steel 48.119 379.24 activation
    10 water 16.840 16.84 cooling
    11 pure copper 89.892 805.43 activation
    12 pure titanium 0.204 0.92 activation
    13 316 stainless steel 382.910 3 017.80 activation
    14 316 stainless steel 1 146.400 9 035.20 activation
    下载: 导出CSV
  • [1] 李德平, 潘自强. 辐射防护手册[M]. 北京: 原子能出版社, 1987.

    Li Deping, Pan Ziqiang. Radiation protection manual. Beijing: Atomic Energy Press, 1987
    [2] Shypailo R J. Stability evaluation and correction of a pulsed neutron generator prompt gamma activation analysis system[J]. J Radioanal Nucl Chem, 2016, 307: 1781-1786. doi: 10.1007/s10967-015-4440-7
    [3] Vuolo M, Bonifetto R, Dulla S, et al. Evaluation of the neutron activation of JET in-vessel components following DT irradiation[J]. Fusion Engineering and Design, 2014, 89: 2071-2075. doi: 10.1016/j.fusengdes.2014.02.050
    [4] Forrest R A. FISPACT-2007: User manual[R]. UKAEA FUS 534, 2007.
    [5] Croff A G. User's manual for Origen2 computer code[R]. ORNL/TM-7175, 1980.
    [6] Liu Xiaoping, Luo Yuetong, Tong Lili, et al. Development and application of MCNP auto-modeling tool: Mcam 3.0[J]. Fusion Engineering and Design, 2005, 75/79: 1275-1279. doi: 10.1016/j.fusengdes.2005.06.236
    [7] Briesmeister J F. MCNP―A general Monte Carlo N-particle transport code[R]. LA-370-M, 2000.
    [8] Rodenas J. Application of the Monte Carlo method to estimate doses due to neutron activation of different materials in a nuclear reactor[J]. Radiation Physics and Chemistry, 2017, 140: 442-446. doi: 10.1016/j.radphyschem.2017.02.015
    [9] Rodenas J, Gallardo S, Weirich F, et al. Application of dosimetry measurements to analyze the neutron activation of a stainless steel sample in a training nuclear reactor[J]. Radiat Phys Chem, 2014, 104: 368-371. doi: 10.1016/j.radphyschem.2014.05.013
    [10] Rodenas J, Abarca A, Gallardo S. Analysis of dose rates received around the storage pool for irradiated control rods in a BWR nuclear power plant[J]. Applied Radiation and Isotopes, 2011, 69: 1104-1107. doi: 10.1016/j.apradiso.2010.10.019
  • 加载中
图(12) / 表(1)
计量
  • 文章访问数:  960
  • HTML全文浏览量:  279
  • PDF下载量:  139
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-01-26
  • 修回日期:  2018-05-18
  • 刊出日期:  2018-09-15

目录

/

返回文章
返回