Matematické Fórum
Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.
Nástěnka
❗22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
❗04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
❗23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.
Nejste přihlášen(a). Přihlásit
#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: 4939
- 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