에탄을 열분해시켜 2.7 million Lb/stream day 의 에틸렌을 생산할 수 있는 반응기를 설계하시오.

반응물의 조성 (몰%): CH4: 3%; C2H6 94%; C3H8 3%

생성물의 조성: C2H4 99% 이상

반응기 압력: 10 psig 이상

반응기는 내부직경 6인치, 길이 32 ft 의 튜브를 활용

반응기온도: 튜브재질을 고려하여 1000oC를 초과하지 않도록 할 것.

반응속도 데이터:

reaction	preexponential factor	activation energy, kcal/mol
C2H6 --> 2CH3	1.0 x 1016, s-1	86
CH3+C2H6 --> CH4 +C2H5	3.16 x 108, L mol-1 s-1	10.8
C2H5 --> C2H4 + H	3.98 x 1013 s-1	38
H + C2H6 --> H2 + C2H5	1.25 x 1011, L mol-1 s-1	9.7
C2H5 + C2H5 --> C4H10	2.511 x 1010, L mol-1 s-1	0
C2H5 + C2H5 --> C2H6 + C2H4	2.511 x 1010, L mol-1 s-1	0
C2H5 + C2H4 --> C3H6 + CH3	3.16 x 109, L mol-1 s-1	19
H + C2H4 --> C2H5	5.0118 x 1010, L mol-1 s-1	-0.8

이라는 이번 설계과제인데요 다음 내용을 토대로 각 Species를 a,b,c,d,e 등으로 나누어서 폴리매스를 통해

값을 구하려고 하는데요 이상하게 에탄의 양이 줄지를 않는 형태로 되네요

아래의 글귀가 폴리매스에 적어놓은 내용입니다. 6.1기준이에요.

d(FA)/d(V) = rA #
d(FB)/d(V) = rB #
d(FC)/d(V) = rC #
d(FD)/d(V) = rD #
d(FE)/d(V) = rE #
d(FF)/d(V) = rF #
d(FG)/d(V) = rG #
d(FH)/d(V) = rH #
d(FI)/d(V) = rI #
r1A = -k1*0.581433338*FA/FT #
r2A = -k2*0.338064726*(FA/FT)*(FB/FT) #
r3D = -k3*0.581433338*FD/FT #
r4A = -k4*0.338064726*(FA/FT)*(FF/FT) #
r5D = -k5*0.338064726*(FD/FT)^2 #
r6D = -k6*0.338064726*(FD/FT)^2 #
r7E = -k7*0.338064726*(FD/FT)*(FE/FT) #
r8F = -k8*0.338064726*(FE/FT)*(FF/FT) #
k1 = 3.321227509 #
k2 = 3596252.41 #
k3 = 5760456.266 #
k4 = 2244173791 #
k5 = 25110000000 #
k6 = 25110000000 #
k7 = 1202192.429 #
k8 = 69820290516 #
FT =   FA+FB+FC+FD+FE+FF+FG+FH+FI   #
rA =   -r1A-r2A-r4A+2*r6D   #
rB =   r1A/2-r2A+r7E   #
rC =   r2A   #
rD =   r2A-r3D-r5D-r6D-r7E+r8F   #
rE =   r3D+r4A-r7E-r8F   #
rF =   r3D-r4A-r8F   #
rG =   r4A   #
rH =   2*r5D   #
rI =   r7E   #
V(0)=0
FA(0)=396000000
FB(0)=0
FC(0)=12638297.87

FD(0)=0
FE(0)=0
FF(0)=0
FG(0)=0
FH(0)=0
FI(0)=0
V(f)=11.67454

지금 어떤게 잘못되었는지를 잘모르겠구요... 도대체 왜이렇게되는지

혹은 이 설계문제와 관련되어서 다른 도움사항 주실 분이 있는지 모르겠습니다 ㅠ

도와주세요