Modeling of thermal phenomena in single laser beam and laser-arc hybrid welding processes using projection method

This paper concerns mathematical and numerical modeling of thermal phenomena accompanying single laser and laser-arc hybrid butt welding of steel sheets. Coupled heat transfer and fluid flow in the fusion zone were described respectively by transient heat transfer equation and Navier–Stokes equation. Laser beam and electric arc heat sources were modeled using different heat source power distributions. Latent heat associated with the material’s state changes, buoyancy forces and liquid material flow through a porous medium were taken into account in considerations. Differential governing equations were numerically solved using projection method combined with finite volume method. Elaborated solution algorithm was implemented into computer solver used for simulation of heat transfer and fluid flow during welding. The geometry of the weld and heat affected zone as well as cooling rates were estimated on the basis of numerically obtained temperature field.     

Use of geometry in finite element thermal radiation combined with ray tracing

In heat transfer for space applications, the exchanges of energy by radiation play a significant role. In this paper, we present a method which combines the geometrical definition of the model with a finite element mesh. The geometrical representation is advantageous for the radiative component of the thermal problem while the finite element mesh is more adapted to the conductive part. Our method naturally combines these two representations of the model. The geometrical primitives are decomposed into cells. The finite element mesh is then projected onto these cells. This results in a ray tracing acceleration technique. Moreover, the ray tracing can be performed on the exact geometry, which is necessary if specular reflectors are present in the model. We explain how the geometrical method can be used with a finite element formulation in order to solve thermal situation including conduction and radiation. We illustrate the method with the model of a satellite.     

微蜂窝环境中无线接收信号的特性分析

以微蜂窝系统为背景,研究了以测试射线法和虚拟源射线跟踪法为代表的射线跟踪技术,并用这两种方法对射线跟踪进行建模,分析比较两种模型的优缺点.利用这两种模型对实际情况进行计算,精确地找出所有给定情况的传播路径并可视化;对多波干涉的振幅进行了统计学分析,对宽带多波干涉现象进行了数学建模,并分析了合成波的包络统计特性.     

Using Multibody Systems for the Investigation of Dynamic Aberrations in High Precision Optics

The coupling of multibody system dynamics and optical simulations by using ray tracing through moving lens systems will be summarized for both rigid and flexible lens systems. Furthermore, a method will be introduced for efficient simulations in the time domain. This method provides a direct integration of the optical simulation into the equations of motion of a multibody system. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Thermomechanical Simulation of the Electron Beam Melting Process for TiAl6V4

In the present contribution a nonlinear thermomechanical finite element model with temperature dependent material parameters is used to simulate the electron beam melting process for TiAl6V4. The beam is modeled as a moving heat source and the isentropic split solution scheme is applied. The temperature and stress distribution during the process were simulated and showed sound results from a qualitative point of view. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

3D finite element analysis of evaporative laser cutting

A three-dimensional computational model of evaporative laser-cutting process has been developed using a finite element method. Steady heat transfer equation is used to model the laser-cutting process with a moving laser. The laser is assumed continuous wave Gaussian beam. The finite element surfaces on evaporation side are nonplanar and approximated by bilinear polynomial surfaces. Semi-infinite elements are introduced to approximate the semi-infinite domain. An iterative scheme is used to handle the geometric nonlinearity due to the unknown groove shape. The convergence studies are performed for various meshes. Numerical results about groove shapes and temperature distributions are presented and also compared with those by semi-analytical methods.     

Finite element analysis of evaporative cutting with a moving high energy pulsed laser

A computational model for simulation of pulsed laser-cutting process has been developed using a finite element method. An unsteady heat transfer model is considered that deals with the material-cutting process using a Gaussian wave laser beam in a pulsed mode. An iterative scheme is used to handle the geometric nonlinearity due to the melting region. The convergence study with mesh refinements and time steps first identifies optimal mesh and time step for the present analyses. Numerical analyses are carried out on the amount of material removal and groove smoothness with laser power (LP) and number of pulses (NPS) while other laser cutting parameters are fixed. The results show that there exist threshold values in number of pulses and laser power in order to achieve two predetermined conditions: (1) amount of material removal and (2) smoothness of groove shape. These values form an envelope called threshold curve that separate the acceptable region from unacceptable one for quality pulsed laser cutting. The effect of velocity also leads to another threshold curve which is determined from both number of pulses and velocity. Finally, the convergence of results in error domain is shown oscillating due to geometric nonlinearity.     

基于分数阶热传导方程激光加热瞬态温度场研究

基于分数阶Taylor(泰勒)级数展开原理,建立单相延迟一阶分数阶近似方程,获得分数阶热传导方程.针对短脉冲激光加热问题建立分数阶热传导方程组,并运用Laplace(拉普拉斯)变换方法进行求解,给出非Gauss(高斯)时间分布的激光内热源温度场解析解.针对具体算例数值研究温度波传播特性.结果表明热传播速度与分数阶阶次有关,分数阶阶次增加,热传播速度减小,温度变化幅度增加.分数阶方程可以用于描述介于扩散方程和热波方程间的热传输过程,且对热传播机制与分数阶热传导方程中分数阶项的关系做了深入剖析.     

MODELING THE EFFECTS OF URBAN HEAT ISLANDS GENERATED MESOSCALE WIND ON AIR POLLUTION DISPERSION IN A PATCHY ATMOSPHERE

Abstract The paper presents a steady‐state two‐dimensional mathematical model to study the dispersion of air pollutants continuously emitting from a point source in a patchy atmosphere. The first patch, which extends to some downwind distance from the source (present in first patch), is the urban center, whereas the second patch is its adjacent countryside. For simulating the effects of urban heat islands on pollutant distribution, a local mesoscale wind is introduced in the first patch along with the usual large‐scale wind. Keeping in mind the real‐life situations, the meteorological, deposition, and removal parameters are taken to be different in different patches. The concentration of pollutants in both the patches has been computed under different stability conditions of the atmosphere through finite difference scheme. The results of the model show that the urban heat island (mesoscale wind) has a tremendous impact on air pollution distribution, both over the urban center as well as its adjacent countryside (rural). The model solutions also reveal about the concentration distribution of pollutants in different parts of the atmosphere in different atmospheric conditions.     

尖劈吸波体的研究和微波暗室的模拟

对尖劈形状吸波体的吸波性能进行了研究,并对导弹导引仿真实验用的微波暗室的性能进行分析和仿真.首先对尖劈体的二维反射性能进行了研究,从单条波线反射的原理出发,得到波束平面反射的统计模型.单条波线的反射通过数值模拟得到;波束反射模型则通过对数值模拟的结果进行统计和拟合得到,最终用多项式表示.对于一些简单或特殊的情况,也给出了解析解.通过分析发现,三维反射和二维反射之间有明确的关系.这种关系可以由三维入射角和反射次数决定,而反射次数可以通过二维模型得到.据此将平面反射模型扩展为三维反射模型,从而得到尖劈形状吸波体的三维反射模型.无回波暗室用于模拟没有背景微波辐射的环境,其关键在于选择合适的吸波材料.基于微波反射通量平衡原理,建立了考虑暗室墙面各点之间的相互影响的耦合模型,从而可以求解出在指定的发射源照射之下墙面各点的辐射强度分布.对模型的求解精度和收敛性进行了验证.基于此模型,对一个导弹引导试验进行了数值模拟,推算出了使用两种不同吸波材料时静区接收到的微波信号的信噪比.     

Different Optimization Models for Crack-Free Laser Welding

Conventional laser beam welding of aluminum alloys often leads to hot cracking. This is caused by a complex process where thermo-mechanical and metallurgical aspects are involved; cf. [3], [2]. A possibility to prevent hot crack initiation yields the multi-beam welding technique (cf. [2]), where additional laser beams are led parallelly besides the main laser beam. There by optimal positions, sizes, and powers of the additional laser beams play an important role otherwise hot cracking can even be enhanced. In [1], [4], resp., a mechanical 1D and thermal 2D model of hot cracking was derived. It provides the basis for different formulations of constrained nonlinear programming problems to identify the optimal parameters of the additional laser beams. In the present paper a comparison between these formulations and between two different optimizers for the so far best formulation are presented. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Fractional order theory of Cattaneo-type thermoelasticity using new fractional derivatives

Transient thermoelastic interactions between materials and the moving heat sources, i.e. Laser additive manufacturing, Laser-assisted thermotherapy, high speed sliding and rolling contacts, are becoming increasingly important. In this work, a unified fractional thermoelastic theory is developed, and applied to study transient responses caused by a moving heat source. Theoretically, new insights on fractional thermoelasticity are provided by introducing new definitions of fractional derivative, i.e. Caputo-Fabrizio, Atangana–Baleanu and Tempered-Caputo type. Numerically, a semi-infinite medium subjected to a source of heat moving with constant velocity is considered within the present model under two different sets of boundary conditions: stress free and temperature given for the first, displacement fixed and thermally adiabatic for the second. Analytical solutions to all responses are firstly formulated in Laplace domain, and then transformed into time domain through numerical method. The numerical results show that Caputo-Fabrizio and Atangana–Baleanu type models predict smaller transient responses than Caputo type theory, while Tempered-Caputo model may give larger results by increasing the tempered parameter. Meanwhile, the effect of fractional order, tempered parameter of Tempered-Caputo model, and the velocity of heat source on all responses is discussed in detail. The time history of responses shows that: for long-term process, the exponential function of TC definition will make sense, and the temperature from TC model is greatly different from that of C model. This work may provide comprehensive understanding for thermoelastic interactions due to moving heat source, and open up possibly wide applications of such new fractional derivatives.     

Monte Carlo simulation of radiative transfer in a medium with varying refractive index specified at discrete points

A Monte Carlo method with discrete ray tracing is developed to simulate radiative transfer in a medium with a spatially varying refractive-index distribution available merely at a set of arbitrary discrete points. We solve the ray equation by a Runge−Kutta Dormand−Prince method to carry out the numerical ray tracing. To retrieve the refractive-index values and gradients needed in the discrete ray tracing, we apply cubic spline interpolation for one-dimensional simulation and a moving least square (MLS) method for two-dimensional simulation. The influence of the basis vectors and the numbers of sampled data used by the MLS method on ray tracing based on the retrieved refractive-index values and gradients has been examined. The results of radiative equilibrium in a planar medium and radiative transfer in two-dimensional media with different geometries and conditions obtained by the present methods are compared with those obtained by solving the integral equations of radiative transfer and the discrete ordinates method. The comparisons show that the present methods generate accurate results for radiative transfer with various geometries, parameters and refractive-index distributions specified at discrete points.     

A Green’s function model for the analysis of laser heating of materials

An analytical model based on Green’s function method is developed to analyze the temperature distribution and heated regions in a material irradiated by a high-energy laser beam. The model is multi-dimensional, transient and incorporates different types of beam characteristics and boundary conditions. The multi-dimensional integration formulas in the Green’s function solution equation are evaluated using an adaptive numerical integration algorithm. A parametric study is conducted to show the effect of various laser beam parameters and material properties on the laser heating process.     

Skew ray tracing and sensitivity analysis of ellipsoidal optical boundary surfaces

The 4 × 4 homogeneous transformation matrix is one of the most commonly applied mathematical tools in the fields of robotics, mechanisms and computer graphics. Here we extend further this mathematical tool to geometrical optics by addressing the following two topics: (1) skew ray tracing to determine the paths of reflected/refracted skew rays crossing ellipsoidal boundary surfaces; and (2) sensitivity analysis to determine via direct mathematical analysis the differential changes of the incident point and the reflected/refracted vector with respect to changes in the incident light source. The proposed ray tracing and sensitivity analysis are projected as the nucleus of other geometrical optical computations.     

Mechanical Analysis of Metallic SLM-Lattices on Small Scales: Finite Element Simulations versus Experiments

The present contribution deals with the mechanics of metallic microlattices from selective laser melting (SLM). Finite element analyses with elasto-plastic material parameters identified in experiments investigate the structural load bearing behavior of different unit cell topologies. Typical failure modes like local buckling as well as global localization in shear bands are analyzed in simulations and experiments for compression tests. Ashby diagrams for the scaling behavior of stiffness and strength at various densities are determined for both bending- and stretch-dominated lattice types. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The gravitational field of a laser beam and the Generalized Uncertainty Principle

Light carries energy and therefore is source of a gravitational field. In the paper it is proven, in particular, that a laser beam is source of non-linear gravitational waves corresponding, from a quantum point of view, to spin-1 massless particles. This fact suggests both a possible solution to the old problem on the lack of gravitational attraction between two laser beams moving parallel and a new derivation of the Generalized Uncertainty Principle.     

On the Motion of an Intensely Heated Evaporating Boundary

In deep penetration welding and cutting (by a laser or an electronbeam) the surface of a metal is strongly evaporated and a holedevelops. The maximum rate of penetration into a given materialis simply shown to be proportional to beam power density. Thereis an initial pre-heating time before the surface starts toevaporate which for any given material depends upon power density(according to an inverse square law). The motion of the evaporatingboundary is dominated and complicated by the discontinuity whichlatent heat represents. However, solutions can be found by usinga perturbation expansion and these all tend to the evaporation-controlledlimit. The maximum velocity is effectively achieved in aboutten times the pre-heat time (depending upon the latent heatin relation to the normal heat capacity). For a stationary source(i.e. not traversing the surface) it is shown that energy considerationsalone imply that the depth of the hole increases without limit.     

微蜂窝环境中无线接收信号的特性分析研究

基于镜像原理,建立了反向射线追踪模型,同时基于有效反射面的概念,给出了辨别有效像点的充要判据.在包含绕射点路径的处理上,提出绕射路径分解方法,使得模型可以处理多次绕射问题.证明了二维情形下能量与场强之间存在正比关系.通过分析视距传播中的反射情形和非视距传播中的绕射情形,给出计算给定传播途径上接收点处场强的离散化方法.利用单频问题中Rayleigh分布规律,提出通过半经验模型得到的合场强幅值估算概率密度函数的方法.