Pengantar Hysys
Pengantar Hysys
Pengantar Hysys
Published
data
meassurements
Thermodynamic
Property
Activity coefficient
(Wilson, NRTL,
UNIQUAC, UNIFAC)
EoS: virial, vdW,
RK, SRK, PR, Gen.
Corrr.)
Process
simulation
Primary choice
Secondary choice
Air Separation
Peng-Robinson (PR)
Soave Redlich-Kwong
(SRK)
Corresponding
States
Gas Processing
PR, SRK
BWRS
Gas Treating
Kent-Eisenberg
Electrolyte NRTL
Petrolium Refining
BK10, Chao-Seader,
Grayson Streed, PR,
SRK, Lee-Kessler-Plocker
Heavy crude
characterization
Petrochemicals-VLE
Petrochemicals-LLE
NRTL, UNIQUAC
Chemicals
UNIFAC
Data, parameters
Electrolytes
Pitzer
Data, parameters,
databank, model
polyelectrolytes
Polymer
Sanchez-Lacombe
EoS,
UNIFAC-FV
Data, parameters,
databanks,models for
polar polymer and block
copolymer
NRTL, UNIQUAC,
UNIFAC
Problem area
Data, parameters
Data, parameters, models
for VLLE
vdW Eos
Empirical Eos
Wong-Sandler
PSRK
Activity
Coef.
2000
1990
MathIas
Huron-Vidal
Peng-Robinson
Soave
Theoritical
EoS
UNIFAC
UNIQUAC
1975
MBWR
NRTL
ASOG
Wilson
1950
Redlich-Kwong
BWR
Beatie-Bridgemen
1900
Van Laar
Meyer
Margules
Teori
PC-Saft
Simulasi
Eksperimen
Saft
PVT behaviour
Diagram PV
Pc
SCF
Diagram PT
Padat
Cair
Triple
point
Uap/Gas
Tc
Cair+Uap
VL
VV
Tc
T<Tc
T>Tc
f ( P, V , T ) 0
Gas ideal (simplest EoS)
-volume individual = 0
PV RT
- no interaction
valid pada tekanan rendah
PV ZRT
Untuk gas ideal, Z = 1
Generalized correlation
2-parameter corresponding state principle (CSP)
Z Z Tr , Pr ,
Pitzer and Curl correlation (1955, 1957)
Z Z 0 Z 1
Dimana Z0 dan Z1 fungsi Tr dan Pr
Tc
2P
0;
2
0
Tc
RT
a
P
2
V b V
Cair+Uap
VL
volume
Tc
T<Tc
VV
27 R 2Tc2
RT
a
;b c
64 Pc
8Pc
T>Tc
Intermolecular attraction
RT
a
P
2
V b V
1949, Redlich-Kwong
(RK)
RT
a
P
12
V b T V b V
RT
a T
P
V b V V b
1976, Peng-Robinson
(PR)
RT
a T
P
2
V b V 2bV b 2
Persamaan umum
Entropy:
Enthalpy:
V
dH C p dT V T
dP
T P
dS Cp
dT V
dP
T T P
EoS
EoS
Gas ideal
EoS gas ideal:
PV RT
dH C p dT
dT
dP
dS Cp
R
T
P
R
P
Equation of State
Ideal gas
Van der Waals
Redlick-Kwong (RK)
Peng-Robinson (PR): khusus dikembangkan
untuk industri petrolium
PR-Stryjek and Vera (PRSV): extended dari PR
untuk menangani non-ideal systems
Kabadi Danner: modifikasi dari SRK untuk VLE
H2O-HC in the dilute region
Sour PR
Sour SRK
Equation of State
Activity Models
Untuk sistim polar atau non-ideal
Jauh lebih empiris dari EoS
Tidak akurat jika diekstrapolasi ke kondisi
operasi diluar range yang ditetapkan (sesuai
data experimen dalam evaluasi parameter)
Parameter-parameter di fitting dari data
eksperimen
Aplikasi sangat terbatas sampai tekanan
sedang.
Untuk keadaan tertentu bisa dilakukan perhitungan
Asimetris (Dual Approach):
Activity model for the liquid
EoS for the gas
Gas Non-polar
Elektrolit
NRTL, Pitzer,
sepanjang BIP
tersedia
Try PRSV,
Wilson, NRTL,
UNIQUAC
Cek BIP
BIP bisa
diestimasi
dari UNIFAC
RESUME APLIKASI PP
Tipe Sistim
Rekomendasi FP
TEG Dehydration
PR
Sour Water
PR, Sour PR
Cryogenic Gas
Processing
PR, PRSV
Air Separation
PR, PRSV
PR
Vacuum Tower
RESUME APLIKASI PP
Tipe Sistim
Rekomendasi FP
Ethylene Tower
High H2 Systems
PR, ZJ GS*
Chemical Systems
T V T L ; PV P L
f V f L
i
i
Isofugacitykomposisi
Koef. Aktifitas:
xiiV yiiL
yi P i xi Pi sat
Gas nyata
RK, SRK dst.
Campuran
gas nyata (i dan j)
Gas ideal
PV=RT
(No intermolecular interaction
P=0 atau V=tak berhingga)
Campuran
gas ideal (i dan j)
(No intermolecular interaction
berarti tidak ada interaksi antar
komponen)
Qm yi y j Qij
i 1 j 1
Qij
Qii Q jj
2
Qij Qii Q jj
1 k
1 l
ij
12
ij
i 1
f id x f
i
i i
Jika aij = 0
Larutan nyata
(Campuran)
Aktifitas koefisien dapat dihitung dengan:
-polynomial expansion: Van Laar, Margules
-semi teoritis model local composition:
Wilson, NRTL, UNIQUAC dst.
Membutuhkan: parameter
interaksi aij yang difitting dari
data eksperimen
Ki
y
xi
iL
Ki V
iV T ,V V , y ;iL T ,V L , x
P P T ,V V , y P T ,V L , x
yiiV P xi i f i 0
i T , x
Kerugian
Umumnya
tidak cukup
akurat untuk fasa uap
dan cair dari komponen
polar, kecuali ada
modifikasi empiris
Mixing rule yang
sederhana tidak mampu
merepresentasikan
kelakuan fasa dari
komponen polar
Keuntungan
Kerugian
Harus
memperhatikan
kondisi reference
Memberikan problem
yang signifikan pada
komponen superkritis
dan pada daerah
mendekati titik kritis
Enthalpy, entropy dan
density harus dihitung
dengan persamaan yang
terpisah
Liquid
Phase
model
Vapor
Phase
model
Vapor
Phase
model
Virial EoS
Fugasitas fasa uap yang menunjukkan interaksi kuat
Contoh: Carboxylic acids
Zat yang punya kecendrungan ikatan hidrogen
kuat pada fasa uap
Virial koefisien bisa dihitung (Hypothetical) dengan
HYSYS jika data komponen murni berikut tersedia:
Tekanan dan suhu kritis
Dipole momen
Rata-rata gyration radius
Association parameter
Association parameter untuk masing-masing
pasangan biner
Range Tekanan:
yP
T
2
i 1
m
i ci
yT
i 1
i ci
Lee-Kesler (1975)
Termasuk generalized correlation, pengembangan dari
metode Pitzer
Z Z 0 Z R Z 0
Dimana:
Z0=compressibility factor untuk simple fluid
ZR=compressibility factor untuk reference fluid
Direkomendasi untuk Tr<0.8
Lee-Kesler Plocker
Termasuk generalized correlation, pengembangan dari
metode Pitzer
Z Z R Z R Z 0
RT
a
P
12
V b T V b V
9 a
3Tr
Pr
0.5
Vr 3 b Tr Vr Vr 3 b
dimana
a 0.42748
b 0.08664
Tc2.5
Tc
a a R
; b b R
Pc
Pc
2
Mixing rule:
b xi bi
i
aij ij ai a j
Binary interaction parameter, ij
Mixing rule 1:
Mixing rule 2:
didifinisikan
RT
a T ,
P
V b V V b
Cukup akurat terutama untuk Hydrocarbon
a T , ac
1 S 1 Tr0.5
RT
a T ,
P
2
V b V 2bV b 2
dimana
a ac
R 2Tc2
ac 0.45724
Pc
RTc
b 0.077480
Pc
1 1 Tr0.5
0.37464 1.5422 0.26992 2
Akurasi dari PR dan SRK
hampir sama namun PR lebih
akurat dalam perhitungan
densitas cairan
HYSYS PR EoS
Sama dengan PR EoS dengan pengembangan pada
range aplikasi untuk perhitungan VLE sistim yang
lebih komplek
PRSV EoS
RT
a T ,
P
2
V b V 2bV b 2
Modifikasi pada
i 1 i 1 T
0.5
r
a ac
R 2Tc2
ac 0.45724
Pc
RTc
b 0.077480
Pc
Gas
Cair
EOS
Gas
Cair
yi i f i 0
Ki
xi
Pi
Activity Coeff.
models
HypNRTL
Sama dengan persamaan original NRTL yang berbeda
hanya parameter sebagai fungsi temperatur
lnGij ij ij
ij
aij bijT
RT
General NRTL
Menggunakan binary interaction parameter lebih
banyak
Digunakan untuk menghitung VLE dan LLE secara
simultan dan range boiling point atau konsentrasi
besar
CH2
ethanol
H2O
ethanol
H2O
a
ij exp ij
RT
CH3
OH
OH
CH2
CH3
H2O
H2O
amk
RT
imk exp
MODEL SELECTION
EoS
Diskripsi
Aplikasi
Ideal Gas
PengRobinson PR
Digunakan
HysysPR
EOS
PengRobinson
Stryjek-Vera
(PRSV)
Braun K 10
SoaveRedlichKwong
(SRK)
Diskripsi
Modifikasi
Aplikasi
Fasa Uap dan
cair
EOS
Diskripsi
Aplikasi
Kabadi
Danner
(KD)
LeeKeslerPlocker
RedlichKwong
Hanya fasa
uap
EOS
Diskripsi
Aplikasi
ZudkevitchJoffee
Virial
Activity Model
Diskripsi
Ideal Solution
Regular Solution
Model
NRTL
Activity Model
General NRTL
Description
Merupakan variasi parameter dari model NRTL
menggunakan 5 parameter dan lebih fleksibel
dari model NRTL.
B C
ij Aij ij ij2 FijT Gij ln T
T T
Aplikasi:
Untuk komponen yang mempunyai boiling point
dengan range lebar.
Jika dibutuhkan penyelesaian yang serempak
dari VLE dan LLE dan ada range yang lebar
terhadap boiling point atau konsentrasi.
UNIQUAC
Menggunakan
Activity Model
Wilson
Chien-Null
Margules
Description
Model
Activity Model
Description
Van Laar
Khususnya
UNIFAC
VLE/LLE
Reference
HYPOTHETICALS
HYSYS
T fp ,Tb ,
Corresponding
state principle
Pure Properties
0
0
0
Ideal gas: C p , H f , G f
Real gas :( PVT
Group Contribution
Method
(Based on UNIFAC)
), H ,G , S
HYPOTHETICALS
Dengan HYPO Manager kita dapat membuat
komponen yang tidak terdapat di library atau
komponen-komponen Hypothetical.
Komponen Hypothetical dapat berupa komponen
murni, defined mixtures, undefined mixtures, atau
solids.
Metode estimasi untuk physical properties murni
dan campuran tersedia (CSP, Group contribution)
Metode untuk estimasi binary interaction parameter
antara komponen hyphotetical dan library tersedia.
HYPO MANAGER
Dengan memilih tab Hypotheticals dari Simulation Basis Manager,
tampilan berikut muncul:
Diskripsi
View
Add
Delete
Translocate
Import
Export
Diskripsi
View Hypo
View Group
Move Hypos
Clone Comps
Buka sebuah new case pada HYSYS dan pilih tab Hypotheticals dari
simulation Basis Manager.
Dari Hypothetical group, click tombol Add untuk membuat
Hypothetical Group yang baru.
Enter HypoAlcohol sebagai nama Group baru.
Builder.
Rumus kimia dari ethanol adalah C2H5OH, dan terdiri dari groups
CH3, CH2, dan OH. Pilih CH3 pada daftar Available UNIFAC
Groups.
Click tombol Add Group (s), sesuai dengan struktur kimia dan
sesuaikan jumlah fungsional grup yang menyusunnya (angka 1
defaault jumlahnya).
UNIFAC Base dan Critical Properties. Click Close icon untuk menutup
tampilan dan kembali ke tampilan input Tabular Hypothetical Input.
Diskripsi
Group Name
Component
Class
Estimation
Methods
Estimate
unknown
Props
Option
Diskripssi
Clone Library
Comps
Notes
DISKRIPSI
View
Add Hypo
Button
Description
Add Solid
Delete
UNIFAC
Base Properties
Base Property
Description
NBP
MW
Molecular Weight
Liq. Density
Liquid Density
Tc
Critical Temperature
Pc
Critical Pressure
Vc
Critical Volume
Acentricity
Acentric Factor
Boiling Point
Unknown
Estimation Methods
Group
Diskripsi
Property to
set Methods
For
View
Group
Diskripsi
Estimation
Method For
Selected Property
View
Property
Critical
Temperature
Critical
Pressure
Default Method
Available
Methods
Aspen, Bergman,
Cavett, Chen Hu,
Eaton Poerter,
Edmister, Group
Contribution, Lee
Kesler, Mathur,
Meissner
Redding, Nokay,
Riazi Dauber,
Roess,
PennState,
Standing, Twu
Aspen, Bergman,
Cavett, Edmister,
Group
Contribution, Lee
Kesler, Lydersen,
Mathur,
PennState, Riazi
Dauber, Rowe,
Standing, Twu
Variables
Affected
Critical
Temperature
Standard Liquid
Density
CONSTALD
Variables
Viscosity
Thetas
Critical
Temperature
Standard Liquid
Density
CONSTALD
Variables
Viscosity
Thetas
Property
Default Method
Available Methods
Critical
Volume
Pitzer
Group Contribution.
Pitzer,Twu
Acentricity
Molecular
Weight
if NBP<155 oF
Bergman is used
all other cases, Lee
Kesler is used
API,Aspen,Aspen Leastq,
Bergman,Hariu
Sage,,Katz
Firoozabadi,Katz nokay,
Lee Kesler, PennState,
Riazi Dauber,Robinson
Peng,Twu,Whitson
Variables Affected
Critical Temperature
Standard Liquid Density
CONSTALD Variables
Viscosity Thetas
w
Gs
Standard Liquid Density
CONSTALD Variables
Viscosity Thetas
Molecular Weight
Property
Default Method
Available
Methods
Variables
Affected
Normal Boiling
Point
Hyprotech
proprietary method
Twu
Vapour Pressure
Liquid Density
for Hydrocarbon,
Lee Kesler is
used
all other cases,
Riedel is used
Yen Woods
Bergman,BergmanP
NA,Chueh
Prausnitz,Gunn
Yamada,Hariu
Sage,,Katz
Firoozabadi,Lee
Kesler,
Twu,Whitson,Yarbor
ough,Yen woods
Normal Boiling
Point
Viscosity Thetas
Antoine Coeficient
PRSV_kappa
Standard Liquid
Density
CONSTALD
Variables
Property
Default Method
Available
Methods
Variables
Affected
Ideal Gas
Enthalpy
Cavett
Cavett,Falon
Watson,Group
Contribution,Lee
Kesler,Modified Lee
Kesler
Ideal H Coefficient
Heat of
Formation
for chemical
structure defined
in UNIFAC groups,
Joback is used
all other
cases,this formula
is used:
Group Contribution
Heat of Formation
Heat of
Combustion
Hyprotech
proprietary method
Group Contribution
Gibbs Coefficient
Property
Default Method
Available
Methods
Variables
Affected
Hyprotech
proprietary, Letsou
Stiel
Viscosity Thetas
Surface
Tension
Brock Bird
Tabular Variables
Radius of
Gyration
Hyprotech proprietary
method
Default only
Liquid
Viscosity
Critical
Temperature
Critical Pressure
Normal Boiling
Point
Molecular Weight
Standard Liquid
Density
Property
Liquid Enthalpy
Vapour Enthalpy
Chao Seader
Acentricity
Chao Seader
Solubility Parameter
Cavett Parameter
Dipole Moment
Property
Enthalpy of Combustion
Vapour Viscosity
Chung
PRSV Kappa 1
Kfactor1
UNIFAC Structure
Objects
Description
UNIFAC
Menampilkan type dan jumlah dari Sub Groups
Structure Group pada Struktur UNIFAc
Add Group(s)
Delete Group
Objects
Description
Free Bonds
Status Bar
Available UNIFAC
Groups
UNIFAC Structure field Menampilkan struktur kimia dari molekul yang sedang
dibangun.
UNIFAC Calculated
Base Properties
UNIFAC Calculated
Critical Properties
Sub Group
Description
Highlighting the
Sub Group
Typing in the
UNIFAC
Structure Input
field
ID Tab
Tab pertama pada hypo
Property View
Critical Tab
Menampilkan based propertis dan critical propertis
Point Tab
Menampilkan Additional Point Properties untuk hypothetical
TDep Tab
Menampilkan Temperature Depedent Properties untuk hypothetical
Ada tiga pilihan informasi pada tampilan :
1. Vapour Enthalpy
Perhitungan didasarkan pada basis massa. Poin reference
untuk persamaan adalah gas ideal pada 0 K.
2. Vapour Pressure
Perhitungan menggunakan persamaan Modified Antoine
OIL MANAGER
HYSYS
Mean Avg. BP
sp gr 60F/60F
P
E
T
R
O
L
I
U
M
Input
Data Lab.
HYSYS
Hypocomponent
1.Characterize assay
2.Generate hypocomponent
3.Install oil in flowsheet