No Title  |POLVER05_0 |1
# Example 12-7 part c
d(Fa)/d(V) =  ra  #
d(Fb)/d(V) = rb #
d(Fc)/d(V) = rc #
d(Fd)/d(V) = rd #
d(T)/d(V) = (Qg-Qr)/sumFiCpi #
d(Ta) / d(V) = Ua*(T-Ta)/m/Cpco*0
ra = r1a+r2a #
rb = r1b #
rc = r1c+r2c #
rd = r2d #
r1a = -k1a*Ca*Cb^2 #
r2c = -k2c*Ca^2*Cc^3 #
k1a = 40*exp ((E1/R)*(1/300-1/T)) #
Ca = Cto*(Fa/Ft)*(To/T)*p #l
Cb = Cto*(Fb/Ft)*(To/T)*p #
Cc = Cto*(Fc/Ft)*(To/T)*p #
k2c = 2*exp((E2/R)*(1/300-1/T)) #
Cpco=10
m = 50 
Ta(0) = 325
r1b = 2*r1a #
r1c = -r1a #
r2a = 2/3*r2c #
r2d = -1/3*r2c #
E1 = 8000 #
E2 = 12000 #
DH1b = -15000 #
DH2a = -10000 #
p = 1 #
sumFiCpi = Cpa*Fa+Cpb*Fb+Cpc*Fc+Cpd*Fd #
Cpa = 10 #
Cpb = 12 #
Cpc = 14 #
Cpd = 16 #
To = 300 #
Ua = 80 #
Cto = 0.2 #
Ft = Fa+Fb+Fc+Fd #
R = 1.987 #
Qg = r1b*DH1b+r2a*DH2a #
Qr = Ua*(T-Ta) #
V(0)=0
Fa(0)=5
Fb(0)=10
Fc(0)=0
Fd(0)=0
T(0)=300
V(f)=10
SCD=if(V>0.0001) then (Fc/Fd) else (0)