source: mso/sample/miscellaneous/sample_cstr_simple.mso @ 51

Last change on this file since 51 was 32, checked in by Argimiro Resende Secchi, 16 years ago

clean up.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 2.0 KB
Line 
1#*---------------------------------------------------------------------
2* This file is property of the author and cannot be used, copyed
3* or modified without permission.
4*
5* Copyright (C) 2002-2004  the author
6*----------------------------------------------------------------------
7* Author: Rafael de Pelegrini Soares
8* $Id: sample_cstr_simple.mso 32 2006-09-23 00:16:00Z arge $
9*---------------------------------------------------------------------*#
10
11using "reactors/cstr";
12
13FlowSheet CSTR
14
15        PARAMETERS
16        Fin as flow_vol;
17        Fout as flow_vol;
18        CA0 as conc_mol;
19        CB0 as conc_mol;
20        FA0 as flow_mol;
21        FB0 as flow_mol;
22        # Troca térmica
23        Tw as temperature;
24        U as heat_trans_coeff;
25        A as area;
26        # Reação
27        k0 as frequency;
28        Ea as energy_mol (Brief = "Energia de Ativação");
29        R  as energy_mol;
30       
31        VARIABLES
32        q as heat_rate;
33        HR as energy_mol;
34        T as temperature;
35        k as frequency;
36        CA as conc_mol;
37        CB as conc_mol;
38        V as volume;
39        rA as reaction_mol;
40        rB as reaction_mol;
41        Cps as cp_mol; # (Unit="cal/mol/K");
42        T0 as temperature;
43       
44        EQUATIONS
45       
46        "Troca térmica"
47        q=U*A*(T-Tw); #j/s
48       
49        "Constante de Reação"
50        k=k0*exp(Ea/R*(1/436.15 - 1*"K"/T));
51       
52        "Taxas de Reação"
53        rA=-k*CA; # (mol/m^3/h)  >> reação de primeira ordem <<
54        rB=-rA;
55       
56        "Calor de Reação"
57        HR=-(80*250)*"cal/mol";
58       
59        "Cp da mistura"
60        Cps=(0.5*250)*3;# *"cal/mol/K";
61       
62        "Balanço de Energia"
63        diff(T)=((U*A*(Tw-T)-FA0*Cps*(T-T0)+HR*rA*V)/(CA0*V*Cps));
64       
65        "Bal. Material Global"
66        diff(V)=Fin-Fout;
67       
68        "Bal. Material A"
69        diff(V*CA)=(FA0-Fout*CA+V*rA);
70       
71        "Bal. Material B"
72        diff(V*CB)=(FB0-Fout*CB+V*rB);
73       
74        "Temperatura da corrente de alimentação"
75        T0= (40 + 273.15 + time*200/1000/"s") * "K";
76       
77        SET
78        Fin = 0.125997903*"m^3/h";
79        Fout = 0.125997903*"m^3/h";
80        CA0=9.3652352*"mol/m^3";
81        CB0=0*"mol/m^3";
82        FA0 = Fin*CA0;
83        FB0 = Fin*CB0;
84        Tw = (140+273)*"K";
85        k0=0.8*"1/h"; #constante para T=436.15 K
86        Ea=28960*"cal/mol";
87        R=1.98*"cal/mol/K";
88        U = 0 *"cal/m^2/K/s";#isolado
89        A = 1*"m^2";
90       
91        INITIAL
92        T = (140 + 273) * "K";
93        V = 0.9*6*"m^3";
94        CA = 2*"mol/m^3";
95        CB = 6*"mol/m^3";
96       
97        OPTIONS
98        time = [0:10:1200];
99end
Note: See TracBrowser for help on using the repository browser.