Problem 9.35 in FelderStream Compounds Reaction1 Feed 1 C2H6 1 C2H6 -> C2H4 + H22 Product 2 C2H4 3 H2 n1 _____________ n2-------->| | ------------> n11 | reactor | n21, n22 |_____________| n23Using the Ceng 301 modules and Data BankHere are your compounds' formulae and names:
No. Formula Name
----------------------------------------
1 C2H6 ethane
2 C2H4 ethylene
3 H2 hydrogen
Here are your reactions:
----------------------------------------
1) C2H6 --> C2H4 + H2
Enter the number of streams: 2
The variables for your compounds have now been created,
>>setne('v',1,1273,1,1)
>>global P f
>>[fs,ts]=p935([.01 .05 .1 .5 1 5 10])
fs =
0.4258 0.3867 0.3683 0.3220 0.3005 0.2474 0.2236
ts =
1.0e+03 *
0.7735 0.8244 0.8478 0.9055 0.9318 0.9953 1.0233
P in atm T in K fraction conv Maple Matlab Maple Matlab 0.01 785.5 773.5 0.4746 0.4258 0.1 864.5 847.8 0.4164 0.3683 1.0 954.1 931.8 0.3425 0.3005 10.0 1050.5 1023.3 0.2527 0.2236restart;function [fs,ts]=p935(Ps)
global P f
fs=[];
ts=[];
kmax=length(Ps);
for k=1:kmax
P=Ps(k);
t=ssec2([300 1200],'er935(x)');
fs=[fs,f];
ts=[ts,t];
endfunction e=er935(t)
global P f
Kp=kequil(t);
f=sqrt(Kp/(P+Kp));
e=reacte('v',t,1,2,f);