Last change
on this file since 1 was
1,
checked in by Rafael de Pelegrini Soares, 16 years ago

Initial import of the library


Property svn:eolstyle 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) 20022004 the author 

6  * 

7  * Author: Rafael de Pelegrini Soares 

8  * $Id: sample_cstr.mso 1 20060620 17:33:53Z rafael $ 

9  **# 

10  

11  using "reactors/cstr"; 

12  

13  FlowSheet 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*(TTw); #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*(TwT)FA0*Cps*(TT0)+HR*rA*V)/(CA0*V*Cps)); 

64  

65  "Bal. Material Global" 

66  diff(V)=FinFout; 

67  

68  "Bal. Material A" 

69  diff(V*CA)=(FA0Fout*CA+V*rA); 

70  

71  "Bal. Material B" 

72  diff(V*CB)=(FB0Fout*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]; 

99  outputLevel = "high"; 

100  end 

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