Context Navigation

HeatExchangerSimplified_cost.mso @ 696

Last change on this file since 696 was 639, checked in by Argimiro Resende Secchi, 15 years ago
Fix cost model due to Heatex changes.
File size: 7.9 KB

Rev	Line
[491]	1	#*-------------------------------------------------------------------
	2	* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
	3	*
	4	* This LIBRARY is free software; you can distribute it and/or modify
	5	* it under the therms of the ALSOC FREE LICENSE as available at
	6	* http://www.enq.ufrgs.br/alsoc.
	7	*
	8	* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
	9	* from http://www.rps.eng.br Copyright (C) 2002-2004.
	10	* All rights reserved.
	11	*
	12	* EMSO is distributed under the therms of the ALSOC LICENSE as
	13	* available at http://www.enq.ufrgs.br/alsoc.
	14	*
	15	*--------------------------------------------------------------------
	16	* Author: Núbia do Carmo Ferreira
	17	* $Id: HeatExchangerDetailed.mso 197 2007-03-08 14:31:57Z bicca $
	18	------------------------------------------------------------------#
	19
	20	using "heat_exchangers/Heatex";
	21
	22	# This model is valid if 14m^2 < A < 1100 m^2
	23
	24	Model HeatExchanger_LMTD_cost as Heatex_LMTD
	25	ATTRIBUTES
	26	Pallete = true;
	27	Icon = "icon/HeatExchanger_LMTD";
	28	Brief = "Heat Exchanger Block - LMTD Method";
	29	Info =
	30	"to be documented.";
	31
	32	PARAMETERS
	33	Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Stainless steel 347", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium", "Hastelloy"],
	34	Default = "Stainless steel 316");
	35	Cost(8,3) as Real;
	36
	37	VARIABLES
	38	Ce as currency (Brief="Capital Cost");
	39	Cb as currency (Brief="Basic Cost");
	40	Fd as positive (Brief="Cost Factor based on the type of the heat exchanger");
	41	Fp as positive (Brief="Cost Factor based on the project pressure");
	42	Fm as positive (Brief="Cost Factor based on the construction material");
	43	Pmax as pressure (Brief="Average Pressure");
	44
	45	EQUATIONS
	46	"Average pressure"
[639]	47	Pmax = max( [OutletHot.P , OutletCold.P] );
[491]	48
	49	"Capital Cost"
	50	Ce = CbFdFp*Fm;
	51
	52	"Basic Cost"
[639]	53	Cb = 'US$'exp(Cost(1,1) + Cost(1,2)ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
[491]	54
	55	"Cost Factor based on the type of the heat exchanger"
[639]	56	Fd = exp(Cost(2,1) + Cost(2,2)ln(A/'m^2') + Cost(2,3)ln(A/'m^2'));
[491]	57
	58	if Pmax <= 700 * 'kPa' then
	59	"Cost Factor based on the project pressure" # verificar
[639]	60	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	61	else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
[639]	62	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	63	else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
[639]	64	Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
[491]	65	else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
[639]	66	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
[491]	67	else
[639]	68	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
[491]	69	end
	70	end
	71	end
	72	end
	73
	74	"Cost Factor based on the construction material"
[639]	75	Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
[491]	76
	77	end
	78
	79	Model Shell_and_tubes_LMTD_cost as Heatex_LMTD
	80	ATTRIBUTES
	81	Pallete = true;
	82	Icon = "icon/Shell_and_Tubes_LMTD";
	83	Brief = "Shell and Tubes Heat Exchanger with 1 or 2 shell pass - LMTD Method";
	84	Info =
	85	"to be documented.";
	86
	87	PARAMETERS
	88	Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Stainless steel 347", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium", "Hastelloy"],
	89	Default = "Stainless steel 316");
	90	Cost(8,3) as Real;
	91
	92	VARIABLES
	93	Ce as currency (Brief="Capital Cost");
	94	Cb as currency (Brief="Basic Cost");
	95	Fd as positive (Brief="Cost Factor based on the type of the heat exchanger");
	96	Fp as positive (Brief="Cost Factor based on the project pressure");
	97	Fm as positive (Brief="Cost Factor based on the construction material");
	98	Pmax as pressure (Brief="Average Pressure");
	99
	100	SET
	101
	102	ExchangerType = "Shell and Tube";
	103
	104	EQUATIONS
	105	"Average pressure"
[639]	106	Pmax = max( [OutletHot.P , OutletCold.P] );
[491]	107
	108	"Capital Cost"
	109	Ce = CbFdFp*Fm;
	110
	111	"Basic Cost"
[639]	112	Cb = 'US$'exp(Cost(1,1) + Cost(1,2)ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
[491]	113
	114	"Cost Factor based on the type of the heat exchanger"
[639]	115	Fd = exp(Cost(2,1) + Cost(2,2)ln(A/'m^2') + Cost(2,3)ln(A/'m^2'));
[491]	116
	117	if Pmax <= 700 * 'kPa' then
	118	"Cost Factor based on the project pressure" # verificar
[639]	119	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	120	else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
[639]	121	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	122	else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
[639]	123	Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
[491]	124	else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
[639]	125	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
[491]	126	else
[639]	127	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
[491]	128	end
	129	end
	130	end
	131	end
	132
	133	"Cost Factor based on the construction material"
[639]	134	Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
[491]	135
	136	end
	137
	138	Model HeatExchanger_NTU_cost as Heatex_NTU
	139	ATTRIBUTES
	140	Pallete = true;
	141	Icon = "icon/HeatExchanger_NTU";
	142	Brief = "Heat Exchanger Block - NTU Method";
	143	Info =
	144	"to be documented.";
	145
	146	PARAMETERS
	147	Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Stainless steel 347", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium", "Hastelloy"],
	148	Default = "Stainless steel 316");
	149	Cost(8,3) as Real;
	150
	151	VARIABLES
	152	Ce as currency (Brief="Capital Cost");
	153	Cb as currency (Brief="Basic Cost");
	154	Fd as positive (Brief="Cost Factor based on the type of the heat exchanger");
	155	Fp as positive (Brief="Cost Factor based on the project pressure");
	156	Fm as positive (Brief="Cost Factor based on the construction material");
	157	Pmax as pressure (Brief="Average Pressure");
	158
	159	EQUATIONS
	160	"Average pressure"
[639]	161	Pmax = max( [OutletHot.P , OutletCold.P] );
[491]	162
	163	"Capital Cost"
	164	Ce = CbFdFp*Fm;
	165
	166	"Basic Cost"
[639]	167	Cb = 'US$'exp(Cost(1,1) + Cost(1,2)ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
[491]	168
	169	"Cost Factor based on the type of the heat exchanger"
[639]	170	Fd = exp(Cost(2,1) + Cost(2,2)ln(A/'m^2') + Cost(2,3)ln(A/'m^2'));
[491]	171
	172	if Pmax <= 700 * 'kPa' then
	173	"Cost Factor based on the project pressure" # verificar
[639]	174	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	175	else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
[639]	176	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	177	else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
[639]	178	Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
[491]	179	else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
[639]	180	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
[491]	181	else
[639]	182	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
[491]	183	end
	184	end
	185	end
	186	end
	187
	188	"Cost Factor based on the construction material"
[639]	189	Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
[491]	190
	191	end
	192
	193	Model Shell_and_tubes_NTU_cost as Heatex_NTU
	194	ATTRIBUTES
	195	Pallete = true;
	196	Icon = "icon/Shell_and_Tubes_NTU";
	197	Brief = "Shell and Tubes Heat Exchanger with 1 or 2 shell pass - NTU Method";
	198	Info =
	199	"to be documented.";
	200
	201	PARAMETERS
	202	Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Stainless steel 347", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium", "Hastelloy"],
	203	Default = "Stainless steel 316");
	204	Cost(8,3) as Real;
	205
	206	VARIABLES
	207	Ce as currency (Brief="Capital Cost");
	208	Cb as currency (Brief="Basic Cost");
	209	Fd as positive (Brief="Cost Factor based on the type of the heat exchanger");
	210	Fp as positive (Brief="Cost Factor based on the project pressure");
	211	Fm as positive (Brief="Cost Factor based on the construction material");
	212	Pmax as pressure (Brief="Average Pressure");
	213
	214	SET
	215
	216	ExchangerType = "Shell and Tube";
	217
	218	EQUATIONS
	219	"Average pressure"
[639]	220	Pmax = max( [OutletHot.P , OutletCold.P] );
[491]	221
	222	"Capital Cost"
	223	Ce = CbFdFp*Fm;
	224
	225	"Basic Cost"
[639]	226	Cb = 'US$'exp(Cost(1,1) + Cost(1,2)ln(A/'m^2') + Cost(1,3)*(ln(A/'m^2'))^2);
[491]	227
	228	"Cost Factor based on the type of the heat exchanger"
[639]	229	Fd = exp(Cost(2,1) + Cost(2,2)ln(A/'m^2') + Cost(2,3)ln(A/'m^2'));
[491]	230
	231	if Pmax <= 700 * 'kPa' then
	232	"Cost Factor based on the project pressure" # verificar
[639]	233	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	234	else if 700 * 'kPa' < Pmax and Pmax < 2100 * 'kPa' then
[639]	235	Fp = Cost(5,1) + Cost(5,2)*ln(A/'m^2');
[491]	236	else if 2100 * 'kPa' < Pmax and Pmax < 4200 * 'kPa' then
[639]	237	Fp = Cost(6,1) + Cost(6,2)*ln(A/'m^2');
[491]	238	else if 4200 * 'kPa' < Pmax and Pmax < 6200 * 'kPa' then
[639]	239	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2');
[491]	240	else
[639]	241	Fp = Cost(7,1) + Cost(7,2)*ln(A/'m^2'); # verificar
[491]	242	end
	243	end
	244	end
	245	end
	246
	247	"Cost Factor based on the construction material"
[639]	248	Fm = Cost(8,1) + Cost(8,2)*ln(A/'m^2');
[491]	249
	250	end
	251

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format