models/johnpye/thermo_types.a4c

REQUIRE "atoms.a4l";
REQUIRE "stream_holdup.a4l";

CONSTANT specific_gas_constant
REFINES real_constant DIMENSION L^2/T^2/TMP;

CONSTANT mass_density_constant
REFINES real_constant DIMENSION M/L^3;

ATOM heat_transfer_coefficient REFINES solver_var
		DIMENSION M/T^3/TMP
		DEFAULT 5{W/m^2/K};

	lower_bound := 0{W/m^2/K};
	upper_bound := 1e50{W/m^2/K};
	nominal := 5{W/m^2/K};

END heat_transfer_coefficient;

ATOM specific_enthalpy REFINES solver_var
		DIMENSION L^2/T^2
		DEFAULT 1000{kJ/kg};

	lower_bound := 0{kJ/kg};
	upper_bound := 1e50{kJ/kg};
	nominal := 1000{kJ/kg};

END specific_enthalpy;

ATOM specific_entropy REFINES solver_var
		DIMENSION L^2/T^2/TMP
		DEFAULT 6{kJ/kg/K};
	
	lower_bound := 0{kJ/kg/K};
	upper_bound := 1e50{kJ/kg/K};
	nominal := 6{kJ/kg/K};

END specific_entropy;

ATOM specific_heat_capacity REFINES solver_var
		DIMENSION L^2/T^2/TMP
		DEFAULT 4.0{kJ/kg/K};
	
	lower_bound := 0{kJ/kg/K};
	upper_bound := 1e50{kJ/kg/K};
	nominal := 4.2{kJ/kg/K};
END specific_heat_capacity;

ATOM heat_capacity REFINES solver_var
    DIMENSION M*L^2/T^2/TMP
    DEFAULT 100{kJ/K};
    lower_bound := 0.0{kJ/K};
    upper_bound := 1e50{kJ/K};
    nominal := 100{kJ/K};
END heat_capacity;

ATOM specific_volume REFINES solver_var
		DIMENSION L^3/M
		DEFAULT 0.001{m^3/kg};
	
	lower_bound := 0{m^3/kg};
	upper_bound := 1e50{m^3/kg};
	nominal := 0.001{m^3/kg};

END specific_volume;

ATOM specific_energy REFINES solver_var
		DIMENSION L^2/T^2
		DEFAULT 1000{kJ/kg};

	lower_bound := 0{kJ/kg};
	upper_bound := 1e50{kJ/kg};
	nominal := 1000{kJ/kg};
END specific_energy;

ATOM delta_specific_enthalpy REFINES solver_var
		DIMENSION L^2/T^2
		DEFAULT 1000{kJ/kg};

	lower_bound := -10000{kJ/kg};
	upper_bound := 10000{kJ/kg};
	nominal := 1000{kJ/kg};

END delta_specific_enthalpy;

ATOM specific_power REFINES solver_var
		DIMENSION L^2/T^3
		DEFAULT 100{W/kg};

	lower_bound := 0{W/kg};
	upper_bound := 1e50{W/kg};
	nominal := 100{W/kg};
END specific_power;

ATOM delta_specific_power REFINES solver_var
		DIMENSION L^2/T^3
		DEFAULT 100{W/kg};

	lower_bound := -1e50{W/kg};
	upper_bound := 1e50{W/kg};
	nominal := 100{W/kg};
END delta_specific_power;

ATOM specific_energy_rate REFINES solver_var
		DIMENSION L^2/T^3
		DEFAULT 100 {J/kg/s};

	lower_bound := -1e50 {kJ/kg/s};
	upper_bound := 1e50 {kJ/kg/s};
	nominal := 100{ J/kg/s};
END specific_energy_rate;

ATOM specific_enthalpy_rate REFINES solver_var
		DIMENSION L^2/T^3 
		DEFAULT 100 {J/kg/s};
	lower_bound := -1e8 {kJ/kg/s};
	upper_bound := 1e8 {kJ/kg/s};
	nominal := 100{ J/kg/s};
END specific_enthalpy_rate;

ATOM ua_value REFINES solver_var
		DIMENSION M/T^3/TMP*L^2
		DEFAULT 1{kW/K};

	lower_bound := 0{W/K};
	upper_bound := 1e50{W/K};
	nominal := 1{kW/K}; (* heat up 1kg water by 1deg in 5 s *)
END ua_value;

ATOM pressure_per_temperature REFINES solver_var
    DIMENSION M/L/T^2/TMP
    DEFAULT 1.0{Pa/K};
    lower_bound := -1e50{Pa/K};
    upper_bound := 1e50{Pa/K};
    nominal := 1.0{Pa/K};
END pressure_per_temperature;

MODEL thermo_state;
	T IS_A temperature;
	rho IS_A mass_density;
	p IS_A pressure;
	u IS_A specific_energy;
	h IS_A specific_enthalpy;
	s IS_A specific_entropy;
	cp IS_A specific_heat_capacity;
	cv IS_A specific_heat_capacity;
	(* w IS_A speed; *)
END thermo_state;
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
18	END heat_transfer_coefficient;
19
20	ATOM specific_enthalpy REFINES solver_var
21	DIMENSION L^2/T^2
22	DEFAULT 1000{kJ/kg};
23
24	lower_bound := 0{kJ/kg};
25	upper_bound := 1e50{kJ/kg};
26	nominal := 1000{kJ/kg};
27
28	END specific_enthalpy;
29
30	ATOM specific_entropy REFINES solver_var
31	DIMENSION L^2/T^2/TMP
32	DEFAULT 6{kJ/kg/K};
33
34	lower_bound := 0{kJ/kg/K};
35	upper_bound := 1e50{kJ/kg/K};
36	nominal := 6{kJ/kg/K};
37
38	END specific_entropy;
39
40	ATOM specific_heat_capacity REFINES solver_var
41	DIMENSION L^2/T^2/TMP
42	DEFAULT 4.0{kJ/kg/K};
43
44	lower_bound := 0{kJ/kg/K};
45	upper_bound := 1e50{kJ/kg/K};
46	nominal := 4.2{kJ/kg/K};
47	END specific_heat_capacity;
48
49	ATOM heat_capacity REFINES solver_var
50	DIMENSION M*L^2/T^2/TMP
51	DEFAULT 100{kJ/K};
52	lower_bound := 0.0{kJ/K};
53	upper_bound := 1e50{kJ/K};
54	nominal := 100{kJ/K};
55	END heat_capacity;
56
57	ATOM specific_volume REFINES solver_var
58	DIMENSION L^3/M
59	DEFAULT 0.001{m^3/kg};
60
61	lower_bound := 0{m^3/kg};
62	upper_bound := 1e50{m^3/kg};
63	nominal := 0.001{m^3/kg};
64
65	END specific_volume;
66
67	ATOM specific_energy REFINES solver_var
68	DIMENSION L^2/T^2
69	DEFAULT 1000{kJ/kg};
70
71	lower_bound := 0{kJ/kg};
72	upper_bound := 1e50{kJ/kg};
73	nominal := 1000{kJ/kg};
74	END specific_energy;
75
76	ATOM delta_specific_enthalpy REFINES solver_var
77	DIMENSION L^2/T^2
78	DEFAULT 1000{kJ/kg};
79
80	lower_bound := -10000{kJ/kg};
81	upper_bound := 10000{kJ/kg};
82	nominal := 1000{kJ/kg};
83
84	END delta_specific_enthalpy;
85
86	ATOM specific_power REFINES solver_var
87	DIMENSION L^2/T^3
88	DEFAULT 100{W/kg};
89
90	lower_bound := 0{W/kg};
91	upper_bound := 1e50{W/kg};
92	nominal := 100{W/kg};
93	END specific_power;
94
95	ATOM delta_specific_power REFINES solver_var
96	DIMENSION L^2/T^3
97	DEFAULT 100{W/kg};
98
99	lower_bound := -1e50{W/kg};
100	upper_bound := 1e50{W/kg};
101	nominal := 100{W/kg};
102	END delta_specific_power;
103
104	ATOM specific_energy_rate REFINES solver_var
105	DIMENSION L^2/T^3
106	DEFAULT 100 {J/kg/s};
107
108	lower_bound := -1e50 {kJ/kg/s};
109	upper_bound := 1e50 {kJ/kg/s};
110	nominal := 100{ J/kg/s};
111	END specific_energy_rate;
112
113	ATOM specific_enthalpy_rate REFINES solver_var
114	DIMENSION L^2/T^3
115	DEFAULT 100 {J/kg/s};
116	lower_bound := -1e8 {kJ/kg/s};
117	upper_bound := 1e8 {kJ/kg/s};
118	nominal := 100{ J/kg/s};
119	END specific_enthalpy_rate;
120
121	ATOM ua_value REFINES solver_var
122	DIMENSION M/T^3/TMP*L^2
123	DEFAULT 1{kW/K};
124
125	lower_bound := 0{W/K};
126	upper_bound := 1e50{W/K};
127	nominal := 1{kW/K}; (* heat up 1kg water by 1deg in 5 s *)
128	END ua_value;
129
130	ATOM pressure_per_temperature REFINES solver_var
131	DIMENSION M/L/T^2/TMP
132	DEFAULT 1.0{Pa/K};
133	lower_bound := -1e50{Pa/K};
134	upper_bound := 1e50{Pa/K};
135	nominal := 1.0{Pa/K};
136	END pressure_per_temperature;
137
138	MODEL thermo_state;
139	T IS_A temperature;
140	rho IS_A mass_density;
141	p IS_A pressure;
142	u IS_A specific_energy;
143	h IS_A specific_enthalpy;
144	s IS_A specific_entropy;
145	cp IS_A specific_heat_capacity;
146	cv IS_A specific_heat_capacity;
147	(* w IS_A speed; *)
148	END thermo_state;