Distillation binary calculations

In this chapter consideration is given to the theory of the process, methods of distillation and calculation of the number of stages required for both binary and multicomponent systems, and discussion on design methods is included for plate and packed columns incorporating a variety of column internals. 

For example, a train has only one degree of freedom because only its speed can be varied a boat has two, an airplane three. When looking at industrial processes, the determination of their degrees of freedom becomes more complex and cannot always be determined intuitively. When the process is more complex, as in the case of binary distillation, the calculation of the degrees of freedom also becomes more involved. Figure 2.25 lists 14 variables of this process, but they are not all independent. 

Binary Mixtures.—Since a mixture of constant boiling point behaves like a single substance on distillation, it is possible, if we know the composition of the mixture distilled, to calculate that of the binary mixture. 

This Excel spreadsheet uses the McCabe-Thiele method to calculate the number of theoretical stages needed for binary distillation. Binary distillation is a common unit of operation that separates two liquids (with one being more volatile... 

Zawidski examined 13 binary mixtures using the refractive index of a distillate to calculate the partial pressures and ps, for in his selection both pressures had to be dealt with. His experimental technique is fully described in his paper likewise, I do not need to consider his theoretical arguments it remains only to draw attention to his diagrams. With the exception of the last-named mixture, the mixtures were examined at one temperature only, that given in parentheses ... 

McCabe-Thie/e Example. Assume a binary system E—H that has ideal vapor—Hquid equiHbria and a relative volatiHty of 2.0. The feed is 100 mol of = 0.6 the required distillate is x = 0.95, and the bottoms x = 0.05, with the compositions identified and the lighter component E. The feed is at the boiling point. To calculate the minimum reflux ratio, the minimum number of theoretical stages, the operating reflux ratio, and the number of theoretical stages, assume the operating reflux ratio is 1.5 times the minimum reflux ratio and there is no subcooling of the reflux stream, then ... 

Simple analytical methods are available for determining minimum stages and minimum reflux ratio. Although developed for binary mixtures, they can often be applied to multicomponent mixtures if the two key components are used. These are the components between which the specification separation must be made frequendy the heavy key is the component with a maximum allowable composition in the distillate and the light key is the component with a maximum allowable specification in the bottoms. On this basis, minimum stages may be calculated by means of the Fenske relationship (34) ... 

In order to determine the packed height it is necessary to obtain a value of the overall number of transfer units methods for doing this are available for binary systems in any standard text covering distillation (73) and, in a more complex way, for multicomponent systems (81). However, it is simpler to calculate the number of required theoretical stages and make the conversion ... 

Commercial computer services are available to do rigorous distillation calculations. Perhaps the licensor will provide copies of rigorous computer runs to validate his balances. Alternately, the operating company can make such runs. For highly non-ideal systems, literature data for binary pairs may have to be sought. In some cases, laboratory equilibrium data may have to be obtained in-house or contracted out to one of several organizations or universities that are in this business. 

Relative volatility is the volatility separation factor in a vapor-liquid system, i.e., the volatility of one component divided by the volatility of the other. It is the tendency for one component in a liquid mixture to separate upon distillation from the other. The term is expressed as fhe ratio of vapor pressure of the more volatile to the less volatile in the liquid mixture, and therefore g is always equal to 1.0 or greater, g means the relationship of the more volatile or low boiler to the less volatile or high boiler at a constant specific temperature. The greater the value of a, the easier will be the desired separation. Relative volatility can be calculated between any two components in a mixture, binary or multicomponent. One of the substances is chosen as the reference to which the other component is compared. 

For a binary distillation, a is calculated at top and bottom conditions and a geometric mean used where the differences are relatively small. 

Multicomponent distillation is by far the common requirement for process plants and refineries, rather than the simpler binary systems. There are many computer programs which have been developed to aid in accurately handling the many iterative calculations required when the system involves three to possibly ten individual components. In order to properly solve a multicomponent design, there should be both heat and material balance at every theoretical tray throughout the calculation. 

Although fractional crystallization has always been the most common method for the separation of diastereomers. When it can be used, binary-phase diagrams for the diastereomeric salts have been used to calculate the efficiency of optical resolution. However, its tediousness and the fact that it is limited to solids prompted a search for other methods. Fractional distillation has given only limited separation, but gas chromatography and preparative liquid chromatography have proved more useful and, in many cases, have supplanted fraetional crystallization, especially where the quantities to be resolved are small. 

Bubble Point Calculation Binary Batch Distillation Column... 

Figure 1.8. Information flow, binary flash distillation calculation (a) Information recycle (b) Information flow...

Integral condensation in which the liquid remains in equilibrium with the uncondensed vapour. The condensation curve can be determined using procedures similar to those for multicomponent flash distillation given in Chapter 11. This will be a relatively simple calculation for a binary mixture, but complex and tedious for mixtures of more than two components. 

A distillation calculation is to be performed on a multicomponent mixture. The vapor-liquid equilibrium for this mixture is likely to exhibit significant departures from ideality, but a complete set of binary interaction parameters is not available. What factors would you consider in assessing whether the missing interaction parameters are likely to have an important effect on the calculations ... 

It is desired to separate a binary mixture by simple distillation. If the feed mixture has a composition of 0.5 mole fraction, calculate the fraction which it is necessary to vaporise in order to obtain ... 

Continuous binary distillation is illustrated by the simulation example CON-STILL. Here the dynamic simulation example is seen as a valuable adjunct to steady state design calculations, since with MADONNA the most important column design parameters (total column plate number, feed plate location and reflux ratio) come under the direct control of the simulator as facilitated by the use of sliders. Provided that sufficient simulation time is allowed for the column conditions to reach steady state, the resultant steady state profiles of composition versus plate number are easily obtained. In this way, the effects of changes in reflux ratio or choice of the optimum plate location on the resultant steady state profiles become almost immediately apparent. 

Process Modes for a Bioreactor 538 Binary Batch Distillation Column 490 Bubble Point Calculation for a Batch Distillation Column 504... 

Key concepts of the calculation of distillation are well illustrated by analysis of the distillation of binary mixtures. Moreover, many real systems are essentially binary or can be treated as binaries made up of two pseudo components, for which it is possible to calculate upper and lower limits to the equipment size for a desired separation. 

The calculational base consists of equilibrium relations and material and energy balances. Equilibrium data for many binary systems are available as tabulations of x vs. y at constant temperature or pressure or in graphical form as on Figure 13.4. Often they can be extended to other pressures or temperatures or expressed in mathematical form as explained in Section 13.1. Sources of equilibrium data are listed in the references. Graphical calculation of distillation problems often is the most convenient... 

The concepts NTU and HTU are defined only for binary distillations and the transfer of a single substance in absorption or stripping. Since most processes of industrial interest involve multicomponents, the HETS of packed towers is the more useful concept, and may be evaluated readily from test data and tray calculations. 

Measurements of binary vapor-liquid equilibria can be expressed in terms of activity coefficients, and then correlated by the Wilson or other suitable equation. Data on all possible pairs of components can be combined to represent the vapor-liquid behavior of the complete mixture. For exploratory purposes, several rapid experimental techniques are applicable. For example, differential ebulliometry can obtain data for several systems in one laboratory day, from which infinite dilution activity coefficients can be calculated and then used to evaluate the parameters of correlating equations. Chromatography also is a well-developed rapid technique for vapor-liquid equilibrium measurement of extractive distillation systems. The low-boiling solvent is deposited on an inert carrier to serve as the adsorbent. The mathematics is known from which the relative volatility of a pair of substances can be calculated from the effluent trace of the elutriated stream. Some of the literature of these two techniques is cited by Walas (1985, pp. 216-217). 

Calculation Methods. An often satisfactory approximation is to take the mixture in the presence of the solvent to be a pseudobinary of the keys on a solvent-free basis, and to employ the McCabe-Thiele or other binary distillation method to find tray and reflux demands. Since the relative volatility varies with concentration of the solvent, different equilibrium curves are used for above and below the feed based on average loads in those zones. Figure 13.25 is of such a construction. 

For the synthesis of heterogeneous batch distillation the liquid-liquid envelope at the decanter temperature is considered in addition to the residue curve map. Therefore, the binary interaction parameters used in predicting liquid-liquid equilibrium are estimated from binary heterogeneous azeotrope or liquid-liquid equilibrium data [8,10], Table 3 shows the calculated purity of original components in each phase split at 25 °C for all heterogeneous azeotropes reported in Table 1. The thermodynamic models and binary coefficients used in the calculation of the liquid-liquid-vapour equilibrium, liquid-liquid equilibrium at 25 °C and the separatrices are reported in Table 2. 

For a binary mixture and given that the relative volatility (a) between the components is constant, instantaneous distillate composition y can be calculated using ... 

The one level optimal control formulation proposed by Mujtaba (1989) is found to be much faster than the classical two-level formulation to obtain optimal recycle policies in binary batch distillation. In addition, the one level formulation is also much more robust. The reason for the robustness is that for every function evaluation of the outer loop problem, the two-level method requires to reinitialise the reflux ratio profile for each new value of (Rl, xRI). This was done automatically in Mujtaba (1989) using the reflux ratio profile calculated at the previous function evaluation in the outer loop so that the inner loop problems (specially problem P2) could be solved in a small number of iterations. However, experience has shown that even after this re-initialisation of the reflux profile sometimes no solutions (even sub-optimal) were obtained. This is due to failure to converge within a maximum limit of function evaluations for the inner loop problems. On the other hand the one level formulation does not require such re-initialisation. The reflux profile was set only at the beginning and a solution was always found within the prescribed number of function evaluations. 

Figure 11.5 illustrates the operating sequence for a binary mixture. Successive passes (multi-pass) are used sequentially using the same column to separate only component. The purity of the distillate product remains the same for each pass but the distillate rate varies (this in turn varies the reflux ratio). This strategy is similar to a time sequenced reflux ratio operation for individual cuts in CBD operation. The total recovery of a component (say A) is calculated from the accumulated amount of distillate from all the passes. 

Up to now, the proposed model has been validated using the total reflux distillation data in the column equipped with the Montz-Pak A3-500 and Montz-Pak B1-250 structured packings. A very good agreement between the calculated and experimental data for binary and ternary mixtures is found. 

Last, the point was made early on that the separation could be calculated on the basis of a light key, but all of the examples were based on a binary mixture. Simple distillation of a multicomponent mixture is analyzed on the basis of the light key as follows,... 

---