#*-------------------------------------------------------------------
* 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.
*
*--------------------------------------------------------------------
* Author: Gerson Balbueno Bicca 
* $Id: HeatExchangerDiscretized.mso 139 2007-01-29 15:50:41Z bicca $
*------------------------------------------------------------------*#

using "HEX_Engine";

Model HeatExchangerDiscretized_Basic 

ATTRIBUTES
	Pallete 	= false;
	Brief 		= "Basic Model for Discretized Heat Exchangers";
	Info 		=
	"write some information";
	
PARAMETERS

outer PP 	    as Plugin(Brief="External Physical Properties");
	HE 	    as Plugin(Brief="STHE Calculations",File="heatex");
outer NComp   as Integer   (Brief="Number of Components");
  M(NComp)  as molweight (Brief="Component Mol Weight");
	
VARIABLES

in  Inlet	    as Inlet_Main_Stream; 	# Hot and Cold Inlets
out Outlet      as Outlet_Main_Stream;	# Hot and Cold Outlets
	Properties  as Main_Properties;		# Hot and Cold Properties
	Details     as Details_Main;
	Tubes   	as Tube_Side_Main;	
	Shell   	as Shell_Side_Main; 
	Resistances as Main_Resistances;
	Baffles     as Baffles_Main;
	
SET

M   = PP.MolecularWeight();

EQUATIONS

"Hot Stream Average Temperature"
	Properties.Hot.Average.T = 0.5*Inlet.Hot.T + 0.5*Outlet.Hot.T;
	
"Cold Stream Average Temperature"
	Properties.Cold.Average.T = 0.5*Inlet.Cold.T + 0.5*Outlet.Cold.T;
	
"Hot Stream Average Pressure"
	Properties.Hot.Average.P = 0.5*Inlet.Hot.P+0.5*Outlet.Hot.P;
	
"Cold Stream Average Pressure"
	Properties.Cold.Average.P = 0.5*Inlet.Cold.P+0.5*Outlet.Cold.P;
	
"Hot Stream Average Molecular Weight"
	Properties.Hot.Average.Mw = sum(M*Inlet.Hot.z);

"Cold Stream Average Molecular Weight"
	Properties.Cold.Average.Mw = sum(M*Inlet.Cold.z);
	

if Inlet.Cold.v equal 0
	
	then	
	
"Cold Stream Average Heat Capacity"
	Properties.Cold.Average.Cp 	= 	PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Heat Capacity"
	Properties.Cold.Inlet.Cp 	= 	PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Heat Capacity"
	Properties.Cold.Outlet.Cp 	= 	PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Average Mass Density"
	Properties.Cold.Average.rho = 	PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Mass Density"
	Properties.Cold.Inlet.rho 	= 	PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Mass Density"
	Properties.Cold.Outlet.rho 	= 	PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Average Viscosity"
	Properties.Cold.Average.Mu 	= 	PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream inlet Viscosity"
	Properties.Cold.Inlet.Mu 	= 	PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
	
"Cold Stream Outlet Viscosity"
	Properties.Cold.Outlet.Mu 	= 	PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Average Conductivity"
	Properties.Cold.Average.K 	= 	PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Conductivity"	
	Properties.Cold.Inlet.K 	= 	PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Conductivity"
	Properties.Cold.Outlet.K 	= 	PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Heat Capacity at Wall Temperature"
	Properties.Cold.Wall.Cp 	= 	PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
	
"Cold Stream Viscosity at Wall Temperature"
	Properties.Cold.Wall.Mu 	= 	PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Conductivity at Wall Temperature"
	Properties.Cold.Wall.K 		= 	PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);


	else

"Cold Stream Average Heat Capacity"
	Properties.Cold.Average.Cp 	=	PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Heat Capacity"	
	Properties.Cold.Inlet.Cp 	=	PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Heat Capacity"	
	Properties.Cold.Outlet.Cp 	=	PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Average Mass Density"
	Properties.Cold.Average.rho =	PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Mass Density"
	Properties.Cold.Inlet.rho 	=	PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Mass Density"	
	Properties.Cold.Outlet.rho 	=	PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Average Viscosity "
	Properties.Cold.Average.Mu 	=	PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Viscosity "	
	Properties.Cold.Inlet.Mu 	=	PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Viscosity "	
	Properties.Cold.Outlet.Mu 	=	PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);

"Cold Stream Average Conductivity "
	Properties.Cold.Average.K 	= 	PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Inlet Conductivity "
	Properties.Cold.Inlet.K 	= 	PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);

"Cold Stream Outlet Conductivity "
	Properties.Cold.Outlet.K 	= 	PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
	
"Cold Stream Heat Capacity at Wall Temperature"
	Properties.Cold.Wall.Cp 	=	PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);


"Cold Stream Viscosity at Wall Temperature"
	Properties.Cold.Wall.Mu 	=	PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);

"Cold Stream Conductivity at Wall Temperature"
	Properties.Cold.Wall.K 		= 	PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
	
	
	
end


if Inlet.Hot.v equal 0

	then

"Hot Stream Average Heat Capacity"
	Properties.Hot.Average.Cp 	=		PP.LiquidCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Inlet Heat Capacity"
	Properties.Hot.Inlet.Cp 	=		PP.LiquidCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);

"Hot Stream Outlet Heat Capacity"
	Properties.Hot.Outlet.Cp 	=		PP.LiquidCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);

"Hot Stream Average Mass Density"
	Properties.Hot.Average.rho 	=		PP.LiquidDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Inlet Mass Density"	
	Properties.Hot.Inlet.rho 	=		PP.LiquidDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);

"Hot Stream Outlet Mass Density"	
	Properties.Hot.Outlet.rho 	=		PP.LiquidDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);

"Hot Stream Average Viscosity"
	Properties.Hot.Average.Mu 	=		PP.LiquidViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);	

"Hot Stream Inlet Viscosity"
	Properties.Hot.Inlet.Mu 	=		PP.LiquidViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);	

"Hot Stream Outlet Viscosity"
	Properties.Hot.Outlet.Mu 	=		PP.LiquidViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);	

"Hot Stream Average Conductivity"
	Properties.Hot.Average.K 	=		PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);	

"Hot Stream Inlet Conductivity"
	Properties.Hot.Inlet.K 		=		PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);	

"Hot Stream Outlet Conductivity"
	Properties.Hot.Outlet.K 	=		PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);	

"Hot Stream Heat Capacity at Wall Temperature"
	Properties.Hot.Wall.Cp 		=		PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Viscosity  at Wall Temperature"
	Properties.Hot.Wall.Mu 		=		PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);	

"Hot Stream Conductivity at Wall Temperature"
	Properties.Hot.Wall.K 		=		PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);	
	

	else

"Hot Stream Average Heat Capacity"
	Properties.Hot.Average.Cp 	=		PP.VapourCp(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Inlet Heat Capacity"
	Properties.Hot.Inlet.Cp 	=		PP.VapourCp(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);

"Hot Stream Outlet Heat Capacity"
	Properties.Hot.Outlet.Cp 	=		PP.VapourCp(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);

"Hot Stream Average Mass Density"
	Properties.Hot.Average.rho 	=		PP.VapourDensity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Inlet Mass Density"	
	Properties.Hot.Inlet.rho 	=		PP.VapourDensity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);

"Hot Stream Outlet Mass Density"
	Properties.Hot.Outlet.rho 	=		PP.VapourDensity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);

"Hot Stream Average Viscosity"
	Properties.Hot.Average.Mu 	=		PP.VapourViscosity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Inlet Viscosity"
	Properties.Hot.Inlet.Mu 	=		PP.VapourViscosity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);

"Hot Stream Outlet Viscosity"
	Properties.Hot.Outlet.Mu 	=		PP.VapourViscosity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);

"Hot Stream Average Conductivity"
	Properties.Hot.Average.K 	=		PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);	

"Hot Stream Inlet Conductivity"
	Properties.Hot.Inlet.K 		=		PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);	
	
"Hot Stream Outlet Conductivity"
	Properties.Hot.Outlet.K 	=		PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);	

"Hot Stream Heat Capacity at Wall Temperature"
	Properties.Hot.Wall.Cp 		=		PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Viscosity at Wall Temperature"
	Properties.Hot.Wall.Mu 		=		PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);

"Hot Stream Conductivity at Wall Temperature"
	Properties.Hot.Wall.K 		=		PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);	


end

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Thermal Details
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
"Hot Stream Heat Capacity"
	Details.Ch =Inlet.Hot.F*Properties.Hot.Average.Cp; 
	
"Cold Stream Heat Capacity"
	Details.Cc =Inlet.Cold.F*Properties.Cold.Average.Cp;
	
"Minimum Heat Capacity"
	Details.Cmin  = min([Details.Ch,Details.Cc]);

"Maximum Heat Capacity"
	Details.Cmax  = max([Details.Ch,Details.Cc]);

"Heat Capacity Ratio"
	Details.Cr    = Details.Cmin/Details.Cmax;
	
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Energy Balance
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
"Energy Balance Hot Stream"
	Details.Q = Inlet.Hot.F*(Inlet.Hot.h-Outlet.Hot.h);

"Energy Balance Cold Stream"
	Details.Q =-Inlet.Cold.F*(Inlet.Cold.h-Outlet.Cold.h);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Material Balance
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
"Flow Mass Inlet Cold Stream"
	Properties.Cold.Inlet.Fw	=  sum(M*Inlet.Cold.z)*Inlet.Cold.F;

"Flow Mass Outlet Cold Stream"
	Properties.Cold.Outlet.Fw	=  sum(M*Outlet.Cold.z)*Outlet.Cold.F;

"Flow Mass Inlet Hot Stream"
	Properties.Hot.Inlet.Fw		=  sum(M*Inlet.Hot.z)*Inlet.Hot.F;

"Flow Mass Outlet Hot Stream"	
	Properties.Hot.Outlet.Fw	=  sum(M*Outlet.Hot.z)*Outlet.Hot.F;

"Molar Balance Hot Stream"
	Inlet.Hot.F  = Outlet.Hot.F;
	
"Molar Balance Cold Stream"
	Inlet.Cold.F = Outlet.Cold.F;

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Constraints
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
"Hot Stream Molar Fraction Constraint"
	Outlet.Hot.z=Inlet.Hot.z;
	
"Cold Stream Molar Fraction Constraint"
	Outlet.Cold.z=Inlet.Cold.z;
	
"No Phase Change In Cold Stream"
	Inlet.Cold.v=Outlet.Cold.v;

"No Phase Change In Hot Stream"
	Inlet.Hot.v=Outlet.Hot.v;

	
end

Model Heatex_Discretized_NTU 		as HeatExchangerDiscretized_Basic

ATTRIBUTES
	Pallete 	= false;
	Brief 		= "write some information";
	Info 		=
	"write some information";
	
VARIABLES

Eft	  	as positive (Brief="Effectiveness",Default=0.05,Lower=1e-8);
	
EQUATIONS

"Exchange Surface Area"
	Details.Q	= Eft*Details.Cmin*(Inlet.Hot.T-Inlet.Cold.T);	

"TEMA E Shell Effectiveness"
	Eft = HE.EshellEffectiveness(Details.Cr,Details.NTU);

end

Model Heatex_Discretized_LMTD 		as HeatExchangerDiscretized_Basic

ATTRIBUTES
	Pallete 	= false;
	Brief 		= "write some information";
	Info 		=
	"write some information";
	
VARIABLES
DT0     as temp_delta	(Brief="Temperature Difference at Inlet",Lower=1);
DTL		as temp_delta	(Brief="Temperature Difference at Outlet",Lower=1);
LMTD	as temp_delta	(Brief="Logarithmic Mean Temperature Difference",Lower=5);
Fc		as positive		(Brief="LMTD Correction Factor",Lower=0.5);
MTD		as temp_delta	(Brief="Mean Temperature Difference",Lower=1);

EQUATIONS
"Exchange Surface Area"
	Details.Q   = Details.Ud*Details.A*MTD;
	
"Mean Temperature Difference"	
	MTD   = Fc*LMTD;
	
"LMTD Correction Factor"
	Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);

"Temperature Difference at Inlet"
	DT0 = Inlet.Hot.T - Outlet.Cold.T;

"Temperature Difference at Outlet"
	DTL = Outlet.Hot.T - Inlet.Cold.T;
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 			Log Mean Temperature Difference
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL])) 
	
	then
"Log Mean Temperature Difference"
	LMTD= (DT0-DTL)/ln(DT0/DTL);

	else
	
if DT0*DTL equal 0
	
	then
"Log Mean Temperature Difference"
	LMTD = 0.5*(DT0+DTL);
	
	else
"Log Mean Temperature Difference"
	LMTD = 0.5*(DT0+DTL)*(1-(DT0-DTL)^2/(DT0*DTL)*(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
	
end
	
end

	
end

Model Profiles

ATTRIBUTES
	Pallete 	= false;
	Brief 		= "write some information";
	Info 		=
	"write some information";
	
PARAMETERS

Zones 	as Integer(Brief="Number of Zones");

VARIABLES

Lz(Zones)       	as length			(Brief="Zone Tube Length");
Area       			as area				(Brief="Area Total");
Q          			as power			(Brief="Total Duty");
PdropTubesNozzle 	as pressure			(Brief="Total Tube Nozzles Pressure Drop");
PdropTubes 			as pressure			(Brief="Total Tube Pressure Drop");
PdropWin		   	as pressure			(Brief="Total Shell Side Window Pressure Drop");
PdropCross		 	as pressure			(Brief="Total Shell Side Cross Flow Pressure Drop");
PdropEnds		  	as pressure			(Brief="Total Shell Side Ends Pressure Drop");
PdropShellNozzle 	as pressure			(Brief="Total Shell Side Nozzles Pressure Drop");
PdropShell		 	as pressure			(Brief="Total Shell Side Pressure Drop");
Udaverage 			as heat_trans_coeff (Brief="Average Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
Ucaverage 			as heat_trans_coeff (Brief="Average Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
hshellaverage 		as heat_trans_coeff	(Brief="Average Shell Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
htubeaverage    	as heat_trans_coeff	(Brief="Average Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
Thot(Zones)     	as temperature		(Brief="Hot Temperature",Lower = 300, Upper = 500);
Tcold(Zones)    	as temperature		(Brief="Cold Temperature",Lower = 300, Upper = 500);
Phot(Zones)     	as pressure			(Brief="Hot Pressure",Lower = 0.8, Upper = 30);
Pcold(Zones)    	as pressure			(Brief="Cold Pressure",Lower = 0.8, Upper = 30);

end

Model E_Shell_NTU_Disc      	 

ATTRIBUTES
	Pallete 	= false;
	Brief 		= "Shell and Tubes Heat Exchanger with 1 shell pass - NTU Method";
	Info 		=
	"write some information";
	
PARAMETERS

	HE 	as Plugin 	(Brief="STHE Calculations",File="heatex");
outer PP 	as Plugin	(Brief="External Physical Properties");
side 	as Integer 		(Brief="Fluid Alocation",Lower=0,Upper=1);
Pi		as constant     (Brief="Pi Number",Default=3.14159265);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Shell Geometrical Parameters
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
Tpass   		as Integer		(Brief="Number of Tube Passes",Lower=1);
Nss				as Integer		(Brief="Number of Sealing Strips pairs",Lower=1);
Dishell 		as length		(Brief="Inside Shell Diameter",Lower=10e-6);
Donozzle_Shell  as length		(Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
Dinozzle_Shell  as length		(Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
Aonozzle_Shell  as area			(Brief="Shell Outlet Nozzle Area",Lower=10e-6);
Ainozzle_Shell  as area			(Brief="Shell Inlet Nozzle Area",Lower=10e-6);
Aeonozzle_Shell as area			(Brief="Shell Outlet Escape Area Under Nozzle",Lower=10e-6);
Aeinozzle_Shell as area			(Brief="Shell Inlet Escape Area Under Nozzle",Lower=10e-6);
Hinozzle_Shell  as length		(Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
Honozzle_Shell  as length		(Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
Lcf 			as length		(Brief="Bundle-to-Shell Clearance",Lower=10e-8);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
#				Tubes Geometrical Parameters						#
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
Ntt			  as Integer		(Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
Pattern       as Integer		(Brief="Tube Layout Characteristic Angle",Lower=30);
Ltube		  as length			(Brief="Effective Tube Length",Lower=0.1);
pitch 		  as length			(Brief="Tube Pitch",Lower=1e-8);
Kwall		  as conductivity 	(Brief="Tube Wall Material Thermal Conductivity");
Dotube  	  as length			(Brief="Tube Outside Diameter",Lower=10e-6);
Ditube 		  as length			(Brief="Tube Inside Diameter",Lower=10e-6);
Donozzle_Tube as length			(Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
Dinozzle_Tube as length			(Brief="Tube Inlet Nozzle Diameter",Lower=10e-6);
Aonozzle_Tube as area			(Brief="Tube Outlet Nozzle Area",Lower=10e-6);
Ainozzle_Tube as area			(Brief="Tube Inlet Nozzle Area",Lower=10e-6);
Kinlet_Tube   as positive		(Brief="Tube Inlet Nozzle Pressure Loss Coeff",Default=1.1);
Koutlet_Tube  as positive		(Brief="Tube Outlet Nozzle Pressure Loss Coeff",Default=0.7);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Baffles Geometrical Parameters
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
Bc      		as Integer		(Brief="Baffle Cut",Default=25,Lower=25);
Nb          	as Integer		(Brief="Number of Baffles",Default=4);
Lcd     		as length		(Brief="Baffle-to-Shell Clearance",Lower=10e-8);
Ltd				as length		(Brief="Tube-to-Bafflehole Clearance",Lower=10e-8);
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
VARIABLES

Unity(Nb+1) 	as Heatex_Discretized_NTU; 
Sumary          as Profiles;

CONNECTIONS

Unity([1:Nb]).Outlet.Hot  to Unity([2:Nb+1]).Inlet.Hot;
Unity([2:Nb+1]).Outlet.Cold to Unity([1:Nb]).Inlet.Cold;

EQUATIONS

"Hot Temperatures"
	Sumary.Thot 	= Unity.Outlet.Hot.T;
	
"Cold Temperatures"
	Sumary.Tcold 	= Unity.Outlet.Cold.T ;

"Hot Pressures"
	Sumary.Phot 	= Unity.Outlet.Hot.P ;
	
"Cold Pressures"
	Sumary.Pcold 	= Unity.Outlet.Cold.P ;

"Average Shell Side Film Coefficient"
	Sumary.hshellaverage       	= sum(Unity.Shell.HeatTransfer.hshell)/Sumary.Zones;

"Average Tube Side Film Coefficient"
	Sumary.htubeaverage        	= sum(Unity.Tubes.HeatTransfer.htube)/Sumary.Zones;

"Average Overall Heat Transfer Coefficient Dirty"
	Sumary.Udaverage       		= sum(Unity.Details.Ud)/Sumary.Zones;
	
"Average Overall Heat Transfer Coefficient Clean"
	Sumary.Ucaverage       		= sum(Unity.Details.Uc)/Sumary.Zones;

"Area Total"
	Sumary.Area 				= sum(Unity.Details.A);

"Duty Total"
	Sumary.Q    				= sum(Unity.Details.Q);

"Length Inlet zone"
	Sumary.Lz(1)    			= Unity(1).Baffles.Lsi;

"Length Outlet zone"
	Sumary.Lz(Nb+1) 			= Unity(1).Baffles.Lso;

"Total Shell Side Pressure Drop"
	Sumary.PdropShell 			= sum(Unity.Shell.PressureDrop.Pdtotal);

"Total Tube Side Pressure Drop"
	Sumary.PdropTubes 			= sum(Unity.Tubes.PressureDrop.Pdtotal);

"Total Tube Side Nozzles Pressure Drop"
	Sumary.PdropTubesNozzle		= sum(Unity.Tubes.PressureDrop.Pdnozzle_in + Unity.Tubes.PressureDrop.Pdnozzle_out);

"Total Shell Side Nozzles Pressure Drop"
	Sumary.PdropShellNozzle		= sum(Unity.Shell.PressureDrop.Pdnozzle_in + Unity.Shell.PressureDrop.Pdnozzle_out);

"Total Shell Side Window Pressure Drop"
	Sumary.PdropWin				= sum(Unity.Shell.PressureDrop.Pdwindow);

"Total Shell Side Cross Flow Pressure Drop"
	Sumary.PdropCross 			= sum(Unity.Shell.PressureDrop.PdCross);

"Total Shell Side Ends Pressure Drop"
	Sumary.PdropEnds			= sum(Unity.Shell.PressureDrop.PdEndZones);

"Shell Side Cross Flow Area" 
	Unity(1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(1).Baffles.Lsi);
	
"Shell Side Cross Flow Area" 
	Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso);


if side equal 1
	
	then
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(1).Shell.PressureDrop.RVsquare_in = Unity(1).Properties.Hot.Inlet.rho*(Unity(1).Shell.PressureDrop.Vnozzle_in)^2;

"Shell Side Outlet Nozzle rho-V^2"
	Unity(Nb+1).Shell.PressureDrop.RVsquare_out = Unity(Nb+1).Properties.Hot.Outlet.rho*(Unity(Nb+1).Shell.PressureDrop.Vnozzle_out)^2;

"Shell Pressure End Zones"
	Unity(1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(1).Shell.HeatTransfer.Re,Unity(1).Baffles.Ls,Unity(1).Baffles.Lso,
	Unity(1).Baffles.Lsi,Unity(1).Properties.Hot.Inlet.Fw,Unity(1).Shell.HeatTransfer.Phi,Unity(1).Properties.Hot.Average.rho);

"Shell Pressure End Zones"
	Unity(Nb+1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(Nb+1).Shell.HeatTransfer.Re,Unity(Nb+1).Baffles.Ls,Unity(Nb+1).Baffles.Lso,
	Unity(Nb+1).Baffles.Lsi,Unity(Nb+1).Properties.Hot.Inlet.Fw,Unity(Nb+1).Shell.HeatTransfer.Phi,Unity(Nb+1).Properties.Hot.Average.rho);

"Pressure Drop Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Pdnozzle_in	= 0.5*Kinlet_Tube*Unity(1).Properties.Cold.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;

"Velocity Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Vnozzle_in	= Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Tube);

"Pressure Drop Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out	= 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Cold.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;

"Velocity Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out	= Unity(Nb+1).Properties.Cold.Inlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Aonozzle_Tube);

"Shell Pressure Drop Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Pdnozzle_in	= (0.5*Unity(1).Properties.Hot.Inlet.Fw^2/Unity(1).Properties.Hot.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Vnozzle_in	= Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Shell);

"Shell Pressure Drop Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out	= (0.5*Unity(Nb+1).Properties.Hot.Outlet.Fw^2/Unity(Nb+1).Properties.Hot.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Vnozzle_out	= Unity(Nb+1).Properties.Hot.Outlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Shell);


	else
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(1).Shell.PressureDrop.RVsquare_in = Unity(1).Properties.Cold.Inlet.rho*(Unity(1).Shell.PressureDrop.Vnozzle_in)^2;

"Shell Side Outlet Nozzle rho-V^2"
	Unity(Nb+1).Shell.PressureDrop.RVsquare_out = Unity(Nb+1).Properties.Cold.Outlet.rho*(Unity(Nb+1).Shell.PressureDrop.Vnozzle_out)^2;

"Shell Pressure End Zones"
	Unity(1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(1).Shell.HeatTransfer.Re,Unity(1).Baffles.Ls,Unity(1).Baffles.Lso,
	Unity(1).Baffles.Lsi,Unity(1).Properties.Cold.Inlet.Fw,Unity(1).Shell.HeatTransfer.Phi,Unity(1).Properties.Cold.Average.rho);

"Shell Pressure End Zones"
	Unity(Nb+1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(Nb+1).Shell.HeatTransfer.Re,Unity(Nb+1).Baffles.Ls,Unity(Nb+1).Baffles.Lso,
	Unity(Nb+1).Baffles.Lsi,Unity(Nb+1).Properties.Cold.Inlet.Fw,Unity(Nb+1).Shell.HeatTransfer.Phi,Unity(Nb+1).Properties.Cold.Average.rho);

"Pressure Drop Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Pdnozzle_in	= 0.5*Kinlet_Tube*Unity(1).Properties.Hot.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;
	
"Velocity Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Vnozzle_in 	= Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Tube);
	
"Pressure Drop Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out	= 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Hot.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;
	
"Velocity Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out	= Unity(Nb+1).Properties.Hot.Inlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Tube);	

"Shell Pressure Drop Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Pdnozzle_in	= (0.5*Unity(1).Properties.Cold.Inlet.Fw^2/Unity(1).Properties.Cold.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Vnozzle_in	= Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Shell);

"Shell Pressure Drop Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out = (0.5*Unity(Nb+1).Properties.Cold.Outlet.Fw^2/Unity(Nb+1).Properties.Cold.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Vnozzle_out  = Unity(Nb+1).Properties.Cold.Outlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Ainozzle_Shell);

end


for i in [2:Nb]

"Pressure Drop Tube Side Inlet Nozzle"
	Unity(i).Tubes.PressureDrop.Pdnozzle_in	= 0;

"Velocity Tube Side Inlet Nozzle"
	Unity(i).Tubes.PressureDrop.Vnozzle_in 	= 0;

"Pressure Drop Tube Side Outlet Nozzle"
	Unity(i).Tubes.PressureDrop.Pdnozzle_out	= 0;	

"Velocity Tube Side Outlet Nozzle"
	Unity(i).Tubes.PressureDrop.Vnozzle_out	= 0;	

"Shell Pressure Drop Inlet Nozzle"
	Unity(i).Shell.PressureDrop.Pdnozzle_in	= 0;

"Velocity Shell Side Inlet Nozzle"
	Unity(i).Shell.PressureDrop.Vnozzle_in	= 0;

"Shell Pressure Drop Outlet Nozzle"
	Unity(i).Shell.PressureDrop.Pdnozzle_out = 0;

"Velocity Shell Side Outlet Nozzle"
	Unity(i).Shell.PressureDrop.Vnozzle_out  = 0;

"Shell Pressure End Zones"
	Unity(i).Shell.PressureDrop.PdEndZones  	= 	0;
	
"Shell Side Outlet Nozzle rho-V^2"
	Unity(i).Shell.PressureDrop.RVsquare_out = 0;
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(i).Shell.PressureDrop.RVsquare_in = 0;

if side equal 1
	
	then
	
"Shell Pressure Drop Cross Flow"
	Unity(i).Shell.PressureDrop.PdCross 		= 	HE.DeltaPcrossIncremental(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Hot.Average.rho);


	else

"Shell Pressure Drop Cross Flow"
	Unity(i).Shell.PressureDrop.PdCross     = 	HE.DeltaPcrossIncremental(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Cold.Average.rho);


end

"Zone Length"
	Sumary.Lz(i) = Unity(1).Baffles.Ls;

"Shell Side Cross Flow Area" 
	Unity(i).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(i).Baffles.Ls);
	
end


for i in [1:Nb+1]
	
if side equal 1	
	
	then 
"Pressure Drop Hot Stream"
	Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P - Unity(i).Shell.PressureDrop.Pdtotal;	

"Pressure Drop Cold Stream"
	Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Tubes.PressureDrop.Pdtotal;
	
"Shell Side Reynolds Number"
	Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;

"Shell Heat Transfer Coefficient"
	Unity(i).Shell.HeatTransfer.hshell =Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;


"Shell Pressure Baffle Window"
	Unity(i).Shell.PressureDrop.Pdwindow		= 	HE.DeltaPwindowIncremental(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.rho,Unity(i).Properties.Hot.Average.Mu,Unity(i).Baffles.Ls);

"Hot Wall Temperature"
	Unity(i).Properties.Hot.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;

"Cold Wall Temperature"	
	Unity(i).Properties.Cold.Wall.Twall = 	(Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;

"Tube Side Velocity"
	Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);

"Tube Side Reynolds Number"
	Unity(i).Tubes.HeatTransfer.Re		=	(Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
	
"Tube Side Prandtl Number"
	Unity(i).Tubes.HeatTransfer.PR 		= ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;

"Tube Side Prandtl Number at Wall"
	Unity(i).Tubes.HeatTransfer.PRw 	= 	((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;

"Tube Side Film Coefficient"
	Unity(i).Tubes.HeatTransfer.htube	= 	HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
	
"Shell Side Prandtl Number"
	Unity(i).Shell.HeatTransfer.PR		=	((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;	

"Shell Side Prandtl Number at Wall"
	Unity(i).Shell.HeatTransfer.PRw		=	((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;	

"Tube Side Pressure Drop"
	Unity(i).Tubes.PressureDrop.PdTube 	= 	HE.DeltaPtubeIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,
	Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));
	
"Shell Side Phi correction for viscosity"
	Unity(i).Shell.HeatTransfer.Phi 	= 	HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);
	
"Tube Side Phi correction for viscosity"
	Unity(i).Tubes.HeatTransfer.Phi  	= 	HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
	
	else

"Pressure Drop Hot Stream"
	Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P- Unity(i).Tubes.PressureDrop.Pdtotal;	
	
"Pressure Drop Cold Stream"
	Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Shell.PressureDrop.Pdtotal;
	
"Shell Side Reynolds Number"
	Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;

"Shell Heat Transfer Coefficient"
	Unity(i).Shell.HeatTransfer.hshell	=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*
	(Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;


"Shell Pressure Baffle Window"
	Unity(i).Shell.PressureDrop.Pdwindow    = 	HE.DeltaPwindowIncremental(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.rho,Unity(i).Properties.Cold.Average.Mu,Unity(i).Baffles.Ls);

"Hot Wall Temperature"
	Unity(i).Properties.Hot.Wall.Twall  	= 	(Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;

"Cold Wall Temperature"	
	Unity(i).Properties.Cold.Wall.Twall   	= 	(Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;

"Tube Side Velocity"
	Unity(i).Tubes.HeatTransfer.Vtube  		= Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);

"Tube Side Reynolds Number"
	Unity(i).Tubes.HeatTransfer.Re			=	(Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
	
"Tube Side Prandtl Number"
	Unity(i).Tubes.HeatTransfer.PR			=	((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
	
"Tube Side Prandtl Number at Wall"
	Unity(i).Tubes.HeatTransfer.PRw			=	((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
	
"Tube Side Film Coefficient"
	Unity(i).Tubes.HeatTransfer.htube		= 	HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
	
"Shell Side Prandtl Number"
	Unity(i).Shell.HeatTransfer.PR			=	((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;

"Shell Side Prandtl Number at Wall"
	Unity(i).Shell.HeatTransfer.PRw			=	((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;

"Tube Side Pressure Drop"
	Unity(i).Tubes.PressureDrop.PdTube 		= 	HE.DeltaPtubeIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,
	Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));

"Shell Side Phi correction for viscosity"
	Unity(i).Shell.HeatTransfer.Phi 		= 	HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
	
"Tube Side Phi correction for viscosity"
	Unity(i).Tubes.HeatTransfer.Phi  		= 	HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);

end

"Tube Resistance"	
	Unity(i).Resistances.Rtube*(Unity(i).Tubes.HeatTransfer.htube*Ditube) = Dotube;
	
"Wall Resistance"
	Unity(i).Resistances.Rwall	=	Dotube*ln(Dotube/Ditube)/(2*Kwall);
	
"Shell Resistance"
	Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)	=	1;
	
"Overall Heat Transfer Coefficient Dirty"
	Unity(i).Details.Ud*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;

"Overall Heat Transfer Coefficient Clean"
	(1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);

"Exchange Surface Area"
	Unity(i).Details.A	=	Pi*Dotube*Ntt*Sumary.Lz(i);

"Baffles Spacing"
	Ltube = Unity(i).Baffles.Lsi+Unity(i).Baffles.Lso+Unity(i).Baffles.Ls*(Nb-1);
	
"Js Factor"	
	Unity(i).Shell.HeatTransfer.Js 			= 	1;
	
"Ji Factor"
	Unity(i).Shell.HeatTransfer.Ji 			= 	HE.JiFactor(Unity(i).Shell.HeatTransfer.Re);

"Jc Factor"
	Unity(i).Shell.HeatTransfer.Jc 			= 	HE.JcFactor();
	
"Jl Factor"
	Unity(i).Shell.HeatTransfer.Jl 			= 	HE.JlFactor(Unity(i).Shell.HeatTransfer.Sm);

"Jb Factor"
	Unity(i).Shell.HeatTransfer.Jb 			= 	HE.JbFactor(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Shell.HeatTransfer.Sm);

"Jr Factor"	
	Unity(i).Shell.HeatTransfer.Jr 			= 	HE.JrFactor(Unity(i).Shell.HeatTransfer.Re);
	
"Total J Factor"
	Unity(i).Shell.HeatTransfer.Jtotal		=	Unity(i).Shell.HeatTransfer.Jc*Unity(i).Shell.HeatTransfer.Jl*Unity(i).Shell.HeatTransfer.Jb*Unity(i).Shell.HeatTransfer.Jr*Unity(i).Shell.HeatTransfer.Js;
end

"Velocity Tube Side Inlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Vnozzle_in	= 0;

"Velocity Tube Side Outlet Nozzle"
	Unity(1).Tubes.PressureDrop.Vnozzle_out	= 0;
	
"Tube Pressure Drop Inlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_in	= 0;

"Tube Pressure Drop Outlet Nozzle"
	Unity(1).Tubes.PressureDrop.Pdnozzle_out	= 0;

"Velocity Shell Side Inlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Vnozzle_in	= 0;

"Velocity Shell Side Outlet Nozzle"
	Unity(1).Shell.PressureDrop.Vnozzle_out	= 0;
	
"Shell Pressure Drop Inlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Pdnozzle_in	= 0;
	
"Shell Pressure Drop Outlet Nozzle"
	Unity(1).Shell.PressureDrop.Pdnozzle_out	= 0;

"Shell Cross Flow Pressure Drop"
	Unity(1).Shell.PressureDrop.PdCross 		= 0;
	
"Shell Cross Flow Pressure Drop"
	Unity(Nb+1).Shell.PressureDrop.PdCross 		= 0;
	
"Shell Side Outlet Nozzle rho-V^2"
	Unity(1).Shell.PressureDrop.RVsquare_out = 0;
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(Nb+1).Shell.PressureDrop.RVsquare_in = 0;

SET
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
#  Set Parameters for heatex Calculation
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
	Sumary.Zones        = Nb+1;
	Pi  				= 3.14159265;
	HE.Tpass 			= Tpass;
	HE.Nss 				= Nss;
	HE.Ntt 				= Ntt;
	HE.Pattern 			= Pattern;
	HE.Bc 				= Bc;
	HE.Donozzle_Shell 	= Donozzle_Shell;
	HE.Dinozzle_Shell 	= Dinozzle_Shell;
	HE.Honozzle_Shell 	= Honozzle_Shell;
	HE.Hinozzle_Shell 	= Hinozzle_Shell;
	HE.Donozzle_Tube 	= Donozzle_Tube;
	HE.Dinozzle_Tube 	= Dinozzle_Tube;
	HE.Nb 			    = Nb;
	HE.Dishell 			= Dishell; 		
	HE.Lcf 				= Lcf;	
	HE.pitch 			= pitch;		
	HE.Dotube  			= Dotube; 	
	HE.Ditube  			= Ditube; 	
	HE.Lcd     			= Lcd; 
	HE.Ltd				= Ltd;
	side				= HE.FluidAlocation();
	
#"Tube Side Inlet Nozzle Area"
	Ainozzle_Tube = (Pi*Dinozzle_Tube*Dinozzle_Tube)/4;
	
#"Tube Side Outlet Nozzle Area"
	Aonozzle_Tube = (Pi*Donozzle_Tube*Donozzle_Tube)/4;
	
#"Tube Inlet Nozzle Pressure Loss Coeff"
	Kinlet_Tube   = 1.1;

#"Tube Outlet Nozzle Pressure Loss Coeff"
	Koutlet_Tube  = 0.7;
	
#"Shell Outlet Nozzle Area"
	Aonozzle_Shell  = (Pi*Donozzle_Shell*Donozzle_Shell)/4;
	
#"Shell Inlet Nozzle Area"
	Ainozzle_Shell  = (Pi*Dinozzle_Shell*Dinozzle_Shell)/4;
	
#"Shell Outlet Escape Area Under Nozzle"
	Aeonozzle_Shell = Pi*Donozzle_Shell*Honozzle_Shell + 0.6*Aonozzle_Shell*(1-Dotube/pitch);
	
#"Shell Inlet Escape Area Under Nozzle"
	Aeinozzle_Shell = Pi*Dinozzle_Shell*Hinozzle_Shell + 0.6*Ainozzle_Shell*(1-Dotube/pitch);
	
	
end

Model E_Shell_LMTD_Disc  
	
ATTRIBUTES
	Pallete 	= false;
	Brief 		= "Shell and Tubes Heat Exchanger with 1 shell pass - LMTD Method";
	Info 		=
	"write some information";
	
PARAMETERS

	HE 	as Plugin 	(Brief="STHE Calculations",File="heatex");
outer PP 	as Plugin	(Brief="External Physical Properties");
side 	as Integer 		(Brief="Fluid Alocation",Lower=0,Upper=1);
Pi		as constant     (Brief="Pi Number",Default=3.14159265);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Shell Geometrical Parameters
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
Tpass   		as Integer		(Brief="Number of Tube Passes",Lower=1);
Nss				as Integer		(Brief="Number of Sealing Strips pairs",Lower=1);
Dishell 		as length		(Brief="Inside Shell Diameter",Lower=10e-6);
Donozzle_Shell  as length		(Brief="Shell Outlet Nozzle Diameter",Lower=10e-6);
Dinozzle_Shell  as length		(Brief="Shell Inlet Nozzle Diameter",Lower=10e-6);
Aonozzle_Shell  as area			(Brief="Shell Outlet Nozzle Area",Lower=10e-6);
Ainozzle_Shell  as area			(Brief="Shell Inlet Nozzle Area",Lower=10e-6);
Aeonozzle_Shell as area			(Brief="Shell Outlet Escape Area Under Nozzle",Lower=10e-6);
Aeinozzle_Shell as area			(Brief="Shell Inlet Escape Area Under Nozzle",Lower=10e-6);
Hinozzle_Shell  as length		(Brief="Height Under Shell Inlet Nozzle",Lower=10e-6);
Honozzle_Shell  as length		(Brief="Height Under Shell Outlet Nozzle",Lower=10e-6);
Lcf 			as length		(Brief="Bundle-to-Shell Clearance",Lower=10e-8);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
#				Tubes Geometrical Parameters						#
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
Ntt			  as Integer		(Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
Pattern       as Integer		(Brief="Tube Layout Characteristic Angle",Lower=30);
Ltube		  as length			(Brief="Effective Tube Length",Lower=0.1);
pitch 		  as length			(Brief="Tube Pitch",Lower=1e-8);
Kwall		  as conductivity 	(Brief="Tube Wall Material Thermal Conductivity");
Dotube  	  as length			(Brief="Tube Outside Diameter",Lower=10e-6);
Ditube 		  as length			(Brief="Tube Inside Diameter",Lower=10e-6);
Donozzle_Tube as length			(Brief="Tube Outlet Nozzle Diameter",Lower=10e-6);
Dinozzle_Tube as length			(Brief="Tube Inlet Nozzle Diameter",Lower=10e-6);
Aonozzle_Tube as area			(Brief="Tube Outlet Nozzle Area",Lower=10e-6);
Ainozzle_Tube as area			(Brief="Tube Inlet Nozzle Area",Lower=10e-6);
Kinlet_Tube   as positive		(Brief="Tube Inlet Nozzle Pressure Loss Coeff",Default=1.1);
Koutlet_Tube  as positive		(Brief="Tube Outlet Nozzle Pressure Loss Coeff",Default=0.7);

#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
# 	Baffles Geometrical Parameters
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
Bc      		as Integer		(Brief="Baffle Cut",Default=25,Lower=25);
Nb          	as Integer		(Brief="Number of Baffles",Default=4);
Lcd     		as length		(Brief="Baffle-to-Shell Clearance",Lower=10e-8);
Ltd				as length		(Brief="Tube-to-Bafflehole Clearance",Lower=10e-8);
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
VARIABLES

Unity(Nb+1) 	as Heatex_Discretized_LMTD; 
Sumary          as Profiles;

CONNECTIONS

Unity([1:Nb]).Outlet.Hot  	to Unity([2:Nb+1]).Inlet.Hot;
Unity([2:Nb+1]).Outlet.Cold to Unity([1:Nb]).Inlet.Cold;

EQUATIONS

"Hot Temperatures"
	Sumary.Thot 	= Unity.Outlet.Hot.T;
	
"Cold Temperatures"
	Sumary.Tcold 	= Unity.Outlet.Cold.T ;

"Hot Pressures"
	Sumary.Phot 	= Unity.Outlet.Hot.P ;
	
"Cold Pressures"
	Sumary.Pcold 	= Unity.Outlet.Cold.P ;

"Average Shell Side Film Coefficient"
	Sumary.hshellaverage       	= sum(Unity.Shell.HeatTransfer.hshell)/Sumary.Zones;

"Average Tube Side Film Coefficient"
	Sumary.htubeaverage        	= sum(Unity.Tubes.HeatTransfer.htube)/Sumary.Zones;

"Average Overall Heat Transfer Coefficient Dirty"
	Sumary.Udaverage       		= sum(Unity.Details.Ud)/Sumary.Zones;
	
"Average Overall Heat Transfer Coefficient Clean"
	Sumary.Ucaverage       		= sum(Unity.Details.Uc)/Sumary.Zones;	

"Area Total"
	Sumary.Area 				= sum(Unity.Details.A);

"Duty Total"
	Sumary.Q    				= sum(Unity.Details.Q);

"Length Inlet zone"
	Sumary.Lz(1)    			= Unity(1).Baffles.Lsi;

"Length Outlet zone"
	Sumary.Lz(Nb+1) 			= Unity(1).Baffles.Lso;

"Total Shell Side Pressure Drop"
	Sumary.PdropShell 			= sum(Unity.Shell.PressureDrop.Pdtotal);

"Total Tube Side Pressure Drop"
	Sumary.PdropTubes 			= sum(Unity.Tubes.PressureDrop.Pdtotal);

"Total Tube Side Nozzles Pressure Drop"
	Sumary.PdropTubesNozzle		= sum(Unity.Tubes.PressureDrop.Pdnozzle_in + Unity.Tubes.PressureDrop.Pdnozzle_out);

"Total Shell Side Nozzles Pressure Drop"
	Sumary.PdropShellNozzle		= sum(Unity.Shell.PressureDrop.Pdnozzle_in + Unity.Shell.PressureDrop.Pdnozzle_out);

"Total Shell Side Window Pressure Drop"
	Sumary.PdropWin				= sum(Unity.Shell.PressureDrop.Pdwindow);

"Total Shell Side Cross Flow Pressure Drop"
	Sumary.PdropCross 			= sum(Unity.Shell.PressureDrop.PdCross);

"Total Shell Side Ends Pressure Drop"
	Sumary.PdropEnds			= sum(Unity.Shell.PressureDrop.PdEndZones);

"Shell Side Cross Flow Area" 
	Unity(1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(1).Baffles.Lsi);
	
"Shell Side Cross Flow Area" 
	Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso);

if side equal 1
	
	then
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(1).Shell.PressureDrop.RVsquare_in = Unity(1).Properties.Hot.Inlet.rho*(Unity(1).Shell.PressureDrop.Vnozzle_in)^2;

"Shell Side Outlet Nozzle rho-V^2"
	Unity(Nb+1).Shell.PressureDrop.RVsquare_out = Unity(Nb+1).Properties.Hot.Outlet.rho*(Unity(Nb+1).Shell.PressureDrop.Vnozzle_out)^2;

"Shell Pressure End Zones"
	Unity(1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(1).Shell.HeatTransfer.Re,Unity(1).Baffles.Ls,Unity(1).Baffles.Lso,
	Unity(1).Baffles.Lsi,Unity(1).Properties.Hot.Inlet.Fw,Unity(1).Shell.HeatTransfer.Phi,Unity(1).Properties.Hot.Average.rho);

"Shell Pressure End Zones"
	Unity(Nb+1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(Nb+1).Shell.HeatTransfer.Re,Unity(Nb+1).Baffles.Ls,Unity(Nb+1).Baffles.Lso,
	Unity(Nb+1).Baffles.Lsi,Unity(Nb+1).Properties.Hot.Inlet.Fw,Unity(Nb+1).Shell.HeatTransfer.Phi,Unity(Nb+1).Properties.Hot.Average.rho);

"Pressure Drop Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Pdnozzle_in	= 0.5*Kinlet_Tube*Unity(1).Properties.Cold.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;

"Velocity Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Vnozzle_in	= Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Tube);

"Pressure Drop Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out	= 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Cold.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;

"Velocity Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out	= Unity(Nb+1).Properties.Cold.Inlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Aonozzle_Tube);

"Shell Pressure Drop Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Pdnozzle_in	= (0.5*Unity(1).Properties.Hot.Inlet.Fw^2/Unity(1).Properties.Hot.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Vnozzle_in	= Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Shell);

"Shell Pressure Drop Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out	= (0.5*Unity(Nb+1).Properties.Hot.Outlet.Fw^2/Unity(Nb+1).Properties.Hot.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Vnozzle_out	= Unity(Nb+1).Properties.Hot.Outlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Shell);


	else
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(1).Shell.PressureDrop.RVsquare_in = Unity(1).Properties.Cold.Inlet.rho*(Unity(1).Shell.PressureDrop.Vnozzle_in)^2;

"Shell Side Outlet Nozzle rho-V^2"
	Unity(Nb+1).Shell.PressureDrop.RVsquare_out = Unity(Nb+1).Properties.Cold.Outlet.rho*(Unity(Nb+1).Shell.PressureDrop.Vnozzle_out)^2;

"Shell Pressure End Zones"
	Unity(1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(1).Shell.HeatTransfer.Re,Unity(1).Baffles.Ls,Unity(1).Baffles.Lso,
	Unity(1).Baffles.Lsi,Unity(1).Properties.Cold.Inlet.Fw,Unity(1).Shell.HeatTransfer.Phi,Unity(1).Properties.Cold.Average.rho);

"Shell Pressure End Zones"
	Unity(Nb+1).Shell.PressureDrop.PdEndZones  	= 	HE.DeltaPendZonesIncremental(Unity(Nb+1).Shell.HeatTransfer.Re,Unity(Nb+1).Baffles.Ls,Unity(Nb+1).Baffles.Lso,
	Unity(Nb+1).Baffles.Lsi,Unity(Nb+1).Properties.Cold.Inlet.Fw,Unity(Nb+1).Shell.HeatTransfer.Phi,Unity(Nb+1).Properties.Cold.Average.rho);

"Pressure Drop Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Pdnozzle_in	= 0.5*Kinlet_Tube*Unity(1).Properties.Hot.Inlet.rho*Unity(1).Tubes.PressureDrop.Vnozzle_in^2;
	
"Velocity Tube Side Inlet Nozzle"
	Unity(1).Tubes.PressureDrop.Vnozzle_in 	= Unity(1).Properties.Hot.Inlet.Fw/(Unity(1).Properties.Hot.Inlet.rho*Ainozzle_Tube);
	
"Pressure Drop Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_out	= 0.5*Koutlet_Tube*Unity(Nb+1).Properties.Hot.Outlet.rho*Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out^2;
	
"Velocity Tube Side Outlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Vnozzle_out	= Unity(Nb+1).Properties.Hot.Inlet.Fw/(Unity(Nb+1).Properties.Hot.Outlet.rho*Aonozzle_Tube);	

"Shell Pressure Drop Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Pdnozzle_in	= (0.5*Unity(1).Properties.Cold.Inlet.Fw^2/Unity(1).Properties.Cold.Inlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Inlet Nozzle"
	Unity(1).Shell.PressureDrop.Vnozzle_in	= Unity(1).Properties.Cold.Inlet.Fw/(Unity(1).Properties.Cold.Inlet.rho*Ainozzle_Shell);

"Shell Pressure Drop Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Pdnozzle_out = (0.5*Unity(Nb+1).Properties.Cold.Outlet.Fw^2/Unity(Nb+1).Properties.Cold.Outlet.rho)*((1/Ainozzle_Shell^2)+(1/Aeinozzle_Shell^2));

"Velocity Shell Side Outlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Vnozzle_out  = Unity(Nb+1).Properties.Cold.Outlet.Fw/(Unity(Nb+1).Properties.Cold.Outlet.rho*Ainozzle_Shell);


end


for i in [2:Nb]

"Pressure Drop Tube Side Inlet Nozzle"
	Unity(i).Tubes.PressureDrop.Pdnozzle_in	= 0;

"Velocity Tube Side Inlet Nozzle"
	Unity(i).Tubes.PressureDrop.Vnozzle_in 	= 0;

"Pressure Drop Tube Side Outlet Nozzle"
	Unity(i).Tubes.PressureDrop.Pdnozzle_out	= 0;	

"Velocity Tube Side Outlet Nozzle"
	Unity(i).Tubes.PressureDrop.Vnozzle_out	= 0;	

"Shell Pressure Drop Inlet Nozzle"
	Unity(i).Shell.PressureDrop.Pdnozzle_in	= 0;

"Velocity Shell Side Inlet Nozzle"
	Unity(i).Shell.PressureDrop.Vnozzle_in	= 0;

"Shell Pressure Drop Outlet Nozzle"
	Unity(i).Shell.PressureDrop.Pdnozzle_out = 0;

"Velocity Shell Side Outlet Nozzle"
	Unity(i).Shell.PressureDrop.Vnozzle_out  = 0;

"Shell Pressure End Zones"
	Unity(i).Shell.PressureDrop.PdEndZones  	= 	0;

"Shell Side Outlet Nozzle rho-V^2"
	Unity(i).Shell.PressureDrop.RVsquare_out = 0;
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(i).Shell.PressureDrop.RVsquare_in = 0;
	
if side equal 1
	
	then
	
"Shell Pressure Drop Cross Flow"
	Unity(i).Shell.PressureDrop.PdCross 		= 	HE.DeltaPcrossIncremental(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Hot.Average.rho);


	else

"Shell Pressure Drop Cross Flow"
	Unity(i).Shell.PressureDrop.PdCross     = 	HE.DeltaPcrossIncremental(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Baffles.Lso,Unity(i).Baffles.Lsi,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Phi,Unity(i).Properties.Cold.Average.rho);


end

"Zone Length"
	Sumary.Lz(i) = Unity(1).Baffles.Ls;

"Shell Side Cross Flow Area" 
	Unity(i).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(i).Baffles.Ls);
	
end


for i in [1:Nb+1]


if side equal 1	
	
	then 
"Pressure Drop Hot Stream"
	Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P - Unity(i).Shell.PressureDrop.Pdtotal;	

"Pressure Drop Cold Stream"
	Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Tubes.PressureDrop.Pdtotal;
	
"Shell Side Reynolds Number"
	Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;

"Shell Heat Transfer Coefficient"
	Unity(i).Shell.HeatTransfer.hshell =Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;

"Shell Pressure Baffle Window"
	Unity(i).Shell.PressureDrop.Pdwindow		= 	HE.DeltaPwindowIncremental(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.rho,Unity(i).Properties.Hot.Average.Mu,Unity(i).Baffles.Ls);

"Hot Wall Temperature"
	Unity(i).Properties.Hot.Wall.Twall  = (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
	
"Cold Wall Temperature"
	Unity(i).Properties.Cold.Wall.Twall = 	(Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;

"Tube Side Velocity"
	Unity(i).Tubes.HeatTransfer.Vtube  	= 	Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);

"Tube Side Reynolds Number"
	Unity(i).Tubes.HeatTransfer.Re		=	(Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
	
"Tube Side Prandtl Number"
	Unity(i).Tubes.HeatTransfer.PR 		= 	((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;

"Tube Side Prandtl Number at Wall"
	Unity(i).Tubes.HeatTransfer.PRw 	= 	((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;

"Tube Side Film Coefficient"
	Unity(i).Tubes.HeatTransfer.htube	= 	HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
	
"Shell Side Prandtl Number"
	Unity(i).Shell.HeatTransfer.PR		=	((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;	

"Shell Side Prandtl Number at Wall"
	Unity(i).Shell.HeatTransfer.PRw		=	((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;	

"Tube Side Pressure Drop"
	Unity(i).Tubes.PressureDrop.PdTube 	= 	HE.DeltaPtubeIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,
	Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));
	
"Shell Side Phi correction for viscosity"
	Unity(i).Shell.HeatTransfer.Phi 	= 	HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);
	
"Tube Side Phi correction for viscosity"
	Unity(i).Tubes.HeatTransfer.Phi  	= 	HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
	
	else

"Pressure Drop Hot Stream"
	Unity(i).Outlet.Hot.P  = Unity(i).Inlet.Hot.P- Unity(i).Tubes.PressureDrop.Pdtotal;	
	
"Pressure Drop Cold Stream"
	Unity(i).Outlet.Cold.P  = Unity(i).Inlet.Cold.P - Unity(i).Shell.PressureDrop.Pdtotal;
	
"Shell Side Reynolds Number"
	Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;

"Shell Heat Transfer Coefficient"
	Unity(i).Shell.HeatTransfer.hshell	=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*
	(Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;


"Shell Pressure Baffle Window"
	Unity(i).Shell.PressureDrop.Pdwindow    = 	HE.DeltaPwindowIncremental(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.rho,Unity(i).Properties.Cold.Average.Mu,Unity(i).Baffles.Ls);


"Hot Wall Temperature"
	Unity(i).Properties.Hot.Wall.Twall  	= 	(Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
	
"Cold Wall Temperature"
	Unity(i).Properties.Cold.Wall.Twall   	= 	(Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;

"Tube Side Velocity"
	Unity(i).Tubes.HeatTransfer.Vtube  		= Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);

"Tube Side Reynolds Number"
	Unity(i).Tubes.HeatTransfer.Re			=	(Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
	
"Tube Side Prandtl Number"
	Unity(i).Tubes.HeatTransfer.PR			=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
	
"Tube Side Prandtl Number at Wall"
	Unity(i).Tubes.HeatTransfer.PRw			=	((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
	
	
"Tube Side Film Coefficient"
	Unity(i).Tubes.HeatTransfer.htube		= 	HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
	
"Shell Side Prandtl Number"
	Unity(i).Shell.HeatTransfer.PR			=	((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;

"Shell Side Prandtl Number at Wall"
	Unity(i).Shell.HeatTransfer.PRw			=	((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;

"Tube Side Pressure Drop"
	Unity(i).Tubes.PressureDrop.PdTube 		= 	HE.DeltaPtubeIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,
	Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));

"Shell Side Phi correction for viscosity"
	Unity(i).Shell.HeatTransfer.Phi 		= 	HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
	
"Tube Side Phi correction for viscosity"
	Unity(i).Tubes.HeatTransfer.Phi  		= 	HE.PhiCorrection(Unity(i).Properties.Hot.Average.Mu,Unity(i).Properties.Hot.Wall.Mu);

end

"Tube Resistance"	
	Unity(i).Resistances.Rtube*(Unity(i).Tubes.HeatTransfer.htube*Ditube) = Dotube;
	
"Wall Resistance"
	Unity(i).Resistances.Rwall	=	Dotube*ln(Dotube/Ditube)/(2*Kwall);
	
"Shell Resistance"
	Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)	=	1;
	
"Overall Heat Transfer Coefficient Dirty"
	Unity(i).Details.Ud*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;

"Overall Heat Transfer Coefficient Clean"
	(1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);

"Exchange Surface Area"
	Unity(i).Details.A	=	Pi*Dotube*Ntt*Sumary.Lz(i);

"Baffles Spacing"
	Ltube = Unity(i).Baffles.Lsi+Unity(i).Baffles.Lso+Unity(i).Baffles.Ls*(Nb-1);

"Js Factor"	
	Unity(i).Shell.HeatTransfer.Js 			= 	1;

"Ji Factor"
	Unity(i).Shell.HeatTransfer.Ji 			= 	HE.JiFactor(Unity(i).Shell.HeatTransfer.Re);

"Jc Factor"
	Unity(i).Shell.HeatTransfer.Jc 			= 	HE.JcFactor();
	
"Jl Factor"
	Unity(i).Shell.HeatTransfer.Jl 			= 	HE.JlFactor(Unity(i).Shell.HeatTransfer.Sm);

"Jb Factor"
	Unity(i).Shell.HeatTransfer.Jb 			= 	HE.JbFactor(Unity(i).Shell.HeatTransfer.Re,Unity(i).Baffles.Ls,Unity(i).Shell.HeatTransfer.Sm);

"Jr Factor"	
	Unity(i).Shell.HeatTransfer.Jr 			= 	HE.JrFactor(Unity(i).Shell.HeatTransfer.Re);
	
"Total J Factor"
	Unity(i).Shell.HeatTransfer.Jtotal		=	Unity(i).Shell.HeatTransfer.Jc*Unity(i).Shell.HeatTransfer.Jl*Unity(i).Shell.HeatTransfer.Jb*Unity(i).Shell.HeatTransfer.Jr*Unity(i).Shell.HeatTransfer.Js;
end

"Velocity Tube Side Inlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Vnozzle_in	= 0;

"Velocity Tube Side Outlet Nozzle"
	Unity(1).Tubes.PressureDrop.Vnozzle_out	= 0;
	
"Tube Pressure Drop Inlet Nozzle"
	Unity(Nb+1).Tubes.PressureDrop.Pdnozzle_in	= 0;

"Tube Pressure Drop Outlet Nozzle"
	Unity(1).Tubes.PressureDrop.Pdnozzle_out	= 0;

"Velocity Shell Side Inlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Vnozzle_in	= 0;

"Velocity Shell Side Outlet Nozzle"
	Unity(1).Shell.PressureDrop.Vnozzle_out	= 0;
	
"Shell Pressure Drop Inlet Nozzle"
	Unity(Nb+1).Shell.PressureDrop.Pdnozzle_in	= 0;
	
"Shell Pressure Drop Outlet Nozzle"
	Unity(1).Shell.PressureDrop.Pdnozzle_out	= 0;

"Shell Cross Flow Pressure Drop"
	Unity(1).Shell.PressureDrop.PdCross 		= 0;
	
"Shell Cross Flow Pressure Drop"
	Unity(Nb+1).Shell.PressureDrop.PdCross 		= 0;
	
"Shell Side Outlet Nozzle rho-V^2"
	Unity(1).Shell.PressureDrop.RVsquare_out = 0;
	
"Shell Side inlet Nozzle rho-V^2"
	Unity(Nb+1).Shell.PressureDrop.RVsquare_in = 0;

SET
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
#  Set Parameters for heatex Calculation
#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
	Sumary.Zones        = Nb+1;
	Pi  				= 3.14159265;
	HE.Tpass 			= Tpass;
	HE.Nss 				= Nss;
	HE.Ntt 				= Ntt;
	HE.Pattern 			= Pattern;
	HE.Bc 				= Bc;
	HE.Donozzle_Shell 	= Donozzle_Shell;
	HE.Dinozzle_Shell 	= Dinozzle_Shell;
	HE.Honozzle_Shell 	= Honozzle_Shell;
	HE.Hinozzle_Shell 	= Hinozzle_Shell;
	HE.Donozzle_Tube 	= Donozzle_Tube;
	HE.Dinozzle_Tube 	= Dinozzle_Tube;
	HE.Nb 			    = Nb;
	HE.Dishell 			= Dishell; 		
	HE.Lcf 				= Lcf;	
	HE.pitch 			= pitch;		
	HE.Dotube  			= Dotube; 	
	HE.Ditube  			= Ditube; 	
	HE.Lcd     			= Lcd; 
	HE.Ltd				= Ltd;
	side				= HE.FluidAlocation();
	
#"Tube Side Inlet Nozzle Area"
	Ainozzle_Tube = (Pi*Dinozzle_Tube*Dinozzle_Tube)/4;
	
#"Tube Side Outlet Nozzle Area"
	Aonozzle_Tube = (Pi*Donozzle_Tube*Donozzle_Tube)/4;
	
#"Tube Inlet Nozzle Pressure Loss Coeff"
	Kinlet_Tube   = 1.1;

#"Tube Outlet Nozzle Pressure Loss Coeff"
	Koutlet_Tube  = 0.7;
	
#"Shell Outlet Nozzle Area"
	Aonozzle_Shell  = (Pi*Donozzle_Shell*Donozzle_Shell)/4;
	
#"Shell Inlet Nozzle Area"
	Ainozzle_Shell  = (Pi*Dinozzle_Shell*Dinozzle_Shell)/4;
	
#"Shell Outlet Escape Area Under Nozzle"
	Aeonozzle_Shell = Pi*Donozzle_Shell*Honozzle_Shell + 0.6*Aonozzle_Shell*(1-Dotube/pitch);
	
#"Shell Inlet Escape Area Under Nozzle"
	Aeinozzle_Shell = Pi*Dinozzle_Shell*Hinozzle_Shell + 0.6*Ainozzle_Shell*(1-Dotube/pitch);
	
end