The diffusive scaling of many finite-velocity kinetic models leads to
a small-relaxation time behavior governed by reduced systems which are
parabolic in nature. Here we demonstrate that standard numerical meth
ods for hyperbolic conservation laws with stiff relaxation fail to cap
ture the right asymptotic behavior. We show how to design numerical sc
hemes for the study of the diffusive limit that possess the discrete a
nalogue of the continuous asymptotic limit. Numerical results for a mo
del of relaxing heat flow and for a model of nonlinear diffusion are p
resented.