水束导引高功率激光的热效应分析.pdf

第 47 卷 第 5 期2023 年 9 月激 光 技 术LASER TECHNOLOGYVol.47,No.5September,2023 文章编号:1001-3806(2023)05-0672-06水束导引高功率激光的热效应分析赵 臻,龙芋宏,黄宇星,张光辉,周 辽,蔺 泽,焦 辉,黄 平(桂林电子科技大学 机电工程学院,桂林 541004)摘要:为了探究水束导引高功率激光的热损失问题,调控水导高功率激光中水-光束耦合的热稳定性,保证激光能量有效地传输至工件表面,采用有限元方法建立水-光束耦合的数值模型,利用射线追踪模拟不同入射激光功率、射流流速和直径时的水射流温度分布,并通过实验数据验证模型的有效性,对比了相同功率下连续激光与脉冲激光的温度分布。结果表明,水射流的温度随着水-光束耦合长度增加而升高,当入射激光功率增大、耦合腔压力降低以及射流直径减小时,都会导致射流的温度升高,严重影响射流的稳定性;由于水射流的冷却作用,平均功率 300 W 的脉冲激光相较于连续激光,温度相差约 5。此研究结果为控制水导高功率激光中水-光束耦合的热效应、提高激光能量的传输效率提供了一定的参考。关键词:激光技术;水导高功率激光;射线追踪;热效应;传输效率中图分类号:TN249;TG665 文献标志码:A doi:10.7510/jgjs.issn.1001-3806.2023.05.015Thermal effect analysis of water-jet guided high-power laserZHAO Zhen,LONG Yuhong,HUANG Yuxing,ZHANG Guanghui,ZHOU Liao,LIN Ze,JIAO Hui,HUANG Ping(School of Mechanical and Electrical Engineering,Guilin University of Electronic Technology,Guilin 541004,China)Abstract:In order to control the thermal stability of water-laser coupling in water-jet guided high-power laser and to ensure the effective transmission of laser energy to the workpiece surface,the thermal loss problem of water-jet guided high-power laser was investigated.The finite element method was used to establish a numerical model of water-laser coupling,and ray tracing was used to simulate the temperature distribution of water jets with different incident laser power,jet velocity,and diameter.The models effectiveness was verified by experimental data,and continuous and pulsed laser temperature distribution was compared under the same power.The results show that the water jet s temperature increases with the water-laser coupling length.When the incident laser power increases,the pressure of the coupling cavity decreases,the jet diameter decreases,and the temperature of the jet increases,which seriously affects the stability of the jet.Due to the water jet s cooling effect,the pulse laser s temperature difference with an average power of 300 W is about 5 compared with that of the continuous laser.The research results provide certain reference significance for controlling the thermal effect of water-laser coupling in water-jet guided high-power laser and improving the transmission efficiency of laser energy.Key words:laser technique;water-jet guided high-power laser;ray tracing;thermal effect;transmission efficiency基金项目:国家自然科学基金资助项目(62004050;52165056);广西自然科学基金资助项目(2019JJD160010;2020JJB170048);广西科技基地和人才专项基金资助项目(2021AC18026);广西制造系统与先进制造技术重点实验室项目(17-259-05-018Z);广 西 青 年 教 师 教 育 项 目(2020KY05020);桂林电子科技大学研究生优秀学位论文培育项目(YCBZ2021073);桂林电子科技大学研究生教育创新计划资助项目(2021YCXS001);广西研究生教育创新计划资助项目(YCSW2022287;YCBZ2022114;YCBZ2021073)作者简介:赵臻(2000-),男,硕士研究生,现主要从事水导激光加工技术的研究。通讯联系人。E-mail:1985251027 收稿日期:2022-08-15;收到修改稿日期:2022-10-24引 言水导激光加工技术在 1994 年被首次提出,此后被应用于半导体、太阳能、电子和生物医疗等各种领域1-4。这项技术利用全反射的原理,将穿过聚焦透镜聚焦后的高功率激光约束在微米级的水射流中,经水射流导引至工件表面,高功率激光对材料烧蚀加工,水射流可冲刷工件表面带走熔渣并冷却工件5。该技术的优点包括:(1)水导激光加工技术不需要对焦,加工长度取决于水射流的长度,有效加工距离可达50 mm 以上;(2)水射流可将加工过程中产生的熔渣带走,清洁工件表面;(3)得益于水射流的冷却作用,第 47 卷 第 5 期赵 臻 水束导引高功率激光的热效应分析加工过程中产生的热影响区可大大减少;(4)水射流类似于多模光纤的作用,改善激光的能量分布,降低加工锥度6。随着制造业对加工效率、加工质量提出更高的需求,水射流耦合的激光功率不断提升,水导激光中常常采用 1064 nm、532 nm 波长激光器。由于 1064 nm 波长的激光价格便宜、稳定性好、通用性强、高功率激光更容易实现等特点,因此在水导激光技术中更常用1064 nm 作为激光源,但是水对 1064 nm 波长激光的吸收系数较大,导致了水-光束耦合中的热效应不容忽视。当高功率激光经过光路传输后,聚焦于耦合腔中的喷嘴表面,会在焦点附近产生热聚集现象,引起耦合腔水温的升高,导致水的折射率变化,使激光的焦点发生偏移7-9。高功率激光与水射流成功耦合后,由于水对激光能量的吸收作用,导致水束的温度升高,同时温度会影响水对激光能量的吸收系数,从而在水束的上、下游形成温度梯度,影响射流表面扰动的生长10-12。扰动生长到一定程度,造成射流的破碎,此时激光从水束中散射出来,造成激光能量的极大损失,无法形成有效加工13。激光功率密度达到一定阈值,光子将与水分子发生非弹性碰撞,造成入射光子频率的变化,这种现象称为受激 喇 曼 散 射14。SPIEGEL 等 人15将532 nm 的激光耦合进入水束,观察到了红光从水束的稳定区发射出来,测得能量的传输效率降低了 30%。目前国内外对于激光与物质相互作用的研究主要集中在激光与材料方面,对水导激光技术中高功率激光与射流相互作用的研究极少,而水-光束耦合的热稳定性又是高功率激光能量传输的关键,因此,有必要对水-光束耦合中的热效应进行系统的研究。基于以上问题,本文中建立了高功率激光与水射流相互作用温度场的数值模型,利用 COMSOL 多物理场仿真软件模拟了不同激光功率、射流流速和直径下水-光束耦合的温度场分布。对比了相同功率下连续激光与脉冲激光的温度分布,同时,求解了常用耦合腔压力及射流直径参数下的最大容忍激光功率。保证水-光束耦合的热稳定性,提高了激光能量在水射流中的传输效率,为水导高功率激光加工的热效应控制提供一定借鉴。1 水束耦合高功率激光束的耦合传输特性1.1 激光在水射流中的耦合传输特性激光与水射流耦合时,激光的光斑直径需要小于水射流的直径,激光的发散全角需要小于水射流所能接收光的最大角度。耦合成功后激光在水射流中传输,由于分子的吸收、微粒的散射等作用,会产生激光的传输损耗,当传输的距离较短时,与在大气中传输一样,激光能量的衰减符合比尔-朗伯定律:P=P0exp(-L)(1)式中,P0、P 分别是水-光束耦合长度为 0 和 L 时的激光功率,为包括吸收和衰减在内的衰减系数。1.2 水的热力学性质为了研究高功率激光与水射流之间的相互作用,有关水的热力学性质将不可忽略,水射流吸收激光的能量,造成温度的升高。不同波长的激光在水中的吸收系数随温度的变化不同,当激光波长在 0.9 m 1.0 m 和 1.07 m1.18 m 范围内,吸收系数随着温度的升高而升高,在 1.0 m1.07 m 波长范围内呈相反的变化趋势。本文中采用 1064 nm 波长的激光,温度 T 与衰减系数 之间的关系由下式给出16:(T)=CT1+exp(-B/T)(2)式中,常数 B 与 C 的值分别为 360 K 和 52.8 K/cm。2 水束耦合高功率激光束的数值模型激光的能量被水吸收从而引起射流温度的变化,为了得到整个射流域上的温度分布,将水射流假设为一圆柱体,划分为有限元单元,如图 1 所示。每个有限元单元都会对传输经过的激光产生吸收作用,从而在单元内生成焦耳热。由于水射流从高压耦合腔的喷嘴射出,具有上百米每秒的速度,因此单元与周围空气之间的对流换热损失也应该被计算,解析每个单元上的温度变化,通过迭代的方式得到整个射流域上的温度变化。图 1 圆柱水射流有限元划分Fig.1 Finite element division of cylindrical water jet如图 1 所示,在射流的传播方向上将圆柱(长度为 ljet)射流划分为有限个厚度为 dz、直径为 djet的单元,每个单元都遵循热力学第一定律,T(z)和 T(z+dz)分别表示 P(z)和 P(z+dz)对应的温度。由于在单元376激 光 技 术2023 年 9 月上并无外界系统做功,内能的变化只能通过热能的形式在单元的边界上传递,根据热力学第一定律,单元内能的增加 dU 可以表示为:dU=dQ+dW=dQ(3)式中,dW 为外界系统对单元所做的功,此处 dW=0;dQ 为外界与单元之间传递的热量,包括激光通过厚度为 dz的单元相互作用损失的激光功率 dP 和对流换热损失的热通量,而 dP 可以由比尔-朗伯定律给出,结合(1)式、(2)式得:dP=P(z)-P(z+dz)=P(z)1-exp-CdzT1+exp(-B/T)(4)射流经喷嘴喷出后被划分为紧密段、过渡段和破碎段,紧密段的水射流通常被用来导引激光。在水射流的紧密段内,可以认为该段截面上的流速处处相同,均为水射流在喷嘴处的出口速度17,并由下式给出:vjet=(2p/)cv(5)式中,p 为耦合腔内外压强差,为水的密度,cv为速度系数,用于补偿摩擦和紊流引起的能量耗散所造成的理论射流速度的降低。假设射流周围的空气静止,根据牛顿冷却公式,当水射流与周围空气发生对流换热时,单元边界上的热通量由下式描述:=qS=h(T-Ta)S(6)式中,q 为边界上的热流密度,S=djetdz为单元的对流换热面积,h 为对流换热系数,T 为单元的温度,Ta为周围空气的温度。由于水射流的缩流效应18-19,水射流的直径小于喷嘴的直径,可得到水射流直径与喷嘴直径之间的关系:djet=kdn(7)式中,dn为喷嘴直径,k 为缩流系数。对流换热系数 h描述了水射流的温度下降,通过努塞尔数 Nu、水射流直径 djet及液体导热系数 计算得到:h=Nudjet(8)因此,对流换热系数主要由努塞尔数 Nu 决定,在水射流从喷嘴喷出流向下流的过程中,Nu 的经验公式为:Nu=0.018Re0.8Pr0.3(9)式中,雷诺数 Re=vjetdjet/,普朗特数 Pr=cpvjet/,为水的动力粘度,cp为水的质量定压热容,为水的导热系数。根据(3)式,传递到单元上的热量等于激光与单元相互作用期间的功率损失与对流换热的热量损失之差:dU=cpmdT=(dP-)dt(10)式中,m=djet2dz/4 为单元的质量,激光与单元的相互作用时间 dt=dz/vjet定义为以射流速度通过 dz的时间,于是单元上下界面的温度差 dT 为:dT=dP-cpdjet2vjet/4(11)在本文模型中,水的定压比热容及密度被认为是常数。3 仿真结果与模型验证本文中均使用 COMSOL 多物理场仿真软件模拟不同参数下的射流截面温度分布。为了简化起见,不考虑水射流表面的波动,水射流的模型被建立为 3 维圆柱体,激光由射线的形式释放,利用射线追踪的方法,经聚焦透镜聚焦至射流顶面,由于射流与周围空气的折射率存在差值,射流被约束在射流内,并在水-气界面发生全反射。在传输的过程中,激光束由于水的吸收,能量沉积在射流中并转换为焦耳热,导致射流温度的升高。同时,由于水压的存在,经喷嘴射出的射流具有较快的速度,因此,本文作者将高速水射流与周围空气之间的对流换热考虑在内。本文中,假设射流周围的空气静止,且射流以恒定的速度流动。根据上述的数值模型可知,影响温度变化的主要因素为激光的入射功率、射流流速以及射流直径的大小。本文作者首先对入射激光功率为 750 W、1000 W、1250 W 及1500 W 下的温度分布进行了模拟,固定耦合腔压力为5 MPa,喷嘴直径 250 m,射流长度 50 mm。图 2 曲线为稳态时射流中心轴线上的温度。可以看到,不管入射激光功率如何变化,随着水-光束耦合长度的增加,射流的温度总是呈现上升的趋势,这是由于在耦合长 图 2 不同入射激光功率的水射流温度分布Fig.2Temperature distribution of water jets with different incident laser powers476第 47 卷 第 5 期赵 臻 水束导引高功率激光的热效应分析度增长的过程中总的热量传递增加,导致射流温度的升高。当入射激光功率为 750 W 时,被水射流吸收的激光功率较少,因此射流的温度变化较小,在稳态时,最高温度为 46。随着激光功率的增大,沉积在射流内的功率增多,生成了更多的焦耳热,射流的温度也逐渐增大,当激光功率为 1500 W 时,射流最高温度达到了 67 左右,温度的升高导致了射流表面扰动的增长速率加快,这对于射流的稳定性是极为不利的。为了验证模拟的有效性,在参数完全一致的情况下,将模拟结果与 CHRISTIAN 等人12的实验结果进行了对比。由图 3 可以看到,随着入射激光功率的增大,模拟值(线条)与实验值(散点)都呈现增长的趋势,并且模拟得到的射流温度与实验测得的温度数值较为符合,经计算平均相对误差约为 3%。图 3 模拟值与实验值的对比Fig.3 Data comparison of the simulated values and the experimental values同样地,本文中也模拟了不同射流流速和不同射流直径的温度分布,根据(5)式,当喷嘴直径不变时耦合腔压力决定了射流的流速,通过改变耦合腔的压力来改变流速,对耦合腔压力为 5 MPa、10 MPa、15 MPa和 20 MPa 时的射流温度分布进行了模拟。如图 4 所示,随着耦合腔压力的增大,水射流的温度逐渐降低,这是因为射流速度的增加使得射流与周围空气之间的换热加快,带走了更多的热量,当耦合腔压力为20 MPa 时,射流的最高温度为 47 左右,温度的降低 图 4 不同耦合腔压力的水射流温度分布Fig.4Temperature distribution of water jets different coupling chamber pressures能有效保证水-光束耦合的热稳定性,有利于激光能量的传输,这一变化趋势也与 DENG 等人20的相关研究符合。在本文中,研究了射流直径对温度分布的影响,作者对直径为50 m、100 m、150 m 及200 m 的射流进行 温 度 模 拟,如 图 5 所 示。当 射 流 直 径 为50 m 时,射流的最高温度远远超过了水的蒸发温度,导致射流的稳定性降低甚至破碎,从而无法约束水射流内的激光,造成激光能量的衰减,传输到工件上的激光能量过低无法实现有效加工。并且最高温度出现在水-光束耦合长度为 25 mm 处,在此之后,随着耦合长度的增加,射流的温度逐渐降低,这是因为当温度上升时,射流与周围空气之间的换热速度加快,而激光功率却随着耦合长度的增加而降低,当在射流内生成的焦耳热小于对流换热量时,射流温度下降。随着射流直径的增加,射流的温度呈降低的趋势,当射流直径较大时,对高功率激光能量的传输过程中温度控制较为有利。图 5 1064 nm 波长激光在不同直径水射流的温度分布Fig.5 Temperature distribution of 1064 nm wavelength laser in different di-ameter water jet此外,为了对 1064 nm 波长激光的热效应进行对比分析,本文作者还对 532 nm 波长激光与水射流耦合的热效应进行了模拟,除了波长不同,其余参数保持不变,模拟结果如图 6 所示。随着射流直径的增加,射流温度分布呈现降低趋势,与 1064 nm 波长激光相比,图 6 532 nm 波长激光在不同直径水射流的温度分布Fig.6 Temperature distribution of 532 nm wavelength laser in different di-ameter water jet576激 光 技 术2023 年 9 月532 nm 波长激光与水射流耦合的热效应要小得多,这是由于 532 nm 波长激光比 1064 nm 波长激光在水中的吸收系数要低得多(532 nm 为 4.2810-4 cm-1,1064 nm 为 0.114 cm-1)15,只有很少的激光能量被水射流吸收转换为焦耳热,大部分的激光能量随射流传输至工件表面,正是因为这一特点,532 nm 波长的激光也常被应用到水导激光加工技术之中。但是,由于 532 nm 激光器价格昂贵、性能稳定性差、高功率不容易实现等特点,通常在 1064 nm 波长吸收率不好的材料才会考虑 532 nm 波长激光。为了探究连续激光和脉冲激光与水射流相互作用的区别,本文作者对平均功率 300 W、1000 W 的连续激光及不同脉冲宽度的脉冲激光温度分布进行求解,结果如图 7 所示。在相同的射流位置处,300 W 纳秒脉冲激光的温度比连续激光的温度低了约 5,1000 W 微秒脉冲激光的温度相较于连续激光低了约7,说明在水射流的循环冷却作用下,脉冲激光相较于连续激光耦合水射流后的热效应明显下降。并且可以看到,在平均功率、重复频率一致的前提下,随着激光脉宽的变化,水射流的温度未出现较大的变化,这说明脉宽并不是影响热效应的主要参数。图 7 连续激光与脉冲激光温度对比Fig.7 Temperature comparison between continuous laser and pulsed laser本文作者还对水导激光常用的参数,即射流直径为 50 m、100 m 以及耦合腔压力为 20 MPa、35 MPa时射流所能容忍的最大激光功率进行了求解,选取射流长度为 20 mm,认为当水射流的最高温度达到100 时的功率为最大容忍功率,此时射流破碎无法进行加工,结果如表 1 所示。表 1 水-光束耦合最大容忍功率Table 1 Maximum tolerance power of water-laser couplingpressure p0/MPajet diameter djet/mmaximum power P/W20501000201002500355012003510035004 结 论建立了水-光束耦合的数值模型,利用仿真软件模拟了不同入射激光功率等参数下的水射流温度分布,通过与实验数据对比,验证了模型的有效性。研究结果表明,激光功率、射流流速和射流直径是影响水射流温度的主要因素,入射激光功率的增大会导致射流温度的升高,入射激光功率为 1500 W 时,射流的最高温度达到了 67,温度的升高不利于射流的稳定。随着耦合腔压力的增大,射流速度加快,带走了更多的热量,可以有效降低水射流的温度。此外,当射流直径为50 m 时,射流的最高温度超过了 100,水射流的热效应严重,这造成了射流的不稳定性,甚至破碎,随着射流直径的增加,射流的温度呈降低的趋势。532 nm激光的热效应分析结果表明,在水导高功率激光技术中,采用 532 nm 激光是一个避免水导耦合高功率激光热效应的选择。入射激光功率为 300 W 时,脉冲激光与连续激光相比较,温度低了约 5,能较好地保证水-光束耦合的稳定性,并且脉宽对于热效应的影响几乎可以忽略。因此,为了保证水-光束耦合的热稳定性,提高激光能量在水射流中的传输效率,可以缩短水-光束耦合长度,在保证水射流稳定的前提下适当提高耦合腔的压力;同时,在满足加工效率和加工质量的前提下,尽量选择合适的激光功率和射流直径,削弱水-光束耦合的热效应,保证激光能量在水射流中的稳定传输。此外,本文中还求解了常用耦合腔压力及射流直径下射流的最大容忍功率,为水导高功率激光加工的功率选择提供了一定借鉴。参考文献1 KENNEDY P K,HAMMER D X,ROCKWELL B A.Laser-induced breakdown in aqueous mediaJ.Progress in Quantum Electronics,1997,21(3):155-248.2 NITIN S.The first coupling of a laser beam to a water jetJ.Pho-tonics Views,2021,18(1):72-76.3 GOBET M,OBI S,PAVIUS M,et al.Implementation of short-pulse lasers for wafer scribing and grooving applicationsJ.Journal of La-ser Micro Nanoengineering,2010,5(1):16-20.676第 47 卷 第 5 期赵 臻 水束导引高功率激光的热效应分析4 QIAO H Ch,CAO Zh H,CUI J F,et al.Experimental study on wa-ter jet guided laser micro-machining of mono-crystalline siliconJ.Optics and Laser Technology,2021,140(4):107057.5PAUCHARD,LEE K,VG N,et al.Advanced micromachining combining nanosecond lasers with water jet-guided laser technologyJ.Proceedings of the SPIE,2009,7201:72010A.6 ZHANG G Y,ZHANG Z,WANG Y F,et al.Gas shrinking laminar flow for robust high-power waterjet laser processing technologyJ.Optics Express,2019,27(26):38635-38644.7 TAN S N.Study on water guided laser processing technology of SiCp/Al compositesD.Harbin:Harbin Institute of Technology,2014:10-21(in Chinese).谭淞年.SiCp/Al 复合材料的水导激光加工技术研究D.哈尔滨:哈尔滨工业大学,2014:10-21.8 BOTYGINA N N,BUKATYI V I,LEVITSKII M E.The influence of thermal defocusing on the measurement of water transmissivity in a la-ser radiation fieldJ.Soviet Physics Journal,1978,21(8):1092-1094.9RICHERZHAGEN B.Interferometer for measuring the absolute re-fractive index of liquid water as a function of temperature at 1.064 mJ.Applied Optics,1996,35(10):1650-1653.10 XU J J,DAVIS S H.Instability of capillary jets with thermocapilla-rityJ.Journal of Fluid Mechanics,1985,161(1):1-25.11 MASHAYE F,ASHGRIZ N.Nonlinear instability of liquid jets with thermocapillarityJ.Journal of Fluid Mechanics,1995,283:97-123.12 CHRISTIAN B,HENNING J,MARKUS E,et al.Thermal investi-gation of interaction between high-power CW-laser radiation and a water-jetJ.Physics Procedia,2016,83:317-327.13 COUTY P,SPIEGEL,VG N,et al.Laser-induced break-up of water jet waveguideJ.Experiments in Fluids,2004,36:919-927.14 PINNICK R G,GILBERT L,FERNANDEZ J,et al.Stimulated Raman scattering and lasing in micrometer-sized cylindrical liquid jets time and spectral dependenceJ.Journal of the Optical Society of America,1992,B9(6):865-870.15SPIEGEL,VG N,WAGNER F R.High efficiency Raman scattering in micrometer-sized water jetsJ.Optical Engineering,2004,43(2):450-454.16MULLICK S,MADHUKAR Y K,KUMAR S,et al.Temperature and intensity dependence of Yb-fiber laser light absorption in waterJ.Applied Optics,2011,50(34):6319-6326.17 WANG Sh W,DING Y,CHENG B,et al.Research progress and mechanism of water-guided laser micromachiningJ.Chinese Jour-nal of Lasers,2022,49(10):1002404(in Chinese).王水旺,丁烨,程柏,等.水导激光微加工机理与研究进展J.中国激光,2022,49(10):1002404.18 YANG L F,JIAO H,HUANG Y X,et al.Simulation study of the flow field in the plane convergent nozzle based on the water guide la-serJ.Laser Technology,2020,44(6):754-761(in Chinese).杨林帆,焦辉,黄宇星,等.基于水导激光平面缩流喷嘴内流场仿真研究J.激光技术,2020,44(6):754-761.19 ZHANG G H,HUANG Y X,HUANG P,et al.Study on energy transmission law of water-laser coupling in water-jet guided laser technologyJ.Laser Technology,2022,46(6):749-754(in Ch-inese).张光辉,黄宇星,黄平,等.水导激光技术中水光耦合传能规律研究J.激光技术,2022,46(6):749-754.20 DENG C,YEO H R,KI H.Electrodynamic simulation of laser beam propagation in waterjet-guided laser processingJ.Optical Express,2020,28(8):11128-11143.776