A phenomenon based model for chlorine dioxide delignification of chemical pulp is introduced. The pulp suspension environment is modeled using the concept of two liquid phases, one inside and the other external to the fiber wall. Physico-chemical processes taking place during delignification are implemented with thermodynamic, mass transfer and reaction kinetic models. A broad library of chemical reactions is introduced. Inclusion of each reaction is justified. The model response is tested against experimental laboratory delignification results (o-delignified birch pulp). The experimental data consists of kappa number, hexenuronic acid, inorganic oxy-chlorine compound, and organochlorine (AOX, OX) measurements at several time points during five delignification experiments. The model predictions are mainly in good agreement with the experimental results. The predictions regarding hypochlorous acid driven processes (HexA removal, organochlorine formation, chlorite and chlorate concentration) are somewhat incoherent, indicating that knowledge regarding the intermediately formed hypochlorous acid is presently insufficient.