Ontola > diskuze
| Nahlásit

CSTR reaktor - ustálený stav

Ahoj, můžu se zeptat, kde je chyba ve vykreslení a zjisštění ustáleného stavu v CSTR reaktoru ? Program MATLAB.
clear all;clc;
[t,x]=ode45(@soustava,[0 500],[0 360]);
plot(t,x(1))
function [pr]=soustava(t,x)
V=0.001;
dH=-83740;
A=0.025;
alfa=70;
T0=273;
cp=2720000;
ca0=6670;
ro=1000;
F=1.67e-6;
R=8.314;
T0=290;
Tc=360;
pr=zeros(2,1);
k=0.55*10^6*exp(-10000/x(1));
r=k*6670*6670*(1-x(2))^2;
pr(1)=(1.67*10^-6*(290-x(1))/0.001)+(0.025*70*(360-x(1))/(0.001*2720000))-r*(-83740)/(2720000);
pr(2)=1.67*10^-6*x(2)/0.001-k*6670^2*(1-x(2));
end
Témata: Nezařazené

Reakce

 Anonym
Odpovídat lze i bez registrace. Dodržujte pravidla Ontoly
Vložit: Obrázek