summaryrefslogtreecommitdiff
path: root/23/CH9/EX9.3/Example_9_3.sce
blob: fdb1682a7e85b36f214b82fb8104e3c668e3d767 (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
38
39
40
41
42
43
44
45
46
47
48
clear;
clc;

//To find Approx Value
function[A]=approx(V,n)
  A=round(V*10^n)/10^n;//V-Value  n-To what place
  funcprot(0)
endfunction  

//Example 9.3
//Caption : Program to Find the Temperature of the High Pressure steam 

//Given Values

x=0.25;
//For Superheated Methane
//By Perry and Green  by linear interpolation
H4=1140;//[KJ/Kg]  @ 300K and 60 Bar
H15=1188.9;//[KJ/Kg]  @ 295K and 1 Bar
//By interpolation based on lnP
T_sat=111.5;//[K]
H9=285.4;//[KJ/Kg]  Saturated Liquid
H12=769.9;//[KJ/Kg] Saturated Vapor
S12=9.521;//[KJ/Kg/K] Saturated vapor

T5=253.6;//[K]
H5=1009.8;//[KJ/Kg]  @ 60 Bar

//From Eqn(9.7)
z=((x*(H12-H5))+H4-H15)/(H9-H15);

H14=((H5-H4)/(1-z))+H15;//[KJ/Kg]
//Whence
T14=227.2;//[K]  @ 1Bar

H7=H5-((1-z)/(1-x)*(H14-H12));//[KJ/Kg]
T7=197.6;//[K]  @ 60Bar

//From Eqn(9.8)
z=approx((H4-H15)/(H9-H15),4);

H7=H4-((1-z)*(H15-H12));
T7=206.6;//[K]

disp('%',z*100,'Fraction of methane liquefied')
disp('K',T7,'Temperature of High Pressure steam entering the throttle valve')

//End