Laser-ablation-induced refractive index fields studied using pulsed digital holographic interferometry

Pulsed digital holographic interferometry has been used to investigate the plume and the shock wave generated in the ablation process of a Q-switched Nd-YAG (λ=1064 nm and pulse duration=12 ns) laser pulse on a polycrystalline boron nitride (PCBN) target under atmospheric air pressure. A special setup based on two synchronised wavelengths from the same laser for simultaneous processing and measurement has been used. Digital holograms were recorded for different time delays using collimated laser light (λ=532 nm) passed through the volume along the target. Numerical data of the integrated refractive index field were calculated and presented as phase maps showing the propagation of the shock wave and the plume generated by the process. Radon inversion has been used to estimate the 3D refractive index fields measured from the projections assuming rotational symmetry. The shock wave density has been calculated using the point explosion model and the shock wave condition equation and its behaviour with time at different power densities ranging from 1.4 to 9.1 GW/cm² is presented. Shock front densities have been calculated from the reconstructed refractive index fields using the Gladstone–Dale equation. A comparison of the shock front density calculated from the reconstructed data and that calculated using the point explosion model at different time delays has been done. The comparison shows quite good agreement between the model and the experimental data. Finally the reconstructed refractive index field has been used to estimate the electron number density distribution within the laser-induced plasma. The electron number density behaviour with distance from the target at different power densities and its behaviour with time are shown. The electron number densities are found to be in the order of 10¹⁸ cm^?3 and decay at a rate of 3×10¹⁵ electrons/cm³ ns.     

大气环境下飞秒激光对铝靶烧蚀过程的研究

利用时间分辨的光阴影成像技术研究了在大气环境下飞秒激光烧蚀铝靶的动态过程. 在入射激光能量为4 mJ, 激光光斑超过1 mm时, 激光烧蚀区表面物质以近似平面冲击波形式向外喷射; 在同样激光能量下、激光光斑较小时(约0.6 mm), 激光烧蚀区以近似半球型冲击波形式向外喷射. 当激光能量比较大时(7 mJ), 发现空气的电离对于激光烧蚀靶材有着重要影响. 在光轴附近烧蚀产生的喷射物具有额外的柱状和半圆型的结构, 叠加在平面冲击波结构上.     

Focusing of shock waves induced by optical breakdown in water

The focusing of laser-generated shock waves by a truncated ellipsoidal reflector was experimentally and numerically investigated. Pressure waveform and distribution around the first (F(1)) and second foci (F(2)) of the ellipsoidal reflector were measured. A neodymium doped yttrium aluminum garnet laser of 1046 nm wavelength and 5 ns pulse duration was used to create an optical breakdown at F(1), which generates a spherically diverging shock wave with a peak pressure of 2.1-5.9 MPa at 1.1 mm stand-off distance and a pulse width at half maximum of 36-65 ns. Upon reflection, a converging shock wave is produced which, upon arriving at F(2), has a leading compressive wave with a peak pressure of 26 MPa and a zero-crossing pulse duration of 0.1 mus, followed by a trailing tensile wave of -3.3 MPa peak pressure and 0.2 mus pulse duration. The -6 dB beam size of the focused shock wave field is 1.6 x 0.2 mm(2) along and transverse to the shock wave propagation direction. Formation of elongated plasmas at high laser energy levels limits the increase in the peak pressure at F(2). General features in the waveform profile of the converging shock wave are in qualitative agreement with numerical simulations based on the Hamilton model.     

连续波激光辐照下金属材料的二维非线性温度场研究

 建立了连续波激光作用下金属材料动态温升的物理模型和有限元计算模型。并分别计算了在平均功率密度1 kW/cm²、光斑直径28 mm的连续波激光辐照下，30CrMnSiA钢及LD10-CS铝合金前、后表面及内部温度分布。得到了给定条件下，材料动态升温过程，靶前表面达到熔点的时刻，以及不同时刻靶前、后表面的温度场分布。并比较了两种材料不同的温度响应特性和不同反射率对温升的影响。     

Comparison of the laser ablation process on Zn and Ti using pulsed digital holographic interferometry

Pulsed digital holographic interferometry has been used to compare the laser ablation process of a Q-switched Nd-YAG laser pulse (wavelength 1064 nm, pulse duration 12 ns) on two different metals (Zn and Ti) under atmospheric air pressure. Digital holograms were recorded for different time delays using collimated laser light (532 nm) passed through the volume along the target. Numerical data of the integrated refractive index field were calculated and presented as phase maps. Intensity maps were calculated from the recorded digital holograms and are used to calculate the attenuation of the probing laser beam by the ablated plume. The different structures of the plume, namely streaks normal to the surface for Zn in contrast to absorbing regions for Ti, indicates that different mechanisms of laser ablation could happen for different metals for the same laser settings and surrounding gas. At a laser fluence of 5 J/cm², phase explosion appears to be the ablation mechanism in case of Zn, while for Ti normal vaporization seems to be the dominant mechanism.     

激光引致激波演化数值研究

在高重复频率激光推进的研究中,激波的合并是发生在激波演化后期的,同时由于脉冲间隔短,脉冲宽度对流场演化的影响也需要详细研究。考虑了激光辐照过程对流场演化的影响,通过数值计算对激波演化特性进行了研究。结果表明,初期波阵面的椭球形状逐渐转化为一个球形,球心与击穿点的距离随时间逐渐减小并最终趋于稳定。基于激波合并的应用,当激波马赫数在1~2之间时,给出了激波波阵面半径随时间的变化规律,以及激波高压区长度和波峰压强随激波波阵面半径变化规律的经验公式。     

Interaction of Shock Waves in Cement Mortar Plate Investigated by the Digital Speckle Correlation Method

Interaction of shock waves in cement mortar plate is studied by digital speckle correlation method and digital high-speed photography technique. When the plates were destroyed by two detonators exploding at the same time, variation of shock wave field is obtained. Experimental results show that the interaction of shock waves will result in a nonlinear huge increase of local normal strain, leading to large deformation and serious destruction. However, the occurrence of this strongly nonlinear phenomenon sensitively depends on the interval between detonators, and it will only appear when the interval is smaller than the diameter of the region where shock waves exist.     

飞秒激光诱导水光学击穿阈值

为了研究飞秒激光诱导水光学击穿阈值随激光脉冲参数的变化关系,采用四阶RungeKutta方法对飞秒激光诱导水光学击穿的椭球体模型进行了不同脉宽(40～540fs)、波长(400～1200nm)和光斑尺寸(0～200μm)下的数值模拟。通过控制变量法得出阈值光强与这些激光脉冲参数的关系曲线图,据此定性分析了阈值光强与激光脉冲参数的变化特征趋势。应用光强与功率、能量、辐照曝光量和电场强度之间的关系,得到了它们随激光脉冲参数脉宽、波长和光斑尺寸的动态关系,这为进一步研究飞秒激光与水和含水介质的相互作用提供了理论依据。     

Laser forming of a bowl shaped surface with a stationary laser beam

Despite a lot of research done in the field of laser forming, generation of a symmetric bowl shaped surface by this process is still a challenge mainly because only a portion of the sheet is momentarily deformed in this process, unlike conventional sheet metal forming like deep drawing where the entire blank undergoes forming simultaneously reducing asymmetry to a minimum. The motion of laser beam also makes the process asymmetric. To counter these limitations this work proposes a new approach for laser forming of a bowl shaped surface by irradiating the centre of a flat circular blank with a stationary laser beam. With high power lasers, power density sufficient for laser forming, can be availed at reasonably large spot sizes. This advantage is exploited in this technique. Effects of duration of laser irradiation and beam spot diameter on the amount of bending and asymmetry in the formed surface were investigated. Laser power was kept constant while varying irradiation time. While varying laser spot diameter laser power was chosen so as to keep the surface temperature nearly constant at just below melting. Experimental conditions promoted almost uniform heating through sheet thickness. The amount of bending increased with irradiation time and spot diameter. It was interesting to observe that blanks bent towards the laser beam for smaller laser beam diameters and the reverse happened for larger spot diameters (~10 times of the sheet thickness). Effect of spot diameter variation has been explained with the help of coupled thermal-structural finite element simulations.     

Acoustic characteristics of pulse detonation engine with ellipsoidal reflector

Acoustic characteristics of a pulse detonation engine(PDE) with and without an ellipsoidal reflector are numerically and experimentally investigated. A two-dimensional(2 D) non-splitting unstructured triangular mesh Euler solver based on the space-time conservation element and solution element(CE/SE) method is employed to simulate the flow field of a PDE.The numerical results clearly demonstrate the external flow field of the PDE. The effect of an ellipsoidal reflector on the flow field characteristic near the PDE exit is investigated. The formation process of reflected shock wave and reflected jet shock are reported in detail. An acoustic measurement system is established for the PDE acoustic testing. The experimental results show that the ellipsoidal reflector changes the sound waveform and directivity of PDE sound. The reflected shock wave and reflected jet shock result in two more positive pressure peaks in the sound waveform. The ellipsoidal reflector changes the directivity of PDE sound from 20 to 0. It is found that the peak sound pressure level(PSPL) and overall sound pressure level(OASPL) each obtain an increment when the PDE is installed with a reflector. The maximum relative increase ratio of PSPL and OASPL are obtained at the focus point F2, whose values are 6.1% and 6.84% respectively. The results of the duration of the PDE sound indicate that the reflecting and focusing wave generated by the reflector result in the increment of A duration and B duration before and near focus point F2. Results show that the ellipsoidal reflector has a great influence on the acoustic characteristic of PDE sound. The research is helpful for understanding the influence of an ellipsoidal reflector on the formation and propagation process of PDE sound.     

Laser driven shock wave studies in gold coated Plexiglas targets

Abstract

Laser-driven shock wave propagation in a transparent material such as Plexiglas coated with a thin overlayer of gold is studied using the technique of high speed optical shadowgraphy. A Nd: glass laser was focussed to produce intensities in the range of 10′²-10′⁴W/cm² on the target, within an irradiation spot diameter of 160 pm, optical shadowgrams were recorded bya second harmonic (0.53 pm wavelength) pulse. Shock pressures and scaling of pressure with laser intensity was studied. Shock pressures in gold-coated Plexiglas target was observed to be considerably higher compared to those in uncoated targets. This enhancement of shock pressure has been explained on the basis of contribution of an X-ray driven ablative heat wave in the gold plasma. Shock pressure values show a close agreement with those obtained from a one-dimensional Langrangian hydrodynamic simulation. Shadowgrams of shock fronts produced by non-uniform spatial laser beam irradiation profiles have shown complete smoothing when a gold layer is used on a Plexiglas target.     

Modeling of an electrohydraulic lithotripter with the KZK equation.

The acoustic pressure field of an electrohydraulic extracorporeal shock wave lithotripter is modeled with a nonlinear parabolic wave equation (the KZK equation). The model accounts for diffraction, nonlinearity, and thermoviscous absorption. A numerical algorithm for solving the KZK equation in the time domain is used to model sound propagation from the mouth of the ellipsoidal reflector of the lithotripter. Propagation within the reflector is modeled with geometrical acoustics. It is shown that nonlinear distortion within the ellipsoidal reflector can play an important role for certain parameters. Calculated waveforms are compared with waveforms measured in a clinical lithotripter and good agreement is found. It is shown that the spatial location of the maximum negative pressure occurs pre-focally which suggests that the strongest cavitation activity will also be in front of the focus. Propagation of shock waves from a lithotripter with a pressure release reflector is considered and because of nonlinear propagation the focal waveform is not the inverse of the rigid reflector. Results from propagation through tissue are presented; waveforms are similar to those predicted in water except that the higher absorption in the tissue decreases the peak amplitude and lengthens the rise time of the shock.     

飞秒脉冲烧蚀单晶硅的超快动力学

利用基于Pump-probe系统的超快时间分辨阴影图的方法,研究了空气中飞秒激光烧蚀单晶硅的动力学过程。实验采用脉宽为50 fs、平均能量密度约35 J/cm2的单脉冲激光烧蚀单晶硅,获取飞秒激光烧蚀单晶硅过程中等离子体和冲击波的形成和发展过程的时间分辨阴影图。实验结果表明:飞秒激光烧蚀单晶硅导致其表面物质喷发的过程是不连续的,分为明显的两次喷发过程。这表明飞秒激光与单晶硅作用的过程中,在不同的时间段可能由不同的机制主导,在前期可能是多光子电离为主,在后期可能是由多光子效应和雪崩效应共同作用。研究还发现,延迟时间较长时,冲击波形状发生畸变。     

针对高压物质密度诊断的激光X射线照相优化设计

激光驱动冲击波能提供高压的物质状态, 是状态方程研究的重要工具. 超短超强激光与固体靶相互作用产生的X射线源, 具有短脉冲、微焦点、高产额、能量可调的特点, 是高压物质密度测量的首选背光源. 本文基于蒙特卡罗程序Geant4建立了X射线照相模型, 客体密度分布由流体力学程序Multi-1D模拟激光冲击加载高压物质获得. 在一维长方体形密度客体情况下, 定义了均方根、峰值偏差与上升沿斜率比三种指标, 对照相图像求解的密度结果进行评价, 开展了照相结果信噪比、分辨率与对比度等参数优化. 并开展了一维圆柱形密度客体的照相模拟, 建立了基于Radon逆变换法的Abel反演算法. 反演结果与模拟设计密度分布符合良好, 要求X射线源半径在5 μm以内; 反演结果与模拟设计密度分布较为一致, 要求X射线源半径在15 μm以内.     

激光聚焦扰动作用下高超声速边界层稳定性实验研究

在马赫数6、单位雷诺数3.1×10⁶/m的条件下对半锥角7°直圆锥边界层稳定性开展了实验研究.以激光聚焦于流场中局部空间而产生的膨胀冲击波作为人工添加的小扰动,分析了该扰动对高超声速圆锥边界层流动稳定性的影响.实验中利用响应频率达到兆赫兹量级的高频压力传感器对圆锥壁面脉动压力进行测量,通过对压力数据进行短时傅里叶分析和功率谱分析发现,相比于不添加激光聚焦扰动的结果,添加激光聚焦扰动使边界层中第二模态波的出现位置提前,且扰动波的幅值大幅度地增加,在相同的流向范围内,激光聚焦扰动将边界层中的扰动波从线性发展阶段推进到非线性发展阶段,其对边界层中扰动波发展的促进效果明显.同时,激光聚焦位置的不同对边界层中扰动波的发展也具有不同的影响.当激光直接聚焦于圆锥壁面X=100 mm位置时,边界层中频率为90 kHz的扰动波幅值增长最快,在X=500 mm的位置处其幅值放大倍数为3.81,相比而言当激光聚焦位置位于圆锥前方自由来流中时,边界层幅值增长最快的扰动波频率大幅减小为73 kHz,相同范围内,其幅值放大倍数为4.51倍.由此可见,当激光聚焦位置位于圆锥上游的自由来流中时,其对边界层中扰动波的影响更为显著.     

激光诱导2024铝合金表面动态应变检测

采用PVDF贴片传感器对脉冲激光作用下2024铝合金表面的动态应变进行了测量,分析了动态应变曲线的特性。结果表明,PVDF贴片传感器在动态应变测量中动态响应快,灵敏度高,可有效应用于脉冲激光诱导材料表面动态应变的实时测量。脉冲激光作用过程中,2024铝合金冲击光斑周围材料先受挤压,后压应变减小。脉冲激光作用结束后,2024铝合金冲击光斑周围材料表面粒子在卸载稀疏波和表面稀疏波的作用下不断往复运动,冲击光斑周围材料甚至受到了拉应变的作用。最后随着时间的推移,材料表面粒子的动态响应经反复震荡后逐渐衰弱形成最终的稳定状态。     

纳秒激光烧蚀铝材料的二维数值模拟

为了探索纳秒脉冲强激光与材料的相互作用机理,建立了二维数值模型,利用有限差分法对纳秒激光脉冲烧蚀金属铝的温度场进行了数值模拟.通过对比不同脉宽、光斑和能量下激光引起的温度场随时间的演化,发现脉冲的前期温度升高比后期快.等温图显示中心温度升高最快,烧蚀轮廓与激光束形状相似,烧蚀深度达1—5 μm.脉宽越长,烧蚀越窄和越深,光斑越大,烧蚀越宽和越浅.数值研究表明,1)激光的脉冲形状、脉宽和功率密度直接影响烧蚀的形状和深度,2)激光功率密度在10⁹ W/cm²量级烧蚀     

Acoustic emission monitoring during laser shock cleaning of silicon wafers

A laser shock cleaning is a new dry cleaning methodology for the effective removal of submicron sized particles from solid surfaces. This technique uses a plasma shock wave produced by laser-induced air breakdown, which has applied to remove nano-scale silica particles from silicon wafer surfaces in this work. In order to characterize the laser shock cleaning process, acoustic waves generated during the shock process are measured in real time by a wide-band microphone and analyzed in the change of process parameters such as laser power density and gas species. It was found that the acoustic intensity is closely correlated with the shock wave intensity. From acoustic analysis, it is seen that acoustic intensity became stronger as incident laser power density increased. In addition, Ar gas has been found to be more effective to enhance the acoustic intensity, which allows higher cleaning performance compared with air or N₂ gas.     

Distribution characteristic of scattering field for an ellipsoidal target irradiated by an electromagnetic wave from an arbitrary direction

It is of great importance for engineering applications to obtain the expression of scattering field for an ellipsoidal target irradiated by an electromagnetic wave from an arbitrary direction. Literature relevant to this problem is seldom found. In this paper, the scattering field for an ellipsoidal target is presented by utilizing the scale transformation of electromagnetic field and the rotation of coordinate system, with an electromagnetic wave projecting on the target from an arbitrary direction. The obtained result is in good agreement with the solution available from the literature if we consider the scale factors to be unity. Taking a conducting ellipsoidal target for sample, we perform the partial simulations of the ellipsoidal model and a plant leaf model by choosing different scale factors. The obtained results show that the distribution characteristic of scattering field is sensitively affected by the polarization of the incident wave and varies not much with the incident wave angle but changes with the observation point. At some points the scattering energy arrives at its maximum.     

激光冲击强化提高K403铸造高温合金疲劳性能

采用YLSS-M60U型高能Nd:YAG激光器,对发动机高压涡轮叶片材料K403/K3铸造高温合金试片进行激光冲击强化处理,强化工艺参数为:激光能量3J,光斑直径2.6mm,脉宽20ns,波长1064nm,吸收保护层为铝箔,约束层为水,搭接率50%,冲击3次。强化后,在420MPa应力水平下进行了室温高周振动疲劳测试,并进行了扫描电镜观察和X射线衍射仪物相分析。研究结果表明:激光冲击强化后,试片疲劳寿命是原始状态试片寿命的2.4倍,激光冲击强化的强冲击波作用使金属发生高应变率塑性变形,以及随之产生的较大较深残余压应力,是金属疲劳性能提高的主要原因。