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;
END values;
METHOD self_test;
	ASSERT x == 1.0;
END self_test;
END testiapwssat1;
1	johnpye	1051	REQUIRE "steam/iapwssatprops.a4c";
2	johnpye	1042
3			MODEL iapwssat REFINES iapwssatprops;
4	johnpye	1043	x IS_A fraction; (* ensures we stay in the required region *)
5	johnpye	1042	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	johnpye	1043
27			MODEL testiapwssat1 REFINES iapwssat;
28			METHODS
29			METHOD values;
30			RUN iapwssatprops::values;
31			x := 1.01;
32			END values;
33	johnpye	1051	METHOD self_test;
34			ASSERT x == 1.0;
35			END self_test;
36	johnpye	1043	END testiapwssat1;