Context Navigation

HEX_Engine.mso @ 490

Last change on this file since 490 was 490, checked in by gerson bicca, 16 years ago
updated
Property svn:eol-style set to `native` Property svn:keywords set to `Id`
File size: 16.7 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 490 2008-04-03 23:38:09Z 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 Tube_Pdrop
86
87	ATTRIBUTES
88	Pallete = false;
89	Brief = "to be documented";
90	Info =
91	"to be documented";
92
93	VARIABLES
94	PdTube as press_delta (Brief="Tube Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa');
95	Pdtotal as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa');
96	Pdnozzle_in as press_delta (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
97	Pdnozzle_out as press_delta (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
98	Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
99	Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
100	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
101
102	EQUATIONS
103	"Total Pressure Drop"
104	Pdtotal = PdTube + Pdnozzle_in + Pdnozzle_out;
105
106	end
107
108	Model Tube_Heat_Transfer
109
110	ATTRIBUTES
111	Pallete = false;
112	Brief = "to be documented";
113	Info =
114	"to be documented";
115
116	VARIABLES
117	Re as positive (Brief="Tube Side Reynolds Number",Default=1000,Lower=1);
118	Nu as positive (Brief="Nusselt Number",Default=0.5,Lower=1e-8);
119	htube as heat_trans_coeff (Brief="Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
120	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
121	PR as positive (Brief="Tube Side Prandtl Number",Default=0.5,Lower=1e-8);
122	Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3);
123	Vtube as velocity (Brief="Tube Side Velocity",Lower=1e-8);
124
125	end
126
127	Model Shell_Pdrop
128
129	ATTRIBUTES
130	Pallete = false;
131	Brief = "to be documented";
132	Info =
133	"to be documented";
134
135	VARIABLES
136	Pideal as press_delta (Brief="Ideal Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
137	PdCross as press_delta (Brief="Cross Flow Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
138	PdEndZones as press_delta (Brief="End Zones Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
139	Pdwindow as press_delta (Brief="Window Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa');
140	Pdtotal as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
141	Pdnozzle_in as press_delta (Brief="Inlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
142	Pdnozzle_out as press_delta (Brief="Outlet Nozzle Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa');
143	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
144	Vnozzle_in as velocity (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
145	Vnozzle_out as velocity (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
146	RVsquare_out as positive (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
147	RVsquare_in as positive (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = 'kg/s^2/m');
148
149	EQUATIONS
150	"Shell Side Total Pressure Drop"
151	Pdtotal = PdCross + PdEndZones + Pdnozzle_in + Pdnozzle_out + Pdwindow;
152
153	end
154
155	Model Shell_Heat_Transfer
156
157	ATTRIBUTES
158	Pallete = false;
159	Brief = "to be documented";
160	Info =
161	"to be documented";
162
163	VARIABLES
164	Re as positive (Brief="Shell Side Reynolds Number",Default=100,Lower=1);
165	PR as positive (Brief="Shell Side Prandtl Number",Default=0.7,Lower=1e-6);
166	hshell as heat_trans_coeff (Brief="Shell Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
167	Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3);
168	end
169
170	Model Baffles_Main
171
172	ATTRIBUTES
173	Pallete = false;
174	Brief = "Main variables in the Baffle section of a shell and tube heat exchanger.";
175
176	PARAMETERS
177
178	BaffleCut as Integer (Brief="Baffle Cut",Default=25,Lower=25);
179	NumberOfBaffles as Real (Brief="Number of Baffles", Symbol = "N_{baffles}", Lower=1);
180
181	VARIABLES
182
183	Inlet_Spacing as length (Brief="Inlet Baffle Spacing",Lower=1e-8, Symbol = "L_{si}", DisplayUnit ='mm' );
184	Central_Spacing as length (Brief="Central Baffle Spacing",Lower=1e-8, Symbol = "L_s", DisplayUnit ='mm' );
185	Outlet_Spacing as length (Brief="Outlet Baffle Spacing",Lower=1e-8, Symbol = "L_{so}", DisplayUnit ='mm' );
186
187	end
188
189	Model Clearances_Main
190
191	ATTRIBUTES
192	Pallete = false;
193	Brief = "Main parameters for diametral clearances in a shell and tube heat exchanger.";
194
195	PARAMETERS
196
197	SealStrip as Integer (Brief="Number of Sealing Strips pairs",Lower=1);
198	Hinozzle_Shell as length (Brief="Height Under Shell Inlet Nozzle",Lower=1E-6);
199	Honozzle_Shell as length (Brief="Height Under Shell Outlet Nozzle",Lower=1E-6);
200	BundleToShell as length (Brief="Bundle-to-Shell Clearance", Symbol = "L_{cf}", Lower=1E-8);
201	BaffleToShell as length (Brief="Baffle-to-Shell Clearance", Symbol = "L_{cd}", Lower=1E-8);
202	TubeToBaffle as length (Brief="Tube-to-Baffle Clearance", Symbol = "L_{td}", Lower=1E-8);
203
204	end
205
206	Model NTU_Basic
207
208	ATTRIBUTES
209	Pallete = false;
210	Brief = "Number of Units Transference Method.";
211	Info =
212	"to be documented";
213
214	VARIABLES
215
216	Ch as positive (Brief="Hot Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
217	Cc as positive (Brief="Cold Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
218	Cr as positive (Brief="Heat Capacity Ratio",Default=0.5,Lower=1e-6);
219	Cmin as positive (Brief="Minimum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
220	Cmax as positive (Brief="Maximum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
221	NTU as positive (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
222	Eft as positive (Brief="Effectiveness",Default=0.5,Lower=1e-8,Upper=1, Symbol ="\varepsilon");
223	Eft1 as positive (Brief="Effectiveness Correction",Lower=1e-8,Default=0.5, Symbol ="\hat {\varepsilon}");
224
225	end
226
227	Model LMTD_Basic
228
229	ATTRIBUTES
230	Pallete = false;
231	Brief = "Log Mean Temperature Difference Method.";
232	Info =
233	"to be documented";
234
235	VARIABLES
236
237	DT0 as temp_delta (Brief="Temperature Difference at Inlet",Lower=1e-6, Symbol ="\Delta T_0");
238	DTL as temp_delta (Brief="Temperature Difference at Outlet",Lower=1e-6, Symbol ="\Delta T_L");
239	LMTD as temp_delta (Brief="Logarithmic Mean Temperature Difference",Lower=1e-6);
240	Fc as positive (Brief="LMTD Correction Factor",Lower=0.1);
241
242	EQUATIONS
243
244	if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL]))
245
246	then
247	"Log Mean Temperature Difference"
248	LMTD= (DT0-DTL)/ln(DT0/DTL);
249
250	else
251
252	if DT0*DTL equal 0
253
254	then
255	"Log Mean Temperature Difference"
256	LMTD = 0.5*(DT0+DTL);
257
258	else
259	"Log Mean Temperature Difference"
260	LMTD = 0.5(DT0+DTL)(1-(DT0-DTL)^2/(DT0DTL)(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
261
262	end
263
264	end
265
266	end
267
268	Model Details_Main
269
270	ATTRIBUTES
271	Pallete = false;
272	Brief = "to be documented";
273	Info =
274	"to be documented";
275
276	VARIABLES
277	A as area (Brief="Exchange Surface Area");
278	Q as power (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
279	Uc as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
280	Ud as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
281
282	end
283
284	Model Tube_Side_Main
285
286	ATTRIBUTES
287	Pallete = false;
288	Brief = "Main variables in the Tube Side section of a shell and tube heat exchanger.";
289
290	PARAMETERS
291
292	NumberOfTubes as Integer (Brief="Total Number of Tubes in Shell",Default=100,Lower=1);
293	Tubepasses as Integer (Brief="Number of Tube Passes", Lower=1);
294	TubeLength as length (Brief="Effective Tube Length",Lower=0.1);
295	TubePitch as length (Brief="Tube Pitch",Lower=1E-8);
296	Kwall as conductivity (Brief="Tube Wall Material Thermal Conductivity");
297	TubeOD as length (Brief="Tube Outside Diameter",Lower=1E-6);
298	TubeID as length (Brief="Tube Inside Diameter",Lower=1E-6);
299	Fouling as positive (Brief="Tubeside Fouling Resistance",Unit='m^2*K/kW', Symbol = "Rf_{tube}", Default=1E-6 , Lower=0);
300	InletNozzleID as length (Brief="Inlet Nozzle Inside Diameter", Lower=1E-6);
301	OutletNozzleID as length (Brief="Outlet Nozzle Inside Diameter", Lower=1E-6);
302
303	VARIABLES
304
305	PressureDrop as Tube_Pdrop (Brief="Tube Side Pressure Drop");
306	HeatTransfer as Tube_Heat_Transfer (Brief="Tube Side Heat Transfer");
307	Properties as Physical_Properties (Brief="Tube Side Properties");
308
309	end
310
311	Model Shell_Side_Main
312
313	ATTRIBUTES
314	Pallete = false;
315	Brief = "Main variables in the Shell Side section of a shell and tube heat exchanger.";
316
317	PARAMETERS
318
319	ShellID as length (Brief="Inside Shell Diameter",Lower=1E-6);
320	Fouling as positive (Brief="Shellside Fouling Resistance",Unit='m^2*K/kW', Symbol = "Rf_{shell}", Default=1E-6 , Lower=0);
321	InletNozzleID as length (Brief="Inlet Nozzle Inside Diameter", Lower=1E-6);
322	OutletNozzleID as length (Brief="Outlet Nozzle Inside Diameter", Lower=1E-6);
323
324	VARIABLES
325
326	PressureDrop as Shell_Pdrop (Brief="Shell Side Pressure Drop");
327	HeatTransfer as Shell_Heat_Transfer (Brief= "Shell Side Heat Transfer");
328	Properties as Physical_Properties (Brief="ShellSide Properties");
329
330	end
331
332	Model Basic_Pdrop
333
334	ATTRIBUTES
335	Pallete = false;
336	Brief = "to be documented";
337	Info =
338	"to be documented";
339
340	VARIABLES
341
342	Pdrop as press_delta (Brief="Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa' , Symbol ="\Delta P");
343	FPdrop as Real (Brief="Pressure Drop : Fraction of Inlet",Lower=0,Upper=0.8);
344
345	end
346
347	Model DoublePipe_HeatTransfer
348
349	ATTRIBUTES
350	Pallete = false;
351	Brief = "to be documented";
352	Info =
353	"to be documented";
354
355	PARAMETERS
356	As as area (Brief="Cross Sectional Area for Flow",Default=0.05,Lower=1e-8);
357	Dh as length (Brief="Hydraulic Diameter of Pipe for Heat Transfer",Lower=1e-8);
358
359	VARIABLES
360	Re as positive (Brief="Reynolds Number",Default=100,Lower=1);
361	hcoeff as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
362	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
363	Nu as positive (Brief="Nusselt Number",Default=0.5,Lower=1e-8);
364	PR as positive (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
365	Phi as positive (Brief="Phi Correction",Default=1,Lower=1e-3);
366	Vmean as velocity (Brief="Tube Velocity",Lower=1e-8);
367
368	end
369
370	Model DoublePipe_PressureDrop
371
372	ATTRIBUTES
373	Pallete = false;
374	Brief = "to be documented";
375	Info =
376	"to be documented";
377
378	PARAMETERS
379
380	Dh as length (Brief="Hydraulic Diameter of Pipe for Pressure Drop",Lower=1e-6);
381
382	VARIABLES
383
384	Pdrop as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=0,DisplayUnit='kPa', Symbol ="\Delta P");
385	Pd_fric as press_delta (Brief="Pressure Drop for friction",Default=0.01, Lower=0,DisplayUnit='kPa', Symbol ="\Delta P_{fric}");
386	Pd_ret as press_delta (Brief="Pressure Drop due to return",Default=0.01, Lower=0,DisplayUnit='kPa', Symbol ="\Delta P_{return}");
387	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
388	Re as positive (Brief="Reynolds Number",Default=100,Lower=1);
389
390	end
391
392	Model Main_DoublePipe
393
394	ATTRIBUTES
395	Pallete = false;
396	Brief = "to be documented";
397	Info =
398	"to be documented";
399
400	VARIABLES
401
402	HeatTransfer as DoublePipe_HeatTransfer (Brief="Double Pipe Heat Transfer",Symbol=" ");
403	PressureDrop as DoublePipe_PressureDrop (Brief="Double Pipe Pressure Drop", Symbol=" ");
404	Properties as Physical_Properties (Brief="Double Pipe Properties",Symbol=" " );
405
406	end
407
408	Model Main_Simplified
409
410	ATTRIBUTES
411	Pallete = false;
412	Brief = "to be documented";
413	Info =
414	"to be documented";
415
416	VARIABLES
417
418	PressureDrop as Basic_Pdrop (Brief="Pressure Drop", Symbol=" ");
419	Properties as Physical_Properties (Brief="Physical Properties", Symbol=" ");
420
421	end
422
423	Model PHE_HeatTransfer
424
425	ATTRIBUTES
426	Pallete = false;
427	Brief = "to be documented";
428	Info =
429	"to be documented";
430
431	VARIABLES
432
433	Re as positive (Brief="Reynolds Number",Default=100,Lower=1);
434	PR as positive (Brief="Prandtl Number",Default=0.5,Lower=1e-8);
435	NTU as positive (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
436	WCp as positive (Brief="Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit='W/K');
437	hcoeff as heat_trans_coeff (Brief="Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
438	Gchannel as flux_mass (Brief ="Channel Mass Flux", Default=1, Lower=1e-6, Symbol ="G^{channel}");
439	Gports as flux_mass (Brief ="Ports Mass Flux", Default=1, Lower=1e-6, Symbol ="G^{ports}");
440	Phi as positive (Brief="Viscosity Correction",Default=1,Lower=1e-6, Symbol="\phi");
441
442	end
443
444	Model PHE_PressureDrop
445
446	ATTRIBUTES
447	Pallete = false;
448	Brief = "to be documented";
449	Info =
450	"to be documented";
451
452	VARIABLES
453
454	DPchannel as press_delta (Brief="Channel Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa', Symbol ="\Delta P^{channel}");
455	DPports as press_delta (Brief="Ports Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa', Symbol ="\Delta P^{ports}");
456	Pdrop as press_delta (Brief="Total Pressure Drop",Default=0.01, Lower=1e-10,DisplayUnit='kPa', Symbol ="\Delta P");
457	fi as fricfactor (Brief="Friction Factor", Default=0.05, Lower=1e-10, Upper=2000);
458	Vchannel as velocity (Brief="Stream Velocity in Channel",Lower=1e-8, Symbol ="V^{channel}");
459	Vports as velocity (Brief="Stream Velocity in Ports",Lower=1e-8, Symbol ="V^{ports}");
460	Npassage as Real (Brief="Number of Channels per Pass", Symbol ="N^{passage}");
461	end
462
463	Model Main_PHE
464
465	ATTRIBUTES
466	Pallete = false;
467	Brief = "to be documented";
468	Info =
469	"to be documented";
470
471	VARIABLES
472
473	HeatTransfer as PHE_HeatTransfer (Brief="PHE Heat Transfer", Symbol = " ");
474	PressureDrop as PHE_PressureDrop (Brief="PHE Pressure Drop", Symbol = " ");
475	Properties as Physical_Properties (Brief="PHE Properties", Symbol = " ");
476
477	end
478
479	Model Thermal_PHE
480
481	ATTRIBUTES
482	Pallete = false;
483	Brief = "to be documented";
484	Info =
485	"to be documented";
486
487	VARIABLES
488	Cr as positive (Brief="Heat Capacity Ratio",Default=0.5,Lower=1e-6);
489	Cmin as positive (Brief="Minimum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
490	Cmax as positive (Brief="Maximum Heat Capacity",Lower=1e-10,Default=1e3,Unit='W/K');
491	NTU as positive (Brief="Number of Units Transference",Default=0.05,Lower=1e-10);
492	Eft as positive (Brief="Effectiveness",Default=0.5,Lower=0.4,Upper=1, Symbol = "\varepsilon");
493	Q as power (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
494	Uc as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
495	Ud as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
496	end

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

Download in other formats:

Original Format