Context Navigation

HEX_Engine.mso @ 659

Last change on this file since 659 was 529, checked in by gerson bicca, 15 years ago
updated Heatex simplified model
Property svn:eol-style set to `native` Property svn:keywords set to `Id`
File size: 14.4 KB

Line
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	* Author: Gerson Balbueno Bicca
16	* $Id: HEX_Engine.mso 529 2008-05-30 01:21:25Z bicca $
17	--------------------------------------------------------------------#
18
19	using "streams";
20
21	Model Properties_Average
22
23	ATTRIBUTES
24	Pallete = false;
25	Brief = "Average physical properties of the streams.";
26	Info =
27	"to be documented.";
28
29	VARIABLES
30	Mw as molweight (Brief="Average Mol Weight",Default=75, Lower=1, Upper=1e8);
31	T as temperature (Brief="Average Temperature",Lower=50);
32	P as pressure (Brief="Average Pressure",Default=1, Lower=1e-10, Upper=2e4, DisplayUnit='kPa');
33	rho as dens_mass (Brief="Stream Density" ,Default=1000, Lower=1e-3, Upper=5e5, Symbol = "\rho");
34	Mu as viscosity (Brief="Stream Viscosity",Lower=0.0001, Symbol = "\mu");
35	Cp as cp_mol (Brief="Stream Molar Heat Capacity", Upper=1e10);
36	K as conductivity (Brief="Stream Thermal Conductivity", Default=1.0, Lower=1e-5, Upper=500);
37
38	end
39
40	Model Properties_In_Out
41
42	ATTRIBUTES
43	Pallete = false;
44	Brief = "Inlet and outlet physical properties of the streams.";
45	Info =
46	"to be documented.";
47
48	VARIABLES
49	Fw as flow_mass (Brief="Stream Mass Flow");
50	rho as dens_mass (Brief="Stream Density" ,Default=1000, Lower=1e-3, Upper=5e5, Symbol = "\rho");
51	end
52
53	Model Properties_Wall
54
55	ATTRIBUTES
56	Pallete = false;
57	Brief = "Physical properties of the streams at wall temperature.";
58	Info =
59	"to be documented.";
60
61	VARIABLES
62
63	Mu as viscosity (Brief="Stream Viscosity",Default=1, Lower=1e-5, Upper=1e5, Symbol = "\mu");
64	Twall as temperature (Brief="Wall Temperature",Lower=50);
65
66	end
67
68	Model Physical_Properties
69
70	ATTRIBUTES
71	Pallete = false;
72	Brief = "to be documented";
73	Info =
74	"to be documented";
75
76	VARIABLES
77
78	Inlet as Properties_In_Out (Brief="Properties at Inlet Stream", Symbol = "_{in}");
79	Average as Properties_Average (Brief="Properties at Average Temperature", Symbol = "_{avg}");
80	Outlet as Properties_In_Out (Brief="Properties at Outlet Stream", Symbol = "_{out}");
81	Wall as Properties_Wall (Brief="Properties at Wall Temperature", Symbol = "_{wall}");
82
83	end
84
85	Model Physical_Properties_Heatex
86
87	ATTRIBUTES
88	Pallete = false;
89	Brief = "to be documented";
90	Info =
91	"to be documented";
92
93	VARIABLES
94
95	Mw as molweight (Brief="Average Mol Weight",Default=75, Lower=1, Upper=1e8);
96	Cp as cp_mol (Brief="Average Molar Heat Capacity", Upper=1e10);
97	Inlet as Properties_In_Out (Brief="Properties at Inlet Stream", Symbol = "_{in}");
98	Outlet as Properties_In_Out (Brief="Properties at Outlet Stream", Symbol = "_{out}");
99
100	end
101
102	Model Tube_Pdrop
103
104	ATTRIBUTES
105	Pallete = false;
106	Brief = "Pressure drop and velocities in the tube side section of a shell and tube heat exchanger.";
107
108	VARIABLES
109	TubeFriction as press_delta (Brief="Tube Pressure Drop due to friction", Symbol = "\Delta P_{tube}", Default=0.01, Lower=1E-10,DisplayUnit='kPa');
110	InletNozzle as press_delta (Brief="Inlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_In}", Default=0.01, Lower=0,DisplayUnit='kPa');
111	OutletNozzle as press_delta (Brief="Outlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_Out}", Default=0.01, Lower=0,DisplayUnit='kPa');
112	Total as press_delta (Brief="Total Pressure Drop", Symbol = "\Delta P_{total}", Default=0.01, Lower=1E-10,DisplayUnit='kPa');
113	Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity", Symbol = "V_{Nozzle\_In}", Default=1, Upper=1e5, Lower=0);
114	Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity", Symbol = "V_{Nozzle\_Out}", Default=1, Upper=1E5, Lower=0);
115	FricFactor as fricfactor (Brief="Friction Factor", Symbol = "f_i", Default=0.05, Lower=1e-10, Upper=2000);
116
117	EQUATIONS
118	"Total Pressure Drop"
119	Total = TubeFriction + InletNozzle + OutletNozzle;
120
121	end
122
123	Model Tube_Heat_Transfer
124
125	ATTRIBUTES
126	Pallete = false;
127	Brief = "to be documented";
128	Info =
129	"to be documented";
130
131	VARIABLES
132	Re as positive (Brief="Tube Side Reynolds Number",Default=1000,Lower=1);
133	Nu as positive (Brief="Nusselt Number",Default=0.5,Lower=1e-8);
134	htube as heat_trans_coeff (Brief="Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
135	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
136	PR as positive (Brief="Tube Side Prandtl Number",Default=0.5,Lower=1e-8);
137	Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3);
138	Vtube as velocity (Brief="Tube Side Velocity",Lower=1e-8);
139
140	end
141
142	Model Shell_Pdrop
143
144	ATTRIBUTES
145	Pallete = false;
146	Brief = "Pressure drop and velocities in the shell side section of a shell and tube heat exchanger.";
147
148	VARIABLES
149	Ideal as press_delta (Brief="Ideal Pressure Drop", Symbol = "\Delta P_{ideal}", Default=0.01, Lower=0,DisplayUnit='kPa');
150	CrossFlow as press_delta (Brief="Cross Flow Pressure Drop", Symbol = "\Delta P_{CrossFlow}", Default=0.01, Lower=0,DisplayUnit='kPa');
151	EndZones as press_delta (Brief="End Zones Pressure Drop", Symbol = "\Delta P_{EndZones}", Default=0.01, Lower=0,DisplayUnit='kPa');
152	Window as press_delta (Brief="Window Pressure Drop", Symbol = "\Delta P_{Window}", Default=0.01, Lower=1e-10,DisplayUnit='kPa');
153	InletNozzle as press_delta (Brief="Inlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_In}", Default=0.01, Lower=0,DisplayUnit='kPa');
154	OutletNozzle as press_delta (Brief="Outlet Nozzle Pressure Drop", Symbol = "\Delta P_{Nozzle\_Out}", Default=0.01, Lower=0,DisplayUnit='kPa');
155	Total as press_delta (Brief="Total Pressure Drop", Symbol = "\Delta P_{Total}", Default=0.01, Lower=0,DisplayUnit='kPa');
156	FricFactor as fricfactor (Brief="Friction Factor", Symbol = "f_i", Default=0.05, Lower=1e-10, Upper=2000);
157	Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity", Symbol = "V_{Nozzle\_In}", Default=1, Upper=1e5, Lower=0);
158	Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity", Symbol = "V_{Nozzle\_Out}", Default=1, Upper=1e5, Lower=0);
159	RVsquare_out as positive (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
160	RVsquare_in as positive (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
161
162	EQUATIONS
163	"Shell Side Total Pressure Drop"
164	Total = CrossFlow+ EndZones + InletNozzle + OutletNozzle + Window;
165
166	end
167
168	Model Shell_Heat_Transfer
169
170	ATTRIBUTES
171	Pallete = false;
172	Brief = "to be documented";
173	Info =
174	"to be documented";
175
176	VARIABLES
177	Re as positive (Brief="Shell Side Reynolds Number",Default=100,Lower=1);
178	PR as positive (Brief="Shell Side Prandtl Number",Default=0.7,Lower=1e-6);
179	hshell as heat_trans_coeff (Brief="Shell Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
180	Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3);
181	end
182
183	Model Baffles_Main
184
185	ATTRIBUTES
186	Pallete = false;
187	Brief = "Main variables in the Baffle section of a shell and tube heat exchanger.";
188
189	PARAMETERS
190
191	BaffleCut as Integer (Brief="Baffle Cut",Default=25,Lower=25);
192	NumberOfBaffles as Real (Brief="Number of Baffles", Symbol = "N_{baffles}", Lower=1);
193
194	VARIABLES
195
196	Inlet_Spacing as length (Brief="Inlet Baffle Spacing",Lower=1e-8, Symbol = "L_{si}", DisplayUnit ='mm' );
197	Central_Spacing as length (Brief="Central Baffle Spacing",Lower=1e-8, Symbol = "L_s", DisplayUnit ='mm' );
198	Outlet_Spacing as length (Brief="Outlet Baffle Spacing",Lower=1e-8, Symbol = "L_{so}", DisplayUnit ='mm' );
199
200	end
201
202	Model Clearances_Main
203
204	ATTRIBUTES
205	Pallete = false;
206	Brief = "Main parameters for diametral clearances in a shell and tube heat exchanger.";
207
208	PARAMETERS
209
210	SealStrip as Integer (Brief="Number of Sealing Strips pairs",Lower=1);
211	Hinozzle_Shell as length (Brief="Height Under Shell Inlet Nozzle",Lower=1E-6);
212	Honozzle_Shell as length (Brief="Height Under Shell Outlet Nozzle",Lower=1E-6);
213	BundleToShell as length (Brief="Bundle-to-Shell Clearance", Symbol = "L_{cf}", Lower=1E-8);
214	BaffleToShell as length (Brief="Baffle-to-Shell Clearance", Symbol = "L_{cd}", Lower=1E-8);
215	TubeToBaffle as length (Brief="Tube-to-Baffle Clearance", Symbol = "L_{td}", Lower=1E-8);
216
217	end
218
219	Model NTU_Basic
220
221	ATTRIBUTES
222	Pallete = false;
223	Brief = "Number of Units Transference Method.";
224	Info =
225	"to be documented";
226
227	VARIABLES
228
229	Ch as positive (Brief="Hot Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
230	Cc as positive (Brief="Cold Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
231	Cr as positive (Brief="Heat Capacity Ratio",Default=0.5,Lower=1e-6);
232	Cmin as positive (Brief="Minimum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
233	Cmax as positive (Brief="Maximum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
234	NTU as positive (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
235	Eft as positive (Brief="Effectiveness",Default=0.5,Lower=1e-8,Upper=1, Symbol ="\varepsilon");
236	Eft1 as positive (Brief="Effectiveness Correction",Lower=1e-8,Default=0.5, Symbol ="\hat {\varepsilon}");
237
238	end
239
240	Model LMTD_Basic
241
242	ATTRIBUTES
243	Pallete = false;
244	Brief = "Log Mean Temperature Difference Method.";
245	Info =
246	"to be documented";
247
248	VARIABLES
249
250	DT0 as temp_delta (Brief="Temperature Difference at Inlet",Lower=1e-6, Symbol ="\Delta T_0");
251	DTL as temp_delta (Brief="Temperature Difference at Outlet",Lower=1e-6, Symbol ="\Delta T_L");
252	LMTD as temp_delta (Brief="Logarithmic Mean Temperature Difference",Lower=1e-6);
253	Fc as positive (Brief="LMTD Correction Factor",Lower=0.1);
254
255	EQUATIONS
256
257	if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL]))
258
259	then
260	"Log Mean Temperature Difference"
261	LMTD= (DT0-DTL)/ln(DT0/DTL);
262
263	else
264
265	if DT0*DTL equal 0
266
267	then
268	"Log Mean Temperature Difference"
269	LMTD = 0.5*(DT0+DTL);
270
271	else
272	"Log Mean Temperature Difference"
273	LMTD = 0.5(DT0+DTL)(1-(DT0-DTL)^2/(DT0DTL)(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
274
275	end
276
277	end
278
279	end
280
281	Model Details_Main
282
283	ATTRIBUTES
284	Pallete = false;
285	Brief = "to be documented";
286	Info =
287	"to be documented";
288
289	VARIABLES
290	A as area (Brief="Exchange Surface Area");
291	Q as power (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
292	Uc as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
293	Ud as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
294
295	end
296
297	Model Tube_Side_Main
298
299	ATTRIBUTES
300	Pallete = false;
301	Brief = "Main variables in the Tube Side section of a shell and tube heat exchanger.";
302
303	PARAMETERS
304
305	NumberOfTubes as Integer (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
306	Tubepasses as Integer (Brief="Number of Tube Passes", Lower=1);
307	TubeLength as length (Brief="Effective Tube Length",Lower=0.1);
308	TubePitch as length (Brief="Tube Pitch",Lower=1E-8);
309	Kwall as conductivity (Brief="Tube Wall Material Thermal Conductivity");
310	TubeOD as length (Brief="Tube Outside Diameter",Lower=1E-6);
311	TubeID as length (Brief="Tube Inside Diameter",Lower=1E-6);
312	Fouling as positive (Brief="Tubeside Fouling Resistance",Unit='m^2*K/kW', Symbol = "Rf_{tube}", Default=1E-6 , Lower=0);
313	InletNozzleID as length (Brief="Inlet Nozzle Inside Diameter", Lower=1E-6);
314	OutletNozzleID as length (Brief="Outlet Nozzle Inside Diameter", Lower=1E-6);
315
316	VARIABLES
317
318	PressureDrop as Tube_Pdrop (Brief="Tube Side Pressure Drop", Symbol = " ");
319	HeatTransfer as Tube_Heat_Transfer (Brief="Tube Side Heat Transfer", Symbol = " ");
320	Properties as Physical_Properties (Brief="Tube Side Properties", Symbol = " ");
321
322	end
323
324	Model Shell_Side_Main
325
326	ATTRIBUTES
327	Pallete = false;
328	Brief = "Main variables in the Shell Side section of a shell and tube heat exchanger.";
329
330	PARAMETERS
331
332	ShellID as length (Brief="Inside Shell Diameter",Lower=1E-6);
333	Fouling as positive (Brief="Shellside Fouling Resistance",Unit='m^2*K/kW', Symbol = "Rf_{shell}", Default=1E-6 , Lower=0);
334	InletNozzleID as length (Brief="Inlet Nozzle Inside Diameter", Lower=1E-6);
335	OutletNozzleID as length (Brief="Outlet Nozzle Inside Diameter", Lower=1E-6);
336
337	VARIABLES
338
339	PressureDrop as Shell_Pdrop (Brief="Shell Side Pressure Drop", Symbol = " ");
340	HeatTransfer as Shell_Heat_Transfer (Brief= "Shell Side Heat Transfer", Symbol = " ");
341	Properties as Physical_Properties (Brief="ShellSide Properties", Symbol = " ");
342
343	end
344
345	Model DoublePipe_HeatTransfer
346
347	ATTRIBUTES
348	Pallete = false;
349	Brief = "to be documented";
350	Info =
351	"to be documented";
352
353	PARAMETERS
354	As as area (Brief="Cross Sectional Area for Flow",Default=0.05,Lower=1e-8);
355	Dh as length (Brief="Hydraulic Diameter of Pipe for Heat Transfer",Lower=1e-8);
356
357	VARIABLES
358	Re as positive (Brief="Reynolds Number",Default=100,Lower=1);
359	hcoeff as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
360	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
361	Nu as positive (Brief="Nusselt Number",Default=0.5,Lower=1e-8);
362	PR as positive (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
363	Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3);
364	Vmean as velocity (Brief="Tube Velocity",Lower=1e-8);
365
366	end
367
368	Model DoublePipe_PressureDrop
369
370	ATTRIBUTES
371	Pallete = false;
372	Brief = "to be documented";
373	Info =
374	"to be documented";
375
376	PARAMETERS
377
378	Dh as length (Brief="Hydraulic Diameter of Pipe for Pressure Drop",Lower=1e-6);
379
380	VARIABLES
381
382	Pdrop as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa', Symbol ="\Delta P");
383	Pd_fric as press_delta (Brief="Pressure Drop for friction",Default=0.01, Lower=0,DisplayUnit='kPa', Symbol ="\Delta P_{fric}");
384	Pd_ret as press_delta (Brief="Pressure Drop due to return",Default=0.01, Lower=0,DisplayUnit='kPa', Symbol ="\Delta P_{return}");
385	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
386	Re as positive (Brief="Reynolds Number",Default=100,Lower=1);
387
388	end
389
390	Model Main_DoublePipe
391
392	ATTRIBUTES
393	Pallete = false;
394	Brief = "to be documented";
395	Info =
396	"to be documented";
397
398	VARIABLES
399
400	HeatTransfer as DoublePipe_HeatTransfer (Brief="Double Pipe Heat Transfer",Symbol=" ");
401	PressureDrop as DoublePipe_PressureDrop (Brief="Double Pipe Pressure Drop", Symbol=" ");
402	Properties as Physical_Properties (Brief="Double Pipe Properties",Symbol=" " );
403
404	end
405

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

Download in other formats:

Original Format