ENTROPY AND EQUILIBRIUM 189

of state for both phases (for example, the SRK equation), and the VLE-condition
µg,i = µl,i gives
φVi yi = φL
i xi (7.44)
where the fugacity coefficients φVi and φLi are determined from the equation of state.
The K value is then Ki = φLi /φV
i . Note that (7.44) can also be used for supercritical
components.

7.5 Flash calculations

V
yi

p, T
F
zi

L
xi

Figure 7.4: Flash tank

Flash calculations are used for processes with vapor/liquid-equilibrium (VLE). A

typical process that requires flash calculations, is when a feed stream (F ) is separated
into a vapor (V ) and liquid (L) product; see Figure 7.4.
In principle, flash calculations are straightforward and involve combining the VLE-
equations with the component mass balances, and in some cases the energy balance.
Some flash calculations are (with a comment on their typical numerical solution or
usage):
1. Bubble point at given T (easy)
2. Bubble point at given p (need to iterate on T )
3. Dew point at given T (easy)
4. Dew point at given p (need to iterate on T )
5. Flash at given p and T (relatively easy)
6. Flash at given p and H (“standard” flash, e.g., for a flash tank after a valve)
7. Flash at given p and S (e.g., for condensing turbine)
8. Flash at given U and V (e.g., for dynamic simulation of an adiabatic flash drum)
The last three flashes are a bit more complicated as they require the use of the energy
balance and relationships for computing H, S, etc. The use of flash calculations is best
illustrated by some examples. Here, we assume that the VLE is given on K-value form,
that is,
yi = Ki xi
190 CHEMICAL AND ENERGY PROCESS ENGINEERING

Table 7.2: Data for flash examples and exercises: Antoine parameters for psat (T ), normal
boiling temperature (Tb ) and heat of vaporization ∆vap H(Tb ) for selected components. Data:
Poling, Prausnitz and O’Connell, The properties of gases and liquids, 5th Ed., McGraw-Hill (2001).

% log10(psat[bar])=A-B/(T[K]+C) Tb[K] dvapHb [J/mol]

A1=3.97786; B1=1064.840; C1=-41.136; Tb1=309.22; dvapHb1=25790; % pentane C5H12
A2=4.00139; B2=1170.875; C2=-48.833; Tb2=341.88; dvapHb2=28850; % hexane C6H14
A3=3.93002; B3=1182.774; C3=-52.532; Tb3=353.93; dvapHb3=29970; % cyclohex C6H12
A4=5.20277; B4=1580.080; C4=-33.650; Tb4=337.69; dvapHb4=35210; % methanol CH3OH
A5=5.11564; B5=1687.537; C5=-42.98; Tb5=373.15; dvapHb5=40660; % water H2O
A6=4.48540; B6= 926.132; C6=-32.98; Tb6=239.82; dvapHb6=23350; % ammonia NH3
A7=3.92828; B7= 803.997; C7=-26.11; Tb7=231.02; dvapHb7=19040; % propane C3H8
A8=4.05075; B8=1356.360; C8=-63.515; Tb8=398.82; dvapHb8=34410; % octane C8H18
A9=4.12285; B9=1639.270; C9=-91.310; Tb9=489.48; dvapHb9=43400; % dodecane C12H26
A10=3.98523; B10=1184.24; C10=-55.578; Tb10=353.24; dvapHb11=30720; % benzene C6H6
A11=4.05043; B11=1327.62; C11=-55.525; Tb11=383.79; dvapHb11=33180; % toluene C7H8

where yi is the vapor phase mole fraction and xi the liquid phase mole fraction for
component i. In general, the “K-value” Ki depends on temperature T , pressure p and
composition (both xi and yi ). We mostly assume ideal mixtures, and use Raoult’s law.
In this case Ki depends on T and p only:

Raoult′ s law : Ki = psat

i (T )/p

In the examples, we compute the vapor pressure psat (T ) using the Antoine parameters
given in Table 7.2.

7.5.1 Bubble point calculations

Let us first consider bubble point calculations, In this case the liquid-phase
composition xi is given (it corresponds to the case where V is very small (V ! 0)
and xi = zi in Figure 7.4). The bubble point of a liquid is the point where the liquid
just starts to evaporate (boil), that is, when the first vapor bubble is formed. If the
temperature is given, then we must lower the pressure until the first bubble is formed.
If the pressure is given, then we must increase the temperature until the first bubble
is formed. In both cases, this corresponds to adjusting T or p until the computed sum
of vapor fractions is just 1, Σyi = 1 or

Σi Ki xi = 1 (7.45)

where xi is given. For the ideal case where Raoult’s law holds this gives

Σi xi psat (T ) = p (7.46)
! i"# $
pi

Example 7.17 Bubble point at given temperature T . A liquid mixture contains 50%
pentane (1), 30% hexane (2) and 20% cyclohexane (3) (all in mol-%), i.e.,

x1 = 0.5; x2 = 0.3; x3 = 0.2

ENTROPY AND EQUILIBRIUM 191

At T = 400 K, the pressure is gradually decreased. What is the bubble pressure and
composition of the first vapor that is formed? Assume ideal liquid mixture and ideal gas
(Raoult’s law).
Solution. The task is to find a p that satisfies (7.46). Since T is given, this is trivial; we
can simply calculate p from (7.46). We start by computing the vapor pressures for the three
components at T = 400K. Using the Antoine data in Table 7.2, we get:

psat
1 (400K) = 10.248 bar

psat
2 (400K) = 4.647 bar

psat
3 (400K) = 3.358 bar

At the bubble point, the liquid phase composition is given, so the partial pressure of each
component is
p1 = x1 psat
1 = 5.124 bar

p2 = x2 psat
2 = 1.394 bar

p3 = x3 psat
3 = 0.672 bar

Thus, from (7.46) the bubble pressure is

p = p1 + p2 + p3 = 7.189 bar

Finally, the vapor composition (composition of the first vapor bubble) is

p1 p2 p3
y1 = = 0.713; y2 = = 0.194; y3 = = 0.093
p p p
For calculation details see the MATLAB code:
T=400; x1=0.5; x2=0.3; x3=0.2
psat1=10^(A1-B1/(T+C1)), psat2=10^(A2-B2/(T+C2)), psat3=10^(A3-B3/(T+C3))
p1=x1*psat1, p2=x2*psat2, p3=x3*psat3, p=p1+p2+p3
y1=p1/p, y2=p2/p, y3=p3/p

Example 7.18 Bubble point at given pressure p. Consider the same liquid mixture
with 50% pentane (1), 30% hexane (2) and 20% cyclohexane (3) (all in mol-%). A p = 5
bar, the temperature is gradually increased. What is the bubble temperature and composition
of the first vapor that is formed?
Solution. In this case, p and xi are given, and (7.46) provides an implicit equation for T
which needs to be solved numerically, for example, by iteration. A straightforward approach
is to use the method from the previous example, and iterate on T until the bubble pressure is
5 bar (for example, using the MATLAB code below). We find T = 382.64 K, and
p1 p2 p3
y1 = = 0.724; y2 = = 0.187; y3 = = 0.089
p p p
% MATLAB:
x1=0.5; x2=0.3; x3=0.2; p=5;
T=fzero(@(T) p-x1*10^(A1-B1/(T+C1))-x2*10^(A2-B2/(T+C2))-x3*10^(A3-B3/(T+C3)) , 400)

7.5.2 Dew point calculations

Let us next consider dew point calculations. In this case the vapor-phase composition
yi is given (it corresponds to the case where L is very small (L ! 0) and yi = zi in
Figure 7.4). The dew point of a vapor (gas) is the point where the vapor just begins
192 CHEMICAL AND ENERGY PROCESS ENGINEERING

to condense, that is, when the first liquid drop is formed. If the temperature is given,
then we must increase the pressure until the first liquid is formed. If the pressure is
given, then we must decrease the temperature until the first liquid is formed. In both
cases, this corresponds to adjusting T or p until Σxi = 1 or

Σi yi /Ki = 1 (7.47)

where yi is given. For an ideal mixture where Raoult’s law holds this gives
yi 1
Σi = (7.48)
psat
i (T ) p

Example 7.19 Dew point at given temperature T . A vapor mixture contains 50%
pentane (1), 30% hexane (2) and 20% cyclohexane (3) (all in mol-%), i.e.,

y1 = 0.5; y2 = 0.3; y3 = 0.2

At T = 400 K, the pressure is gradually increased. What is the dew point pressure and
the composition of the first liquid that is formed? Assume ideal liquid mixture and ideal gas
(Raoult’s law).
Solution. The task is to find the value of p that satisfies (7.48). Since T is given, this is
trivial; we can simply calculate 1/p from (7.48). With the data from Example 7.17 we get:
1 0.5 0.3 0.2
= + = = 0.1729bar−1
p 10.248 4.647 3.358

and we find p = 5.78 bar. The liquid phase composition is xi = yi p/psat

i (T ) and we find

0.5 · 5.78 0.3 · 5.78 0.2 · 5.78

x1 = = 0.282, x2 = = 0.373, x3 = = 0.345
10.248 4.647 3.749
% MATLAB:
T=400; y1=0.5; y2=0.3; y3=0.2
psat1=10^(A1-B1/(T+C1)), psat2=10^(A2-B2/(T+C2)), psat3=10^(A3-B3/(T+C3))
p=1/(y1/psat1 + y2/psat2 + y3/psat3)
x1=y1*p/psat1, x2=y2*p/psat2, x3=y3*p/psat3

Example 7.20 Dew point at given pressure p. Consider the same vapor mixture with
50% pentane (1), 30% hexane (2) and 20% cyclohexane (3). At p = 5 bar, the temperature is
gradually decreased. What is the dew point temperature and the composition of the first liquid
that is formed?
Solution. In this case, p and yi are given, and (7.48) provides an implicit equation for
T which needs to be solved numerically (e.g., using the MATLAB code below). We find
T = 393.30 K, and from xi = yi p/psat
i (T ) we find

x1 = 0.278; x2 = 0.375; x3 = 0.347

% MATLAB:
y1=0.5; y2=0.3; y3=0.2; p=5;
T=fzero(@(T) 1/p-y1/10^(A1-B1/(T+C1))-y2/10^(A2-B2/(T+C2))-y3/10^(A3-B3/(T+C3)) , 400)

Example 7.21 Dew point with non-condensable components. Calculate the

temperature and composition of a liquid in equilibrium with a gas mixture containing 10%
pentane (1), 10% hexane and 80% nitrogen (3) at 3 bar. Nitrogen is far above its critical
point and may be considered non-condensable.
ENTROPY AND EQUILIBRIUM 193

Solution. To find the dew-point we use Σi xi = 1. However, nitrogen is assumed non-

condensable so x3 = 0. Thus, this component should not be included in (7.48), which becomes
y1 y2 1
+ sat =
psat
1 (T ) p2 (T ) p
Solving this implicit equation in T numerically (e.g., using the MATLAB code below) gives
T = 314.82K and from xi = yi p/psati (T ) the liquid composition is

x1 = 0.245; x2 = 0.755; x3 = 0

7.5.3 Flash with liquid and vapor products

Next, consider a flash where a feed F (with composition zi ) is split into a vapor
product V (with composition yi ) and a liquid product (with composition xi ); see
Figure 7.4 on page 189. For each of the Nc components, we can write a material
balance
F zi = Lxi + V yi (7.49)
In addition, the vapor and liquid is assumed to be in equilibrium,
yi = Ki xi
The K-values Ki = Ki (T, P, xi , yi ) are computed from the VLE model. In addition,
we have the two relationships Σi xi = 1 and Σi yi = 1. With a given feed (F, zi ), we
then have 3Nc + 2 equations in 3Nc + 4 unknowns (xi , yi , Ki , L, V, T, p). Thus, we need
two additional specifications, and with these the equation set should be solvable.

pT -flash
The simplest flash is usually to specify p and T (pT -flash), because Ki depends
mainly on p and T . Let us show one common approach for solving the resulting
equations, which has good numerical properties. Substituting yi = Ki xi into the
mass balance (7.49) gives F zi = Lxi + V Ki xi , and solving with respect to xi gives
xi = (F zi /(L + V Ki ). Here, introduce L = F − L (total mass balance) to derive
zi
xi = V
1+ F (Ki − 1)
Here, we cannot directly calculate xi because the vapor split V /F is not known. To
find V /F we may use the relationship Σi xi = 1 or alternatively Σi yi = Σi Ki xi = 1.
However, it has been found that the combination Σi (yi −xi ) = 0 results in an equation
with good numerical properties; this is the so-called Rachford-Rice flash equation5
zi (Ki − 1)
Σi =0 (7.50)
1 + VF (Ki − 1)
which is a monotonic function in V /F and is thus easy to solve numerically. A physical
solution must satisfy 0 ≤ V /F ≤ 1. If we assume that Raoult’s holds, then Ki depends
5 Rachford, H.H. and Rice, J.D.: “Procedure for Use of Electrical Digital Computers in Calculating
Flash Vaporization Hydrocarbon Equilibrium,” Journal of Petroleum Technology, Sec. 1, p. 19,
Oct. 1952.
194 CHEMICAL AND ENERGY PROCESS ENGINEERING

on p and T only: Ki = psati (T )/p. Then, with T and p specified, we know Ki and the
Rachford-Rice equation (7.50) can be solved for V /F . For non-ideal cases, Ki depends
also on xi and yi , so one approach is add an outer iteration loop on Ki .

Example 7.22 pT -flash. A feed F is split into a vapor product V and a liquid product L
in a flash tank (see Figure 7.4 on page 189). The feed is 50% pentane, 30% hexane and 20%
cyclohexane (all in mol-%). In the tank, T = 390K and p = 5 bar. For example, we may have
a heat exchanger that keeps constant temperature and a valve on the vapor product stream
that keeps constant pressure. We want to find the product split and product compositions.
Assume ideal liquid mixture and ideal gas (Raoult’s law).
Comment. This is a quite close-boiling mixture and we have already found that at 5 bar the
bubble point temperature is 382.64 K (Example 7.18) and the dew point temperature is 393.30
K (Example 7.20). The temperature in the flash tank must be between these temperatures for
a two-phase solution to exist (which it does in our case since T = 390 K).
Solution. The feed mixture of pentane (1), hexane (2) and cyclohexane (3) is

z1 = 0.5; z2 = 0.3; z3 = 0.2

We have Ki = psat
i (T )/p and at T = 390K and p= 5 bar, we find with the Antoine parameters
in Table 7.2:
K1 = 1.685, K2 = 0.742, K3 = 0.532
Now, zi and Ki are known, and the Rachford-Rice equation (7.50) is solved numerically to
find the vapor split V /F = 0.6915. The resulting liquid and vapor compositions are (for details
see the MATLAB code below):

x1 = 0.3393, x2 = 0.3651, x3 = 0.2956

y1 = 0.5717, y2 = 0.2709, y3 = 0.1574

% MATLAB:
z1=0.5; z2=0.3; z3=0.2; p=5; T=390;
psat1=10^(A1-B1/(T+C1)); psat2=10^(A2-B2/(T+C2)); psat3=10^(A3-B3/(T+C3));
K1=psat1/p; K2=psat2/p; K3=psat3/p; k1=1/(K1-1); k2=1/(K2-1); k3=1/(K3-1);
% Solve Rachford-Rice equation numerically to find a=V/F:
a=fzero(@(a) z1/(k1+a) + z2/(k2+a) + z3/(k3+a) , 0.5)
x1=z1/(1+a*(K1-1)), x2=z2/(1+a*(K2-1)), x3=z3/(1+a*(K3-1))
y1=K1*x1, y2=K2*x2, y3=K3*x3

Example 7.23 Condenser and flash drum for ammonia synthesis. The exit gas from
an ammonia reactor is at 250 bar and contains 61.5% H2 , 20.5% N2 and 18% N H3 . The gas
is cooled to 25o C (partly condensed), and is then separated in a flash drum into a recycled
vapor stream V and a liquid product L containing most of the ammonia. We want to calculate
the product compositions (L and V ) from the flash drum.
Data. In spite of the high pressure, we assume for simplicity ideal gas. Use vapor pressure
data for ammonia from Table 7.2 and Henry’s law coefficients for N2 and H2 from page 187.
For ammonia, we assume ideal liquid mixture, i.e., γNH3 = 1 (which is reasonable since the
liquid phase is almost pure ammonia).
Solution. The feed mixture of H2 (1), N2 (2) and N H3 (3) is

z1 = 0.615, z2 = 0.205, z3 = 0.18

For ammonia, we have at T = 298.15 K and p = 250 bar (Raoult’s law):

psat
3 (T ) 9.83 bar
K3 = = = 0.0393
p 250 bar
ENTROPY AND EQUILIBRIUM 195

For H2 and N2 , we have from the given data for Henry’s coefficient at 25o C (298.15 K):
H1 (T ) 15200 bar
K1 = = = 60.8
p 250 bar
H2 (T ) 8900 bar
K2 = = = 35.6
p 250 bar
Now, zi and Ki are known, and the Rachford-Rice equation (7.50) is solved numerically to
find the vapor split V /F = 0.8500. The resulting liquid and vapor compositions of the products
are
x1 = 0.0119, x2 = 0.0067, x3 = 0.9814
y1 = 0.7214, y2 = 0.2400, y3 = 0.0386
This agrees well with flow sheet data from a commercial ammonia plant.

Other flashes
For other flashes, like the pH-flash (which is relevant for an adiabatic flash tank), one
must include also the energy balance. For example, for an adiabatic flash tank, the
steady-state energy balance gives that the enthalpy H is constant. That is, Hin = Hout ,
and we get
F hF = V hV + LhL (7.51)
! "# $
H
where hV and hL [kJ/mol; kJ/kg] depend primarily on T , but in general also on xi , yi
and p. One solution approach is to use the pT -flash described above, and iterate on T
in an outer loop until the requirement on H is satisfied. Another approach is to solve
the equations simultaneously, as shown for the dynamic adiabatic flash of methanol
and ethanol in Example 11.18 (page 317).

7.5.4 Flash exercises

Exercise 7.8 ∗ Bubble and dew point at given temperature. A hydrocarbon mixture
contains 10% propane, 80% hexane and 10% dodecane. (a) Find the bubble point pressure at
300 K. (b) Find the dew point pressure at 300 K.

Exercise 7.9 ∗ Bubble and dew point at given pressure. A hydrocarbon mixture
contains 10 mol-% propane, 80% hexane and 10% dodecane. (a) Find the bubble point
temperature at 1 bar. (b) Find the dew point temperature at 1 bar.

Exercise 7.10 Bubble point at given pressure. A liquid mixture contains 4 mol-%
hexane and the rest is octane. What is the composition of the first vapor formed if the total
pressure is 1 atm?

Exercise 7.11 ∗ Flash at given p and T . A feed to a flash tank is 100 mol/s and contains
10% propane, 80% hexane and 10% dodecane. Find the amount of vapor product and the
compositions when T = 350K and p = 2bar.

Exercise 7.12 Flash calculation for binary mixture. Calculate the amount of liquid
that will remain at equilibrium when a mixture of 7 kg hexane and 3 kg toluene is vaporized
at 95o C and 1.5 bar.
Data: Molecular weights are 86.17 and 92.13.
196 CHEMICAL AND ENERGY PROCESS ENGINEERING

Exercise 7.13 ∗ Bubble and dew point calculations. (a) A gas mixture of 15 mol-%
benzene, 5 mol-% toluene and the rest nitrogen is compressed isothermally at 100o C until
condensation occurs. What will be the composition of the initial condensate?
(b) Calculate the temperature and composition of a vapor in equilibrium with a liquid that
is 25 mol-% benzene and 75 mol-% toluene at 1 atm. Is this a bubble point or a dew point?
(c) Calculate the temperature and composition of a liquid in equilibrium with a gas
mixture containing 15 mol-% benzene, 25 mol-% toluene and the rest nitrogen (which may be
considered non-condensable) at 1 atm. Is this a bubble point or a dew point?

Exercise 7.14 Condenser for exhaust gas. The exhaust gas from a natural gas power
plant is at 1 bar and contains 76% N2 (1), 12% O2 (2), 4% CO2 (3) and 8% H2 O (4). The
gas is cooled to 25o C (partly condensed), and is then separated in a flash drum into a gas
product V and a liquid product L containing most of the water. Find the compositions of the
product streams. Are we able to remove any significant amount of CO2 in the water?
Data: Use pure component vapor pressure data for water and Henry’s law coefficients for
the gas components (see page 187).