Example RE13-1 Falling Off the Upper Steady State |POLVER05_0 |1 d(Ca)/d(t) = 1/tau*(Ca0-Ca)+ra # d(Cb)/d(t) = 1/tau*(Cb0-Cb)+rb # d(Cc)/d(t) = 1/tau*(0-Cc)+rc # d(Cm)/d(t) = 1/tau*(Cm0-Cm) # d(T)/d(t) = (Q-Fa0*ThetaCp*(T-T0)+(-36000)*ra*V)/NCp # Fa0 = 80 # T0 = 70 # V = (1/7.484)*500 # UA = 16000 # Ta1 = 60 # k = 16.96e12*exp(-32400/1.987/(T+460)) # Fb0 = 1000 # Fm0 = 100 # mc = 1000 # ra = -k*Ca # rb = -k*Ca # rc = k*Ca # Nm = Cm*V # Na = Ca*V # Nb = Cb*V # Nc = Cc*V # ThetaCp = 35+Fb0/Fa0*18+Fm0/Fa0*19.5 # v0 = Fa0/0.923+Fb0/3.45+Fm0/1.54 # Ta2 = T-(T-Ta1)*exp(-UA/(18*mc)) # Ca0 = Fa0/v0 # Cb0 = Fb0/v0 # Cm0 = Fm0/v0 # Q = mc*18*(Ta1-Ta2) # tau = V/v0 # NCp = Na*35+Nb*18+Nc*46+Nm*19.5 # X = 1-Ca/(80/v0) # t(0)=0 Ca(0)=0.03789 Cb(0)=2.12 Cc(0)=0.143 Cm(0)=0.2265 T(0)=138.53 t(f)=4