Membrane Separation of a Ternary Mixture​​Authors: Housam Binous and Ahmed Bellagi

reference : H . Binous, Gas Permeation Computations with Mathematica®, Chemical Engineering Education, 40 (2), pp . 140 - 144, 2006.
We compute the permeate and residue compositions obtained using a membrane separation module. Membrane area is also computed. The feed is a ternary mixture and the complete mixing model is used. This problem was solved using an iterative technique by C. J. Geankoplis in Transport Processes and Unit Operations, 3rd edition, Prentice Hall, 1993. In the present work, we use FindRoot to solve seven simultaneous equations with seven unknowns. Results obtained using both methods are in agreement.
Reference:
C. J. Geankoplis, Transport Processes and Unit Operations, 3rd edition, Prentice Hall, 1993.
Feed is a ternary mixture containing components A, B and C. Feed flow rate is
q
f
=1.0
4
10
3
cm
(STP)/s.
Feed composition is
x
fA
=0.25,
x
fB
=0.55,
x
fC
=0.20.
Membrane thickness is t=2.54
-3
10
cm.
The fraction of the feed that permeates through the membrane is known as the stage cut, θ=0.25, and is given by:
Permeabilities are expressed in
3
cm
(STP)cm(s
2
cm
cmHg).
Permeabiilties of components A, B and C are equal to
P'
A
=200
-10
10
,
P'
B
=50
-10
10
and
P'
C
=25
-10
10
.
The pressure of the feed side is
p
h
=300cmHg
while
p
l
=
30 cmHg is the pressure on the permeate side.
The three rate of permeation equations are called rate1, rate2 and rate3 where
A
m
is the membrane area.
In[]:=
rate1=
q
P
y
PA

P'
A
t
A
m
(
p
h
x
OA
-
p
l
y
PA
)
Out[]=
q
P
y
PA

A
m
(
p
h
x
OA
-
p
l
y
PA
)
′
P
A
t
In[]:=
rate2=
q
P
y
PB

P'
B
t
A
m
(
p
h
x
OB
-
p
l
y
PB
)
Out[]=
q
P
y
PB

A
m
(
p
h
x
OB
-
p
l
y
PB
)
′
P
B
t
In[]:=
rate3=
q
P
y
PC

P'
C
t
A
m
(
p
h
x
OC
-
p
l
y
PC
)
Out[]=
q
P
y
PC

A
m
(
p
h
x
OC
-
p
l
y
PC
)
′
P
C
t
The three material balance equations are called matbalance1, matbalance2 and matbalance3.
In[]:=
matbalance1=
x
OA

x
fA
-θ
y
PA
1-θ
Out[]=
x
OA

x
fA
-θ
y
PA
1-θ
In[]:=
matbalance2=
x
OB

x
fB
-θ
y
PB
1-θ
Out[]=
x
OB

x
fB
-θ
y
PB
1-θ
In[]:=
matbalance3=
x
OC

x
fC
-θ
y
PC
1-θ
Out[]=
x
OC

x
fC
-θ
y
PC
1-θ
The two summation rules are called summation1 and summation2.
In[]:=
summation1=
y
PA
+
y
PB
+
y
PC
1
Out[]=
y
PA
+
y
PB
+
y
PC
1
In[]:=
summation2=
x
OA
+
x
OB
+
x
OC
1
Out[]=
x
OA
+
x
OB
+
x
OC
1
In[]:=
x
fA
=0.25;
x
fB
=0.55;
x
fC
=0.20;
In[]:=
θ=0.25;
In[]:=
P'
A
=20010^-10;
P'
B
=5010^-10;
P'
C
=2510^-10;
In[]:=
t=2.54 10^-3;
In[]:=
q
P
=θ
q
f
;
In[]:=
q
f
=1.010^4;
In[]:=
p
h
=300;
p
l
=30;
Using FindRoot, we solve the seven simultaneous equations where the seven unknowns are
A
m
,
y
PA
,
y
PB
,
y
PC
,
x
OA
,
x
OB
and
​
​
x
OC
.
Results in agreement with those obtained using the iteration method presented by C. J. Geankoplis, Transport Processes and Unit Operations, 3rd edition, Prentice Hall, 1993. Permeate mole fraction of components A, B and C are 0.4552, 0.4502 and 0.0944. Residue mole fraction of components A, B and C are 0.1815, 0.5832 and 0.2351. The membrane area is 3.5417
6
10
2
cm
.
In[]:=
FindRoot[{rate1,rate2,rate3,matbalance1,matbalance2,matbalance3,summation1},{
y
PA
,0.2},{
y
PB
,0.2},{
y
PC
,0.2},{
A
m
,10^6},{
x
OA
,0.2},{
x
OB
,0.2},{
x
OC
,0.2}]
Out[]=

y
PA
0.455281,
y
PB
0.450286,
y
PC
0.0944335,
A
m
3.54176×
6
10
,
x
OA
0.181573,
x
OB
0.583238,
x
OC
0.235189