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$
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Parameters
c_p = 1005.0;
c_v = 718.0;
gamma_0 = c_p/c_v;
alpha = (gamma_0-1)/gamma_0;
|
>
>
>
|
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 18:08:28 peterg
# %% To RCS
# %%
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Parameters
c_p = 1005.0;
c_v = 718.0;
gamma_0 = c_p/c_v;
alpha = (gamma_0-1)/gamma_0;
|
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
|
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;
|
|
|
|
|
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
|
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;
|