Page 1
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 1/13
T-xydiagramforMETHA-01/
CHLOR-01
Liquid/vapormolefraction,METHA-0
1
Temperature, C
0.000
0.02
0.00
0.0!
0.100
0.12
0.10
0.1!
0.200
0.22
0
.20
0.2!
0."00
0."2
0."0
0."!
0.#00
0.#2
0.#0
0.#!
0.00
0.2
0.0
0.!
0.$00
0.$2
0.$0
0.$!
0.!00
0.!2
0.!0
0.!!
0.%00
0.%2
0.%0
0.%!
0.&00
0.&2
0.&0
0.&!
1.000
1##
1#$
1#%
10
12
1#
1$
1%
1$0
1$2
1$#
1$$
'12.0atm
(12.0atm
)ri*c+er, Ariel
A110$2!!%
CE 122, rofe**or pat3ie4ic5
ro6lem 7et 1
18 a8 i8
Page 2
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 2/13
a8 ii8
y-x diagram for METHA-01/CHLOR-01
Liquid/vapor mole fraction, METHA-01
9 a p o r m o l e f r a c t i o n ,
M E T
H A - 0
1
0.00 0.0 0.10 0.1 0.20 0.2 0."0 0." 0.#0 0.# 0.0 0. 0.$0 0.$ 0.!0 0.! 0.%0 0.% 0.&0 0.& 1.000.00
0.0
0.10
0.1
0.20
0.2
0."0
0."
0.#0
0.#
0.0
0.
0.$0
0.$
0.!0
0.!
0.%0
0.%
0.&0
0.&
1.00
12.0 atm
Page 3
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 3/13
a8 iii8
Activity coecientsfor METHA-01/CHLOR-01
Liquid/vapor molefraction, METHA-01
A c t i v i t ( c o e : c i e n t *
0. 000 0. 02 0. 00 0. 0! 0. 100 0.12 0. 10 0. 1! 0. 200 0. 22 0. 20 0. 2! 0. "00 0. "2 0. "0 0. "! 0. #00 0. #2 0. #0 0. #! 0. 00 0.2 0. 0 0. ! 0. $00 0. $2 0. $0 0. $! 0. !00 0. !2 0. !0 0. !! 0. %00 0. %2 0. %0 0. %! 0. &00 0.&2 0. &0 0. &! 1. 0000.&0
0.&2
0.&#
0.&$
0.&%
1.00
1.02
1.0#
1.0$
1.0%
1.10
Liquid1METHA-01 12.0atm
Liquid1CHL;-01 12.0atm
68 i8
Page 4
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 4/13
T-xydiagramforMETHA-0
1/CHLOR-01
Liquid/vapormolefraction,MET H
A-01
Temperature, C
0.000
0.02
0.00
0.0!
0.100
0.12
0.10
0.1!
0.200
0.22
0.20
0.2!
0."00
0."2
0."0
0."!
0.#00
0.#2
0.#0
0.#!
0.00
0.2
0.0
0.!
0.$00
0.$2
0.$0
0.$!
0.!00
0.!2
0.!0
0.!!
0.%00
0.%2
0.%0
0.%!
0.&00
0.&2
0.&0
0.&!
1.000
1#0
1#2
1##
1#$
1#%
10
12
1#
1$
1%
1$0
1$2
1$#
1$$
'12.0atm
(12.0atm
68 ii8
Page 5
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 5/13
y-x diagram for METHA-01/CHLOR-01
Liquid/vapor mole fraction, METHA-01
9 a p o r m o l e f r a
c t i o n ,
M E T H A - 0
1
0.00 0.0 0.10 0.1 0.20 0.2 0."0 0." 0.#0 0.# 0.0 0. 0.$0 0.$ 0.!0 0.! 0.%0 0.% 0.&0 0.& 1.000.00
0.0
0.10
0.1
0.20
0.2
0."0
0."
0.#0
0.#
0.0
0.
0.$0
0.$
0.!0
0.!
0.%0
0.%
0.&0
0.&
1.00
12.0 atm
68 iii8
Page 6
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 6/13
ActivitycoecientsforME
THA-01/CHLOR-01
Liquid/vapormolefraction,MET
HA-01
Activit(coe:cient*0.000
0.02
0.00
0.0!
0.100
0.12
0.10
0.1!
0.200
0.22
0.20
0.2!
0."00
0."2
0."0
0."!
0.#00
0.#2
0.#0
0.#!
0.00
0.2
0.0
0.!
0.$00
0.$2
0.$0
0.$!
0.!00
0.!2
0.!0
0.!!
0.%00
0.%2
0.%0
0.%!
0.&00
0.&2
0.&0
0.&!
1.000
1.01.1
1.21."
1.#1.
1.$1.!
1.%1.&
2.02.1
2.22."
2.#2.
2.$2.!
2.%2.&
".0".1
".2"."
".#".
".$".!
".%".&
#.0#.1
Liquid1METHA-0112.0atm
Liquid1CHL;-0112.0atm
c8 i8
Page 7
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 7/13
T-xydiagramforMETHA-0
1/CHLOR-01
Liquid/vapormolefraction,MET
HA-01
Temperature, < 0.000
0.02
0.00
0.0!
0.100
0.12
0.10
0.1!
0.200
0.22
0.20
0.2!
0."00
0."2
0."0
0."!
0.#00
0.#2
0.#0
0.#!
0.00
0.2
0.0
0.!
0.$00
0.$2
0.$0
0.$!
0.!00
0.!2
0.!0
0.!!
0.%00
0.%2
0.%0
0.%
!
0.&00
0.&2
0.&0
0.&!
1.000
#12
#1#
#1$
#1%
#20
#22
#2#
#2$
#2%
#"0
#"2
#"#
#"$
#"%
'12.0atm
(12.0atm
c8 ii8
Page 8
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 8/13
y-x diagram for METHA-01/CHLOR-01
Liquid/vapor mole fraction, METHA-01
9 a p o r m o l e f r a
c t i o n ,
M E T H A - 0
1
0.00 0.0 0.10 0.1 0.20 0.2 0."0 0." 0.#0 0.# 0.0 0. 0.$0 0.$ 0.!0 0.! 0.%0 0.% 0.&0 0.& 1.000.00
0.0
0.10
0.1
0.20
0.2
0."0
0."
0.#0
0.#
0.0
0.
0.$0
0.$
0.!0
0.!
0.%0
0.%
0.&0
0.&
1.00
12.0 atm
c8 iii8
Page 9
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 9/13
ActivitycoecientsforME
THA-01/CHLOR-01
Liquid/vapormolefraction,MET
HA-01
Activit(coe:cient*0.000
0.02
0.00
0.0!
0.100
0.12
0.10
0.1!
0.200
0.22
0.20
0.2!
0."00
0."2
0."0
0."!
0.#00
0.#2
0.#0
0.#!
0.00
0.2
0.0
0.!
0.$00
0.$2
0.$0
0.$!
0.!00
0.!2
0.!0
0.!!
0.%00
0.%2
0.%0
0.%!
0.&00
0.&2
0.&0
0.&!
1.000
1.01.1
1.21."
1.#1.
1.$1.!
1.%1.&
2.02.1
2.22."
2.#2.
2.$2.!
2.%2.&
".0".1
".2"."
".#".
".$".!
".%".&
#.0#.1
Liquid1METHA-0112.0atm
Liquid1CHL;-0112.0atm
Page 10
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 10/13
"8 a8
68
Page 11
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 11/13
c8
d8 i8 T+e mi'ture *tream temperature remain* t+e *ame a* t+e inlet* 4+en t+e deal
model i* applied in =a8 and =68 6ut 4+en t+e ;TL model i* applied t+e mi'ture
temperature ri*e* 6( 2.&>C. T+e ;TL model do doe* not +ave t+e *ame outlet
*tream temperature a* t+e t4o inlet temperature* 6ecau*e it ta3e* into account
molecular interaction* 4+ile t+e ideal model doe*n?t.
ii8 T+e volume rate of t+e mi'ed *tream i* equal to t+e *um of t+e t4o inlet *tream*
for t+e ideal model* 6ut not equal for t+e ;TL model. t i* not equal in t+e ;TL
model 6ecau*e t+i* model ta3e* into account additional molecular interaction* 4+ile
t+e ideal model doe*n?t.
#. @tili5ed Matla6 to complete t+e*e computation*
Matla6 Code
clc;
clear all;
syms x; % x is the quality in this case
e = 1e-5; % Setting the tolerance
dx = e + 1;
x = [.5, .5, .5, .5, .5, .5, .5] ; % nitially assumed quality
! = [.""#$, .1&1, ."'($, .&)$", ."5'$, .1&1, .151]; % *ale ata
= [1).&, 5.&, &.), 1.$', .$1, .#&, .&']; % *ale ata
= sum/ /!.0/1- .2 /1+x.0/-1 ; % 3achord43ice lash unction
1 = sum/ /!.0/1-.6& .2 /1+x.0/-1.6& ; % eri7ati7e o the unction
count = "; % 8ill count ho9 many iterations are ta:en
9hile /as/dx e % 8ill continue iterating until the error is less than
tolerance
dx = .21; % <art o equation (-1"
Page 12
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 12/13
x = x - dx; % >dating the 7alue o x
= sum/ /!.0/1- .2 /1+x.0/-1 ; % >date the 3achord43ice lash
unction
1 = sum/ /!.0/1-.6& .2 /1+x.0/-1.6& ; % >date the deri7ati7e o
unction
count = count + 1; % ncrement the counter
i/as/dx ? e % 8ill out>ut liquid and 7a>or com>ositions i the
quality change is less than the tolerance
liquid@om>osition = !.2/1-x+x.0;
7a>or@om>osition = .0liquid@om>osition;
i/xor/sum/liquid@om>osition A= 1,sum/7a>or@om>osition A= 1
%@hec: to see i liquid and 7a>or com>ositions add u> to equal 1
>rint/BCrrorD Eiquid or Fa>or @om>osition Sum is not equal to 1D
GnB;
else
[ro9s, columns] = si!e/x;
% they do, >rint 7a>or2liquid com>ositions
>rint/BEiquid Cxit Stream @om>ositionH GnB;
or I = 1Hcolumns,
>rint/Bx%d = %). GnB, I, liquid@om>osition/I;
end
>rint/BGnFa>or Cxit Stream @om>ositionH GnB;
or I = 1Hcolumns,
>rint/By%d = %). GnB, I, 7a>or@om>osition/I;
end
end
end
i /count &""
>rint/BCrrorD Solution is not con7ergingD GnB; % Sho9 error
message i the unction isnBt con7erging
rea:;
endend
i /count ? &""
>rint/BGn*he Juality isH %). GnB, x/1; %Sho9 the ans9er i the
unction con7erged
>rint/BKumer o iterations = %d GnB,count;
end
utput in Command indo4
Liquid E'it 7tream Compo*ition
'1 B 0.001
'2 B 0.0#0
'" B 0.0#
'# B 0.1!
' B 0.0$2
Page 13
7/17/2019 PS1
http://slidepdf.com/reader/full/ps1563db83b550346aa9a91c5cf 13/13
'$ B 0.1$
'! B 0.21
9apor E'it 7tream Compo*ition
(1 B 0.01#
(2 B 0.20%
(" B 0.11%
(# B 0."#$
( B 0.0$
($ B 0.112
(! B 0.1#$
T+e ualit( i* 0.#&
um6er of iteration* B 2