烧蚀对强脉冲离子束在高分子材料中能量沉积的影响*

张世健 喻晓 钟昊玟 梁国营 许莫非 张楠 任建慧 匡仕成 颜莎 Gennady Efimovich Remnev 乐小云†

1) (北京航空航天大学物理学院， 北京 100191)2) (北京航空航天大学， 大数据精准医疗高精尖创新中心， 北京 100191)3) (北京航空航天大学， 先进核能材料与物理北京市重点实验室， 北京 100191)4) (托木斯克理工大学， 托木斯克 634050， 俄罗斯)5) (北京大学重离子物理研究所， 北京 100871)(2020 年2 月12日收到; 2020 年3 月23日收到修改稿)

高能量密度纳秒量级强脉冲离子束辐照材料表面的烧蚀产物和束流的相互作用， 可能对束流在靶中的能量沉积产生影响， 进而影响烧蚀情况下的束流分析和相关应用的优化. 本文采用红外成像方法对横截面能量密度1.5—1.8 J/cm2的强脉冲离子束在304不锈钢和高分子材料上的能量沉积进行了测量分析. 结果表明在高分子材料上， 在超过一定能量密度后， 束流引发材料表面烧蚀产物的屏蔽效应使得大部分束流能量不能沉积在靶上. 采用有限元方法对束流引发的温度场分布进行了计算， 验证了高分子材料的低热导率以及低分解温度使其在脉冲辐照早期即开始热解， 烧蚀产物对后续束流能量的进一步沉积产生屏蔽. 此类效应在金属上存在的可能性和对束流诊断等应用的影响， 亦进行了讨论.

1 引 言

强脉冲离子束(intense pulsed ion beam，IPIB)源于20世纪60年代对惯性约束核聚变(inertial confinement fusion， ICF)点火技术的研究[1]. 当IPIB作用于材料表面时， 能够在材料表面微米尺度深度内形成极高的功率密度沉积， 导致材料表面温度剧烈地上升和下降， 伴随产生快速熔化、汽化、重凝. 在这个过程中， 材料表面的硬度和韧性等参数能够得到显著的提高[2，3]. IPIB的这种特性使其在材料处理及表面改性方面获得应用并展现出良好的发展前景[4，5].

当IPIB能量密度较高时， 靶表面在剧烈的辐照热效应下会产生烧蚀等离子体并向外扩散[6，7].利用IPIB产生的烧蚀等离子体可以以较高的沉积速率进行薄膜制备[8，9]. 在一定反应气体条件下， 利用IPIB产生的烧蚀等离子体可以制备纳米粉末[10，11].在束流辐照过程中， 烧蚀产物可能与束流发生相互作用从而对束流的能量沉积产生影响. 尤其是束流中的离子可能被烧蚀产物阻止， 使其能量耗散在烧蚀产物中而不能充分沉积在靶上， 即烧蚀产物可以对离子在靶上的能量沉积产生屏蔽. 这对较高能量密度下束流的诊断和在烧蚀条件下束流辐照效应的分析， 都会产生影响. 以往的研究对IPIB烧蚀效应的探索主要集中在对烧蚀质量损失和烧蚀产物的研究上[12-14]， 在这些研究中， 辐照材料主要为金属， 烧蚀效应相对较弱， 烧蚀产物密度较低， 对束流的屏蔽效应并不明显. 而在采用IPIB对高分子材料进行的改性研究中， 由于靶材导热率和分解温度均较低， 在较弱的束流辐照下靶材表面即可能产生烧蚀， 并和束流相互作用对其能量沉积产生影响. 故在有较为稠密烧蚀产物的情况下， 研究辐照过程中烧蚀产物和束的相互作用对于认识辐照参数、理解辐照机制具有重要的意义.

本文选用具有较低的热导率以及较低的分解温度的高分子材料聚碳酸酯(polycarbonate， PC)、聚氯乙烯(Polyvinyl chloride， PVC)、聚甲基丙烯酸甲酯(polymethyl methacrylate， PMMA)作为靶材， 选用304不锈钢作为对比材料研究了较为稠密的烧蚀产物对IPIB在靶上能量沉积的影响. 为了对于烧蚀的程度进行预测， 采用蒙特卡罗和有限元方法对辐照产生的温度场分布进行了计算， 并结合实验数据进行了分析.

2 实 验

IPIB辐照实验在北京航空航天大学物理学院BIPPAB-450强脉冲粒子加速器上进行. 该加速器加速电压最高达450 kV， 最大束流密度150 A/cm2，束流横截面能量密度达1.8 J/cm2， 脉冲长度(半高宽)约80 ns. IPIB通过磁绝缘二极管产生， 离子成分为70%的H+和30%的Cn+， 为了提高束流密度， 采用15 cm的圆锥型铜束流器辅助束流聚焦.

烧蚀实验选用尺寸分别为150 mm × 150 mm ×0.125 mm， 150 mm × 150 mm × 0.15 mm， 150 mm ×150 mm × 0.25 mm和150 mm × 150 mm × 0.1 mm的PC， PVC， PMMA和304不锈钢作为靶材. 靶材垂直于束流的传输方向放置于聚束铜罩的出口.IPIB在靶材上形成的热斑通过机械臂控制的FLUKE Ti25红外相机在脉冲发射后0.1 s内获得， 并由此计算束流沉积于靶上的横截面能量密度分布[15].

3 模 拟

为了计算IPIB辐照材料产生的热场分布， 采用傅里叶传热方程描述IPIB辐照材料引发的热场分布:

其中 ρ (T)，C (T)，λ (T) 分别 代 表材料的 密 度， 比热容和热导率; P 为热源， 即IPIB在材料中产生的功率密度分布[16].

其中 U (x，y) 为离子束横截面能量密度分布， 可通过红外成像或量热器获得; d(z)为将离子阻止本领dE/dz(z) 深度归一化的离子能损分布函数， 表征IPIB能量沉积沿深度的概率密度分布; g(t)为时间归一化的IPIB功率演化函数， 表征IPIB能量沉积随时间的概率密度分布.

初始条件为

其中T0为靶的初始温度， 取环境温度298 K.

对于边界条件， 取斯特藩-玻尔兹曼边界条件

其中j为辐射功率密度通量， ε为材料辐射系数，σ为斯特藩-玻尔兹曼常数， T为靶表面的温度，T0为环境温度.

4 结果与讨论

离子束二极管采用几何聚焦， 同时采用锥形聚束器对束流横向分布进行约束以提高束流密度[17].

图1为IPIB辐照前后304不锈钢和高分子材料背面温度分布图. 束流在不锈钢上产生的热斑近似圆形， 而且束流辐照区域和未辐照区域存在较为清晰的边界， 如图1(b)所示. 根据不锈钢上的热斑， 束流中心区域能量密度最高， 但对于高分子材料， 如图1(c)—图1(e)所示， 束斑中心区域存在较为明显的低温区域， 即在束流能量较高的位置， 沉积的能量密度反而较低.

图 1 IPIB辐照前后靶背面温度分布图 (a) 辐照前304不锈钢; (b) 辐照后304不锈钢; (c) 辐照后PC; (d) 辐照后PVC; (e) 辐照后PMMAFig. 1. Distribution of temperature on rear face before and after IPIB irradiation: (a) 304 stainless steel， before irradiation;(b) 304 stainless steel， after irradiation; (c) PC， after irradiation; (d) PVC， after irradiation; (e) PMMA， after irradiation.

图2 为IPIB在304不锈钢和高分子材料上束流中心部位横截面能量密度分布. 根据温度变化的范围可知， IPIB在不锈钢和高分子材料上产生的能量沉积区域是近似的， 但是当能量密度超过一定阈值之后， 在高分子材料上， 沉积能量密度会产生明显的下降， 且随着束流能量密度的提高， 沉积能量的下降更为显著. 如图2(b)所示， 当能量密度大于0.28 J/cm2时， IPIB沉积在PC上的能量密度开始明显降低， 直到束流焦点位置， 能量沉积降低至0.16 J/cm2. 对于PVC和PMMC， 束流密度在不超过0.2 J/cm2时会引发类似效应， 如图2(c)—图2(d)所示.

图 2 沿x方向能量密度分布图 (a) 304不锈钢; (b) PC; (c) PVC; (d) PMMAFig. 2. Distribution of energy density along x direction: (a) 304 stainless steel; (b) PC; (c) PVC; (d) PMMA.

图 3 能量密度为1 J/cm2的IPIB产生的功率密度 (a) 304不锈钢; (b) PCFig. 3. IPIB power density distribution with cross-sectional energy density 1 J/cm2 in (a) 304 stainless steel; (b) PC.

为了更准确地分析出现这种现象的原因， 使用有限元方法对IPIB辐照材料产生的热场分布进行分析. 图3为IPIB辐照304不锈钢和PC产生的功率密度， 由于离子射程的差异， 能量密度为1 J/cm2的IPIB辐照与304不锈钢和PC表面产生的最大功率密度分别为1.7 × 1017和6.6 ×1016W/m3. 以现有参数对PC进行热场模拟， 如图4(a)所示， 在能量密度为1 J/cm2的IPIB辐照下， 其表面数微米深度内在脉冲前期约100 ns时已达到热解温度(约580 K)[18-20]， 这意味着在束流能量沉积的初期， 材料表面区域会发生剧烈的烧蚀， 表面数微米深度范围内都会由于烧蚀而脱离材料表面， 使得烧蚀产物足够稠密， 可以屏蔽束流中后续的离子， 从而对束流的能量沉积产生显著影响. 对于更高能量密度的IPIB， 材料发生烧蚀的时间更早， 程度也更强烈， 产生的对能量的屏蔽效应也更明显， 因此在能量密度更高的辐照中心区域沉积在材料上的能量反而更少， 如图2(b)—图2(d)所示. 由于PC的热导率很低， 束流辐照会产生急遽的温度上升， 且在表面达到最高温度之后， 温度降低的速率很小， 会有利于烧蚀的产生. 对于该能量密度的IPIB产生的温度场， 如图4(b)所示， 在304不锈钢上最高的温度约为1800 K， 由于该温度显著低于沸点， 可以认为没有剧烈的烧蚀产生，不会发生能量屏蔽现象， 如图2(a)所示. 由于304不锈钢的热导率较高， 在达到温度峰值之后，热场向靶内部扩展， 并使近表面区域温度降低的速率显著高于PC， 进一步降低了烧蚀发生的可能性.

图 4 1 J/cm2的IPIB作用下的热场变化 (a) 304不锈钢; (b) PCFig. 4. Thermal field distribution after irradiation of IPIB with cross-sectional energy density of 1 J/cm2 in (a) 304 stainless steel; (b) PC.

5 结 论

为了研究烧蚀过程对IPIB能量沉积的影响，使用高分子材料和304不锈钢作为靶材， 通过红外成像系统测量IPIB在靶上的能量密度沉积分布，对存在烧蚀情况下束流的能量沉积进行了实验和理论上的分析. 研究表明， 高分子材料的低热导率使其在接受IPIB辐照时在近表面区域快速升温，较低的分解温度使其在IPIB脉冲的初级阶段发生较为剧烈的烧蚀， 烧蚀产物会在高分子材料的表面形成一个屏障对束流中的后续离子产生阻止作用.当束流能量密度提高时， 由于烧蚀更加剧烈， 束流能量的屏蔽效应更加显著. 在束流能量密度较高时， 在金属等材料表面也有可能产生类似的效应.这种烧蚀屏蔽效应的存在， 使得束流能量较强， 足以在辐照初期形成稠密烧蚀产物影响能量沉积的情况下的束流能量沉积的分析存在困难， 基于量热原理的方法可能会产生较大的误差， 分析烧蚀过程相关参数或分析束流的等离子体参数可能是束流参数表征的重要补充.