Scribd
Upload
16 views6 pages

Cre Assignment 2

The document demonstrates the use of a recycle reactor with a focus on the van de Vusse reaction kinetics. It includes Python code for plotting concentration trajectories and determining optimal parameters for maximizing the concentration of product B. The optimal recycle ratio and residence time are identified as 0.68 and 0.3, respectively, with further improvements possible using a CSTR in series with a PFR.

Uploaded by

pratik Kewalramani
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
16 views6 pages

Cre Assignment 2

The document demonstrates the use of a recycle reactor with a focus on the van de Vusse reaction kinetics. It includes Python code for plotting concentration trajectories and determining optimal parameters for maximizing the concentration of product B. The optimal recycle ratio and residence time are identified as 0.68 and 0.3, respectively, with further improvements possible using a CSTR in series with a PFR.

Uploaded by

pratik Kewalramani
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

Recycle reactor

January 19, 2025

1 PFR trajectory convex hull demonstration


1.1 Import necessary python packages
[2]: # artools
import sys
sys.path.append('../../../artools')
from importlib import reload
import artools
artools = reload(artools)
import numpy as np
# anaconda
from ipywidgets import interact

import scipy as sp
from scipy.spatial import ConvexHull

import matplotlib.pyplot as plt


from matplotlib.tri import Triangulation
%matplotlib inline

plt.style.use("ggplot")
from mpl_toolkits.mplot3d import Axes3D

1.2 Define kinetics


We use the van de vusse system to demonstrate:

A→B→C

2A → D

[84]: # 2D van de Vusse kinetics


# A -> B -> C

def rate_fn(C,t):
cA = C[0]
cB = C[1]

1
#rate constants
k1 = 1.0
k2 = 1.0
k3 = 10.0

#r = [rA, rB]
return np.array([-k1*cA - 2*k3*cA**2,
k1*cA - k2*cB])

[320]: def plot_fn(L=0.85, pfr_tend=0.21, N=25):


Cf = np.array([1.0, 0.0])

pfr_cs, pfr_ts = artools.pfrTrajectory(Cf, rate_fn, 10)


cstr_cs, cstr_ts = artools.cstrLocus_fast(Cf, rate_fn, 1000, 200)

fig = plt.figure(figsize=(5,5))
ax = fig.gca()

ax.plot(cstr_cs[:, 0], cstr_cs[:, 1], "bx")


ax.plot(pfr_cs[:, 0], pfr_cs[:, 1], "r-")

pfr_cs, pfr_ts = artools.pfrTrajectory(Cf, rate_fn, pfr_tend)


ax.plot(pfr_cs[:, 0], pfr_cs[:, 1], "k-")

for i in range(N):
c_mix = L*pfr_cs[-1, :] + (1.0 - L)*Cf
ax.plot(c_mix[0],c_mix[1],"*")
pfr_cs, pfr_ts = artools.pfrTrajectory(c_mix, rate_fn, pfr_tend)

ax.plot(pfr_cs[:, 0], pfr_cs[:, 1], "k-")

ax.plot(pfr_cs[-1, 0], pfr_cs[-1, 1], "ko")

#plt.show(fig)

index=np.linspace(0,len(cstr_cs)-1,10)
for j in index:
j=int(j)

c_cstr_output=cstr_cs[j,:]

pfr_cs_from_cstr, pfr_ts = artools.pfrTrajectory(c_cstr_output, rate_fn,␣


15)
,→

ax.plot(pfr_cs_from_cstr[:, 0], pfr_cs_from_cstr[:, 1])


plt.text(c_mix[0],c_mix[1],"starting for recycle *")

2
ax.set_xlabel("C_A")
ax.set_ylabel("C_B")

plt.show(fig)
interact(plot_fn, L=(0,1,0.01), pfr_tend=(0,10,0.05), N=(1,25,1))
plot_fn(L=0.68,pfr_tend=0.3)
print("The maximum conc of B is seen by hit and trial of parameters, the␣
,→parameters found potimal given that we can only use recycle reactor are given␣

,→by Recycle ratio=0.68, residance time=0.3, with maximum value of cb ~~ 0.095.␣

,→other wise more better configuration can be found by using cstr in series with␣

,→pfr as shown in figure.")

interactive(children=(FloatSlider(value=0.85, description='L', max=1.0, step=0.


,→01), FloatSlider(value=0.21, de...

The maximum conc of B is seen by hit and trial of parameters, the parameters
found potimal given that we can only use recycle reactor are given by Recycle
ratio=0.68, residance time=0.3, with maximum value of cb ~~ 0.095. other wise
more better configuration can be found by using cstr in series with pfr as shown
in figure.

3
[315]: def feild_plot(c,t):
for k in [0.1,1]:

ca=np.linspace(0.001,k,50)

cb=np.linspace(0.001,0.1,10)
fig = plt.figure(figsize=(10,5))
ax = fig.gca()
for i in ca:
for j in cb:
normal_rate=rate_fn([i,j],t)
normal_rate=normal_rate/np.linalg.norm(normal_rate)
ax.quiver(i, j, normal_rate[0],␣
,→normal_rate[1],scale=90,color="b")

plt.show()

[316]: feild_plot(1,1)

4
[ ]:

[ ]:

OPTIMAL PLOT

5
6

You might also like