Autothermal Reactor

This autothermal reactor consists of two concentric tubes. The feed enters an annular region between the tubes at temperature and flows from right to left. The reactants are at temperature when they enter the inner tube, which is a packed catalyst bed, where they flow left to right. An exothermic reaction in the catalyst bed raises the gas temperature and heat transfers from the catalyst bed to preheat the gas in the annular region. To solve the ordinary differential equations (ODEs) that model this reactor, a shooting method is used; a value of (the blue circle) is guessed and then changed until the correct value of is obtained. This is done by moving the slider for until the red circle is inside the green circle. Two reactor conditions are possible (two solutions to the mass and energy balances), one of which corresponds to essentially no reaction. Reactor startup determines which condition is obtained.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


Material balance on packed bed:
where and are the molar flow rate (mol/s) and concentration () of , is distance down the reactor (m), is the cross-sectional radium of the reactor (m), is the rate constant (1/s), is a pre-exponential factor (1/s), is activation energy (J/mol), is the ideal gas constant (J/[mol K]), is temperature (K), is volumetric flow rate (), is the inlet molar flow rate of (mol/s), is the ideal gas constant (), and is pressure (bar).
Energy balance on packed bed:
where is heat of reaction (J/mol), is heat transfer rate (J/s), is mass heat capacity (J/[g K]), is the overall heat transfer coefficient (), is the temperature of the annular region (K), is the feed temperature to the reactor (K), and is the feed temperature to the packed bed (K).
Energy balance on annular region:
The sign of the heat transfer term in the annular region is negative because integration is opposite to the direction of flow.
For the solution in which the reaction rate is nonzero, the feed temperature in the annular region increases as the feed flows to the left and is heated by the packed bed. The temperature increases in the packed bed due to the exothermic reaction and some of the energy transfers to the feed gas in the annular region. As the reactant is used up, the temperature in the packed bed reaches a maximum and decreases as heat is transferred to the feed gas.
The screencast video at [2] shows how to use this Demonstration.
[1] J. G. Ekerdt and J. B. Rawlings, Chemical Reactor Analysis and Design Fundamentals, Madison, WI: Nob Hill Publishing, 2002 pp. 326–332.


(University of Colorado Boulder, Department of Chemical and Biological Engineering)
"Autothermal Reactor"
 Wolfram Demonstrations Project
 Published: July 10, 2013
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2017 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+