A Monte Carlo simulation in 3D voxel space, extended with a finite element simulation for temperature increase, heat diffusion and thermal damage in the voxel space. Primarily targeted for tissue optics, but can be used in any environment.