Abstrak —Listrik  meringankan  pekerjaan  manusia.  Dalam mengalirkan listrik pembangkit listrik tenaga uap menggunakan  motor  listrik  sebagai starting awal.  Untuk menjaga  motor  listrik  tetap  bekerja  maksimal  diperlukan proses pendinginan. Selama ini motor listrik hanya didinginkan menggunakan  fan,  namun  masalah  lain  muncul  ketika  udara luar lembab atau saat udara mengandung debu/partikulat. Hal ini  menyebabkan  dinding  motor  listrik  ditempeli  oleh debu/partikulat  yang  akan  menyebabkan  panas  tidak  bisa keluar  dari  motor,  sehingga  terjadi overheating dan  motor mudah mengalami kerusakan. Tujuan dari penelitian ini adalah merancang sebuah alat pendingin motor dengan menggunakan alat  penukar  panas compact  tipe  plate  fin .  Pada  penelitian  ini data operasi diperoleh dari  pembangkit listrik tenaga uap dan perancangan desain alat penukar panas berdasarkan spesifikasi yang ada pada buku Kays and London.  Setelah dihitung akan dibandingkan  nilai  UA LMTD dengan  UA desain .  Variasi  pada penelitian  ini  adalah  beban  yang  harus  didinginkan  oleh heat exchanger fungsi  temperatur  masuk  dan  keluar  dari heat exchanger .  Hasil dari penelitian diperoleh dimensi plate fin heat exchanger yang  memiliki  spesifikasi  sesuai  desain  permukaan 5.3 dari buku Kays and London dengan panjang 1,556 m, lebar 0,897  m,  dan  tinggi  1,299  m.  Sedangkan  nilai  UA LMTD sebesar 17,364 kW/K dan UA desain sebesar 17,599 kW/K. Kesimpulannya pada  analisa  performa  terhadap  variasi  beban  yang  harus didinginkan  yaitu  semakin  tinggi  beban  yang  didinginkan semakin  tinggi  pula  nilai effectiveness .  Hal  ini  terlihat  pada beban 105 % dengan efektivitas 0,8845.
Kata  Kunci —desain,  plate  fin,  ∆T LMTD ,  performa  heat exchanger

## I. PENDAHULUAN
EBUTUHAN  listrik  terus  meningkat  dari  tahun  ke tahun.  Oleh  karena  itu  dibutuhkan  dibutuhkan  suatu pembangkit  tenaga yang  memadai.  Berbagai  jenis pembangkit  tenaga  dikembangkan  di  Indonesia.  Salah satunya pembangkit listrik tenaga uap.
Gambar  1.  merupakan  siklus  sederhana  dari  pembangkit listrik tenaga uap yang terdiri pesawat penting seperti: boiler, turbin, condenser, boiler feed pump , dan lain-lain.
Siklus  sistem  pembangkit  listrik  tenaga  uap  (PLTU) bermula dari air yang dipompa lalu  masuk ke dalam boiler. Di  dalam  boiler  terjadi  pembakaran  lalu  air  menjadi  uap. Kemudian uap panas digunakan untuk menggerakkan turbin serta generator yang menghasilkan listrik.  Jenis motor yang digunakan  adalah  BFP  ( Boiler  Feed  Pump )  untuk mengalirkan  fluida,  dimana  pompa  ini  digerakkan  oleh sebuah motor listrik. Dalam pengoperasiannya motor bekerja terus-menerus.
Motor  yang  dioperasikan  tiap  hari  membutuhkan pemeliharaan  ( maintenance ).  Caranya  dengan  melakukan pendinginan.  Selama  ini  pendinginan  adalah  dengan mengalirkan udara bebas  yang ada di sekeliling ruangan ke dalam dinding motor listrik melalui winding yang dipasang di bagian kanan dan kiri dari motor listrik.
Gambar 1. Instalasi Pembangkit Listrik Tenaga Uap (Moran, 2005)

## Gambar 2. Skema pola piker
Namun  ada  keterbatasan  pada  proses  pendinginan  ini  yaitu pada  saat  udara  bebas  yang  dialirkan  mengandung debu/kotoran.  Debu/kotoran  akan  menempel  pada  dinding motor  dan  akan  mengakibatkan overheating karena  panas dari  motor  listrik  tidak  bisa  terbuang  dengan  baik. Keterbatasan  lain  saat  udara  lembab,  udara  bebas  akan mengandung uap air. Dan akan mengakibatkan motor listrik cepat rusak.
Pada penelitian ini akan dirancang sebuah heat exchanger dengan menggunakan tipe plate fin heat exchanger , dimana fluida yang mengalir pada sisi plate adalah udara dan fluida yang mengalir pada sisi fin adalah udara ( air to air ). Besarnya losses motor  dapat dihitung  dengan  perumusan power  input  dikurangi  power  output,  dimana losses motor merupakan  beban  yang  harus  didinginkan  oleh heat exchanger yang  akan  dirancang.  Untuk  satu  motor  akan dirancang dua buah heat exchanger .
Dari perhitungan losses didapatkan persamaan-persamaan sebagai berikut:
Power Input = 𝑉 × 𝐼 × 𝑐𝑜𝑠 Ø √3 (1) 𝑄 = 𝑃 − 𝑊 out motor (2) dimana:
V = Voltase (kV) I = Arus listrik (Ampere) Q = beban yang harus didinginkan oleh heat exchanger
(kW)
P = daya motor (kW)
Persamaan untuk mendapatkan dimensi heat exchanger yang sesuai dengan spesifikasi yang ada adalah sebagai berikut: 𝑄 𝑚𝑜𝑡𝑜𝑟 = 𝑚̇ 𝑢𝑑𝑎𝑟𝑎 . 𝑐 𝑝 . ∆𝑇 (3)
dimana:
Q motor = laju perpindahan panas (kW) 𝑚̇ 𝑢𝑑𝑎𝑟𝑎 = laju alir massa fluida (kg/s)
Persamaan untuk  menghitung temperatur udara yang keluar dari heat exchanger :

## Desain Compact Heat Exchanger Tipe Plate Fin sebagai Pendingin Motor pada Boiler Feed Pump
Nur Sapti Marsheliyana dan Djatmiko Ichsani Jurusan Teknik Mesin, Fakultas Teknologi Industri, Institut Teknologi Sepuluh Nopember (ITS)
Jl. Arief Rahman Hakim, Surabaya 60111 Indonesia
e-mail : djatmiko@me.its.ac.id
K
𝑇 𝑐,𝑜 = 𝑚̇ ℎ . 𝑐 𝑝ℎ . (𝑇 ℎ,𝑖 − 𝑇 ℎ,𝑜 ) 𝑚̇ 𝑐. 𝑐 𝑝𝑐 + 𝑇 𝑐,𝑖 (4)
dimana: 𝑚̇ ℎ = laju alir massa fluida panas (kg/s) 𝑚̇ 𝑐 = laju alir massa fluida dingin (kg/s)
𝐶 𝑝,ℎ = panas spesifik fluida panas (kJ/kg K)
𝐶 𝑝,𝑐 = panas spesifik fluida dingin (kJ/kg K)
𝑇 ℎ,𝑖 = temperature masuk dari fluida panas (K)
𝑇 ℎ,𝑜 = temperature keluar dari fluida panas (K)
𝑇 𝑐,𝑖 = temperature masuk dari fluida dingin (K)
𝑇 𝑐,𝑜 = temperature keluar dari fluida dingin (K) Persamaan  untuk  mendapatkan  nilai  UA  berdasarkan  nilai ∆T LMTD adalah sebagai berikut:
∆𝑇 𝐿𝑀𝑇𝐷 = (𝑇 ℎ,𝑖 − 𝑇 𝑐,𝑜 )(𝑇 ℎ,𝑜 − 𝑇 𝑐,𝑖 ) (𝑙𝑛 (𝑇 ℎ,𝑖 − 𝑇 𝑐,𝑜 ) (𝑇 ℎ,𝑜 − 𝑇 𝑐,𝑖 )) (5)
Nilai F dapat dicari dengan menggunakan grafik (Lihat Gambar 3)
Nilai  UA  digunakan  untuk  memperhitungkan  efektivitas dari alat penukar panas.
𝑈𝐴 = 𝑞 𝐹. ∆𝑇 𝐿𝑀𝑇𝐷 (6)
dimana: 𝑞 = laju perpindahan panas, (kW) 𝑈 = overall heat transfer coefficient , (kJ/s.m 2 .K) A = luas bidang perpindahan panas (m 2 ) 𝐹 = faktor koreksi
Perancangan  dimensi  ini  disesuaikan  dengan  nilai  UA LMTD , dan  dari  data  spesifikasi  desain  permukaan plate  fin  heat exchanger yang  ada  di  buku  Kays  and  London,  adalah sebagai berikut:
𝑉 HE = 𝑝. 𝑙. 𝑡 (7)
dimana:
𝑉 HE = volume heat exchanger , (m 3 ) 𝑝 = panjang, (m)
𝑙 = lebar, (m)
𝑡 = tinggi, (m)
Persamaan  untuk  menghitung  nilai  luasan frontal dari  satu sisi heat exchanger adalah sebagai berikut:
𝐴 𝑓𝑟 = 𝑝. 𝑙 (8)
dimana:
𝐴 𝑓𝑟 = luas frontal dari penukar panas, (m 2 ) Persamaan untuk  mengitung  nilai  rasio  perpindahan  panas dari alat penukar panas adalah sebagai berikut:
α = 𝐴 𝑉 total (9) = 𝑏. 𝛽 𝑏 1 + 𝑏 2 + 2𝑎 (10)
dimana: α = rasio perpindahan daerah dari satu sisi untuk total volume  dari  alat  penukar  panas  (hanya  diberikan untuk crossed-rod matrices , tubular and finned and type surface , m 2 /m 3 ) 𝑏 1 = 𝑏 2 = jarak plate (m) 𝛽 =  perbanding  antara  luas  perpindahan  panas  total terhadap volume antar plat (m 2 /m 3 ) 𝑎 = tebal plate (m)
Gambar 3. Faktor koreksi untuk cross flow heat exchanger dengan kedua fluida tidak bercampur (incropera, 2007)
Persamaan  untuk  menghitung  nilai  luas  area  bebas  dari sebuah alat penukar panas adalah sebagai berikut:
𝐴 𝑓𝑓 = 𝜎. 𝐴 𝑓𝑟,ℎ dimana, 𝜎 = 𝛼. 𝑟 ℎ (12) dimana:
A ff = luas area bebas, (m 2 ) δ    = ketebalan fin , (m) 𝑟 ℎ = jari-jari hidraulik → 𝑟 ℎ = 𝐴𝑐.𝐿 𝐴 (m)               (13) 𝐿 = panjang aliran di satu sisi (m)
Persamaan  untuk  menghitung  nilai  kecepatan  maksimum fluida dalam sebuah heat exchanger , adalah:
𝐺 = 𝑚̇ 𝐴 𝑓𝑓 (14)
dimana:
G = kecepatan  maksimum, (kg/m 2 .s)
Nilai reynolds digunakan untuk mengetahui jenis aliran yang ada dalam alat penukar panas. Adapun persamaannya adalah sebagai berikut:
𝑅𝑒 = 4. 𝑟 ℎ 𝜇 . 𝐺                                                                             (15)
dimana:
𝑅𝑒 = bilangan reynolds 𝐷 ℎ = diameter hidrolis, (m) 𝜇 = viskositas absolut udara, (N s/m 2 )
Persamaan  untuk  menghitung  bilangan  Prandalt  adalah sebagai berikut:
𝑃𝑟 = 𝑐 𝑝 . 𝜇 𝑘 (16)
dimana:
𝑃𝑟 = bilangan prandalt 𝑘      = konduktivitas termal fluida, (W/m K) Persamaan untuk menghitung bilangan nusselt adalah sebagai berikut:
𝑁𝑢 ℎ = 0,023 𝑅𝑒 ℎ 4
5 𝑃𝑟 0,3 (17)
𝑁𝑢 𝑐 = 0,023 𝑅𝑒 𝑐 4/5 𝑃𝑟 0,4 (18) dimana:
𝑁𝑢 = bilangan nusselt
Persamaan untuk menghitung nilai koefisien konveksi adalah sebagai berikut:
ℎ = 𝑁𝑢. 𝑘 𝐷 ℎ (19)
Persamaan  untuk  menghitung  effisiensi  dari  alat  penukar panas adalah sebagai berikut:
𝑚 = √ 2. ℎ 𝑘. 𝛿 (20)
dimana: m = fin effectiveness parameter , (/m) h  = koefisien konveksi , (kW/m 2 .K)
Gambar 4. Basic components of a plate fin heat exchanger (Shah and
Webb, 1983)
Gambar 5. Skema proses pendinginan motor listrik
𝑙 ≈ 𝑏 2 (21)
dimana: l = panjang fin efektif (m)
Maka,  persamaan  untuk  menghitung  effisiensi  fin  pada  alat penukar panas:
𝜂 𝑓 = 𝑡𝑎𝑛ℎ𝑚 𝑙 𝑚𝑙 (22)
dimana: 𝜂 𝑓 = efisiensi fin
𝜂 𝑜 = 1 − 𝐴 𝑓 𝐴 (1 − 𝜂 𝑓 )                                                             (23)
dimana: 𝜂 𝑜 = effisiensi overall 𝐴𝑓 𝐴 = diperoleh dari tabel 9.3 buku Kays and London Sehingga, UA dapat dihitung dengan persamaan berikut:
1 𝑈𝐴 = 1 (𝜂 𝑜 . ℎ. 𝐴) ℎ + 1 (𝜂 𝑜 . ℎ. 𝐴) 𝑐 (24)
Kapasitas panas ( C ) merupakan kemampuan dari alat penukar panas untuk menyerap panas. Adapun persamaannya adalah sebagai berikut:
𝐶 ℎ = 𝑚̇ ℎ . 𝑐 𝑝,ℎ (25) 𝐶 𝑐 = 𝑚̇ 𝑐 . 𝑐 𝑝,𝑐 (26) dimana:
𝐶 ℎ = kapasitas panas fluida panas, (kJ/s. K)
𝐶 𝑐 = kapasitas panas fluida dingin, (kJ/s. K)
Gambar 6. Skema sistem heat exchanger
Gambar 7. Skema sistem heat exchanger dalam sebuah ruang
Keterangan gambar: 1  : Udara dingin masuk ( T c,i )
2  : Blower
3  : Heat Exchanger
4  : Blower
5  : Motor listrik
6  : Udara keluar ( T c,o ) 7  : Ducting
Dalam  sebuah  perancangan  alat  penukar  panas  biasanya menggunakan perhitungan NTU ( Number of Transfer Unit ). Persamaannya adalah sebagai berikut:
NTU = 𝑈𝐴 𝐶 𝑚𝑖𝑛 (27)
dimana:
C min = kapasitas panas diantara dua fluida, (kJ/s. K)
Effectiveness ( ε ) dapat digunakan untuk menentukan unjuk kerja  dari  sebuah heat  exchanger .  Adapun  persamaannya adalah sebagai berikut:
𝜀 𝑡𝑒𝑜𝑟𝑖 = 1 − 𝑒𝑥𝑝 [( 1 𝐶 𝑟 ) . (𝑁𝑇𝑈) 0,22 . {𝑒𝑥𝑝[−𝐶 𝑟 (𝑁𝑇𝑈) 0,78 ]
−    1}]                                                                                     (28)
Gambar  4  merupakan  bagian-bagian  dari plate  fin  heat exchanger
Gambar 8. Winding motor Gambar 9. Distribusi temperature

## II. METODE PENELITIAN
Pada penelitian ini dilakukan dalam tiga tahap. Tahap awal dari penelitian ini adalah dilakukan studi literatur mengenai jenis-jenis alat penukar panas yang berasal dari buku, jurnal, maupun penelitian terdahulu.
Tahap  yang  kedua  adalah  penentuan  data  teknis  yang berasal dari data operasi motor yang berupa W out motor = 5683 kW, arus listrik= 386 Ampere, voltage= 10 kV dan cos Ø= 0,86. Serta T ambient sekitar plant = 34° C dan T ambient di luat plant = 30° C Tahap  ketiga  merupakan  tahap  perhitungan  matematis. Setelah  didapatkan  data  operasi  motor  dan  data  pendukung lain,  dilakukan  perencanaan  dimensi  yang  sesuai  dengan keinginan kemudian dianalisa performa dari heat exchanger hasil  perancangan  dengan  hasil  teori  yang  ada.  Hasil performa  akan  divariasikan  berdasarkan  perubahan  beban yang harus didinginkan ( Q = 50 %, 75 %, 100 % dan 105 %). Sistem pendinginan motor listrik bermula dari udara  yang dialirkan  secara  konveksi  paksa  melalui fan/blower. Lalu udara panas yang keluar dari motor listrik akan didinginkan dengan udara yang dilewatkan heat exchanger , supaya udara tersebut  dapat  digunakan  lagi  untuk  mendinginkan  motor listrik.  Pada  proses  pendinginan  motor  listrik  ini  digunakan sistem tertutup.

## III. ANALISA DAN PEMBAHASAN
A. Skema sistem heat exchanger yang akan dirancang (Lihat Gambar 6):
Mula-mula udara bebas di luar plant akan dialirkan secara konveksi  melalui blower masuk  kedalam  heat  exchanger, kemudian  udara  panas  yang  keluar  dari  motor  akan didinginkan dengan udara bebas tadi. Proses pendinginan ini dilakukan dengan sistem tertutup. Pemasangan blower kedua dekat dengan winding dari motor untuk mencegah udara tidak bias masuk ke dalam motor.
G ambar 10. Konstruksi plate fin heat exchanger (Kays and London, 1964)
Gambar 11. Grafik beban fungsi temperatur udara

## B. Perhitungan laju massa alir udara (𝑚̇)
Setelah  dilakukan  pengukuran  kecepatan  udara  yang masuk  ke  dalam winding dari  motor  diperoleh  kecepatan udara pada masing-masing area sebesar v 1 = 14 m/s; v 2 = 10,5 m/s; v 3 = 14,7 m/s. Kemudian kita akan menghitung nilai dari laju  alir  massa  udara  yang  melewati  masing-masing  laluan, sehingga diperoleh nilai laju alir massa udara total. Laju alir massa  udara  total  inilah  yang  digunakan  untuk  menghitung temperatur  udara  yang  keluar  dari  motor.  Perhitungan temperatur udara yan masuk dan keluar dari heat exchanger (Gambar 9) C. Perhitungan dimensi heat exchanger yang sesuai dengan data  operasi  pembangkit  listrik  tenaga  uap  dan spesifikasi untuk plate fin heat exchanger buku Kays and London (Gambar 10)
Setelah dilakukan perhitungan matematis diperoleh dimensi heat exchanger yang sesuai dengan data operasi pembangkit listrik tenaga uap adalah sebagai berikut:
1. Panjang HE   : 1,556 m
2. Lebar HE     : 0,897 m
3. Tinggi HE    : 1,299 m
4. Tebal plat ( a )  : 0,001 m
Berdasarkan  perhitungan  temperatur  yang  masuk  dan keluar  dari heat  exchanger ,  dapat  dibuat  grafik  temperatur fungsi variasi beban yang harus didinginkan ( Q = 50 %, 75 %, 100 % dan 105 %).
Variasi  pembebanan  ini  akan  meningkat  seiring  dengan meningkatnya temperatur udara panas yang masuk ke dalam Heat Exchanger . Begitu juga dengan temperatur udara dingin yang  keluar  dari heat  exchanger akan  meningkat.    Sesuai denan  perumusan,  apabila  beban  ( q )  meningkat,  untuk nilai 𝑚̇ ℎ , c p,h dan 𝑇 ℎ,𝑜 yang  tetap,  maka  nilai 𝑇 ℎ,𝑖𝑛 akan meningkat.  Begitu  pula  dengan  perubahan  beban  yang divariasikan  untuk  fluida  dingin,  dimana  saat  beban  ( q ) meningkat,  untuk  nilai 𝑚̇ 𝑐 , 𝑐 𝑝,𝑐 dan 𝑇 𝑐,𝑖𝑛 yang  tetap,  maka nilai 𝑇 𝑐,𝑜 juga akan meningkat.
300 315 330 345 360 0 50 100 150 200 T em pera tur (K ) Beban yang harus didinginkan (kW) Grafik beban vs temperatur Temperatur udara panas Temperatur udara dingin
Gambar 12. Grafik beban fungsi effectiveness
Apabila nilai beban ( q ) meningkat maka nilai effectiveness (ɛ)  akan  meningkat  juga.    Sesuai  perumusan  yang  ada = 𝑈𝐴. ∆𝑇 𝐿𝑀𝑇𝐷 , untuk setiap kenaikan q maka nilai dari ∆𝑇 𝐿𝑀𝑇𝐷 juga  akan  meningkat  (seperti  pembahasan  Grafik  4.11), sehingga nilai UA juga akan meningkat. Dari NTU = 𝑈𝐴 𝐶𝑚𝑖𝑛 , untuk nilai 𝐶 𝑚𝑖𝑛 yang tetap apabila nilai UA bertambah besar maka  nilai  NTU  juga  akan  bertambah  besar  juga.  Nilai effectiveness (ɛ)  merupakan  fungsi  dari  NTU  dan Cr ,  jika nilai NTU meningkat maka nilai effectiveness (ɛ) juga akan meningkat.  Sehingga  dapat  disimpulkan  bahwa  semakin meningkat nilai pembebanan ( q ), maka nilai effectiveness (ɛ) juga akan semakin meningkat. Dapat dilihat saat beban yang didinginkan  105  %,  maka  diperoleh  nilai  efektivitasnya adalah 0,8845.

## IV. KESIMPULAN
Hasil penelitian yang dilakukan, kemudian perhitungan  dan  analisa  performa heat  exchanger ,  maka penulis dapat memberikan beberapa kesimpulan, diantaranya:
1. Dengan beban pendinginan sebesar 167,5 kW konstruksi heat exchanger yang sesuai adalah sesuai dengan surface designation 5.3  pada  buku  Kays  and  London  dengan panjang HE 1.556 m, lebar HE 0.897 m dan tinggi HE 1.229 m yang memiliki spesifikasi sebagai berikut:
a. Material heat exchanger adalah aluminium b. Dengan Plate  spacing ( b )  =  0.0119  m, radius hydraulic ( r h )  =  0.0015  m,  ketebalan fin ( δ )  = 0.00015  m, β =  617.1449  m 2 /m 3 , fin  area/total  fin ( Af/A ) = 0.719
c. Dan Plate thickness ( a ) = 0.001 m
2. Performa heat  exchanger dapat  ditinjau  dari  nilai effectiveness. Semakin  tinggi  beban  pendinginan  pada heat  exchanger maka  nilai effectiveness juga  semakin tinggi.  Hal  ini  terlihat  pada  beban  yang  harus didinginkan 105 % dengan efektivitas 0,8845.

## UCAPAN TERIMA KASIH
Penulis Nur Sapti Marsheliyana mengucapkan terima kasih kepada Bapak Prof. Dr. Ir. Djatmiko Ichsani, M.Eng, selaku dosen  pembimbing  tugas  akhir,  yang  telah  membimbing dalam mata kuliah tugas akhir.
Tabel 1. Surface Geometry, Plate Fin Surface dari tabel 9.3 “ Compact HE , Kays-London”
Sisi  Panas  dan Sisi Dingin Surface designation 5.3 Plate spacing , b (mm) 11,938 Hydraulic radius , r h (mm) 1,537 Fin thickness , δ (mm) 0,152 Transfer  area/volume  betwees  plates ,β, m²/m³ 617,1449 Fin or extended area/total area 0,719 Tabel 2. Tabel NTU dan effectiveness % Beban (kW) UA (kW/K) NTU ɛ 50 % 83,750 8,8619 1,7009 0,6853 75 % 125,625 13,0229 2,5515 0,7821 100 % 167,5 17,3638 3,4019 0,8369 105 % 175,875 18,2319 3,5702 0,8845

## DAFTAR PUSTAKA
[1] Fox,  Robert  W.  Pritchard,  Philip  J.  McDonald,  Alan  T.  2003. Introduction To Fluid Mechanics. John Wiley & Sons (Asia) Pte Ltd: Asia
[2] Incropera,  Frank  P.  Dewitt,  David  P.  Bergman,  Theodore  L.  Lavine, Adrienne  S.  2007. Fundamentals  of  Heat  and  Mass  Transfer .  John Wiley & Sons (Asia) Pte Ltd: Hoboken.
[3] International  site  for  spirax  sarco.  2013: Steam  Consumtion  of  Pipes and Air Heater , <URL: http://www.spiraxsarco.com/resources/steam- engineering-tutoria/steam-engineering-principle-and-heat- transfer/steam-consumtion-of pipes and air heater.asp>.
[4] Maiti,  D.K  and  Sarangi  S.  K. Heat  Transfer  and  Flow  Friction Characteristic  of  Plate  Fin  Heat  Exchanger  surface- A  Numerical Study  P,hD.  Disertation,  Indian  Institut  of  Technologi  Kharagpur (2002)
[5] Moran,  Michael  J.  Shapiro,  Howard  N.  2004. Fundamentals  of Engineering  Thermodynamics .  John Wiley  &  Sons  Ltd,  The  Atrium, Shouthern Gate, Chichester, West Sussex PO19 8SQ: England
[6] Natasha  fitri.  Susanti.  2010. Analisis  Unjuk  Kerja  Alat  Pendingin Oksigen Cair (HE 38) dengan Variasi Desain Permukaan Studi Kasus pada Unit Separasi Udara PT. Gresik Gases Indonesia (BOC-Gases) .
Surabaya: Teknik Mesin
[7] Wen  Kang,  Shung  &  Chau  Tseng.Shin. Chuan Rod.  Ying. 5-6  April 2007.  “ Analysis of Effectiveness  and Pressure  Drop  in Micro Cross- Flow  Heat  Exchanger ”,  Applied  Thermal  Engineering;  Volume  27: 877-885
[8] W.M.  Kays  &  London.A.L.  1964. Compact  Heat  Exchanger  second edition . Mc Graw Hill: United States of America. 0.00 0.30 0.60 0.90 0 30 60 90 120 150 180 210 e ff ec tiveness (ɛ ) Beban yang harus didinginkan (kW) Grafik beban vs effectiveness