1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
package process
model Inlet
extends Simulator.Streams.Material_Stream;
extends Simulator.Files.Thermodynamic_Packages.Raoults_Law;
end Inlet;
model Condensor
extends Simulator.Unit_Operations.Distillation_Column.Cond;
extends Simulator.Files.Thermodynamic_Packages.Raoults_Law;
end Condensor;
model Tray
extends Simulator.Unit_Operations.Distillation_Column.DistTray;
extends Simulator.Files.Thermodynamic_Packages.Raoults_Law;
end Tray;
model Reboiler
extends Simulator.Unit_Operations.Distillation_Column.Reb;
extends Simulator.Files.Thermodynamic_Packages.Raoults_Law;
end Reboiler;
model DistColumn
extends Simulator.Unit_Operations.Distillation_Column.DistCol;
Condensor condensor(NOC = NOC, comp = comp, condType = condType, boolFeed = boolFeed[1], T(start = 300));
Reboiler reboiler(NOC = NOC, comp = comp, boolFeed = boolFeed[noOfStages]);
Tray tray[noOfStages - 2](each NOC = NOC, each comp = comp, boolFeed = boolFeed[2:noOfStages -1], each liqMolFlo(each start = 150), each vapMolFlo(each start = 150));
end DistColumn;
model main
import data =Simulator.Files.Chemsep_Database;
parameter data.Methanol meth;
parameter data.Water wat;
parameter Integer NOC=2;
parameter data.General_Properties comp[NOC]= {meth,wat};
process.Inlet feed(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {-68, 46}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
process.DistColumn DC(NOC = NOC, comp = comp, condType = "Total", feedStages = {5}, noOfFeeds = 1, noOfHeatLoads = 0, noOfSideDraws = 0, noOfStages = 10) annotation(
Placement(visible = true, transformation(origin = {46, 52}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
process.Inlet top(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {82, 64}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
process.Inlet bottom(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {84, -18}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
Simulator.Streams.Energy_Stream CD annotation(
Placement(visible = true, transformation(origin = {48, 82}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
Simulator.Streams.Energy_Stream RD annotation(
Placement(visible = true, transformation(origin = {52, 8}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
Simulator.Unit_Operations.Rigorous_HX HE(Calculation_Method = "Hot_Fluid_Outlet_Temperature", Heat_Loss = 0, NOC = NOC, comp = comp, deltap_cold = 0, deltap_hot = 0) annotation(
Placement(visible = true, transformation(origin = {-14, -18}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
process.Inlet coldfeed(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {-66, -28}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
process.Inlet residue(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {28, -74}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
equation
connect(bottom.outlet, HE.Hot_In) annotation(
Line(points = {{94, -18}, {98, -18}, {98, -42}, {-42, -42}, {-42, -14}, {-22, -14}, {-22, -14}}));
connect(feed.inlet, HE.Cold_Out) annotation(
Line(points = {{-78, 46}, {-90, 46}, {-90, -2}, {22, -2}, {22, -24}, {-6, -24}, {-6, -20}}));
connect(HE.Hot_Out, residue.inlet) annotation(
Line(points = {{-6, -14}, {4, -14}, {4, -74}, {18, -74}}));
connect(DC.reboiler_duty, RD.inlet) annotation(
Line(points = {{53, 42}, {53, 8}, {42, 8}}));
connect(DC.condensor_duty, CD.inlet) annotation(
Line(points = {{52, 62}, {28, 62}, {28, 82}, {38, 82}}));
connect(DC.distillate, top.inlet) annotation(
Line(points = {{56, 59}, {58, 59}, {58, 64}, {72, 64}}));
connect(DC.bottoms, bottom.inlet) annotation(
Line(points = {{56, 45}, {68, 45}, {68, -18}, {74, -18}}));
connect(feed.outlet, DC.feed[1]) annotation(
Line(points = {{-58, 46}, {-11, 46}, {-11, 52}, {36, 52}}));
connect(coldfeed.outlet, HE.Cold_In) annotation(
Line(points = {{-56, -28}, {-28, -28}, {-28, -20}, {-22, -20}, {-22, -20}}));
feed.T = 324.84;
coldfeed.P = 101325;
coldfeed.compMolFrac[1, :] = {0.35, 0.65};
coldfeed.totMolFlo[1] = 60.07;
DC.condensor.P = 101325;
DC.reboiler.P = 101325;
DC.refluxRatio = 2;
bottom.totMolFlo[1] = 36.74;
residue.T = 327.84;
HE.U = 1000;
end main;
end process;
|
