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)