We derive diffusive macroscopic equations for the particle and energy density of a system whose time evolution is described by a kinetic equation for the one particle position and velocity function f(r,v,t) that consists of a part that conserves energy and momentum such as the Boltzmann equation and an external randomization of the particle velocity directions that breaks the momentum conservation. Rescaling space and time by epsilon and epsilon square respectively and carrying out a Hilbert expansion in epsilon around a local equilibrium Maxwellian yields coupled diffusion equations with specified Onsager coefficients for the particle and energy density. Our analysis includes a system of hard disks at intermediate densities by using the Enskog equation for the collision kernel.