Abstract

A short-cut method for batch distillation columns working at constant reflux was applied to solve a problem of 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 the minimum. 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).

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

2. 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 o-xylene (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.

2.2. 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 short-cut 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 (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 () 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 on-specification 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.

2.3. 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 to obtain and ()/() to obtain for each point. From these values we calculated the relative volatility for each reported point. Higher values of ,  , 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 (), while the simpler separation was benzene/toluene ().

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] 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 increased from 6.33 to 7.16.

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.

2.4. Quality Indexes for Batch Distillation Column Operation
2.4.1. 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 (, , , and ) that are collected on the main cut receivers, the time required to eliminate a volatile component (, , and ) 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 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.

2.4.2. Total Annual Cost for Comparison

Total annual cost given by 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  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: = Number of theoretical stages, = Reflux ratio used in the first separation: distillation of benzene (i) from (ii), (iii), and (iv) (4 components). For this separation we used , , , , , = Reflux ratio used in the second separation: distillation of (ii) from (iii) and (iv) (3 components), = Reflux ratio used in the third and last separation of (iii) from (iv) (2 components). For this separation the reflux ratios 2.0 and 5.0 are smaller than the minimum reflux and the distillation column cannot operate, then we used , , .

As an example, 20 ABC means a distillation column with 20 theoretical stages, (A) for the first separation (four components), (B) for the second separation (three components), and (C) for the third separation (two components).

2.4.3. 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.

3. Results and Discussion

3.1. 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 and compare it with the result of (h) in Table 2 for minimum reflux ratio with the Underwood correlation . Once (g) and (h) in Table 2 were solved, the results were used to find the number of minimum stages () that should be applied in other steps of the method.

In the first iterations we began with , minimum reflux from Gilliland-Eduljee provided a value of and . Using successive substitution we calculated the new values and for this value we obtained and and a new value for . The last values for the first iteration were ,  , and .

In the process, time was increasing and decreased from 6.77 to 2.25 at the last time reported (2.2 hours), but minimum reflux () 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 () 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) and 3(b) show the results for a column with and while Figures 4(a) and 4(b) show the same results but using .

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 and , respectively. The difference between the original 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 and .

For reflux ratios of 5 and 10 and the same number of theoretical stages, more distillate product with 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 , we extrapolated the data of versus time in Figure 2(a) obtaining first equation: And then solving for , 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 , we get 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: .

Applying the same procedure for and provided and 11 hours, respectively to eliminate the benzene from the bottom of the column (). 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 and low number of theoretical stages distillate composition drops below 0.90, but for the combination of and the punctual composition for benzene in distillate is always close to unity for the three values of reflux ratio. Similar results were performed for 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 ) 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 (P4-CAP) for the four components separation by dividing the collected benzene product by the total time for the run.

3.2. Three Components Separation

Beginning with  kmol we run the program first for and 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 () 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 , 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 , 30, 40, and 50 and then the partial Luyben factor P3-CAP for the three component separation is also very close for the same reflux ratio.

3.3. 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 and . 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 versus time, until it is zero, or we may extrapolate versus time until it is zero. The first one will assure that there are no ethylbenzene at the bottom and gives hours. The other means that there is no ethylbenzene at the bottom and at the plates or packing of the column and will be hours; we get 73 kmol of o-xylene with .

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.

3.4. Luyben Capacity Factor, (CAP)

With the information on Tables 46 the Luyben Capacity factor may be calculated for any combination. By example 20 ABC is the batch distillation column with 20 theoretical stages, for the first separation, for the second separation, and for the third separation. Application of (10) gives the capacity factor en kmol/h as follows: 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.

3.5. 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 m3.

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.

3.6. 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.

4. 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 .

There were also good combinations with low (20ABE) and high (50ABC) theoretical stages.

5. 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.2 times 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 m2 at an approximated cost of $119,100 USD. The required stainless steel reboiler should have 119.74 m2 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.

Nomenclature

,  , , , :Constants for Antoine equation (1) with   in (K) and   in mm Hg
: Area of downcomer  m2
: Area of condenser  m2
AP: Annual profit  USD
: Area of reboiler  m2
CAP: Luyben capacity factor  kmol/h, or kg/h
: Heat capacity of cooling water  K/kg-K
D: Distillate  kmol
: Diameter of column  m
:Fractionoperational vapor velocity  dimensionless
(i),(ii), (iii), (iv): Benzene,  toluene,  ethylbenzene,  xylene
: Equilibriumconstant   dimensionless
: Liquid flow  kmol/s
: Flow rate of cooling water  kg/s
: Molecular mass  kg/kmol
: Flow rate of steam  kg/s
: Number of theoretical  stages  dimensionless
: Number of components  dimensionless
: Minimum number of theoretical stages  dimensionless
: Pressure of saturation for component   bar
: Total pressure  bar
: Heat at condenser  J/h
: Heat at reboiler  J/h
Reflux ratio dimensionless
: Minimum  reflux ratio  dimensionless
Time  h
: Time of equilibration  h
: Time of finalization  h
: Time of process  h
TAC: Total annual costUSDollar
: Temperature of  cooling water  K
: Temperature of heating steam  K
: Superficial velocity  m/s
: Global heat coefficient  J/s-m2-K
: Flooding velocity of vapor  m/s
: Boilup vapor  kmol/h
: Mole fraction  on liquid  dimensionless
: Mole fraction in vapor  dimensionless
: Liquid at bottom of column  kmol.
Greek Letters
: Relative volatility ()  dimensionless
: Density  kg/m3
: Viscosity  kg/m-s
Heat of condensation or vaporization  J/kg-K.
Subscripts
0: Starting time
: Component
: Feed
: Liquid phase
: Distillate
: Bottom.

Acknowledgments

The authors acknowledge David Caliz Rodriguez (Student of Instituto Tecnologico de VillaHermosa, Mexico) for the implementation of the first Excel computer program in 2011 Scientific Summer Mexican Program and Jessica Ortiz Palafox (Student of FIQ-UADY) for the implementation of the computer program in Mathcad to estimate distillation column parameters and cost functions.