# Ticket #124: ColunaDissertação.mso

File ColunaDissertação.mso, 2.7 KB (added by , 13 years ago) |
---|

Line | |
---|---|

1 | |

2 | using "stage_separators/column"; |

3 | |

4 | # column section with 8 trays |

5 | FlowSheet Distillation_kettle_cond_Test |

6 | PARAMETERS |

7 | PP as Plugin(Brief="External Physical Properties", |

8 | Type="PP", |

9 | Components = [ "o-xylene", "m-xylene", "p-xylene","n-propylbenzene"], |

10 | LiquidModel = "PR", |

11 | VapourModel = "PR" |

12 | ); |

13 | NComp as Integer; |

14 | |

15 | VARIABLES |

16 | Qc as energy_source (Brief="Heat rate removed from condenser"); |

17 | Qr as energy_source (Brief="Heat rate supplied to reboiler"); |

18 | |

19 | SET |

20 | NComp = PP.NumberOfComponents; |

21 | |

22 | DEVICES |

23 | col as Distillation_kettle_cond; |

24 | feed as source; |

25 | zero as stream; |

26 | |

27 | CONNECTIONS |

28 | feed.Outlet to col.trays(65).Inlet; |

29 | zero to col.reb.Inlet; |

30 | zero to col.trays([1:64]).Inlet; |

31 | zero to col.trays([66:col.NTrays]).Inlet; |

32 | Qc.OutletQ to col.cond.InletQ; |

33 | Qr.OutletQ to col.reb.InletQ; |

34 | |

35 | SPECIFY |

36 | feed.Fw = 71500 * 'kg/h'; #fornecida pelo artigo de Fabio, 71500 kg/h - deve ser colocada em mol/s - Calcular novo valor |

37 | feed.T = (172+273.15) * 'K'; #fornecida pelo artigo de Fabio |

38 | feed.P = 2.4 * 'kgf/cm^2'; #fornecida pelo artigo de Fabio |

39 | feed.Composition =[0.242, 0.170, 0.585, 0.003]; #estimada por mim com base em informações do artigo |

40 | |

41 | |

42 | zero.F = 0 * 'kmol/h'; #corrente necessária para conexão |

43 | zero.T = 300 * 'K'; |

44 | zero.P = 1 * 'atm'; |

45 | zero.z = 1/NComp; |

46 | zero.v = 0; |

47 | zero.h = 0 * 'J/mol'; |

48 | |

49 | #INSERIR MAIS DADOS DA DISSERTAÇÃO |

50 | |

51 | col.sptop.Outlet2.F = 280000 * 'kmol/h'; |

52 | col.reb.OutletL.F = 6000 * 'kmol/h'; |

53 | col.sptop.frac = 0.2; |

54 | col.cond.OutletV.F = 0 * 'kmol/h'; |

55 | Qr.OutletQ.Q = 3.7743e6 * 'kJ/h'; |

56 | Qc.OutletQ.Q = -3.71e6 * 'kJ/h'; |

57 | col.pump1.dP = 16 * 'kPa'; |

58 | col.trays.Emv = 1; |

59 | #col.cond.OutletL.z(1)=0.98; #fornecida pelo artigo |

60 | |

61 | col.alfaTopo = 2; |

62 | |

63 | SET |

64 | col.NTrays = 130; |

65 | col.cond.V = 2 * 'm^3'; |

66 | col.cond.Across = 1 * 'm^2'; |

67 | col.trays.V = 4 * 'ft^3'; |

68 | col.trays.Ah = 0.394 * 'ft^2'; |

69 | col.trays.lw = 20.94 * 'in'; |

70 | col.trays.hw = 0.125 * 'ft'; |

71 | col.trays.Q = 0 * 'kW'; |

72 | col.trays.beta = 0.8; |

73 | col.trays.alfa = 4; |

74 | col.trays.Ap = 3.94 * 'ft^2'; |

75 | col.reb.V = 2 * 'm^3'; |

76 | col.reb.Across = 1 * 'm^2'; |

77 | |

78 | INITIAL |

79 | # condenser |

80 | col.cond.OutletL.T = 260 *'K'; |

81 | col.cond.Level = 1 * 'm'; |

82 | col.cond.OutletL.z([1:3]) = [0.98, 0.01, 0.01]; #0.98 do artigo |

83 | |

84 | # reboiler |

85 | col.reb.OutletL.T = 350 *'K'; |

86 | col.reb.Level = 1 * 'm'; |

87 | col.reb.OutletL.z([1:3]) = [0.1, 0.7, 0.15]; |

88 | |

89 | # column trays |

90 | col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K'; |

91 | col.trays.Level = 1.2 * col.trays.hw; |

92 | col.trays.OutletL.z([1:3]) = [0.5, 0.05, 0.01]; |

93 | |

94 | OPTIONS |

95 | TimeStep = 0.1; |

96 | TimeEnd = 2; |

97 | NLASolver(File="nlasolver"); |

98 | #GuessFile="Distillation_kettle_cond_Test.rlt"; |

99 | #Dynamic = false; |

100 | end |