PR#1 Root Finding

TM3101 Aplikasi Metode Numerik Perminyakan

f ( zi )

f ' ( zi )

z 4 + 2.5zi 1.3
zi +1 = zi i
4zi3 + 2.5
Mulai dengan z0 = 1.0 hingga a < 0.0001
Contoh perhitungan untuk z0 dan z1
f ( 1) = 14 + 2.5 ( 1) 1.3 = 2.2
3

f ' ( 1) = 4 ( 1) + 2.5 = 6.5

f ( z0 )

f ' ( z0 )

2.2
= 0.66154
6.5

n = z n z n 1
0 = 0

1 = 0.66154 1 = 0.33846

Karena 1 > 0.0001, iterasi kita lanjutkan

Berikut adalah hasil perhitungan selengkapnya

Zi

F(zi)

f'(zi)

zi+1

ei

0
1
2
3
4
5

1
0,661538462
0,512450909
0,495964945
0,495824675
0,495824665

2,2
0,545368916
0,050089151
0,000419127
2,90337E-08
0

6,5
3,658044606
3,0382906
2,987992263
2,987578333
2,987578305

0,661538462
0,512450909
0,495964945
0,495824675
0,495824665
0,495824665

0,33846154
0,14908755
0,01648596
0,00014027
9,7182E-09
0

Dari table di atas didapat harga z = 0.495824665 yang merupakan akar dari persamaan di
atas dengan = 0 < 0.0001
No.2 Solution of the Peng-Robinson EOS: Newton Method

The Peng-Robinson EOS

at
+ b ( Vm b ) ( Vm b ) = RT
p +

Vm ( Vm + b )

Dengan menggunakan variable dimensionless maka persamaan di atas dapat ditulis

atb
2b 3b 2p 2
b 2 b3p3
bp
2 at
2
Z +
1Z + p

Z + p 3 3 + 2 2 + 3 3 = 0
2 2 RT

RT
R T R T
R 2T 2
R T
R T
3

atau
Z3 + c 2 Z2 + c1Z + c0 = 0
b = 0.07780

Dengan

RTc
pc

a c = 0.45724

R 2Tc 2
pc

)(

a t = a c 1 + 0.37464 + 1.54226 0.269922 1 T Tc

Data-data yang diketahui:

T = 373.15 K

Pc = 3.72 MPa

P =1.5 MPa

= 0.193

R = 0.0083145 (untk P dlm KPa)

Perhitungan:
b=

0.07780 ( 0.0083145 ) ( 425.2 )

a c = 0.45724

= 7.33462176 x105

( 0.0083145 ) 2 ( 425.2 ) 2
3750

= 1.523954 x103

a t = 1.523954e3 1 + 0.37464 + 1.5422 ( 0.193) 0.26992 ( 0.193)

a t = 1.654196 x103
Sehingga
c2 =

bp
1
RT

c2 =

7.33462176e 5 x1500
1
0.0083145 x 373.15

) ( 1

373.15 425.2

C2 = -0.9645
a
2b 3b 2 P 2
c1 = P t

R 2T 2 RT R 2T 2

)(
2

3 7.34e 5 15002
2 7.34e5
1.654196e 3

c1 = 1500

2
2 0.0083145 x 373.15
2
2
0.0083145 x 373.15
0.0083145 x 373.15

C1 = 0.183079662
atb
b 2 b3 p3
c0 = p
+
+
R 3T3 R 2T 2 R 3T3

5
( 1500 )

7.34e
1.65e 3 7.34e5
7.34e 5
2
c0 = 1500
+
+
0.00831453 373.153 0.00831452 373.152 0.00831453 373.153

C0 = -0.00783882
Sehingga kita dapatkan persamaan sbb:
f(Z) = Z3 - 0.9645Z2 + 0.18308Z - 0.00784
f(Z) = 3Z2 1.929Z + 0.18308
Zi +1 = Zi

Z3 0.9645Z2 + 0.18308Z 0.00784

Menentukan selang:
Berikut nilai trial untuk z
Z
-2
-1
0
1
f(Z)
-12.232
-2.15542
-7.84e-3
0.21074
Karena f(a) < 0 dan f(b) > 0 maka didapatkan selang 0 < Z <1

2
4.50032

3
18.8609

Kita coba dengan Z0 = 1

zi

F(zi)

f'(zi)

zi+1

ei

0
1
2
3
4
5
6

1
0,831956494
0,751407938
0,729386881
0,727734827
0,727725766
0,727725766

0,21074
0,052734394
0,009413027
0,000614743
3,33521E-06
1,00063E-10
-1,38778E-17

1,25408
0,654690747
0,427455756
0,372108374
0,368073454
0,368051369
0,368051368

0,831956494
0,751407938
0,729386881
0,727734827
0,727725766
0,727725766
0,727725766

0,16804351
0,08054856
0,02202106
0,00165205
9,0613E-06
2,7187E-10
0

Dari table di atas didapatkan harga Z = 0.727725766 sebagai solusi dari fungsi factor
deviasi butane pada T = 373.15 K dan P = 1.5 MPa. Harga error yang didapat
= 0 < 0.001 sehingga hasil Z dapat diterima.

Note:
Jika didapatkan tiga nilai akar maka kita pilih nilai yang paling besar. Hal ini
dikarenakan nilai akar yang paling besar menunjukkan solusi jika substansi berada pada
fasa gas. Dan seperti yang kita tahu bahwa pada P = 1.5 MPa butane berada pada fasa gas.
No.3 Heat Balance: false Position Method
h = 4662 + 1.983T 0.001031T 2 + 4.243 106 T3 2.946 109 T 4 + 7.218 1013 T5
h = f(T)
T1 = 300 K f ( 300 ) = 4607.438 KJ Kg
hh
2= h+ 1
12000 KJ
h2 =
+ ( 4067.438 KJ Kg )
5 Kg
h 2 = 1667.438 KJ Kg
Sehingga didapat:
No.4 Well deliverability: Secant Method

Asumsi p = 5000 psi

4.38 p p wf

= 93.68 p wf 2 106
0.62

f ( p wf ) = 4.38 p p wf

0.5

93.68 p wf 2 106

f ( p wf ) = 4.38 50002 p wf 2
Secant Method:

0.62

x i +1 = x i

0.5

93.68 p wf 2 106

=0

0.5

=0

f ( x i ) ( x i 1 x i )
f ( x i 1 ) f ( x i )

Kita mulai dengan x0 = Pwf0 = 0

Pwfi

f(Pwfi)

0
1
2
3

1000
1100
1376,250498
1788,322381

164894,09
121068,47
72478,994
16470,727

Pwfi-1

f(Pwfi-1)

Pwfi+1

1000
164894,086 1376,2505
1100
121068,473 1788,32238
1376,2505 72478,9942 1909,50312

e
100
276,250498
412,071883

4
5
6
7
8

1909,503125
1917,188143
1917,241768
1917,241787
1917,241787

982,24563
6,8065316
0,0023758
5,995E-09
-3,201E-10

1788,3224
1909,5031
1917,1881
1917,2418
1917,2418

16470,7269
982,245634
6,80653163
0,00237575
5,9954E-09

1917,18814
1917,24177
1917,24179
1917,24179
1917,24179

121,180743
7,68501791
0,05362541
1,8724E-05
4,7294E-11

Dari table di atas dengan menggunakan Secant Method dapat diketahui bahwa salah satu
akar dari persamaan f(Pwf) di atas adalah pwf = 1917.241787 dengan a = 4.7294E-11
No.5 Isotherm Flash: Newton-Raphson Method
Persamaan untuk flash calculation adalah sebagai berikut

( Ki 1) zi
i 1) V + 1
i =1
( K1 1) z1 + ( K 2 1) z 2 + ( K3 1) z3
f ( Vi ) =
( K1 1) V + 1 ( K 2 1) V + 1 ( K 3 1) V + 1
( 61 1) 0.3396 ( 9 1) 0.0646 ( 0.035 1) 0.6068
f ( Vi ) =
+
+
( 61 1) V + 1 ( 9 1) V + 1
( 0.035 1) + 1
f ( Vi ) =

h =3

(K

20.376 0.5168
0.585562
+

60V + 1 8V + 1 0.965V + 1
1222.56
4.1344
0.5651
f ' ( Vi ) =

( 60V + 1) 2 ( 8V + 1) 2 ( 0.965V + 1) 2
f ( Vi ) =

Newton-Raphson Method
Vi +1 = Vi

f ( Vi )

f ' ( Vi )

Menentukan selang:
Menentukan selang:
Berikut adalah hasil perhitungan V dan f(V) beserta grafiknya
V

f(V)

-2
-1
0
1
2
3
4
5

-0.405530736
-0.717180432
20.307238
-16.33888785
0.828433253
0.44225028
0.304950282
0.233387335

25
20
15
10
5

Series1

0
-3

-2

-1
0
-5
-10
-15
-20

10

11

Dari tabel di atas kita ambil selang 0 < V < 1

Vi

F(Vi)

F'(Vi)

0
1

0
-1226.1293

20.3072
-0.0008

1246.44
-0.0008

1.34089E-07
1226.109619

20.3071

-1226.1293
1.34089E-07
1226.109619

-0.0008

20.3071

1.34093E-07
1226.109619

1.34093E-07
1226.109619

1.34093E-07

20.3071

3
4

-0.0008

-0.0008

1.34093E-07
1226.109619

1246.42

1246.42
-0.0008
1246.42

Dapat kita lihat disini bahwa akar dari persamaan di atas tidak dapat dicari dengan cara
iterasi pada Newton-Raphson Method (bentuk fungsinya tidak konvergence) sehingga
digunakan Exel dengan menggunakan Goal Seek. Hasil dari Goal Seek yang digunakan
adalah sebagai berikut:
V

f(V)

-0.8
-0.9
-0.95
0.8
0.6
0.4
0.3
0.35
0.38
0.39
0.395
0.398
0.3954
0.39545
0.39546
0.3955
0.3956
0.3957
0.3958
0.39578
0.395785
0.3957854
0.3957855
0.395786

-0.859688215
-0.781193808
-0.747657499
-2.082579813
-0.751077459
-0.01559642
0.400266232
0.177980988
0.059435003
0.021591577
0.002921909
-0.00820647
0.001435035
0.001249244
0.001212087
0.001063467
0.000691959
0.000320512
-5.08765E-05
2.33964E-05
4.82792E-06
3.34246E-06
2.97109E-06
1.11426E-06

Dengan menggunakan goal seek didapatkan akar yaitu V = 0.395786

Untuk menentukan komposisi fasa vapour (yi) dan komposisi dari fasa liquid (xi)
digunakan persamaan sbb:

yi =

K i zi
( K i 1)V + 1

xi =

zi
( K i 1)V + 1

Berikut adalah hasil perhitungan selengkapnya:

Ki

zi

yi

xi

61
9
0.035

0.3396
0.0646
0.6068

0.005431287
0.006136824
0.388571361

0.013723
0.015505
0.981771