Numerical diffusion solver