Steady States in a Catalyst Pellet - T vs

Example 3.2.13. Multiple Steady States in a Catalyst Pellet - T] vs. 3  

The non-isothermal reaction in a rectangular pellet (example 3.2.2) is again considered. The dimensionless concentration is governed by  

For a given value of I , equation (3.63) can be solved for the boundary conditions ((3.65) and (3.66)). Once the numerical solution is obtained the effectiveness factor can be calculated using equation (3.67). It is of interest to plot the effectiveness factor as a function of O. For this purpose, one can solve the boundary value problem for different values of O and predict the effectiveness factor. However, since multiple steady states occur different initial guesses have to be used to capture all of the steady states. This boundary value problem is difficult to solve because for every value of O, there can 3 different qs and, hence. 

We know that for low values of O, the concentration at the center y(0) is close to 1. For high values of O, y(0) approaches zero. We take the concentration at the center as 

Warning, cannot evaluate the solution further right of. 20038679, stop condition 1 violated 

---