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

using "heat_exchangers/DoublePipe.mso";

FlowSheet SeriesOfPipes

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

	exchanger1 	as DoublePipe_NTU;
	exchanger2 	as DoublePipe_LMTD;
	exchanger3 	as DoublePipe_NTU;
	
	InletHot   	as simple_source; 
	InletCold   	as simple_source; 
	OutletHot 	as sink;
	OutletCold 	as sink;
	
	
CONNECTIONS

	InletHot.Outlet					to exchanger1.InletOuter;
	InletCold.Outlet					to exchanger3.InletInner;
	
	exchanger3.OutletInner 	to exchanger2.InletInner;
	exchanger1.OutletOuter 	to exchanger2.InletOuter;
	
	exchanger2.OutletInner 	to exchanger1.InletInner;
	exchanger2.OutletOuter 	to exchanger3.InletOuter;
	
	exchanger1.OutletInner 	to OutletCold.Inlet;
	exchanger3.OutletOuter 	to OutletHot.Inlet;
	
SET

NComp   = PP.NumberOfComponents;

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#   Options for convergence
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

exchanger1.CalculationApproach = "Simplified";
exchanger1.Qestimated                   = 50*'kW';
exchanger1.Eftestimated                 = 0.64;

exchanger2.CalculationApproach = "Full";
exchanger2.Qestimated                   = 40*'kW';

exchanger3.CalculationApproach = "Full";
exchanger3.Qestimated                    = 30*'kW';
exchanger3.Eftestimated                  = 0.4;

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

exchanger1.FlowDirection  							= "counter";
exchanger1.InnerTurbulentCorrelation  	= "SiederTate";
exchanger1.OuterTurbulentCorrelation  	= "SiederTate";
exchanger1.OuterLaminarCorrelation     	= "Schlunder";

exchanger2.FlowDirection  							= "counter";
exchanger2.InnerTurbulentCorrelation  	= "SiederTate";
exchanger2.OuterTurbulentCorrelation  	= "SiederTate";
exchanger2.OuterLaminarCorrelation     	= "Schlunder";

exchanger3.FlowDirection  							= "counter";
exchanger3.InnerTurbulentCorrelation  	= "SiederTate";
exchanger3.OuterTurbulentCorrelation  	= "SiederTate";
exchanger3.OuterLaminarCorrelation     	= "Schlunder";
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# Double Pipe Geometrical Parameters and Alocation
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

exchanger1.Geometry.DoInner  	=	42.16*'mm';
exchanger1.Geometry.DiInner		=	35.05*'mm';
exchanger1.Geometry.DiOuter		=	52.50*'mm';
exchanger1.Geometry.Kwall   		= 0.057 *'kW/m/K';
exchanger1.Geometry.Lpipe			=	9*'m';

exchanger2.Geometry.DoInner  	=	42.16*'mm';
exchanger2.Geometry.DiInner		=	35.05*'mm';
exchanger2.Geometry.DiOuter		=	52.50*'mm';
exchanger2.Geometry.Kwall   		= 0.057 *'kW/m/K';
exchanger2.Geometry.Lpipe			=	9*'m';

exchanger3.Geometry.DoInner  	=	42.16*'mm';
exchanger3.Geometry.DiInner		=	35.05*'mm';
exchanger3.Geometry.DiOuter		=	52.50*'mm';
exchanger3.Geometry.Kwall   		= 0.057 *'kW/m/K';
exchanger3.Geometry.Lpipe			= 9*'m';

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 	Fouling
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
exchanger1.Geometry.Rfi = 0.00217*'m^2*K/kW';
exchanger1.Geometry.Rfo = 0.0008*'m^2*K/kW';

exchanger2.Geometry.Rfi = 0.00217*'m^2*K/kW';
exchanger2.Geometry.Rfo = 0.0008*'m^2*K/kW';

exchanger3.Geometry.Rfi = 0.00217*'m^2*K/kW';
exchanger3.Geometry.Rfo = 0.0008*'m^2*K/kW';

SPECIFY
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 		Hot Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	InletHot.Outlet.F 			= 51   			* 'kmol/h';
	InletHot.Outlet.T 	= (90+273.15)   * 'K';
	InletHot.Outlet.P 	= 2 			* 'atm';
	InletHot.Outlet.z 	= [1,0];
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#    	Cold Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	InletCold.Outlet.F 			= 251   			* 'kmol/h';
	InletCold.Outlet.P = 2     		* 'atm';
	InletCold.Outlet.T = (33+273.15)	* 'K';
    InletCold.Outlet.z = [0,1];
	
OPTIONS

Dynamic  = false;
GuessFile = "SeriesOfPipes";

end