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

Initial import of the library


Property svn:eolstyle set to
native

Property svn:keywords set to
Id

File size:
2.0 KB

Rev  Line  

[1]  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.