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

gag
Příspěvky: 44
Reputace:   
 

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: 3275
Reputace:   90 
 

Re: CSTR - ustálené stavy

Ahoj, proč myslíš, že je tam chyba?


Popelka - pohádka o neprosté funkci nabývající minima v jediném bodě

Offline

 

#3 19. 01. 2021 18:14

gag
Příspěvky: 44
Reputace:   
 

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

 

Zápatí

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson