A Computationally Efficient Two-Stage Two-Temperature Model for Multi-pulse Femtosecond Laser Heat Accumulation in Tungsten
摘要
Ultrafast laser processing at megahertz repetition rates produces cumulative inter-pulse heat buildup that governs damage thresholds, melt onset, and surface quality. The two-temperature model (TTM) resolves the nonequilibrium electron–lattice physics of each pulse, but becomes impractical for multi-pulse sequences typical of MHz-rate processing. We present a two-stage decoupling method that separates each pulse cycle into a stiff TTM integration during the pulse and relaxation period, followed by Crank–Nicolson diffusion through the inter-pulse gap using temperature-dependent thermal conductivity from tabulated data. Applied to tungsten, the solver completes a 10,000-pulse simulation in under 17 min on a single workstation, approximately 600 times faster than the most efficient previously reported spatially resolved multi-pulse TTM. Multi-pulse simulations predict progressive inter-pulse temperature buildup and a growing surface temperature inversion (