Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor

Numerical solutions for the conversion of an isothermal axially dispersed plug flow reactor accomplishing first and second order reactions are obtained. The Danckwerts boundary conditions for the so-called "closed" configuration are used. We use the shooting method since it is a split-boundary-value problem. Damkohler number (Da) and Inverse Peclet number (1/Pe) must have same order of magnitude unless one gets a stiff ODE which cannot be solved easily with ode45.