Pulsed laser heating – Fourier and electron kinetic theory approaches

Laser surface pulse heating of engineering metals is in demand in the metal industry and investigation into laser pulse heating becomes fruitful in this regard. Application of Fourier theory to heat conduction due to high power laser irradiation may give closed form solution to the problem. On the other hand, the heat flux through a given plane depends on the electron energy distribution through the material and at the scale of distance required to examine the problem, the material can no longer be considered as being homogeneous continuum, therefore, errors may occur when considering the Fourier theory in laser heating process. The problem requires to be examined in the quantum field. The present study examines the pulse laser heating process when considering both Fourier conduction and electron-kinetic theory approaches. Analytical solution to Fourier conduction equation is obtained for intensity exponential pulses while numerical scheme is introduced to solve the heat transfer equation resulted from kinetic theory approach. It is found that both Fourier and electron kinetic theory approaches result in similar temperature profiles for the pulses having the same energy content. In the case of electron kinetic theory approach the rise time for surface temperature to reach the melting point is shorter than that obtained from the analytical solution. Received on 23 February 1998     

Numerical approach to pulsed laser heating of semi-infinite aluminum substance

Processing of the reflective materials, such as aluminum, with a pulsed CO₂ laser beam depends largely on laser output power and pulse form. To enhance the understanding of the effect of pulse parameters on laser machining a modeling of laser induced heating is essential. The present study develops the heat transfer model allowing temporal variation of CO₂ laser output pulse, phase change process and temperature dependent thermal properties. A numerical technique is introduced to solve the resulting heat transfer equation. Aluminum is selected as workpiece and its surface reflectivity is taken into account in the computation. Thermal integration due to repetitive pulsing is also discussed. It is found that time corresponding to maximum temperature can be predicted by proper selection of pulsed parameters and the ability of the material to follow the laser pulse profile depends upon the pulse shape.     

Laser repetitive pulse heating influence of pulse duty on temperature rise

Laser repetitive pulse heating of the solid surface is considered and influence of duty cycle in the repetitive pulses on the temperature rise and temperature difference at the surface is examined. Laser power distribution with Gaussian profile is assumed and 3D heating model is accommodated in the heating analysis. A numerical method is employed to solve governing equations of heat transfer. It is found that the rate of surface temperature rise is high and temperature difference is large for high duty cycle during the consecutive pulse irradiation. This in turn enhances the thermal cyclic loading of the heated surface.     

非傅里叶热传导研究进展

傅里叶定律能够精确描述大多数的热传导问题,但对于超短脉冲激光加热等热作用的周期时间极短的超急速、超常规热传导等问题,非傅里叶效应将会显得至关重要．对非傅里叶热传导的实质、模型、模型的求解及应用与实验等几个方面的研究进展做了一个较详尽的概括与评述,并指出了今后需要着重研究的方向．      

电磁加热条件下皮肤组织的生物热力学行为

研究目的是开发一种数学方法来计算传热过程、热引起的力学响应以及相应的疼痛等级, 从而对临床上应用的各种加热疗法之间的差别进行定量评估. 采用基于有限差分法的数值模拟方法, 基于无限大和均匀化假设, 分析了各种热疗法中皮肤组织的温度、烧伤和热应力分布. 研究发现: 充血对热损伤的影响很小, 但对皮肤的温度分布影响很大, 而这又反过来显著影响由此产生的热应力场; 对于激光加热, 光波越短则峰值温度越高, 但峰值更接近皮肤表面温度; 激光和微波加热所产生的热应力集中于表皮顶层, 因为发热量沿皮肤深度方向呈指数衰减; 薄角质层(常常被忽略)对皮肤组织的热力学响应起主导作用.      

飞秒激光辐射纳米薄膜脱附的能量输运机理研究

从微观能量传输角度,对飞秒脉冲激光的辐射传热及其产生的金属薄膜吸附层分子脱附过程的能量输运机理进行了研究。以非平衡态传热和随机过程理论为基础,建立了描述吸附分子脱附过程的数学模型;并分别对单脉冲与双脉冲飞秒激光辐射下,金属Ru基体表面层的电子和声子的温度分布、吸附层分子CO的脱附概率分布进行了数值模拟。其结果与实验规律相符,为飞秒脉冲激光应用于微机电系统中杂质脱附提供了重要信息。     

短脉冲激光加热分数阶导热及其热应力研究

短脉冲激光加热引起材料内部复杂的传热过程及热变形,现有的以Fourier定律或Cattaneo-Vernotte松弛方程结合弹性理论为框架建立起来热应力理论在刻画其热物理过程存在严重缺陷. 本文基于分数阶微积分理论, 以半空间为研究对象, 建立了分数阶Cattaneo热传导方程和相应的热应力方程, 给出了问题的初始条件和边界条件, 采用拉普拉斯变换方法, 给出了非高斯时间分布激光热源辐射下温度场和热应力场的解析解, 研究了短脉冲激光加热的温度场及热应力场的热物理行为. 数值计算中, 首先对理论解进行数值验证, 然后取分数阶变量$p=0.5$研究温度场和热应力场的变化特点及激光参数对温度和热应力的影响,最后数值计算分数阶参数对温度和热应力场的影响. 计算结果表明, 分数阶Cattaneo传热方程和热应力方程描述的温度和热应力任然具有波动特性,与经典的Fourier传热模型和标准的Cattaneo传热模型相比, 分数阶阶次越大, 热波波速越小, 热波波动性越明显; 反之, 则热波波速越大, 热扩散性越强.激光加热和冷却的速度越快, 温度上升和下降的速度越快, 压应力和拉应力交替变化越快, 温度变化幅值越小, 热应力幅值影响不明显.      

Applications of Crank-Nicolson method with ADI in laser transformation hardening

A two-dimensional numerical solution for pulsed laser transformation hardening is developed using the finite difference method (FDM). The FDM has been developed using Crank-Nicolson scheme which solved by using alternating-direction implicit method. If this present model was compared to the analytical solution, then the Crank-Nicolson scheme showed better results in terms of accuracy, consistency, stability, convergence, and performance than to the explicit scheme. The longer heating duration, higher laser beam intensity, and greater number of pulse had influences on increasing the maximum temperature. The repetitive heating had influences on extending the heat duration and increasing the initial temperature of domain. The shorter cooling duration in repetitive pulse produced higher maximum temperature. The thinner material’s thickness increased the cooling rate, which finally increased the possibility of austenite to transform into martensite phase. In addition, it was also found that the higher maximum temperature always reduced the cooling rate value when temperature cools down toward to the starting temperature of martensite formation. It reduced the possibility of martensite formation. It was also seen that the heat was conducted more effective to the axial direction than to the radial direction.     

热-力联合作用下柱壳结构变形的计算分析

本文研究了激光连续加热和恒内压作用下柱壳结构的变形规律。以4340钢材料作为研究对象，其本构方程选用热粘塑性本构模型（Johnson-Cook模型）；选取某一应变率为临界值（本文中临界应变率取为1s^-1)。考察了激光功率密度，预载荷大小，激光作用时间等对结构变形规律的影响。主要结论有：结构失稳是导致激光辐照下充内压柱壳破坏的重要原因；激光功率密度和内压越高，结构破坏时间越短；根据辐照时间的长短，预测到了结构急速破坏，延迟破坏和不破坏三种模式。本工作对于深入认识激光作用下预载结构的热－力联合破坏有一定意义。     

Analysis of nano channel formation in quartz cubes by laser-induced process

A novel laser processing technique was developed for making channels in the nano regime in this paper. A Nd:YAG laser was used to dry fabricate micro channels (25μm～100μm di ameter) in a 1 cm^3 fused silica substrate by thermal-induced processing. By controlling the locations of these initiating micro channels on a silica cube, 1D-controllable self-connecting nano fractures can be formed as rectangular channels. These nano channels are smooth and with extremely high aspect ratio (～10^4 depth to width ratio). A possible mechanism is proposed to explain the formation of the nano channels. This laser-based nano channel fabrication technique is fast and inexpensive, and with potential applications in capillary electrophoresis and electro-osmosis driven nano-filtration.     

Transient behavior of boiling bubbles generated on the small heater of a thermal ink jet printhead

The fundamental behavior of boiling bubbles generated on a small film heater used for thermal ink jet (TIJ) printers is investigated experimentally under the condition of a single pulse heating in a pool of water. The pulse power and the pulse width are varied in wide ranges that include the printing conditions. As the pulse power is increased or the pulse width increased at a fixed high pulse power, numerous fine bubbles appear simultaneously on the heater and then coalesce into a thin vapor film to grow to a vapor bubble, before collapsing at the center of the heater. For a long pulse width sequence, the coalesced bubble repeats the growth and collapse. Bubble behavior is also studied in the same heat flux range using a platinum film heater enabling surface temperature measurement. From a comparison of the two heaters, the dominant mechanism of nucleation on the TIJ heater is believed to be spontaneous nucleation at around the heating rate for printing. The dependence of the size and lifetime of the coalesced bubble on pulse power and pulse width are examined. Based on the analytical model presented by Asai [J. Heat Transfer 113 (1991) 973], the pressure impulse arising during the rapid evaporation of the superheated liquid, presumed to dominate the subsequent growth of the coalesced bubble, is estimated from the measured size of the coalesced bubble. The relationship between the pressure impulse and the superheat energy in the liquid is discussed.     

Thermal analysis of thin multi-layer metal films during femtosecond laser heating

Multi-layer metals films are widely used in modern engineering applications such as gold-coated metal mirrors used in high power laser systems. A transient heat flux model is derived to analyze multi-layer metal films under laser heating. The two separate system composed of electrons and the lattice is considered to take into account the electron–lattice interaction. The present model predicted the effects of underlying chromium’s thermal properties on temperature rise of the top gold layer. The effects of two adjacent and different metals with different electron–lattice coupling factors are analyzed for the heating mechanism of different lattices. The derived transient model combined with the two different conservation equations for the lattice and electrons are applied for the ultra short-pulse laser heating of a multi-layer film composed of gold and chromium.     

Analytical analysis of the dual-phase-lag model of bioheat transfer equation during transient heating of skin tissue

Dual-phase-lag model of bioheat transfer equation is utilized in treating the transient heat transfer problems in skin tissue considering prevalent heating conditions in thermal therapy applications, namely, pulse train and periodic heat flux. Comparisons between the presented analytical results for limiting cases and previous studies display an excellent agreement. The effects of temperature gradient relaxation time on the tissue temperature, damage, and also on the blood perfusion in skin tissue are studied.     

A GENERALIZED THERMOELASTIC DIFFUSION PROBLEM OF THIN PLATE HEATED BY THE ULTRASHORT LASER PULSES 1)

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势, 近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域. 本文基于L-S型广义热弹扩散理论, 建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论, 它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应. 推导了所建理论的控制方程, 并基于拉普拉斯积分变换获得了控制方程的解. 作为算例, 利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题, 得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律. 结果表明: 传热对传质影响显著, 传质对传热影响甚微; 非局部参数对位移、应力影响显著, 对温度、化学势和浓度几乎没有影响. 该理论及求解方法的建立, 旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.     

Heat transfer analysis of a semi-infinite solid heated by a laser beam

Laser heating of semi-infinite solid substance is introduced and analytical solution appropriate to the laser machining power intensities is obtained. Both conduction limited and non-conduction limited cases are considered. The limiting cases in non-conduction limited equation are investigated and the equation governing the conduction limited heating is deduced. In the case of non-conduction limited heating process, the location of maximum temperature and its functional relation with the substance and laser properties are developed. Steady state analysis of the heating process is taken into account and limiting laser power intensity is concluded. In addition, the applicability of Fourier heat conduction equation in laser heating is discussed.     

超短激光脉冲加热薄板的广义热弹扩散问题

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势, 近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域. 本文基于L-S型广义热弹扩散理论, 建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论, 它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应. 推导了所建理论的控制方程, 并基于拉普拉斯积分变换获得了控制方程的解. 作为算例, 利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题, 得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律. 结果表明: 传热对传质影响显著, 传质对传热影响甚微; 非局部参数对位移、应力影响显著, 对温度、化学势和浓度几乎没有影响. 该理论及求解方法的建立, 旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.      

基于时域间断有限元方法的生物组织非傅立叶热行为数值分析

在热传导分析中,当热流与温度梯度存在时间延迟时,需采用非傅立叶热传导模型进行分析。生物组织具有较强的热松弛时间系数,承受激光、微波及烧烫等作用时,其呈现出较强的非傅立叶行为。本文对脉冲热源作用下生物组织的非傅立叶热传导进行研究,针对强脉冲引起的温度场在空间域的高梯度变化、波阵面的间断行为以及通用传统时域数值方法会带来虚假数值振荡的特点,提出采用所发展的时域间断Galerkin有限元法（DG-FEM ）进行求解计算。对多种脉冲热源作用下的非傅立叶热传导过程进行数值模拟,通过考量强脉冲作用下温度场分布和热致生物组织损伤行为的影响,表明了本文所发展的DGFEM 能够有效、准确地描述温度场空间分布和热传导过程以及非傅立叶行为下的生物热损伤更为明显,在生物组织热行为分析中应该受到重视。     

Analytical solution for temperature field in thin film initially heated by a short-pulse laser source

Analytical solution for electron and lattice temperature distribution in the solid initially heated by a laser short-pulse is presented. Strained parameters method is introduced when formulating electron and lattice temperature distributions. Laser short pulse heating of gold film is simulated numerically and temperature data at the end of the heating pulse are adopted as initial condition to the governing equations of energy transport for analytical solutions. This enables to solve the governing equations of energy analytically in the cooling period. It is found that electron temperature decays sharply while lattice site temperature increases gradually in the surface regions during the cooling cycle. As the depth from the surface increases change in both temperatures become gradual.     

Influence of laser pulse shape both on temperature profile and hardened layer depth

This paper presents a mathematical modelling and numerical calculations of heat conduction problems where laser generated heat is assumed as a surface heat source. Also the effect of a laser time structure on a hardened layer depth is examined. Temperature profiles for different laser pulse shapes are determined from the solution of a linear one-dimensional heat conduction equation for semi-infinite medium and discussed in terms of the parameters evolution such as dimensionless: temperature, heat flux, hardening depth, laser impulse duration and increasing time of triangular pulse shape.     

Pulse laser ablation of ground glass

We have studied pulse laser ablation of ground glass surface. Ablation process in air and in vacuum has been observed by high speed framing camera. Burst of small fragments of glass has been observed in the present experiment, when ground glass surface is laser ablated through glass plate from rear side. Production of macro particles by laser ablation is an inherent characteristic of ground glass, and no similar phenomena have been observed in case of ablation of other transparent materials. By using ns-duration Nd:YAG laser of 100–400 mJ/pulse, observed maximum particle velocity ranges from 0.5 to 1.5 km/s. In order to understand the particle generation process, scanning electron micrograph observation of the ground surface in the region of no laser irradiation and also in the damaged region, has been made. Cleavage surface structure has been evidenced in the damaged surface area, which stems from plenty of micro cracks covering the virgin glass surface. Effects of surface roughness on the particle generation were studied by using ground glass of quartz with different surface roughness. Produced glass particles were captured in vacuum by a polymethylmetacrylate (PMMA) plate. Size distribution was obtained by analysing the trapped particles on the PMMA plate and revealed that most probable particles size has almost no dependence on the initial surface roughness. Finally, we applied the phenomena to ignite pentaerithritoltetranitrate (PETN) powder explosive, and succeeded in igniting (PETN) powder only by laser ablation of ground glass.