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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
package BTX
model ms"material stream"
extends Simulator.Streams.Material_Stream;
extends Simulator.Files.Thermodynamic_Packages.NRTL;
end ms;
model Condensor
extends Simulator.Unit_Operations.Distillation_Column.Cond;
extends Simulator.Files.Thermodynamic_Packages.NRTL;
end Condensor;
model Tray
extends Simulator.Unit_Operations.Distillation_Column.DistTray;
extends Simulator.Files.Thermodynamic_Packages.NRTL;
end Tray;
model Reboiler
extends Simulator.Unit_Operations.Distillation_Column.Reb;
extends Simulator.Files.Thermodynamic_Packages.NRTL;
end Reboiler;
model DistColumn
extends Simulator.Unit_Operations.Distillation_Column.DistCol;
Condensor condensor(NOC = NOC, comp = comp, condType = condType,boolFeed = boolFeed[1], T(start = 400));
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 = 10), each vapMolFlo(each start = 10));
end DistColumn;
model flowsheet
import feed_data = Simulator.Files.Chemsep_Database;
parameter feed_data.Benzene benzene;
parameter feed_data.Toluene toluene;
parameter feed_data.Pxylene xylene;
parameter Integer NOC = 3;
parameter feed_data.General_Properties comp[NOC] = {benzene, toluene, xylene};
BTX.ms feed(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {-76, 12}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
BTX.DistColumn distColumn1(NOC = NOC, comp = comp, feedStages = {10}, noOfFeeds = 1, noOfStages = 20) annotation(
Placement(visible = true, transformation(origin = {-26, 14}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
BTX.ms ben(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {26, 62}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
BTX.ms bott(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {42, -22}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
Simulator.Streams.Energy_Stream energy_Stream1 annotation(
Placement(visible = true, transformation(origin = {-27, 35}, extent = {{-5, -5}, {5, 5}}, rotation = 0)));
Simulator.Streams.Energy_Stream energy_Stream2 annotation(
Placement(visible = true, transformation(origin = {-25, -7}, extent = {{-5, -5}, {5, 5}}, rotation = 0)));
BTX.DistColumn distColumn2(NOC = NOC, comp = comp, feedStages = {15}, noOfStages = 30) annotation(
Placement(visible = true, transformation(origin = {102, 18}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
BTX.ms tol(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {160, 50}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
ms xyl(NOC = NOC, comp = comp) annotation(
Placement(visible = true, transformation(origin = {158, -22}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
Simulator.Streams.Energy_Stream energy_Stream3 annotation(
Placement(visible = true, transformation(origin = {102, -2.22045e-16}, extent = {{-4, -4}, {4, 4}}, rotation = 0)));
Simulator.Streams.Energy_Stream energy_Stream4 annotation(
Placement(visible = true, transformation(origin = {102, 36}, extent = {{-4, -4}, {4, 4}}, rotation = 0)));
equation
connect(distColumn1.distillate, ben.inlet) annotation(
Line(points = {{-16, 22}, {-10, 22}, {-10, 62}, {16, 62}}));
connect(energy_Stream1.outlet, distColumn1.condensor_duty) annotation(
Line(points = {{-22, 35}, {-22, 29.5}, {-20, 29.5}, {-20, 24}}));
connect(distColumn2.reboiler_duty, energy_Stream3.outlet) annotation(
Line(points = {{110, 8}, {110, 3.5}, {106, 3.5}, {106, 0}}));
connect(distColumn2.condensor_duty, energy_Stream4.outlet) annotation(
Line(points = {{108, 28}, {108, 32}, {106, 32}, {106, 36}}));
connect(distColumn1.reboiler_duty, energy_Stream2.outlet) annotation(
Line(points = {{-18, 4}, {-18, -1.5}, {-20, -1.5}, {-20, -7}}));
connect(distColumn2.distillate, tol.inlet) annotation(
Line(points = {{112, 26}, {138, 26}, {138, 50}, {150, 50}}));
connect(distColumn2.bottoms, xyl.inlet) annotation(
Line(points = {{112, 10}, {138, 10}, {138, -22}, {148, -22}, {148, -22}}));
connect(bott.outlet, distColumn2.feed[1]) annotation(
Line(points = {{52, -22}, {68, -22}, {68, 18}, {92, 18}, {92, 18}}));
connect(distColumn1.bottoms, bott.inlet) annotation(
Line(points = {{-16, 6}, {-10, 6}, {-10, -22}, {32, -22}}));
connect(feed.outlet, distColumn1.feed[1]) annotation(
Line(points = {{-66, 12}, {-36, 12}, {-36, 14}, {-36, 14}}));
feed.P = 101325;
feed.T = 298.15;
feed.compMolFrac[1, :] = {0.333, 0.333, 0.333};
feed.totMolFlo[1] = 10.049;
distColumn1.condensor.P = 101325;
distColumn1.reboiler.P = 101325;
distColumn1.refluxRatio = 2;
bott.totMolFlo[1] = 6.7;
distColumn2.condensor.P = 101325;
distColumn2.reboiler.P = 101325;
distColumn2.refluxRatio = 3;
xyl.totMolFlo[1] = 3.349;
annotation(
Diagram(coordinateSystem(extent = {{-100, -100}, {200, 100}}, initialScale = 0.1), graphics = {Text(origin = {14, 106}, extent = {{-44, 20}, {96, -50}}, textString = "BTX separation columns", fontName = "MV Boli"), Text(origin = {125, 76}, extent = {{-69, 14}, {-1, -4}}, textString = "-MehulKumar Sutariya", fontName = "Segoe Script"), Text(origin = {63, -93}, extent = {{-91, -1}, {87, -7}}, textString = "Total time = 1500 seconds"), Text(origin = {19, 55}, extent = {{-7, 5}, {25, -17}}, textString = "Benzene"), Text(origin = {157, 36}, extent = {{-15, 8}, {19, -10}}, textString = "Toluene"), Text(origin = {166, -43}, extent = {{-18, 9}, {8, -1}}, textString = "Xylene"), Text(origin = {-78, -9}, extent = {{-10, 7}, {10, -7}}, textString = "FEED")}),
Icon(coordinateSystem(extent = {{-100, -100}, {200, 100}})),
__OpenModelica_commandLineOptions = "");
end flowsheet;
end BTX;
|
