/[ascend]/trunk/models/johnpye/thermo_types.a4c
ViewVC logotype

Contents of /trunk/models/johnpye/thermo_types.a4c

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2067 - (show annotations) (download) (as text)
Thu Aug 13 10:39:43 2009 UTC (11 years ago) by jpye
File MIME type: text/x-ascend
File size: 4144 byte(s)
Add thermal_resistance type.
Formatting of thermodynamics.a4l.
Changing to 'molar_gas_constant' in eos-vanderwaals.
1 REQUIRE "atoms.a4l";
2 REQUIRE "stream_holdup.a4l";
3
4 CONSTANT specific_gas_constant
5 REFINES real_constant DIMENSION L^2/T^2/TMP;
6
7 CONSTANT mass_density_constant
8 REFINES real_constant DIMENSION M/L^3;
9
10 ATOM heat_transfer_coefficient REFINES solver_var
11 DIMENSION M/T^3/TMP
12 DEFAULT 5{W/m^2/K};
13
14 lower_bound := 0{W/m^2/K};
15 upper_bound := 1e50{W/m^2/K};
16 nominal := 5{W/m^2/K};
17 END heat_transfer_coefficient;
18
19 ATOM specific_enthalpy REFINES solver_var
20 DIMENSION L^2/T^2
21 DEFAULT 1000{kJ/kg};
22
23 lower_bound := 0{kJ/kg};
24 upper_bound := 1e50{kJ/kg};
25 nominal := 1000{kJ/kg};
26
27 END specific_enthalpy;
28
29 ATOM specific_entropy REFINES solver_var
30 DIMENSION L^2/T^2/TMP
31 DEFAULT 6{kJ/kg/K};
32
33 lower_bound := 0{kJ/kg/K};
34 upper_bound := 1e50{kJ/kg/K};
35 nominal := 6{kJ/kg/K};
36
37 END specific_entropy;
38
39 ATOM specific_heat_capacity REFINES solver_var
40 DIMENSION L^2/T^2/TMP
41 DEFAULT 4.0{kJ/kg/K};
42
43 lower_bound := 0{kJ/kg/K};
44 upper_bound := 1e50{kJ/kg/K};
45 nominal := 4.2{kJ/kg/K};
46 END specific_heat_capacity;
47
48 ATOM heat_capacity REFINES solver_var
49 DIMENSION M*L^2/T^2/TMP
50 DEFAULT 100{kJ/K};
51 lower_bound := 0.0{kJ/K};
52 upper_bound := 1e50{kJ/K};
53 nominal := 100{kJ/K};
54 END heat_capacity;
55
56 ATOM specific_volume REFINES solver_var
57 DIMENSION L^3/M
58 DEFAULT 0.001{m^3/kg};
59
60 lower_bound := 0{m^3/kg};
61 upper_bound := 1e50{m^3/kg};
62 nominal := 0.001{m^3/kg};
63
64 END specific_volume;
65
66 ATOM specific_energy REFINES solver_var
67 DIMENSION L^2/T^2
68 DEFAULT 1000{kJ/kg};
69
70 lower_bound := 0{kJ/kg};
71 upper_bound := 1e50{kJ/kg};
72 nominal := 1000{kJ/kg};
73 END specific_energy;
74
75 ATOM delta_specific_enthalpy REFINES solver_var
76 DIMENSION L^2/T^2
77 DEFAULT 1000{kJ/kg};
78
79 lower_bound := -10000{kJ/kg};
80 upper_bound := 10000{kJ/kg};
81 nominal := 1000{kJ/kg};
82
83 END delta_specific_enthalpy;
84
85 ATOM specific_power REFINES solver_var
86 DIMENSION L^2/T^3
87 DEFAULT 100{W/kg};
88
89 lower_bound := 0{W/kg};
90 upper_bound := 1e50{W/kg};
91 nominal := 100{W/kg};
92 END specific_power;
93
94 ATOM delta_specific_power REFINES solver_var
95 DIMENSION L^2/T^3
96 DEFAULT 100{W/kg};
97
98 lower_bound := -1e50{W/kg};
99 upper_bound := 1e50{W/kg};
100 nominal := 100{W/kg};
101 END delta_specific_power;
102
103 ATOM specific_energy_rate REFINES solver_var
104 DIMENSION L^2/T^3
105 DEFAULT 100 {J/kg/s};
106
107 lower_bound := -1e50 {kJ/kg/s};
108 upper_bound := 1e50 {kJ/kg/s};
109 nominal := 100{ J/kg/s};
110 END specific_energy_rate;
111
112 ATOM specific_enthalpy_rate REFINES solver_var
113 DIMENSION L^2/T^3
114 DEFAULT 100 {J/kg/s};
115 lower_bound := -1e8 {kJ/kg/s};
116 upper_bound := 1e8 {kJ/kg/s};
117 nominal := 100{ J/kg/s};
118 END specific_enthalpy_rate;
119
120 ATOM ua_value REFINES solver_var
121 DIMENSION M/T^3/TMP*L^2
122 DEFAULT 1{kW/K};
123
124 lower_bound := 0{W/K};
125 upper_bound := 1e50{W/K};
126 nominal := 1{kW/K}; (* heat up 1kg water by 1deg in 5 s *)
127 END ua_value;
128
129 ATOM thermal_resistance REFINES solver_var
130 DIMENSION TMP*T^3/M/L^2
131 DEFAULT 0.5 {K/W};
132 lower_bound := 0 {K/W};
133 upper_bound := 1e12 {K/W};
134 nominal := 0.5 {K/W};
135 END thermal_resistance;
136
137 ATOM R_value REFINES solver_var
138 DIMENSION TMP*T^3/M
139 DEFAULT 1{K*m^2/W};
140 lower_bound := 0{K*m^2/W};
141 upper_bound := 1e6 {K*m^2/W};
142 nominal := 1 {K*m^2/W};
143 END R_value;
144
145 ATOM pressure_per_temperature REFINES solver_var
146 DIMENSION M/L/T^2/TMP
147 DEFAULT 1.0{Pa/K};
148 lower_bound := -1e50{Pa/K};
149 upper_bound := 1e50{Pa/K};
150 nominal := 1.0{Pa/K};
151 END pressure_per_temperature;
152
153 ATOM energy_rate_per_length REFINES solver_var
154 DIMENSION M*L/T^3
155 DEFAULT 1000{W/m};
156 lower_bound := -1e50{W/m};
157 upper_bound := 1e50{W/m};
158 nominal := 1000{W/m};
159 END energy_rate_per_length;
160
161 ATOM energy_flux REFINES solver_var
162 DIMENSION M/T^3
163 DEFAULT 1000{W/m^2};
164 lower_bound := -1e50{W/m^2};
165 upper_bound := 1e50{W/m^2};
166 nominal := 1000{W/m^2};
167 END energy_flux;
168
169 (* for use in heat exchangers: C_c, C_h, C_min, etc. *)
170 ATOM capacity_rate REFINES power_per_temperature;
171 END capacity_rate;
172
173 MODEL thermo_state;
174 T IS_A temperature;
175 rho IS_A mass_density;
176 p IS_A pressure;
177 u IS_A specific_energy;
178 h IS_A specific_enthalpy;
179 s IS_A specific_entropy;
180 cp IS_A specific_heat_capacity;
181 cv IS_A specific_heat_capacity;
182 (* w IS_A speed; *)
183 END thermo_state;

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22