#Startup of a CSTR 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) = (Qg-Qr)/NCp # Fa0 = 36.3 # T0 = 297 # V = 1.89 # UA = 7262 # dh=-20013 Ta1 = 288.7 # k = 16.96e12*exp(-18012/1.987/(T)) # Fb0 = 453.6 # Fm0 = 45.4 # mc = 453.6# 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/14.8+Fb0/55.3+Fm0/24.7 # Ta2 = T-(T-Ta1)*exp(-UA/(18*mc)) # Ca0 = Fa0/v0 # Cb0 = Fb0/v0 # Cm0 = Fm0/v0 # Qr2 = mc*18*(Ta2-Ta1) # tau = V/v0 # NCp = Na*35+Nb*18+Nc*46+Nm*19.5 # Qr1=Fa0*ThetaCp*(T-T0)# Qr=Qr1+Qr2 Qg=ra*V*dh t(0)=0 Ca(0)=0 Cb(0)=55.3 Cc(0)=0 Cm(0)=0 T(0)=297 t(f)=4