1. PENDAHULUAN

Sampai pertengahan tahun 1976 refrigeran halokarbon menjadi pilihan utama sebagai fluida kerja untuk mesin pendingin [1]. Dari segi teknik memang refrigeran jenis ini mempunyai keunggulan dibanding refrigeran lainnya, ditambah juga sifatnya yang tidak beracun dan tidak mudah terbakar. Penggunaan refrigeran halokarbon disepakati untuk dikurangi setelah terbukti bahwa

refrigeran jenis ini termasuk dalam kategori jenis zat yang berpotensi merusak lapisan ozon [2].

Salah satu refrigeran yang sekarang ini dikembangkan adalah refrigeran hidrokarbon. Refrigeran hidrokarbon adalah refrigeran campuran berbagai zat yang termasuk dalam golongan hidrokarbon. Zat-zat yang terlibat dalam pencampuran ini tidak dapat bereaksi satu sama lain karena kuatnya ikatan atom karbon pada tiap zat

tersebut. Sifat-sifat campuran hidrokarbon dapat didekati dari sifat-sifat komponen murni yang membentuknya dengan menyertakan parameter-parameter tambahan seperti fugasitas dan faktor asentrik.

2. TERMODINAMIKA CAMPURAN HIDROKARBON

Sesuai dengan tingkat keadaannya, campuran hidrokarbon dapat berada dalam satu fasa dan dapat pula berada dalam dua fasa. Campuran yang terdiri atas satu fasa disebut campuran homogen, sedangkan campuran yang terdiri atas dua fasa atau lebih disebut campuran heterogen. Persamaan-persamaan berikut ini diambil dari buku Applied Hydrocarbon Thermodynamics karangan Wayne C. Edmister dan Byung Ik Lee.

Persamaan-persamaan pada campuran homogen dapat diterapkan pada campuran heterogen jika semua fasa berada dalam tingkat kesetimbangan. Semua persamaan termodinamika untuk zat murni dapat diterapkan langsung pada campuran homogen. Sifat-sifat termodinamika sebagai fungsi dari temperatur dan volume adalah sebagai berikut:

1. Energi dalam

\[\frac{U - U^*}{RT} = \frac{1}{RT} \int \left\{ T \left( \frac{\partial P}{\partial T} \right)_V - P \right\} dV \tag{1}\]

2. Entalpi

\[\frac{H - H^*}{RT} = Z - 1 + \frac{1}{RT} \int_{-\infty}^{V} \left\{ T \left( \frac{\partial P}{\partial T} \right)_{V} - P \right\} dV \tag{2}\]

3. Entropi

\[\frac{S - S_0^*}{R} + \ln \frac{P}{P_0} = \ln Z + \frac{1}{R} \int_{\infty}^{V} \left\{ \left( \frac{\partial P}{\partial T} \right)_{V} - \frac{R}{V} \right\} dV \quad (3)\]

4. Fugasitas

\[ln\left(\frac{f}{P}\right) = \frac{1}{RT} \int_{0}^{P} \left(V - \frac{RT}{P}\right) dP \tag{4}\]

Sembarang sifat termodinamika pada campuran homogen dapat didefinisikan oleh dua sifat termodinamika lain. Sedangkan pada campuran heterogen kuantitas tiap komponen juga mempengaruhi sifat-sifat tersebut. Persamaan-persamaan termodinamika pada campuran heterogen adalah sebagai berikut:

1. \[dU = (TdS) - (PdV) + \sum_{I}^{N} \left( \frac{\partial U}{\partial n_i} \right)_{S,V,n_i} dn_i\] (5)

2. \[dH = (TdS) - (VdP) + \sum_{i=1}^{N} \left(\frac{\partial H}{\partial n_i}\right)_{S,P,n_i} dn_i\] (6)

3. \[dA = -(SdT) - (PdV) + \sum_{i=1}^{N} \left(\frac{\partial A}{\partial n_i}\right)_{T,V,n_i} dn_i \qquad (7)\]

4. \[dG = -(SdT) + (VdP) + \sum_{i=1}^{N} \left(\frac{\partial G}{\partial n_i}\right)_{T,P,n_i} dn_i\] (8)

Suku-suku terakhir pada keempat persamaan merupakan perubahan sifat terhadap perubahan komposisi. Turunanturunan ini bernilai sama dan disebut potensial kimia.

3. KESETIMBANGAN FASA

Setiap sistem terbuka pada campuran heterogen dapat ditinjau sebagai sistem homogen dari suatu sistem tertutup jika keseluruhan sistem berada pada tingkat keadaan kesetimbangan. Pada tingkat keadaan ini tidak ada perbedaan potensial yang dapat menyebabkan aliran massa atau energi melewati batas fasa.

Karena tidak ada perubahan sifat-sifat sistem pada tingkat keadaan kesetimbangan, maka perubahan entropi total, perubahan volume total, perubahan energi dalam total, dan jumlah mol total sistem juga bernilai nol. Berangkat dari kenyataan ini, tingkat keadaan kesetimbangan sistem yang terdiri atas N komponen tercapai jika temperatur dan tekanan tiap fasa adalah sama serta potensial kimia komponen di tiap fasa adalah sama. Dalam bentuk persamaan dapat ditulis sebagai berikut:

\[T^{(1)} = T^{(2)} = \dots = T^{(M)}\] \[P^{(1)} = P^{(2)} = \dots = P^{(M)}\] \[\mu_i^1 = \mu_i^2 = \dots = \mu_i^M \quad (i = 1, 2, \dots, N)\] (9)

Parameter lain yang digunakan dalam kesetimbangan fasa adalah konstanta kesetimbangan. Konstanta kesetimbangan komponen i dalam campuran \((K_i)\) didefinisikan sebagai berikut

\[K_i = \frac{Y_i}{X_i} \tag{10}\]

Konstanta kesetimbangan dapat pula dinyatakan dalam sifat termodinamika, yaitu:

\[K_i = \frac{\phi_i^L}{\phi_i^V} \tag{11}\]

4. PERSAMAAN TINGKAT KEADAAN

Persamaan sifat termodinamika yang telah diturunkan di atas dapat dipecahkan jika diketahui hubungan PVT dan sifat-sifat gas ideal. Untuk sistem campuran diperlukan satu data tambahan, yaitu komposisi campuran. Umumnya hubugan PVT diambil dari persamaan tingkat keadaan yang telah ada. Dari penelitian yang telah dilakukan sebelumnya, umumnya persamaan tingkat keadaan Soave memberikan hasil yang lebih akurat dalam memprediksi sifat-sifat campuran hidrokarbon [2]. Persamaan Soave diturunkan dari persamaan Redlich-Kwong dan berbentuk

\[P = \frac{RT}{V - b} - \frac{a}{V(V + b)} \tag{12}\]

Koefisien a merupakan perkalian dari parameter fungsi temperatur dan parameter bukan fungsi temperatur. Koefisien a dapat ditulis dalam bentuk

\[a = a' f(T) \tag{13}\]

Persamaan (12) dapat pula ditulis dalam bentuk persamaan pangkat tiga dari volume V sebagai berikut:

\[V^{3} - \left(\frac{RT}{P}\right)V^{2} + \left(\frac{a}{P} - \frac{bRT}{P} - b^{2}\right)V - \frac{ab}{P} = 0\] (14)

Dalam diagram PV bentuk persamaan di atas diperlihatkan pada Gambar (1). Seperti yang ditunjukkan pada gambar tersebut, semua garis isotermal di bawah

garis isotermal kritis mempunyai tiga akar, yaitu A, B, dan D. Titik A dan D menyatakan volume cair jenuh dan volume uap jenuh, sedangkan titik B tidak mempunyai arti fisik. Ketika temperatur mendekati temperatur kritis, tiga akar ini menjadi lebih dekat satu sama lain dan akhirnya menjadi identik di titik kritis. Kondisi ini secara aljabar dapat dijabarkan sebagai berikut

\[(V - V_c)^3 = V^3 - 3V_cV^2 + 3V_c^2V - V_c^3 = 0\] (15)

Gambar 1 Diagram PV Zat Murni

Pada titik kritis, persamaan (14) dan (15) harus identik. Dengan demikian nilai koefisien kedua persamaan harus sama pula, yaitu

\[3V_c = \frac{RT_c}{P_c} \tag{16}\]

\[3V_c^2 = \frac{a_c}{P_c} - \frac{bRT_c}{P_c} - b^2 \tag{17}\]

\[3V_c^3 = \frac{a_c b}{P_c} \tag{18}\]

\[a_c = a' f(T_c) \tag{19}\]

Jika persamaan (16), (17), dan (18) digabung maka akan diperoleh

\[b^{3} + (3V_{c})b^{2} + (3V_{c}^{2})b - V_{c}^{3} = 0\] (20)

Persamaan (20) ini dapat juga ditulis dalam bentuk

\[(b^3 + 3b^2V_c + 3bV_c^2 + V_c^3) = 2V_c^3\] (21)

atau

\[(b + V_c)^3 = 2V_c^3\] \[b = (2^{V/3} - 1)V_c\] (22)

Jika persamaan (16) dan persamaan (22) digabung maka akan diperoleh

\[b = \frac{(2^{1/3} - 1)RT_c}{3P_c} = 0,08664 \frac{RT_c}{P_c}\] (23)

Menggabungkan persamaan (16), (18), dan (23) akan memberikan

\[a_c = \frac{(RT_c)^2}{9(2^{1/3} - 1)P_c} = 0.42748 \frac{(RT_c)^2}{P_c}\] (24)

Dari persamaan (13) dan (19) akan diperoleh

\[a = a_c \frac{f(T)}{f(T_c)} \tag{25}\]

Nilai \(f(T)/f(T_c)\) pada persamaan di atas biasanya diganti dengan \(\alpha\), sehingga persamaan tersebut dapat ditulis dalam bentuk

\[a = a_c \alpha \tag{26}\]

Pada saat sistem berada pada titik kritis, \(\alpha\) akan bernilai satu.

Jika parameter V pada persamaan (12) diganti dengan ZRT/P maka akan diperoleh persamaan baru dalam bentuk derajat tiga, yaitu

\[Z^{3} - Z^{2} + (A - B - B^{2})Z - AB = 0\] (27)

dimana

\[A = \frac{aP}{\left(RT\right)^2} \tag{28}\]

\[B = \frac{bP}{RT} \tag{29}\]

Faktor penting yang membedakan persamaan Soave dengan persamaan-persamaan sebelumnya adalah dalam mendefinisikan \(\alpha\). Persamaan (25) dan (26) menyatakan bahwa \(\alpha\) tergantung pada T dan Tc sistem tersebut. Soave menambahkan satu parameter tambahan dalam menentukan \(\alpha\) dan disebut sebagai faktor asentrik. Faktor asentrik digunakan sebagai ukuran bentuk molekul sistem dan nilainya telah ditabelkan.

Soave mendefinisikan nilai \(\alpha\) sebagai berikut

\[\alpha = \left[1 + m\left(1 - T_r^{0.5}\right)\right]^2 \tag{30}\]

\[m = 0.48 + 1.574\omega - 0.176\omega^2 \tag{31}\]

5. ATURAN PENCAMPURAN

Persamaan tingkat keadaan yang telah dijabarkan di atas berlaku untuk zat murni. Persamaan tersebut dapat juga digunakan untuk zat campuran dengan menerapkan aturan pencampuran. Aturan pencampuran ini dimaksudkan untuk merefleksikan pengaruh komposisi terhadap sifat termodinamika zat campuran.

Untuk campuran yang terdiri atas N komponen, koefisien persamaan Soave adalah

\[a = \sum_{i}^{N} \sum_{j}^{N} x_{i} x_{j} \left( a_{i} a_{j} \right)^{0.5} \left( 1 - k_{ij} \right)\] (32)

\[b = \sum_{i=1}^{N} x_i b_i \tag{33}\]

Koefisien interaksi dapat diperoleh dari data kesetimbangan uap-cair secara eksperimen. Untuk campuran hidrokarbon, koefisien interaksi ini bernilai nol. Jika koefisien interaksi bernilai nol maka persamaan (32) dapat disederhanakan menjadi

\[a = \left[\sum_{i}^{N} x_{i} a_{i}^{0.5}\right]^{2} \tag{34}\]

6. FUNGSI TERMODINAMIKA CAMPURAN DARI PERSAMAAN TINGKAT KEADAAN SOAVE

engan menggunakan hubungan PVT yang telah peroleh dari persamaan tingkat keadaan Soave, peroleh persamaan-persamaan sifat termodinamika bagai berikut:

Entalpi

\[\frac{H - H^*}{RT} = Z - 1 - \frac{A}{B} \left[ 1 + \frac{1}{a} \sum_{i}^{N} \sum_{j}^{N} x_i x_j m_j \right]\] \[\left( a_i a_{cj} T_{rj} \right)^{0.5} \left( 1 - k_{ij} \right) \ln \left( 1 + \frac{B}{Z} \right)\] (35)

\[\frac{S - S_0^*}{R} + \ln \frac{P}{P_0} = \ln (Z - B) - \frac{A}{B}\] \[\left[ \sum_{i}^{N} \sum_{j}^{N} x_i x_j m_j (a_i a_{cj} T_{rj}) (1 - k_{ij}) \right] + \ln \left( 1 + \frac{B}{Z} \right)\] (36)

Koefisien fugasitas

\[\ln \phi_i = -\ln(Z - B) + (Z - 1)B_i^*\] \[-\frac{A}{B} \left(A_i^* - B_i^*\right) \ln\left(1 + \frac{B}{Z}\right)\] (37)

dimana:

\[A_{i}^{\bullet} = \frac{1}{an} \left[ \frac{\partial (n^{2} a)}{\partial n_{i}} \right]_{T,n_{i}}\] \[B_{i}^{\bullet} = \frac{1}{b} \left[ \frac{\partial (nb)}{\partial n_{i}} \right]_{T,n_{i}}\]

7. PERANGKAT LUNAK PERHITUNGAN SIFAT-SIFAT TERMODINAMIKA

Carena terlalu banyaknya parameter yang terlibat dan proses iterasi dalam orde ribuan, diperlukan bantuan computer untuk melakukan perhitungan sifat-sifat ermodinamika campuran hidrokarbon. Perangkat lunak yang dibuat dapat digunakan untuk membantu nenganalisis fenomena lain sepanjang berkaitan dengan sifat-sifat termodinamika campuran hidrokarbon.

Perhitungan sifat termodinamika selalu dimulai dengan perhitungan tingkat keadaan jenuh campuran berdasarkan tekanan atau temperatur yang diketahui. Diagram alir perhitungan temperatur uap jenuh ditampilkan pada Gambar 2. Pada kasus ini parameter yang diketahui adalah tekanan dan fraksi mol campuran. Pola perhitungan yang sama dapat diterapkan untuk perhitungan temperatur cair jenuh. Pada kasus ini, tebakan nilai \(K_i\) digunakan untuk menghitung \(X_i\) (fraksi mol bagian cair) dan harus selalu diperbaharui sampai jumlah totalnya satu.

Diagram alir untuk perhitungan tekanan jenuh sama dengan diagram alir untuk perhitungan temperatur jenuh. Untuk temperatur dan fraksi mol campuran sebagai parameter yang diketahui, akan diperoleh nilai tekanan

jenuh campuran. Perhitungan dapat dilanjutkan dengan perhitungan entalpi jenuh dan entropi jenuh dengan persamaan (35) dan persamaan (36).

Gambar 2 Diagram Alir perhitungan Temperatur Uap Jenuh

Jika tekanan, temperatur, dan fraksi mol suatu campuran hidrokarbon diketahui, dengan melakukan perhitungan tingkat keadaan jenuh maka fasa campuran akan diketahui. Untuk kasus temperatur sistem lebih kecil dibandingkan temperatur cair jenuh, fasa campuran adalah cair. Sebaliknya untuk kasus temperatur sistem lebih besar dibandingkan temperatur uap jenuh, fasa campuran adalah uap. Sedangkan untuk kasus temperatur sistem diantara temperatur cair jenuh dan temperatur uap jenuh, fasa campuran adalah cair dan uap. Saat memecahkan persamaan (27) selalu akan diperoleh tiga nilai Z. Nilai Z terkecil adalah nilai Z fasa cair, nilai Z terbesar adalah nilai Z fasa uap. Sedangkan nilai Z yang lain tidak mempunyai arti apa-apa. Metode yang digunakan dalam pemecahan akar persamaan (27) adalah Metode Newton-Rapshon. Jika fasa sistem adalah cair, maka untuk mendapatkan nilai Z terkecil adalah dengan memberikan tebakan Z sebesar nol. Sebaliknya jika fasa sistem adalah uap maka tebakan Z adalah satu. Perhitungan tingkat keadaan kesetimbangan akan dilakukan jika fasa sistem adalah campuran.

Metode perhitungan yang sama dapat dilakukan untuk kasus tiga parameter lainnya yang diketahui, misalnya tekanan, entalpi, dan fraksi mol campuran. Dalam kasus ini entalpi sistem yang akan dijadikan landasan untuk memeriksa fasa sistem. Contoh diagram alir perhitungan untuk kasus tekanan, temperatur, dan fraksi mol sebagai parameter yang diketahui ditampilkan pada Gambar 3.

Gambar 3 Diagram Alir Perhitungan Sifat termodinamika

8. HASIL PERHITUNGAN DAN PEMBAHASAN

Dengan bantuan perangkat lunak yang telah dibuat aka dilakukan analisis pengaruh etana terhadap refrigera campuran hidrokarbon. Sebelum perangkat luna digunakan, terlebih dahulu dilakukan pemeriksaan has keluaran program dengan data eksperimen dan data databel sifat termodinamika. Dari hasil pemeriksaa tekanan jenuh dan temperatur jenuh kesalahan yan terjadi adalah sekitar 2%. Sedangkan untuk entalp spesifik, volume spesifik, dan entropi spesifik, kesalaha yang terjadi sekitar 0,1% [3]

Analisis ini akan dilakukan terhadap dua jenis refrigera campuran hidrokarbon, yaitu HCR-12 dan HCR-22 HCR-12 adalah refrigeran untuk mengganti R-12 sedangkan HCR-22 adalah refrigeran untuk menggan R-22. Semua parameter yang ditinjau akan selah dibandingkan dengan R-12 dan R-22.

Tinjauan pertama adalah tekanan dan temperatur jenul campuran. Gambar (4) dan Gambar (5) menampilkan grafik hubungan tekanan dan temperatur jenuh ini untul berbagai persentase massa kandungan etana pada HCR 12 dan HCR-22.

Gambar 4 Kurva Tekanan-Temperatur Jenuh R-12 dan HCR-12 dengan Berbagai Persentase massa Kandungan Etana

Gambar 5 Kurva Tekanan-Temperatur Jenuh R-22 dan HCR-22 dengan Berbagai Persentase Massa Kandungan Etana

engan memperhatikan Gambar (4) diketahui bahwa makin besar persentase massa etana yang terkandung da HCR-12, semakin besar pula penympangan kurva T jenuh terhadap R-12. Keadaan ini terjadi pada kanan rendah maupun tekanan tinggi. Kurva HCR-12 engan kandungan etana 0% merupakan kurva yang aling mendekati kurva R-12.

ntuk HCR-22, sampai tekanan 4 bar terlihat enyimpangan yang semakin besar jika persentase massa tana semakin besar. Namun semakin besar tekanan urva R-22 semakin mendekati kurva HCR-22 dengan andungan etana 4%. Pada tekanan tinggi kurva HCR-22 engan kandungan etana 4% merupakan kurva yang aling mendekati kurva R-22. Sebaliknya pada tekanan endah, kurva HCR-22 dengan kandungan etana 0% nerupakan kurva yang paling mendekati kurva R-22.

Finjauan kedua adalah nilai COP masing-masing efrigeran. Sistem pendingin yang akan ditinjau litampilkan pada Gambar (6) dan dengan ketentuan sebagai berikut:

• Refrigeran masuk kompresor pada tingkat keadaan uap jenuh.
• Refrigeran keluar kondensor pada tingkat keadaan cair jenuh.
• Proses isentropis pada kompresor dan isoentalpi pada katup ekspansi.
• d. Untuk HCR-12 tekanan sistem adalah 2 bar dan 10 bar. Untuk HCR-22 tekanan sistem adalah 5 bar dan 16 bar.

Gambar 6 Sistem Pendingin

Hasil perhitungan COP untuk berbagai persentase kandungan etana dalam HCR-12 ditampilkan pada Gambar (7).

Terlihat jelas bahwa COP sistem pendingin menurun dengan bertambahnya kandungan etana dalam HCR-12. Sistem pendingin yang menggunakan HCR-12 dengan 0% etana mempunyai COP sebesar 3,73. Sementara sistem pendingin yang menggunakan HCR-12 dengan 10% etana mempunyai COP sebesar 2,65. Dalam kasus ini terjadi penurunan COP sebesar 28,9% jika kandungan etana dalam HCR-12 sebesar 10%. Sebagai bahan

perbandingan, untuk sistem pendingin yang sama dengan menggunakan R-12 mempunyai COP sebesar 3,78.

Gambar 7 Kurva COP sistem pendingin untuk berbagai nilai persentase kandungan etana dalam HCR-12

Kurva COP sistem pendingin terhadap persentase kandungan etana dalam HCR-22 ditampilkan pada Gambar (8). Terlihat bahwa COP sistem pendingin juga menurun dengan bertambahnya kandungan etana dalam HCR-22. Terjadi penurunan COP sebesar 50,1% jika kandungan etana dalam HCR-22 sebesar 10%. Untuk sistem pendingin yang sama dengan menggunakan R-22 mempunyai COP sebesar 5,27.

Gambar 8 Kurva COP sistem pendingin untuk berbagai nilai persentase kandungan etana dalam HCR-

Temperatur refrigeran pada evaporator dan kondensor juga menunjukkan kecenderungan menurun dengan bertambahnya kandungan etana pada HCR-12 dan HCR-22. Kurva temperatur evaporator dan temperatur kondensor terhadap persentase kandungan etana dalam HCR-12 dan HCR-22 ditampilkan pada Gambar (9) dan Gambar (10).

Gambar 9 Kurva temperatur evaporator dan temperatur kondensor sistem pendingin untuk berbagai nilai persentase kandungan etana dalam HCR-12

Gambar 10 Kurva Temperatur Evaporator dan Temperatur Kondensor sistem pendingin\nuntuk berbagai nilai persentase kandungan\netana dalam HCR-12

9. KESIMPULAN

A. Pengaruh Etana pada HCR-12

Semakin besar kandungan etana semakin besar pula penyimpangan kurva P-T jenuh refrigeran terhadap kurva P-T jenuh pada R-12. Kondisi ini juga diikuti dengan menurunnya COP sistem pendingin. Kandungan 10% etana dapat menurunkan COP sistem sebesar 28,9%. Sementara temperatur evaporator dan temperatur kondensor juga menurun dengan bertambahnya kandungan etana dalam HCR-12.

B. Pengaruh Etana pada HCR-22

Pada tekanan dibawah 4 bar, besar penyimpangan kurva P-T jenuh HCR-22 terhadap kurva P-T jenuh R-22 sebanding dengan kandungan etana dalam HCR-22 tersebut. Pada tekanan di atas 4 bar, semakin besar tekanan kurva P-T jenuh R-22 semakin mendekati kurva P-T jenuh HCR-22 dengan kandungan etana sebesar 4%. Sementara COP sistem pendingin menurun dengan bertambahnya kandungan etana dalam HCR-22. Kandungan 10% etana dapat menurunkan COP sistem pendingin sebesar 50,1%. Temperatur evaporator dan temperatur kondensor juga menurun dengan bertambahnya kandungan etana dalam HCR-22.

SUMBER DANA

Penelitian ini didana dari Proyek Penelitian Hibah Bersaing IV Tahun Anggaran 1998/1999.

DAFTAR DAN ARTI LAMBANG

• a Konstanta persamaan Soave
• A Konstanta persamaan Soave
• b Konstanta persamaan Soave
• B Konstanta persamaan Soave
• a_c Konstanta persamaan Soave pada tingkat keadaan kritik

• f Fugasitas
• H Entalpi (kJ)
• H Entalpi pada tingkat keadaan standar (kJ)
• kii Koefisien interaksi
• P Tekanan (bar)
• P_c Tekanan pada tingkat keadaan kritis (bar)
• P* Tekanan pada tingkat keadaan standar (bar)
• Q Kalor (Watt)
• S Entropi (kJ/K)
• S Entropi (kJ/K)
• T Temperatur (K)
• T_c Temperatur pada tingkat keadaan kritis (K)
• T* Tekanan pada tingkat keadaan standar (K)
• U Energi dalam (kJ)
• U* Energi dalam pada tingkat keadaan standar (kJ)
• V Volume (m³)
• W. Daya kompresor (Watt)
• X Fraksi mol fasa cair
• Y Fraksi mol fasa uap
• Z Faktor kompresibilitas
• ∝ Konstanta fungsi temperatur
• μ Potensial kimia (kJ)
• φ Koefisien fugasitas