No Title  |POLVER05_0
d(Fa)/d(V) = ra
d(Fb)/d(V) = rb
d(Fc)/d(V) = rc  -Rc
d(Fd)/d(V) = rd
d(Fcsg)/d(V) = Rc*Cc
d(y)/d(V) = -alpharho/2/y*(Ft/Fto)
Fosg = 0.1
vosg = 5
vsg = vosg*(Fosg+Fcsg)/Fosg
P = y*100
Ccsg = Fcsg/vsg
Ft = Fa+Fb+Fc+Fd
r1a = -k1a*Ca*Cb^2
r2c = -k2c*Ca^2*Cc^3
r1b = 2*r1a
r1c = -r1a
r2a = 2/3*r2c
r2d = -1/3*r2c
ra = r1a+r2a
rb = r1b
rc = r1c+r2c
rd = r2d
Ca = Cto*(Fa/Ft)*y
Cb = Cto*(Fb/Ft)*y
Cc = Cto*(Fc/Ft)*y
Cd = Cto*(Fd/Ft)*y
k1a = 1000
k2c = 60000
Fto = 20
alpharho = .0405
kcc = 2
Rc = kcc*(Cc -Ccsg)
Cto = 0.2
Scd = if(V>0.0001)then((Fc+Fcsg)/Fd)else(0)
V(0)=0
Fa(0)=10
Fb(0)=10
Fc(0)=0
Fd(0)=0
Fcsg(0)=0
y(0)=1
V(f)=50