Skip to main content

Mole balance for a Plug-Flow reactor

In principle any balance, be it in baking a cake, accounting, electrical engineering, mechanical engineering or chemical engineering is suppose to be simple, but for some reason maths always find a way to complicate things a little. So lets get the balance going.

The general balance for any thing in life is as follows:
EQ1: General Balance equation
Before trying to visualize the math, lets first try and understand the hardware we are trying to describe mathematically. This is what a PFR (Plug Flow Reactor) looks like.

FIG1: Plug Flow Reactor
And this is how we are going to graphically illustrate a plug flow reactor.

FIG2: Mole balance of species i in dv
with these images in mind, lets try and explain the mathematics that describes it, by taking the mole balance and substituting values that is of importance to the system we obtain:

EQ2: General Mole Balance
with Fi0 and Fi being the initial and final feed rates respectively in [moles/time], Gi being gain [mole/time] and the last term is your accumulation term [moles/time]. We must firstly consider the Gi term because it is a bit of a arbitrary concept, especially for the mole balance on a PFR.

Gi or gain is a function of ri, the rate of reaction, this term is obtained experimentally Gi is defined as :
EQ3: Gain 
 with V the reactor volume and ri, rate of reaction, defined as:
EQ4: Rate of reaction 
We will not really be using EQ4 for the explanation of the mole balance on the PFR because dCi can be of zeroth, first, second etc order and thus changes for various systems.

Substitution of EQ3 into EQ2 we get:
EQ5
the above equation is still just a general mole balance, to make it applicable to our PFR we need to make some assumptions that is suitable for using when designing a PFR.

The assumptions are:

  • There is no radial variation within the volume 
  • Steady state conditions
  • The flow inside the pipe to be Plug-Flow
After the assumptions EQ5 becomes;
EQ6
Lets now consider a small enough section of pipe, i.e. a differential volume (the orange part of FIG2), so that ri can be assumed constant, the gain term, EQ3, can then be described as follows;
EQ7

Doing a mole balance on this differential section of the pipe, i.e. substituting EQ7 into EQ6, and taking the limit of that equation we get;

EQ8
EQ8 is equivalent to the definition of the derivative, this yields:
EQ9: Differential  form of  a steady state mole balance of a PFR
Integration yields:
EQ10: Integral form of PFR design  equation

The equation can be written in terms of conversion as well, and can easily be done for a batch reactor and a Continuous Stirred Tank Reactor, CSTR.

Comments

Popular posts from this blog

How to plot an XY diagram for McCabe Thiele method in Aspen Plus

Its been a really long time since I made any posts, sorry for that but work sometimes gets in the way of the things I like doing, like updating my blog. This post might seem a bit stupid but I tried to find out how to do this on the web and it was quite impossible to find a easy step by step guide for making a XY plot in Aspen Plus for use in the McCabe Thiele method for designing a distillation coulomb, so I had to figure it out myself after a lot of fellow class mates had the same problem as me. So without any further a dew here is my 13 steps to making a XY plot in Aspen Plus for a binary mixture. For this example I used the old classic ethanol and water problem. Step1:  Open Aspen Plus and this is the first screen to pop up. Choose "Blank Simulation" an press "Ok" Step 1: Begin new Aspen simulation Step2:  You are now in the simulation environment and to start the simulation you need to enter the "Data Browser" (Usual simulations you would

Designing a C5/C6 Isomerization unit

As a final year student we need to design a catalytic isomerization unit that effectively rises the octane number by isomerization of n-butane, n-pentane and n-hexane into their respective isomers. I was responsible for the design of the reactor, this was probably the most difficult thing I have done in engineering since I started the course. It took a lot of hard work, but in return I am most certainly a much better engineer.... (P.S. The most difficult part of the project was finding the correct reaction kinetics, and just a note to the reader, I do not fully trust the kinetics that was used in this design.) 3.1        Reactions and Kinetics n-Butane In a process such as the Penex process where carbon strings from C4-C6 are present a large amount of possible reactions can take place in the form of isomerisation and cracking. We will begin by considering the short straight chains and their reaction and kinetics first then move to the cyclic molecules and their reactions. Con

What is fugacity really?

Fugacity is supposed to be easily understood, especially because that was the main purpose of its invention. Lets just get some history straight first. In the beginning there was chemical potential, usually denoted as  μ, and defined as: Definition of Chemical Potential with G and n i being the Gibbs free energy and moles of compound i respectively. As you see  μ is used to describe the Gibbs free energy, and if you know thermodynamics you should know that Gibbs is used to describe chemical equilibrium in pure and mixed systems, if you didn't know, this is very important. The problem however with chemical potential is that it is an abstract concept, like the amount of honest bankers, in other words, it can not be directly measured. In its stead fugacity was invented by one, G.N. Lewis (Also famous for his Lewis dot structure of molecules and his mustache). G.N. Lewis Partial fugacity, is the pseudo-pressure exerted by the mixture of molecules in the vapor phase, and it i