Thermal Diffusion & Heat Accumulation

How far heat spreads during a laser pulse, whether it stays within the optical absorption depth, and how fast you can pulse before heat builds up from one pulse to the next. Includes room-temperature diffusivities for common metals, glasses, and polymers.

Inputs

δ = 1/α = λ/(4πκ). Metals in the near-IR are roughly 10–30 nm. Transparent materials range from micrometres to metres.

Results

Thermal diffusion lengthLth = 2√(Dτ)—
Alternative convention√(Dτ)—
Lateral spread relative to beam radiusLth/w₀—
Thermal confinementLth vs. δ—
Heat-spreading time across spottw = w₀²/4D—
Heat-accumulation onsetf ≈ 4D/w₀²—
Conservative onset (full diameter)4D/d²—
Pulse period vs. tw—

Diffusion length vs. pulse duration

How it is calculated

Heat released at a surface spreads with the diffusion equation. After a time \(t\), the temperature profile from an instantaneous planar source falls as \(\exp(-x^2/4Dt)\), so it drops to 1/e at

\[ L_\text{th} = 2\sqrt{D\tau}, \qquad D = \frac{k}{\rho\, c_p}, \]

where \(k\) is the thermal conductivity, \(\rho\) the density, and \(c_p\) the specific heat. Authors differ on the prefactor: \(\sqrt{D\tau}\), \(\sqrt{2D\tau}\), and \(2\sqrt{D\tau}\) all appear. Both common conventions are shown, so check which one a source uses before comparing numbers. Diffusion length scales only with \(\sqrt{\tau}\). Going from 10 ns to 10 ps shrinks it by a factor of about 32, not 1000.

Thermal confinement

If \(L_\text{th} \lesssim \delta\), the optical penetration depth, heat stays in the volume where it was absorbed during the pulse. The pulse is thermally confined. Energy density then builds up efficiently and ablation is relatively clean, with a small heat-affected zone. If \(L_\text{th} \gg \delta\), heat leaks into the bulk during the pulse. More energy is needed to reach ablation, and melt and recast appear. For metals, with \(\delta \approx\) 10–30 nm, confinement requires picosecond or shorter pulses.

Heat accumulation between pulses

After each pulse, heat must spread out of the irradiated spot before the next pulse arrives. Spreading over the beam radius takes about

\[ t_w \approx \frac{w_0^2}{4D} \quad\Longrightarrow\quad f_\text{acc} \approx \frac{4D}{w_0^2}. \]

Above roughly \(f_\text{acc}\), residual heat from earlier pulses raises the baseline temperature. The process then behaves more like a quasi-CW heat source: the heat-affected zone grows, melt and burr formation increase, and the effective threshold drops. Using the full diameter instead of the radius gives a 4× lower, more conservative onset. The transition is gradual because residual heat decays slowly, so treat \(f_\text{acc}\) as an order-of-magnitude guide. Pulse overlap raises the number of pulses per point and makes accumulation worse (see the overlap tool).

Material properties

Materialk (W/m·K)ρ (kg/m³)cp (J/kg·K)D (mm²/s)
Copper4018933385117
Aluminum (pure)237270290397.1
Silicon148233071289.2
Stainless steel 30414.979004773.95
Ti-6Al-4V6.744305262.88
Fused silica1.3822037400.85
Borosilicate (Borofloat 33)1.222308300.65
Soda-lime glass1.025007500.53
Water (25 °C)0.60799741810.146
PMMA0.19119014200.11

These are typical room-temperature values. Alloys, glasses, and polymers vary with composition and supplier, and every material's diffusivity changes with temperature. Most metals' diffusivity falls as they heat toward melting. Use measured data for your material when accuracy matters.

Worked example

Stainless steel 304 (D = 3.95 mm²/s) with 10 ns pulses: \(L_\text{th} = 2\sqrt{3.95\ \text{mm}^2/\text{s}\times 10\ \text{ns}} \approx 0.40\ \mu\text{m}\). That is about 20× the optical penetration depth, so the process is not thermally confined. The same steel with 10 ps pulses gives 12.6 nm, which is comparable to δ. For a 30 µm spot, \(t_w = (15\ \mu\text{m})^2/(4D) = 14\ \mu\text{s}\), so heat accumulation sets in around 70 kHz. At 100 kHz with significant overlap, expect accumulation. In fused silica the same spot accumulates heat from about 15 kHz. In copper the onset is about 2 MHz.

Assumptions and limits

References

  1. T. L. Bergman, A. S. Lavine, F. P. Incropera, and D. P. DeWitt, Fundamentals of Heat and Mass Transfer, 7th ed., Wiley (2011), Appendix A (property data).
  2. H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids, 2nd ed., Oxford University Press (1959).
  3. D. Bäuerle, Laser Processing and Chemistry, 4th ed., Springer (2011).