Bulletin Board

# 26 Non-isothermal CSTR part (c)

Calculated values of NLE variables

	Variable	Value	f(x)	Initial Guess
1	FA	0.2487268	-5.193E-10	0.3
2	FB	0.4598077	5.456E-10	0.3
3	FC	0.0731988	-1.83E-11	0.3
4	FD	0.0625495	-1.562E-11	0.3

	Variable	Value
1	CA	0.035927
2	CA0	0.1626016
3	CB	0.0664163
4	CC	0.0105731
5	CD	0.0090349
6	CPA	16.
7	CPB	16.
8	CPC	16.
9	CPD	10.
10	dCP	4.
11	E1	5600.
12	E2	6000.
13	E3	7000.
14	FA0	0.8130081
15	FT	0.8442829
16	HA	1.46E+04
17	HB	5500.
18	HC	4000.
19	HD	3000.
20	k1	0.1570634
21	k2	0.0110212
22	k3	0.0047089
23	PA0	4.
24	Q	-3915.548
25	R	0.082
26	R2	1.987
27	rA	-0.0056428
28	rB	0.0045981
29	rC	0.000732
30	rD	0.0006255
31	SCD	1.170254
32	sumrdH	52.16361
33	T	400.
34	T0	300.
35	v	6.92312
36	V	100.
37	v0	5.
38	X	0.694066

Nonlinear equations

1	f(FA) = FA0+rA*V-FA = 0
2	f(FB) = rB*V-FB = 0
3	f(FC) = rC*V-FC = 0
4	f(FD) = rD*V-FD = 0

Explicit equations

1	PA0 = 4
2	v0 = 5
3	T0 = 300
4	R = 0.082
5	CA0 = PA0/(R*T0)
6	FA0 = v0*CA0
7	T = 400
8	V = 100
9	FT = FA+FB+FC+FD
10	v = v0*FT/FA0*T/T0
11	CA = FA/v
12	CB = FB/v
13	CC = FC/v
14	CD = FD/v
15	R2 = 1.987
16	E1 = 5600
17	E2 = 6000
18	E3 = 7000
19	k1 = 1.5E-2*exp(-E1/R2*(1/T-1/T0))
20	k2 = 8.9E-4*exp(-E2/R2*(1/T-1/T0))
21	k3 = 2.5E-4*exp(-E3/R2*(1/T-1/T0))
22	rA = -k1*CA
23	rB = k1*CA-k2*CB-k3*CB
24	rC = k2*CB
25	rD = 2*k3*CB
26	SCD = FC/FD
27	X = (FA0-FA)/FA0
28	HC = 4000
29	HD = 3000
30	CPD = 10
31	CPA = 16
32	CPB = 16
33	CPC = 16
34	dCP = 2*CPD-CPB
35	HA = 14600
36	HB = 5500
37	sumrdH = -k1*CA*(HB-HA)-k2*CB*(HC-HB)-k3*CB*((2*HD-HB)+dCP*(T-298))
38	Q = FA0*CPA*(T-T0)-V*sumrdH