Laser short pulse heating: Influence of pulse intensity on temperature and stress fields

Laser short-pulse heating of gold surface is considered and the influence of laser pulse intensity on the temperature and stress fields is investigated. Laser step input pulses with different pulse lengths and the same energy content are employed in the simulations. The electron kinetic theory approach employing thermomechanical coupling is introduced to model the non-equilibrium energy transport in the electron and lattice sub-systems. Thermal stress development in the lattice sub-system and temperature rise in the lattice and electron sub-systems are computed. It is found that electron temperature rises rapidly while lattice site temperature rise is gradual in the early heating period, which is more pronounced for high intensity pulses. Thermal stress component in the axial direction is compressive and its magnitude is considerably less than the yielding limit of the substrate material.     

Laser short-pulse heating of metallic surface: Consideration of Seebeck effect

In the present study, laser short-pulse heating is formulated using an electron kinetic theory approach. Temperature predictions are compared with that obtained from two-equation model. Seebeck effect is considered during the heating process. The predicted Seebeck coefficients are compared with the results based on the early formulation. Electron excess energy loss due to Seebeck effect is compared with electron mean energy. It is found that Seebeck coefficient decays sharply in the surface region due to sharp decay of electron temperature in this region. Seebeck coefficient obtained from the present study is in agreement with the predictions based on the early formulation, provided electron temperature is used in the previous formulation. However, Seebeck coefficient differs significantly once the lattice site temperature is used in the previous formulation. Electron excess energy loss due to Seebeck effect is considerably less than electron mean energy, i.e. the ratio is in the order of 10⁻⁵.     

Effect of temperature oscillation on thermal characteristics of an aluminum thin film

Energy transport in aluminum thin film is examined due to temperature disturbance at the film edge. Thermal separation of electron and lattice systems is considered in the analysis, and temperature variation in each sub-system is formulated. The transient analysis of frequency-dependent and frequency-independent phonon radiative transport incorporating electron–phonon coupling is carried out in the thin film. The dispersion relations of aluminum are used in the frequency-dependent analysis. Temperature at one edge of the film is oscillated at various frequencies, and temporal response of phonon intensity distribution in the film is predicted numerically using the discrete ordinate method. To assess the phonon transport characteristics, equivalent equilibrium temperature is introduced. It is found that equivalent equilibrium temperature in the electron and lattice sub-systems oscillates due to temperature oscillation at the film edge. The amplitude of temperature oscillation reduces as the distance along the film thickness increases toward the low-temperature edge of the film. Equivalent equilibrium temperature attains lower values for the frequency-dependent solution of the phonon transport equation than that corresponding to frequency-independent solution.     

Laser cutting of rectangular geometry into alumina tiles

Laser cutting of rectangular geometry into the 5 mm thick alumina tiles is carried out. Temperature and stress fields, which are developed during the cutting process, are simulated in line with the experimental conditions. The morphological changes in the cutting sections are examined using optical and electron scanning microscopes, energy dispersive spectroscopy, and X-ray diffraction technique. The predictions of surface temperature and the residual stress are validated through the experimental data. It is found that von Mises stress attains high values in the region of the mid-thickness of the workpiece. The laser cut sections are free from major cracks and large scale sideways burning. The predictions of surface temperature and residual stress agree well with their counterparts obtained from the experiment.     

Elastic displacement of surface due to laser picosecond pulse heating of gold

Laser picosecond pulse heating initiates nonequilibrium energy transport in the surface vicinity of the metallic substrates. In this case, electron temperature rises rapidly while lattice site temperature is slow during the heating period. Although the rise of lattice site temperature is low, the temperature gradient is high. This results in elastic displacement of the surface. To model the heating process thermomechanical coupling needs to be introduced in the analysis. In the present study, picosecond pulse heating of gold surface is considered. Three-dimensional analysis of electron kinetic theory is introduced when formulating the electron and the lattice site temperatures. The analysis is extended to include the thermomechanical coupling due to mechanical response of the substrate material. The lattice site temperatures obtained from the electron kinetic theory are compared with the predictions of the two-equation model. It is found that both models predict almost identical temperature profiles in the surface vicinity of the substrate material. In addition, the surface displacement on the order of 10⁻¹¹ m is predicted.     

Laser short-pulse heating of silicon film with the presence of metallic substrate

Laser interaction of silicon film located at he top of metallic substrate is examined and energy transport in electron and lattice sub-systems are formulated using the electron kinetic theory approach. The simulations are repeated for different substrate materials, namely gold, silver, and copper. It is found that electron temperature in the silicon film rises in the vicinity of the silicon–metallic substrate interface, despite the fact that energy absorption from the irradiated filed is significantly low in the silicon film. Lattice site temperature rises rapidly in the early heating period at the interface. In addition, lattice site temperature increase is higher in the silicon film than that corresponding to the metallic substrate.     

Exact solution for temerature field due to non-equilibrium heating of solid substrate

Non-equilibrium energy transfer between electron and lattice sub-systems due to short-pulse heating is formulated and the closed form solution for electron and lattice site temperatures is presented. The electron kinetic theory approach is incorporated to formulate non-equilibrium energy transfer in the electron and lattice sub-systems. The method of Lie point symmetries is used in the exact solution of governing energy equation. In the analysis, the volumetric heat source, representing the laser heating pulse, and surface heat source, corresponding to short thermal contact of the surface, are incorporated and the analytical solutions for each heating source are presented. Electron temperature distribution obtained from the closed form solution is compared with its counterpart predicted from the numerical simulation. It is found that the results obtained from the closed form agree well with electron temperature predictions obtained from numerical simulation.     

Cemented carbide cutting tool: Laser processing and thermal stress analysis

Laser treatment of cemented carbide tool surface consisting of W, C, TiC, TaC is examined and thermal stress developed due to temperature gradients in the laser treated region is predicted numerically. Temperature rise in the substrate material is computed numerically using the Fourier heating model. Experiment is carried out to treat the tool surfaces using a CO₂ laser while SEM, XRD and EDS are carried out for morphological and structural characterization of the treated surface. Laser parameters were selected include the laser output power, duty cycle, assisting gas pressure, scanning speed, and nominal focus setting of the focusing lens. It is found that temperature gradient attains significantly high values below the surface particularly for titanium and tantalum carbides, which in turn, results in high thermal stress generation in this region. SEM examination of laser treated surface and its cross section reveals that crack initiation below the surface occurs and crack extends over the depth of the laser treated region.     

Laser bending of AISI 304 steel sheets: Thermal stress analysis

Laser induced bending of steel sheet is carried out and thermal stress developed in the heated region is examined. Temperature and stress fields are predicted using the finite element model. The microstructural changes in the melted region are investigated through scanning electron microscope, energy dispersive spectroscopy and X-ray diffraction. The residual stress developed at the surface vicinity of the laser treated region is measured using the X-ray diffraction technique, which is then compared with its counterpart predicted from the simulations. It is found that the residual stress at the surface vicinity is compressive and the prediction of the residual stress agrees well with that obtained from the X-ray diffraction technique. In addition, surface temperature predictions are in good agreement with the thermocouple data. The laser treated region is free from major cracks and large cavities.     

Plasma kinetic control in a tokamak

Plasma kinetic (temperature and density) control is developed in terms of optimal control theory. This involved solving a distributed parameter control problem for which no previous general theory or techniques could be practically implemented. The methodology used involves reducing the one-dimensional particle and energy transport partial differential equations into a system of ordinary differential equations. The latter are linearized and put into standard control theory format. This technique can be generalized to any distributed parameter system whose variables can be modeled as simple analytic functions of the spatial coordinates. This distributed parameter control technique shows excellent control characteristics when applied to realistic plasma temperature and density profiles. Tests were made on temperatures and densities which had been perturbed about 10% below their desired value and profiles that were significantly more peaked than the required shape. Results were obtained for a simplified model problem with specific empirical transport coefficients and a one species plasma     

稠密气固流动二阶矩摩擦应力模型和模拟

基于统一二阶矩两相湍流应力模型和稠密颗粒动力学理论,建立了欧拉-欧拉双流体二阶矩颗粒摩擦应力模型。模型充分反映各向异性的气固两相相间雷诺应力相互作用,并引入有效颗粒弹性恢复系数,考虑了因颗粒表面不光滑产生的摩擦力对湍流流动结构和颗粒弥散特性的影响,对于下降管的计算结果与实验吻合较好。因颗粒摩擦产生能量耗散降低了颗粒温度和导致颗粒雷诺应力再分布。在入口和出口区域因颗粒碰撞频率较高而产生的能量耗散对流动结构影响明显。     

Numerical Investigation Into the Highly Nonlinear Heat Transfer Equation with Bremsstrahlung Emission in the Inertial Confinement Fusion Plasmas

A highly nonlinear parabolic partial differential equation that models the electron heat transfer process in laser inertial fusion has been solved numerically. The strong temperature dependence of the electron thermal conductivity and heat loss term (Bremsstrahlung emission) makes this a highly nonlinear process. In this case, an efficient numerical method is developed for the energy transport mechanism from the region of energy deposition into the ablation surface by a combination of the Crank‐Nicolson scheme and the Newton‐Raphson method. The quantitative behavior of the electron temperature and the comparison between analytic and numerical solutions are also investigated. For more clarification, the accuracy and conservation of energy in the computations are tested. The numerical results can be used to evaluate the nonlinear electron heat conduction, considering the released energy of the laser pulse at the Deuterium‐Tritium (DT) targets and preheating by heat conduction ahead of a compression shock in the inertial confinement fusion (ICF) approach. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of electron kinetic,two-temperature,and one-temperature models for laser short-pulse heating

The laser short-pulse heating of metallic workpieces initiates the non-equilibrium heating in the surface vicinity of the substrate. The material response to the non-equilibrium heating cannot be predicted accurately by the one-temperature model. Consequently, new models pertinent to laser short-pulse heating are needed. In the present study, laser short-pulse heating of gold, copper, and lead is considered. The material responses to the laser short-pulse due to the electron kinetic theory and the two-temperature and the one-temperature models are examined in detail. The differences between the collisional and diffusional heating mechanisms are presented. The conditions for the convergence of conduction mechanisms are discussed. The electron kinetic theory, the two-temperature, and the one-temperature predictions are compared for three substrates. It is found that the electron kinetic theory predictions differ from the predictions of the one-temperature model in the surface vicinity of the substrate during the early heating duration. As the heating progresses, both models predict similar temperature profiles. The electron kinetic theory and the two-temperature model predictions are in good agreement. PACS 44.10.+i; 42.62.-b     

Determination of pyrolysis temperature for charring materials

An energy and mass balanced method of determining the pyrolysis temperature is proposed. The concept is to find the pyrolysis temperature that consumes the same amount of energy to produce the same amount of mass when using the pyrolysis front model as when using finite rate kinetics models for the entire charring process. The resulting pyrolysis temperature has the form of pyrolysis rate weighted average temperature. Comparisons between finite rate kinetics and pyrolysis front models for various boundary conditions, geometries, heats of decomposition, kinetic parameters and assumptions used in the literature were made to assess the proposed method. Models using energy and mass balanced pyrolysis temperature show good agreement with finite rate models and the experiments. Extensive numerical studies on various factors influencing the charring material pyrolysis show that heat flux, sample size, heat of decomposition and kinetic parameters are the most important factors for determining an appropriate pyrolysis temperature. Thermal conductivity, specific heat and density have a lesser effect on the pyrolysis temperature. For practical application, a non-dimensional correlation is developed to determine the appropriate pyrolysis temperature without solving the problem by using finite rate models. With this correlation the energy and mass balanced pyrolysis temperature can be determined with a standard deviation of 7.6 K. These predictions are validated by comparison with measurements of wood cylinder pyrolysis. A good agreement suggests that simpler pyrolysis front models yield practically useful and accurate results given an appropriate pyrolysis temperature.     

Fluid turbulence in quantum plasmas

Nonlinear fluid simulations are developed by us to investigate the properties of fully developed two-dimensional (2D) electron fluid turbulence in a very dense Fermi (quantum) plasma. We find that a 2D quantum electron plasma exhibits dual cascades, in which the electron number density cascades towards smaller turbulent scales, while the electrostatic potential forms larger scale eddies. The characteristic turbulent spectrum associated with the nonlinear electron plasma oscillations (EPO) is determined critically by a ratio of the energy density of the EPOs and the electron kinetic energy density of quantum plasmas. The turbulent transport corresponding to the large-scale potential distribution is predominant in comparison with the small-scale electron number density variation, a result that is consistent with the classical diffusion theory.     

Linking high- and low-temperature plasticity in bulk metallic glasses: thermal activation,extreme value statistics and kinetic freezing

Abstract

At temperatures well below their glass transition, the deformation properties of bulk metallic glasses are characterized by a sharp transition from elasticity to plasticity, a reproducible yield stress and an approximately linear decrease of this stress with increasing temperature. In the present work, it is shown that when the well-known properties of the undercooled liquid regime, in terms of the underlying potential energy landscape, are assumed to be also valid at low temperature, a thermal activation model is able to reproduce the observed onset of macroscopic yield. At these temperatures, the thermal accessibility of the complex potential energy landscape is drastically reduced, and the statistics of extreme value and the phenomenon of kinetic freezing become important, affecting the spatial heterogeneity of the irreversible structural transitions mediating the elastic-to-plastic transition. As the temperature increases and approaches the glass transition temperature, the theory is able to smoothly transit to the high-temperature deformation regime where plasticity is known to be well described by thermally activated viscoplastic models.     

n‐type conductivity in sublimation‐grown AlN bulk crystals

AlN crystals grown by physical vapour‐phase transport in the presence of a SiC doping source possess n‐type conductivity. The net donor concentration attains up to mid 10¹⁷ cm^–3. The investigation reveals shallow donors forming an impurity band and acceptor‐like electron traps at about 0.5 eV below the conduction band edge. Thermal electron emission from these traps is responsible for the observed n‐type conductivity. The shallow donors are suggested to be due to Si atoms on Al sites. The majority of them is assumed to be compensated by deep acceptors in the lower half of the band gap. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

太赫兹量子级联激光器中有源区上激发态电子向高能级泄漏的研究

利用热力学统计理论和激光器输出特性理论,建立了太赫兹量子级联激光器(THz QCL)有源区中上激发态电子往更高能级电子态泄漏的计算模型,以输出功率度量电子泄漏程度研究分析了晶格温度和量子阱势垒高度对电子泄漏的影响.数值仿真结果表明,晶格温度上升会加剧电子泄漏,并且电子从上激发态泄漏到束缚态的数量大于泄漏到阱外连续态,同时温度的上升也会降低激光输出功率.增加量子阱势垒高度能抑制电子泄漏,并且有源区量子阱结构中存在一个最优量子阱势垒高度. THz QCL经过最优量子阱势垒高度优化后,工作温度得到提升,其输出功率相比于以往的结果也有所提高.研究结果对优化THz QCL有源区结构、抑制电子泄漏和改善激光器输出特性有指导作用.     

双温度氩-氮等离子体热力学和输运性质计算

获得覆盖较宽温度和压力范围内的等离子体热力学和输运性质是开展等离子体传热和流动过程数值模拟的必要条件.本文通过联立Saha方程、道尔顿分压定律以及电荷准中性条件求解等离子体组分;采用理想气体动力学理论计算等离子体热力学性质;基于Chapman-Enskog方法求解等离子体输运性质.利用上述方法计算了压力为0.1, 1.0和10.0 atm (1 atm=101325 Pa),电子温度在300—30000 K范围内,非局域热力学平衡(电子温度不等于重粒子温度)条件下氩-氮等离子体的热力学和输运性质.结果表明压力和非平衡度会影响等离子体中各化学反应过程,从而对氩-氮等离子体的热力学及输运性质有较大的影响.在局域热力学平衡条件下,计算获得的氩-氮等离子体输运性质和文献报道的数据符合良好.