We study the propagation of radiative heat (Marshak) waves, using modified P 1 -approximation equations.In relatively optically-thin media the heat propagation is supersonic, i.e. hydrodynamic motion is negligible, and thus can be described by the radiative transfer Boltzmann equation, coupled with the material energy equation. However, the exact thermal radiative transfer problem is still difficult to solve and requires massive simulation capabilities. Hence, there still exists a need for adequate approximations that are comparatively easy to carry out. Classic approximations, such as the classic diffusion and classic P 1 , fail to describe the correct heat wave velocity, when the optical depth is not sufficiently high. Therefore, we use the recently developed discontinuous asymptotic P 1 approximation, which is a time-dependent analogy for the adjustment of the discontinuous asymptotic diffusion for two different zones. This approximation was tested via several benchmarks, showing better results than other common approximations, and has also demonstrated a good agreement with a main Marshak wave experiment and its Monte-Carlo gray simulation.Here we derive energy expansion of the discontinuous asymptotic P 1 approximation in slab geometry, and test it with numerous experimental results for propagating Marshak waves inside low density foams. The new approximation describes the heat wave propagation with good agreement. Furthermore, a comparison of the simulations to exact IMC slab-geometry multi-group simulations, in this wide range of experimental conditions, demonstrates the superiority of this approximation to others.