No Title |POLVER05_0 |1
#Adiabatic semibatch reactor
#mole balances and initial conditions
d(NA)/d(t)= rA*V # mol/min
NA(0)= 500 # reactor contains only A initially, 500 mol
d(NB)/d(t)= FBi + rA*V # B is added in feed stream with flow rate FBi
NB(0)= 0
d(NC)/d(t)= - rA*V
NC(0)= 0
d(V)/d(t)= vo #reactor volume increases because B is added
V(0) = 50 #initial reactor volume (L)
rA= - k*cA*cB
t(0)=0
t(f)=120 # stop integration after 120 min
#Energy balance and initial conditions
d(temp)/d(t)=(FBi*CPB*(tempin-temp)+delH*rA*V)/(NA*CPA+NB*CPB+NC*CPC)
temp(0)=298 # initial reactor temperature (K)
#Parameters
CPA=15 #heat capacity (cal/(mol K))
CPB=15
CPC=30
delH=-6000 #heat of reaction (cal/mol)
R=1.98 #gas constant (cal/(mol K))
tempin=273+50 # feed temperature of B (K)
vo=1.0 #volumetric feed rate of B (L/min)
NAo=500 # Initial number of moles of A in reactor (mol)
#Other equations
#molar flow rate (mol/min) of B into reactor; flow stops after 50 min
FBi= if (t<50) then (10) else (0)
#rate constant (L/(mol min); activation energy in cal/mol
k=1.7*10^4*exp(-10000/R/temp)
cA=NA/V # concentration of A in reactor
cB=NB/V # concentration of B in reactor
X=(NAo-NA)/NAo # fractional conversion of A