C04 10/04/2010 Page 158

158 Chapter 4 Single Equilibrium Stages and Flash Calculations

EXAMPLE 4.8 Liquid–Liquid Equilibrium for a §4.7 SOLID–LIQUID SYSTEMS

Four-Component Mixture.
Solid–liquid separations include leaching, crystallization,
An azeotropic mixture of isopropanol, acetone, and water is dehy- and adsorption. In leaching (solid–liquid extraction), a multi-
drated with ethyl acetate in a system of two distillation columns. component solid mixture is separated by contacting the solid
Benzene was previously used as the dehydrating agent, but legisla- with a solvent that selectively dissolves some of the solid
tion has made benzene undesirable because it is carcinogenic. Ethyl species. Although this operation is quite similar to liquid–
acetate is far less toxic. The overhead vapor from the first column, liquid extraction, leaching is a much more difficult operation
with the composition below, at 20 psia and 80 C, is condensed and
in practice in that diffusion in solids is very slow compared to
cooled to 35 C, without significant pressure drop, causing the for-
mation of two liquid phases assumed to be in equilibrium. Estimate
diffusion in liquids, thus making it difficult to achieve equili-
the amounts of the phases in kg/h and the equilibrium phase compo- brium. Also, it is impossible to completely separate a solid
sitions in wt%. phase from a liquid phase. A solids-free liquid phase can be
obtained, but the solids will always be accompanied by some
liquid. In comparison, complete separation of two liquid
Component kg/h
phases is fairly easy to achieve.
Isopropanol 4,250 Crystallization or precipitation of a component from a liq-
Acetone 850 uid mixture is an operation in which equilibrium can be
Water 2,300
achieved, but a sharp phase separation is again impossible. A
Ethyl acetate
drying step is always needed because crystals occlude liquid.
43,700
A third application of solid–liquid systems, adsorption, in-
volves use of a porous solid agent that does not undergo
phase or composition change. Instead, it selectively adsorbs
Note that the specification of this problem conforms
liquid species, on its exterior and interior surfaces. Adsorbed
with the degrees of freedom predicted by (4-4), which for
C ¼ 4 is 9. species are then desorbed and the solid adsorbing agent is
regenerated for repeated use. Variations of adsorption include
ion exchange and chromatography. A solid–liquid system is
Solution also utilized in membrane-separation operations, where the
This example was solved with the CHEMCAD program using the solid is a membrane that selectively absorbs and transports
UNIFAC method to estimate liquid-phase activity coefficients. The selected species.
results are: Solid–liquid separation processes, such as leaching and cry-
stallization, almost always involve phase-separation operations
such as gravity sedimentation, filtration, and centrifugation.
Weight Fraction

Component Organic-Rich Phase Water-Rich Phase §4.7.1 Leaching

Isopropanol 0.0843 0.0615 In Figure 4.19, the solid feed consists of particles of compo-
Acetone 0.0169 0.0115 nents A and B. The solvent, C, selectively dissolves B. Over-
Water 0.0019 0.8888 flow from the stage is a solids-free solvent C and dissolved B.
Ethyl acetate 0.8969 0.0382 The underflow is a slurry of liquid and solid A. In an ideal
1.0000 1.0000 leaching stage, all of the solute is dissolved by the solvent,
Flow rate, kg/h 48,617 2,483 whereas A is not dissolved. Also, the composition of the re-
tained liquid phase in the underflow slurry is identical to the
composition of the liquid overflow, and that overflow is free
It is of interest to compare the distribution coefficients from the of solids. The mass ratio of solid to liquid in the underflow
UNIFAC method to values given in Perry’s Handbook [1]: depends on the properties of the phases and the type of

Distribution Coefficient (wt% Basis)

Solid feed, F Overflow, V
Component UNIFAC Perry’s Handbook Insoluble A Liquid
Solute B B, C
Isopropanol 1.37 1.205 (20 C) Liquid
Acetone 1.47 1.50 (30 C) solvent, S
C MIXER-SETTLER
Water 0.0021 —
Ethyl acetate 23.5 —

Underflow, U
Results for isopropanol and acetone are in agreement at these dilute Liquid Solid
conditions, considering the temperature differences. B, C A
Figure 4.19 Leaching stage.
C04 10/04/2010 Page 159

§4.7 Solid–Liquid Systems 159

1.0
F constant-solution underflow. Figure 4.20b depicts ideal
leaching conditions when XA varies with XB. This is varia-
xB, Mass solute/mass of liquid

0.8
ble-solution underflow. In both cases, the assumptions are:
Overflow
(1) an entering feed, F, free of solvent such that XB ¼ 1; (2) a

Underflow
0.6
solids-free and solute-free solvent, S, such that YA ¼ 0 and YB
¼ 0; and (3) equilibrium between exiting liquid solutions in
0.4
Tie line underflow, U, and overflow, V, such that XB ¼ YB; and (4) a
V M U solids-free overflow, V, such that YA ¼ 0.
0.2 A mixing point, M, can be defined for (F þ S), equal to
that for the sum of the products of the leaching stage, (U þ
0 V). Typical mixing points, and inlet and outlet compositions,
S xA, Mass of solid/mass of liquid 1.0 are included in Figures 4.20a and b. In both cases, as shown
(a) in the next example, the inverse-lever-arm rule can be applied
to line UMV to obtain flow rates of U and V.
F
1.0
EXAMPLE 4.9 Leaching of Soybeans to Recover Oil.
xB, Mass solute/mass of liquid

Overflow

0.8
Soybeans are a predominant oilseed crop, followed by cottonseed,
peanuts, and sunflower seed. While soybeans are not consumed
Un

0.6 directly by humans, they can be processed to produce valuable prod-

de
rfl

ucts. Production of soybeans in the United States began after World

o

Tie line
w

0.4 War II, increasing in recent years to more than 140 billion lb/yr.
V M U Most soybeans are converted to soy oil and vitamins like niacin and
0.2 lecithin for humans, and defatted meal for livestock. Compared to
other vegetable oils, soy oil is more economical and healthier. Typi-
0 cally, 100 pounds of soybeans yields 18 lb of soy oil and 79 lb of
S xA, Mass of solid/mass of liquid 1.0 defatted meal.
To recover their oil, soybeans are first cleaned, cracked to loosen
(b) the seeds from the hulls, dehulled, and dried to 10–11% moisture.
Figure 4.20 Underflow–overflow conditions for ideal leaching: Before leaching, the soybeans are flaked to increase the mass-trans-
(a) constant-solution underflow; (b) variable-solution underflow. fer rate of the oil out of the bean. They are leached with hexane to
recover the oil. Following leaching, the hexane overflow is separated
from the soy oil and recovered for recycle by evaporation, while the
equipment, and is best determined from experience or tests underflow is treated to remove residual hexane, and toasted with hot
with prototype equipment. In general, if the viscosity of the air to produce defatted meal. Modern soybean extraction plants
liquid phases increases with increasing solute concentration, crush up to 3,000 tons of soybeans per day.
the mass ratio of solid to liquid in the underflow decreases Oil is to be leached from 100,000 kg/h of soybean flakes, con-
because the solid retains more liquid. taining 19 wt% oil, in a single equilibrium stage by 100,000 kg/h of
Ideal leaching calculations can be done algebraically or a hexane solvent. Experimental data indicate that the oil content of
with diagrams like Figure 4.20. Let: the flakes will be reduced to 0.5 wt%. For the type of equipment to
be used, the expected contents of the underflows is as follows:
F ¼ total mass flow rate of feed to be leached
b, Mass fraction of 0.68 0.67 0.65 0.62 0.58 0.53
S ¼ total mass flow rate of entering solvent
solids in underflow
U ¼ total mass flow rate of the underflow, including solids Mass ratio of solute in 0.0 0.2 0.4 0.6 0.8 1.0
V ¼ total mass flow rate of the overflow underflow liquid, XB
XA ¼ mass ratio of insoluble solid A to (solute B þ solvent
C) in the feed flow, F, or underflow, U Calculate, both graphically and analytically, compositions and
YA ¼ mass ratio of insoluble solid A to (solute B þ solvent flow rates of the underflow and overflow, assuming an ideal leaching
stage. What % of oil in the feed is recovered?
C) in the entering solvent flow, S, or overflow, V
XB ¼ mass ratio of solute B to (solute B þ solvent C) in the Solution
feed flow, F, or underflow, U
The flakes contain (0.19)(100,000) ¼ 19,000 kg/h of oil and
YB ¼ mass ratio of solute B to (solute B þ solvent C) in the
(100,000  19,000) ¼ 81,000 kg/h of insolubles. However, all of
solvent flow, S, or overflow, V the oil is not leached. For convenience in the calculations, lump the
Figure 4.20a depicts ideal leaching conditions where, in the unleached oil with the insolubles to give an effective A. The flow
underflow, the mass ratio of insoluble solid to liquid, XA, is a rate of unleached oil ¼ (81,000)(0.5/99.5) ¼ 407 kg/h. Therefore,
constant, independent of the concentration, XB, of solute in the flow rate of A is taken as (81,000 þ 407) ¼ 81,407 kg/h and the
oil in the feed is just the amount leached, or (19,000  407) ¼
the solids-free liquid. The resulting tie line is vertical. This is
C04 10/04/2010 Page 160

160 Chapter 4 Single Equilibrium Stages and Flash Calculations

18,593 kg/h of B. Therefore, in the feed, F, YA ¼ (81,407/18,593) ¼ Oil flow rate in the feed is 19,000 kg/h. The oil flow rate in the
4.38, and XB ¼ 1.0. overflow ¼ YBV ¼ 0.157(79,000) ¼ 12,400 kg/h. Thus, the oil in the
The sum of the liquid solutions in the underflow and overflow feed that is recovered in the overflow ¼ 12,400/19,000 ¼ 0.653 or
includes 100,000 kg/h of hexane and 18,593 kg/h of leached oil. 65.3%. Adding washing stages, as described in §5.2, can increase
Therefore, for the underflow and overflow, XB ¼ YB ¼ [18,593/ the oil recovery.
(100,000 þ 18,593)] ¼ 0.157.
This is a case of variable-solution underflow. Using data in the Algebraic Method
above table, convert values of b to values of XA, As with the graphical method, XB ¼ 0.157, giving a value from
kg=h A bU b the previous table of XA ¼ 2.05. Then, since the flow rate of solids in
XA ¼ ¼ ¼ ð1Þ the underflow ¼ 81,407 kg/h, the flow rate of liquid in the underflow
kg=h ðB þ CÞ ð1  bÞU ð1  bÞ
¼ 81,407/2.05 ¼ 39,711 kg/h. The total flow rate of underflow is U
Using (1), the following values of XA are computed from the pre- ¼ 81,407 þ 39,711 ¼ 121,118 kg/h. By mass balance, the flow rate
vious table. of overflow ¼ 200,000  121,118 ¼ 78,882 kg/h. These values are
close to those obtained graphically. The percentage recovery of oil,
XA 2.13 2.03 1.86 1.63 1.38 1.13 and the underflow and overflow, are computed as before.
XB 0.0 0.2 0.4 0.6 0.8 1.0

§4.7.2 Crystallization
Graphical Method
Figure 4.21 is a plot of XA as a function of XB. Because no solids Crystallization takes place from aqueous or nonaqueous solu-
leave in the overflow, that line is horizontal at XA ¼ 0. Plotted are the tions. Consider a binary mixture of two organic chemicals
feeds, F, and hexane, S, with a straight line between them. A point for such as naphthalene and benzene, whose solid–liquid equili-
the overflow, V, is plotted at XA ¼ 0 and, from above, XB ¼ 0.157. brium diagram at 1 atm is shown in Figure 4.22. Points A and
Since YB ¼ XB ¼ 0.157, the value of XA in the underflow is at the B are melting (freezing) points of pure benzene (5.5 C) and
intersection of a vertical line from overflow, V, to the underflow line.
pure naphthalene (80.2 C). When benzene is dissolved in liq-
This value is XA ¼ 2.05. Lines FS and UV intersect at point M.
uid naphthalene or vice versa, the freezing point is depressed.
In the overflow, from XB ¼ 0.157, mass fractions of solute B and
solvent C are, respectively, 0.157 and (1  0.157) ¼ 0.843. In the Point E is the eutectic point, corresponding to a eutectic tem-
underflow, using XA ¼ 2.05 and XB ¼ 0.157, mass fractions of solids perature (3 C) and composition (80 wt% benzene). ‘‘Eutec-
B and C are [2.05/(1 þ 2.05)] ¼ 0.672, 0.157(1  0.672) ¼ 0.0515, tic’’ is derived from a Greek word meaning ‘‘easily fused,’’
and (1  0.672  0.0515) ¼ 0.2765, respectively. and in Figure 4.22 it represents the binary mixture of naph-
thalene and benzene with the lowest freezing (melting) point.
5.0
Points located above the curve AEB correspond to a ho-
4.5
F mogeneous liquid phase. Curve AE is the solubility curve for
benzene in naphthalene. For example, at 0 C solubility is
XA, Mass of solid/mass of liquid

4.0
very high, 87 wt% benzene. Curve EB is the solubility curve
3.5
for naphthalene in benzene. At 25 C, solubility is 41 wt%
3.0 naphthalene and at 50 C, it is much higher. For most mix-
2.5
tures, solubility increases with temperature.
U
2.0 Under 80 B
flow
1.5
70
1.0
M 60
0.5
S V Overflow 50
0
Temperature, °C

0 0.2 0.4 0.6 0.8 1 P

XB, Mass solute/mass of liquid 40 Homogeneous solution
F
Figure 4.21 Constructions for Example 4.9.
30 Solid C10H8 + solution

G I
The inverse-lever-arm rule is used to compute the underflow and 20 Solid C6H6 H
+ solution
overflow. The rule applies only to the liquid phases in the two exit-
10
ing streams because Figure 4.21 is on a solids-free basis. The mass A
ratio of liquid flow rate in the underflow to liquid flow rate in the
0
overflow is the ratio of line MV to line MU. With M located at XA C J D
E Mixture solid C6H6 and solid C10H8
¼ 0.69, this ratio ¼ (0.69  0.0)/(2.05  0.69) ¼ 0.51. Thus, the –10
liquid flow rate in the underflow ¼ (100,000 þ 18,593)(0.51)/(1 þ 0 10 20 30 40 50 60 70 80 90 100
0.51) ¼ 40,054 kg/h. Adding the flow rates of carrier and unex- Weight percent C10H8 in solution
tracted oil gives U ¼ 40,054 þ 81,407 ¼ 121,461 kg/h or, say, Figure 4.22 Solubility of naphthalene in benzene.
121,000 kg/h. The overflow rate ¼ V ¼ 200,000  121,000 ¼
[Adapted from O.A. Hougen, K.M. Watson, and R.A. Ragatz, Chemical Pro-
79,000 kg/h. cess Principles. Part I, 2nd ed., John Wiley & Sons, New York (1954).]
C04 10/04/2010 Page 161

§4.7 Solid–Liquid Systems 161

If a liquid solution represented by point P is cooled along Crystallization of a salt from an aqueous solution can
the vertical dashed line, it remains liquid until the line inter- be complicated by the formation of water hydrates. These
sects the solubility curve at point F. are stable, solid compounds that exist within certain
If the temperature is lowered further, crystals of naphtha- temperature ranges. For example, MgSO4 forms the stable
lene form and the remaining liquid, the mother liquor, hydrates MgSO412H2O, MgSO47H2O, MgSO46H2O, and
becomes richer in benzene. When point G is reached, pure MgSO4H2O. The high hydrate exists at low temperatures; the
naphthalene crystals and a mother liquor, point H on solubil- low hydrate exists at higher temperatures.
ity curve EB, coexist, the solution composition being 37 wt% A simpler example is that of Na2SO4 and water. As seen in
naphthalene. By the Gibbs phase rule, (4-1), with C ¼ 2 and the phase diagram in Figure 4.23, only one stable hydrate is
P ¼ 2, F ¼ 2. Thus for fixed T and P, compositions are formed, Na2SO410H2O, known as Glauber’s salt. Since the
fixed. The fraction of solution crystallized can be determined molecular weights are 142.05 for Na2SO4 and 18.016 for H2O,
by the inverse-lever-arm rule. In Figure 4.22, the fraction is the weight percent Na2SO4 in the decahydrate is 44.1, which is
kg naphthalene crystals/kg original solution ¼ length of line the vertical line BFG.
GH/length of line HI ¼ (52  37)/(100  37) ¼ 0.238. The water freezing point, 0 C, is at A, but the melting
As the temperature is lowered, line CED, corresponding to point of Na2SO4, 884 C, is not on the diagram. The decahy-
the eutectic temperature, is reached at point J, where the two- drate melts at 32.4 C, point B, to form solid Na2SO4 and a
phase system consists of naphthalene crystals and a mother mother liquor, point C, of 32.5 wt% Na2SO4. As Na2SO4 dis-
liquor of eutectic composition E. Any further removal of heat solves in water, the freezing point is depressed slightly along
causes the eutectic solution to solidify. curve AE until the eutectic, point E, is reached. Curves EC
and CD represent solubilities of decahydrate crystals and
anhydrous sodium sulfate in water. The solubility of Na2SO4
EXAMPLE 4.10 Crystallization of Naphthalene from a decreases slightly with increasing temperature, which is
Solution with Benzene. unusual. In the region below GFBHI, a solid solution of anhy-
Eight thousand kg/h of a solution of 80 wt% naphthalene and 20 wt
drous and decahydrate forms exist. The amounts of coexisting
% benzene at 70 C is cooled to 30 C to form naphthalene crystals. phases can be found by the inverse-lever-arm rule.
If equilibrium is achieved, determine the kg of crystals formed and
the composition in wt% of the mother liquor.
EXAMPLE 4.11 Crystallization of Na2SO4 from Water.
Solution
A 30 wt% aqueous Na2SO4 solution of 5,000 lb/h enters a cooling-
From Figure 4.22, at 30 C, the solubility of naphthalene is 45 wt%. type crystallizer at 50 C. At what temperature will crystallization
By the inverse-lever-arm rule, for an original 80 wt% solution, begin? Will the crystals be decahydrate or the anhydrous form? At
kg naphthalene crystals ð80  45Þ what temperature will the mixture crystallize 50% of the Na2SO4?
¼ ¼ 0:636
kg original mixture ð100  45Þ
Solution
The flow rate of crystals ¼ 0.636 (8,000) ¼ 5,090 kg/h.
The remaining 2,910 kg/h of mother liquor is 55 wt% benzene. From Figure 4.23, the 30 wt% Na2SO4 solution at 50 C corresponds
to a point in the homogeneous liquid solution region. If a vertical

60
D

50
Homogeneous Solid Na2SO4 + solution
solution

40 C
B H
Temperature, °C

30
Ice + solution

20

Solids – Na2SO4 + Na2SO4⋅ 10H2O

Na2SO4⋅ 10H2O
10 + solution
A
0 E
Na2SO4⋅ 10H2O
F
Ice +
eutectic + eutectic
G I
–10
0 10 20 30 40 50 60 70 80 90 100
Weight percent Na2SO4

Figure 4.23 Solubility of sodium sulfate in water.

[Adapted from O.A. Hougen, K.M. Watson, and R.A. Ragatz, Chemical Process Principles. Part I, 2nd ed., John Wiley & Sons, New York (1954).]
C04 10/04/2010 Page 162

162 Chapter 4 Single Equilibrium Stages and Flash Calculations

line is dropped from that point, it intersects solubility curve EC at Solid adsorbent, C,
of mass amount S
31 C. Below this temperature, the crystals are the decahydrate.
The feed contains (0.30)(5,000) ¼ 1,500 lb/h of Na2SO4 and Liquid, Q
(5,000  1,500) ¼ 3,500 lb/h of H2O. Thus, (0.5)(1,500) ¼ 750 lb/h cB
are to be crystallized. The decahydrate crystals include water of hy- Equilibrium
Liquid mixture Solid, S
dration in an amount given by a ratio of molecular weights or q*B
  Carrier, A
ð10Þð18:016Þ Solute, B, of concentration cB,
750 ¼ 950 lb=h of total volume amount Q
ð142:05Þ
Figure 4.25 Equilibrium stage for liquid adsorption.
The total amount of decahydrate is 750 þ 950 ¼ 1,700 lb/h. The
water remaining in the mother liquor is 3,500  950 ¼ 2,550 lb/h.
The composition of the mother liquor is 750/(2,550 þ 750) (100%) mol/g of carbon or 0.203 g phenol/g of carbon). Thus, the
¼ 22.7 wt% Na4SO4. From Figure 4.23, the temperature corre- affinity of this adsorbent for phenol is high. The extent of
sponding to 22.7 wt% Na2SO4 on the solubility curve EC is 26 C. adsorption depends on the process used to produce the acti-
vated carbon. Adsorption isotherms can be used to determine
the amount of adsorbent required to selectively remove a
§4.7.3 Liquid Adsorption given amount of solute from a liquid.
When a liquid contacts a microporous solid, adsorption takes Consider the ideal, single-stage adsorption process of Fig-
place on the external and internal solid surfaces until equil- ure 4.25, where A is the carrier liquid, B is the solute, and C
brium is reached. The solid adsorbent is essentially insoluble is the solid adsorbent. Let: cB ¼ concentration of solute in the
in the liquid. The component(s) adsorbed are called solutes carrier liquid, mol/unit volume; qB ¼ concentration of
when in the liquid and adsorbates upon adsorption. The adsorbate, mol/unit mass of adsorbent; Q ¼ volume of liquid
higher the concentration of solute, the higher the adsorbate (assumed to remain constant during adsorption); and S ¼
concentration on the adsorbent. Component(s) of the liquid mass of adsorbent (solute-free basis).
other than the solute(s) are called the solvent or carrier and A solute material balance, assuming that the entering
are assumed not to adsorb. adsorbent is free of solute and that equilibrium is achieved,
No theory for predicting adsorption-equilibrium curves, as designated by the asterisk superscript on q, gives
ðFÞ
based on molecular properties of the solute and solid, is uni- c B Q ¼ c B Q ¼ qB S ð4-28Þ
versally embraced, so laboratory measurements must be per- This equation can be rearranged in the form of a straight line
formed to provide data for plotting curves, called adsorption that can be plotted on a graph of the type in Figure 4.24 to
isotherms. Figure 4.24, taken from the data of Fritz and obtain a graphical solution for cB and qB . Solving (4-28) for qB,
Schuluender [13], is an isotherm for the adsorption of phenol
Q ðFÞ Q
from an aqueous solution onto activated carbon at 20 C. q B ¼  c B þ cB ð4-29Þ
Activated carbon is a microcrystalline, nongraphitic form of S S
ðFÞ
carbon, whose microporous structure gives it a high internal The intercept on the cB axis is cB Q=S, and the slope is
surface area per unit mass of carbon, and therefore a high (Q=S). The intersection of (4-29) with the adsorption
capacity for adsorption. Activated carbon preferentially isotherm is the equilibrium condition cB and qB .
adsorbs organic compounds when contacted with water con- Alternatively, an algebraic solution can be obtained if the
taining dissolved organics. adsorption isotherm for equilibrium-liquid adsorption of a
As shown in Figure 4.24, as the concentration of phenol in species i can be fitted to an equation. For example, the
water increases, adsorption increases rapidly at first, then Freundlich equation discussed in Chapter 15 is of the form
increases slowly. When the concentration of phenol is 1.0 ð1=nÞ
qi ¼ Aci ð4-30Þ
mmol/L (0.001 mol/L of aqueous solution or 0.000001 mol/g
of aqueous solution), the concentration of phenol on the acti- where A and n depend on the solute, carrier, and adsorbent.
vated carbon is somewhat more than 2.16 mmol/g (0.00216 Constant, n, is greater than 1, and A is a function of tempera-
ture. Freundlich developed his equation from data on the
4 adsorption on charcoal of organic solutes from aqueous
solutions. Substitution of (4-30) into (4-29) gives
mmole
Adsorption, q*, gram

3 ð1=nÞ Q ðFÞ Q
AcB ¼  cB þ cB ð4-31Þ
S S
2
which is a nonlinear equation in cB that is solved numerically
by an iterative method, as illustrated in the following example.
1

0 1 2 3 4 5 6 7 EXAMPLE 4.12 Adsorption of Phenol on Activated

Equilibrium concentration, c,
mmole Carbon.
liter
Figure 4.24 Adsorption isotherm for phenol from an aqueous A 1.0-liter solution of 0.010 mol of phenol in water is brought to
solution in the presence of activated carbon at 20 C. equilibrium at 20 C with 5 g of activated carbon having the
C04 10/04/2010 Page 163

§4.8 Gas–Liquid Systems 163

adsorption isotherm shown in Figure 4.24. Determine the percent temperature is above the critical temperatures of most or all
adsorption and equilibrium concentration of phenol on carbon by of the species. Thus, in gas–liquid systems, the components
(a) a graphical method, and (b) a numerical algebraic method. For of the gas are not easily condensed.
the latter case, the curve of Figure 4.24 is fitted with the Freundlich Even when components of a gas mixture are at a tempera-
equation, (4-30), giving ture above critical, they dissolve in a liquid solvent to an
ð1=4:35Þ
qB ¼ 2:16cB ð1Þ extent that depends on temperature and their partial pressure
in the gas mixture. With good mixing, equilibrium between
the two phases can be achieved in a short time unless the liq-
Solution
uid is very viscous.
From the data, cB(F) ¼ 10 mmol/L, Q ¼ 1 L, and S ¼ 5 g. No widely accepted theory for gas–liquid solubilities
exists. Instead, plots of experimental data, or empirical corre-
(a) Graphical method.
From (4-29), qB ¼ ð15ÞcB þ 10ð15Þ ¼ 0:2cB þ 2.
lations, are used. Experimental data for 13 pure gases dis-
Plot this equation, with a slope of 0.2 and an intercept of 2, on solved in water are plotted in Figure 4.26 over a range of
Figure 4.24. An intersection with the equilibrium curve will temperatures from 0 to 100 C. The ordinate is the gas mole
occur at qB ¼ 1:9 mmol/g and cB ¼ 0.57 mmol/liter. Thus, the fraction in the liquid when gas pressure is 1 atm. The curves
adsorption of phenol is of Figure 4.26 can be used to estimate the solubility in water
ðFÞ
cB  cB 10  0:57
ðFÞ
¼ ¼ 0:94 or 94%
cB 10
6
(b) Numerical algebraic method.
4
Applying Eq. (1) from the problem statement and (4-31),
NH3
2:16c0:23
B ¼ 0:2cB þ 2 ð2Þ 2

or f fcB g ¼ 2:16c0:23
B þ 0:2cB  2 ¼ 0 ð3Þ 10–1
8
This nonlinear equation for cB can be solved by an iterative numeri- 6
cal technique. For example, Newton’s method [14], applied to 4
Eq. (3), uses the iteration rule:
ðkþ1Þ ðkÞ 0 ðkÞ
cB ¼ cB  f ðkÞ fcB g=f fcB g ð4Þ 2
, mole fraction per atmosphere

where k is the iteration index. For this example, f{cB} is given by SO2
10–2
Eq. (3) and f 0 {cB} is obtained by differentiating with respect to cB: 8
0 ðkÞ
f fcB g ¼ 0:497c0:77
B þ 0:2 6

4 Cl2
A convenient initial guess for cB is 100% adsorption of
ð0Þ
phenol to give qB ¼ 2 mmol=g. Then, from (4-30), cB ¼ ðqB =AÞn ¼ Br2
4:35 2
ð2=2:16Þ ¼ 0:72 mmol=L, where the (0) superscript designates H2S
the starting guess. The Newton iteration rule of Eq. (4-1) is now Cl2
used, giving the following results: 10–3
8
H
1

6
Values of

ðkÞ ðkþ1Þ
k cB f (k){cB} f 0 (k){cB} cB H2S
4
0 0.72 0.1468 0.8400 0.545 CO2
1 0.545 –0.0122 0.9928 0.558 2
2 0.558 –0.00009 0.9793 0.558
10–4
8
These results indicate convergence to f{cB} ¼ 0 for a value of 6
C2H4
cB ¼ 0.558 after only three iterations. From Eq. (1), qB ¼
2:16ð0:558Þð1=4:35Þ ¼ 1:89 mmol=g. Numerical and graphical meth- 4
C2H6
ods are in agreement. CH4
2
CO

O2 H2
10–5
N2
§4.8 GAS–LIQUID SYSTEMS
0 10 20 30 40 50 60 70 80 90 100
Vapor–liquid systems were covered in § 4.2, 4.3, and 4.4, Temperature, °C
wherein the vapor was mostly condensable. Although the Figure 4.26 Henry’s law constant for solubility of gases in water.
terms vapor and gas are often used interchangeably, the term [Adapted from O.A. Hougen, K.M. Watson, and R.A. Ragatz, Chemical Pro-
gas often designates a mixture for which the ambient cess Principles. Part I, 2nd ed., John Wiley & Sons, New York (1954).]