File Exchange

image thumbnail

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

version 1.0.0.0 (39.2 KB) by Housam Binous
Prediction of the conversion of an isothermal axially dispersed plug flow reactor.

13 Downloads

Updated 10 Nov 2004

No License

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.

Cite As

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 .

Comments and Ratings (13)

David Han

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.

Jerry Che

Ethan Ma

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.

Cynthia Chu

Joel

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.

OCFE?

aymen mohamed

thank you for this efforts

hussein AHMED

sarma ravi krishna

hello
its good but not converging for high peclet numbers. by the way is this shooting point method??????? i am actually looking for OCFE

V. Laura Barrio

MATLAB Release Compatibility
Created with R13
Compatible with any release
Platform Compatibility
Windows macOS Linux