Charbel et al. (1999) have investigated the flow stability in conical spouted beds by modeling the effective solid stress distribution in the annular region as a function of the failure state under which the spout-annulus interface is formed. However, due to the complexity in solving the momentum balance equations for the solid phase, the model application has been restricted to one case. This work is aimed at developing a new methodology to solve this set of algebraic differential equations, using optimization techniques to determine the model parameter. This methodology has improved the model solution extending its application for any other case.