Distillation
Distillation
where pv and pV1 are the vapor pressures in Pascal at absolute temperature T and T1 in K. λ is
the molar latent heat of vaporization which is independent of temperature.
Antoine Equation:
‣ Typical representative values of the constants A, B and C are given in the following Table 2.1 (Ghosal et al., 1993).
‣ Table 2.1: Typical representative values of the constants A, B and C
Phase Diagram
‣ For binary mixture phase diagram only two-component mixture, (e.g. A (more volatile) and B
(less volatile)) are considered.
‣ There are two types of phase diagram: constant pressure and constant temperature.
Constant Pressure Phase Diagram
‣ The Figure 2.1 shows a constant pressure phase diagram for an ideal solution (one that obeys
Raoult's Law).
‣ At constant pressure, depending on relative concentrations of each component in the liquid,
many boiling point temperatures are possible for mixture of liquids (solutions) as shown in
phase diagram (Figure 2.1).
‣ For mixture, the temperature is called bubble point temperature when the liquid starts to boil
and dew point when the vapor starts to condense.
‣ Boiling of a liquid mixture takes place over a range of boiling points. Likewise, condensation
of a vapor mixture takes place over a range of condensation points.
‣ The upper curve in the boiling point diagram is called the dew-point curve (DPC) while the
‣ At each temperature, the vapor and the liquid are in equilibrium.
‣ The constant pressure phase diagram is more commonly used in the analysis of vapor-liquid
equilibrium.
‣ where, yA = mole fraction of component ‘A’ in the vapor, xA = mole fraction of component ‘A’
in the liquid.
‣ In general, relative volatility of a mixture changes with the mixture composition. For binary
mixture, xB = 1-xA.
‣ where PA is the partial pressure of component A in the vapor, PB is the partial
pressure of component B in the vapor and P is the total pressure of the system.
‣ Thus if the relative volatility between two components is equal to one, separation
is not possible by distillation.
‣ The larger the value of α, above 1.0, the greater the degree of separability, i.e. the
easier the separation.
‣ Some examples of optimal relative volatility that are used for distillation process design are
given in Table 2.2.
Table 2.2: Some optimal relative volatility that are used for distillation process
Vapor-Liquid Equilibrium (VLE)
‣ It is useful for graphical design in determining the
number of theoretical stages required for a distillation
column.
‣ A typical equilibrium curve for a binary mixture on x-y
plot is shown in Figure 2.3.
‣ It contains less information than the phase diagram (i.e.
temperature is not included), but it is most commonly
used.
‣ The VLE plot expresses the bubble-point and the dew-
point of a binary mixture at constant pressure.
‣ The curved line in Figure 2.3 is called the equilibrium
Figure 2.3: Equilibrium line and it describes the compositions of the liquid and
diagram for a benzene-toluene vapor in equilibrium at some fixed pressure.
mixture at 1 atmosphere
Vapor-Liquid Equilibrium (VLE)
‣ This particular VLE plot (Figure 2.3) shows a binary
ideal mixture that has a uniform vapor-liquid
equilibrium that is relatively easy to separate.
‣ On the other hand, the VLE plots shown in Figure 2.4
represented for non-ideal systems.
‣ These non-ideal VLE systems will present more
difficult separation.
‣ The most challenging VLE curves are generated by
azeotropic systems. An azeotrope is a liquid mixture
which when vaporized, produces the same
composition as the liquid.
‣ Two types of azeotropes are known: minimum-boiling
Figure2.4:Vapor-liquid and maximum-boiling (less common).
equilibrium curve for non-ideal
systems
Vapor-Liquid Equilibrium (VLE)
‣ The most challenging VLE curves are generated by azeotropic systems.
‣ An azeotrope is a liquid mixture which when vaporized, produces the same composition as the
liquid.
‣ Two types of azeotropes are known: minimum-boiling and maximum-boiling (less common).
Ethanol-water system (at 1 atm, 89.4 mole %, 78.2 oC; Carbon-disulfide - acetone (61.0 mole
% CS2, 39.25 oC, 1 atm) and Benzene - water (29.6 mole% H2O, 69.25 oC, 1 atm) are
minimum-boiling azeotropes. Hydrochloric acid - water (11.1 mole% HCl, 110 oC, 1 atm);
‣ Acetone - chloroform (65.5 mole% chloroform, 64.5 oC, 1 atm) are the examples of
Maximum-boiling azeotropes.
‣ The Figure 2.5 shows two different azeotropic systems, one with a minimum boiling point
(Figure 2.5a) and one with a maximum boiling point (Figure 2.5b).
‣ The points of intersections of the equilibrium curves with the diagonal lines are called
azeotropic points.
‣ An azeotrope cannot be separated by conventional distillation. However, vacuum distillation
may be used as the lower pressures which can shift the azeotropic point.
Figure 2.5: VLE curves for azeotropic systems: (a) for maximum boiling point, (b) for minimum boiling point
‣ Although most distillations are carried out at atmospheric or near atmospheric pressure, it is
not uncommon to distill at other pressures.
‣ High pressure distillation (typically 3 - 20 atm) usually occurs in thermally integrated
processes.
‣ In those cases the equilibrium curve becomes narrower at higher pressures as shown in Figure
2.6. Separability becomes less at higher pressures.
Distillation columns and their process calculations
‣ There are many types of distillation columns each of which is designed to perform specific
types of separations.
‣ One way of classifying distillation column type is to look at how they are operated. Based on
operation, they are of two types: batch or differential and continuous columns.
Batch or differential distillation columns and their process calculation
‣ In batch operation, the feed is introduced batch-wise to the column. That is, the column is
charged with a 'batch' and then the distillation process is carried out.
‣ When the desired task is achieved, a next batch of feed is introduced. Consider a binary
mixture of components A (more volatile) and B (less volatile).
‣ The system consists of a batch of liquid (fixed quantity) inside a kettle (or still) fitted with
heating element and a condenser to condense the vapor produced as shown in Figure 2.7.
‣ The condensed vapor is known as the distillate.
‣ The distillate is collected in a condensate receiver. The liquid remaining in the still is known as
the residual.
‣ The process is unsteady state.
‣ The concentration changes can be analyzed using the phase diagram, and detailed
mathematical calculations carried out using the Rayleigh Equation.
‣ As the process is unsteady state, the derivation is based on a differential approach to changes in
concentration with time.
‣ Let L1 = initial moles of liquid originally in still, L2 = final moles of liquid remained in still,
x1 = initial liquid composition in still (mole fraction of A), x2 = final liquid composition in
still (mole fraction A).
‣ At any time t, the amount of liquid in the still is L, with mole fraction of A in the liquid being
x. After a small differential time (t + dt), a small amount of vapor dL is produced, and the
composition of A in the vapor is y (mole fraction).
‣ The vapor is assumed to be in equilibrium with the residue liquid.
‣ The amount of liquid in the still is thus reduced from L to (L - dL), while the liquid
composition changed from x to (x - dx).
‣ Then the material balance on A can be written as: Initial amount in still = Amount left in still +
Amount vaporized
‣ where HF, HV and HL are enthalpies of feed, vapor and liquid respectively which can be
obtained from enthalpy-concentration diagram for the mixture.
2. When feed is cold liquid or superheated vapor:
‣ q can be alternatively defined as the heat required to convert 1 mole of feed from its entering
condition to a saturated vapor; divided by the molar latent heat of vaporization.
‣ Based on this definition, one can calculate the q-value from the following Equations for the
case whereby q > 1 (cold liquid feed) and q < 0 (superheated vapor feed) as
2.1.For cold liquid feed:
‣ Component balances for the more volatile component in the rectifying and stripping sections
are:
‣ For rectifying section:
‣ For stripping section:
‣ At the feed point where the two operating lines intersect can be written as:
Figure 2.14: Feed tray with fraction of liquid and vapor of feed
‣ From component balance around the entire column, it can be written as
‣ For a given feed condition, xF and q are fixed, therefore the q-line is a straight line with slope -
q / (1-q) and intercept xF/(1-q). If x = xF , then from Equation above y = xF.
‣ At this condition the q-line passes through the point (xF, xF) on the 45 diagonal.
‣ Different values of q will result in different slope of the q-line.
‣ Different q-lines for different feed conditions are shown in Figure 2.15.
Figure 2.15: Different q-lines for different feed conditions
Determination of the stripping section operating line (SOL)
‣ The stripping section operating line (SOL) can be obtained from the ROL and qline without
doing any material balance.
‣ The SOL can be drawn by connecting point xB on the diagonal to the point of intersection
between the ROL and q-line.
‣ The SOL will change if q-line is changed at fixed ROL.
‣ The change of SOL with different q-lines for a given ROL at constant R and xD is shown in
Figure 2.16.
‣ It is applied for the whole column. Every tray is assumed to have the same efficiency.
‣ The overall efficiency depends on the (i) geometry and design of the contacting trays, (ii) flow
rates and flow paths of vapor and liquid streams, (iii) Compositions and properties of vapor
and liquid streams (Treybal, 1981; Seader and Henley, 1998).
‣ The overall efficiency can be calculated from the following correlations:
‣ The Drickamer-Bradford empirical correlation:
Murphree Efficiency
‣ The efficiency of the tray can also be calculated based on semi-theoretical models which can
be interpreted by the Murphree Tray Efficiency EM.
‣ In this case it is assumed that the vapor and liquid between trays are well-mixed and have
uniform composition.
‣ It is defined for each tray according to the separation achieved on each tray based on either the
liquid phase or the vapor phase. For a given component, it can be expressed as:
‣ Based on vapor phase:
‣ The points represented by the feed, distillate and bottom streams can be plotted on the
enthalpy-concentration diagram as shown in Figure 2.21.
‣
‣
Figure 2.21: Representation of feed, distillate and bottom streams on the enthalpy-concentration
Analysis of binary distillation by Ponchon-Savarit Method
‣ Comparing the last Equation with the points plotted on Figure (2.21), it is found that the left hand side of the last
Equation represents the slope of the straight line between the points (xB, HB-QR/B) and (xF, HF).
‣ The right hand side of the Equation represents the slope of the straight line passing through the points (xF, HF)
and (xD, HD-QC/D).
‣ From this it can be said that all three points are on a same straight line.
‣ The amount of distillate as per Equation is proportional to the horizontal distance and the amount of bottoms is
proportional to the horizontal distance FD xx , then from the overall material balance it can be interpreted that the
amount of feed is proportional to the horizontal distance BD xx .
‣ This leads to the inverse lever rule. Based on this principle, the analysis of distillation column is called Enthalpy-
composition analysis or Ponchon-Savarit analysis of distillation column.
‣ The sum of the liquid and vapor leaving the plate must equal the total flow to the plate. So the
tie line must pass through the point P which represents the addition point of the vapor Vn and
liquid Ln leaving the tray.
‣ The intersection of the tie line with the enthalpy-composition curves will represent the
compositions of these streams.
‣ The distance is proportional to the quantity Ln-1 and the distance is proportional to
the quantity Vn+1 as per lever rule.
‣ So
‣ The sum of the liquid and vapor leaving the plate must equal the total flow to the plate. So the
tie line must pass through the point P which represents the addition point of the vapor Vn and
liquid Ln leaving the tray.
‣ The intersection of the tie line with the enthalpy-composition curves will represent the
compositions of these streams.
Figure 2.22: Operation principle of stage-wise binary distillation on enthalpy-concentration
diagram
Introduction to Multicomponent Distillation
‣ In industry, most of the distillation processes involve with more than two components. The
multicomponent separations are carried out by using the same type of distillation columns,
reboilers, condensers, heat exchangers and so on.
‣ However some fundamental differences are there which is to be thoroughly understood by the
designer. These differences are from phase rule to specify the thermodynamic conditions of a
stream at equilibrium.
‣ In multicomponent systems, the same degree of freedom is not achieved because of the
presence of other components.
‣ Neither the distillate nor the bottoms composition is completely specified. The components
that have their distillate and bottoms fractional recoveries specified are called key components.
‣ The most volatile of the keys is called the light key (LK) and the least volatile is called the
heavy key (HK).
‣ The other components are called non-keys (NK). Light no-key (LNK) is referred when non-
key is more volatile than the light key whereas heavy non-key (HNK) is less volatile than the
heavy key. Proper selection of key components is important if a multicomponent separation is
‣ Several short-cut methods are used for carrying out calculations in multicomponent systems.
‣ These involve generally an estimation of the minimum number of trays, the estimation of
minimum reflux rate and number of stages at finite reflux for simple fractionators.
‣ Although rigorous computer methods are available to solve multicomponent separation
problems, approximate methods are used in practice.
‣ A widely used approximate method is commonly referred to as the Fenske-Underwood-
Gilliland method.
Estimation of Minimum number of trays: Fenske Equation
‣ Fenske (1932) was the first to derive an Equation to calculate minimum number of trays for multicomponent distillation at
total reflux.
‣ The derivation was based on the assumptions that the stages are equilibrium stages.
‣ Consider a multicomponent distillation column operating at total reflux as shown in Figure 5.24.
‣ Equilibrium relation for the light key component on the top tray is
Figure 2.24: Multicomponent column at minimum trays
Minimum Reflux: Underwood Equations
‣ For multi-component systems, if one or more of the components appear in only one of the
products, there occur separate pinch points in both the stripping and rectifying sections.
‣ In this case, Underwood developed an alternative analysis to find the minimum reflux ratio
(Wankat, 1988).
‣ The presence of non-distributing heavy non-keys results a pinch point of constant composition
at minimum reflux in the rectifying section whereas the presence of non-distributing light non-
keys, a pinch point will occur in the stripping section.
‣ Let us consider the pinch point is in the rectifying section.
‣ The mass balance for component i around the top portion of the rectifying section as illustrated
in Figure 5.24 is
Analysis of binary distillation by Ponchon-Savarit Method
Analysis of binary distillation by Ponchon-Savarit Method
Analysis of binary distillation by Ponchon-Savarit Method