International efforts have focused recently on the development of tungsten surfaces that can intercept energetic ionized and neutral atoms, and heat fluxes in the divertor region of magnetic fusion confinement devices. The combination of transient heating and local swelling due to implanted helium and hydrogen atoms has been experimentally shown to lead to severe surface and sub-surface damage. We present here a computational model to determine the relationship between the thermo-mechanical loading conditions, and the onset of damage and failure of tungsten surfaces. The model is based on thermo-elasticity, coupled with a grain boundary damage mode that includes contact cohesive elements for grain boundary sliding and fracture. This mechanics model is also coupled with a transient heat conduction model for temperature distributions following rapid thermal pulses. Results of the computational model are compared to experiments on tungsten bombarded with energetic helium and deuterium particle fluxes.