#*-------------------------------------------------------------------
* 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 - LMTD Method
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Gerson B. Bicca
* $Id: Sample_DoublePipe_convergence.mso $
*--------------------------------------------------------------------*#
using "heat_exchangers/DoublePipe";


#* Note: This sample file explain how to obtain a solution using a double pipe simplified model at conditions with better convergence characteristics

step 1:
	- use the option below to obtain a good approximation of the final solution:

	CalculationApproach = "Simplified";		setting the simplified model
	Qestimated = 70*'kW';							estimate heat duty for the simplified model (a good value to heat duty)

step 2:
	- Once you got convergence, use the EMSO graphical interface and save the result to a file (eg. result.rlt)

step 3:
	- CalculationApproach = "Full"; 				setting the full model
	- guessFile="result.rlt" 								use the option guessFile

*#

FlowSheet TESTE_1

PARAMETERS
	PP 			as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
									Components = ["aniline","benzene"]);
	NComp 		as Integer;
	
DEVICES

	exchanger 	as DoublePipe_LMTD;
	InletHot   		as simple_source; 
	InletCold   	as simple_source; 
	OutletHot 		as simple_sink;
	OutletCold 	as simple_sink;
	
	
CONNECTIONS

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

NComp   		= PP.NumberOfComponents;

#exchanger.CalculationApproach = "Simplified";	
exchanger.Qestimated = 30*'kW';

exchanger.CalculationApproach = "Full";

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

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

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

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

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 	Fouling
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
exchanger.Geometry.Rfi = 0.00018*'m^2*K/W';
exchanger.Geometry.Rfo = 0*'m^2*K/W';

SPECIFY
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 		Hot Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	#InletHot.Fw 			= 4396   	 * 'kg/h';
	InletHot.Outlet.F 			= 47.2038   	 * 'kmol/h';
	InletHot.Outlet.T 	= 338.71   * 'K';
	InletHot.Outlet.P 	= 413.685 * 'kPa';
	InletHot.Outlet.z 	= [1,0];
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#    	Cold Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	#InletCold.Fw 			= 4536   	* 'kg/h';
	InletCold.Outlet.F 	= 58.0694   	* 'kmol/h';
	InletCold.Outlet.T 	= 288.71     	* 'K';
	InletCold.Outlet.P 	= 413.685 * 'kPa';
    InletCold.Outlet.z 	= [0,1];
	
OPTIONS

Dynamic  = false;
GuessFile = "TESTE_1";

end

FlowSheet TESTE_2

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

	exchanger 	as DoublePipe_LMTD;
	InletHot   		as simple_source; 
	InletCold   	as simple_source; 
	OutletHot 		as simple_sink;
	OutletCold 	as simple_sink;
	
	
CONNECTIONS

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

NComp   		= PP.NumberOfComponents;

#exchanger.CalculationApproach = "Simplified";
exchanger.Qestimated = 40*'kW';

exchanger.CalculationApproach = "Full";

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

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

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

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

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 	Fouling
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
exchanger.Geometry.Rfi = 0.00018*'m^2*K/W';
exchanger.Geometry.Rfo = 0*'m^2*K/W';

SPECIFY
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 		Hot Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	#InletHot.Fw 			= 4396   	 * 'kg/h';
	InletHot.Outlet.F 			= 47.2038   	 * 'kmol/h';
	InletHot.Outlet.T 	= 338.71   * 'K';
	InletHot.Outlet.P 	= 413.685 * 'kPa';
	InletHot.Outlet.z 	= [1,0];
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#    	Cold Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	#InletCold.Fw 			= 4536   	* 'kg/h';
	InletCold.Outlet.F 	= 58.0694   	* 'kmol/h';
	InletCold.Outlet.T 	= 288.71     	* 'K';
	InletCold.Outlet.P 	= 413.685 * 'kPa';
    InletCold.Outlet.z 	= [0,1];
	
OPTIONS

Dynamic  = false;
GuessFile = "TESTE_2";

end

FlowSheet TESTE_3

PARAMETERS
	PP 			as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
									Components = [ "aniline", "toluene"]);
	NComp 		as Integer;
	
DEVICES

	exchanger 	as DoublePipe_LMTD;
	InletHot   		as simple_source; 
	InletCold   	as simple_source; 
	OutletHot 		as simple_sink;
	OutletCold 	as simple_sink;
	
	
CONNECTIONS

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

NComp   		= PP.NumberOfComponents;

#exchanger.CalculationApproach = "Simplified";
exchanger.Qestimated = 30*'kW';

exchanger.CalculationApproach = "Full";

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

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

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

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

#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 	Fouling
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
exchanger.Geometry.Rfi = 0.00018*'m^2*K/W';
exchanger.Geometry.Rfo = 0*'m^2*K/W';

SPECIFY
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# 		Hot Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	#InletHot.Fw 			= 4396   	 * 'kg/h';
	InletHot.Outlet.F 			= 47.2038   	 * 'kmol/h';
	InletHot.Outlet.T 	= 413   * 'K';
	InletHot.Outlet.P 	= 413.685 * 'kPa';
	InletHot.Outlet.z 	= [1,0];
	
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#    	Cold Stream
#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	#InletCold.Fw 			= 4536   	* 'kg/h';
	InletCold.Outlet.F 	= 58.0694   	* 'kmol/h';
	InletCold.Outlet.T 	= 330    	* 'K';
	InletCold.Outlet.P 	= 413.685 * 'kPa';
    InletCold.Outlet.z 	= [0,1];
	
OPTIONS

Dynamic  = false;
GuessFile = "TESTE_3";
end