#Example RE13-3 Proportional-Integral Control of a CSTR |POLVER05_0 |1 d(Ca)/d(t) = 1/tau*(Ca0-Ca)+ra # d(Cb)/d(t) = 1/tau*(Cb0-Cb)+ra # d(Cc)/d(t) = 1/tau*(0-Cc)-ra # d(Cm)/d(t) = 1/tau*(Cm0-Cm) # d(T)/d(t) = (-Q-Fa0*ThetaCp*(T-T0)+(-36000)*ra*V)/NCp # d(I)/d(t) = T-Tsp # Fa0 = 80 # T0 = 70 # V = (1/7.484)*500 # Tsp = 138 # UA = 16000 # Ta1 = 60 # kc = 8.5 # k = 16.96e12*exp(-32400/1.987/(T+460)) # Fb0 = 1000 # Fm0 = 100 # mc0 = 1000 # ra = -k*Ca # NCp = Ca*V*35+Cb*V*18+Cc*V*46+Cm*V*19.5 # ThetaCp = 35+Fb0/Fa0*18+Fm0/Fa0*19.5 # v0 = Fa0/0.923+Fb0/3.45+Fm0/1.54 # Ca0 = Fa0/v0 # Cb0 = Fb0/v0 # Cm0 = Fm0/v0 # tau = V/v0 # X = (Ca0-Ca)/Ca0 # mc = mc0+kc*(T-Tsp)+kc/tau*I # Ta2 = T-(T-Ta1)*exp(-UA/(18*mc)) # Q = mc*18*(Ta2-Ta1) # t(0)=0 Ca(0)=0.03789 Cb(0)=2.12 Cc(0)=0.143 Cm(0)=0.2265 T(0)=138.53 I(0)=0 t(f)=4