Laser ablation in an ambient environment is becoming increasingly
important in science and technology. It is used in applications ranging
from chemical analysis via mass spectroscopy, to pulsed laser deposition
and nanoparticle manufacturing. We describe numerical schemes for a
multiphase hydrodynamic model of nanosecond laser ablation expressing
energy, momentum, and mass conservation in the target material, as well
as in the expanding plasma plume, along with collisional and radiative
processes for laser-induced breakdown (plasma formation). Numerical
simulations for copper in a helium background gas are presented and the
efficiency of various ODE integrators is compared.