blob: 48486f8b99b74ff56dc05011252e7366c2604d34 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
//Harriot P.,2003,Chemical Reactor Design (I-Edition) Marcel Dekker,Inc.,USA,pp 436.
//Chapter-5 Ex5.2 Pg No. 194
//Title: Maximum internal temperature difference
//=============================================================================================================
clear
format(16)
clc
//INPUT
T_C=200;//Temperature(°C)
P=1.2;//Pressure (atm)
f_ethylene=0.05;//fraction of ethylene
k_s=8*10^(-4);//Solid conductivity (cal/sec cm°C)
D_e=0.02;//Diffusivity for ethylene (cm2/s)
del_H= -32.7*10^(3);//Heat of reaction (cal)
V_ref=22400;// reference volume(cm3)
T_ref=273;//Reference Temperature (K)
P_ref=1;//Reference Pressure (atm)
T_K=T_C+273;//Reaction Temperature (K)
//CALCULATION
C_s=f_ethylene*P*T_ref/(V_ref*T_K*P_ref);
Tc_minus_Ts=D_e*C_s*(-del_H)/k_s;//Refer equation 5.51 Pg No. 194
//OUTPUT
mprintf('\n\tThe maximum internal temperature difference %0.3f °C',Tc_minus_Ts);
//FILE OUTPUT
fid= mopen('.\Chapter5-Ex2-Output.txt','w');
mfprintf(fid,'\n\tThe maximum internal temperature difference %0.3f °C',Tc_minus_Ts);
mclose(fid);
//=====================================================END OF PROGRAM=================================================
|
