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
33
34
35
|
//Hougen O.A., Watson K.M., Ragatz R.A., 2004. Chemical process principles Part-1: Material and Energy Balances(II Edition). CBS Publishers & Distributors, New Delhi, pp 504
//Chapter-8, Illustration 6, Page 282
//Title: Calculation of total enthalpy
//=============================================================================
clear
clc
//INPUT
AW = 65.4; //Atomic weight of zinc
T = [1000 0 419 907]; //Given temperature, solid state temperature, melting point and boiling point of zinc in degree C
CP = [0.105 0.109]; //Mean specific heat of solid from 0-419 degree C and liquid from 419-907 degree C in cal per gram degree C obtained from Fig 63, Page 260
lamda1 = 1660; //Heat of fusion in cal per g-atom obtained from Table 24, Page 272
CP1 = 4.97; //Molal heat capacity of zinc vapor at constant preesure in cal per g-mole
//CALCULATION
T1 = T+273; //Given temperature, solid state temperature, melting point and boiling point of zinc in K
lamda2 = T1(4)*(8.75+4.571*log10(T1(4))); //Heat of vaporization at normal boiling point in cal per g-mole
Lamda1 = CP(1)*(T(3)-T(2)); //Heat absorbed by solid in cal per gram
Lamda2 = lamda1/AW; //Heat of fusion in cal per gram
Lamda3 = CP(2)*(T(4)-T(3)); //Heat absorbed by liquid in cal per gram
Lamda4 = lamda2/AW; //Heat of vaporization in cal per gram
Lamda5 = CP1*(T(1)-T(4))/AW; //Heat absorbed by vapor in cal per gram
Lamda = Lamda1+Lamda2+Lamda3+Lamda4+Lamda5; //Total enthalpy in cal per gram
//OUTPUT
// Console Output
mprintf('\n Total enthalpy of zinc vapor at 1000 degree C = %3.0f cal per gram',Lamda);
// File Output
fd= mopen('.\Chapter8_Example6_Output.txt','w');
mfprintf(fd,'\n Total enthalpy of zinc vapor at 1000 degree C = %3.0f cal per gram',Lamda);
mclose(fd);
//=============================END OF PROGRMAM=================================
|
