Study on thermal response of aluminum alloy plate subjected to airflow and laser irradiation
-
摘要: 为了获得激光加工等应用中气流环境激光辐照材料的热响应特性和熔穿规律,采用面热源近似和不定常热传导分析方法构建了包括热传导、对流换热、熔化烧蚀等主要机制的物理模型,并通过实验校核后进行了气流环境激光辐照铝合金热学响应定量评估,获得气流环境不同激光功率密度下铝合金靶板温升和烧蚀熔穿规律。结果表明,相同靶参数下,熔穿时间随着激光功率密度增加而急剧减小,并且熔穿时间变化规律与热平衡积分方法气化烧蚀模型结果一致;不同条件下熔穿时间显示,在激光功率密度较小(500 W/cm2附近)时,气流引起对流换热因素对激光烧蚀熔穿影响较大,而在激光功率密度较高(1500 W/cm2以上)时,气流引起对流换热因素对激光烧蚀熔穿影响较小。建立的气流条件下激光辐照热响应模型与实际物理过程更接近,计算获得的熔穿规律为激光辐照模型合理简化提供了定量的参考依据。Abstract: In order to obtain the thermal response and melting through characteristics of materials irradiated by high energy laser with airflow environment in field of laser processing and other applications, a physical model considering the mechanisms of heat conduction, convective heat transfer, melting and ablation was constructed by using the approximation of plane heat source and unsteady method. The model was validated by experiments and was used for evaluating quantitatively the thermal responses of material. Then, the temperature rise and penetration law of aluminum alloy plates subjected to different laser power density in airflow were obtained. The results indicate that the melting through time decreases sharply with the increase of laser power density under the same target parameters. And the calculated melting through time is consistent with the results by using the gasification ablation model of thermal balance integration method in the literature. The analysis of perforation time under different conditions shows that the airflow factor has a significant influence on laser ablation melting through when the laser power density is small (approximately around 500 W/cm2). When the laser power density is large (above 1500 W/cm2), the effect of airflow through convective heat transfer is insignificant. In addition, the physical model considered in calculation is close to the actual physical process. Above calculation results and analytical understanding provide a quantitative basis for simplifying reasonably the laser ablation models.
-
Keywords:
- laser technique /
- laser irradiation /
- airflow /
- thermal response /
- laser ablation
-
0. 引言
奇异光学是现代光学的重要分支[1-2]。自1973年NYE等人提出波在传输过程中波前存在缺陷[3],光学中的奇异现象就引起了研究人员的兴趣。1992年,ALLEN及其同事基于拉盖尔-高斯光束提出了光的轨道角动量[4],引起了研究奇异光学的热潮。在奇异光学中,相位奇异是指标量光场相位为涡旋相位,表现为光场中存在强度为0的点,围绕该点的封闭路径的相位变化是2π的整数倍,该倍数即为拓扑荷。对于相位奇异,其拓扑荷通常为整数,但是也存在为非整数的情况。BERRY从理论上预测了分数涡旋相位的相位阶梯不连续将导致奇异线,并在其传播过程中会产生一系列复杂的演化[5]。LEACH等人在实验中观察到了该现象[6]。GBUR则发现这种复杂的演化与希尔伯特旅馆悖论密切相关[7]。XIONG等人利用连续镜面反射镜实现了对分数涡旋光束的动态调控[8]。
此外,对于矢量光场,存在着偏振奇异现象[9]。对于2维非均匀椭圆偏振光场,偏振奇异是指光场中无法定义其偏振方向的圆偏振点(C-points)或无法定义其偏振旋向的线偏振线(L-lines)。对于C点,其拓扑荷为±1/2。此外,还有光场中均为线偏振,但偏振方向成径向或旋向分布的柱矢量光束(V-points)和高阶偏振奇点。目前,偏振奇异光场可以由基于拉盖尔-高斯光束[10]、贝塞尔-高斯光束[11]、马修-高斯光束[12]、圆艾里光束[13]、分数涡旋光束[14]等新型光场得到。研究人员已经将偏振奇异应用于光学显微操作[15]、光学微加工[16]和光学通信等方面[17]。同时,偏振奇异还带来了很多有趣的物理现象,比如伪拓扑性[18]、奇点结19]、奇点爆炸[20]以及C点偶极子的产生与湮灭等[21]。对于V点,其拓扑荷通常为整数,然而ZHANG等人研究了V点拓扑荷为分数的情形,可以观察到偏振奇异线[22]。此后,他们又创造实现了Julia分数矢量光场,并验证了其伪拓扑性[18]。GU等人从理论和实验两方面研究了弱聚焦分数矢量光束的自旋角动量分离和传输特性[23]。
另一方面,无衍射光束因其独特的传输特性近些年来也越来越受到研究人员的关注,而其中艾里光束更是由于其加速性和自愈性吸引了大量的研究[24]。基于艾里光束,EFREMIDIS等人又提出了圆艾里光束[25],而后研究人员将其与涡旋相位结合,提出了艾里涡旋光束,并研究了其传播性质[26-28]。虽然C点的拓扑荷为分数,但是其通常基于整数拓扑荷的涡旋光束生成得到,对基于分数涡旋光束的偏振奇异光场虽已提出,但尚未研究其传输特性。本文中拟用携带分数涡旋相位的圆艾里光束来进行偏振奇异光场的传播特性研究。利用两束正交圆偏振的圆艾里光束的叠加,其中一束拓扑荷为0,另一束拓扑荷为分数,可以得到基于分数涡旋光束的偏振奇异光场,并且可以控制偏振奇异光场的拓扑荷并非传统的C点的±1/2,圆艾里光束的自聚焦特性会影响该光场在自由空间传播中的拓扑结构。本文中的研究展示了基于圆艾里光束的分数偏振奇点在传播过程中的拓扑结构演变,为其应用于自由空间光通信做好了理论基础,从而拓宽了圆艾里光束的应用场景。
1. 基于分数圆艾里涡旋光束的偏振奇异光场
为生成得到偏振奇异光场,采用了生成全庞加莱光束的方法[10],即叠加两束正交偏振的圆偏振光。位于初始平面z=0的叠加光场可以表示为[10]:
\boldsymbol{E}(\boldsymbol{r} ; \gamma)=\cos \gamma \cdot \boldsymbol{e}_1 U_1(\boldsymbol{r})+\sin \gamma \cdot \boldsymbol{e}_{\mathrm{r}} U_2(\boldsymbol{r}) (1) 式中: r为空间向量; γ为控制光束强度的常量参数; el和er分别为左旋和右旋圆偏振基矢量; U1(r)和U2(r)分别为携带或不携带涡旋相位的圆艾里光束。位于初始平面的携带涡旋相位的圆艾里光束,其电场在柱坐标系(r, φ, z)中可表示[26]:
\begin{aligned} U(r, \varphi, z=0)= & A_0 \cdot A\left(\frac{r_0-r}{w_0}\right) \exp \left(a \frac{r_0-r}{w_0}\right) \cdot \\ & {[r \exp (\mathrm{i} \varphi)]^l } \end{aligned} (2) 式中: A0为振幅常量; A(·)为艾里函数; r0为主艾里光环的半径; w0为径向缩放; a为决定传播距离的衰变参数; l为涡旋相位的拓扑荷。对于l=0,即代表圆艾里光束不携带涡旋相位。
为生成分数偏振奇异光场,将U1设置为携带涡旋相位且拓扑荷l为分数数值的圆艾里光束,U2为不携带涡旋相位的圆艾里光束。设定光场U1中l的数值为0.1、0.3、0.5、0.7、0.9和1.0,设定γ=π/4,r0=1 mm,w0=0.08 mm,a=0.1,光场尺寸为4 mm×4 mm。将上述参数结合式(1)和式(2),即U1为左旋圆偏振、U2为右旋圆偏振时,初始平面光场的偏振态分布如图 1所示。图中偏振椭圆的不同颜色代表了不同的旋向,绿色为右旋偏振,红色为左旋偏振。图 1展示了左旋分量U1中的拓扑荷不同时,光场偏振态分布的演化过程。当l=0.1时,光场偏振态接近于线偏振;随着l数值的不断变大,光场中偏振椭圆的椭圆率也随之变大;当l=1.0时,光场中出现了一个圆偏振点,即C点,此时光场的偏振态分布是一个典型的“柠檬”结构;当l取其它数值时,虽然其偏振态分布与“柠檬”结构相似,但是其光场中并没有圆偏振点,可将其称为“准柠檬”结构。
与之类似,将U1设置为右旋圆偏振、U2为左旋圆偏振时,初始平面光场的偏振态分布如图 2所示。此时初始平面的偏振态从接近于线偏振逐渐变化到“星”结构;当l取其它数值时,其偏振态分布可称为“准星”结构。
l>1.0时的情形同样值得关注,以U1为左旋圆偏振为例,当拓扑荷l取值为1.1、1.3、1.5、1.7、1.9和2.0时,初始平面的偏振态分布如图 3所示。结合图 1中l=1.0的情形,图 3展示了光场从“柠檬”到高阶偏振奇异光场的演化过程。
2. 偏振奇异光场的动态传输
对于圆艾里光束,无论是否携带涡旋相位,其在自由空间传输时的光场都无法用解析式U(r, φ, z)来表示,因此只能用数值模拟的方式来展现其在某特定传输距离z时的光场分布。在此采用角谱的方法来进行计算,即在初始平面通过快速傅里叶变换将其变换到频域,通过乘以传递函数,得到传输距离z处的频谱,再经过快速傅里叶逆变换,得到其在z处空域的光场分布[29]。对于传输距离z,分别选取为100 mm,400 mm,464 mm和500 mm,可以数值模拟U1和U2在各传输距离的光强分布。其中对于U1,分别选取l=0.5和l=1.5为例,如图 4所示。可以看到,无论是否携带涡旋相位,也无论拓扑荷l的数值,随着传播距离z的不断变大,U1和U2主光环的尺寸在不断减小,这印证了圆艾里光束的自聚焦特性。在传输距离z=464 mm时,光场尺寸达到最小值,而后又随着传输距离不断变大,即二者的焦距均为464 mm。
根据图 4中两个正交圆偏振分量的光场分布,可以计算得出在传输距离z处叠加光场的偏振态分布,如图 5所示。需注意的是, 图 5中第1列光场尺寸是3 mm×3 mm,而后3列的光场尺寸为1 mm×1 mm。从图示来看,光场中偏振椭圆的分布规律并不明显,将在下一节中讨论这一问题。
3. 分析与讨论
为研究光场的拓扑结构,斯托克斯相位是一个非常重要的参量[30]。斯托克斯相位场中的相位奇异点对应矢量光场中的偏振奇异点,其表达式为[30]:
S_{12}=S_1+\mathrm{i} S_2 (3) 式中: S1、S2分别为斯托克斯参量; i为虚数单位。斯托克斯相位ϕ=arctan(S2/S1),对于U1, l=0.5+U2, l=0, 在不同传输距离z处的斯托克斯相位场计算如图 6所示。在初始平面,斯托克斯相位从-π/2变化到π/2;而在其传输过程中,其相位的变化均为从-π到π。在传输过程的初段,可以看到其中心涡旋结构遭到破坏,经过一段距离之后,如z为464 mm和500 mm处,发现相位场中又至少出现了一个相位奇异点,如图中红色圆圈所示。
如果将传输距离z=464 mm和z=500 mm两处的光场尺寸缩放到0.2 mm×0.2 mm,这个局部区域的光场偏振态分布如图 7所示。图中黄线表示S1=0,蓝线表示S2=0,两者交点,即图中实心红圆即为偏振奇异点。也就是说在初始平面不存在偏振奇异点的情况下,经过一定距离的传输,光场中出现了偏振奇异点。
对于U1, l=1.5+U2, l=0的情形,不同传输距离的斯托克斯相位场如图 8所示。由图可以看出,在初始平面,斯托克斯相位场仅包含一个相位奇异点,而在传输距离z=464 mm和z=500 mm处,至少都包含两个相位奇异点。
同样地,可以得出在z=464 mm和z=500 mm处中心0.2 mm×0.2 mm区域的偏振态分布,如图 9所示。结合图 8可以看出,在传输初期拓扑结构遭到破坏的情况下,经过一定的传输距离之后依旧可以重建,并且出现了更多的偏振奇异点。结合图 7可以看出,基于圆艾里光束的分数偏振奇点在空间中传播时拓扑结构的可恢复性,利用这一特性可将该偏振奇点用于自由空间光通信中。
4. 结论
由于圆艾里光束的自聚焦特性,叠加光场在传播过程之初,主瓣半径不断减小直至其焦点;由于分数涡旋光束的传输特性,叠加光场在传输之初存在缺口,然后得以闭合。无论其在初始平面有无C点,当叠加光场传输到焦点时,光场中出现了C点,并在焦点之后的传输过程中将出现更多的C点。与基于整数圆艾里涡旋光束的偏振奇异光场相比[13],两者均具有奇异拓扑结构的自恢复特性,然而基于分数的涡旋光束会给恢复后的光场中带来更多的C点。该工作从理论和数值模拟研究了对基于分数涡旋光束的偏振奇异光场在自由空间传输时拓扑结构的演变规律,为其应用提供理论基础。
-
表 1 模型校验用算例参数
Table 1 Parameters of simulation case for model verification
number power density/(W·cm-2) irradiation time/s material material dimension/mm flow velocity/(m·s-1) case 1 1190 10 LY12 \varnothing30×2 0 case 2 1190 10 LY12 \varnothing30×1 0 shot 1 800 10 LY12 100×100×3 120 shot 2 1400 10 LY12 100×100×3 120 -
[1] 孙承纬. 激光辐照效应[M]. 北京: 国防工业出版社, 2002: 1-50. SUN Ch W. Laser irradiation effect[M]. Beijing: National Defense Industry Press, 2002: 1-50(in Chinese).
[2] WANG P, ZHANG Z, HAO B, et al. Investigation on heat transfer and ablation mechanism of CFRP by different laser scanning directions[J]. Composites Part, 2023, B262: 110827.
[3] WAN H L, LIN J P, MIN J Y. Effect of laser ablation treatment on corrosion resistance of adhesive-bonded Al alloy joints[J]. Surface and Coatings Technology, 2018, 345: 13-21. DOI: 10.1016/j.surfcoat.2018.03.087
[4] WANG Y, ZHANG M Y, DONG Y H, et al. Morphology modeling and validation in nanosecond pulsed laser ablation of metallic materials[J]. Precision Engineering, 2023, 79: 34-42. DOI: 10.1016/j.precisioneng.2022.08.010
[5] DOU H Q, YAO C Zh, LIU H, et al. Femtosecond laser ablation of Al-Mg alloy in vacuum and air[J]. Applied Surface Science, 2018, 447: 388-392. DOI: 10.1016/j.apsusc.2018.04.003
[6] LI P F, WANG Y H, LI L L, et al. Ablation oxidation and surface quality during laser polishing of TA15 aviation titanium alloy[J]. Journal of Materials Research and Technology, 2023, 23: 6101-6114. DOI: 10.1016/j.jmrt.2023.02.209
[7] YIN J H, CAO Y Zh, CUI Y W, et al. Nd ∶YAG laser ablation of aluminum alloy 6061 before and after silicon dioxide coating[J]. Journal of Alloys and Compounds, 2021, 877: 160329. DOI: 10.1016/j.jallcom.2021.160329
[8] 焦路光, 赵国民, 江厚满. 切向气流作用下激光对典型金属靶的辐照效应[J]. 中国光学, 2011, 4(1): 77-81. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGA201101015.htm JIAO L G, ZHAO G M, JIANG H M. Irradiation effects of laser on typical metal targets under tangential airflow[J]. Chinese Optics, 2011, 4(1): 77-81(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGA201101015.htm
[9] 张潇允, 张巍, 夏盛强, 等. 高功率激光辐照CFRP的温度场和应力场的数值分析[J]. 激光技术, 2021, 45(5): 636-641. DOI: 10.7510/jgjs.issn.1001-3806.2021.05.017 ZHANG X Y, ZHANG W, XIA Sh Q, et al. Study on numerical analysis of temperature field and stress field of carbon fiber reinforced polymers irradiated by high power laser[J]. Laser Technology, 2021, 45(5): 636-641(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2021.05.017
[10] XING X D, MA T, WANG R X, et al. Dynamic rupture of metal sheet subjected to laser irradiation and tangential subsonic airflow[J]. Theoretical & Applied Mechanics Letters, 2018, 8(4): 272-276. DOI: 10.3969/j.issn.2095-0349.2018.04.008
[11] WANG J, MA Y, LIU Y, et al. Experimental investigation on laser ablation of C/SiC composites subjected to supersonic airflow[J]. Optics & Laser Technology, 2019, 113: 399-406. DOI: 10.5768/JAO201940.0301007
[12] 伍俊英, 杨利军, 吴宝, 等. 强激光烧蚀铝靶实验及数值模拟[J]. 北京理工大学学报, 2018, 38(10): 1018-1024. https://www.cnki.com.cn/Article/CJFDTOTAL-BJLG201810005.htm WU J Y, YANG L J, WU B, et al. Experiment and numerical simulation of ablation on aluminum target with high-intensity laser[J]. Transactions of Beijing institute of Technology, 2018, 38(10): 1018-1024(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-BJLG201810005.htm
[13] BAEK W K, LEE K C, AN S I, et al. Melt-through characteristics in continuous beam irradiation of flying metal samples in flow speeds up to 85 m/s[J]. Optics & Laser Technology, 2013, 45: 250-255.
[14] 张萌, 方英武, 张广鹏. 纳秒脉冲激光辐照铝靶碎片动态响应的仿真研究[J]. 激光技术, 2023, 47(4): 541-546. DOI: 10.7510/jgjs.issn.1001-3806.2023.04.015 ZHANG M, FANG Y W, ZHANG G P. Numerical simulation of dynamic response for aluminum target debris irradiated by nanosecond pulse laser[J]. Laser Technology, 2023, 47(4): 541-546(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2023.04.015
[15] ZHAO W N, MA T, SONG H W, et al. Effects of tangential supersonic airflow on the laser ablation of laminated CFRP[J]. Journal of Materials Research and Technology, 2021, 14: 1985-1997. DOI: 10.1016/j.jmrt.2021.07.101
[16] 张永强, 陶彦辉, 张黎, 等. 激光辐照超音速气流下TA15钛合金和LY12铝合金的热响应[J]. 强激光与粒子束, 2014, 26(8): 081005. DOI: 10.11884/HPLPB201426.081005 ZHANG Y Q, TAO Y H, ZHANG L, et al. Thermal response of TA15 titanium alloy and LY12 Aluminium alloy irradiated by laser under supersonic tangential flow[J]. High Power Laser and Particle Beams, 2014, 26(8): 081005 (in Chinese). DOI: 10.11884/HPLPB201426.081005
[17] HARRACH R J. Analytical solutions for laser heating and burnthrough of opaque solid slabs[J]. Journal of Applied Physics, 1977, 48(6): 2370-2383. DOI: 10.1063/1.323999
[18] BOLEY C D, CUTTER K P, FOCHS S N, et al. Interaction of a high-power laser beam with metal sheets[J]. Journal of Applied Physics, 2010, 107(4): 043106. DOI: 10.1063/1.3284204
[19] 梅丽芳, 秦建红, 严东兵. 活性激光焊接304不锈钢温度场的数值与试验研究[J]. 激光技术, 2020, 44(4): 492-496. DOI: 10.7510/jgjs.issn.1001-3806.2020.04.016 MEI L F, QIN J H, YAN D B. Numerical and experimental study on temperature field of activated laser welding 304 stainless steel[J]. Laser Technology, 2020, 44(4): 492-496(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2020.04.016
[20] 王伟平, 刘仓理, 王春彦, 等. 切向气流对激光加热材料的影响[J]. 强激光与粒子束, 1996, 8(3): 373-377. https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY603.010.htm WANG W P, LIU C L, WANG Ch Y, et al. Tangential airflow influence on laser heating materials[J]. High Power Laser and Particle Beams, 1996, 8(3): 373-377(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY603.010.htm
[21] 张家雷, 谭福利, 刘仓理, 等. 激光熔蚀金属平板的三维有限体积法数值分析[J]. 强激光与粒子束, 2015, 27(8): 081002. https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY201508003.htm ZHANG J L, TAN F L, LIU C L, et al. 3D finite volume method simulation on melt removal effects of metal plate under laser irradiation[J]. High Power Laser and Particle Beams, 2015, 27(8): 081002(in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-QJGY201508003.htm
-
期刊类型引用(1)
1. 冯姣姣,段美玲,单晶,王灵辉,薛婷. 部分相干混合位错光束在生物组织传输中的偏振特性. 物理学报. 2024(18): 137-150 . 百度学术
其他类型引用(0)