blob: 465c5ce82b4bebce3b1c7f6ea7e311da1921230a (
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
33
34
35
36
37
|
// Y.V.C.Rao ,1997.Chemical Engineering Thermodynamics.Universities Press,Hyderabad,India.
//Chapter-4,Example 10,Page 105
//Title:Final temperature
//================================================================================================================
clear
clc
//INPUT
N=100;//number of moles of carbon dioxide in mol
T1=298;//initial temperature in K
Q=1;//energy added as heat in MJ
a=45.369;//coefficient in the specific heat capacity expression where Cp=a+bT+eT^-2
b=8.688*10^-3;//coefficient in the specific heat capacity expression where Cp=a+bT+eT^-2
e=-9.619*10^5;//coefficient in the specific heat capacity expression where Cp=a+bT+eT^-2
//Where Cp is in J/molK
//CALCULATION
delh=Q*10^6/N;//calculation of enthalpy in J
Tguess=520;//The final temperature guess value in K used for solving the system of equations
//The system of equations are defined by :
//T2=T1+(delh/Cpm)--->A
//Cpm=a+(b*Tm)+(e/T1T2)--->B
Cpm_guess=a+(b*((T1+Tguess)/2))+(e/(T1*Tguess));//calculation of Cpm guess (in J/molK) to be used to determine T2 from Equation A
T2_guess=T1+(delh/Cpm_guess);//calculation of T2 using Equation A (in K) based on the value of Cpm guess computed using Equation B
tolerance=1e-6;//defining the tolerance limit to obtain convergence
while abs(T2_guess-Tguess)>tolerance
Tguess=T2_guess;
Cpm_guess=a+(b*((T1+Tguess)/2))+(e/(T1*Tguess));
T2_guess=T1+(delh/Cpm_guess);//the iteration process to solve the system of equations
end
T2=T2_guess;//value of the final temperature of CO2 obtained upon solving the system of equations(A and B) in K
//OUTPUT
mprintf('\n The final temperature= %0.1f K\n',T2);
//===============================================END OF PROGRAM===================================================
|
