#*-------------------------------------------------------------------
* 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 of HeatExchangerSimplified_cost
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Núbia do Carmo Ferreira
* $Id: Sample_Simplified.mso 145 2007-01-30 18:20:00Z rafael $
*--------------------------------------------------------------------*#

using "costs/HeatExchangerSimplified_cost";


FlowSheet sample_HeatExchanger_LMTD_cost

DEVICES

	exchanger 	as HeatExchanger_LMTD_cost;
	streamhot 	as source;
	streamcold 	as source;


CONNECTIONS

	streamhot.Outlet 	to exchanger.InletHot;
	streamcold.Outlet 	to exchanger.InletCold;
	
PARAMETERS
	PP 			as Plugin 	(Type="PP",LiquidModel    = "RK",VapourModel   = "RK",
										Components 	= ["water"]);
	NComp 		as Integer;

SET

	NComp   			= PP.NumberOfComponents;
	
	exchanger.FlowDirection      = "counter";

	# costs
	exchanger.Material = "Stainless steel 316";
	
	exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
	exchanger.Cost(2,:) = [-0.9003,0.0906,0];
	exchanger.Cost(3,:) = [1.35,0,0];
	exchanger.Cost(4,:) = [-0.7844,0.0830,0];
	exchanger.Cost(5,:) = [0.8955,0.04981,0];
	exchanger.Cost(6,:) = [1.2002,0.07140,0];
	exchanger.Cost(7,:) = [1.4272,0.12088,0];
	exchanger.Cost(8,:) = [0.8608,0.23296,0];
	

SPECIFY

	exchanger.Method.Fc 					 = 1;
	
	exchanger.Details.A 					 = 139.35*'m^2';
	exchanger.Details.Ud 					 = 0.75*'W/(m^2*K)';
	exchanger.Details.Uc 					 = 0.75*'W/(m^2*K)';
	
	exchanger.HotSide.PressureDrop.Pdrop  = 0*'atm';
	exchanger.ColdSide.PressureDrop.Pdrop = 0*'atm';
	
	streamhot.Outlet.F 					= 1000    * 'kmol/h';
	streamhot.Outlet.T 					= 450     * 'K';
	streamhot.Outlet.P 					= 1.5     * 'atm';
	streamhot.Outlet.z 					= [1];

	streamcold.Outlet.F  					= 500   * 'kmol/h';
	streamcold.Outlet.P  					= 1     * 'atm';
	streamcold.Outlet.T  					= 315   * 'K';
    streamcold.Outlet.z  					= [1];
	
	
OPTIONS

 Dynamic = false;

end


FlowSheet sample_HeatExchanger_NTU_cost
	
DEVICES

	exchanger 	as HeatExchanger_NTU_cost;
	streamhot 	as source;
	streamcold 	as source;

CONNECTIONS

	streamhot.Outlet 	to exchanger.InletHot;
	streamcold.Outlet 	to exchanger.InletCold;
	
PARAMETERS
	PP 			as Plugin 	(Type="PP",LiquidModel    = "RK",VapourModel   = "RK",
										Components 	= ["water"]);
	NComp 		as Integer;
	
SET


	NComp   			= PP.NumberOfComponents;
	
	exchanger.FlowDirection      = "counter";

	# costs
	exchanger.Material = "Stainless steel 316";
	
	exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
	exchanger.Cost(2,:) = [-0.9003,0.0906,0];
	exchanger.Cost(3,:) = [1.35,0,0];
	exchanger.Cost(4,:) = [-0.7844,0.0830,0];
	exchanger.Cost(5,:) = [0.8955,0.04981,0];
	exchanger.Cost(6,:) = [1.2002,0.07140,0];
	exchanger.Cost(7,:) = [1.4272,0.12088,0];
	exchanger.Cost(8,:) = [0.8608,0.23296,0];

SPECIFY

	exchanger.Details.A 					 = 65.031*'m^2';
	exchanger.Details.Ud 					 = 0.75*'W/(m^2*K)';
	exchanger.Details.Uc 					 = 0.75*'W/(m^2*K)';
	
	exchanger.HotSide.PressureDrop.Pdrop  = 0*'atm';
	exchanger.ColdSide.PressureDrop.Pdrop = 0*'atm';
	
	streamhot.Outlet.F 					= 1000    * 'kmol/h';
	streamhot.Outlet.T 					= 450     * 'K';
	streamhot.Outlet.P 					= 1.5     * 'atm';
	streamhot.Outlet.z 					= [1];

	streamcold.Outlet.F  					= 500   * 'kmol/h';
	streamcold.Outlet.P  					= 1     * 'atm';
	streamcold.Outlet.T  					= 315   * 'K';
    streamcold.Outlet.z  					= [1];
	
OPTIONS

Dynamic = false;
GuessFile = "GuessSimplified";

end

FlowSheet sample_Shell_and_tubes_LMTD_cost

PARAMETERS

		PP 		as Plugin (Brief="Physical Properties",Type="PP",LiquidModel = "PR",VapourModel = "PR",
								Components 	= ["water","n-butane", "benzene" ]);
	NComp 	as Integer;
	
DEVICES

	exchanger 	 as Shell_and_tubes_LMTD_cost;
	streamhot 	 as source;
	streamcold   as source;
	OutCold      as sink;
	OutHot       as sink;

SET

	NComp   			= PP.NumberOfComponents;

	exchanger.LMTDcorrection = "Bowmann";
	exchanger.ShellType	        = "Eshell";
	
	# costs
	exchanger.Material = "Stainless steel 316";
	
	exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
	exchanger.Cost(2,:) = [-0.9003,0.0906,0];
	exchanger.Cost(3,:) = [1.35,0,0];
	exchanger.Cost(4,:) = [-0.7844,0.0830,0];
	exchanger.Cost(5,:) = [0.8955,0.04981,0];
	exchanger.Cost(6,:) = [1.2002,0.07140,0];
	exchanger.Cost(7,:) = [1.4272,0.12088,0];
	exchanger.Cost(8,:) = [0.8608,0.23296,0];


CONNECTIONS

	streamhot.Outlet 	to exchanger.InletHot;
	streamcold.Outlet 	to exchanger.InletCold;

	exchanger.OutletHot  to OutHot.Inlet;
	exchanger.OutletCold to OutCold.Inlet;

SPECIFY

exchanger.Details.Ud 						= 210*'W/(m^2*K)';
exchanger.Details.Uc 						= 210*'W/(m^2*K)';

exchanger.HotSide.PressureDrop.Pdrop 	= 0.1*'kPa';
exchanger.ColdSide.PressureDrop.Pdrop   	= 0.1*'kPa';

streamhot.Outlet.F		= 36    * 'kmol/h';
streamhot.Outlet.T 		= 363.15   * 'K';
streamhot.Outlet.P 		= 120   * 'kPa';
streamhot.Outlet.z 		= [1,0,0];

streamcold.Outlet.F 		= 10 * 'kmol/h';
streamcold.Outlet.P 		= 120 * 'kPa';
streamcold.Outlet.T 		= 303.15 * 'K';	
streamcold.Outlet.z 		= [0.1,0.7, 0.2];

exchanger.Details.A     = 15*'m^2';
	
OPTIONS

 Dynamic	= false;
 GuessFile   = "GuessShellandTubes";

end


FlowSheet sample_Shell_and_tubes_NTU_cost

PARAMETERS

	PP 		as Plugin (Brief="Physical Properties",Type="PP",LiquidModel = "PR",VapourModel = "PR",
								Components 	= ["water","n-butane", "benzene" ]);
	NComp 	as Integer;
	
DEVICES

	exchanger 	 as Shell_and_tubes_NTU_cost;
	streamhot 	 as source;
	streamcold  as source;
	OutCold      as sink;
	OutHot       as sink;

SET

	NComp   			= PP.NumberOfComponents;

	exchanger.ShellType	        = "Eshell";

	# costs
	exchanger.Material = "Stainless steel 316";
	
	exchanger.Cost(1,:) = [8.202,0.01506,0.06811];
	exchanger.Cost(2,:) = [-0.9003,0.0906,0];
	exchanger.Cost(3,:) = [1.35,0,0];
	exchanger.Cost(4,:) = [-0.7844,0.0830,0];
	exchanger.Cost(5,:) = [0.8955,0.04981,0];
	exchanger.Cost(6,:) = [1.2002,0.07140,0];
	exchanger.Cost(7,:) = [1.4272,0.12088,0];
	exchanger.Cost(8,:) = [0.8608,0.23296,0];

CONNECTIONS

	streamhot.Outlet 	to exchanger.InletHot;
	streamcold.Outlet 	to exchanger.InletCold;

	exchanger.OutletHot  to OutHot.Inlet;
	exchanger.OutletCold to OutCold.Inlet;

SPECIFY

exchanger.Details.Ud 						= 210*'W/(m^2*K)';
exchanger.Details.Uc 						= 210*'W/(m^2*K)';

exchanger.HotSide.PressureDrop.Pdrop 	= 0.1*'kPa';
exchanger.ColdSide.PressureDrop.Pdrop   	= 0.1*'kPa';

streamhot.Outlet.F		= 36    * 'kmol/h';
streamhot.Outlet.T 		= 363.15   * 'K';
streamhot.Outlet.P 		= 120   * 'kPa';
streamhot.Outlet.z 		= [1,0,0];

streamcold.Outlet.F 		= 10 * 'kmol/h';
streamcold.Outlet.P 		= 120 * 'kPa';
streamcold.Outlet.T 		= 303.15 * 'K';	
streamcold.Outlet.z 		= [0.1,0.7, 0.2];

exchanger.Details.A     = 15*'m^2';
	
OPTIONS

 Dynamic	= false;

end