Plasma pressure over time-space evolution law for femtosecond pulses laser shock peening
-
摘要: 为研究飞秒脉冲激光冲击强化中等离子体压力时空演化规律，利用考虑电子态密度（DOS）效应的模型计算了电子热容和电声耦合系数随电子温度的演化规律，并与采用QEOS（quotidian equation of state）模型计算结果进行了对比；提出DOS飞秒脉冲激光冲击强化模型，计算得到电子温度、晶格温度、等离子体羽位置时间演化规律和等离子体压力时空演化规律，并与QEOS飞秒脉冲激光冲击强化模型结果进行了对比。结果表明：DOS飞秒脉冲激光冲击强化模型计算得到的等离子体羽位置随时间的演化规律与实验结果吻合程度更好；增加激光能量或功率密度、考虑电子DOS效应会增加电子、晶格温度和等离子体压力。Abstract: The purpose of this research work is to look into the time-space evolution of plasma pressure for femtosecond pulse laser shock peening (fs-LSP). In this study, propose a model to understand plasma pressure over time-space process in fs-LSP based on the first principle, improved two temperature equations, and plasma hydrodynamic equations. Firstly analyze the plasma plume front location with respect to time by solving the plasma hydrodynamic equations. The simulated results by the electron DOS (density of state) femtosecond pulse laser shock peening model are in better agreement with the experiment results than the QEOS (quotidian equation of state) femtosecond pulse laser shock peening model. The DOS femtosecond pulse laser shock peening model was shown to be effective and superior. Then use the DOS model to calculate how the electron heat capacity and electron-phonon coefficient with respect to electron temperature. Electron heat capacity calculated by the QEOS model is larger than calculated by the electron DOS model, whereas the electron-phonon coefficient is the reverse. Moreover, the electron-phonon coefficient calculated by the QEOS model shows linear variation with respect to the electron temperature, which is the reverse of that calculated by the electron DOS model. Therefore, the electron DOS effect should be considered in two-temperature equations. Next see a graph of electron and lattice temperature with respect to time using the modified two-temperature equations to calculate. Increasing laser energy, decreasing pulse width, and considering the electron DOS effect will increase the electron’s peak temperature, and equilibrium temperature of electron and lattice systems, and reduce the electron-phonon relaxation time. Finally, we utilize the results of the two temperature equations as the initial condition to substitute into the plasma hydrodynamic equations to compute the plasma pressure. plasma peak pressure will rise as laser energy is increased, the pulse width is decreased, and the electron DOS effect is taken into account.
-
在众多的油气资源中，气态的油气资源，例如天然气、瓦斯和原油、成品油挥发气体等，属于典型的易燃易爆气体，稍有不慎就可能引发火灾爆炸安全事故，从而造成重大人员伤亡和巨额经济损失。近年来，油气火灾爆炸事故仍时有发生，造成大量人员伤亡和难以估量的经济损失。为了避免受限空间内油气爆炸事故的发生，减少人员伤亡与财产损失，尚需加强对油气爆炸抑制过程的研究，从而提出有效的油气爆炸抑制技术和措施。
可燃气体爆炸抑制过程是一个非常复杂的过程，从抑制手段上看可以分为主动抑制和被动抑制2类[1]。然而，无论采用何种抑制手段，其目的均是控制可燃气体爆炸过程中火焰的发展，从而使火焰在传播过程中逐渐减速衰减并熄灭。然而，从掌握的文献来看，涉及可燃气体爆炸抑制剂[2-3]和抑制技术[4-5]的文献较多，而涉及抑爆机理的文献较少，特别是涉及油气爆炸抑制过程火焰行为的文献十分少见。
氮气是一种非常稳定的惰性阻燃气体，广泛应用于石油化工、冶金等领域中有防火防爆要求的各种场合。因此，许多文献中对氮气的惰化特性[6]、稀释特性[7]以及防火防爆性能[8]开展了较为深入的探究。然而在这些文献中，氮气常常被用作预混抑制剂以惰化、稀释可燃气体，而将氮气作为非预混抑制剂来抑制受限空间内可燃气体爆炸还鲜有报道。
本文中拟在激波管上通过可视化实验手段，对油气爆炸氮气非预混抑制过程进行可视化实验研究，从而获得油气爆炸氮气非预混抑制过程的火焰行为，以期为油料洞库或坑道等涉油气受限空间内油气爆炸的氮气非预混抑制技术和装备的研究提供参考。
1. 实验装置与实验方法
实验的主要装置包括可组装式激波管(200 mm×200 mm×5500 mm)、油气循环系统、数据采集系统、点火系统、高速摄影仪、计算机等。数据采集系统主要包括压力采集系统和浓度测试系统等。图 1是氮气非预混抑制油气爆炸时的主要实验装置布置示意图。
为了直观地分析氮气非预混抑爆过程的火焰行为和抑爆机理，在激波管上安装了可视化段(280 mm×200 mm)，并利用高速摄影仪捕捉火焰遭遇氮气后的火焰行为，以揭示火焰在氮气中的行为规律和氮气非预混抑爆过程的机理。
实验中，用FASTCAM-Ultima 512型高速摄影仪捕捉火焰在氮气非预混段内的行为，根据高速摄影仪的快门速度与采样频率就可以获得火焰前锋速度[9-10]。具体来说就是根据相邻照片中火焰前锋的位置差和相邻照片的时间差来计算火焰速度。
油气循环系统包括真空循环泵、油气蒸发装置(具体结构及工作原理见文献[11-12]。压力采集系统主要由压力传感器、采集卡和计算机组成。沿激波管共布置了3个压力传感器，实验中以这3个压力传感器采集到数据的平均值作为激波管内的压力值。浓度采集系统主要由GXH-1050型红外分析仪和NHA-502型汽车尾气分析仪组成。实验中油气浓度由GXH-1050型红外分析仪测试，而其他气体组分，如O2、CO、CO2等，均由NHA-502型汽车尾气分析仪测量。
实验用的氮气由氮气制备系统提供，氮气纯度约为95%(其余气体组分为主要为氧气)。非预混氮气段用厚度为0.2 mm的塑料薄膜隔离形成，由于激波管可组装，因此非预混氮气段长度可调。实验时，先由氮气发生装置制备好所需氮气，之后设置好氮气抑爆段长度、位置并充入氮气，然后向激波管内充入一定浓度的油气，最后点火并采集实验数据。
为了使实验结果具有代表性，实验中油气当量比(按异辛烷计算)设置为1，经计算，此时的初始油气体积分数约为1.61%。本文实验的初始条件设置为：氮气段长度为226 cm，氮气段近端与点火源的距离为55 cm，氮气段氮气和氧气的体积分数分别为95%和5%，油气初始体积分数为1.61%，氧气的初始体积分数为20.6%，点火能为5 J，初始温度为300 K，初始相对压力为0 Pa。
2. 实验结果分析
2.1 氮气非预混抑制过程超压与主要组分体积分数变化
图 2是激波管内油气爆炸有无氮气非预混抑制的超压(p)曲线。由图 2可以看出，无氮气非预混抑制时油气爆炸最大超压为0.430 MPa，达到最大超压时的耗时为0.342 s，超压上升速率为1.26 MPa/s。而当有氮气非预混抑制时油气爆炸最大超压仅为0.126 MPa，达到最大超压时的耗时为0.210 s，超压上升速率约为0.60 MPa/s。有氮气非预混抑制时的最大超压比无氮气非预混抑制时下降了约70.7%，平均超压上升速率比无氮气非预混抑制时下降了约52.4%。
图 2 有无氮气非预混抑制时的超压曲线Figure 2. Curves of overpressure of the explosion with and without suppression by non-premixed nitrogen实验中，当氮气段上游的油气被点燃后，火焰面后的已燃高温高压气体会推动火焰面前的未燃气体向下游运动，此时氮气段上游的薄膜会先在压力作用下破裂，这会导致火焰前方的未燃气体和氮气段内的氮气混合。因此，火焰必然会进入氮气段传播。此时，如果火焰能持续穿过氮气段，则抑制实验失败，如果火焰未穿过氮气段，则抑制实验成功。抑制实验能否成功，与氮气非预混段长度、氮气段上游的油气段长度、氮气段内氮气的纯度、初始油气浓度等因素有关[13]。但是仅以油气爆炸最大超压和平均超压上升速率来说明油气爆炸氮气非预混抑制成功与否还不够充分，还必须提供其他证据。表 1是油气爆炸氮气非预混抑制实验前后激波管内主要气体组分的体积分数(φ)。实验后激波管内油气、O2、CO2和CO的体积分数分别为1.19%、15.24%、2.86%和0.89%。这表明，实验后激波管内仍有大量的油气和O2剩余，并且均处于油气爆炸所需的极限[14]内。这说明，本次实验中氮气非预混段后的油气混合物并未发生燃烧，即火焰并未穿过氮气非预混段，本次抑制实验是成功的。
表 1 实验前后激波管内主要气体组分的体积分数Table 1. Volume fraction of main gas components before and after the experiment测量条件 φ/% 油气 O2 CO2 CO 实验前 1.61 20.42 0.01 0 实验后 1.19 15.24 2.86 0.89 2.2 氮气非预混抑制过程火焰行为与火焰速度
图 3是高速摄影仪捕捉到的火焰前锋在氮气非预混段中熄灭的过程。其中高速摄影仪采样频率为1000 s-1，快门时间周期为0.001 s。
图 3 高速摄影仪捕捉到的火焰前锋在氮气中熄灭的过程Figure 3. Photos of flame front in extinction, captured by high speed camera由图 3可以看出：
(1) 油气爆炸火焰前锋进入氮气段后，在初期仍保持上游的发展趋势向前发展，此时火焰前锋仍由较为清晰的轮廓，火焰前锋的形状略呈向下游凸出的抛物线(228~232 ms)。这主要是由于火焰受上游膨胀压力波作用，从上游进入氮气段时，流场速度大于火焰的燃烧速度，因此火焰仍具有一定的“惯性”，持续发生燃烧反应并放出热量，这也导致火焰前锋边界有能量支撑，呈现较为清晰的轮廓。这一阶段火焰最主要的特征是惯性相持。
(2) 油气爆炸火焰前锋从234 ms开始由抛物线状逐渐演化为布满激波管横断面的矩形，火焰的颜色也从外围逐渐由亮黄色变成暗红色，火焰前锋的轮廓也变得模糊(234~244 ms)。这主要是由于从234 ms开始，氮气对油气爆炸火焰的抑制效果开始显现，氮气分子在火焰放热产生的对流机制作用下，开始由火焰外围向火焰内部扩散。由于氮气分子是惰性的，它在火焰高温下作为第3体[15]参与油气爆炸的化学反应，其主要反应过程可用以下方程式表示：
CH+2O+N2→OH+CO+N2 (1) OH+OH+N2→H2O+O+N2 (2) O+O+N2→O2+N2 (3) 式(1)~(3)中左边的氮气分子自身能量较低，右边氮气分子能量较高，它们参与反应主要是携带走CH、O、OH、CO等高能自由基的能量，在碰撞反应过程中高能自由基将内能传递给氮气分子，成为氮气分子的动能，氮气分子促使链式反应向中止链大量发展。正是由于这一机制的存在，使得链式反应的方向发生了变化，链式反应逐渐中止。相应的，火焰温度随着化学反应历程的改变逐渐降低，这导致碳烟颗粒激发态的辐射强度下降，相应的辐射光颜色也由亮黄色逐渐向暗红色过渡。火焰在这一阶段最主要的特征是抑制衰减。图 4显示了处于抑制衰减阶段的火焰结构。从图 4可以看出，此时的火焰由衰减抑制区和核心区火焰构成，火焰与氮气的相互作用主要发生在衰减抑制区内。
(3) 火焰前锋在246~276 ms之间持续衰减，颜色由黄红色逐渐变为暗红色，火焰区形状不能自持并逐渐扩散消失。这是由于此时火焰前锋在氮气的抑制作用下，化学反应速度持续衰减，化学反应放热已不能为火焰的继续传播提供能量，火焰区在流场作用下已不能维持自身形状，逐渐向四周扩散，这又导致热量损失进一步的加剧，火焰区温度迅速下降，火焰厚度迅速变薄，最终火焰完全在氮气中熄灭。火焰在这一阶段最主要的特征是扩散熄灭。
从上述分析可以看出，氮气非预混抑爆过程大致经历了3个阶段，即惯性相持期、抑制衰减期、扩散熄灭期。这3个阶段分别耗时约6、12和32 ms。
图 5是根据火焰行为的高速摄影照片得到的氮气非预混抑爆过程中火焰前锋速度(v)随时间的变化曲线。由图 5可以看出：(1)氮气非预混抑爆过程中火焰前锋的速度随时间波动，这主要是由于火焰在传播过程中，放热与热损失、混合气体的速度场、边界条件(如开口或闭口)以及油气爆炸压力波相互作用的共同作用导致了受限空间火焰速度的波动。这一现象与机理在文献中[11]已进行了详细的讨论。(2)油气爆炸火焰进入氮气抑制段时速度约为20 m/s，在上游惯性作用下火焰速度上升至43 m/s，之后氮气的抑制作用开始显现，火焰开始减速。在约238 ms时火焰速度降低至零，此时的火焰为驻留火焰，之后，火焰的减速并没有停止，火焰开始向后方运动，火焰速度出现了负值，在约240 ms时火焰速度出现最小值，约为-12.92 m/s。之后为扩散熄灭期，火焰在振荡流场影响下，速度又开始增加，在244 ms时火焰速度又上升至10.77 m/s。然而之后，油气爆炸火焰随流场开始扩散，火焰放热、火焰温度、火焰面积等迅速减小，直至火焰完全熄灭。(3)在抑制衰减期内，火焰速度的衰减可用以下线性拟合公式描述：
图 5 氮气非预混抑爆过程中火焰速度随时间的变化Figure 5. Variation of the flame speed with time in suppression processv=−0.5756t+1368 (4) 式中：火焰速度v的单位为m/s，时间t的单位为ms。拟合因子R2为0.983。可以看出火焰速度的衰减强度约为576 m/s，因此氮气对火焰的抑制效果非常显著。
3. 结论
本文中在激波管内完成了油气爆炸的氮气非预混抑制实验。通过可视化手段，得到了油气爆炸火焰在氮气段内衰减熄灭的照片。实验数据分析表明，氮气非预混抑制手段能显著降低油气爆炸过程的超压和平均超压上升速率；油气爆炸的氮气非预混抑爆过程经历了惯性相持阶段、抑制衰减阶段和扩散熄灭阶段3个阶段；惰性氮气分子在火焰高温下作为第3体参与化学反应并携带走高能自由基的能量，促使链式反应向中止链大量发展，这是油气爆炸氮气非预混抑制过程的主要机理；处于抑制衰减阶段的火焰结构由衰减抑制区和核心区火焰构成，火焰与氮气的相互作用主要发生在衰减抑制区内；在抑制衰减期内，火焰速度的衰减可用线性拟合公式描述。
-
图 1 飞秒脉冲激光冲击强化示意
Figure 1. Schematic illustration of femtosecond pulses laser shock peening
图 3 电子热容（Ce）和电声耦合系数（G）随电子温度演化规律
Figure 3. Electron heat capacity (Ce) and electron-phonon coefficient (G) versus electron temperature
-
[1] HU Y X, YAO Z Q, HU J. 3-D FEM simulation of laser shock processing [J]. Surface and Coatings Technology, 2006, 201(3/4): 1426–1435. DOI: 10.1016/j.surfcoat.2006.02.018. [2] 吴先前, 段祝平, 黄晨光, 等. 激光冲击强化过程中蒸气等离子体压力计算的耦合模型 [J]. 爆炸与冲击, 2012, 32(1): 1–7. DOI: 10.11883/1001-1455(2012)01-0001-07.WU X Q, DUAN Z P, HUANG C G, et al. A coupling model for computing plasma pressure induced by laser shock peening [J]. Explosion and Shock Waves, 2012, 32(1): 1–7. DOI: 10.11883/1001-1455(2012)01-0001-07. [3] DEVAUX D, FABBRO R, TOLLIER L, et al. Generation of shock waves by laser-induced plasma in confined geometry [J]. Journal of Applied Physics, 1993, 74(4): 2268–2273. DOI: 10.1063/1.354710. [4] WU B X, SHIN Y C. Laser pulse transmission through the water breakdown plasma in laser shock peening [J]. Applied Physics Letters, 2006, 88(4): 041116. DOI: 10.1063/1.2168022. [5] WU X Q, DUAN Z P, SONG H W, et al. Shock pressure induced by glass-confined laser shock peening: experiments, modeling and simulation [J]. Journal of Applied Physics, 2011, 110(5): 053112. DOI: 10.1063/1.3633266. [6] XIONG Q L, SHIMADA T, KITAMURA T, et al. Atomic investigation of effects of coating and confinement layer on laser shock peening [J]. Optics & Laser Technology, 2020, 131: 106409. DOI: 10.1016/j.optlastec.2020.106409. [7] RONDEPIERRE A, ÜNALDI S, ROUCHAUSSE Y, et al. Beam size dependency of a laser-induced plasma in confined regime: shortening of the plasma release. Influence on pressure and thermal loading [J]. Optics & Laser Technology, 2021, 135(1): 106689. DOI: 10.1016/j.optlastec.2020.106689. [8] WU B X, TAO S, LEI S T. Numerical modeling of laser shock peening with femtosecond laser pulses and comparisons to experiments [J]. Applied Surface Science, 2010, 256(13): 4376–4382. DOI: 10.1016/j.apsusc.2010.02.034. [9] NAKANO H, MIYAUTI S, BUTANI N, et al. Femtosecond laser peening of stainless steel [J]. Journal of Laser Micro/Nanoengineering, 2009, 4(1): 35–38. DOI: 10.2961/jlmn.2009.01.0007. [10] LEE D, KANNATEY-ASIBU JR E. Experimental investigation of laser shock peening using femtosecond laser pulses [J]. Journal of Laser Applications, 2011, 23(2): 022004. DOI: 10.2351/1.3573370. [11] AGEEV E I, BYCHENKOV V Y, IONIN A A, et al. Double-pulse femtosecond laser peening of aluminum alloy AA5038: effect of inter-pulse delay on transient optical plume emission and final surface micro-hardness [J]. Applied Physics Letters, 2016, 109(21): 211902. DOI: 10.1063/1.4968594. [12] SANO T, EIMURA T, KASHIWABARA R, et al. Femtosecond laser peening of 2024 aluminum alloy without a sacrificial overlay under atmospheric conditions [J]. Journal of Laser Applications, 2017, 29(1): 012005. DOI: 10.2351/1.4967013. [13] HOPPIUS J S, KUKREJA L M, KNYAZEVA M, et al. On femtosecond laser shock peening of stainless steel AISI 316 [J]. Applied Surface Science, 2018, 435: 1120–1124. DOI: 10.1016/j.apsusc.2017.11.145. [14] WANG H, PÖHL F, YAN K, et al. Effects of femtosecond laser shock peening in distilled water on the surface characterizations of NiTi shape memory alloy [J]. Applied Surface Science, 2019, 471: 869–877. DOI: 10.1016/j.apsusc.2018.12.087. [15] AGEEV E I, ANDREEVA Y M, IONIN A A, et al. Single-shot femtosecond laser processing of Al-alloy surface: an interplay between Mbar shock waves, enhanced microhardness, residual stresses, and chemical modification [J]. Optics & Laser Technology, 2020, 126: 106131. DOI: 10.1016/j.optlastec.2020.106131. [16] CHEN L, WANG Z S, GAO S, et al. Investigation on femtosecond laser shock peening of commercially pure copper without ablative layer and confinement layer in air [J]. Optics & Laser Technology, 2022, 153: 108207. DOI: 10.1016/j.optlastec.2022.108207. [17] TAN S, WU J J, ZHANG Y, et al. A model of ultra-short pulsed laser ablation of metal with considering plasma shielding and non-fourier effect [J]. Energies, 2018, 11(11): 3163. DOI: 10.3390/en11113163. [18] KIRAN KUMAR K, SAMUEL G L, SHUNMUGAM M S. Theoretical and experimental investigations of ultra-short pulse laser interaction on Ti6Al4V alloy [J]. Journal of Materials Processing Technology, 2019, 263: 266–275. DOI: 10.1016/j.jmatprotec.2018.08.028. [19] FAIRAND B P, CLAUER A H. Laser generation of high-amplitude stress waves in materials [J]. Journal of Applied Physics, 1979, 50(3): 1497–1502. DOI: 10.1063/1.326137. [20] FABBRO R, FOURNIER J, BALLARD P, et al. Physical study of laser-produced plasma in confined geometry [J]. Journal of Applied Physics, 1990, 68(2): 775–784. DOI: 10.1063/1.346783. [21] PEYRE P, SOLLIER A, CHAIEB I, et al. FEM simulation of residual stresses induced by laser peening [J]. The European Physical Journal Applied Physics, 2003, 23(2): 83–98. DOI: 10.1051/epjap:2003037. [22] WU B X, SHIN Y C. A self-closed thermal model for laser shock peening under the water confinement regime configuration and comparisons to experiments [J]. Journal of Applied Physics, 2005, 97(11): 113517. DOI: 10.1063/1.1915537. [23] LAVILLE S, VIDAL F, JOHNSTON T W, et al. Fluid modeling of the laser ablation depth as a function of the pulse duration for conductors [J]. Physical Review E, 2002, 66(6): 066415. DOI: 10.1103/PhysRevE.66.066415. [24] RETHFELD B, SOKOLOWSKI-TINTEN K, VON DER LINDE D, et al. Timescales in the response of materials to femtosecond laser excitation [J]. Applied Physics A, 2004, 79(4): 767–769. DOI: 10.1007/s00339-004-2805-9. [25] WU B X, SHIN Y C. A simple model for high fluence ultra-short pulsed laser metal ablation [J]. Applied Surface Science, 2007, 253(8): 4079–4084. DOI: 10.1016/j.apsusc.2006.09.007. [26] ALEXOPOULOU V E, MARKOPOULOS A P. A critical assessment regarding two-temperature models: an investigation of the different forms of two-temperature models, the various ultrashort pulsed laser models and computational methods [J]. Archives of Computational Methods in Engineering, 2023, 31(6). DOI: 10.1007/s11831-023-09974-1. [27] LIN Z B, ZHIGILEI L V, CELLI V. Electron-phonon coupling and electron heat capacity of metals under conditions of strong electron-phonon nonequilibrium [J]. Physical Review B, 2008, 77(7): 075133. DOI: 10.1103/PhysRevB.77.075133. [28] ZHANG Z Y, NIAN Q, DOUMANIDIS C C, et al. First-principles modeling of laser-matter interaction and plasma dynamics in nanosecond pulsed laser shock processing [J]. Journal of Applied Physics, 2018, 123(5): 054901. DOI: 10.1063/1.5021894. [29] WU B X, SHIN Y C. A one-dimensional hydrodynamic model for pressures induced near the coating-water interface during laser shock peening [J]. Journal of Applied Physics, 2007, 101(2): 023510. DOI: 10.1063/1.2426981. [30] ZHANG N, ZHU X N, YANG J J, et al. Time-resolved shadowgraphs of material ejection in intense femtosecond laser ablation of aluminum [J]. Physical Review Letters, 2007, 99(16): 167602. DOI: 10.1103/PhysRevLett.99.167602. 期刊类型引用(7)
1. 慕雪蒙，丛海勇，毕明树. 全氟己酮抑制乙醇/汽油替代物爆炸特性及机理研究. 火灾科学. 2023(04): 216-224 . 
百度学术2. 何昆，石英杰，李孝斌，田国兵. 惰性气体抑制丙酮蒸气爆炸实验对比研究. 消防科学与技术. 2021(01): 16-19 . 百度学术3. 胡文超，杜扬，梁建军，张培理，刘冲. 不同管径管道油气惰化抑爆实验研究. 当代化工. 2020(09): 1910-1914 . 百度学术4. 胡文超，杜扬，张培理. 油气惰化抑爆研究进展. 辽宁化工. 2019(08): 808-810 . 百度学术5. 郑立刚，李刚，王亚磊，朱小超，窦增果，杜德朋，余明高. 开口阻塞比对粉体抑制甲烷爆炸的影响研究. 爆炸与冲击. 2019(11): 144-153 . 
本站查看6. 姚丽蓉，赵海洋，赵德银，何金蓬. 塔一联含氮天然气对加热炉燃烧的影响. 石油与天然气化工. 2018(05): 51-57 . 百度学术7. 张培理，齐圣，王世茂，程顺国. 93号汽油蒸气绝热燃烧温度的计算分析. 后勤工程学院学报. 2017(02): 33-39 . 百度学术其他类型引用(5)
-
下载:
下载:
百度学术
计量
- 文章访问数: 385
- HTML全文浏览量: 128
- PDF下载量: 80
- 被引次数: 12