Lecture Notes 3 - Mass Transfer 2
Lecture Notes 3 - Mass Transfer 2
Lecture Notes 3 - Mass Transfer 2
SEPARATION PROCESSES I
Diffusivities are generally higher for gases than liquids, both of which are higher
than for solids. For example, the molecules in liquids are closer together than in
gases, and hence, diffusing molecules will collide more often and diffuse more
slowly than in gases. This explains why diffusion coefficients in liquids are
relatively smaller than in gases.
Diffusivity values
Gases: 5×10‒6 to 1×10‒5 m2 /s
Liquids: 1×10‒10 to 1×10‒9 m2 /s
Solids: 1×10‒14 to 1×10‒10 m2 /s
Measurement of diffusivity can be made in the laboratory, but there are also
semi-theoretical equations (correlated from experimental data) that can be used
to predict the values of diffusivities for different systems.
1
(i) Fuller, Schettler, and Giddings Equation:
1⁄2
1.00 × 10−7 𝑇 1.75 (1⁄𝑀 + 1⁄𝑀 )
𝐴 𝐵
𝐷𝐴𝐵 = 1⁄ 3 1⁄ 3 2
𝑃𝑇 [(∑ 𝑣𝐴 ) + (∑ 𝑣𝐵 ) ]
𝑇 is absolute temperature (K); 𝑃 is total gas pressure (atm); 𝑀A and 𝑀B are molar
masses (g/mol); 𝑣𝐴, 𝑣𝐵 are atomic diffusion volumes (or structural volume
increments). Values of atomic diffusion volumes for common substances are
given in Table 18.2-2 (Geankoplis, 5th Ed.). Here, 𝐷AB is proportional to 1/𝑃 and
to 𝑇1.75.
𝑇
𝐷𝐴𝐵 𝛼 ( ), 𝜇 = viscosity
𝜇
9.96 × 10−16 𝑇
𝐷𝐴𝐵 = 1/3
𝜇𝑉𝐴
2
(ii) Wilke-Chang equation:
𝑇
𝐷𝐴𝐵 = 1.173 × 10−16 (∅𝑀𝐵 )1/2
𝜇𝐵 𝑉𝐴0.6
Solvent 𝜱𝑩
Water 2.6
Methanol 1.9
Ethanol 1.5
Others 1.0
Example 2.4
Estimate the diffusivity of methanol vapour in carbon monoxide gas at 1 atm and
100 °C, using the Fuller et al method.
Species B: CO
3
From table 18.2-2:
1
1.00 × 10−7 (373)1.75 (1/32.04 + 1/28.01)2
DAB =
1 1 2
P [(29.9)3 + (18.9)3 ]
In many cases, to have a fluid in convective flow, the fluid velocity is increased
until turbulent diffusion or eddy current occurs. Convective mass transfer from
the surface to a fluid, or vice versa, can be represented as follows:
4
The simplest model is to assume a linear profile, then an equation like Fick’s law
can be written for the flux:
𝑑𝑐𝐴
𝐽𝐴 = −(𝐷𝐴𝐵 + 𝜀𝑀 )
𝑑𝑧
where DAB is the molecular diffusivity in m2/s and εM is eddy diffusivity in m2/s.
In practice 𝑑𝑧 (𝑖. 𝑒. 𝑧1 − 𝑧1 ) is very small ((𝜇𝑚). It cannot be predicted, hence
∆z is included in the proportionality constant referred to as convective mass-
transfer coefficient. The flux, denoted by NA is thus written as:
𝑁𝐴 = 𝑘𝑐 (𝑐𝐴1 − 𝑐𝐴2 )
Because there are many measures of composition (other than 𝑐𝐴 ) there are
equally many forms of the MTC.
𝑝𝐴 𝑅𝑇
For perfect gases: 𝑝𝐴 = 𝑦𝐴 𝑃; 𝑐𝐴 = ; 𝑘𝑐 = 𝑘𝐺 𝑅𝑇 = 𝑘𝑐 and 𝑘𝑦 = 𝑘𝐺 𝑃
𝑅𝑇 𝑃
𝑐𝐴 𝑘𝑥
𝑥𝐴 = ; 𝑘𝑐 = 𝑘𝐿 =
𝑐 𝑐
Variables needed to describe the system include those descriptive of the system
geometry, the flow and fluid properties and the quantity of importance, 𝑘𝑐 .
5
Analysis by the Buckingham π-method:
There are six variables and three fundamental dimensions (M, L, T). Hence three
dimensionless groups (6 – 3 = 3) are needed to describe the system. Using the
core variables 𝐷, and 𝑑 gives the following three dimensionless groups:
𝑘𝑐 𝐿
1. Sherwood number: 𝑆ℎ =
𝐷
𝐿υ𝜌
2. Reynolds number: 𝑅𝑒 =
𝜇
𝜇
3. Schmidt number: 𝑆𝑐 =
𝜌𝐷
The three relations show that 𝑆ℎ = 𝑓(𝑅𝑒, 𝑆𝑐). The general equation is of the
form:
6
(Gases) (Liquids)
The interface:
• Assumed to have zero thickness, hence offers zero resistance to mass transfer.
• The required driving force for any finite flux is zero. (A system with zero
driving force is an equilibrium system).
7
We consider a two-film or two-resistance model, involving two mass transfer
films (gas phase and liquid phase films), each is characterized by a film mass
transfer coefficient, MTC (denoted by lower-case k).
where 𝑘y and 𝑘x are the locally applicable mass transfer coefficients for the gas
and liquid films.
At steady state, the flux of A is everywhere the same (through each film and
through the interface). Thus, the rate of transfer to the interface is equal to the
rate of transfer from the interface.
(𝑦𝐴 − 𝑦𝐴𝑖 ) 𝑘𝑥
→ =−
(𝑥𝐴 − 𝑥𝐴𝑖 ) 𝑘𝑦
The ratio −𝑘x/𝑘y is the slope of the line connecting the equilibrium composition
with the bulk composition. The equation can be represented on a graph as
follows:
8
The equilibrium concentrations, 𝑦𝐴𝑖 and 𝑥𝐴𝑖 are related by equilibrium relations.
Note that these compositions do not exist physically. They are hypothetical
mixtures that would be in equilibrium with the bulk of the other phase. To
calculate the rate of mass transfer, a more appropriate expression for 𝑁A is used
which uses an overall mass transfer coefficient, normally indicated with an
upper-case letter, 𝐾.
𝑝𝐴 = 𝑦𝐴 𝑃
where 𝑦𝐴 is the mole fraction of component A in the gas phase, and P the total
pressure.
9
Henry’s Law: For dilute solutions, the equilibrium relation between the partial
pressure of component A in the gas phase, 𝑝A, and mole fraction in the liquid
phase, 𝑥 A, is expressed by a straight-line Henry’s law equation.
𝑝𝐴 = 𝑥𝐴 𝐻
where 𝐻 is Henry’s law constant in atm/mole fraction for the given system.
𝑦𝐴 = 𝑝𝐴 ⁄𝑃 = 𝐻′ 𝑥𝐴 , 𝐻 ′ = 𝐻 ⁄𝑃
where 𝐻′ is the Henry’s law constant in mole frac gas/mole frac liquid and is
equal to 𝐻/𝑃. (𝐻′ depends on total pressure, whereas 𝐻 does not).
10
As can be seen from the graph, 𝑥𝐴 is uniquely related to 𝑦𝐴∗ for a given condition
of 𝑇 and 𝑃. From the graph:
Therefore
𝑁𝐴 𝑁𝐴 𝑁𝐴 1 1 𝑚′
= + 𝑚′ ⟹ = +
𝐾𝑦 𝑘𝑦 𝑘𝑥 𝐾𝑦 𝑘𝑦 𝑘𝑥
1 1 1
= ′′ +
𝐾𝑥 𝑚 𝑘𝑦 𝑘𝑥
For dilute solutions, with 𝑦𝐴 = 𝐻′ 𝑥𝐴∗ (i.e., assuming the equilibrium curve is
linear), the slope of the equilibrium curve is the Henry’s law constant. Hence,
1 1 𝐻′ 1 1 1
= + and = ′ +
𝐾𝑦 𝑘𝑦 𝑘𝑥 𝐾𝑥 𝐻 𝑘𝑦 𝑘𝑥
Limiting cases
• If 𝑘 y and 𝑘 x are nearly equal and the slope is so small so that solute A is
relatively insoluble in liquid, the term 𝑚’/𝑘x will be negligible compared to
1/𝑘y. Hence,
1 1
≈
𝐾𝑦 𝑘𝑦
We say the rate of mass transfer is gas-phase controlled (i.e., the overall
resistance lies in the gas phase only).
1 1
≈
𝐾𝑥 𝑘𝑥
• For cases where 𝑘y and 𝑘x are not nearly equal, then it will be the relative
size of the ratio 𝑘𝑥 ⁄𝑘𝑦 and of slope which will determine the location of the
controlling mass transfer resistance.
Example 2.4
For a system in which component (A) is transferring from the liquid to the gas
phase, the equilibrium is given by 𝑦A∗= 0.75 𝑥A. At one point in the apparatus
the liquid contains 90 mol% of (A) and the gas contains 45 mol% of (A). The
individual gas film mass transfer coefficient at this point in the apparatus is
12
0.02716 kmol/m2.s, and 70% of the overall resistance to mass transfer is known
to be encountered in the gas film. Determine:
Solution:
𝑦𝐴∗
𝑥𝐴∗ 𝑥𝐴𝑖
𝑦𝐴𝑖
𝑦𝐴 = 0.45
Distance, z
Part (a)
1 1 1 1 𝑘𝑚𝑜𝑙
= 0.7 ( ) → = 0.7 ( ) ⟹ 𝐾𝑦 = 0.019
𝑘𝑦 𝐾𝑦 0.02716 𝐾𝑦 𝑚2 . 𝑠
𝑁𝐴 = 𝐾𝑦 (𝑦𝐴∗ − 𝑦𝐴 )
Part (b)
13
Part (c)
1 1 𝐻′ 1 1 0.75
= + → = + ⟹ 𝑘𝑥 = 0.0476
𝐾𝑦 𝑘𝑦 𝑘𝑥 0.019 0.02716 𝑘𝑥
1 1 1 1 1 1
= ′ + = +
𝐾𝑥 𝐻 𝑘𝑦 𝑘𝑥 𝐾𝑥 (0.75)(0.02716) 0.0476
𝑘𝑚𝑜𝑙
⟹ 𝐾𝑥 = 0.0142
𝑚2 . 𝑠
14