source: trunk/BioModel/heat_exchangers/heater_m.mso @ 1008

Last change on this file since 1008 was 1008, checked in by Argimiro Resende Secchi, 21 months ago

Adding BioModel? to the MSO library.
File size: 27.9 KB
Line 
1#*-------------------------------------------------------------------
2* Biorrefinaria Petrobras
3*--------------------------------------------------------------------
4* Nome do arquivo: heater_m.mso
5* Projeto: Modelo integrado de producao de etanol 1G/2G
6* Conteudo: trocador de calor
7*--------------------------------------------------------------------*#
8
9#*-------------------------------------------------------------------
10*
11* Versao 2.2
12* Data:    03/2016
13* Autores:   Anderson R. A. Lino e Gabriel C. Fonseca
14*
15*--------------------------------------------------------------------
16* Descricao: modelo simplificado do trocador de calor que sera 
17* empregado na biorrefinaria
18*--------------------------------------------------------------------
19
20*--------------------------------------------------------------------
21*Hipoteses assumidas: 1 - operacao em estado estacionario
22*                     2 - sem perdas de calor para o ambiente
23*--------------------------------------------------------------------
24
25*--------------------------------------------------------------------
26*Notas: Foi feito o flowsheet teste para averiguar o modelo.
27*As siglas R, L, U, D sao indicativos da posicao da corrente de saida
28*sendo R=right, L=legth, U=up e D=down
29*--------------------------------------------------------------------*#
30
31using "energy_stream";
32using "main_stream";
33
34Model heatex_mL
35        ATTRIBUTES
36        Pallete = false;
37        Brief = "Model of a Heat Exchanger";
38        Info =
39"== GENERAL ==
40* Model of a simplified heat exchanger;
41* The model performes mass and energy balances only.
42       
43== ASSUMPTIONS ==
44* Steady-State operation;
45* No heat loss to the surroundings.
46";
47
48
49#*-------------------------------------------------------------------
50* Declaracao de variaveis
51*--------------------------------------------------------------------*#
52       
53VARIABLES
54
55in  Inlet               as main_stream          (Brief="Inlet Stream", PosX=1.0, PosY=0.5, Protected = false, Symbol="_{in}");
56out Outlet              as main_stream_PH       (Brief="Outlet Stream", PosX=0.0, PosY=0.5, Protected = false, Symbol="_{out}");
57        Q                       as power                        (Brief="Heat Duty", Default=0, Lower=-1e10, Upper=1e10);
58       
59        Pdrop           as press_delta          (Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
60        lmtd            as temp_delta           (Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
61        U                       as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
62        A                       as area                         (Brief = "Exchange Surface Area", Lower = 1e-6);
63       
64
65#*-------------------------------------------------------------------
66* Equacoes do modelo
67*--------------------------------------------------------------------*#
68
69        EQUATIONS
70
71        "Energy Balance"
72                Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);
73
74        "Molar Balance (Fluid Phase)"
75                Inlet.Fluid.F  = Outlet.Fluid.F;
76               
77        "Molar Balance (Solid Phase)"
78                Inlet.Solid.F = Outlet.Solid.F;
79
80        "Fraction Constraint (Fluid Phase)"
81                Outlet.Fluid.z = Inlet.Fluid.z;
82               
83        "Fraction Constraint (Solid Phase)"
84                Outlet.Solid.z = Inlet.Solid.z;
85               
86        "Pressure Drop"
87                Outlet.P = Inlet.P - Pdrop;
88       
89end
90
91Model heater_mLD as heatex_mL
92        ATTRIBUTES
93        Pallete         = true;
94        Icon            = "icon/heaterD";
95        Brief           = "Model of a Heat Exchanger";
96        Info =
97"== GENERAL ==
98*Model of a simplified heat exchanger;
99*The model performes mass and energy balances only.
100
101== ASSUMPTIONS ==
102* Steady-State operation;
103* No heat loss to the surroundings.
104
105== SPECIFY ==
106* The inlet streams:
107  flow rate
108  temperature
109  pressure
110  stream composition;
111* Outlet stream temperature;
112* Overall Heat Transfer Coefficient;
113* Exchange Surface Area;
114* Pressure drop.
115
116== SET ==
117* Phase of the fluid entering the heater;
118* Number of stream components(Ncomp/NcompS);
119* Basic composition (mass or molar).
120";
121
122        VARIABLES
123in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{in}");
124
125        EQUATIONS
126       
127        "Equate Heat Stream"
128        Inlet_q.Q = Q;
129       
130        "General equation of heat exchange"
131        Inlet_q.Q = U * A * lmtd;
132       
133end
134
135Model cooler_mLD as heatex_mL
136        ATTRIBUTES
137        Pallete         = true;
138        Icon            = "icon/heaterD";
139        Brief           = "Model of a Heat Exchanger";
140        Info =
141"== GENERAL ==
142*Model of a simplified heat exchanger;
143*The model performes mass and energy balances only.
144
145== ASSUMPTIONS ==
146* Steady-State operation;
147* No heat loss to the surroundings.
148
149== SPECIFY ==
150* The inlet streams:
151  flow rate
152  temperature
153  pressure
154  stream composition;
155* Outlet vapour fraction;
156* Overall Heat Transfer Coefficient;
157* Exchange Surface Area;
158* Pressure drop.
159
160== SET ==
161* Phase of the fluid entering the heater;
162* Number of stream components(Ncomp/NcompS);
163* Basic composition (mass or molar).
164";
165       
166        VARIABLES
167out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{out}");
168
169        EQUATIONS
170        "Equate Heat Stream"
171        Outlet_q.Q = -Q;
172       
173        "General equation of heat exchange"
174        Outlet_q.Q = U * A * lmtd;
175       
176end
177
178Model heater_mLU as heatex_mL
179        ATTRIBUTES
180        Pallete         = true;
181        Icon            = "icon/heaterU";
182        Brief           = "Model of a Heat Exchanger";
183        Info =
184"== GENERAL ==
185*Model of a simplified heat exchanger;
186*The model performes mass and energy balances only.
187
188== ASSUMPTIONS ==
189* Steady-State operation;
190* No heat loss to the surroundings.
191
192== SPECIFY ==
193* The inlet streams:
194  flow rate
195  temperature
196  pressure
197  stream composition;
198* Outlet stream temperature;
199* Overall Heat Transfer Coefficient;
200* Exchange Surface Area;
201* Pressure drop.
202
203== SET ==
204* Phase of the fluid entering the heater;
205* Number of stream components(Ncomp/NcompS);
206* Basic composition (mass or molar).
207";
208       
209        VARIABLES
210in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{in}");
211
212        EQUATIONS
213        "Equate Heat Stream"
214        Inlet_q.Q = Q;
215       
216        "General equation of heat exchange"
217        Inlet_q.Q = U * A * lmtd;
218       
219end
220
221Model cooler_mLU as heatex_mL
222        ATTRIBUTES
223        Pallete         = true;
224        Icon            = "icon/heaterU";
225        Brief           = "Model of a Heat Exchanger";
226        Info =
227"== GENERAL ==
228*Model of a simplified heat exchanger;
229*The model performes mass and energy balances only.
230
231== ASSUMPTIONS ==
232* Steady-State operation;
233* No heat loss to the surroundings.
234
235== SPECIFY ==
236* The inlet streams:
237  flow rate
238  temperature
239  pressure
240  stream composition;
241* Outlet vapour fraction;
242* Overall Heat Transfer Coefficient;
243* Exchange Surface Area;
244* Pressure drop.
245
246== SET ==
247* Phase of the fluid entering the heater;
248* Number of stream components(Ncomp/NcompS);
249* Basic composition (mass or molar).
250";
251       
252        VARIABLES
253out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{out}");
254
255        EQUATIONS
256        "Equate Heat Stream"
257        Outlet_q.Q = -Q;
258       
259        "General equation of heat exchange"
260        Outlet_q.Q = U * A * lmtd;
261       
262end
263
264Model heatex_mR
265
266        ATTRIBUTES
267        Pallete         = false;
268        Brief           = "Model of a Heat Exchanger";
269        Info =
270"== GENERAL ==
271* Model of a simplified heat exchanger;
272* The model performes mass and energy balances only.
273       
274== ASSUMPTIONS ==
275* Steady-State operation;
276* No heat loss to the surroundings.
277";
278
279#*-------------------------------------------------------------------
280* Declaracao de variaveis
281*--------------------------------------------------------------------*#
282       
283VARIABLES
284
285in  Inlet               as main_stream          (Brief="Inlet Stream", PosX=0.0, PosY=0.5, Protected = false, Symbol="_{in}");
286out Outlet              as main_stream_PH       (Brief="Outlet Stream", PosX=1.0, PosY=0.5, Protected = false, Symbol="_{out}");
287        Q                       as power                        (Brief="Heat Duty", Default=0, Lower=-1e10, Upper=1e10);
288       
289        Pdrop           as press_delta          (Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
290        lmtd            as temp_delta           (Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
291        U                       as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
292        A                       as area                         (Brief = "Exchange Surface Area", Lower = 1e-6);
293       
294#*-------------------------------------------------------------------
295* Equacoes do modelo
296*--------------------------------------------------------------------*#
297
298        EQUATIONS
299
300        "Energy Balance"
301                Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);
302
303        "Molar Balance (Fluid Phase)"
304                Inlet.Fluid.F  = Outlet.Fluid.F;
305               
306        "Molar Balance (Solid Phase)"
307                Inlet.Solid.F = Outlet.Solid.F;
308
309        "Fraction Constraint (Fluid Phase)"
310                Outlet.Fluid.z  = Inlet.Fluid.z;
311               
312        "Fraction Constraint (Solid Phase)"
313                Outlet.Solid.z = Inlet.Solid.z;
314               
315        "Pressure Drop"
316                Outlet.P = Inlet.P - Pdrop;
317end
318
319Model heater_mRD as heatex_mR
320        ATTRIBUTES
321        Pallete         = true;
322        Icon            = "icon/heaterD";
323        Brief           = "Model of a Heat Exchanger";
324        Info =
325"== GENERAL ==
326*Model of a simplified heat exchanger;
327*The model performes mass and energy balances only.
328
329== ASSUMPTIONS ==
330* Steady-State operation;
331* No heat loss to the surroundings.
332
333== SPECIFY ==
334* The inlet streams:
335  flow rate
336  temperature
337  pressure
338  stream composition;
339* Outlet stream temperature;
340* Overall Heat Transfer Coefficient;
341* Exchange Surface Area;
342* Pressure drop.
343
344== SET ==
345* Phase of the fluid entering the heater;
346* Number of stream components(Ncomp/NcompS);
347* Basic composition (mass or molar).
348";
349       
350        VARIABLES
351in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{in}");
352
353        EQUATIONS
354        "Equate Heat Stream"
355        Inlet_q.Q = Q;
356       
357        "General equation of heat exchange"
358        Inlet_q.Q = U * A * lmtd;
359       
360end
361
362Model cooler_mRD as heatex_mR
363        ATTRIBUTES
364        Pallete         = true;
365        Icon            = "icon/heaterD";
366        Brief           = "Model of a Heat Exchanger";
367        Info =
368"== GENERAL ==
369*Model of a simplified heat exchanger;
370*The model performes mass and energy balances only.
371
372== ASSUMPTIONS ==
373* Steady-State operation;
374* No heat loss to the surroundings.
375
376== SPECIFY ==
377* The inlet streams:
378  flow rate
379  temperature
380  pressure
381  stream composition;
382* Outlet vapour fraction;
383* Overall Heat Transfer Coefficient;
384* Exchange Surface Area;
385* Pressure drop.
386
387== SET ==
388* Phase of the fluid entering the heater;
389* Number of stream components(Ncomp/NcompS);
390* Basic composition (mass or molar).
391";
392       
393        VARIABLES
394out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=0.5, PosY=1.0, Symbol="_{out}");
395
396        EQUATIONS
397        "Equate Heat Stream"
398        Outlet_q.Q = -Q;
399       
400        "General equation of heat exchange"
401        Outlet_q.Q = U * A * lmtd;
402       
403end
404
405Model heater_mRU as heatex_mR
406        ATTRIBUTES
407        Pallete         = true;
408        Icon            = "icon/heaterU";
409        Brief           = "Model of a Heat Exchanger";
410        Info =
411"== GENERAL ==
412*Model of a simplified heat exchanger;
413*The model performes mass and energy balances only.
414
415== ASSUMPTIONS ==
416* Steady-State operation;
417* No heat loss to the surroundings.
418
419== SPECIFY ==
420* The inlet streams:
421  flow rate
422  temperature
423  pressure
424  stream composition;
425* Outlet stream temperature;
426* Overall Heat Transfer Coefficient;
427* Exchange Surface Area;
428* Pressure drop.
429
430== SET ==
431* Phase of the fluid entering the heater;
432* Number of stream components(Ncomp/NcompS);
433* Basic composition (mass or molar).
434";
435       
436        VARIABLES
437in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{in}");
438
439        EQUATIONS
440        "Equate Heat Stream"
441        Inlet_q.Q = Q;
442       
443        "General equation of heat exchange"
444        Inlet_q.Q = U * A * lmtd;
445       
446end
447
448Model cooler_mRU as heatex_mR
449        ATTRIBUTES
450        Pallete         = true;
451        Icon            = "icon/heaterU";
452        Brief           = "Model of a Heat Exchanger";
453        Info =
454"== GENERAL ==
455*Model of a simplified heat exchanger;
456*The model performes mass and energy balances only.
457
458== ASSUMPTIONS ==
459* Steady-State operation;
460* No heat loss to the surroundings.
461
462== SPECIFY ==
463* The inlet streams:
464  flow rate
465  temperature
466  pressure
467  stream composition;
468* Outlet vapour fraction;
469* Overall Heat Transfer Coefficient;
470* Exchange Surface Area;
471* Pressure drop.
472
473== SET ==
474* Phase of the fluid entering the heater;
475* Number of stream components(Ncomp/NcompS);
476* Basic composition (mass or molar).
477";
478       
479        VARIABLES
480out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=0.5, PosY=0.0, Symbol="_{out}");
481
482        EQUATIONS
483        "Equate Heat Stream"
484        Outlet_q.Q = -Q;
485       
486        "General equation of heat exchange"
487        Outlet_q.Q = U * A * lmtd;
488       
489end
490
491Model heatex_mU
492
493        ATTRIBUTES
494        Pallete         = false;
495        Brief           = "Model of a Heat Exchanger";
496        Info =
497"== GENERAL ==
498* Model of a simplified heat exchanger;
499* The model performes mass and energy balances only.
500       
501== ASSUMPTIONS ==
502* Steady-State operation;
503* No heat loss to the surroundings.
504";
505
506#*-------------------------------------------------------------------
507* Declaracao de variaveis
508*--------------------------------------------------------------------*#
509       
510VARIABLES
511
512in  Inlet               as main_stream          (Brief="Inlet Stream", PosX=0.5, PosY=1.0, Protected = false, Symbol="_{in}");
513out Outlet              as main_stream_PH       (Brief="Outlet Stream", PosX=0.5, PosY=0.0, Protected = false, Symbol="_{out}");
514        Q                       as power                        (Brief="Heat Duty", Default=7000, Lower=-1e10, Upper=1e10);
515       
516        Pdrop           as press_delta          (Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
517        lmtd            as temp_delta           (Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
518        U                       as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
519        A                       as area                         (Brief = "Exchange Surface Area", Lower = 1e-6);
520       
521#*-------------------------------------------------------------------
522* Equacoes do modelo
523*--------------------------------------------------------------------*#
524
525        EQUATIONS
526
527        "Energy Balance"
528                Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);
529
530        "Molar Balance (Fluid Phase)"
531                Inlet.Fluid.F  = Outlet.Fluid.F;
532               
533        "Molar Balance (Solid Phase)"
534                Inlet.Solid.F = Outlet.Solid.F;
535
536        "Fraction Constraint (Fluid Phase)"
537                Outlet.Fluid.z = Inlet.Fluid.z;
538               
539        "Fraction Constraint (Solid Phase)"
540                Outlet.Solid.z = Inlet.Solid.z;
541               
542        "Pressure Drop"
543                Outlet.P = Inlet.P - Pdrop;     
544end
545
546Model heater_mUR as heatex_mU
547        ATTRIBUTES
548        Pallete         = true;
549        Icon            = "icon/heaterR";
550        Brief           = "Model of a Heat Exchanger";
551        Info =
552"== GENERAL ==
553*Model of a simplified heat exchanger;
554*The model performes mass and energy balances only.
555
556== ASSUMPTIONS ==
557* Steady-State operation;
558* No heat loss to the surroundings.
559
560== SPECIFY ==
561* The inlet streams:
562  flow rate
563  temperature
564  pressure
565  stream composition;
566* Outlet stream temperature;
567* Overall Heat Transfer Coefficient;
568* Exchange Surface Area;
569* Pressure drop.
570
571== SET ==
572* Phase of the fluid entering the heater;
573* Number of stream components(Ncomp/NcompS);
574* Basic composition (mass or molar).
575";
576
577        VARIABLES
578in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{in}");
579
580        EQUATIONS
581        "Equate Heat Stream"
582        Inlet_q.Q = Q;
583       
584        "General equation of heat exchange"
585        Inlet_q.Q = U * A * lmtd;
586       
587end
588
589Model cooler_mUR as heatex_mU
590        ATTRIBUTES
591        Pallete         = true;
592        Icon            = "icon/heaterR";
593        Brief           = "Model of a Heat Exchanger";
594        Info =
595"== GENERAL ==
596*Model of a simplified heat exchanger;
597*The model performes mass and energy balances only.
598
599== ASSUMPTIONS ==
600* Steady-State operation;
601* No heat loss to the surroundings.
602
603== SPECIFY ==
604* The inlet streams:
605  flow rate
606  temperature
607  pressure
608  stream composition;
609* Outlet vapour fraction;
610* Overall Heat Transfer Coefficient;
611* Exchange Surface Area;
612* Pressure drop.
613
614== SET ==
615* Phase of the fluid entering the heater;
616* Number of stream components(Ncomp/NcompS);
617* Basic composition (mass or molar).
618";
619
620        VARIABLES
621out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{out}");
622
623        EQUATIONS
624        "Equate Heat Stream"
625        Outlet_q.Q = -Q;
626       
627        "General equation of heat exchange"
628        Outlet_q.Q = U * A * lmtd;
629       
630end
631
632Model heater_mUL as heatex_mU
633        ATTRIBUTES
634        Pallete         = true;
635        Icon            = "icon/heaterL";
636        Brief           = "Model of a Heat Exchanger";
637        Info =
638"== GENERAL ==
639*Model of a simplified heat exchanger;
640*The model performes mass and energy balances only.
641
642== ASSUMPTIONS ==
643* Steady-State operation;
644* No heat loss to the surroundings.
645
646== SPECIFY ==
647* The inlet streams:
648  flow rate
649  temperature
650  pressure
651  stream composition;
652* Outlet stream temperature;
653* Overall Heat Transfer Coefficient;
654* Exchange Surface Area;
655* Pressure drop.
656
657== SET ==
658* Phase of the fluid entering the heater;
659* Number of stream components(Ncomp/NcompS);
660* Basic composition (mass or molar).
661";
662       
663        VARIABLES
664in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{in}");
665
666        EQUATIONS
667        "Equate Heat Stream"
668        Inlet_q.Q = Q;
669       
670        "General equation of heat exchange"
671        Inlet_q.Q = U * A * lmtd;
672       
673end
674
675Model cooler_mUL as heatex_mU
676        ATTRIBUTES
677        Pallete         = true;
678        Icon            = "icon/heaterL";
679        Brief           = "Model of a Heat Exchanger";
680        Info =
681"== GENERAL ==
682*Model of a simplified heat exchanger;
683*The model performes mass and energy balances only.
684
685== ASSUMPTIONS ==
686* Steady-State operation;
687* No heat loss to the surroundings.
688
689== SPECIFY ==
690* The inlet streams:
691  flow rate
692  temperature
693  pressure
694  stream composition;
695* Outlet vapour fraction;
696* Overall Heat Transfer Coefficient;
697* Exchange Surface Area;
698* Pressure drop.
699
700== SET ==
701* Phase of the fluid entering the heater;
702* Number of stream components(Ncomp/NcompS);
703* Basic composition (mass or molar).
704";
705       
706        VARIABLES
707out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{out}");
708
709        EQUATIONS
710        "Equate Heat Stream"
711        Outlet_q.Q = -Q;
712       
713        "General equation of heat exchange"
714        Outlet_q.Q = U * A * lmtd;
715       
716end
717
718Model heatex_mD
719
720        ATTRIBUTES
721        Pallete         = false;
722        Brief           = "Model of a Heat Exchanger";
723        Info =
724"== GENERAL ==
725* Model of a simplified heat exchanger;
726* The model performes mass and energy balances only.
727       
728== ASSUMPTIONS ==
729* Steady-State operation;
730* No heat loss to the surroundings.
731";
732
733#*-------------------------------------------------------------------
734* Declaracao de variaveis
735*--------------------------------------------------------------------*#
736       
737VARIABLES
738
739in  Inlet               as main_stream          (Brief="Inlet Stream", PosX=0.5, PosY=0.0, Protected = false, Symbol="_{in}");
740out Outlet              as main_stream_PH       (Brief="Outlet Stream", PosX=0.5, PosY=1.0, Protected = false, Symbol="_{out}");
741        Q                       as power                        (Brief="Heat Duty", Default=0, Lower=-1e10, Upper=1e10);
742       
743        Pdrop           as press_delta          (Brief="Pressure Drop",Default=0.01, Lower=-0.001,DisplayUnit='kPa' , Symbol ="\Delta P");
744        lmtd            as temp_delta           (Brief = "Normalized Mean Log", Lower = 1e-6, Symbol = "LMTD");
745        U                       as heat_trans_coeff (Brief = "Overall Heat Transfer Coefficient", Default=0.69445, Lower=1e-6, Upper=1e10);
746        A                       as area                         (Brief = "Exchange Surface Area", Lower = 1e-6);
747       
748#*-------------------------------------------------------------------
749* Equacoes do modelo
750*--------------------------------------------------------------------*#
751
752        EQUATIONS
753
754        "Energy Balance"
755                Q = Inlet.Fluid.F*(Outlet.Fluid.h - Inlet.Fluid.h) + Inlet.Solid.F*(Outlet.Solid.h - Inlet.Solid.h);
756
757        "Molar Balance (Fluid Phase)"
758                Inlet.Fluid.F  = Outlet.Fluid.F;
759               
760        "Molar Balance (Solid Phase)"
761                Inlet.Solid.F = Outlet.Solid.F;
762
763        "Fraction Constraint (Fluid Phase)"
764                Outlet.Fluid.z = Inlet.Fluid.z;
765               
766        "Fraction Constraint (Solid Phase)"
767                Outlet.Solid.z = Inlet.Solid.z;
768               
769        "Pressure Drop"
770                Outlet.P = Inlet.P - Pdrop;     
771end
772
773Model heater_mDR as heatex_mD
774        ATTRIBUTES
775        Pallete         = true;
776        Icon            = "icon/heaterR";
777        Brief           = "Model of a Heat Exchanger";
778        Info =
779"== GENERAL ==
780*Model of a simplified heat exchanger;
781*The model performes mass and energy balances only.
782
783== ASSUMPTIONS ==
784* Steady-State operation;
785* No heat loss to the surroundings.
786
787== SPECIFY ==
788* The inlet streams:
789  flow rate
790  temperature
791  pressure
792  stream composition;
793* Outlet stream temperature;
794* Overall Heat Transfer Coefficient;
795* Exchange Surface Area;
796* Pressure drop.
797
798== SET ==
799* Phase of the fluid entering the heater;
800* Number of stream components(Ncomp/NcompS);
801* Basic composition (mass or molar).
802";
803       
804        VARIABLES
805in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{in}");
806
807        EQUATIONS
808        "Equate Heat Stream"
809        Inlet_q.Q = Q;
810       
811        "General equation of heat exchange"
812        Inlet_q.Q = U * A * lmtd;
813       
814end
815
816Model cooler_mDR as heatex_mD
817        ATTRIBUTES
818        Pallete         = true;
819        Icon            = "icon/heaterL";
820        Brief           = "Model of a Heat Exchanger";
821        Info =
822"== GENERAL ==
823*Model of a simplified heat exchanger;
824*The model performes mass and energy balances only.
825
826== ASSUMPTIONS ==
827* Steady-State operation;
828* No heat loss to the surroundings.
829
830== SPECIFY ==
831* The inlet streams:
832  flow rate
833  temperature
834  pressure
835  stream composition;
836* Outlet vapour fraction;
837* Overall Heat Transfer Coefficient;
838* Exchange Surface Area;
839* Pressure drop.
840
841== SET ==
842* Phase of the fluid entering the heater;
843* Number of stream components(Ncomp/NcompS);
844* Basic composition (mass or molar).
845";
846       
847        VARIABLES
848out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=1.0, PosY=0.5, Symbol="_{out}");
849
850        EQUATIONS
851        "Equate Heat Stream"
852        Outlet_q.Q = -Q;
853       
854        "General equation of heat exchange"
855        Outlet_q.Q = U * A * lmtd;
856       
857end
858
859Model heater_mDL as heatex_mD
860        ATTRIBUTES
861        Pallete         = true;
862        Icon            = "icon/heaterR";
863        Brief           = "Model of a Heat Exchanger";
864        Info =
865"== GENERAL ==
866*Model of a simplified heat exchanger;
867*The model performes mass and energy balances only.
868
869== ASSUMPTIONS ==
870* Steady-State operation;
871* No heat loss to the surroundings.
872
873== SPECIFY ==
874* The inlet streams:
875  flow rate
876  temperature
877  pressure
878  stream composition;
879* Outlet stream temperature;
880* Overall Heat Transfer Coefficient;
881* Exchange Surface Area;
882* Pressure drop.
883
884== SET ==
885* Phase of the fluid entering the heater;
886* Number of stream components(Ncomp/NcompS);
887* Basic composition (mass or molar).
888";
889       
890        VARIABLES
891in      Inlet_q as heat_stream  (Brief = "Inlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{in}");
892
893        EQUATIONS
894        "Equate Heat Stream"
895        Inlet_q.Q = Q;
896       
897        "General equation of heat exchange"
898        Inlet_q.Q = U * A * lmtd;
899       
900end
901
902Model cooler_mDL as heatex_mD
903        ATTRIBUTES
904        Pallete         = true;
905        Icon            = "icon/heaterR";
906        Brief           = "Model of a Heat Exchanger";
907        Info =
908"== GENERAL ==
909*Model of a simplified heat exchanger;
910*The model performes mass and energy balances only.
911
912== ASSUMPTIONS ==
913* Steady-State operation;
914* No heat loss to the surroundings.
915
916== SPECIFY ==
917* The inlet streams:
918  flow rate
919  temperature
920  pressure
921  stream composition;
922* Outlet vapour fraction;
923* Overall Heat Transfer Coefficient;
924* Exchange Surface Area;
925* Pressure drop.
926
927== SET ==
928* Phase of the fluid entering the heater;
929* Number of stream components(Ncomp/NcompS);
930* Basic composition (mass or molar).
931";
932       
933        VARIABLES
934out     Outlet_q        as heat_stream  (Brief = "Outlet Heat Stream", PosX=0.0, PosY=0.5, Symbol="_{out}");
935
936        EQUATIONS
937        "Equate Heat Stream"
938        Outlet_q.Q = -Q;
939       
940        "General equation of heat exchange"
941        Outlet_q.Q = U * A * lmtd;
942       
943end
944
945FlowSheet teste_heater_m
946       
947#*-------------------------------------------------------------------
948* Declaracao de dispositivos (ou blocos contendo o modelo)
949*--------------------------------------------------------------------*#
950       
951        DEVICES
952        SHot as main_sourceR;
953        H as heater_mLD;
954        SQ as heat_sourceR;
955       
956#*-------------------------------------------------------------------
957* Especifica as conexoes entre os modelos
958*--------------------------------------------------------------------*#
959       
960        CONNECTIONS
961        SHot.Outlet to H.Inlet;
962        SQ.Outlet_q to H.Inlet_q;
963       
964#*-------------------------------------------------------------------
965* Especifica variaveis definidas no modelo
966*--------------------------------------------------------------------*#
967       
968        SPECIFY
969
970        SHot.Fluid.Fw = 52200 * 'kg/h';
971        SHot.Solid.Fw = 48800 * 'kg/h';
972        SHot.T = 340 * 'K';
973        SHot.P = 1 * 'atm';
974        SHot.CompositionOfSolid = [0.5110, 0.0657, 0.2190, 0.0001, 0, 0, 0, 0.0004, 0];
975        SHot.CompositionOfFluid(1:11) = [0.982744, 0.000223825, 0.000463338, 0.00172757, 0, 0.00068594, 0.000144428, 0, 0, 0.000383508, 0.0130415];
976        SHot.CompositionOfFluid(12:NComp) = 0;
977        H.Pdrop = 0 * 'atm';
978        H.Outlet.T = 350.0 * 'K';
979        H.U = 0.69445 * 'kW/m^2/K';
980        H.lmtd = 10 * 'K';
981
982
983
984#*-------------------------------------------------------------------
985#Parametros
986*--------------------------------------------------------------------*#
987         
988        PARAMETERS
989        PP as Plugin    (Brief = "External Physical Properties",
990                Type="PP",
991                Project = "../Flowsheets/v2_2/Fluid_v2_2.vrtherm"
992        );
993        PPS as Plugin   (Brief = "External Physical Properties",
994                Type="PP",
995                Project = "../Flowsheets/v2_2/Solid_v2_2.vrtherm"
996        );
997       
998        NComp   as Integer      (Brief = "Number of chemical components in the fluid phase");
999        NCompS  as Integer      (Brief = "Number of chemical components in the solid phase");
1000       
1001#*-------------------------------------------------------------------
1002* Define o valor dos parametros declarados no modelo
1003*--------------------------------------------------------------------*#
1004       
1005        SET
1006        NComp = PP.NumberOfComponents();
1007        NCompS = PPS.NumberOfComponents();
1008        SHot.ValidPhases = "Liquid-Only";
1009        SHot.CompositionBasis = "Molar";
1010#*-------------------------------------------------------------------
1011* Opcoes de Solver
1012*--------------------------------------------------------------------*#
1013       
1014        OPTIONS
1015        Dynamic = false;
1016        #GuessFile = "teste3.rlt";
1017       
1018end
1019
1020FlowSheet teste_cooler_m
1021       
1022#*-------------------------------------------------------------------
1023* Declaracao de dispositivos (ou blocos contendo o modelo)
1024*--------------------------------------------------------------------*#
1025       
1026        DEVICES
1027        Scold as main_sourceR;
1028        H as cooler_mLD;
1029       
1030       
1031#*-------------------------------------------------------------------
1032* Especifica as conexoes entre os modelos
1033*--------------------------------------------------------------------*#
1034       
1035        CONNECTIONS
1036        Scold.Outlet to H.Inlet;
1037       
1038#*-------------------------------------------------------------------
1039* Especifica variaveis definidas no modelo
1040*--------------------------------------------------------------------*#
1041       
1042        SPECIFY
1043
1044        Scold.Fluid.Fw = 52200 * 'kg/h';
1045        Scold.Solid.Fw = 4800 * 'kg/h';
1046        Scold.T = 345 * 'K';
1047        Scold.P = 1 * 'atm';
1048        Scold.CompositionOfSolid = [0.5110, 0.0657, 0.2190, 0.0001, 0, 0, 0, 0.0004, 0];
1049        Scold.CompositionOfFluid(1:11) = [1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0];
1050        Scold.CompositionOfFluid(12:NComp) = 0;
1051        H.Pdrop = 0 * 'atm';
1052        H.Outlet.T = 335 * 'K';
1053        H.U = 0.69445 * 'kW/m^2/K';
1054        H.lmtd = 10 * 'K';
1055
1056
1057
1058#*-------------------------------------------------------------------
1059#Parametros
1060*--------------------------------------------------------------------*#
1061         
1062        PARAMETERS
1063        PP as Plugin    (Brief = "External Physical Properties",
1064                Type="PP",
1065                Project = "../Flowsheets/v2_2/Fluid_v2_2.vrtherm"
1066        );
1067        PPS as Plugin   (Brief = "External Physical Properties",
1068                Type="PP",
1069                Project = "../Flowsheets/v2_2/Solid_v2_2.vrtherm"
1070        );
1071       
1072        NComp   as Integer      (Brief = "Number of chemical components in the fluid phase");
1073        NCompS  as Integer      (Brief = "Number of chemical components in the solid phase");
1074       
1075#*-------------------------------------------------------------------
1076* Define o valor dos parametros declarados no modelo
1077*--------------------------------------------------------------------*#
1078       
1079       
1080        SET
1081        NComp = PP.NumberOfComponents();
1082        NCompS = PPS.NumberOfComponents();
1083        Scold.ValidPhases = "Liquid-Only";
1084        Scold.CompositionBasis = "Molar";
1085       
1086#*-------------------------------------------------------------------
1087* Opcoes de Solver
1088*--------------------------------------------------------------------*#
1089       
1090        OPTIONS
1091        Dynamic = false;
1092
1093       
1094end
1095
Note: See TracBrowser for help on using the repository browser.