  We consider two models which were both designed to describe the movement of
eukaryotic cells responding to chemical signals. Besides a common standard
parabolic equation for the diffusion of a chemoattractant, like chemokines or
growth factors, the two models differ for the equations describing the movement
of cells. The first model is based on a quasilinear hyperbolic system with
damping, the other one on a degenerate parabolic equation. The two models have
the same stationary solutions, which may contain some regions with vacuum. We
first explain in details how to discretize the quasilinear hyperbolic system
through an upwinding technique, which uses an adapted reconstruction, which is
able to deal with the transitions to vacuum. Then we concentrate on the
analysis of asymptotic preserving properties of the scheme towards a
discretization of the parabolic equation, obtained in the large time and large
damping limit, in order to present a numerical comparison between the
asymptotic behavior of these two models. Finally we perform an accurate
numerical comparison of the two models in the time asymptotic regime, which
shows that the respective solutions have a quite different behavior for large
times.
