No Title  |POLVER05_0 |1
Xe = ((thetaB+1)*Kc- (((thetaB+1)*Kc)^2-4*(Kc-4)*(Kc*thetaB))^0.5)/(2*(Kc-4)) 
d(Ta)/d(W) = Uarho*(T-Ta)/(mc*Cpcool) #
d(p)/d(W) = -alpha/2*(T/To)/p #
d(T)/d(W) = (Qg - Qr)/(Fao*sumcp) #
d(X) / d(W) = -ra/Fao 
X(0) = 0
Qg = (-ra)*(-Hr)#
Qr = Uarho*(T-Ta)#
alpha = .0002 #
To = 330 #
Uarho = 0.5 #
mc = 1000 #
Cpcool = 18 #
Hr = -20000 #
Fao = 5 #
thetaI = 1 #
CpI = 40 #
CpA = 20 #
thetaB = 1 #
CpB = 20 #
Cto = 0.3 #
Ea = 25000 #
Kc = 1000*(exp(Hr/1.987*(1/303-1/T))) #
k = .004*exp(Ea/1.987*(1/310-1/T)) #
yao = 1/(1+thetaB+thetaI) #
Cao = yao*Cto #
sumcp = (thetaI*CpI+CpA+thetaB*CpB) #
Ca = Cao*(1-X)*p*To/T #
Cb = Cao*(thetaB-X)*p*To/T #
Cc = Cao*2*X*p*To/T #
ra = -k*(Ca*Cb-Cc^2/Kc) #
W(0)=0
Ta(0)=320
p(0)=1
T(0)=330
W(f) = 4500