Расчет значений k0 и e в уравнении аррениуса с использованием метода наименьших квадратов и поиска решения excel. Расчет значений константы равновесия с использованием метода наименьших квадратов, страница 7

Твх,°С	Плот ор,м3/мг	Объем, м3	Твых,°С	Степ абс,%
170	13	22	58,9	69,06
180	14	25	52,4	77,10
170	13	30	44,0	87,15
160	18	21	46,1	80,4
188	17	27	43,8	88,45
200	16	24	54,1	76,63
210	19	22	53,3	76,72
150	20	25	32,1	99,12
174	21	26	33,6	98,57
182	21	26	35,3	96,81
190	21	26	37,0	95,16
170	18	26	39,5	92,11
160	17	29	35,4	99,54
170	15	24	49,8	79,19
180	15	24	52,1	77,4
190	15	24	54,3	75,74
210	15	24	58,7	72,77
225	16	22	64,3	68,49
210	18	29	43,0	91,87
150	18	19	47,8	76,54
186	14	25	53,7	76,09
190	14	25	54,6	75,45

приложение 2

Зависимость константы равновесия от температуры

Т,°С	Кр эксп
400	443
420	265
440	112
460	108
480	72
520	35
540	25
560	16
580	12
600	9,5
620	7

приложение 3

Зависимость константы скорости реакции от температуры

Т, °С	К, 1/с
400	0,4
405	0,4
415	0,5
435	0,6
455	0,8
485	1,3
505	1,6
515	1,8
530	2,2
575	3,5
595	4,3
605	4,7
615	5,2

приложение 4

n=11: r=8.315

dim t(n), ke(n), x(n), y(n), kr(n)

s1=0: s2=0: s3=0: s4=0

data 400,443,…

for i = 1 to n

read t(i), ke(i)

x(i)=1/(t(i)+273): y(i)=log(ke(i))

s1 = s1+x(i): s2=s2+y(i): s3=s3+x(i)^2: s4=s4+x(i)*y(i)

next i

z=s3*n-s1^2

a=(s4*n-s2*s1)/z: ea=-a*r

b=(s3*s2-s1*s4)/z: k0=exp(b)

? “k0=”; k0

? “ea=”; ea

for i=1 to n

kr(i)=k0*exp(-ea/(r*(t(i)+273)))

? using “######.###”; t(i), ke(i), kr(i)

next i

приложение 5

n=11: r=8,315

dim kr(n), k(n), t(n),x(n),y(n)

data 400, 0.4, 405, 0.4, …

for i=1 to n

read t(i), k(i)

y(i)=log(k(i))/2.3

x(i)=1/(t(i)+273)

s1=s1+x(i)

s2=s2+y(i)

s3=s3+x(i)*x(i)

s4=s4+x(i)*Y(i)

next i

z=s3*n-s1^2

a=(s4*n-s1*s2)/z

b=(s3*s2-s1*s4)/z

? “a=”; a

? “b=”; b

for i=1 to n

kr(i)=10^(a/(t(i)+273)+b)

? using “####.## “; k(i), t(i), kr(i)

next i

приложение 6

n=10

input t0

input V

input g

data 21200, 0.33, 8.31441

read q, cp, r

input a0

input b0

input c0

t0=t0+273: g=g/3600: q=q/22.4: tau=V/g: d=tau/n: a=a0: b=b0: c=c0:t=t0

for i=1 to n

k=19949.17*exp(-60914.62/(r*t0))

kp=10^(4794.21/t0-4.5248)

w=k*a*b/(a+0.8*c)*(1-(c/(kp*a*b^0.5))^2)

p=q/cp*w

a=a-w*d: t=t+p*d

b=b0-0.5*(a0-a): c=c0-(a-a0)

next I

st=(a0-a)/a0: g=g*(1-a0*st/2)*3600

? “t=”; t-273

? “g=”; g

? “a=”; a

? “b=”; b

? “c=”; c

приложение 7

dim txh(10), txk(10), tgh(10), tgk(10), wx(10), wg(10)

k=0.01745: cena=483: td=20

input nx, ng

for i=1 to nx

input txh(i), txk(i), wx(i)

next i

? “x  txh  txk  g  tgh  tgk  stoim”

for i=1 to nx

for g=1 to ng

tx0=txh(i): tx1=txk(i)

tg0=tgh(g): tg1=tgk(g)

y=(tg0-tx0-td)

if y<0 then 5 else goto 6

5 print using “## “; i, g

print “tepla net”

goto 4

6 if (tg1-tx0)<td then tg1=tx0+td

if (tg0-tx1)<tg then tx1=tg0-td

qx=(tx1-tx0)*wx(i)

qg=(tg0-tg1)*wg(g)

q=qx

if qx>qg then q=qg

tx1=tx0+q/wx(i)

tg1=tg0-q/wg(g)

dt1=tg0-tx1

dt2=tg1-tx0

r=log(dt1/dt2)

dt=(dt1+dt2)/2

if abs(r)>0.693 then dt=(dt1-dt2)/r

f=q/k/dt

stoim=cena*f^0.6

print using “#  ###.#  ###.#  #  ###.#  ###.#  ######## “;i, tx0, tx1, g

4 next: next

список литературы