#*-------------------------------------------------------------------
* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
*
* This LIBRARY is free software; you can distribute it and/or modify
* it under the therms of the ALSOC FREE LICENSE as available at
* http://www.enq.ufrgs.br/alsoc.
*
* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
* from http://www.rps.eng.br Copyright (C) 2002-2004.
* All rights reserved.
*
* EMSO is distributed under the therms of the ALSOC LICENSE as
* available at http://www.enq.ufrgs.br/alsoc.
*
*--------------------------------------------------------------------
* Sample file for models valve_simplified and valve
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Paula B. Staudt, Estefane Horn, Núbia do Carmo Ferreira
* $Id: sample_valve.mso 587 2008-08-04 22:45:14Z bicca $
*--------------------------------------------------------------------*#

using "pressure_changers/valve";

FlowSheet valve_Test_V
# with vapour conditions
	PARAMETERS
	PP	as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = [ "isobutane", "benzene"],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp	as Integer;

	SET
	NComp = PP.NumberOfComponents;
	
	DEVICES
	v1 as valve_simplified;
	s1 as source;
	
	CONNECTIONS
	s1.Outlet to v1.Inlet;

	SPECIFY
	s1.P = 150 * 'kPa';
	s1.T = 281.75 * 'K';
	s1.Composition = [0.664, 0.336];

	v1.x = 0.866;
	v1.Outlet.P = 148 * 'kPa';

	SET
	v1.k = 20 * 'gal/min/psi^0.5';
	
	OPTIONS
	Dynamic = false;
	NLASolver(
		RelativeAccuracy = 1e-6
	);
end

FlowSheet valve_Test_L
#with liquid conditions
	PARAMETERS
	PP	as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = [ "isobutane", "benzene"],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp	as Integer;

	SET
	NComp = PP.NumberOfComponents;
	
	DEVICES
	v1 as valve_simplified;
	s1 as source;
	
	CONNECTIONS
	s1.Outlet to v1.Inlet;

	SPECIFY
	s1.P = 185 * 'kPa';
	s1.T = 328.12 * 'K';
	s1.Composition = [0.001848, 0.9982];

	v1.x = 0.3363;
	v1.Outlet.P = 183 * 'kPa';

	SET
	v1.k = 20 * 'gal/min/psi^0.5';
	
	OPTIONS
	Dynamic = false;
	NLASolver(
		RelativeAccuracy = 1e-6
	);
end