# Three-Dimensional Temperature Distribution Produced by a Moving Laser Beam

Authors: R. Uyhan
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

An axisymmetric laser beam, moving with constant speed, heats a
thin infrared absorbing layer sandwiched between two plastic sheets. We use a simplified theoretical model to study the three-dimensional
unsteady temperature field produced by the moving laser beam.


---

## Body

## 1. Introduction

High quality hard copy of electronically stored images can be produced by thermal dye diffusion printing [1]. Such copies are of near photographic quality and have been used in such diverse areas as medical imaging and credit card personalisation. A theoretical model of the laser induced thermal dye-diffusion process for a stationary axisymmetric laser was given by [2], based upon the ICI Imagedata L2T2 system as described by [3].In the L2T2 system a thin  2μm dye layer is sandwiched between a 6 μm layer of PET (polyethylene-tetrathalate) and a relatively thick receiver sheet of thickness 40 μm. Heat is supplied directly to the dye-coat layer by a laser (see Figure 1). A suitable infrared absorbing material such as substituted copper phthalocyanine is incorporated into the dye coat layer which absorbs radiation strongly at about 800 nm [3]. The dye therefore diffuses rapidly into the receiver sheet producing a small dye dot. The total amount of dye transferred depends on the temperature because the diffusion coefficient is known experimentally to be a sensitive function of temperature [4]. In the L2T2 system there is a thin barrier layer at the PET/dye layer interface which prevents dye diffusing into the PET film.Figure 1
Illustration of the geometrical situation (not to scale).The required image is formed by a large number of small dye dots produced in the receiver sheet as the pulsed laser beam moves relative to the layered system. In [2] a simplified model of the thermal dye diffusion process was studied in which the laser beam was fixed relative to the receiver sheet over the heating period typically 100 μs. The laser beam was then moved to the next pixel and the process was repeated and so on, producing a series of dots. It was assumed that there was no thermal interaction between adjacent pixels, so that we could consider the thermal dye diffusion process for a single typical pixel. In this case the time-dependent temperature and dye diffusion fields produced by the laser beam are axi-symmetric. A review of previous theoretical studies is given in [2, 5].In practice the axisymmetric laser beam is moved continuously relative to the receiver sheet and is switched on and off in a periodic fashion, so that a series of dots is produced. In this situation both the temperature and the dye diffusion fields will depend on three space coordinates and time, so the task of finding the temperature and dye fields is more difficult. In this paper we will use a simple mathematical model to study the temperature field.The temperature field produced by a Gaussian laser beam is of interest in other practical applications. Thus, for example, [6] studied the temperature field when the beam moves over the surface of a metal sheet of finite thickness, a study relevant to metal forming, while [7] studied the temperature field in a thin metal sheet, so that it could be assumed that the temperature was constant across it. In both these studies the beam had constant wattage and moved with constant speed, so that the temperature distribution took the form of a travelling wave.In the present case the polymer layers have relatively low conductivity. As a consequence the receiver sheet of thickness 40μm can be regarded as infinitely thick for all practical purposes over the time scales of interest, but the finite thickness of the top 6 μm PET layer, above the dye layer, has to be taken into account.


## 2. The Mathematical Model

Consider the layered system illustrated in Figure1. The z-coordinate is measured into the layered system from the air interface which is situated at z=0. The top (PET) carrier layer has thickness h1 and the dye-coat layer has thickness h2. The receiver sheet is relatively thick (about 40 μm) and can be regarded as infinitely thick in practice. We take a frame of reference in which the polymer sheet is stationary and the laser beam is moved with constant speed v in the positive x-direction over the surface z=0. We will assume that the axis of the laser is at x=0, y=0 at time t=0, so that at time τ, its axis will be at x=vτ, y=0. It is assumed in this paper that the thermal properties of the three layers shown in Figure 1 have the same constant values.Consider first the absorption of laser radiation by the infrared absorber in the dye-coat layer. Assuming that the distribution of the absorber is uniform, then the laser radiation flux per unit area,I(r,z,τ) at time τ, satisfies the differential equation:
(1)∂I∂z=-λh2I,
for h1≤z≤h1+h2 where λ is a constant. The quantity which is measured experimentally is the optical density which equals λ/ln10. Assuming that the radial dependence of the radiation flux on z=h1 is Gaussian, it follows that, at time τ,
(2)I(r,z,τ)=I0exp(-λ(z-h1)h2-r2r02),
where I0 depends on τ, r0 is a characteristic radius for the moving beam, and r2=(x-vτ)2+y2. The wattage of the beam W(τ) is found by integrating 2πrI(r,z,τ) radially over the beam at z=h1, and we find that W=πr02I0. It follows that the laser energy absorbed per unit volume per unit time in the dye layer is
(3)-∂I∂z=W(τ)λπr02h2exp(-λ(z-h1)h2-r2r02).
This is of course the largest at the top of the dye layer, z=h1. We note that the total energy absorbed per unit area per unit time is
(4)Wπr02(1-e-λ)exp(-r2r02).
This is the largest when λ=∞, in which case all the heat is generated in an infinitely thin layer at z=h1.Since, by assumption, the thermal properties have the same constant values in the three layers, we can construct the temperature field produced by the moving laser by combining instantaneous source solutions.At timeτ and at z=z′, the laser beam generates a heat source distribution in the dye layer q(τ,z′)exp(-r2/r02)dz′dτ over a short time dτ in a thin layer of thickness dz′.Then from (3),
(5)q(τ,z′)=W(τ)λπr02h2exp[-λ(z′-h1)h2].
If the medium is infinite, has uniform properties, and has zero initial temperature, the temperature field for t>τ is given by
(6)G(x,y,z,t,z′,τ)=q(τ,z′)γ(u)2ρcpπαuexp[-(z-z′)24αu-γ(u)r2r02],
where u=t-τ,  γ=r02/(r02+4αu), α is the thermal diffusivity, ρ is the density, and cp is the specific heat.The solution (6) satisfies the heat conduction equation
(7)∂G∂t=α{∂2G∂z2+1r∂∂r(r∂G∂r)}.
As t→τ, we can see that G→0, except near z=z′, where it is large, and then (6) takes the approximate form
(8)G≃q(τ,z′)2ρcpπαuexp[-(z-z′)24αu-r2r02],
which is the corresponding “one-dimensional” heat source solution of
(9)∂G∂t=α∂2G∂z2
(see, e.g., [8, 9]). If instead heat is produced throughout the dye-coat layer h1<z<h1+h2, the corresponding temperature field T can be found by integrating the instantaneous heat source solution G with respect to z′ from h1 to h1+h2 and with respect to τ from -∞ to the current time t.Thus we find that(10)T(x,y,z,t)=∫-∞t‍∫h1h1+h2‍q(τ,z′)γ(u)2ρcpπαu×exp[-(z-z′)24αu-γ(u)r2r02]dz′dτ,
where q(τ,z′) is defined in (5). This is the temperature field produced by a moving laser in an infinite medium initially at zero temperature.In the real problem, the initial temperature isTA and there is no heat flux at z=0, the air-polymer interface. To incorporate this boundary condition into the model, we introduce an appropriate image heat source distribution in z<0. Then by symmetry, there will be no heat flux across z=0. The effect of this is to give, after some manipulation,
(11)T-TA=∫-∞t‍∫h1h1+h2‍q(τ,z′)γ(u)2ρcpπαu×f(z,z′,u)exp-γ(u)r2r02dz′dτ,
where
(12)f(z,z′,u)=exp(-(z-z′)24αu)+exp(-(z+z′)24αu).
In normal operation the laser beam is switched on and off in a periodic fashion. Assuming that the period is τp and the laser is switched on for a time τon and that when it is operating W=W0, a constant, then we can take
(13)W(τ)={W0,a(n)≤τ<b(n),0,otherwise,
where a(n)=nτp, b(n)=min(a(n)+τon,t), and n takes all integer values from -∞ to n=n*, where a(n*)≤τ<a(n*+1). Thus
(14)T-TA=Q0∑n=-∞n*‍∫a(n)b(n)‍∫h1h1+h2‍γ(u)πατf(z,z′,u)×exp[-λ(z′-h1)h2-γ(u)r2r02]dz′dτ,
where u=t-τ, r2=(x-vτ)2+y2,  and Q0=W0λ/2πr02h2ρcp.For the L2T2 system, typical values of physical parameters arev= 0.15 m/s, τp=150 μs,  τon=100 μs, ρ=1.3×103 kg/m3, cp=1.9×103 J/kg·K, α=3.2×10−8 m2 s−1, h1=6 μm, h2=2 μm, W0=75 mW , r0=15 μm, μ=0.7, λ=1, and TA=20°C.


## 3. Results and Discussion

The double integral (14) was evaluated using a standard NAG library routine. For the values of the parameters chosen, we found that the number of terms in the infinite series is small. This is because the axis of the beam passes over any given point in a time of order 2r0/v= 200 μs, which is comparable to the heating time τon=100 μs. Thus only a few laser pulses can affect the temperature at a given point.Consider the temperature field aty=0 near x=0. The axis of the laser passes over x=0 at time τ=0. However, the temperature field near x=0 will depend on the precise times when the laser was switched on and off. In our calculations we have assumed that the laser is switched on at τ=0, ±τp, 2±τp…. In particular this means that the laser was turned off for −50 μs <τ<0 if τon=100 μs, and consequently the maximum temperature in the dye coat layer cannot be at x=0. Numerical calculation shows that the maximum temperature in the dye coat layer is near x= 2.5 μm, and hence our numerical results are shown for this value of x. The position of maximum temperature will depend on the parameter values chosen.Thez-temperature profile at various times for x=2.5 μm, y=0 is shown in Figure 2. It can be seen that the temperature rises between −100 μs and −50 μs but then falls until t=0 because the laser is turned off. The temperature then rises rapidly to about 900°C over the next 100 μs because the laser is switched on at τ=0, and its axis is close to x=2.5 μm. The laser is then turned off again and the temperature falls slowly to 700°C at t=200 μs, as heat conducts upwards and downwards from the dye coat layer, with negligible heat input from the laser. This is shown in Figure 3.Figure 2
Temperature profiles onx= 2.5 μm, y=0 as a function of z at times −100, −50, 0, 50, 100, and 200 μs.Figure 3
Temperature profiles onx= 2.5 μm, y=0 as a function of z at times 250, 300, 350, 400, 500, 550, and 600 μs.Thez-temperature profiles shown in Figures 2 and 3 are complicated mainly because they show that at least two laser pulses contribute significantly. These profiles are much more complicated than these for a stationary laser (see Figures 4 and 5) reproduced from our paper [2] for the same heating time. Note that the maximum temperature is about 1000°C in this case. This suggests that for these standard values, the process of dye diffusion is not well described by our stationary laser model. In general the maximum temperatures will be reduced as v increases from zero and so will be dye diffusion, which is a sensitive function of temperature.Figure 4
Temperature profiles onr=0 at times t=25,50,75, and 100 μs for an uniformly distributed absorber in the dye layer (heating phase). The laser is stationary.Figure 5
Temperature profiles onr=0 at times t=100, 200, 300, and 400 μs for an uniformly distributed absorber in the dye layer (cooling phase). The laser is stationary.The effect of an increase in the speed of the laser to 0.5 m/s is illustrated in Figures6 and 7. As before, the laser is turned off at τ=-50 μs, but now the laser beam axis is at vτ=-25 μm. The heating before this time is therefore unimportant since x=2.5 μm is well outside the laser beam. Note that we have chosen to show the z-temperature profiles at x=2.5 μm, as in the standard case. The maximum temperatures will be at a slightly larger value of x.Figure 6
Temperature profiles onx= 2.5 μm, y=0, and v=0.5 as a function of z at times 0, 50, 100, 150, and 200 μs.Figure 7
Temperature profiles onx= 2.5 μm, y=0, and v=0.5 as a function of z at times 250, 300, 350, 400, 450, 500, 550, and 600 μs.The laser is turned on atτ=0 for 100 μs, but for much of that time, little heating occurs at x=2.5 μm since the laser beam has moved further to the right. Thus much lower temperatures of order 400°C are achieved at t=50 μs. Figure 5 shows the typical cooling curves as heat conducts upwards and downwards from the dye coat layer. In Figure 8 the temperature profiles on y=0, z=7 micro meter and v=0.1 m/s as a function of x for the different of time.Figure 8
Temperature profiles onz=7, y=0, and v=0.1 as a function of x at times 0, 50, 100, 150, 200, and 250 μs.The diffusion coefficient will be much smaller at these lower temperatures, but there will be negligible interaction with other pulses. Higher temperatures could of course be achieved by increasing the laser wattage.


## 4. Conclusions

We have developed a simple theoretical model which allows us to calculate thez-temperature field in the dye coat layer produced by a moving laser beam. Our main conclusion is that, for the standard data, the temperature at a given point is significantly influenced by at least two pulses leading to a complicated z-temperature profile in the dye coat layer as time changes. This implies that the dye-diffusion process is much more complicated than that predicted by a stationary laser model.



---

*Source: 649590-2013-11-27.xml*