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
|
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
|
+
+
+
-
-
+
+
-
-
+
+
-
-
+
+
-
-
+
+
|
# Numerical parameter file (OttoCycle_input.txt)
# Generated by MTT at Thu Dec 4 11:17:09 GMT 1997
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% Version control history
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# %% $Id$
# %% $Log$
# %% Revision 1.2 2000/12/28 18:15:52 peterg
# %% To RCS
# %%
# %% Revision 1.1 1998/08/10 14:42:13 peterg
# %% Initial revision
# %%
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Set the inputs
if ((t>=0.0)&&(t<1.0)) #Compression
u(1) = 0.0; # Entropy flow
u(2) = -0.8; # Volume rate-of-change
OttoCycle_Heat = 0.0; # Entropy flow
OttoCycle_Work = -0.8; # Volume rate-of-change
endif;
if ((t>=1.0)&&(t<2.0)) #Heating
u(1) = 1000; # Entropy flow
u(2) = 0.0; # Volume rate-of-change
OttoCycle_Heat = 1000; # Entropy flow
OttoCycle_Work = 0.0; # Volume rate-of-change
endif;
if ((t>=2.0)&&(t<3.0)) #Expansion
u(1) = 0.0; # Entropy flow
u(2) = 0.8; # Volume rate-of-change
OttoCycle_Heat = 0.0; # Entropy flow
OttoCycle_Work = 0.8; # Volume rate-of-change
endif;
if (t>=3.0) #Cooling
u(1) = -1000; # Entropy flow
u(2) = 0.0; # Volume rate-of-change
OttoCycle_Heat = -1000; # Entropy flow
OttoCycle_Work = 0.0; # Volume rate-of-change
endif;
|