#*-------------------------------------------------------------------
* 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 model E_Shell_NTU_Det
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO (www.vrtech.com.br) to run.
*
*----------------------------------------------------------------------
* Author: Gerson B. Bicca
* $Id: NTU_Method.mso 89 2006-12-11 15:26:11Z paula $
*--------------------------------------------------------------------*#
using "heat_exchangers/HeatExchangerDetailed.mso";

 FlowSheet NTU_Method

	
DEVICES

	exchanger 	as E_Shell_NTU_Det;
	streamhot 	as streamTP;
	streamcold 	as streamTP;

CONNECTIONS

	streamhot 				to exchanger.Inlet.Hot;
	streamcold 				to exchanger.Inlet.Cold;
	
PARAMETERS

	PP 			as CalcObject 	(File="vrpp");
	NComp 		as Integer;
	
SET
#===============================================================
#   Components - State Equation - Correlations
#===============================================================
	
	PP.LiquidModel    	= "RK";
	PP.VapourModel    	= "RK";
	PP.Components 		= ["water"];
	NComp   			= PP.NumberOfComponents;
	exchanger.HE.HotSide          = "Shell";
	
	exchanger.HE.TurbulentFlow     = "Petukhov";
#	exchanger.HE.TurbulentFlow     = "SiederTate";
	
#=====================================================================
# 	Shell Geometrical Parameters
#=====================================================================    
exchanger.Tpass			= 2;
exchanger.Dishell   	= 1.2  	*"m";
exchanger.Lcf       	= 0.007 *"m";
exchanger.Nss       	= 3;


#=====================================================================
# 	Tubes Geometrical Parameters
#=====================================================================
exchanger.Ntt         	= 360;
exchanger.Pattern   	= 30;
exchanger.pitch     	= 0.03048  	*"m";
exchanger.Ltube      	= 5.3  *"m";
exchanger.Ditube      	= 0.01397  *"m";
exchanger.Dotube    	= 0.025146  	*"m";
exchanger.Kwall   		= 0.5588  *"W/m/K";

#=====================================================================
# 	Baffles Geometrical Parameters
#=====================================================================
exchanger.Lcd       	= 5.5e-3 	*"m";
exchanger.Bc        	= 25;
exchanger.Ltd       	= 5e-4   	*"m";
exchanger.Nb   			= 5;
exchanger.Donozzle_Shell    = 0.15405  	*"m";
exchanger.Dinozzle_Shell    = 0.15405  	*"m";
exchanger.Honozzle_Shell    = 0.0225  	*"m";
exchanger.Hinozzle_Shell    = 0.0225	*"m";
exchanger.Donozzle_Tube     = 0.15405  	*"m";
exchanger.Dinozzle_Tube    	= 0.15405  	*"m";

SPECIFY
#============================================
#   Hot Streams
#============================================
	streamhot.F 					= 100    * "kmol/h";
	streamhot.T 					= 419    * "K";
	streamhot.P 					= 1     * "atm";
	streamhot.z 					= [1];

#============================================
#   Cold Streams
#============================================	
	streamcold.F  					= 50   * "kmol/h";
	streamcold.P  					= 1     * "atm";
	streamcold.T  					= 363   * "K";
    streamcold.z  					= [1];
#=====================================================================
# 	Baffles Spacing
#=====================================================================
	
	exchanger.Baffles.Ls	    	= 1   	*"m";
	exchanger.Baffles.Lsi	    	= 1   	*"m";
	

OPTIONS
#============================================
#   Simulation Options
#============================================
 mode 			= "steady";
  relativeAccuracy = 1e-8;

end