In this work we numerically study the diffusive limit of run & tumble kinetic models for cell motion due to chemotaxis by means of asymptotic preserving schemes. It is well-known that the diffusive limit of these models leads to the classical Patlak-Keller-Segel macroscopic model for chemotaxis. We will show that the proposed scheme is able to accurately approximate the solutions before blow-up time for small parameter. Moreover, the numerical results indicate that the global solutions of the kinetic models stabilize for long times to steady states for all the analyzed parameter range. We also generalize these asymptotic preserving schemes to two dimensional kinetic models in the radial case. The blow-up of solutions is numerically investigated in all these cases.Hereis the density of chemoattractant, whose governing equation will be described later.A very interesting consequence is that even in one dimension, this system blows-up over this critical mass. This makes possible to numerically analyse this highly subtle phenomena. This modified Keller-Segel system was numerically solved in [3] by using optimal transportation ideas and the scheme proven to be convergent in the subcritical case.Alt in [1, 2] derived (1.4) from a transport equation similar to (1.1) via a stochastic model. Later the kinetic system (1.1)-(1.2) was formulated in [26]. Othmer and Hillen in [16,27] studied the formal diffusion limit of this system by moment expansions. In [11] Chalub et al. proved that, in three dimensions with suitable assumptions on the turning kernel T ε , the solution to the kinetic system (1.1)-(1.2) globally exists for any initial total mass. Moreover, they gave a rigorous proof of the convergence to the macroscopic limit for all time invervals of existence of the limiting Keller-Segel system. In the following work [18,19] the results are extended to more general cases.