Design of Batch Distillation Columns Using Short-Cut Method at Constant Reflux

A short-cutmethod for batch distillation columnsworking at constant refluxwas applied to solve a problemof four components that needed to be separated and purified to a mole fraction of 0.97 or better. Distillation columns with 10, 20, 30, 40, and 50 theoretical stages were used; reflux ratio was varied between 2 and 20.Three quality indexes were used and compared: Luyben’s capacity factor, total annual cost, and annual profit. The best combinations of theoretical stages and reflux ratio were obtained for each method. It was found that the best combinations always required reflux ratios close to theminimum. Overall, annual profit was the best quality index, while the best combination was a distillation column with 30 stages, and reflux ratio’s of 2.0 for separation of benzene (i), 5.0 for the separation of toluene (ii), and 20 for the separation of ethylbenzene (iii) and purification of o-xylene (iv).


Introduction
Distillation is the most commonly used separation process in the chemistry and petrochemical industry, mostly employed in continuous processes but also used in discontinuous processes.Batch distillation is widely used for the separation of specialty and fine chemicals and for the recovery of small quantities of solvent during the production of high purity and added value products.Batch processing is the main feature of the pharmaceutical, biochemical, and specialty chemical industries.
One advantage of batch distillation is the possibility to separate a mixture of several components with only one column and, although a batch distillation column uses more energy than a continuous column, many times it is convenient to operate in batch mode.For these reasons the design, analysis, and optimization of batch distillation columns have obtained attention by several researchers.
Figure 1 shows a conventional batch distillation column consisting of a reboiler or pot, a column with plates or packing, a condenser, a reflux tank, accumulative recipients for products (main cuts), and byproducts (off cuts).
In batch distillation columns the feed may be charged into the reboiler at the beginning of the operation and then heat is provided at the reboiler to evaporate part of the liquid to generate a vapor that rises through the column until it reaches the condenser where it is converted to liquid and collected at the reflux tank.From this reflux tank a part is returned to the column as liquid reflux that descends through the column contacting with the vapor counter-currently if the column is full of packing, or contacting the vapor in a crosscurrent pattern if the column contains plates.
Mass transfer occurs when vapor and liquid interact and as a result the light components rise with the vapor and the heavy compounds flow within the liquid.In this way a profile of concentration, temperature, and pressure is drawn throughout the column.At the bottom, pressure and temperature are higher because the heavy compounds have higher boiling points, and because the vapor has a loss of pressure when passing through the packing or plates.The more volatile components concentrate at the top of the column and the less volatile or heavy compounds concentrate at the bottom.These steps at the beginning of the operation are required in order to put the column at total reflux and this operation is used in continuous and batch processes.In normal operation of a batch process, liquid distillate is withdrawn from the reflux tank, taking away the more volatile compounds first and increasing at the same time the concentration of the less volatile or heavy components at the bottom.
The main characteristic of a batch distillation column is that concentrations and temperatures are changing with time at any part of the column.
There are two main operating methods for batch distillation columns: constant reflux (with variable product composition at distillate) and variable reflux (with constant product composition, for one component, at distillated product).Another operating method is the optimal reflux policy, which tries to get an economical combination of the two main operating methods, according to established criteria like: minimum time, maximum product, or maximum profit.
This dynamic process may be modeled by mass and energy balances and coupling them with equilibrium ratios obtained from thermodynamic or experimental vapor-liquid equilibrium data.
The models may vary from McCabe Thiele binary distillation with constant or variable reflux, to short-cut or approximate to rigorous method.Each one has advantages and disadvantages.In this work we use a constant reflux short-cut method to design a batch distillation system (column, condenser, and reboiler) and test three quality indexes: Luyven's capacity factor, total annual cost, and annual profit.

Materials and Methods
2.1.The Problem.Design (size) a distillation column with condenser and reboiler operating at atmospheric pressure with constant reflux, to separate and purify an equimolar mixture of benzene (i), toluene (ii), ethylbenzene (iii), and oxylene (iv), containing 100 kmol of each, at molar fractions of 0.97 for each component.For vapor liquid equilibrium compare reported data for binary pairs, predicted for multicomponent mixture, versus predicted by Raoult law and ideal gases.
If we receive 400 kmol of the mixture per day, then we would need to size the distillation column filled with sieve trays; we would also need to calculate the area of the condenser and reboiler, as well as the steam and cooling water needed in the condenser.In addition, report Luyben's capacity factor, total annual cost, and annual profit for several alternatives for number of stages, for example, 10, 20, 30, 40, and 50 theoretical plates, and for several values of reflux ratios.

Previous Work and Model Description. It was probably
Diwekar [1] in her 1988 Ph.D. dissertation who first reported the short-cut method for batch distillation columns, but we think that Diwekar and Madhavan paper [2] is considered the first in a known journal.Their paper described both constant and variable reflux.In 1993, Sundaram and Evans [3] reported their method for simulating batch distillation operation.It only covered constant reflux ratio.Both papers assumed that, at any time, the batch column was identical to the rectifying section of a continuous column and applied the short-cut (FUG) method: Fenske [4], Underwood [5,6], and Gilliland [7] assuming that the relative volatility of the components remained constant.Both methods began with setting Barolo and Guarice [8] extended their continuous shortcut method to batch distillation columns and used an infinite stage column.
Seader and Henley [9,10] applied Sundaram and Evans [3] method and showed how the FUG equations were used to solve a simulation problem with constant reflux.They presented global (1) and component (i) (2) mass balances: Solving (1) for  (using  for time increment index) and for time (using 0 for initial): Solving (2) for   i (using  for time increment index) and (i) for any component: Equations ( 3) and ( 5) form the basis of the calculations.Equation (3) gets a new value for total molar mass at the new time at the bottom of the column, and (5) gets the concentration of the component (i) at the bottom at the new time.
Luyben [11] did not use the short-cut method but studied and reported simulation results for binary separations covering the following practical aspects of optimal design: number of trays, reflux ratio, tray holdup, and amount of initial charge to the still.For a complete batch cycle he included startup time or equilibration time (  ), time needed to bring the column to steady state condition on total reflux; overhead product withdrawal or process time (  ), period during which on-specification distillate product is produced; final bottoms product purification time (  ), time required to remove enough low boiler from the still pot and column trays so that the bottoms product will be on-specification when the liquid in the column drains down into the still pot.He established that the design objective was to maximize capacity and defined capacity as the amount of total onspecification products produced per unit time (i.e., the sum of distillate and bottoms).He concluded that, as is usually the case in continuous distillation columns, the optimum combination of reflux ratio and number of trays is usually in the direction of several trays and low reflux ratio.
Al-Tuwaim and Luyben [12] used simulation results of a rigorous model for binary and ternary mixtures with the following assumptions: theoretical trays, equimolar overflows, constant relative volatilities, constant tray holdup (1 lbmol), constant reflux drum holdup (10 lbmol), and total slop recycle.Their approach was to find the optimal reflux ratios for columns with various numbers of trays and specify the capacity factor proposed by Luyben [11], which applied for a multicomponent mixture, resulting in the following: For a given separation we need to specify relative volatility of the system, product purities, energy cost, and material of construction of the column.The authors also used the total annual cost that involves capital cost (column, condenser, and reboiler) and operation cost (steam and cooling water) for a one-year operation period.In this work, to solve the problem, the strategy was to simulate the operation of a batch distillation column with 10 to 50 theoretical stages and reflux ratios of 2, 5, 10, 15 and 20, and rank the runs according with three qualification parameters: Luyben's capacity factor, total annual cost, and annual profit.It was expected that the qualification parameters may point to the most economical combinations of tray and reflux ratio.
From the papers of Diwekar and Madhavan [2] and Sundaram and Evans [3] for constant reflux ratio, we found that both methods were similar and differed only in small details.Seader and Henley [9,10] used the sequence of equations shown in Table 2 that apply FUG method.Gilliland in 1940 [7] presented his results about the relationship between the number of minimum stages and the number of theoretical stages with minimum reflux ratio and operational reflux ratio in graphical form.In 1972, Gilliland's results were correlated by Molokanov et al. [13] and in 1975 by Eduljee [14] with the simpler expression of (g) in Table 2.The book of Mujtaba [15] also covers the short-cut method for batch distillation.

Equilibrium Ratios and Relative
Volatility Used for the Short-Cut Method.For liquid vapor equilibrium we used the binary experimental data reported in Gmehling et al. [16].We used the reported  1 / 1 to obtain  1 and (1 −  1 )/(1 −  1 ) to obtain  2 for each point.From these values we calculated the relative volatility  =  1 / 2 for each reported point.Higher values of  1 ,  2 , and volatilities were obtained at the bottom of the distillation column where the temperature was higher.Table 3 shows the values obtained for the four components in the three separations.It is observed that the most difficult separation was ethylbenzene/o-xylene ( = 1.04), while the simpler separation was benzene/toluene ( = 1.34).
In order to get relative volatilities for the mixture, we ran simulations of continuous distillation columns using Aspen plus with Peng Robinson and NRTL models that have low deviation (0.10 averaged) between experimental and predicted binary data.Table 3 also shows the values predicted by (7) that represents Raoult's Law, at the beginning of the separation at 80 and 110 ∘ C at the top and bottom of the column, respectively, The vapor pressure of each component was predicted by Antoine's equation (8).Constants , , , , and  were taken from Table 23-IV from Yaws [17] log It is observed that all thermodynamic models, such as activity coefficient model (NRTL), equation of state (Peng Robinson), and Raoult Law, predict the same ranking in terms of ease of separation for benzene/toluene and the most difficult for ETB/o-xylene.Also, it is observed that the prediction of ideal gas and ideal liquid behavior predicted by Raoult law is very good, probably because the four components have similar structures as they form part of the aromatic family.From the results of the simulation with Aspen plus we also obtained physical properties at several points of the column.We decided to take the bottom of the column as the reference for physical properties and for relative volatilities to be used in the short-cut method.When averaged values for relative volatilities were used, the results had the same tendency and the separations improved because  i-ii increased from 6.33 to 7.16. Hengstebek-Geddes: Stepping forward in time and calculating new molar mass at bottom and new molar fraction of all components With the new value of  min , estimate the molar fraction of the reference component at distillate With the new value of  min , calculate molar fraction of nonreference component at distillate Repeat steps 3 to 6 for several cycles.

S t o p
Computer programs in Excel, Mathcad, and Visual Basic were created with the equations of Table 2 and were applied to simulate the batch columns that fixed the value of the boilup vapor at 100 kmol/hr and varied the number of theoretical plates (10,20,30,40,50) and external reflux ratio (/ = 2, 5, 10, 15, 20).They were initially applied for the separation of component (i) (benzene), the lightest component with a purity of 0.97 in mass fraction.

Luyben Capacity Factor for Comparison.
Application of (6), neglecting equilibrium time (time needed to bring the column to steady state conditions on total reflux), and assigning 1 hour for feed and cleaning result in the following: The numerator shows the molar or mass amounts of products that satisfy the desired concentration, while the denominator is made up of three different parameters: the harvest time for each product ( i ,  ii ,  iii , and  iv ) that are collected on the main cut receivers, the time required to eliminate a volatile component ( i-ii ,  ii-iii , and  iii-iv ) whose mixtures are discharged on the off cut tanks, and the time needed for feeding and discharging the batch distillation column ( ⋅ , which equals 1 hour).It should be noted that  iii-iv is the time elapsed after the concentration of ethylbenzene drops below 0.97 and the time the concentration of xylene at bottom reaches 0.97.

Total Annual Cost for Comparison. Total annual cost
given by TAC = Annual cost of equipment + Annual cost of utilities (10) is an index used for comparison of alternatives.It takes into account the annual capital cost of equipment: distillation column (plates and shell), reboiler, and condenser, and also  annual operating cost (steam for reboiler and cooling water for condenser).With the boilup fixed at  = 100 kmol/h and external reflux ratio  = / specified, liquid flow rate ( =  * ) and internal ratio L/V may be calculated as With the molar flow rates of liquid and vapor, the diameter of the column may be calculated using physical properties (, , ) and geometric parameters to predict flooding velocity (  ) and fixing operational vapor velocity as a fraction () of it, the column diameter is With the molar flow of the vapor  that needed to be condensed and the heat of condensation of distillate vapor,   was calculated, and from it   was determined with the heat transfer equation.Additionally, the flow rate of cooling water was obtained At the bottom of the column heat was needed to vaporize  kmol/h of the liquid and with the heat of vaporization of the liquid at the bottom we obtained   , as well as the flow rate of the steam needed as main utility as follows: In order to identify the runs, the following nomenclature was used: N IJK, where: As an example, 20 ABC means a distillation column with 20 theoretical stages,  = 2 (A) for the first separation (four components),  = 5 (B) for the second separation (three components), and  = 10 (C) for the third separation (two components).

Annual Profit for Comparison.
A complete economical analysis would require the calculation of annual gain or profit: Annual income is the amount of money collected when the four harvested products are sold.Table 6.4 of Turton et al. [18] provides averaged values from the Chemical Market Reporter for January-April, 2001 see Table 1.
The selling price for ethylbenzene was assumed to be equal to that of ortho-xylene.

Four Components Separation. Figures 2(a) and 2(b) show
the results for a batch distillation column operating at constant reflux ratio of 2.0 with ten theoretical stages.In Figure 2(a) it is observed that at the bottom of the column, the composition of benzene decreased with time while the composition of the heavy components increased.In Figure 2(b) it is observed that the instant distillate composition of benzene was 1.0 at the beginning and then decreased to 0.86 at 2.2 hours, which was the last time reported, but the averaged composition of benzene remained above 0.97 all this time.We call this harvest time.It is observed that the larger changes in distillate composition occurred close to the final time.As observed in step four of Table 2, the computer program was needed to solve (g) in Table 2 for minimum reflux ratio with the Gilliland Eduljee correlation  min- and compare it with the result of (h) in Table 2 for minimum reflux ratio with the Underwood correlation  min- .Once (g) and (h) in Table 2 were solved, the results were used to find the number of minimum stages ( min ) that should be applied in other steps of the method.
In the first iterations we began with  min = /2 = 5, minimum reflux from Gilliland-Eduljee provided a value of  min- = 1.42 and  min- = 0.74.Using successive substitution we calculated the new values  min = 6.7956 and for this value we obtained  min- = 0.7400 and  min- = 0.7484 and a new value for  min = 6.7766.The last values for the first iteration were  min = 6.7766,  min- = 0.7483, and  min- = 0.7483.
In the process, time was increasing and  min decreased from 6.77 to 2.25 at the last time reported (2.2 hours), but minimum reflux ( min ) of Gilliland and Underwood was gradually increasing starting from about 0.75.Once the value of the Underwood minimum reflux or Guilliland-Eduljee equals or surpasses the value of the fixed operational reflux (2.0 in this case), the term (- min ) in Gilliland correlation turns negative and hence the program cannot continue with the calculations.Fortunately, the calculations performed before this point provided all the useful information we needed.
Figures 3 and 4 help us compare the use of different reflux ratios.Figures 3(a For reflux ratios of 5 and 10 and the same number of theoretical stages, more distillate product with    ≥ 0.97 was obtained, but also more time was needed.Since the beginning of this first distillation the harvested benzene product goes to main cut 1 tank of Figure 1 and at harvest time the distillate is deviated to the off cut 1 tank where a mixture of benzene, toluene and a little of ethylbenzene was kept. In order to know the time needed for all the benzene to disappear in the bottom of the column with  = 2.0, we extrapolated the data of  W-i versus time in Figure 2(a) obtaining first equation: And then solving for  W-i = 0.0, we find that about 3.0 hours are needed to eliminate all the benzene from the bottom.Nonetheless, it is likely that some benzene will be at the plates of packing of the column.Assuming that we were able to separate in the first distillation the 100 kmol of benzene, solving (4) for   = 300, we get  = 3.0 hours.The time needed for the elimination of benzene from the four-component original mixture is a rough approximation.When comparing a rigorous method versus a short-cut method, usually the first uses more time than the second.Also, when benzene is being depleting from the reboiler the other light components, such as toluene and ethylbenzene are being separated in the later stages of the distillation column.
Also in the approximation taken, we found that the composition of the three compounds that remained at the bottom (toluene, ethylbenzene, and o-xylene) had the same composition:  W-i ≈ 0.333.
Applying the same procedure for  = 5 and  = 10 provided  = 6 and 11 hours, respectively to eliminate the benzene from the bottom of the column ( W-i = 0).Also   the molar fraction of the other components tended to be equimolar at 0.333.
Figures 3(b) and 4(b) show that with the combination of low reflux  = 2 and low number of theoretical stages  = 10 distillate composition drops below 0.90, but for the combination of  = 40 and  = 2 the punctual composition for benzene in distillate is always close to unity for the three values of reflux ratio.Similar results were performed for  = 30 and 50.
Table 4 show harvest time and collected products for the different alternatives in this four-components separation.It is observed that with more reflux, better product (with  W-i ≥ 0.97) is collected but harvested and total time increase with respect to the lowest reflux of 2.0.It is possible to calculate a partial Luyben index (P 4 -CAP) for the four components separation by dividing the collected benzene product by the total time for the run.

Three Components Separation.
Beginning with  0 = 300 kmol we run the program first for  = 10 and  = 2 and observe that the first molar fraction of component ii (toluene) is 0.9995, the second molar fraction drops to 0.9538, and then the average composition is lower than the required 0.97.We see that this combination of stages and reflux is not enough to get the desired molar fraction of toluene distillate.This happens because the relative volatility for three components is lower (2.09) than the one for four components (2.38), and we may anticipate also difficulties for the two component separation with relative volatility of only 1.28.We are going to see, when we run the separation for the binary mixture of ethylbenzene o-xylene, that the reflux ratios of 2.0 and 5.0 will not be enough to get good product ( W-ii ≥ 0.97) and only with the reflux ratio of 10.0 we get the desired molar fraction above 0.97.In order to get more points to analyze we add reflux ratios of 15 and 20 and use the columns with  = 20, 30, 40, and 50.
Using a distillation column with 20 stages with reflux ratio of 2.0 we get good compositions (above 0.97) and we got a harvest time of 0.6 hours and collected 20 kmol of toluene.Extrapolating the composition of toluene at bottom we need 3.3 hours, and using (4) we need 3.0 hours to distillate other 100 kmol that represents all the toluene feed.
In order to simplify the comparison between all runs we are going to assume that (4) predicts well the total time to distillate a product and that the components will remain in equimolar proportion.
All combinations with 20, 30, 40, and 50 stages and reflux ratios of 2, 5, 10, 15, and 20 (A, B, C, D, and E) work fine and the results are shown in Table 5.
Again, for a batch distillation at constant reflux with a fixed number of theoretical stages, when the reflux ratio increases the harvested product increases, but also the harvested and total time increase.The results on harvested product, harvest time, and total time are very similar for all the cases with  = 20, 30, 40, and 50 and then the partial Luyben factor P 3 -CAP for the three component separation is also very close for the same reflux ratio.

Two-Component
Separation.Now only ethylbenzene and o-xylene remain and we are going to assume that we have 200 kmol of an equimolar mixture.
Figure 5 shows the results for the simulation with  = 20 and  = 10.It is observed that molar fraction of ethylbenzene at the distillate remains above 0.97 only 0.5 hours.We only collect 4.55 kmol of ethylbenzene.We need to turn the following distillate to off cut 3 that will contain ethylbenzene and o-xylene.We need to do this until we get o-xylene with molar fraction above 0.97 in the bottom of distillation column.We may extrapolate  W-iii versus time, until it is zero, or we may extrapolate  D-iii versus time until it is zero.The first one will assure that there are no ethylbenzene at the bottom and gives  = 13.69 hours.The other means that there is no ethylbenzene at the bottom and at the plates or packing of the column and will be  = 13.87 hours; we get 73 kmol of o-xylene with  W-iv ≥ 0.97.
Table 6 shows the results for this binary separation of ethylbenzene and xylene including the two harvested products (iii) and (iv), the time for harvest component (iii) at distillate and the time to harvest (iv) at bottom.

Luyben Capacity Factor, (CAP). With the information on
Tables 4-6 the Luyben Capacity factor may be calculated for any combination.By example 20 ABC is the batch distillation column with 20 theoretical stages,  = 2 for the first separation,  = 5 for the second separation, and  = 10 for the third separation.Application of (10) Table 7 provides Luyben capacity factor (CAP) for most of the combination studied and Figure 6 shows best and worst cases.Best Luyben capacity factor uses low reflux ratio and high number of theoretical stages.

Total Annual Cost, (TAC).
The CAPCOST computer program from Turton et al. book [18] was used to estimate the cost of distillation column with sieve trays, condenser, and reboiler.It was assumed that the equipment will last for 10 years and then the contribution for one year was calculated by dividing the cost between the depreciation period.From Table 6.3 the cost for low pressure steam at 5 bars and 160 ∘ C was fixed at 16.22 USD/1000 kg, and the cost for cooling water was 14.8 USD/1000 m 3 .Table 7 shows calculated values for TAC for all the runs and Figure 6 presents best and worst cases.
It is observed that best total annual cost uses low reflux ratio (20AAD, 30AAE) and low number of theoretical stages, although some good combinations use 40 and 50 theoretical stages (40AAE and 50AAC).These alternatives require low energy at the reboiler and offer low TACs.
Luyben capacity factor only takes into account harvest product and processing time.It does not consider cost of the column, reboiler, and condenser, or cost of steam and cooling water.Meanwhile total annual cost it has into account cost of equipment and operational cost of utilities, but ignore harvested product.
This finding of best cases for TAC and Luyben CAP using low reflux ratio and a given number of theoretical stages is analogous to continuous distillation columns that locate the optimum operational reflux ratio at values given by ( 18), Peters and Timmerhaus [19].
As the reflux ratio is increased from the minimum value, the number of plates decreased, the column diameter increases, and the reboiler steam and condenser cooling-water requirements increase.Although the numerical value for reboiler and condenser duty is similar, the annual cost of steam is close to ten times the cost of cooling water.Then, the total annual cost is dominated by the cost of steam.

Annual Profit, (AP).
With harvested data in kmol for each compound, converted to kg and to money with the selling price, then with the total amount of time needed for one run and assuming that in one year the plant will operate 7872 hours (about 90% of the time) the total annual income is calculated.The total annual cost is subtracted to this number and we get the annual profit.Table 7 resumes annual profit for all the runs and Figure 6 shows best and worst cases for the three qualifying parameters.

Analysis
At the beginning of this project only two qualifying parameters: Luyben's capacity factor (CAP) and total annual cost (TAC) were used.The best combinations of theoretical stages and reflux ratio only coincided when using low reflux ratio, probably a reflux slightly above minimum reflux, but it was not clear which number of stages was optimum.Using CAP, batch distillation columns with high number of trays were recommended for eight of the ten best cases, probably because for high number of plates N, the reflux ratio could be small.
Using total annual cost, the best runs used 20 and 30 theoretical stages, although half of the best options used 40 and 50 theoretical stages.
Reading the new Diwekar's book [20] and Diwekar and Madhavan paper [21] we see that optimization is described as a three-step process, and that in the second step a measure for system effectiveness is needed.They propose annual profit given by ( 15), but calculated in a different form.Because we already had the total annual cost, we only needed to calculate the total annual income that was related to Luyben's capacity factor.
We observed that CAP recommends many and TAC recommends few theoretical stages.The combination with highest annual profit AP combined a regular number of theoretical stages (30) and a reflux slightly higher than the minimum, close to 1.2  min .
There were also good combinations with low (20ABE) and high (50ABC) theoretical stages.

Conclusions
The equations required to simulate the short-cut method are few and are easy to program and use in a short computer program that uses relative volatilities as the parameter to predict both the liquid vapor equilibrium and the difficulty of separation.With the short-cut method many cases can be solved by varying the number of theoretical stages and the value of reflux ratio.
To select the best combination, three quality indexes were compared: Luyben's capacity factor, total annual cost, and annual profit.The best combination always required a reflux ratio close to the minimum, which is analogous to the optimum operational reflux equal to 1.1-1.2times the minimum reflux ratio for continuous distillation columns.
The best qualifying parameter was annual profit and, as it used the other two for its calculation, we believe it is representative of the process.
The best combination was case 30ABE.It used a distillation column with 30 plates and a column diameter of 0.91 m, with reflux ratio's of 2.0 (A) for the first separation of four components, 5.0 (B) for the second separation of three components, and 20 (E) for the last separation of two components.The cost of this stainless steel column was estimated at $358,520 USD in 2012.It requires a condenser with an area of 35.85 m 2 at an approximated cost of $119,100 USD.The required stainless steel reboiler should have 119.74m 2 of area and an estimated cost of $810,400 USD.With a depreciation of 0.10 per year the total annual cost is $128,802 USD.
The mass flow rate of cooling water and steam required are 80,087 and 2,443,200 kg/year, respectively, providing an annual operational cost of $2,523,287 USD.Overall, the total annual cost is $2,652,089 USD.
The annual income of the four components was calculated at $17,988,162 USD and subtracting to this value TAC, it should provide an annual profit of $15,336,073 USD before taxes.
From the three quality indexes used only annual profit allowed the selection of the best economical alternative.

Figure 1 :
Figure 1: Batch distillation column with three receivers for products (main cut) and three receivers for mixtures (off cut).

Figure 2 :
Figure 2: Results for four components batch distillation with short-cut method at constant reflux of 2.0; (a) bottom composition, (b) distillate composition.
) and 3(b) show the results for a column with  = 10 and  = 2, 5, 10 while Figures 4(a) and 4(b) show the same results but using  = 40.Figures 3(a) and 4(a) are very similar.They show that for a constant reflux of 2.0 the original bottom is depleted from 400 kmol to 327 and 330 kmol for  = 10 and  = 40, respectively.The difference between the original   = 400 kmol and amount remaining at the bottom appears as distillate harvested = 73 and 70 kmol of a mixture of benzene with a little of toluene for  = 10 and  = 40.

Figure 3 :
Figure 3: Results for four components batch distillation at reflux ratios of 2, 5, and 10 for a column with 10 theoretical stages.

Figure 4 :
Figure 4: Results for four components batch distillation at reflux ratios of 2, 5, and 10 for a column with 40 theoretical stages.

Figure 6 :
Figure 6: Best and worst cases for luyben capacity factor (CAP), total annual cost (TAC), and annual profit (AP).

Table 3 :
Equilibrium constant and relative volatilities for the mixture: benzene, toluene, ethylbenzene, and o-xylene at atmospheric pressure.

Table 7 :
Luyben capacity factor (kmol/h), total annual cost (US dollars), and annual profit (US dollars) for all run tested.
, , , , : Constants for Antoine equation (1) with  in (K) and  in mm Hg   : A r e a o f d o w n c o m e rm 2 : Mole fraction on liquid dimensionless : Mole fraction in vapor dimensionless : Liquid at bottom of column kmol.