The BFKL pomeron in perturbative QCD is plagued by the lack of unitarity and diffusion into the infra-red region of gluon virtualities. These two problems are intimately related. We perform numerical studies of the evolution equation proposed by Balitsky and Kovchegov which unitarizes the BFKL pomeron. We show how diffusion into the infra-red region is suppressed due to the emergence of a saturation scale and scaling behaviour. We study universality of this phenomenon as well as its dependence on subleading corrections to the BFKL pomeron such as the running coupling and kinematic constraint. These corrections are very important for phenomenological applications. (orig.)Available from TIB Hannover: RA 2999(2001-172) / FIZ - Fachinformations...