models/steam/iapwssat.a4c

REQUIRE "steam/iapwssatprops.a4c";

MODEL iapwssat REFINES iapwssatprops;
	x IS_A fraction; (* ensures we stay in the required region *)
	h IS_A specific_enthalpy;
	u IS_A specific_energy;
	v IS_A specific_volume;

	h_expr: h = (h_g - h_f)*x + h_f;
	u_expr: u = (u_g - u_f)*x + u_f;
	v_expr: v = (rho_g + (rho_f - rho_g)*x ) / ( rho_f * rho_g );
METHODS
METHOD specify;
	RUN iapwssatprops::specify;
	FIX x;
END specify;
METHOD values;
	RUN iapwssatprops::values;
	x := 0.5;
END values;
METHOD on_load;
	RUN reset; RUN values;
	RUN bound_self;
END on_load;
END iapwssat;

MODEL testiapwssat1 REFINES iapwssat;
METHODS
METHOD values;
	RUN iapwssatprops::values;
	x := 1.01;
	x.upper_bound := 1.0;
	x.lower_bound := 0.0;
END values;
METHOD self_test;
	ASSERT x <= x.upper_bound;
	ASSERT x >= x.lower_bound;
END self_test;
END testiapwssat1;
1	REQUIRE "steam/iapwssatprops.a4c";
2
3	MODEL iapwssat REFINES iapwssatprops;
4	x IS_A fraction; (* ensures we stay in the required region *)
5	h IS_A specific_enthalpy;
6	u IS_A specific_energy;
7	v IS_A specific_volume;
8
9	h_expr: h = (h_g - h_f)*x + h_f;
10	u_expr: u = (u_g - u_f)*x + u_f;
11	v_expr: v = (rho_g + (rho_f - rho_g)x ) / ( rho_f rho_g );
12	METHODS
13	METHOD specify;
14	RUN iapwssatprops::specify;
15	FIX x;
16	END specify;
17	METHOD values;
18	RUN iapwssatprops::values;
19	x := 0.5;
20	END values;
21	METHOD on_load;
22	RUN reset; RUN values;
23	RUN bound_self;
24	END on_load;
25	END iapwssat;
26
27	MODEL testiapwssat1 REFINES iapwssat;
28	METHODS
29	METHOD values;
30	RUN iapwssatprops::values;
31	x := 1.01;
32	x.upper_bound := 1.0;
33	x.lower_bound := 0.0;
34	END values;
35	METHOD self_test;
36	ASSERT x <= x.upper_bound;
37	ASSERT x >= x.lower_bound;
38	END self_test;
39	END testiapwssat1;