激光系统中元件热效应时空特性分析 下载: 860次
1 引言
热效应是光学系统设计中要考虑的重要因素之一,固体激光[1-2]、合束光栅[3]、测量系统[4]的性能均受热效应的影响。随着激光系统向高能量、高功率发展,内部光路复杂度增加,引入了更多的光学元件,光学元件受到激光辐照升温、变形,使得激光波面发生畸变,降低了激光束光束质量,影响了远场激光强度分布,因此,采用自适应光学(AO)对元件热效应进行校正,可以得到更好的光束质量。
元件热效应变化规律以及对热相差的补偿方法,国内外有很多丰富的研究成果。Peñano 等[5]研究了凹透镜、凸透镜热效应对光束评价因子M2的影响,Sun等[6]研究了主次镜在激光辐照下温升、热变形以及热相差Zernike系数变化,盛朝霞等[7-8]研究了窗口镜热效应对光束质量的影响,An等[9-10]研究了窗口镜、介质气体综合热效应分析,刘文广等[11-12]等研究了元件热效应对激光光束远场传输特性的影响,Mao等[13]研究了补偿热相差方法,周琼等[14-15]研究了反射镜、分光镜热效应相差特性与波前补偿问题。
Si反射镜、SiO2窗口镜是激光系统中常用元件,在激光辐照过程中,两者热相差的时间、空间变化规律不同,对光学系统自适应校正、传输特性也产生不同的影响。本文基于光束控制多物理耦合仿真分析软件EasyLaser[16],分析了Si反射镜、SiO2窗口镜温升和热相差分布的时间、空间分布特征,对比分析了激光辐照强度分布对元件热相差的影响;搭建仿真光路,分析不同分布激光辐照Si反射镜、SiO2窗口镜热效应的AO效果的影响,为光学系统设计与实验结果诊断、AO校正提供参考。
2 元件热相差理论
激光辐照到光学元件上,若不考虑元件材料物性随温度的变化,并忽略辐射、对流散热对元件温度分布的影响,则温度分布随时间变化的热传导方程为
忽略激光穿过透射元件的能量损失,则
元件表面吸收的激光能量作为热边界条件,即
假设元件初始温度分布均匀,则
式中:ρ、c、κ分别为材料密度、热容与热导率;q(x,y,z,t)为元件体内吸收激光能量;β(x,y,z)为元件体吸收率;η(x,y)为元件表面吸收率;I(x,y,t)为激光强度分布;T0为元件初始温度分布;L为元件厚度。
得到元件温度分布后,求解热弹性方程[17],即可得到元件变形。
元件热相差主要是由于激光辐照元件升温,导致元件变形、折射率发射变化,使得激光波前发生改变,影响激光传输特性。反射镜的热相差Δϕ1正比于元件表面变形,即
窗口镜的热相差由前后表面变形以及由折射率变化引起的透射相差引起。折射率变化主要有两个原因:1)热光效应,即温度变化引起折射率变化;2)弹光效应,即应力引起折射率变化,一般后者可以忽略[10]。因此,窗口镜热相差Δϕ2表达式为
式中:k为波矢;uz为元件表面热变形;θ为激光入射角度;n0为元件折射率;
3 元件热效应特性对比分析
Si反射镜、SiO2窗口镜是光学系统中常用元件,两种材料属性[10-11]如
高能激光采用非稳腔输出的光斑一般为环形束,如化学氧碘激光。计算中,假设环形区域内激光强度分布均匀,即平顶环形束,表示为
表 1. 元件材料物理参数
Table 1. Physical parameters of the materials
|
式中:rin和rout分别是光束的内、外半径。计算中取rout=50 mm,rin=15 mm,I0=1000 W/cm2,辐照时间为20 s。为了方便比较分析,假设Si反射镜、SiO2窗口镜厚度均为30 mm,直径为300 mm,元件对辐照激光总的吸收率均为200×10-6。Si反射镜只有迎光面受激光辐照,因此计算中取辐照面吸收率为200×10-6,入射角度为22.5°。SiO2窗口镜为透射元件,可分为前、后表面吸收与体吸收,考虑了两种情形:1)只考虑面吸收,即前、后表面吸收率均为100×10-6,体吸收率为0;2)窗口镜以体吸收为主,体吸收率为60×10-6 cm-1,前、后表面吸收率为10×10-6。
在激光辐照期间,Si镜温升随辐照时间增加而增加,且温度分布呈环状,而镜面中心凹区温升比辐照环形区域略小,如
图 1. Si反射镜温升分布。(a)表面最大温升随时间变化;(b)表面温升径向分布
Fig. 1. Temperature increase distribution of Si reflective mirror. (a) Maximal temperature increase change with time; (b) radial distribution of surface temperature increase
图 2. Si反射镜热相差RMS随时间变化
Fig. 2. RMS of thermal aberration of Si reflective mirror change with time
SiO2材料的热导率比Si材料约小2个量级,当受到激光辐照时,SiO2热扩散很慢,20 s的扩散长度约为0.8 cm,导致SiO2温升随辐照时间增长近似线性增加,且温升分布与辐照激光光斑分布相似。即使激光停止辐照后,在很长时间内,温升分布仍然与辐照激光分布相似,如
图 3. SiO2窗口镜不同时刻温升分布。(a)面吸收;(b)体吸收
Fig. 3. Temperature increase distribution of SiO2 window mirror at different moments. (a) Face absorption; (b) body absorption
SiO2窗口镜表面、体吸收的不同处主要表现在表面最大温升随时间变化不同:只考虑面吸收时,表面温升较高,但是激光停止辐照后,表面温度下降很快;以体吸收为主时,表面最大温升较低,激光停止辐照后,表面温度下降较慢。SiO2窗口镜的热相差包括表面变形与热光效应两部分,当激光的总吸收率相同时,热相差RMS相近,如
图 4. SiO2窗口镜热效应的时间变化。(a)表面最大温升;(b)热相差RMS
Fig. 4. Thermal effect of SiO2 window mirror change with time. (a) Maximal temperature increase; (b) RMS of thermal aberration
4 自适应光学校正
采用光束控制多物理耦合仿真分析软件EasyLaser搭建元件热分析及自适应光学校正仿真光路,并设置光路中元件的总吸收率相同,仿真中主要分析Si反射镜、SiO2窗口镜热效应对光束传输与自适应光学校正的影响,不考虑其他光学元件和激光传输光路中气体热效应。
AO校正原理如下:激光束通过变形镜、反射镜或窗口镜,主激光由分光镜透射后进行光束远场诊断,反射光进入哈特曼探测器对热畸变波面进行探测和重构,根据哈特曼探测器与变形镜间的共轭关系,计算变形镜驱动器的行程量,基于控制传函驱动变形镜产生与热畸变波面共轭的面形分布,达到畸变波面的自适应光学校正。若断开图中“AO控制传函”向“变形镜”的反馈信号,仿真光路表示AO开环。两个光路中AO采用同样的计算参数:哈特曼探测器为14×14子孔径,子孔径像素数为31,像素视场角为30 μrad;变形镜单元数为127;AO控制传函采用比例积分算法。
图 5. 元件热效应自适应校正仿真光路。(a) Si反射镜;(b) SiO2窗口镜
Fig. 5. Adaptive correction simulation optical path of thermal effect. (a) Reflective mirrors Si; (b) window mirror SiO2
首先考虑均匀环形光斑辐照情形。SiO2窗口镜与Si反射镜表现出不同的特性:随着时间增加,SiO2窗口镜热相差RMS值随辐照时间线性增加,而5块Si反射镜总热相差的RMS在辐照前5 s内增加很快,随后增加速度减慢;在激光辐照时间内,Si反射镜总的热相差比SiO2镜大,如
图 6. 均匀光斑辐照时,元件热相差RMS随时间变化
Fig. 6. RMS of thermal aberration under uniform intensity laser irradiation change with time
以不考虑元件热效应、激光传输到远场的光斑峰值强度为参考,得到光学自适应校正前后归一化激光远场峰值强度的变化,对比5块Si反射镜、1块SiO2窗口镜热效应激光远场传输的影响,如
图 7. 均匀光斑辐照时,远场光斑归一化峰值强度随时间变化。(a) AO开环;(b) AO闭环
Fig. 7. Normalized peak intensity of far-field under uniform laser irradiation change with time. (a) AO open; (b) AO closed
实验中,光斑强度分布不均匀,取某次实验测量光斑分布,如
图 9. 20 s时元件热相差分布。(a) Si反射镜;(b) SiO2窗口镜
Fig. 9. Distribution of thermal aberration at 20 s. (a) Reflective mirrors Si; (b) window mirror SiO2
在实验光斑辐照下,激光强度空间分布引起两种元件热相差RMS大于均匀分布激光辐照,SiO2窗口镜的增幅尤其明显,如
图 10. 实验光斑辐照时,元件热相差RMS随时间变化
Fig. 10. RMS of thermal aberration under experimental laser irradiation change with time
图 11. 实验光斑辐照时,远场光斑归一化峰值强度校正前后随时间变化。(a) AO开环;(b) AO闭环
Fig. 11. Normalized peak intensity of far-field under experimental laser irradiation change with time before and after correction. (a) AO open; (b) AO closed
模拟分析结果表明,采用AO校正元件热效应同时,采取以下措施也有助于减小元件热效应对激光传输的影响:1)采用导热性能好的元件;2)减少元件对激光的吸收率;3)提高辐照激光光斑强度均匀性。
5 结论
光学元件热效应与激光光斑分布、材料特性等密切相关,具有不同的时空特性。对于Si反射镜与SiO2窗口镜,温升与热相差的时空分布主要区别有:1)SiO2窗口镜热导率小,表面温升分布与激光光斑强度分布相似,且在激光停止辐照后仍然能保持较长时间;Si镜热导率大,温度分布平衡很快达到动平衡,因此,当吸收激光总能量相同时,SiO2窗口镜表面最大温升以及热相差RMS值均比Si反射镜高许多;2)当吸收激光总能量相同时,SiO2窗口镜面吸收、体吸收特性与表面最大温升有关,但对窗口镜热相差RMS值影响很小;3)SiO2窗口镜热相差RMS值随着激光辐照时间增加而增加,且在辐照停止后影响仍然维持较长时间;Si反射镜热相差RMS值随激光辐照时间增加先快速增加,然后增加幅度减小,且在辐照停止后影响快速下降;4)当辐照激光光斑不均匀时,光斑强度空间分布会增大元件热相差RMS,SiO2窗口镜热相差AO校正效果变差,而Si反射镜热相差的校正效果几乎不影响。
由于这两类元件温升、热相差的时空变化特征区别较大,当实验中出现由于元件热效应引起激光远场强度变化时,可以通过测量激光辐照期间、辐照停止后波面的变化,反推出影响激光系统性能的元件类型,并采取适当措施降低元件热效应的影响。
[2] Sikocinski P, Novak O, Smrz M, et al. Time-resolved measurement of thermally induced aberrations in a cryogenically cooled Yb∶YAG slab with a wavefront sensor[J]. Applied Physics B, 2016, 122(4): 1-10.
[3] 徐娇, 陈俊明, 陈鹏, 等. 960线光谱合束光栅的热畸变分析[J]. 光学学报, 2018, 38(5): 0505002.
[4] Ramette J, Kasprzack M, Brooks A, et al. Analytical model for ring heater thermal compensation in the advanced laser interferometer gravitational-wave observatory[J]. Applied Optics, 2016, 55(10): 2619-2625.
[5] Peñano J, Sprangle P, Ting A, et al. Optical quality of high-power laser beams in lenses[J]. Journal of the Optical Society of America B, 2009, 26(3): 503-510.
[6] Sun Y, Li F, Yang W Q, et al. Thermal effects of optical antenna under the irradiation of laser[J]. Proceedings of SPIE, 2017, 10463: 104631F.
[7] 盛朝霞, 王再军. 强激光输出窗口热行为对光束质量的影响[J]. 激光技术, 2008, 32(3): 278-280.
[9] An J Z, Li Y K, Du X W. Thermal effects of a laser beam tube consisting of a window and nonflowing gas[J]. Optics Letters, 2004, 29(24): 2899-2901.
[10] 安建祝, 李有宽, 杜祥琬. 激光窗口热透镜效应对光束质量的影响[J]. 强激光与粒子束, 2004, 16(4): 429-433.
[11] 刘文广, 饶鹏, 华卫红. 非均匀激光辐照下硅镜热变形对光束传输特性的影响[J]. 强激光与粒子束, 2008, 20(10): 1615-1619.
Liu W G, Rao P, Hua W H. Effects of thermal distortion of Si mirror irradiated by non-uniformity laser intensity on laser propagation[J]. High Power Laser and Particle Beams, 2008, 20(10): 1615-1619.
[12] Guo R H, Chen N, Zhuang X Y, et al. Simulations of far-field optical beam quality influenced by the thermal distortion of the secondary mirror for high-power laser system[J]. Proceedings of SPIE, 2015, 9255: 92553G.
[13] Mao Y J, Li S K, Sun G, et al. The thermal aberration analysis of a lithography projection lens[J]. Proceedings of SPIE, 2017, 10147: 101471P.
[14] 周琼, 姜宗福, 习锋杰. 内光路系统的缩放模型及其热变形像差的波前预补偿[J]. 中国激光, 2012, 39(4): 0402008.
[15] 周琼, 姜宗福, 习锋杰. 光学系统中白宝石分光镜的热变形像差分析[J]. 中国激光, 2012, 39(10): 1002001.
[16] 谢晓钢, 张建柱, 岳玉芳, 等. 激光系统组件化仿真软件EasyLaser[J]. 强激光与粒子束, 2013, 25(10): 2536-2540.
[17] 严宗达, 王洪礼. 热应力[M]. 北京: 高等教育出版社, 1993: 100.
Yan ZD, Wang HL. Heat stress[M]. Beijing: Higher Education Press, 1993: 100.
Article Outline
胡鹏, 张建柱. 激光系统中元件热效应时空特性分析[J]. 光学学报, 2020, 40(20): 2014001. Peng Hu, Jianzhu Zhang. Analysis of Spatio-Temporal Characters of Thermal Effects of Optical Components in Laser System[J]. Acta Optica Sinica, 2020, 40(20): 2014001.