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 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 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;

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 specific_power REFINES solver_var
77	DIMENSION L^2/T^3
78	DEFAULT 100{W/kg};
79
80	lower_bound := 0{W/kg};
81	upper_bound := 1e50{W/kg};
82	nominal := 100{W/kg};
83	END specific_power;
84
85	ATOM delta_specific_power REFINES solver_var
86	DIMENSION L^2/T^3
87	DEFAULT 100{W/kg};
88
89	lower_bound := -1e50{W/kg};
90	upper_bound := 1e50{W/kg};
91	nominal := 100{W/kg};
92	END delta_specific_power;
93
94	ATOM ua_value REFINES solver_var
95	DIMENSION M/T^3/TMP*L^2
96	DEFAULT 1{kW/K};
97
98	lower_bound := 0{W/K};
99	upper_bound := 1e50{W/K};
100	nominal := 1{kW/K}; (* heat up 1kg water by 1deg in 5 s *)
101	END ua_value;
102
103	MODEL thermo_state;
104	T IS_A temperature;
105	rho IS_A mass_density;
106	p IS_A pressure;
107	u IS_A specific_energy;
108	h IS_A specific_enthalpy;
109	s IS_A specific_entropy;
110	cp IS_A specific_heat_capacity;
111	cv IS_A specific_heat_capacity;
112	(* w IS_A speed; *)
113	END thermo_state;