1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
# Numerical parameter file (SimpleGasTurbine_numpar.txt)
# Generated by MTT at Tue Mar 31 12:15:00 BST 1998
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% Version control history
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% $Id$
# %% $Log$
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#Dummies
alpha = 1;
c_v = 1;
density = 1;
ideal_gas = 1;
|
>
>
>
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
# Numerical parameter file (SimpleGasTurbine_numpar.txt)
# Generated by MTT at Tue Mar 31 12:15:00 BST 1998
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% Version control history
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% $Id$
# %% $Log$
# %% Revision 1.1 2000/12/28 16:55:29 peterg
# %% To RCS
# %%
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#Dummies
alpha = 1;
c_v = 1;
density = 1;
ideal_gas = 1;
|
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
m_c = (p_3*v_c)/(t_3*r);
%Equate pressures
p_4 = p_1;
p_2 = p_3;
%Compute ss temperatures (isentropic)
t_2 = t_1*(p_2/p_1)^alpha;
t_4 = t_3*(p_4/p_3)^alpha;
%Find the steady-state work output
w_0 = c_p*(t_3-t_4) - c_p*(t_2-t_1);
%Unit mass flow
mdot = 1;
%Corresponding shaft speed
omega_0 = mdot/k;
%Compute the corresponding load resistance (to absorb that work)
r_l = w_0/(omega_0)^2;
%Compute shaft inertia to give unit time constant (j_s*r_l)
j_s = r_l;
%Find angular momentum to give shaft speed omega_0
mom_0 = omega_0*j_s;
|
|
|
|
|
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
|
m_c = (p_3*v_c)/(t_3*r);
%Equate pressures
p_4 = p_1;
p_2 = p_3;
%Compute ss temperatures (isentropic)
t_2 = t_1*pow((p_2/p_1),alpha);
t_4 = t_3*pow((p_4/p_3),alpha);
%Find the steady-state work output
w_0 = c_p*(t_3-t_4) - c_p*(t_2-t_1);
%Unit mass flow
mdot = 1;
%Corresponding shaft speed
omega_0 = mdot/k;
%Compute the corresponding load resistance (to absorb that work)
r_l = w_0/pow((omega_0),2);
%Compute shaft inertia to give unit time constant (j_s*r_l)
j_s = r_l;
%Find angular momentum to give shaft speed omega_0
mom_0 = omega_0*j_s;
|