05.05.2013 Views

HEMISKI Zbirka re{eni zada~i REAKTORI 3

HEMISKI Zbirka re{eni zada~i REAKTORI 3

HEMISKI Zbirka re{eni zada~i REAKTORI 3

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Petti del. Kataliti~ki <strong>re</strong>aktor so fiksen sloj katalizator (PBR)<br />

Programata vo POLYMATH gi vklu~uva ravenkata za<br />

dizajn (4) i toplinskiot bilans (5) zaedno so brzinskiot izraz<br />

(2) i izrazite za molskite koncentracii (1). Isto taka, vo<br />

programata ja pi{uvame i kvadratnata ravenka (6). Se dodavaat<br />

i izrazite za temperaturnite zavisnosti na konstantite. Se<br />

dobivaat slednive <strong>re</strong>zultati (programa, <strong>re</strong>zultati i grafi~ki<br />

prikaz):<br />

POLYMATH Results<br />

Calculated values of the DEQ variables<br />

Variable initial value minimal value maximal value final value<br />

W 0 0 15 15<br />

X 0 0 0.7370304 0.7370304<br />

CAo 0.271 0.271 0.271 0.271<br />

FAo 5 5 5 5<br />

To 450 450 450 450<br />

T 450 450 1187.0304 1187.0304<br />

k 0.5010696 0.5010696 515.51332 515.51332<br />

K 2.5E+04 32.759106 2.5E+04 32.759106<br />

CB 0 0 0.0599532 0.0599532<br />

CA 0.271 0.0427821 0.271 0.0427821<br />

R 0.0367991 9.705E-05 1.2217073 9.705E-05<br />

A 6.098E+06 7989.9459 6.098E+06 7989.9459<br />

B 1.22E+07 1.717E+04 1.22E+07 1.717E+04<br />

Xe 0.9939255 0.737043 0.9939255 0.737043<br />

ODE Report (RKF45)<br />

Diffe<strong>re</strong>ntial equations as ente<strong>re</strong>d by the user<br />

[1] d(X)/d(W) = R/FAo<br />

Explicit equations as ente<strong>re</strong>d by the user<br />

[1] CAo = 0.271<br />

[2] FAo = 5<br />

[3] To = 450<br />

[4] T = To+1000*X<br />

[5] k = 3.56*(10^4)*exp(-5027/T)<br />

[6] K = 25000*exp(-10.69*(T-To)/T)<br />

[7] CB = CAo*X*To/T/2/(1-0.5*X)<br />

[8] CA = CAo*(1-X)*To/T/(1-0.5*X)<br />

[9] R = k*((CA^2)-(CB/K))<br />

[10] A = 2*K*To*CAo<br />

[11] B = (2*A)+T<br />

[12] Xe = (B-(((B^2)-(2*A*B))^0.5))/B<br />

351

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!