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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
%SUMMARY RG CR for 2 port gas dymanics R: isentropic nozzle
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % $Id$
% % $Log$
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
OPERATOR RG;
% Port 1 : Enthalpy flow
FOR ALL gamma,R,A,T_u,P_u,T_d,P_d
LET RG(gamma,R,A, flow, 1,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) )
*(gamma/(gamma-1))*R*T_u;
% Port 2 : Mass flow
FOR ALL gamma,R,A,T_u,P_u,T_d,P_d
LET RG(gamma,R,A, flow, 2,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) );
% Port 3 : Enthalpy flow
FOR ALL gamma,R,A,T_u,P_u,T_d,P_d
LET RG(gamma,R,A, flow, 3,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) )
*(gamma/(gamma-1))*R*T_u;
% Port 4 : Mass flow
FOR ALL gamma,R,A,T_u,P_u,T_d,P_d
LET RG(gamma,R,A, flow, 4,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) );
END;;
|
>
>
>
|
|
|
|
|
|
|
|
|
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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
%SUMMARY RG CR for 2 port gas dymanics R: isentropic nozzle
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % $Id$
% % $Log$
% % Revision 1.1 1998/03/04 15:37:48 peterg
% % Initial revision
% %
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
OPERATOR RG;
% Port 1 : Enthalpy flow
FOR ALL COMP, gamma,R,A,T_u,P_u,T_d,P_d
LET RG(COMP, gamma,R,A, flow, 1,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) )
*(gamma/(gamma-1))*R*T_u;
% Port 2 : Mass flow
FOR ALL COMP, gamma,R,A,T_u,P_u,T_d,P_d
LET RG(COMP, gamma,R,A, flow, 2,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) );
% Port 3 : Enthalpy flow
FOR ALL COMP, gamma,R,A,T_u,P_u,T_d,P_d
LET RG(COMP, gamma,R,A, flow, 3,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) )
*(gamma/(gamma-1))*R*T_u;
% Port 4 : Mass flow
FOR ALL COMP, gamma,R,A,T_u,P_u,T_d,P_d
LET RG(COMP, gamma,R,A, flow, 4,
T_u,effort,1,
P_u,effort,2,
T_d,effort,3,
P_d,effort,4)
= A*( P_u/sqrt(T_u) )
*sqrt( (2*gamma)/(R*(gamma-1)) )
*sqrt( (P_d/P_u)^(2/gamma) - (P_d/P_u)^((1+gamma)/gamma) );
END;;
|