Use of orthogonal collocation methods for the modeling of catalyst particles—II analysis of stability