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.
Housam Binous (2020). Prediction of Conversion with an Isothermal Axially-Dispersed Plug Flow Reactor (https://www.mathworks.com/matlabcentral/fileexchange/6206-prediction-of-conversion-with-an-isothermal-axially-dispersed-plug-flow-reactor), MATLAB Central File Exchange. Retrieved .
The clarity of the code could be improved with the addition of comments. Additionally, the code would be more versatile with the addition of different solvers to accommodate large Peclet numbers.
Code works but is not intuitive as there are no comments.
This code is sufficient for calculating the fractional conversion for a range of values. The code converges very efficiently with a high level of precision that is retained when various ODE solvers are implemented.
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