In this contribution, the kinetics of the relaxation process during the Hebb-Wagner polarization of a nominally undoped perovskite-type oxide ABO₃ is simulated by means of numerical methods. The advantages of the theoretical approach are the prediction of the expected relaxation time for experiments and the exact reproduction of the fluxes of all involved mobile species. Thus, it is possible to investigate the dependence of the relaxation time on different experimental variables, such as the applied voltage, acceptor concentration, temperature, and external oxygen partial pressure.