#*-------------------------------------------------------------------
* 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.
*
*--------------------------------------------------------------------
* Hairpin Heat Exchanger - NTU and LMTD Method
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Gerson B. Bicca
* $Id: Sample_hairpin.mso  $
*--------------------------------------------------------------------*#

using "heat_exchangers/Hairpin.mso";

FlowSheet LMTD_Method

PARAMETERS
	PP 			as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
									Components = ["n-hexane","water"]);
	NComp 		as Integer;
	
DEVICES

	exchanger 	as Hairpin_LMTD;
	InletHot   	as simple_source; 
	InletCold   	as simple_source; 
	OutletHot 	as sink;
	OutletCold 	as sink;
	
	
CONNECTIONS

	InletHot.Outlet			to exchanger.InletInner;
	InletCold.Outlet			to exchanger.InletOuter;
	exchanger.OutletOuter 	to OutletCold.Inlet;
	exchanger.OutletInner 	to OutletHot.Inlet;
	
SET

NComp   = PP.NumberOfComponents;

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#   Options
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

exchanger.FlowDirection  					= "counter";
exchanger.InnerTurbulentCorrelation  	= "SiederTate";

exchanger.OuterTurbulentCorrelation  	= "SiederTate";
exchanger.OuterLaminarCorrelation     	= "Schlunder";
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# Hairpin Geometrical Parameters and Alocation
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

exchanger.DoInner  		=	0.04826*'m';
exchanger.DiInner		=	0.04089*'m';
exchanger.DiOuter		=	0.07793*'m';
exchanger.Kwall   		= 	0.057 *'kW/m/K';
exchanger.Lpipe			=	3*'m';

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 	Fouling
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
exchanger.Rfi = 0.0018*'m^2*K/kW';
exchanger.Rfo = 0.0021*'m^2*K/kW';

SPECIFY
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 		Hot Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	InletHot.Outlet.F = 10   			* 'kmol/h';
	InletHot.Outlet.T = (68+273.15)   * 'K';
	InletHot.Outlet.P = 5.4 			* 'bar';
	InletHot.Outlet.z = [1,0];
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#    	Cold Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	InletCold.Outlet.F = 5   			* 'kmol/h';
	InletCold.Outlet.P = 5.1     		* 'bar';
	InletCold.Outlet.T = (23+273.15)	* 'K';
    InletCold.Outlet.z = [0,1];
	
OPTIONS

Dynamic  = false;
GuessFile = "GuessPipe";

end

FlowSheet NTU_Method

PARAMETERS

	PP 			as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
									Components = ["n-hexane","water"]);
	NComp 		as Integer;
	
DEVICES

	exchanger 	as Hairpin_NTU;
	InletHot   	as simple_source; 
	InletCold   	as simple_source; 
	OutletHot 	as sink;
	OutletCold 	as sink;
	
	
CONNECTIONS

	InletHot.Outlet			to exchanger.InletInner;
	InletCold.Outlet			to exchanger.InletOuter;
	exchanger.OutletOuter 	to OutletCold.Inlet;
	exchanger.OutletInner 	to OutletHot.Inlet;
	
SET

NComp   = PP.NumberOfComponents;

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#   Options
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

exchanger.FlowDirection  = "counter";
exchanger.InnerTurbulentCorrelation  = "SiederTate";

exchanger.OuterTurbulentCorrelation  = "SiederTate";
exchanger.OuterLaminarCorrelation     = "Schlunder";
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# Hairpin Geometrical Parameters and Alocation
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

exchanger.DoInner  		=	0.04826*'m';
exchanger.DiInner		=	0.04089*'m';
exchanger.DiOuter		=	0.07793*'m';
exchanger.Kwall   		= 	0.057 *'kW/m/K';
exchanger.Lpipe			=	3*'m';

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 	Fouling
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
exchanger.Rfi = 0.0018*'m^2*K/kW';
exchanger.Rfo = 0.0021*'m^2*K/kW';

SPECIFY
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 		Hot Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	InletHot.Outlet.F = 10   			* 'kmol/h';
	InletHot.Outlet.T = (68+273.15)   * 'K';
	InletHot.Outlet.P = 5.4 			* 'bar';
	InletHot.Outlet.z = [1,0];
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#    	Cold Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	InletCold.Outlet.F = 5   			* 'kmol/h';
	InletCold.Outlet.P = 5.1     		* 'bar';
	InletCold.Outlet.T = (23+273.15)	* 'K';
    InletCold.Outlet.z = [0,1];

OPTIONS

Dynamic = false;

GuessFile = "GuessPipe";

end