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
|
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
|
-
-
-
-
+
+
+
+
-
+
-
-
-
-
+
+
+
+
-
-
+
+
-
-
-
-
+
+
+
+
-
+
-
-
-
-
+
+
+
+
-
+
-
|
%SUMMARY polytrop CR for gas turbine compressor
OPERATOR polytrop;
% Port 1 generates zero flow
FOR ALL deltaP,temperature,pressure,k,deltaT
LET polytrop(k, flow, 1,
deltaP,effort,1,
deltaT,effort,2,
FOR ALL Ipressure,temperature,Fpressure,gamma,enthflow
LET polytrop(gamma, flow, 1,
Fpressure,effort,1,
enthflow,flow,2,
temperature,effort,3,
pressure,effort,4)
Ipressure,effort,4)
= 0;
% Port 2 generates deltaT
FOR ALL deltaP,temperature,pressure,k,deltaT
LET polytrop(k, effort, 2,
deltaP,effort,1,
deltaT,effort,2,
FOR ALL Ipressure,temperature,Fpressure,gamma,enthflow
LET polytrop(gamma, effort, 2,
Fpressure,effort,1,
enthflow,flow,2,
temperature,effort,3,
pressure,effort,4)
= temperature*((1-(deltaP/pressure)^((k-1)/k)-1);
Ipressure,effort,4)
= temperature*((Ipressure/Fpressure)^(gamma)-1);
% Port 3 generates zero flow
FOR ALL deltaP,temperature,pressure,k,deltaT
LET polytrop(k, flow, 3,
deltaP,effort,1,
deltaT,effort,2,
FOR ALL Ipressure,temperature,Fpressure,gamma,enthflow
LET polytrop(gamma, flow, 3,
Fpressure,effort,1,
enthflow,flow,2,
temperature,effort,3,
pressure,effort,4)
Ipressure,effort,4)
= 0;
% Port 4 generates zero flow
FOR ALL deltaP,temperature,pressure,k,deltaT
LET polytrop(k, flow, 4,
deltaP,effort,1,
deltaT,effort,2,
FOR ALL Ipressure,temperature,Fpressure,gamma,enthflow
LET polytrop(gamma, flow, 4,
Fpressure,effort,1,
enthflow,flow,2,
temperature,effort,3,
pressure,effort,4)
Ipressure,effort,4)
= 0;
|