Optimizing Reactor Sizing for Gas-Phase Irreversible Reaction

In summary, the maximum conversion that can be achieved with these two reactors connected in series is 0.445, and the overall conversion is 0.595 if two 400-liter CSTRs are connected in series. The overall conversion is 0.52 if two 400-liter CSTRs are connected in parallel with half of the feed going to each reactor.
  • #1
gfd43tg
Gold Member
950
50

Homework Statement


A 400-liter CSTR and a 100-liter PFR are available to process 1.0 liter of feed per second. The feed contains 41% A, 41% B, and 18% inerts. The irreversible gas-phase reaction

[tex] A + B → C [/tex]

is to be carried out at 10 atm and 227 °C. The rate of reaction in mol/liter min is given below as a function of conversion

[tex]\begin{pmatrix}
-r_{A} & 0.2 & 0.0167 & 0.00488 & 0.00286 & 0.00204 \\
\hline
X & 0.0 & 0.1 & 0.4 & 0.7 & 0.9
\end{pmatrix}
[/tex]

(a) What is the maximum conversion that can be achieved with these two reactors connected in series (Ans. ##X_{C,P}## = 0.445, ##X_{P,C}## = 0.515.)

(b) What would be the overall conversion if two 400-liter CSTRs were connected in series for the same feed and operating conditions? (Ans: X = 0.595)

(c) What would be the overall conversion if two 400-liter CSTRs were connected in parallel with half of the feed going to each reactor? (Ans: X = 0.52)

(d) What is the volume of a single tubular reactor necessary to achieve 60% conversion if the molar feed rate is 2 mol A/min? (Ans: V = 180 liters.)

(e) If the total pressure were reduced by a factor of 10, would the conversion increase, decrease, or remain the same?

(f) Give a critique of the answers to this problem.


Homework Equations


The Attempt at a Solution


Right now I am working on part (a), and have some troubles. First off, is it fair to say that if I use a basis of one second, then the initial concentration of A, which I am using as my limiting reactant, is 0.41 mol/1 L, for 0.41 M?

In my attempt, I start off with a CSTR connected to a PFR to see the overall conversion, which I will then do a PFR connected to a CSTR and compare the overall conversion.

My material balance on the CSTR
[tex] \frac {dN_{A}}{dt} = F_{A0} - F_{A1} + r_{A1}V_{1}[/tex]

Knowing that ##F_{A1} = F_{A0} - F_{A0}X_{1}##, and assuming steady state, the mole balance simplifies to

[tex] 0 = F_{A0} - F_{A0} + F_{A0}X_{1} + r_{A1}V_{1}[/tex]

[tex]\frac {-F_{A0}X_{1}}{V_{1}} = r_{A1}[/tex]

[tex] \frac {-0.41X_{1}}{400} = r_{A1}[/tex]

Here I run into a problem. I don't know what the conversion of the CSTR is, hence I don't know what the reaction rate is, so I have two unknowns and only one equation.

I put the data into matlab, and here is my code
Code:
rA = [-0.2 -0.0167 -0.00488 -0.00286 -0.00204];
X = [0 0.1 0.4 0.7 0.9];
FA0 = 0.41;
plot(X,-FA0./(rA/60))
xlabel('X')
ylabel('-FA0/rA')
title('FA0/-rA vs. Conversion')

This reaction appears to be zero order, based on the nearly straight, which implies a zero order reaction. This may be a clue, but I'm not sure.
ImageUploadedByPhysics Forums1410785249.625731.jpg
 

Attachments

  • 2-10b.jpg
    2-10b.jpg
    12.7 KB · Views: 951
Last edited:
Physics news on Phys.org
  • #2
The reaction rate is a function of conversion r = r(X), so you have one equation in one unknown for the conversion. The equation is non-linear, and you have to interpolate in your table to solve for the conversion in the CSTR. This non-linear equation can be solved in a variety of ways, like secant method or half-interval technique.

chet
 
  • #3
I did a polynomial fit to the data given (finding ##r_{A} = f(X)## is proving not to be so trivial),

Code:
polyfit(X,rA/60,2)

ans =

   -0.0081    0.0097   -0.0024

Hence my ##r_{A} = -0.0081X^2 + 0.0097X - 0.0024##

Then I equate it ##\frac {-0.41X_{1}}{400} = -.0081X_{1}^{2} + 0.0097X_{1} - 0.0024## and find the roots,

Code:
 roots([-.0081 0.010725 -0.0024])

ans =

    1.0389
    0.2852

So I take it that ##X_{1} = 0.285##, well that solves the conversion of the CSTR, but I still need to find the conversion of the PFR, ##X_{2}##.

Is it that case that ##r_{A2}## would not be the same function as ##r_{A1} = f(X)##, or will they be the same function of ##X##?

I then try to find the conversion of the PFR using the integral form to find the volume of a PFR

[tex] V = F_{A0}\int_{X_{in}}^{X_{out}} \frac {dX}{-r_{A}}[/tex]

keep in mind that when I used matlab, I plotted ##X## against ##r_{A}##, and not against ##-r_{A}##, which is why I omit the negative sign in the integral when I do it below.

[tex] 0.41 \int_{0.285}^{X_{2}} \frac {dX}{-.0081X^2 + 0.0097X - .0024} = 100 [/tex]

Then I run with guess and check in matlab

Code:
 f = @(x) 1./(-.0081*x.^2 + .0097*x - .0024)

f = 

    @(x)1./(-.0081*x.^2+.0097*x-.0024)

0.41*quad(f,.285,.445)

ans =

 -107.3562

0.41*quad(f,.285,.446)

ans =

  351.9439

But this is obviously behaving wildly
 
Last edited:
  • #4
It looks like you have the right idea. How does your polynomial fit to the data look if you compare it vs the data. My guess is that, when you make the plot, you will not be happy with what you see, since the function r(X) is monotonic, and your fit is not going to be monotonic. A polynomial fit is going to turn out to be not such a good choice. Try plotting log r vs X and see how it looks. Then fit a polynomial to logr (or lnr) vs X. I think you will be much happier with the fit. Then re-solve the CSTR case and see what you get. In the case of the PFR, might I suggest doing the integration numerically, rather than trying to get an analytic solution. Using an automatic integrator (if you have access to one) might things easier, especially one that will automatically stop the integration when you satisfy the final condition on the total volume.

Chet
 
  • #5
Okay, I will try that. The polynomial fit is not that great. By the way, when I am using the function quad in MATLAB, that is a numeric integrator, I guess it has an input function. I can just use trapz to integrate the data only though.

This problem is really making me dig deep into using matlab, got to figure out how to do curve fitting beyond a polynomial fit, awesome!

EDIT: Wow, after spending all this time with matlab, I realized that Excel can do this a million times more easily, and I got a logarithmic curve fit,

##r_{A} = 0.0001 \ln{X} + 1*10^{-5}##

I solve and find ##X_{1} = 0.1656##.

I can see why you say to numerically integrate. When I try to integrate

[tex] 0.41 \int_{0.1656}^{X_{2}} \frac {dX}{0.0001 \ln{X} + 1*10^{-5}} [/tex],

and
Code:
fhan = @(x) 1./(0.0001*log(x) + 1e-5); 
0.41*quad(fhan,0.1656,0.4)

ans =

 -841.2465
I get extremely erroneous results.

The plot
Code:
x = linspace(0,1,1000);
plot(x,fhan(x))

looks completely different from
Code:
plot(X,fhan(X))

This is so weird, when I plot ##r_{A}## vs. ##X##, it has a logarithmic fit, but when I do ##1/r_{A}## vs. ##X##, it has a linear fit.
 
Last edited:
  • #6
This has become quite irritating. I believe I found a better fit than a logarithm, and this is my equation

[tex] r_{A} = \frac{1}{-31412X - 122.76}[/tex]

When I fit this, it actually fits pretty nicely. I uploaded my plots in the PDF as well as code.

This is how I calculated the conversion through the CSTR, ##X_{1}##.
Code:
fZero = 

    @(x)(1./(-31412*x-122.76))+0.001025*x

EDU>> fzero(fZero,0.5)

ans =

    0.1743

However, when I try to find the area

Code:
fThis

fThis = 

    @(x)-31412*x-122.76

EDU>> FA0*abs(quad(fThis,0.1743,0.445))

ans =

   1.0932e+03
I am off by a factor of 10 for the volume, so I can't explain what is going on!
 

Attachments

  • curvefit.pdf
    135.1 KB · Views: 417
  • #7
Maylis said:
This has become quite irritating. I believe I found a better fit than a logarithm, and this is my equation

[tex] r_{A} = \frac{1}{-31412X - 122.76}[/tex]

When I fit this, it actually fits pretty nicely. I uploaded my plots in the PDF as well as code.

This is how I calculated the conversion through the CSTR, ##X_{1}##.
Code:
fZero = 

    @(x)(1./(-31412*x-122.76))+0.001025*x

EDU>> fzero(fZero,0.5)

ans =

    0.1743

However, when I try to find the area

Code:
fThis

fThis = 

    @(x)-31412*x-122.76

EDU>> FA0*abs(quad(fThis,0.1743,0.445))

ans =

   1.0932e+03
I am off by a factor of 10 for the volume, so I can't explain what is going on!
None of the functionalities you showed look like they fit the data. Why don't you just do the problem with linear interpolation between the data points and see what you get? It's got to be better than what you have been getting so far. This will at least give you a start. You can then try to develop a better functionality later.

Chet
 
  • #8
Try this:

For 0≤x≤0.1, ##r=-0.2e^{-24.829X}##

For 0.1≤X≤0.4, ##r=-0.02517e^{-4.101X}##

For 0.4≤X≤0.7, ##r=-0.0.00995e^{-1.781X}##

For 0.7≤X≤0.9, ##r=-0.00933e^{-1.6895X}##

Chet
 
  • #9
I don't understand how linearly interpolating between the points will help me solve for the conversion. Don't I need to simultaneously know ##r_{A}##?
 
  • #10
Maylis said:
I don't understand how linearly interpolating between the points will help me solve for the conversion. Don't I need to simultaneously know ##r_{A}##?
Interpolating r vs X between the data points is conceptually the same as fitting your polynomial or your other functions to the r vs X data. At least the interpolation is monotonic and passes through all the data points. I prefer the interpolation I provided in my most recent post, involving logarithmic/exponential interpolation. Each segment is even curved in the correct way.

For the CSTR, you are just solving 2 equations in two unknowns (X1 and r(X1)), correct?

I ran the calculation for a CSTR, and got a conversion of about 0.2, indicating that the solution lies in the second interpolation interval.

chet
 

Related to Optimizing Reactor Sizing for Gas-Phase Irreversible Reaction

What is conversion and why is it important in reactor sizing?

Conversion is the fraction of reactant that is converted into product in a chemical reaction. It is important in reactor sizing because it determines the amount of reactant that needs to be fed into the reactor and the amount of product that will be produced.

How is conversion calculated?

Conversion is calculated by dividing the moles of reactant that have been converted into product by the initial moles of reactant, and multiplying by 100% to get a percentage. Alternatively, it can also be calculated by dividing the difference between the initial moles of reactant and the moles of reactant present after the reaction by the initial moles of reactant, and multiplying by 100%.

What factors affect conversion in a chemical reaction?

The factors that affect conversion in a chemical reaction are the reaction kinetics, temperature, pressure, concentration of reactants, and the presence of catalysts or inhibitors.

What is a reactor and what are the different types of reactors?

A reactor is a vessel or tank in which a chemical reaction takes place. There are three main types of reactors: batch reactors, continuous flow reactors, and semi-batch reactors. Batch reactors operate by adding all the reactants at once and allowing the reaction to take place until completion. Continuous flow reactors continuously feed reactants into the reactor and continuously remove products. Semi-batch reactors are a combination of batch and continuous flow reactors, where one or more reactants are fed into the reactor continuously while another reactant is added in batches.

How is reactor sizing determined?

Reactor sizing is determined by considering factors such as the desired conversion, reaction kinetics, and reaction rate. It also takes into account the type of reactor being used and the operating conditions, such as temperature and pressure. Reactor sizing involves mathematical calculations and can also be determined through experiments or simulations.

Similar threads

  • Engineering and Comp Sci Homework Help
Replies
1
Views
2K
  • Engineering and Comp Sci Homework Help
Replies
3
Views
1K
  • Materials and Chemical Engineering
Replies
10
Views
36K
Back
Top