WOLFRAM NOTEBOOK

WOLFRAM|DEMONSTRATIONS PROJECT

Concentration and Temperature Profiles in an Axially Dispersed Adiabatic Tubular Reactor

concentration profile
temperature profile
low Péclet region
high Péclet region
Pe=20
low Péclet regionPéclet number
1.5
high Péclet regionPéclet number
20
This Demonstration computes the concentration and temperature profiles in a tubular reactor for low Péclet numbers (
Pe)
using the shooting method (blue solid curve), as well as the orthogonal collocation method with shifted Legendre polynomials and eight internal collocation points (the red dots). Perfect agreement is obtained (see the first two snapshots).
For large values of the
Pe
number, the reactor behavior approaches that of a plug flow reactor (PFR), and in the limit as
Pe
, the system become singular. Consequently the boundary value problem becomes increasingly stiff with increasing
Pe
. For
Pe>3
, the shooting method fails to converge. Thus for
Pe>3
, only orthogonal collocation results are shown (see red dots). To assess the accuracy of the collocation method, we computed the residual error for the temperature and concentration profiles. The residual error is shown in the plots; it measures how closely the solution satisfies the ordinary differential equations and boundary conditions. We also compared our results with data taken from [1] at
Pe=20
(blue circles), as can be seen in the last two snapshots. The concentration profile is in good agreement with [1], but our temperature profile is at best only in qualitative agreement.
Wolfram Cloud

You are using a browser not supported by the Wolfram Cloud

Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.


I understand and wish to continue anyway »

You are using a browser not supported by the Wolfram Cloud. Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.