2 Mass Transfer in Gas-liquid Systems

In chemical separations, we can use thermodynamic models to predict the composition in each phase at equilibrium.  For example, Raoult’s Law describes the compositions of vapor and liquid phases at equilibrium.  Mass transfer models help us understand how we can manipulate the process to reach equilibrium in a faster or more economical manner.

This chapter will briefly review fundamentals of mass transfer in gas-liquid systems because many of the separation processes in this handbook involve the movement of species between gas and liquid phases.

Diffusion-based Mass Transfer

J_{Az} = molar flux of A relative to the molar-average velocity of the mixture in the z direction

D_{AB} = mutual diffusion coefficient of A in B (=D_{BA})

c_A = molar concentration of A

N_A= molar flux of species A

N = total molar flux

x_A= mole fraction of species A

c = total molar concentration

(1.1)   \begin{equation*} J_{Az}= -D_{AB}\,\frac{dc_A}{dz} \end{equation*}

(1.2)   \begin{equation*} N_A= x_A N- cD_{AB}\, \frac{dx_A}{dz} \end{equation*}

 

Falling Liquid Film with Gaseous Solute A Diffusing into Liquid B

\Gamma = liquid flow rate per unit width of film =\rho \bar{u}_y r_H (mass length^{-1} time^{-1})

\delta = film thickness (length)

\eta = metric used in selecting the appropriate equation for calculating k_c

\mu = dynamic viscosity of the liquid (pressure time)

\delta = density of the liquid (mass volume-1)

A = available area for mass transfer = L*W (length2)

c_A = concentration of species A in liquid B (mol volume-1)

c_{Ai} = concentration of species A in liquid B at the gas/liquid interface (mol volume^{-1})

c_{A0} = concentration of species A in liquid B when the liquid enters the enters (mol volume^{-1})

\Delta c_A = log mean concentration difference driving force (mol volume^{-1})

\overline {c}_{A} = bulk concentration of species A in liquid B at any position y (mol volume^{-1})

\overline {c}_{AL} = bulk concentration of species A in liquid B at any position y = L (mol volume^{-1})

D_{AB} = diffusivity of solution A in liquid B (length^2 time^{-1})

g = gravitational constant (length time^{-2})

H_A = Henry’s Law constant for solute A in our liquid at our system temperature (volume pressure mol^{-1})

k_{c,\,{\rm ave} = average mass transfer coefficient (length time^{-1})

L = film height (length)

n_A = molar flow rate of species A (mol time^{-1})

N_{\rm Re} = Reynolds number (dimensionless)

N_{\rm Sc} = Schmidt number (dimensionless)

N_{\rm Pe,\,M} = Peclet number for mass transfer (dimensionless)

P_A = partial pressure of species A in the gas phase (pressure)

r_H = flow cross section per wetted perimeter (length)

\bar{u}_y = bulk velocity of the falling film in the y direction (length time^{-1})

W = film width (length)

Reynolds number for a falling film

(1.3)   \begin{equation*} N_{\rm Re} = \frac{\rho\bar{u}_y4r_H}{\mu} = \frac{4\Gamma}{\mu} \end{equation*}

Rate

(1.4)   \begin{equation*} n_{A} = \bar{u}_y{\delta}W(\bar{c}_{AL}-c_{A0}) \end{equation*}

(1.5)   \begin{equation*} n_{A} = k_cA{\Delta}c_A \end{equation*}

Film thickness

(1.6)   \begin{equation*} \delta = \left(\frac{3\overline {u}_y\mu}{\rho g}\right)^{1/2} = \left(\frac{3\mu\Gamma}{\rho^{2}g}\right)^{1/3} \end{equation*}

Hydraulic radius for a falling film     Flow cross section: \delta W; wetted perimeter: W

(1.7)   \begin{displaymath} {r_H} = \frac{{\delta}W}{W} ={\delta} \end{displaymath}

Schmidt number

(1.8)   \begin{displaymath} {N_{\rm Sc}} = \frac{{\mu}}{{\rho}D_{AB}} \end{displaymath}

Peclet number for mass transfer

(1.9)   \begin{displaymath} {N_{\rm Pe,M}} ={N_{\rm Re}}{N_{\rm Sc}}= \frac{4{\delta}\overline u_y}{D_{AB}}=\frac{4{\Gamma}}{{\rho}D_{AB}} \end{displaymath}

Mass transfer coefficient k_{c,{\rm ave}}

(1.10)   \begin{equation*} {\eta}=\frac{2D_{AB}L}{3{\delta}^{2}\overline u_y}=\frac{{8}/{3}}{N_{\rm Re}N_{\rm Sc}({\delta}/L)}=\frac{8/3}{({\delta}/L)N_{\rm Pe,M}} \end{equation*}

                       if \eta < 0.001,

(1.11)   \begin{equation*} k_{c,{\rm ave}}=\left(\frac{6D_{AB}{\Gamma}}{{\pi}{\delta}{\rho}L}\right)^{1/2} \end{equation*}

                       if \eta > 0.1,

(1.12)   \begin{equation*} k_{c,{\rm ave}}=\left(\overline u_y\frac{{\delta}}{L}\right)\left(0.241+5.1213{\eta}\right) \end{equation*}

                       if \eta > 1,

(1.13)   \begin{equation*} k_{c,{\rm ave}}=3.414\frac{D_{AB}}{\delta} \end{equation*}

Henry’s Law, applies at the interface

(1.14)   \begin{displaymath} c_A = \frac{P_A}{H_A} \end{displaymath}

Log mean concentration difference

(1.15)   \begin{equation*} {\Delta}c_A = (c_{A_i}-\bar{c}_A)_{\rm LM}=\frac{(c_{A_i}-c_{A_0})-(c_{A_i}-\bar{c}_{A_L})}{\ln\left[{{(c_{A_i}-c_{A_0})/(c_{A_i}-\bar{c}_{A_L})}}\right]} \end{equation*}

Example

10 mL/s of water at 25^\circC flows down a wall that is 1.0 m wide and 3.0 m high. This film is in contact with pure CO2 at 1.0 atm, 25^\circC. Find the rate of absorption of CO2 into the water (kmol/s).

  • \rho_{\rm H2O} = 998 kg/m3
  • \mu_{\rm H2O} = 8.90×10-4 Pa s
  • D_{\rm CO2,H2O} = 1.96×10-9 m2/s
  • H_{\rm CO2,H2O} = 29.41 L atm mol-1

Velocity and Concentration Boundary Layers for Laminar Horizontal Flow Across a Flat Plate

\delta = thickness (height) of the velocity boundary layer (length)

\delta_c = thickness (height) of the concentration boundary layer (length)

\mu = fluid viscosity (pressure time) or (mass length-1 time-1)

\rho = fluid density (mass volume-1)

c_A = concentration of species A in fluid B at some point (x, z) (mol volume-1)

c_{A0} = initial concentration of species A in fluid B (mol volume-1)

c_{Ai} = concentration of species A in fluid B at the interface (mol volume-1)

D_{AB} = diffusivity of solute A in fluid B (length2 time-1)

N_{\rm Re,x} = Reynolds number (dimensionless)

N_{\rm Sc} = Schmidt number (dimensionless)

u_{0} = free-stream velocity of the fluid in the x direction (length time-1)

u_{x} = fluid velocity in the x direction at some point (x, z) (length time-1)

x = horizontal distance from the leading edge of the flat plate (length)

z = vertical distance from the plate surface (length)

 

(2.1)   \begin{displaymath} \frac{\delta}{x}=\frac{4.96}{{N_{\rm Re,x}}^{0.5}} \end{displaymath}

(2.2)   \begin{displaymath} N_{\rm Re,x}=\frac{xu_0{\rho}}{\mu} \end{displaymath}

(2.3)   \begin{equation*} \frac{u_x}{u_0}=1.5\left(\frac{z}{\delta}\right)-0.5\left(\frac{z}{\delta}\right)^{3} \end{equation*}

(2.4)   \begin{equation*} {\delta}_c =\frac{\delta}{{N_{\rm Sc}^{1/3}}} \end{equation*}

(2.5)   \begin{displaymath} N_{\rm Sc}=\frac{\mu}{{\rho}D_{AB}} \end{displaymath}

(2.6)   \begin{displaymath} \frac{(c_{A_i}-c_A)}{(c_{A_i}-c_{A_0})}=1.5\left(\frac{z}{\delta_c}\right)-0.5{\left(\frac{z}{\delta_c}\right)^{3}} \end{displaymath}

Example

Air at 100^\circC and 1.0 atm with a free-stream velocity of 5.0 m/s flows over a 3.0 m long flat plate made of naphthalene.

(a) At what horizontal position does flow become turbulent?

(b) What is the thickness of the velocity boundary layer at that point?

(c) What is the thickness of the concentration boundary layer at that point?

(d) What is the concentration of naphthalene in the air at this horizontal position and vertical position that is half of the height of the concentration boundary layer?

  • \mu_{\rm air} = 2.15×10-5 kg m-1 s-1
  • \rho_{\rm air} = 0.0327 kmol m-3
  • D_{AB} = 0.94×10-5 m2 s-1
  • c_{Ai} = 4.3×10-3 kmol m-3

 

Film Theory and Overall Mass Transfer Coefficients

\delta = thickness of the film in which c_{A} \neq c_{Ab} (length)

c = total concentration of liquid B (mol volume-1)

c_{Ai} = concentration of species A in liquid B at the interface (mol volume-1)

c_{Ab} = bulk concentration of species A in liquid B (mol volume-1)

c^*_{A} = concentration of species A in liquid B at equilibrium with the bulk gas phase (mol volume-1)

D_{AB} = diffusivity of species A in liquid B (area time-1)

H'_{A} = Henry’s Law constant for equation of the form c_{A}=H'_{A}P_{A};
(mol volume-1 pressure-1)

k_{c} = liquid-phase mass transfer coefficient, with respect to concentration driving force (length time-1)

K_{G} = overall mass transfer coefficient, with respect to pressure driving force (mol time-1 area-1 pressure-1)

K_{L} = overall mass transfer coefficient, with respect to concentration driving force (length time-1)

k_{p} = gas-phase mass transfer coefficient, pressure driving force (mol time-1 area-1 pressure-1)

N_{A} = molar flux (mol area-1 time-1)

P_{Ab}=p_{Ab} = partial pressure of species A in the bulk gaseous phase (pressure)

P^*_A = partial pressure of species A in a gas at equilibrium with the bulk liquid phase (pressure) x_{Ai}   = mole fraction of species A in liquid B at the interface

x_{Ab} = bulk mole fraction of species A in liquid B

 

Film Theory

(3.1)   \begin{equation*} N_A = \left(\frac{D_{AB}}{\delta}\right)(c_{Ai}-c_{Ab}) = \left(\frac{cD_{AB}}{\delta}\right)(x_{Ai}-x_{Ab}) \end{equation*}

 

Example

SO2 is absorbed from air into water using a packed absorption tower. At a specific location in the tower, we know that the pressure of SO2 is 0.15 atm. Measurements of the gas composition above and below this location in the tower have told us that the flux of SO2 into the water is 0.0270 kmol SO2 m-2 hr-1. We were able to sample the bulk liquid phase at this location and found that it contained 3.0×10-4 mol SO2/mol. Assuming that this system fits film theory, find the thickness of the film.

  • D_{\rm SO_2,H_2O} = 1.7*10-9 m2 s-1
  • H_{\rm SO_2,H_2O} = 1.2 L atm mol-1

 

Two-Film Theory, Overall Mass Transfer Coefficients, Gas/Liquid

(3.2)   \begin{equation*} N_A = \left(P_{Ab}H'_A-c_{Ab}\right)\left({\frac{H'_A}{k_p}+\frac{1}{k_c}\right)^{-1} \end{equation*}

(3.3)   \begin{equation*} \frac{1}{K_L}=\frac{H'_A}{k_p}+\frac{1}{k_c} \end{equation*}

(3.4)   \begin{equation*} N_A=K_L(c^{*}_A-c_{Ab}) \end{equation*}

(3.5)   \begin{equation*} N_A=\left(P_{Ab}-P^{*}_A\right)\left({\frac{1}{k_p}+\frac{1}{H'_Ak_c}\right)^{-1} \end{equation*}

(3.6)   \begin{equation*} \frac{1}{K_G}=\frac{1}{k_p}+\frac{1}{H'_Ak_c} \end{equation*}

(3.7)   \begin{equation*} N_A=K_G(P_{Ab}-P^{*}_A) \end{equation*}

 

Example

We intend to use water to absorb SO2 from air. The incoming air is at 50°C and 2.0 atm and contains 0.085 mol SO2/mol. The incoming water, also at 50°C, already contains 0.0010 mol SO2/mol.

(a) Which phase is most limiting to mass transfer?

(b) What is the expected initial flux value?

  • k_c = 0.18 m hr-1
  • k_p = 0.040 kmol hr-1 m-2 kPa-1
  • H'_A (50°C) = 0.76 mol atm-1 L-1

 

 

License

Icon for the Creative Commons Attribution-NonCommercial 4.0 International License

Chemical Engineering Separations: A Handbook for Students Copyright © 2021 by Monica H. Lamm and Laura R. Jarboe is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License, except where otherwise noted.