Numerical simulation of phase transition process of laser drilling on GH4037 nickel-based superalloy
-
摘要: 为了更准确地研究激光打孔相变过程,基于流体传热和流体力学理论,建立了GH4037镍基高温合金激光打孔相变模型。模型中考虑了重力、粘滞力、反冲压力的作用,以及材料的固-液相变和液-气相变过程,通过数值计算得到了激光打孔相变过程的温度场和速度场。结果表明,气化材料的反冲压力可以加快熔池的流动,在激光功率为2000W、脉宽为1.70ms时,材料最大气化蒸发速率可以达到250m/s。该模型为进一步开展激光打孔研究提供了理论基础。
-
关键词:
- 激光技术 /
- 激光打孔 /
- 数值模拟 /
- GH4037镍基高温合金
Abstract: In order to study the phase transformation process of laser drilling more accurately, based on the theories of fluid heat transfer and fluid mechanics, a phase transition model of laser drilling on GH4037 nickel-based superalloy was established. In the model, the effect of gravity, viscous force, recoil pressure, solid-liquid phase transition and liquid-gas phase transition process were considered, and the temperature field and velocity field of laser drilling phase transition process were obtained by numerical calculation. The results show that, the reaction pressure of the gasification material can accelerate the flow of the molten pool. When the laser power is 2000W and the pulse width is 1.70ms, the maximum evaporation rate of the material can reach 250m/s. This model provides a theoretical basis for further research on laser drilling. -
引言
电磁超材料是由用来模拟电/磁偶极子的亚波长结构单元排列而成的等效介质。它的介电常数和磁导率[1-3]可以通过设计结构单元及其排列方式来调控。利用超材料,人们不但可以实现对电磁波传播的有效操控而且可以模拟一些有趣的物理效应,比如电磁诱导透明效应[4]、克尔效应[5-6]、负折射效应[7-8]等,同时还可以用来设计一些特殊的微波器件,如超棱镜[9-10]、隐身斗篷[11-12]等。因此,电磁超材料能够用于设计民用或军用上的微波器件[13]。
在2008年,LANDY等人首次提出了电磁超材料完美吸波器的设计概念[14],随后,TAO[15]等人基于电谐振设计了一种工作在THz波段的超材料吸波器。由此电磁吸波器便成为了研究热点,在理论和实验研究上都取得了许多令人瞩目的成果。目前有关吸波器的研究工作主要集中在多频带吸波[16-18]、宽带吸波[19-21]、极化不敏感[22]、大角度入射[23-24]、超薄结构设计[25]等方面。然而,常规吸波器最大的缺陷是很难得到可调谐的吸收频谱。
等离子体超材料是超材料的一种,它一般由等离子体介质周期结构构成,并可以通过调节等离子体的频率(密度)或激励不同谐振区域来改变其物理特性。用等离子体材料设计吸波器可以得到可调谐的吸收频谱。关于等离子体超材料吸波器的研究尚处于起步阶段,因此,研究等离子体超材料吸波器的设计方法和理论具有重要的意义。
作者设计了一种能工作在微波频段的基于等离子体超材料的吸波器,通过改变等离子体的激励区域来实现对吸收频谱的动态调控。在理论上,用全波仿真法计算得到的吸收频率研究了该吸收器的电磁特性,并探讨了该吸波器吸收效应的物理机理。
1. 理论模型
图 1是作者设计的等离子体材料吸波器的单元结构示意图。图 1a是正视图,图 1b是侧视图。由图 1可知,该结构单元的底层为金属铜板(电导率为σ=5.8×107S/m),中间层为介质基板,上层由被截断的等离子体环形结构和三角形金属铜片构成。介质基板为FR-4(相对介电常数εr=4.3,损耗角正切tanδ=0.025),介质基板厚度h=2.9mm,介质基板的边长l=72mm,上层贴片和底层金属铜板的厚度w=0.04mm,外层等离子体谐振环(ring 1)和内层等离子体谐振环(ring 2)的间距a=1mm,外层环的边长q=50.821mm,表示为3l/5+22x/5,其中x=√3a,内层环的边长r=43.2mm,表示为3l/5,中间三角贴片的边长s=18.6mm,表示为3l10−√3x,外层环的上下谐振单元间距p=14.314mm,表示为24x/5+2c,其中c=0.2mm,内层环的上下谐振单元间距e=2mm,外层环上面谐振单元的宽度b=1.8mm,内层环上下谐振单元的宽度t=1.2mm,外层环下面的“V”型结构中左边谐振单元的宽度d=1.2mm,与内层环下面的“V”型谐振单元的间距g=1mm,外层环下面的“V”型结构中右边谐振单元的宽度u=1.2mm,与内层环下面的“V”型谐振单元的间距v=0.9mm,上层等离子体的介电常数用Drude模型来描述:
εp(ω)=1−ω2pω2+jωωc (1) 式中, ω表示角频率,等离子频率ωp=2.9×1015rad/s, 碰撞频率ωc=1.65×10141/s。电磁波波矢方向为沿着-z方向垂直入射,本文中提及的TE波为电场平行于y轴,磁场平行x轴。吸波器的吸收率A(ω)可以表示为:
A(ω)=1−R(ω)−T(ω) (2) 式中, R(ω)为反射率,T(ω)为透射率。由于本文中设计的吸波器底层为金属铜板,所以透射率T(ω)=0,那么吸收率表示为:
A(ω)=1−R(ω) (3) 由(3)式可知,如果R(ω)越小,那么A(ω)的值就越大。
2. 结果分析与讨论
图 2中给出了激励不同等离子体谐振单元时的吸收频谱。如图 2a所示(ring 1中的等离子体谐振结构被激励,而ring 2中的等离子体谐振结构未被激励),吸收频谱中存在着两个吸收率较高的频点,它们的值分别为89.69%和90.04%,位于9.2GHz和11.71GHz。显然,此时该吸波器的吸波效果较差。如果ring 1和ring 2中等离子体谐振结构同时被激励,其吸收频谱如图 2b所示。由图 2b所示,在8GHz~12GHz的频率范围内存在着两个吸收率大于90%的吸收频域,它们分别为9.16GHz~9.36GHz和11.16GHz~11.38GHz。在9.19GHz, 9.31GHz和11.31GHz处的吸收率分别达到99.03%, 95.92%和99.21%。比较图 2a和图 2b中的结果可知,ring 1和ring 2中等离子体谐振结构同时被激励时,该吸波器的吸波特性得到了明显地改善。显然,人为地改变上层等离子体谐振结构的激励形态,不仅可以改善该吸波器的吸收特性,还能实现对吸收频域的可调谐。换句话说,可以以最为简单的单极化多频吸波器为起点,实现极化不敏感的宽频吸波器的设计。从这个角度出发,采用最简单的电谐振吸波器的设计思路(单极化),通过等离子体超材料的可调谐性实现了由频谱中的单频吸收与窄带吸收间的相互切换(极化敏感),也为下一步设计极化不敏感的宽频吸波器奠定了基础。
为了研究该吸波器工作的物理本质,进一步地理解电磁波损耗的机理,对该吸波器在f1=9.19GHz和f2=9.31GHz两个频点处的表面电流进行了计算。图 3是该吸波器在这个两点频点上的上表面和背面的电流分布图。由图 3可知,当垂直入射的TE波的频率分别为f1=9.19GHz和f2=9.31GHz时,在该吸波器的上下表面都会形成表面电流。对于上表面而言,其表面电流主要集中在上表面的两个“V”形等离子体谐振结构的底部,其电流方向如图中黑色箭头所示。由图 3还可知,在这两个吸收频点该吸波器的背面同样也会产生表面电流,其表面电流也主要集中在正对上表面“V”形谐振结构底部处,其表面电流的方向与上表面的相同。此时,该吸波器的上下表面可以看成一个电偶极子,将在上下层间形成电谐振。这种电谐振将会和外部入射的电磁波进行耦合,使得入射电磁波的能量在介质基板中被损耗,从而实现了该吸波器在这两个频点吸波。因此,该吸波器主要是电谐振的形式来实现对入射电磁波的吸收。
图 4中给出了入射角θ(入射电磁波矢与y-O-z面的夹角)与吸收频谱的关系图。由图 4a可知,当入射角θ分别为0°, 30°和60°时,该吸波器的吸收频谱几乎没有发生变化,在9.16GHz~9.36GHz和11.16GHz~11.38GHz这两个频域范围内其吸收率都可以达到90%以上,其吸收率峰值分别位于9.32GHz, 9.39GHz和11.3GHz,它们的值分别为98.5%, 98.8%和99.6%。显然,该吸波器的角度稳定性较好,随着入射角θ的增大,该吸波器的最大吸收率会有所增加。为了进一步说明这个问题,图 4b中给出了入射角θ与吸收率的关系图。由图 4b可知,当入射角θ由0°连续变化到80°时,该吸波器的吸波特性几乎保存不变。随着入射角θ的增大,吸收频域上下边缘会略微地向高频方向移动。当θ=80°时,改变吸波器的工作频域(吸收率大于90%)变为9.27GHz~9.45GHz和11.16GHz~11.38GHz。因此,由图 4中的结果可知,该吸波器具有较好的角度稳定性。
为了进一步研究该吸波器的特性,在图 5中给出了结构参量c和v对吸收频谱的影响。图 5a给出了其它参量不变的情下,结构参量c为0.2mm,0.4mm,0.6mm和0.8mm时的吸收频谱。由图 5a可知,c=0.2mm时,在频率范围8.7GHz~9.9GHz内该吸波器工作频域为9.16GHz~9.36GHz,且有两个吸收峰值分别为99%和94.9%位于f1=9.19GHz和f2=9.31GHz。随着c的增大,该吸波器的工作频域将逐渐减小。当c为0.4mm,0.6mm和0.8mm时,该吸波器的工作频域分别为9.16GHz~9.23GHz,9.16GHz~9.21GHz和9.16GHz~9.20GHz。但f1=9.19GHz处的第一吸收峰值将几乎保持不变,而第二吸收峰将随着c的增大而减小,并向高频方向移动。为了进一步说明这个问题,图 5b中给出了吸收频谱与结构参量c的关系图。由图 5b可知,f1=9.19GHz处的第一吸收峰值几乎不会随着c的增大而变化,而该吸波器工作频域将会随着c的增大先增大后减小,最后再保持不变。c=0.2mm时,该吸波器的工作带宽最宽,即0.20GHz。当c≥0.6mm时,该吸波器的工作带宽最窄,且几乎不会随着c的增大而发生变化,即0.04GHz,位于9.16GHz~9.20GHz。显然,改变参量c的值不仅能够对该吸波器的工作频域进行调谐,还能对第二吸收峰的大小和位置进行调谐。同样,图 5c中给出了其它参量不变的情况下, 结构参量v为0.9mm, 1.1mm, 1.3mm和1.5mm时的吸收频谱。由图 5c可知,v=0.9mm时,该吸波器工作频域为9.16GHz~9.36GHz,且两个吸收峰值分别位于f1=9.19GHz和f2=9.31GHz。随着参量v的增大,该吸波器的工作带宽也将逐渐地减少。但第一吸收峰的频率将保持f1=9.19GHz不变,而第二吸收峰的频率f2将随着参量v的增大而向低频方向移动。当v为1.1mm, 1.3mm和1.5mm时,该吸波器的工作频域分别为9.16GHz~9.34GHz, 9.16GHz~9.279GHz和9.16GHz~9.277GHz。图 5d中给出了吸收频谱与结构参量v的关系图。由图 5d可知,该吸波器工作频域将会随着v的增大先增大后减小。v=0.9mm时,该吸波器的工作带宽最大且等于0.20GHz。v=1.7mm时,该吸波器的工作带宽最小且等于0.13GHz,位于9.14GHz~9.27GHz。显然,改变参量v的值对该吸波器的工作带宽有明显的调谐作用。综上所述,只要人为地改变等离子体谐振结构的激励区域(如改变结构参量c和v的大小)就可以得到可调谐的吸收频谱,通过相应的参量优化就能拓展该吸波器的工作频域。
3. 结论
设计了一款基于等离子体超材料的吸波器,能够实现在TE波下得到可调谐的吸收频谱。采用了全波仿真方法对该吸波器的吸收频谱和表面电流图进行了计算,并探讨了结构参量c, v和入射角度θ对吸收率的影响。研究结果表明:通过激励不同的等离子体谐振结构不但可以改善该吸波器的吸收特性,而且还能获得可调谐的吸收频谱。该吸波器的工作频谱位于9.16GHz~9.36GHz和11.16GHz~11.38GHz,其吸波机理主要是通过电谐振的方式使得入射电磁波的能量在介质基板中得以损耗。当入射角θ的值发生改变时,该吸波器的工作频域几乎不发生变化,即该吸波器具有较好的角度稳定性。改变参量c和v的大小都能够实现对吸收频谱的调谐。增加参量c的值,该吸波器的工作频域将先增大后减少最后趋于一个定值(9.16GHz~9.20GHz),而增加参量v的值,该吸波器的工作频域将先增大后减少,其最小值为9.14GHz~9.27GHz(v=1.7mm)。吸收频谱中的第一峰的位置(f1=9.19GHz)将不会随着c和v的变化而变化,而第二吸收峰的位置f2将随着c的增大而向高频方向移动,而随着v的增大而向低频方向移动。显然,只要人为地调控激励等离子体谐振单元的区域,不但能够都到可调谐的吸收频谱, 还能拓展吸波器的工作带宽。该研究也为设计新型波器吸提供了思路。
-
Table 1 Physical properties of air
T/K cp/(J·kg-1·K-1) κ/(W·m·-1K-1) μ /(Pa·s) 300 1.005×103 0.0262 1.983×10-5 500 1.029×103 0.0403 2.671×10-5 1000 1.141×103 0.0675 4.152×10-5 1500 1.230×103 0.0946 5.400×10-5 2000 1.338×103 0.1240 6.500×10-5 2500 1.688×103 0.1750 7.670×10-5 Table 2 Physical properties of GH4037 Nickel-based superalloy
propertiy value specific heat of solid phase cp, s 440J/(kg·K) thermal conductivity of solid phase κs 13.8W/(m·K) thermal conductivity of melting phase κm 23.9W/(m·K) dynamic viscosity μ 0.006Pa·s latent heat of melting Hm 296kJ/kg latent heat of vaporization Hv 6423kJ/kg melting temperature Tm 1585K vaporization temperature Tv 3005K -
[1] WANG Zh P, LIU Y, WANG Y. Experiment study on laser cutting and process simulation[J]. Manufacturing Technology & Machine Tool, 2016(6):118-122(in Chinese).
[2] FU B Y, OUYANG B Sh, LIU W D, et al. Study on fiber laser rotary drilling process of 0.12mm SUS304[J]. Optical and Technique, 2016, 42(2):126-129(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/gxjs201602007
[3] QIAN X Zh, WANG Q Q, REN N F. Optimization of laser drilling processing parameters for SUS304 based on orthogonal experiments[J]. Laser Technology, 2017, 41(4):578-581(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jgjs201704024
[4] HAN F M, XU Sh Zh, SONG W L, et al. Study of nanosecond laser ablation on aluminum and stainless steel targets[J].Chinese Journal of Lasers, 2016, 43(2):203005(in Chinese). DOI: 10.3788/CJL
[5] HANON M M, AKMAN E, OZTOPRAK B G, et al. Experimental and theoretical investigation of the drilling of alumina ceramic using Nd:YAG pulsed laser[J]. Optics & Laser Technology, 2012, 44(4):913-922. http://www.sciencedirect.com/science/article/pii/S0030399211003446
[6] NATH A K, HANSDAH D, ROY S, et al. A study on laser drilling of thin steel sheet in air and underwater[J]. Journal of Applied Physics, 2010, 107(12):123103. DOI: 10.1063/1.3447866
[7] LOW D K Y, LI L, BYRD P J. The effects of process parameters on spatter deposition in laser percussion drilling[J]. Optics & Laser Technology, 2000, 32(5):347-354. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=5190fa0016263b302e17a48f8fa6e4cc
[8] LOU D Y, XIONG H, WU Y G, et al. Damage threshold and drilling mechanism of aluminum plate by nanosecond laser with different pulse widths[J]. Laser Technology, 2017, 41(3):427-432(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jgjs201703024
[9] XU B Zh, QI L J, WANG W, ZHU X, et al. Study on the optimum matching parameters of the combined laser drilling[J]. Laser Technology, 2018, 42(1):5-10(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jgjs201801002
[10] SWIFT-HOOK D T, GICK A E F. Penetration welding with laser[J]. Welding Journal, 1973, 52:492-499. DOI: 10.1063-1.334571/
[11] WU Y M, LI J Ch. Simulation of laser drilling based on semi-analytic solution[J]. Laser Technology, 2009, 33(2):205-208(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jgjs200902024
[12] XU B Q, WANG H, XU G D, et al. Numerical modeling of laser-induced molten pool for laser interaction with metal material[J]. Journal of Jiangsu University (Natural Science Edition), 2010, 31(3):358-362(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=jslgdxxb201003024
[13] ZHANG L, QIN Y, LU J, et al. Analysis of energy efficiency ratio for long-pulse laser drilling based on the heat conduction equation[J]. Journal of Optoelectronics·Laser, 2015, 26(10):1990-1996(in Chinese).
[14] CHU Q C, YU G, LU G Q, et al. Two-dimensional numerical investigation for the effects of laser process parameters on hole type during laser drilling[J]. Chinese Journal of Lasers, 2011, 38(6):0603001(in Chinese). DOI: 10.3788/CJL
[15] BEGIC-HAJDAREVIC D, BIJELONJA I. Experimental and numerical investigation of temperature distribution and hole geometry during laser drilling process[J]. Procedia Engineering, 2015, 100:384-393. DOI: 10.1016/j.proeng.2015.01.382
[16] SONG L L, SHI G Q, LI Zh G. Simulation on laser drilling temperature field by using ANSYS[J]. Journal of Changchun University of Science and Technology, 2006, 29(4):19-21(in Chinese). http://www.cnki.com.cn/Article/CJFDTotal-CGJM200604005.htm
[17] SOWDARI D, MAJUMDAR P. Finite element analysis of laser irradiated metal heating and melting processes[J]. Optics & Laser Technology, 2010, 42(6):855-865. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=e42595431b3082f04bedfc162ab3a22a
[18] JIE Q, ZHANG G H, CHEN Y W. Real fluid dynamics[M]. Beijing:Tsinghua University Press, 1986:383-385(in Chinese).
[19] THE COUNCIL OF EDITORS OF CHINA AERONAUTICAL MATERIALS HANDBOOK. China aeronautical materials handbook[M]. Beijing:Standards Press of China, 2002:260-261(in Chinese).
[20] ZHANG Y, SHEN Z, NI X. Modeling and simulation on long pulse laser drilling processing[J]. International Journal of Heat & Mass Transfer, 2014, 73(4):429-437. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dcea0d2ef8345b96fa4352481cbbbd6d
-
期刊类型引用(1)
1. 张浩,马宇,章海锋,杨靖,刘佳轩. 一种带宽展宽的等离子体超材料吸波体的设计. 激光技术. 2019(02): 256-262 . 本站查看
其他类型引用(2)