Tubular Reactor: Mass Balance

Download as pdf or txt
Download as pdf or txt
You are on page 1of 21

Tubular Reactor

Mass balance:
Reaction

Diffusion Advection/Convection

Boundary conditions:

(Danckwerts)
Numerical Methods for Chemical Engineers 1
Tubular Reactor

Mass balance:

Non-dimensionalization:

Numerical Methods for Chemical Engineers 2


Tubular Reactor

Mass balance:

Non-dimensionalization:

Numerical Methods for Chemical Engineers 3


Tubular Reactor

Mass balance:

Peclet number:

Damkhler number:
Numerical Methods for Chemical Engineers 4
Tubular Reactor

Mass balance:

Steady state tubular reactor

Numerical Methods for Chemical Engineers 5


Tubular Reactor

Discretize in space

Numerical Methods for Chemical Engineers 6


Tubular Reactor

Discretize in space

Numerical Methods for Chemical Engineers 7


Tubular Reactor
Centered Finite Difference Backward Finite Difference

Discretize in space

Numerical Methods for Chemical Engineers 8


Tubular Reactor

Boundary conditions!

Discretize in space

Numerical Methods for Chemical Engineers 9


Tubular Reactor

Discretize in space

Numerical Methods for Chemical Engineers 10


Tubular Reactor

Ghost point Ghost point


Discretize in space

Numerical Methods for Chemical Engineers 11


Tubular Reactor

System of nonlinear equations!!!

Numerical Methods for Chemical Engineers 12


Assignment 1
1. Solve the steady state tubular reactor for 20 different Peclet
numbers (between 0.01 and 100) and for a first (n=1) and a second
(n=2) order reaction. Use a Damkhler number of unity.
Complete the template rhs.m by implementing the non-linear
equations to solve
2. Plot the conversion at the end of the reactor vs. the
Peclet number for both reaction orders
Also plot the ratio between the conversions of the first order and
second order reaction
What is better for these reactions, a lot of back-mixing (Pe small,

CSTR) or ideal plug flow (Pe large, PFR)?


What influence does the reaction order have overall and at low or
high Peclet numbers?
Complete the template TubReact_steady_state.m
Numerical Methods for Chemical Engineers 13
Tubular Reactor

Mass balance:
Reaction

Diffusion Advection/Convection

Boundary conditions:
Numerical Methods for Chemical Engineers 14
Tubular Reactor

Mass balance:

Dynamic tubular reactor


Boundary conditions:

Numerical Methods for Chemical Engineers 15


Tubular Reactor

Discretize in space

Numerical Methods for Chemical Engineers 16


Tubular Reactor

ODEs

Discretize in space

Numerical Methods for Chemical Engineers 17


Tubular Reactor
ODEs

Ghost point Ghost point


Discretize in space

Numerical Methods for Chemical Engineers 18


Time index Tubular Reactor Order of reaction

E.g. Explicit Euler, But in general ver stiff

Discretize in time

Numerical Methods for Chemical Engineers 19


Tubular Reactor

System of nonlinear ODEs!!!


Stiff...
Numerical Methods for Chemical Engineers 20
Assignment 2
1.Solve the dynamic tubular reactor from initial 0
to final time of 5 with MATLABs ode23s
Use the rhs.m from assignment 1 and the
template TubReact_dynamic.m
Consider only a first order reaction with Pe=100
and Da=1
2.Plot the conversion at the end of the reactor vs.
dimensionless time
3.At what time does the solution reach a steady
state, i.e. how many reactor volumes of solvent
will you need?
Numerical Methods for Chemical Engineers 21

You might also like