Diffusive BGK approximations for nonlinear multidimensional parabolic equations