Use of orthogonal collocation methods for the modeling of catalyst particles—I. Analysis of the multiplicity of the solutions