#1 19. 01. 2021 16:16
CSTR - ustálené stavy
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
Offline
#2 19. 01. 2021 18:06
- check_drummer
- Příspěvky: 5578
- Reputace: 106
Re: CSTR - ustálené stavy
Ahoj, proč myslíš, že je tam chyba?
"Máte úhel beta." "No to nemám."
Offline
#3 19. 01. 2021 18:14
Re: CSTR - ustálené stavy
↑ check_drummer:
protože mi mají vyjít 3 ustálené stavy při teplotách Ts1 = 314 K, Ts2 = 451K a Ts3 = 336 K
ale když dám plot, tak se mi tam žádný graf nezobrazuje. Bilance hmoty a energie (pr(2), pr(1)) mám snad dobře vyjádřené.
Offline