Analisis Pengaruh Variasi Kondisi Operasional terhadap Kapasitas Aliran Kompresor Sentrifugal

Pendahuluan

Kompresor sentrifugal merupakan peralatan vital dalam berbagai industri proses, termasuk pada unit pemisah udara (ASU). Dalam pengoperasiannya, kapasitas aliran kompresor dapat dipengaruhi oleh beberapa faktor seperti temperatur inlet, rasio kompresi, dan daya motor. Artikel ini membahas bagaimana variasi kondisi operasional tersebut mempengaruhi kapasitas aliran kompresor sentrifugal berdasarkan pendekatan isothermal dan efisiensi kompresor yang diasumsikan tetap.

Contoh Spesifikasi Desain Kompresor

• Tipe: Sentrifugal multi-stage
• Medium: Nitrogen
• Flow Desain: 96000 Nm³/h
• Jangkauan Kendali: 75% ~ 105%
• Inlet Pressure: 0.522 MPa abs
• Outlet Pressure: 3.1 MPa abs
• Inlet Temperature: 37°C
• Shaft Power: 6850 kW
• Efisiensi Motor: 97.6%

Dasar Teori Koreksi Kapasitas Aliran

Untuk menyederhanakan analisis, proses kompresi dianggap isothermal, dan efisiensi kompresor dianggap tetap. Kerja kompresi isothermal reversibel per satuan massa gas ideal adalah:

\( w = R \cdot T_{in} \cdot \ln(PR) \)

Karena daya poros adalah P = ṁ × w / η, maka pada efisiensi kompresor (η) yang tetap, massa alir sebanding linear dengan daya poros, berbanding terbalik dengan temperatur inlet, dan berbanding terbalik dengan ln(PR). Satuan Nm³/h mewakili massa alir (volume pada kondisi normal), sehingga hubungan antara perubahan parameter dan kapasitas aliran (Q_real) dapat dihitung menggunakan persamaan koreksi relatif terhadap titik desain:

\( Q_{real} = Q_{design} \times \frac{P_{real}}{P_{design}} \times \frac{T_{design}}{T_{real}} \times \frac{\ln(PR_{design})}{\ln(PR_{real})} \)

Rasio kompresi (PR) dihitung berdasarkan tekanan absolut pada outlet dan inlet kompresor:

\( PR = \frac{Press_{outlet}}{Press_{inlet}} \)

di mana:

• Q: kapasitas aliran (misalnya dalam Nm³/h)
• T: temperatur inlet dalam Kelvin
• PR: rasio kompresi (outlet/inlet dalam tekanan absolut)
• P: shaft power (misalnya dalam kW)

Persamaan ini tidak menggunakan akar pangkat dua. Pada efisiensi konstan, jika daya poros naik dua kali lipat (dengan temperatur inlet dan rasio kompresi tetap), massa alir yang dikompresi juga naik dua kali lipat. Ini konsekuensi langsung dari kekekalan energi. Bentuk akar kuadrat memang muncul pada konteks lain yang strukturnya mirip, misalnya kompensasi tekanan/temperatur pada flowmeter differential-pressure (Q ∝ √(ΔP/ρ)), tetapi tidak berlaku untuk estimasi flow dari daya poros. Hukum afinitas (P ∝ N³) juga tidak berlaku di sini, karena kompresor beroperasi pada kecepatan tetap dengan pengaturan IGV, bukan dengan variasi kecepatan.

Faktor rasio kompresi menggunakan ln(PR), bukan PR langsung, karena kerja kompresi isothermal bersifat logaritmik terhadap rasio tekanan.

Catatan tentang Proses Isothermal dan Peran Intercooler/Aftercooler

Dalam perhitungan ini digunakan pendekatan proses isothermal, yaitu mengasumsikan bahwa suhu gas tetap konstan selama proses kompresi. Hal ini menyederhanakan perhitungan karena perubahan energi internal diabaikan, dan hanya kerja tekanan yang dihitung. Formula isothermal juga tidak memerlukan nilai k (Cp/Cv).

Namun secara nyata, kompresor sentrifugal bekerja mendekati proses adiabatik pada tiap stage: gas dikompresi tanpa pertukaran panas sempurna, sehingga suhu gas meningkat. Untuk menurunkan suhu gas kembali dan meningkatkan efisiensi keseluruhan, sistem biasanya dilengkapi dengan intercooler (di antara stage) dan aftercooler (di akhir proses kompresi). Pendinginan ini berfungsi:

• Menurunkan suhu gas sebelum masuk stage berikutnya → mengurangi beban kerja stage berikutnya,
• Menurunkan temperatur akhir sebelum masuk ke sistem proses downstream → menjaga integritas peralatan.

Tiap stage tetap merupakan proses adiabatik/politropik, bukan berubah menjadi isothermal. Yang terjadi adalah kerja kompresi total (seluruh stage plus pendinginan antar-stage) menjadi lebih dekat ke kerja isothermal. Semakin efektif pendinginan antar-stage, semakin dekat kerja aktual ke benchmark tersebut. Karena itu pendekatan isothermal layak dipakai sebagai model pembanding performa antar kondisi operasi, yang dikalibrasi relatif terhadap titik desain.

Pendekatan ini bukan pengganti flowmeter. Ketidakpastian absolutnya realistis di orde ±3–5% (akurasi pengukuran daya, variasi efisiensi motor, perubahan efisiensi kompresor akibat IGV dan fouling, rugi mekanis), dan presisinya menurun semakin jauh titik operasi dari titik desain.

Contoh perhitungan kapasitas jika daya motor 7800 kW, kondisi lain tetap seperti desain

Dengan motor 7800 kW dan efisiensi 97.6%, daya poros (shaft power) adalah:

\( P_{shaft} = 7800 \times 0.976 = 7612.8 \text{ kW} \)

Menggunakan formula:

\( Q_{real} = 96000 \times \frac{7612.8}{6850} = 106690 \text{ Nm}^3/\text{h} \; (\approx 111\% \text{ desain}) \)

Angka ini berada di luar jangkauan kendali (maksimum 105%). Artinya, pada temperatur inlet dan rasio kompresi desain, motor tidak akan pernah terbebani sampai 7800 kW. Batas operasi praktisnya ditentukan oleh bukaan IGV, batas choke, dan kurva kompresor, bukan oleh kapasitas motor.

Contoh perhitungan daya motor jika kapasitas 100800 Nm³/h, kondisi lain tetap seperti desain

Jika diketahui flow aktual sebesar 100800 Nm³/h (105% desain), dengan temperatur inlet dan rasio kompresi tetap sama seperti desain, maka daya motor yang dibutuhkan dapat dihitung dengan membalik formula:

\( P_{real} = P_{design} \times \frac{Q_{real}}{Q_{design}} \)

\( P_{real} = \frac{6850}{0.976} \times \frac{100800}{96000} = 7018.44 \times 1.05 = 7369.36 \text{ kW} \)

Hasil Perhitungan dan Analisis

1. Variasi Daya Motor (kondisi lain tetap seperti desain)

Analisis: Kapasitas aliran meningkat linear terhadap daya poros. Peningkatan flow 5% dari desain membutuhkan kenaikan daya poros sekitar 5% juga (daya motor ≈7369 kW), selama temperatur inlet, rasio kompresi, dan efisiensi kompresor tidak berubah.

2. Variasi Temperatur Inlet (daya dan rasio kompresi tetap)

Analisis: Penurunan temperatur inlet menurunkan kerja kompresi spesifik (kerja per kg gas sebanding dengan T_in). Pada daya poros yang sama, kompresor mampu mengompresi massa gas yang lebih besar, sehingga kapasitas aliran meningkat. Setiap penurunan 1°C memberi tambahan flow sekitar 0.3%. Efisiensi kompresor sendiri diasumsikan tetap dalam model ini. Yang turun adalah konsumsi energi spesifik (kWh per Nm³), bukan efisiensi kompresor yang naik.

3. Variasi Rasio Kompresi (daya dan temperatur inlet tetap)

Analisis: Rasio kompresi yang lebih tinggi membutuhkan kerja spesifik yang lebih besar (sebanding dengan ln(PR)), sehingga kapasitas aliran menurun secara gradual jika daya tetap. Pada PR sekitar 5.94, kenaikan PR sebesar 1% menurunkan flow sekitar 0.56% (= 1/ln(PR)). Sensitivitas ini lebih besar pada PR rendah: pada PR sekitar 2, kenaikan PR 1% menurunkan flow sekitar 1.4%. Tekanan inlet yang lebih rendah dengan target outlet tetap akan meningkatkan rasio kompresi dan konsumsi energi spesifik. Tetapi kondisi proses tetap harus diperhatikan, karena mungkin proses membutuhkan tekanan yang lebih tinggi dibandingkan flow yang lebih besar.

Formula Python

Berikut contoh penerapan kalkulasi dalam bahasa Python:

import math

# Titik desain
Q_DESIGN = 96000.0           # Nm3/h
P_DESIGN = 6850.0            # kW, shaft power
T_DESIGN = 37.0 + 273.15     # K
PR_DESIGN = 3.1 / 0.522      # rasio kompresi desain (~5.94)
ETA_MOTOR = 0.976

def flow_estimate(p_motor_kw, t_in_c, pr):
    """Estimasi flow (Nm3/h) - isothermal, efisiensi kompresor tetap."""
    p_shaft = p_motor_kw * ETA_MOTOR
    t_in = t_in_c + 273.15
    return (Q_DESIGN * (p_shaft / P_DESIGN)
            * (T_DESIGN / t_in)
            * (math.log(PR_DESIGN) / math.log(pr)))

def motor_power_required(q, t_in_c, pr):
    """Daya motor (kW) yang dibutuhkan untuk flow q (Nm3/h)."""
    t_in = t_in_c + 273.15
    p_shaft = (P_DESIGN * (q / Q_DESIGN)
               * (t_in / T_DESIGN)
               * (math.log(pr) / math.log(PR_DESIGN)))
    return p_shaft / ETA_MOTOR

print(f"Flow @ 7800 kW      : {flow_estimate(7800, 37, PR_DESIGN):.0f} Nm3/h")
print(f"Daya @ 100800 Nm3/h : {motor_power_required(100800, 37, PR_DESIGN):.2f} kW")

📎 Bisa langsung dicoba di: 👉 Python Online Compiler

Formula Structured Text

Contoh penerapan kalkulasi dalam bahasa Structured Text (ST) untuk sistem kontrol PLC atau DCS:

FUNCTION FC_FlowEstimateIsothermal : REAL
VAR_INPUT
    P_Motor_kW : REAL;   (* daya motor aktual, kW *)
    T_In_C     : REAL;   (* temperatur inlet aktual, degC *)
    P_In_MPa   : REAL;   (* tekanan inlet, MPa abs *)
    P_Out_MPa  : REAL;   (* tekanan outlet, MPa abs *)
END_VAR
VAR CONSTANT
    Q_DESIGN  : REAL := 96000.0;   (* Nm3/h *)
    P_DESIGN  : REAL := 6850.0;    (* kW shaft *)
    T_DESIGN  : REAL := 310.15;    (* K *)
    PR_DESIGN : REAL := 5.9387;    (* 3.1 / 0.522 *)
    ETA_MOTOR : REAL := 0.976;
END_VAR
VAR
    PR      : REAL;
    P_Shaft : REAL;
END_VAR

(* Proteksi: tekanan inlet harus valid dan PR > 1 agar LN(PR) > 0 *)
IF P_In_MPa <= 0.0 THEN
    FC_FlowEstimateIsothermal := 0.0;
    RETURN;
END_IF;
PR := P_Out_MPa / P_In_MPa;
IF PR <= 1.0 THEN
    FC_FlowEstimateIsothermal := 0.0;
    RETURN;
END_IF;

P_Shaft := P_Motor_kW * ETA_MOTOR;

(* Linear terhadap daya, 1/T, dan 1/LN(PR) - tanpa akar kuadrat *)
FC_FlowEstimateIsothermal := Q_DESIGN
    * (P_Shaft / P_DESIGN)
    * (T_DESIGN / (T_In_C + 273.15))
    * (LN(PR_DESIGN) / LN(PR));
END_FUNCTION

Formula WolframAlpha

Contoh penerapan kalkulasi untuk WolframAlpha (di WolframAlpha, log adalah logaritma natural):

📎 Link siap jalan: 👉 Klik di sini untuk mencari daya motor

Kesimpulan

Analisis variasi ini menunjukkan bahwa kapasitas aliran kompresor sentrifugal dipengaruhi oleh daya poros, temperatur masuk, dan rasio kompresi. Dengan pendekatan isothermal dan efisiensi tetap, flow berbanding lurus dengan daya poros, berbanding terbalik dengan temperatur inlet absolut, dan berbanding terbalik dengan ln(PR). Dalam sistem yang dikendalikan dengan motor tetap, pemahaman pengaruh ini penting untuk optimasi beban dan pemanfaatan maksimum kapasitas tanpa risiko overloading.

Pada kondisi desain, flow 105% hanya membutuhkan daya motor sekitar 7370 kW. Motor 7800 kW (desain ≈7018 kW) memberi margin sekitar 6% di atas kebutuhan batas atas jangkauan kendali. Margin ini berguna untuk menampung kondisi yang lebih berat dari desain, misalnya temperatur inlet lebih tinggi, rasio kompresi lebih besar, atau penurunan efisiensi akibat fouling. Karena itu pemilihannya tetap merupakan keputusan engineering yang logis untuk menjamin fleksibilitas operasi. Model ini bersifat relatif terhadap titik desain dan sebaiknya divalidasi terhadap data flowmeter atau compressor map sebelum dipakai sebagai dasar keputusan operasional.

Catatan revisi: versi sebelumnya artikel ini menggunakan akar kuadrat dan PR langsung dalam persamaan koreksi. Bentuk tersebut keliru: kerja kompresi isothermal sebanding dengan ln(PR), dan pada efisiensi tetap flow berbanding linear dengan daya poros. Akibatnya, versi lama meredam pengaruh variasi daya dan temperatur sekitar separuh dari yang seharusnya. Persamaan, contoh perhitungan, tabel, grafik, dan kode telah diperbarui.

Artikel ini ditulis dengan bantuan kecerdasan buatan dengan arahan, penyesuaian, dan validasi oleh penulis.