Ultimate Mass Transfer Notes- CHE 120 Winter 2017

Welcome to this study guide. Hw uses matlab code. A lot of plug n

chug. Exams are more theoretical. Try to read the book(Fundamentals
of Momentum, Mass, and Heat Transfer Welty Rorrer Foster) and go
through relevant equations. Plug in whatever equation into MATLAB (or
any other calculator to solve). Careful with units. Theyre not all the
same for each equation, e.g H, henrys constant can be in atm or
needs a conversion factor such as density/Mol weight. This does not
have a lot of important equations due to not useful in matlab
such as continuity or other differential equations. Check book.
MATLAB may be a bit difficult to read but most of it should
make sense. Plug into matlab script and see all relevant variables,
equations, constants, etc. If there was a table, I ran the script and then
plugged that into a variable I define, such as phi, in the script. That
way you can use the same equations and simply change the numbers
for the given problem, but be careful with the units!
P.S shoutout to Wyman. Yeah hes boring and what not but I feel like he
actually cares.
General Equations

R= 82.06 cm^3*atm/molK

%linear interpolation:
x=7.7108; x1=7; x2=8; y1= 0.7896; y2= 0.7712;
y=y1+(x-x1)*(y2-y1)/(x2-x1);

Chp 24 Fundamentals of
Mass Transfer
%Chapter 24 Equations
MWavg=sum(yi*MWi); %avg moelcular weight= sum(yi*MWi)
rhoavg=MWavg*C; % average gas density(g/L)=MW(g/mol)*C(mol/L)
PA=yA*P
CA=PA/(R*T)
C=P/(R*T)
BK=1.3*10^-16; %boltzman constant in ergs/K
N=6.022*10^-23;

%knudson coefficient
DKA=4850*dpore*sqrt(T/MA) ; %24-58 b.

%A is H2S
%B is CH4
MA=38.08; %g/gmole
MB=16.05; %g/gmole
N=6.022*10^23;
BK=1.3*10^-16; %boltzman constant in ergs/K
sigmaA=3.623;%given
sigmaB= 3.822; %appendix K.2
sigmaAB=(sigmaA+sigmaB)/2;
MA=34; %MW of H2S
MB=16; %MW of CH4
phiA=301.1; %epsilon A over BK, given\
phiB=136.5;
epsB=BK*136.5; %appendix K.2- eps/BK = 136.5 so you multiply
eps=BK*sqrt(phiA*epsB/BK);% equation 24-45
%or alternatively
eps=BK*sqrt(phiA*phiB); % and they are the same
phi=eps/BK;
AppendixKrowI= T/phi % 1.4946 Appendix K.1
%assume linear for linearInterpolation:
x1=1.45; x2=1.5; y1=1.215; y2= 1.198; x=AppendixKrowI;
y=y1+(x-x1)*(y2-y1)/(x2-x1);
OhmD=y
%OhmD=1.1998;
%Binary Gas-Phase Diffusion Coeff- Hirschfelder Equation 24-33
DABhirschfelder=0.001858*T^(3/2)*(1/MA+1/MB)^.5/(P*sigmaAB^2*OhmD)
%cm^2/
%Fuller,Schettler, and Giddlings Equation
%A=H2S; B=CH4
%use table 24.3
SVA=2*1.98+17; %sum of A from table 24.3 H2S-> H2+S->2*1.98+17
SVB=16.5+4*1.98; %sum of A from table 24.3 CH4-> C+4H->
DABfuller=0.001*T^(1.75)*sqrt(1/MA+1/MB)/(P*(SVA^(1/3)+SVB^(1/3))^2)
DAB=DABfuller;

%solve DprimeAe
%eps=0.5;
Dpore=20*10^-7;
Dka=4850*Dpore*sqrt(T/MA)
DAe=(1/DAB+1/Dka)^-1
DprimeAe=eps^2*DAe;

muB=1.394; %viscocity water in cp

DAB=T/muB*7.4*10^-8*sqrt(phiB*MB)/VA^0.6;

HW:
24.5 Consider the following properties of the atmosphere of
the planet Mars at a particular measurement point on the surface,
as measured by the Mars Rover:
Average surface pressure: 6.1 mbar
Average temperature: 210 K
Atmospheric gas composition (by volume): carbon dioxide
(CO2) 95.32%; nitrogen (N2) 2.7%; argon (Ar) 1.6%; oxygen
(O2) 0.13%; carbon monoxide (CO) 0.08%.
a. What is the partial pressure of CO2 in the Martian atmosphere,
pA, in units of Pa?
b. What is the molar concentration of CO2 in the Martian
atmosphere, cA?
c. What is the total molar concentration of all gases in the
Martian atmosphere?
d. What is the total mass concentration of all gases the Martian
atmosphere?
%masshw1chp24p4.m
clear all
clc

P=6.1*10^-3*10^5; %mbar converted to bar converted to atm

T=210; %K
R=8.314;
yA=0.9532; %mole fraction of gas
PA=yA*P; %P in atm
CA=PA/(R*T);% gmol/cm^3
C=P/(R*T); % gmol/cm^3

SyiMWi= 43.41; %sum of yi*MWi

MWavg=SyiMWi;
rhoavg=MWavg*C; %g/m^3

24.12 Sour natural gas is contaminated with hydrogen sulfide.

The H2S vapors are commonly treated removed by passing
the gas through a packed bed of adsorbent particles. In the
present process, natural gas with a bulk gas composition is
99 mole% methane (CH4) and 1.0 mole% H2S will be treated
with a porous adsorbent material of void fraction of 0.50 and
mean pore diameter of 20 nm at 30C and 15.0 atm total system
pressure. To design the adsorption bed, it is necessary to estimate
the diffusion coefficient of H2S within the porous adsorbent.
The LennardJones parameters for H2S are s . 3.623 and
e/k . 301.1 K.
a. What is molar concentration of H2S in the gas mixture?
b. Estimate the binary gas-phase molecular diffusion
coefficient of methane in H2S by two methods, and compare
the results.
c. Estimate the effective diffusion coefficient of H 2S within the
porous material, assuming methane and hydrogen sulfide fill
the gas void space within the porous material.

%masshw1chp24p12.m
clear all
clc

%part a
yA=0.01;
P=15; %atm
R=82.06;
T=303; %K
PA=yA*P;
CA=PA/(R*T); %6.0328e-06 mol/cm^3

%part b
%A is H2S
%B is CH4
MA=38.08; %g/gmole
MB=16.05; %g/gmole
N=6.022*10^23;
BK=1.3*10^-16; %boltzman constant in ergs/K
sigmaA=3.623;%given
sigmaB= 3.822; %appendix K.2
sigmaAB=(sigmaA+sigmaB)/2;
MA=34; %MW of H2S
MB=16; %MW of CH4
phiA=301.1; %epsilon A over BK, given\
phiB=136.5;
epsB=BK*136.5; %appendix K.2- eps/BK = 136.5 so you multiply
eps=BK*sqrt(phiA*epsB/BK);% equation 24-45
%or alternatively
eps=BK*sqrt(phiA*phiB); % and they are the same
phi=eps/BK;
AppendixKrowI= T/phi % 1.4946 Appendix K.1
%assume linear for linearInterpolation:
x1=1.45; x2=1.5; y1=1.215; y2= 1.198; x=AppendixKrowI;
y=y1+(x-x1)*(y2-y1)/(x2-x1);
OhmD=y
%OhmD=1.1998;
%Binary Gas-Phase Diffusion Coeff- Hirschfelder Equation 24-33
DABhirschfelder=0.001858*T^(3/2)*(1/MA+1/MB)^.5/(P*sigmaAB^2*OhmD)
%cm^2/
%Fuller,Schettler, and Giddlings Equation
%A=H2S; B=CH4
%use table 24.3
SVA=2*1.98+17; %sum of A from table 24.3 H2S-> H2+S->2*1.98+17
SVB=16.5+4*1.98; %sum of A from table 24.3 CH4-> C+4H->
DABfuller=0.001*T^(1.75)*sqrt(1/MA+1/MB)/(P*(SVA^(1/3)+SVB^(1/3))^2)
DAB=DABfuller;

%solve DprimeAe
%eps=0.5; Im just not sure about the difference between eps and epsAB
its for porous material
Dpore=20*10^-7;
Dka=4850*Dpore*sqrt(T/MA)
DAe=(1/DAB+1/Dka)^-1
DprimeAe=eps^2*DAe;

24.18 Estimate the liquid diffusion coefficients of the following

solutes that are transferred through dilute solutions:
a. Methanol in water at 288 K
b. Water in methanol at 288 K
c. n-butanol in water at 288 K
d. Water in n-butanol at 288 K
Compare these estimated values to the available experimental
values reported in Appendix J.2.

%masshw1chp24p18.m
clear all
clc

T=288;
MB=18;% molecular of water
y=y1+(x-x1)*(y2-y1)/(x2-x1);
muB=1.394; %viscocity water in cp
phiB=2.6;%table 24.5 page 450
VA=(4*3.7+14.8+7.4);%volumes of H, C, O respectively
DAB=T/muB*7.4*10^-8*sqrt(phiB*MB)/VA^0.6

Chp 25 Differential
Equations of Mass transfer
Chp 26 Steady State
Molecular Diffusion
Ex 1
Diffusivity

%mass chp26 ex 1
clear all
clc

%equation
%NA=cDAB/delz*ln(1-ya2/(1-ya1))
%WA=NA*S
P=1;%1 atm
R=0.08206; %m^3*atm/kgmole K
T=273+35; %K
c=P/(R*T); %0.0396 kg mole /m^3

PA=(115.1/760);% mm Hg converted to atm

YA1=PA/P; %0.152;
YA2=0; %"the mole fraction of TCE vapor at the exit of the tank is taken as 0"

z2=5; %m
z1=0.2; %m

DAB=0.088/100^2; %cm^2/s converted to m^2/s

NAz=c*DAB/(z2-z1)*log((1-YA2)/(1-YA1))

D=2;
S=pi*D^2/4;
%book: NAz=1.197*10^-8
W=NAz*S*131.4*3600*24% kg/kgmole * s * h/ h / day -total kg TCE/day

W= 0.423 kg TCE / day

Midterm Review
Chp 27 Unsteady-
State Molecular
Diffusion
27.4
One step of the manufacturing of silicon solar cells is the molecular diffusion (doping) of
elemental phosphorous (P) into crystalline silicon to make an n-type semiconductor. This
Pdoped layer needs to be at least 0:467 mm into the 200-mm thick wafer. The present
diffusion process is carried out at 1000 C. Data for the total amount of phosphorous atoms
loaded into the silicon wafer vs. time at 1000 C are presented in the figure below. The
maximum solubility of phosphorous within crystalline silicon is 1:0 10 21 P atoms=cm3 at 1000
C. The square silicon wafer has a surface area of 100 cm 2 (10 cm=side). Initially, there is no
phosphorous impurity in the crystalline silicon. What is the concentration of P atoms (atoms
P/cm3) doped into the silicon at a depth of 0:467 mm after 40 min, based on the data
provided?

%masshw4ch27p4.m
clear all
clc
Cao=0;
mao=0;
S=100; %cm^2

CAS=1*10^21; %P atoms
t=40*60; %mins converted to hour
mat=6*10^18/100; %from slope in units P atoms/ sqrt(t)
DAB=pi*mat^2/(CAS^2*S^2*4) %cm^2/s
z=0.468*10^-4; %micro m conveted to cm
phi = z/(2*sqrt(DAB*t)) %from error function phi= 0.8983
%from appendix L erf(phi)=0.797
x1=0.85; x2=0.9; y1=0.7707; y2=0.797; x=phi;
y=y1+(x-x1)*(y2-y1)/(x2-x1); %y=erf(phi)=erfphi= 0.7961
erfphi= 0.7961;
CA=CAS*(1-erfphi) %atoms/cm^3

27.9 A simple experiment is set up to measure the effective diffusion coefficient of blue
dextran dye into a rigid agarose gel. An aqueous, well-mixed solution of liquid height 10 cm
containing 1:0 g=L (1:0 10_3 g=cm3) of the dye rests over the rigid gel of 2:0 cm thickness, as
shown in the figure in the right-hand column above. The solubility of the dye in both water and
the agarose gel are the samei.e., the concentration of dye on the water side of the water-gel
interface is equal to the concentration of dye on the gel side of the water-gel interface. The
aqueous gel is generally considered a water-hydrated polymer, where solute molecules diffuse
through the hydrated regions. Initially, there is no dye dissolved the gel. After 24 h, a tiny
section of the gel, 2:0 mm from the surface, is very carefully excised with a syringe needle.
The concentration of the dye within the gel sample is measured by a spectrophtometer, and
reads 0:203 g dye=L (2: 03 10_4 g=cm3). Based on the experimental measurements, estimate
the effective diffusion coefficient of blue dextran dye in the agarose gel. Justify by a calculation
that it is reasonable to assume that the concentration of the dye in the well-mixed aqueous
solution does not change significantly over the course of the experiment.

Chp 28 Convective
Mass Transfer

Mass Hw 5 Week 8
Andrew Sanchez

28.4

A thin (1.0-mm-thick) coat of fresh paint has just been

 sprayed over a 1.5-m by 1.5-m square steel body part, which
approximates a flat surface. The paint contains a volatile solvent
that initially constitutes 30 wt% of the wet paint. The initial
density of the wet paint is 1.5 g/cm3. The freshly painted part is
introduced into a drying chamber. Air is blown into the rectangular
drying chamber at a volumetric flow rate of 60 m3/min, as
shown in the figure below, which has dimensions L.1.5 m,
H.1.0 m,W.1.5 m). The temperature of the air stream and the
steel body part are both maintained at 27C, and the total system
pressure is 1.0 atm. The molecular weight of the solvent is 78 g/
gmole, the vapor pressure exerted by the solvent at 27C is
105mm Hg, and the molecular diffusion coefficient of the
solvent vapor in air at 27C and 1.0 atm is 0.097 cm2/s.

What is the Schmidt number and the average Sherwood

number (ShL) for the mass-transfer process?
b. What is the estimated solvent evaporation rate from the
surface of the whole body part in units of g/min? It may
be assumed that convection mass-transfer limits the evaporation
rate, and that the concentration of solvent vapor in the
bulk gas is finite, but can be approximated as cA; 0.
c. Using the results from part (a), how long will it take for the
paint to completely dry?
d. What are the hydrodynamic (d) and concentration boundarylayer
(dc) thicknesses at x.L.1.5 m? How does this
compare to H, the height of the drying chamber?
e. What is the new required air flow rate (m3/

a. Find Shmidt:
a. Sc= v(air)/DAB
i. Why is vair not Volumetric flow rate of air/ cross-sectional Area- vair
from appendix I- 1.5689*10^-5
ii. DAB- given above
iii. Vai
DAB= 0.097; %cm^2/s given
vair= 1.5689*10^-5; %m^2/s
Sc=vair/DAB*100^2= 1.6174
iv. Why is v0 volumetric flowrate
W=1.5;
H=1;
L=1.5;

%solve Shl
v0= 60; %m^3
vinf=v0/(W*H*60)
%need vinf for Re
%need Re for ShL
ReL= vinf*L/vair;
 ShL=0.664*ReL^.5*Sc^(1/3)=196.778

%part b solve WA

kc=ShL*DAB/(L*100)% 0.127cm/s
CA=(105/760)/(82.06*300)%=5.61*10^-6;
W=kc*CA*L^2*100^2*78*60 = 75.1987 g/min

m=1.5*1.5^2*100^2/10*.3

t=m/W

del=5*1.5/sqrt(ReL)
delc=del/Sc^(1/3)

finalans=v0*(150/W)^2

28.7 In a manufacturing process, an organic solvent (methyl

ethyl ketone, MEK) is used to dissolve a thin coating of a
polymer film away from a nonporous flat surface of length
20 cm and width 10 cm, as shown in the figure (above right).
The thickness of the polymer film is initially uniform at l a .
0:20 mm (0.02 cm). In the present process, a volumetric flow
rate of 30 cm3/s of MEK liquid solvent is added to an open flat
pan of length 30 cm and width 10 cm. The depth of the liquid
MEK solvent in the pan is maintained at 2.0 cm. It may be
assumed that the concentration of dissolved polymer in the
bulk solvent is essentially zero (cA; .0) even though in reality
the concentration of dissolved polymer in the solvent increases
very slightly from the entrance to the exit of the pan. It may
also be assumed that the change film thickness during the
dissolution process does not affect the convection masstransfer
process. Let A.polymer (solute), B.MEK (liquid
solvent).
a. What are the Schmidt number and the average Sherwood
number for the mass-transfer process?
b. What is the average flux of dissolved polymer (in units of
g A/cm2 s) from the surface?
c. Plot out the local kc;x vs. position x, for 0 < x L.
Potentially useful data: DAB . 3.0 10_6 cm2/s, the diffusion
coefficient of dissolved polymer (solute A) in liquid MEK
solvent at 25C; rA,solid.1.05 g/cm3, mass density of solid
polymer film material; cAs.cA
*.0.04 g/cm3, maximum solubility
of dissolved polymer (solute A) in MEK solvent at
25C; nB . 6.0 10_3 cm2/s, kinematic viscosity of liquid
MEK at 25C; rB.0.80 g/cm3, mass density of liquid MEK
solvent at 25C.

%masshw5 ch 28.4

clear all
clc

%28.7
%a solve Shmidt
DAB= 3*10^-6; %cm^2/s given
vair= 6*10^-3; %cm^2/s
Sc=vair/DAB
v0=30;
W=10;
L=20;

vinf=v0/(W*L)
Re=vinf*L/vair
ShL=0.664*Re^(1/2)*Sc^(1/3)
ShL=0.664*Sc^(1/3)*(v0/(W*vair))^0.5

kc=ShL*DAB/L
Ca=0.04
NA=kc*Ca

x=[0:0.1:L];
kc=ShL*DAB./x;
plot (x,kc)

28.12 A well-mixed open pond contains wastewater that is

contaminated with a dilute concentration of dissolved methylene
chloride. The pond is rectangular with dimensions of 500 m by
100 m, as shown in the figure (above right). Air at 27C and
1.0 atm blows parallel to the surface of the pond with a bulk
velocity of 7.5 m/s. At 20C and 1.0 atm, for the gas phase
(A.methylene chloride, B.air), the diffusion coefficient (DAB)
is 0.085 cm2/s, and kinematic viscosity (nB) is 0.15 cm2/s. At
27C, for the liquid phase, (A.methylene chloride, B.liquid
water), the diffusion coefficient (DAB) is 1.07 10_5 cm2/s, and
the kinematic viscosity (nB) is 0.010 cm2/s.
a. At what position across the pond is the air flow no longer
laminar? Would it be reasonable to assume that the mean gas
film mass-transfer coefficient for methylene chloride in air is
dominated by turbulent flow mass transfer?
b. As part of an engineering analysis to predict the emissions
rate of methylene chloride (species A) from the pond,
determine the average gas film mass-transfer coefficient
 associated with the mass-transfer methylene chloride from
the liquid surface to the bulk air stream.
c. Compare the Schmidt number for methylene chloride in the
gas phase vs. the liquid phase, and explain why the values
are different.

a. omg I thought it was vA or vB but thats kinematic viscocity- thats why you can use the tables
b. also theyre doing 2*10^5 for Re but why not 2000

Mass Hw 5 Week 8
Andrew Sanchez

28.4

A thin (1.0-mm-thick) coat of fresh paint has just been

sprayed over a 1.5-m by 1.5-m square steel body part, which
approximates a flat surface. The paint contains a volatile solvent
that initially constitutes 30 wt% of the wet paint. The initial
density of the wet paint is 1.5 g/cm3. The freshly painted part is
introduced into a drying chamber. Air is blown into the rectangular
drying chamber at a volumetric flow rate of 60 m3/min, as
shown in the figure below, which has dimensions L.1.5 m,
H.1.0 m,W.1.5 m). The temperature of the air stream and the
steel body part are both maintained at 27C, and the total system
pressure is 1.0 atm. The molecular weight of the solvent is 78 g/
gmole, the vapor pressure exerted by the solvent at 27C is
105mm Hg, and the molecular diffusion coefficient of the
solvent vapor in air at 27C and 1.0 atm is 0.097 cm2/s.

What is the Schmidt number and the average Sherwood

number (ShL) for the mass-transfer process?
b. What is the estimated solvent evaporation rate from the
 surface of the whole body part in units of g/min? It may
be assumed that convection mass-transfer limits the evaporation
rate, and that the concentration of solvent vapor in the
bulk gas is finite, but can be approximated as cA; 0.
c. Using the results from part (a), how long will it take for the
paint to completely dry?
d. What are the hydrodynamic (d) and concentration boundarylayer
(dc) thicknesses at x.L.1.5 m? How does this
compare to H, the height of the drying chamber?
e. What is the new required air flow rate (m3/

b. Find Shmidt:
a. Sc= v(air)/DAB
i. Why is vair not Volumetric flow rate of air/ cross-sectional Area- vair
from appendix I- 1.5689*10^-5
ii. DAB- given above
iii. Vai
DAB= 0.097; %cm^2/s given
vair= 1.5689*10^-5; %m^2/s
Sc=vair/DAB*100^2= 1.6174
iv. Why is v0 volumetric flowrate
W=1.5;
H=1;
L=1.5;

%solve Shl
v0= 60; %m^3
vinf=v0/(W*H*60)
%need vinf for Re
%need Re for ShL
ReL= vinf*L/vair;
ShL=0.664*ReL^.5*Sc^(1/3)=196.778

%part b solve WA

kc=ShL*DAB/(L*100)% 0.127cm/s
CA=(105/760)/(82.06*300)%=5.61*10^-6;
W=kc*CA*L^2*100^2*78*60 = 75.1987 g/min

Part c

%masshw5 ch 28.4

clear all
clc

%28.7
%a solve Shmidt
DAB= 3*10^-6; %cm^2/s given
vair= 6*10^-3; %cm^2/s
Sc=vair/DAB
v0=30;
W=10;
L=20;

vinf=v0/(W*L)
Re=vinf*L/vair
ShL=0.664*Re^(1/2)*Sc^(1/3)
ShL=0.664*Sc^(1/3)*(v0/(W*vair))^0.5

kc=ShL*DAB/L
Ca=0.04
NA=kc*Ca

x=[0:0.1:L];
kc=ShL*DAB./x;
plot (x,kc)

28.12 A well-mixed open pond contains wastewater that is

contaminated with a dilute concentration of dissolved methylene
chloride. The pond is rectangular with dimensions of 500 m by
100 m, as shown in the figure (above right). Air at 27C and
1.0 atm blows parallel to the surface of the pond with a bulk
velocity of 7.5 m/s. At 20C and 1.0 atm, for the gas phase
(A.methylene chloride, B.air), the diffusion coefficient (DAB)
is 0.085 cm2/s, and kinematic viscosity (nB) is 0.15 cm2/s. At
27C, for the liquid phase, (A.methylene chloride, B.liquid
water), the diffusion coefficient (DAB) is 1.07 10_5 cm2/s, and
the kinematic viscosity (nB) is 0.010 cm2/s.
a. At what position across the pond is the air flow no longer
laminar? Would it be reasonable to assume that the mean gas

Chp 29 Convective
Mass Transfer b/w
phases- [Advanced
Thermo part. 1]
Hw 8

29.6 It is desired to recover hexane vapor (solute A) from air

using an absorption process. The absorption solvent is a nonvolatile
mineral oil, which has a mass density of 0.80 g/cm 3 and a
molecular weight of 180 g/gmole. In the dilute concentration
range, the equilibrium relationship for the dissolution of hexane
vapor in the mineral oil at 20C is described by pA;i
. H xA;i,
where H . 0.15 atm. At the present conditions of operation,
the hexane partial pressure in the bulk gas stream is 0.015 atm,
and the dissolved hexane in the bulk absorption solvent is
5.0 mole%. The total system pressure is 1.50 atm, and the
temperature is 20C. The liquid film mass-transfer coefficient
kx is 0.01 kgmole/m2 s, and the gas film mass-transfer
coefficient ky is 0.02 kgmole/m2 s.
a. What is the overall mass-transfer coefficient based on the
liquid phase, KL, and molar flux NA?
b. What is the composition of hexane at the gasliquid interface,
in terms of pA,i and xA,i?

%masshw8ch29p6.m
clear all
clc
format long

H= 0.15; % atm
P= 1.5; %atm
MWB= 180;%g/gmole
rhoB=0.8; %g/cm^3
kx= 0.01; % kgmole/m^2 s
ky= 0.02; % kgmole/m^2 s

kG=ky/P;
CL=rhoB/MWB;
kL= kx*MWB/rhoB*1000/100^3
KL=(100^3/1000*rhoB/(H*MWB*kG)+1/kL)^-1 %0.00105 m/s
%must convert H by multiplying by kG/rhoB due to units given in problem

xA= 0.05;
PA=0.015; %atm
Csat=PA/H*CL;
CAL=xA*CL;
NA=KL*(Csat-CAL)*100^3 % convert to mol / m s

xAi=.05-NA/kx/1000
PAi=H*xAi

Problem 29.19

29.19 Consider the waste treatment process shown in the figure

below. In this process, wastewater containing a dissolved TCE
concentration of 50 gmole/m3 enters a clarifier, which is essentially
a shallow, well-mixed tank with an exposed liquid surface.
The overall diameter of the tank is 20.0 m, and themaximumdepth
of the liquid in the tank is 4.0 m. The clarifier is enclosed to contain
the gases that are emitted from the wastewater. Fresh air is blown
into this enclosure to sweep away the gases emitted from the
clarifier, which are sent to an incinerator. Samples of effluent gas
phase and effluent liquid phase were measured for TCE content.
The TCE composition in the effluent gas is 4.0
mole%, whereas the TCE concentration dissolved in the effluent
liquid phase is 10 gmole TCE/m3 liquid. The clarifier operates at
1.0 atm and constant temperature of 20C. In independent pilot
plant studies, for the transferring species TCE, the liquid film
mass-transfer coefficient for the clarifier was kx . 200 gmole/
m2 s, whereas the gas film mass-transfer coefficient for the
clarifier was ky . 0.1 gmole/m2 s. Equilibrium data for the
air-TCE-water system at 20C is represented by Henrys law of
the form PA;i . H xA;i, with H . 550 atm.
Fresh air Effluent air +
4.0 mole% TCE
4.0 mole% TCE in air
cAL,o = 50 gmole/m3 cAL = 10 gmole TCE/m3
vo = ? m3/h
Effluent liquid
cAL = 10 gmole TCE/m3
D = 20 m
a. What is the overall mass-transfer coefficient based on the
liquid phase, as KL?
b. What is the flux of TCE from the clarifier liquid surface?
c. Develop a well-mixed, steady material balance model for the
process. What is the inlet volumetric flow rate of wastewater,
vo (in units of m3/hr) needed to ensure that the liquid effluent
TCE concentration is cAL . 10 gmole TCE/m3?

%masshw8ch29p19.m
clearall
clc
formatlong

H=550;%atm
P=1;%atm
MWB=18;
rhoB=1000;
kx=200;%gmole/m^2s
ky=0.1;%gmole/m^2s
KG=ky/P/1000;
kL=kx*MWB/rhoB/1000;
KL=(1/(H*KG)+1/kL)^1
kG=ky/P;
CL=rhoB/MWB;
kL=kx*MWB/rhoB*1000/100^3
KL=(100^3/1000*rhoB/(H*MWB*kG)+1/kL)^1%0.00105m/s
%mustconvertHbymultiplyingbykG/rhoBduetounitsgiveninproblem

xA=0.04;
PA=0.015;%atm
Csat=PA/H*CL;
CAL=10;
NA=KL*(CAL1000*Csat)*100^3%converttomol/ms
Final Exam Practice