#*-------------------------------------------------------------------
* 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 column cost model
*--------------------------------------------------------------------
* 
* This sample file needs VRTherm DEMO(www.vrtech.com.br) to run
* SectionColumn_Test and needs VRTherm full to run the distillation
* column flowsheet.
*
*----------------------------------------------------------------------
* Author: Núbia do Carmo Ferreira
* $Id: sample_column.mso 247 2007-04-24 13:44:18Z rafael $
*--------------------------------------------------------------------*#

using "costs/column_cost";

FlowSheet sample_Distillation_kettle_cond_cost
	
	PARAMETERS
	PP	as Plugin(Brief="Physical Properties",
		Type="PP",
		Components = [ "isobutane", "n-pentane", "propylene",
		"benzene", "isobutene" ],
		LiquidModel = "PR",
		VapourModel = "PR"
	);
	NComp	as Integer;
	
	
	SET
	NComp = PP.NumberOfComponents;

	DEVICES
	col as Distillation_kettle_cond_cost;
	feed as source;
	zero as stream;
	Qc as energy_source;
	Qr as energy_source;
	
	CONNECTIONS
	feed.Outlet to col.trays(5).Inlet;
	zero 		to col.reb.Inlet;
	zero 		to col.trays([1:4]).Inlet;
	zero 		to col.trays([6:col.NTrays]).Inlet;
	Qc.OutletQ 	to col.cond.InletQ;
	Qr.OutletQ	to col.reb.InletQ;
	
	SPECIFY
	feed.F = 113.4 * 'kmol/h';
	feed.T = 291 * 'K';
	feed.P = 168.3 * 'kPa';
	feed.Composition = 1/NComp;
	
	zero.F = 0 * 'kmol/h';
	zero.T = 300 * 'K';
	zero.P = 1 * 'atm';
	zero.z = 1/NComp;
	zero.v = 0;
	zero.h = 0 * 'J/mol';
	
	col.sptop.Outlet2.F = 85 * 'kmol/h';
	col.reb.OutletL.F = 28.4 * 'kmol/h';
	col.sptop.frac = 0.444445;
	col.cond.OutletV.F = 0 * 'kmol/h';
	Qr.OutletQ = 3.7743e6 * 'kJ/h';
	Qc.OutletQ = -3.71e6 * 'kJ/h';
	col.pump1.dP = 16 * 'kPa';
	col.trays.Emv = 1;
	col.alfaTopo = 2;
	
	SET
	col.NTrays = 8;
	col.cond.V = 2 * 'm^3';
	col.cond.Across = 1 * 'm^2';
	col.trays.V = 4 * 'ft^3';
	col.trays.Ah = 0.394 * 'ft^2';
	col.trays.lw = 20.94 * 'in';
	col.trays.hw = 0.125 * 'ft';
	col.trays.Q = 0 * 'kW';
	col.trays.beta = 0.6;
	col.trays.alfa = 4;
	col.trays.Ap = 3.94 * 'ft^2';
	col.reb.V = 2 * 'm^3';
	col.reb.Across = 1 * 'm^2';
	
	#cost
	col.Material = "Stainless steel 316";
	col.Tray_Material = "Stainless steel 304";
	col.Tray_Type = "Valve";
	col.Di = 0.9 * 'm';
	col.Lt = 17.5 * 'm';
	col.Tb = 0.01250 * 'm';
	col.Tp = 0.01349 * 'm';
	col.D  = 2.4 * 'm';
	col.dens_mass_material = 8000 * 'kg/m^3';

	col.Cost(1,:) = [6.950,0.1808,0.02468,0.01580];
	col.Cost(2,:) = [834.86,0,0,0];
	col.Cost(3,:) = [6.448,0.21887,0.02297,0];
	col.Cost(4,:) = [1017.0,0,0,0];
	col.Cost(5,:) = [1.7,0,0,0]; #For "Stainless steel 304" - tray material
	col.Cost(6,:) = [278.38,0.1739,0,0];
	col.Cost(7,:) = [1.189,0.1894,0,0];
	col.Cost(8,:) = [1.401,0.2376,0,0];
	col.Cost(9,:) = [1.525,0.2585,0,0];
	col.Cost(10,:) = [2.306,0.3674,0,0];
	col.Cost(11,:) = [2.25,1.0414,0,0];
	col.Cost(12,:) = [1,0,0,0]; #For valve trays
	
	INITIAL
	# condenser
	col.cond.OutletL.T = 260 *'K';
	col.cond.Level = 1 * 'm';
	col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];

	# reboiler
	col.reb.OutletL.T = 300 *'K';
	col.reb.Level = 1 * 'm';
	col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];

	# column trays
	col.trays.OutletL.T = [290:(300-290)/(col.NTrays-1):300] * 'K';
	col.trays.Level = 1.2 * col.trays.hw;
	col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];

	OPTIONS
	TimeStep = 0.1;
	TimeEnd = 50;
	#Dynamic = false;	
end