Laboratoire des matériaux composites céramiques et polymères (LaMaCoP) Faculté des sciences de Sfax BP 805 Sfax 3000 Tunisie. imed_boukhris@yahoo.frAbstract: This paper reports on validation of a numerical model for both transient and steady bipolar charge transport, trapping and recombination in polymeric insulators, such as low density polyethylene. The numerical model is based on the high precision Runge-Kutta method for determining mobile and trapped charge local density within the sample, during the application of a DC voltage. In addition, we have applied the finite element method considered as the most general one, which can be applied to any sample shape and boundary conditions. It is applied to resolve Poisson's equation, thus providing the potential and its gradient within the sample and at the dielectric-electrode interfaces. The principal results show the appearance of charge packets for the first time in modeling works although they have long been reported in experimental studies on polyethylene materials. The conditions for the appearance of such packets are discussed. For the net space charge density and the conduction current, the dynamics are in a good agreement with those observed in some experimental works during the transit time of the bipolar carriers. The model is already numerically validated for the space charge densities and the electric field distribution under low DC voltage.