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.
Thanks, it was very helpful for an initial analysis of my problem (low Peclet numbers). However, as already mentioned by other user, it does not converge when higher Peclet numbers are used. I would advise everyone who might have this problem to use bvp4c matlab function instead.
thank you for this efforts
its good but not converging for high peclet numbers. by the way is this shooting point method??????? i am actually looking for OCFE