## OPTIMALISASI RONGGA TERHADAP VARIASI DERAJAT KEVAKUMAN SEBAGAI ISOLATOR

## MULYONO
Jurusan Teknik Mesin Fakultas Teknik Universitas Muhammadiyah Malang
E-mail: mail: mulyono@umm.ac.id

## AbSTrAK
Material isolasi sering digunakan untuk mempertahankan temperatur zat pada keadaan dingin dan panas. Material isolasi terkadang dipandang kurang praktis karena mahal dan tebal. Material isolasi dapat digantikan dengan cara Material isolasi dapat digantikan dengan cara menggunakan rongga udara bertekanan rendah atau vakum. Membuat derajat kevakuman 100% pada rongga udara tidaklah mudah. Permasalahan yang timbul dalam penelitian ini adalah bagaimana laju perpindahan panas konveksi alamiah yang terjadi pada variasi derajat kevakuman dan aspek ratio rongga yang berbeda agar diperoleh aspek ratio rongga yang optimal dalam penyimpanan. Tujuan penelitian ini adalah untuk mendapatkan aspek rongga dan besarnya tekanan vakum yang paling optimal dan mampu sebagai isolator dalam suatu rongga. Manfaat dari penelitian ini adalah dipergunakan untuk menyimpan bahan/zat/lainnya yang tidak tahan pada kondisi temperatur/tekanan atmosfir (bahan Kriogenik). Dari penelitian ini diperoleh derajat kevakuman –60 cm Hg dan L/ d = 5.96 dengan harga laju perpindahan panas konveksi alamiah lebih rendah dibandingkan dengan derajat kevakuman –20 cm Hg dan –40 cm Hg, sehingga mempunyai arti bahwa derajat kevakuman –60 cm Hg mempunyai hambatan termal lebih besar (isolator) dibandingkan dengan derajat kevakuman –20 cm Hg dan –40 cm Hg.
Kata kunci : derajat kevakuman, rongga, perpindahan panas, aspek ratio rongga

## AbSTrACT
Insulation materials are often utilized to maintain the temperature, either low or high. However, since the insulation materials often need to be thick, thus, costly, they are frequently deemed impractical. Therefore, in the attempt to replace the insulation materials, a cavity with low vacuum pressure is opted for. The purpose of the study is to find answers to the following question: "How do vacuum pressure variation and cavity ratio affect the rate of natural heat convection through the a cavity?" Pertinent to the question, the study was aimed to find the appropriate value of the vacuum pressure which can function well as an insulator. This study is significant in the attempt to lower down the rate of heat transfer taking place in a system vis-à-vis the surrounding media. The study found out that the degree of emptiness of -60 cm Hg and L/ L/ d = 5.96, = 5.96, results in a lower rate of heat transfer compared with -20 cm Hg and –40 cm Hg. This means that the vacuum pressure of –60 cm Hg bears a bigger thermal resistance than the -20 cm Hg and –40 cm Hg do.
Keywords : vacuum pressure, cavity, heat transfer, aspects of cavity ratio

## PENDAHULUAN
Perkembangan teknologi sekarang ini, untuk menyimpan zat-zat cair kriogenik yang bersuhu rendah (sampai kira-kira –250° C), misalnya hidrogen cair untuk waktu yang lama, telah dikembangkannya superisolator (superinsolation). Superisolator yang Superisolator yang paling efektif terdiri dari lapisan-lapisan rangkap yang  terbuat  dari  bahan  yang  berdaya  refleksi tinggi dengan isolator-isolator sebagai pengantara. Keseluruhan sistem ini dihampakan agar konduksi melalui  udara  menjadi  minimum  (Wang,  2002). Material isolasi digunakan untuk mempertahankan temperatur zat pada keadaan dingin dan panas. Material  isolasi  terkadang  dipandang  kurang praktis karena mahal dan tebal. Material isolasi
dapat  digantikan  dengan  menggunakan  rongga udara bertekanan rendah atau vakum. Membuat derajat kevakuman 100% pada rongga udara tidaklah mudah, biasanya masih ada udara didalam rongga yang mengakibatkan terjadinya perpindahan panas konveksi alami (natural convection) antara kedua permukaan pembentuk rongga tersebut (Holman, 1988).
Laju  perpindahan  panas  pada  suatu  rongga dihitung berdasarkan hukum newton pendinginan. Koefisien perpindahan panas konveksi(h) pada rongga dipengaruhi oleh perbedaan temperatur dari kedua dinding pembatas, selain itu juga merupakan fungsi dari geometri rongga, orientasi dari rongga dan sifat- sifat yang dimiliki fluida. Besarnya laju perpindahan
panas  antara  dua  permukaan  yang  membentuk rongga sangat dipengaruhi pula oleh kondisi sifat-sifat (properties) yaitu; tekanan, temperatur, massa jenis, konduktivitas, viskositas dan sebagainya dari fluida yang berada di dalam rongga tersebut (Setterfield, 2002). Laju perpindahan panas pada suatu rongga (cavity) dapat dikurangi dengan cara penghampaan (pemvakuman). Penghampaan merupakan suatu metode  yang  sangat  populer  didalam  kehidupan sehari-hari. Cara ini sering digunakan pada thermos, pelat-pelat absorber pada pemanas matahari, sebagai isolasi untuk mengurangi kehilangan panas pada reaktor nuklir, pendinginan pada tangki sampah radioaktif, ventilasi ruangan dan seterusnya (Wang, 2002).
Permasalahan yang timbul dalam penelitian ini adalah; seberapa besar laju perpindahan panas konveksi alamiah yang terjadi pada rongga vakum dan  aspek  ratio  berapa  yang  memberikan  daya penyimpanan  yang  optimal.  Tujuan  penelitian ini adalah untuk mendapatkan aspek rongga dan besarnya  tekanan  vakum  yang  paling  optimal dan mampu sebagai isolator dalam suatu rongga. Manfaat dari penelitian ini adalah dipergunakan untuk menyimpan bahan/zat/lainnya yang tidak tahan pada kondisi temperatur/tekanan atmosfir (bahan Kriogenik).

## METODE
Model eksperimen akan dijelaskan pada Gambar 1 berikut.
perpindahan  panas  konveksi  alami (natural  convection) antara  kedua  permukaan  pembentuk  rongga tersebut (Holman JP, 1988) .
Laju  perpindahan  panas  pada  suatu  rongga  dihitung  berdasarkan  hukum  newton  pendinginan. Koefisien perpindahan panas konveksi(h) pada rongga dipengaruhi oleh perbedaan temperatur dari kedua dinding pembatas, selain itu juga merupakan fungsi dari  geometri rongga, orientasi dari rongga dan sifat- sifat  yang  dimiliki  fluida. Besarnya  laju  perpindahan  panas  antara  dua  permukaan yang  membentuk rongga  sangat dipengaruhi pula oleh  kondisi  sifat-sifat (properties) yaitu ; tekanan, temperatur, massa jenis,  konduktivitas,  viskositas  dan  sebagainya  dari  fluida  yang  berada  di  dalam  rongga  tersebut (Setterfield Barry, 2002) . Laju perpindahan panas pada suatu rongga (cavity) dapat dikurangi dengan cara penghampaan  (pemvakuman).  Penghampaan merupakan  suatu  metode  yang  sangat  populer  didalam kehidupan  sehari-hari. Cara  ini  sering  digunakan  pada thermos,  pelat-pelat  absorber  pada  pemanas matahari,  sebagai  isolasi  untuk  mengurangi  kehilangan  panas  pada  reaktor  nuklir,  pendinginan  pada tangki sampah radioaktif, ventilasi ruangan dan seterusnya (Wang Q.W, 2002).
Permasalahan  yang  timbul  dalam  penelitian  ini  adalah; seberapa  besar    laju  perpindahan  panas konveksi  alamiah  yang  terjadi  pada  rongga  vakum  dan  aspek  ratio  berapa  yang  memberikan  daya penyimpanan yang optimal. Tujuan penelitian ini adalah untuk mendapatkan aspek rongga dan  besarnya tekanan  vakum  yang  paling  optimal  dan    mampu    sebagai  isolator  dalam  suatu  rongga. Manfaat    dari penelitian  ini  adalah  dipergunakan  untuk  menyimpan  bahan/zat/lainnya  yang  tidak  tahan  pada  kondisi temperatur/tekanan atmosfir (bahan Kriogenik).
.

## METODE
Model eksperimen akan dijelaskan pada Gambar 1 berikut.
Analisa proses perpindahan panas konveksi alamiah dalam tabung annulus yang diberikan perlakuan berupa  variasi  derajat  kevakuman,  dilakukan  analisa  dimensional  dengan  pendekatan “DALIL BUCKINGHAM” atau Buckingham  method ,  yang  mana  akan  didapatkan  suatu  fungsi-fungsi  tanpa dimensi. Persamaan  fungsi  koeffisien  perpindahan  panas  konveksi  alamiah  yang  menyatakan  hubungan antara variabel-variabel adalah sebagai berikut ;
h  = f { , k, Cp, , , p, f s T T g } atau Pompa Vacuum P Aquisi data Monitor Computer Katup 300 mm r w1 80 mm Gambar 1.  Model Eksperimen r w2
Gambar 1. Model eksperimen
Analisis  proses  perpindahan  panas  konveksi alamiah  dalam  tabung  annulus  yang  diberikan perlakuan  berupa  variasi  derajat  kevakuman, dilakukan Analisis dimensional dengan pendekatan "DALIL p BUCKINGHAM"  atau  Buckingham method, yang mana akan didapatkan suatu fungsi- fungsi tanpa dimensi. Persamaan fungsi koeffisien
perpindahan  panas  konveksi  alamiah  yang menyatakan  hubungan  antara  variabel-variabel adalah sebagai berikut.
h = f {d, k, Cp, m , r , p, b g (T – Tf)}  .................... (1)
atau
f {h, d, k, Cp, m, r , p, b g (T s – Tj)}= 0 ................ (2)
Buckingham methode dengan cara menggunakan dimensi  MLtT  (massa,  panjang,  waktu  dan temperatur) akan didapat grup parameter tanpa dimensi dari fungsi diatas, yaitu:
p 1 = f (h, d a1 , k b1 , m cl , p dl , r el ) » bilangan Nussel
p 2 = f (C p , d a2 , k b2 , m c2 , p d2 , r e2 ) » bilangan Prandtl
p 3 = f ( b g, (T s –T f ), d a3 , k b4 , m c4 , p d4 , r e4 ) » bil. Grashoft
p 4 = ( L d ) » aspek ratio rongga ........................ (3)
sehingga  didapatkan  persamaan  koeffisien perpindahan panas konveksi alamiah dalam bentuk variabel tanpa dimensi sebagai berikut:
p 1 = f ( p 2 , p 3 , p 4 ), atau:
p 4 = f ( Pr, Gr s L d ) atau:
Nu = f ( Ra s, L d ) .................................................. (4)
Analisis dimensi diperoleh koeffisien perpindahan panas (variabel tetap) konveksi alamiah sebagai fungsi dari bilangan Rayleigh, aspek ratio rongga. Eksperimen  ini  dilakukan  untuk  mendapatkan hubungan variabel antara koeffisien perpindahan panas konveksi alamiah (bilangan Nusselt) terhadap bilangan  Prandtl,  bilangan  Grashoft  dengan melakukan  perlakuan  berupa  variasi  derajat kevakuman  (tekanan  rendah)  dan  aspek  ratio diameter tabung/rongga.
Berdasarkan  Gambar  1,  peralatan  yang digunakan dalam eksperimen ini adalah sebagai berikut: 1) Barometer vacuum, untuk menggukur
tekanan  vacuum;  2)  Thermometer  biasa,  untuk mengukur  temperatur  ruangan;  3)  Tabung anulus  dari  bahan  baja  carbon  chrom,  sebagai objek penelitian; 4) Thermokopel/sensor tipe lm- Thermokopel/sensor tipe lm- 35, untuk mengukur temperatur dinding tabung; 5) Pompa vakum (yang digunakan untuk AC), untuk memvakumkan tabung annulus/rongga; 6) Akuisisi data, difungsikan untuk merekam input data dari sensor atau termokopel yang merubah input panas menjadi volt selanjutnya menjadi data temperature; 7) Seperangkat komputer (PC), untuk menampilkan data hasil proses dari data aquisisi; 8) Katup, untuk membuka dan menutup aliran udara ke rongga; 9) Fluida kerja, tabung dalam berupa es dan tabung Fluida kerja, tabung dalam berupa es dan tabung luar berupa udara. Mengukur temperatur sejumlah titik pada dinding tabung digunakan thermokopel tipe lm-35 yang dihubungkan dengan data akuisisi yang telah dikalibrasi, selanjutnya dihubungkan ke  seperangkat  komputer  untuk  membaca  data temperatur.
Berdasarkan Gambar 1 yaitu skema penelitian, maka untuk melakukan penelitian dan pengambilan data, langkah-langkahnya adalah sebagai berikut: 1) Mempersiapkan tabung annulus, thermokople, barometer,  es,  pompa  vakum  dan  komputer selanjutnya alat diset; 2) Operasikan pompa vakum Operasikan pompa vakum hingga  mencapai  tekanan  vakum  yang  telah ditentukan; 3) Bila tekanan vakum (tekanan yang direncanakan) telah tercapai, selanjutnya katup ditutup agar tidak berhubungan dengan lingkungan; 4)  Mencatat  seluruh  data  yaitu:  temperatur dinding, tekanan vakum, temperatur udara vakum, temperatur es, temperatur lingkungan semuanya dicatat  oleh  data  aquisisi  yang  dihubungkan  ke computer; 5) Point 3 diulang dengan tekanan vacuum dengan tekanan vacuum yang lain dan seterusnya.

## HASIL DAN PEMBAHASAN
Gambar  2  menunjukkan  hubungan  antara derajat/tingkat  kevakuman  terhadap  bilangan Nusselt adanya variasi aspek ratio rongga  L d pada posisi  90  derajat  dan  temperatur  permukaan: 45° C. Pada Gambar 2 dapat dijelaskan bahwa tingkat/ derajat kevakumanan semakin besar maka harga
bilangan Nusselt semakin menurun, karena semakin besar derajat kevakuman akan memberikan harga kerapatan jenis udara dalam rongga semakin besar. Kerapatan jenis semakin besar berarti pergerakan fluida semakin lambat, sehingga akan memengaruhi selisih  kerapatan  jenis  free  stream  dengan kerapatan jenis mula-mula ( r ¥ – r). Pada aspek ratio rongga  L d yang besar(rongga tipis) akan memberikan selisih kerapatan yang kecil sehingga mengakibatkan pergerakan fluida(gaya apung) dan harga bilangan Nusselt yang kecil. Pada Gambar 2 dengan diberikannya perlakuan kevakuman untuk aspek ratio rongga ( L d = 5.90) mempunyai harga bilangan Nusselt lebih kecil dibandingkan aspek ratio yang lain.
0 50000 100000 150000 200000 250000 300000 350000 10 20 30 40 50 60 70 R a Tingkat Kevakuman (cm Hg) L/dha= L/dha=7.87
L/dha=11.8
Gambar 3. Tingkat Kevakuman vs Rayleigh pada Posisi 90 derajat dengan Temperatur Permukaan 45° C

## Tabel 1. Variabel Perlakuan
Diameter tabung dalam (D1) … Es Diameter tabung luar (D2)... udara Temperatur permukaan (° C) Tekanan kevakuman (cm Hg) 50,8 mm 76,2 mm 101,6 mm 35 ° C 45 ° C 55 ° C 65 ° C -20 cm Hg -40 cm Hg -60 cm Hg 5 0 0 5 5 0 0 1 0 5 0 0 1 5 5 0 0 2 0 5 0 0 1 0 2 0 3 0 4 0 5 0 6 0 7 0 N u T e ka n an  K ev a ku m a n  ( cm  H g ) L /d h a = 5 .9 L /d h a = 7 .8 7 L /d h a = 1 1 .8
Gambar 2. Hubungan Tingkat Kevakuman terhadap Hubungan Tingkat Kevakuman terhadap Bilangan Nusselt
Gambar 3 merupakan grafik hubungan derajat kevakuman terhadap bilangan Rayleigh pada posisi 90 derajat dengan temperatur permukaan 45° C. Pada Gambar 3 dapat dijelaskan bahwa semakin besar  tingkat  kevakuman  pada  rongga,  harga bilangan Nusselt cenderung naik. Bagi aspek ratio rongga semakin kecil ( L d = 5.90 ) atau rongga besar mempunyai bilangan Rayleigh yang paling besar, artinya bila tekanan kevakuman dinaikkan maka kerapatan akan naik, sehingga kecepatan fluida/ udara dalam rongga menjadi sangat kecil. Kecepatan fluida  bergerak  dipengaruhi  oleh  perbedaan kerapatan  free  stream  dengan  fluida  di  dinding. Terjadinya perbedaan kerapatan disebabkan oleh kondisi temperatur permukaan dinding (T w2-in ).
Pada aspek ratio rongga semakin besar ( L d ) atau ketebalan rongga adalah tipis, maka harga bilangan Rayleigh semakin kecil, karena lapisan batas termal yang terbentuk adalah tipis sehingga perbedaan temperatur free stream dengan temperatur dinding (T w2–n – T ¥ ) adalah kecil. Akibat perbedaan temperatur yang kecil maka harga bilangan grashoft juga kecil, sehingga bilangan Rayleigh juga kecil/ rendah, karena bilangan Rayleigh merupakan fungsi dari bilangan Grashof dan Prandtl.
Guna memberikan satu argumentasi terhadap pengaruh derajat kevakuman dalam rongga, validasi yang digunakan adalah mengamati temperatur yang terjadi selama es mencair (waktu yang dibutuhkan). Gambar 4 menjelaskan seberapa besar waktu yang digunakan untuk mencair dengan adanya variasi derajat kevakuman. Tekanan kevakuman = –60 cmHg, sehingga mempunyai waktu yang lebih lama untuk es mencair, karena semakin tekanan kevakuman meningkat  tingkat  hambatan  termal  di  rongga semakin besar, sehingga arus temperatur menuju ruang es terhambat. Derajat kevakuman yang besar akan mengakibatkan perubahan kerapatan jenis menjadi besar dan kecepatan aliran fluida menjadi lambatakibatnya gaya apung yang terjadi dirongga lemah, artinya arus konveksi(gaya bodi) menjadi kecil.
-6 -4 -2 0 2 4 6 8 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 Waktu (detik)
Te m pe ra tu r ( C P=-20 cm hg P = -40 cm hg P = -60  cm hg
Gambar 4. Validasi Eksperimen
SIMPULAN Hasil penelitian yang telah dilakukan dapat disimpulkan bahwa suatu rongga yang diberi  perlakuan variasi derajat kevakuman akan mampu memberikan kontribusi sebagai isolator. Derajat kevakuman –60 cm Hg mempunyai laju perpindahan panas konveksi alamiah lebih rendah dibandingkan dengan derajat kevakuman  –20  cm  Hg  dan  –40  cm  Hg,  artinya  derajat  kevakuman  –60  cm  Hg  mempunyai  harga hambatan termal lebih besar (cocok sebagai isolator) khusunya dalam proses penyimpanan atau tingkat penyimpanannya lebih optimal.
Aspek ratio rongga  ( G L ) = 5.90 mempunyai harga bilangan Nusselt yang lebih kecil dibandingkan dengan  aspek  ratio  rongga  ( G L )  =  11.8,  artinya  aspek  ratio  rongga  5.90  mempunyai  kecepatan pergerakan molekul  membawa energi lebih rendah. Rongga yang kecil (aspek ratio yang besar = 11.8), kecepatan aliran udara di dalam rongga relatif lebih besar karena bilangan Nusselt berbanding langsung dengan bilangan Reynolds. Kecepatan  yang  tinggi  mengakibatkan  bilangan Reynolds tinggi,  akibatnya  bilangan Nusselt juga tinggi, sehingga nilai koefisien perpindahan panas konveksi alamiah (h) pada celah kecil adalah tinggi. Sehingga perpindahan panas konveksi pada celah kecil mempunyai nilai yang besar.
Hasil uji validasi eksperimen, dapat disampaikan bahwa semakin besar derajat kevakuman (-60 cm Hg) maka waktu mencair es lebih lama, bila dibandingkan dengan derajat kevakuman yang lebih rendah (-20 cm  Hg  atau  –40  cm  Hg).  Tekanan  kevakuman  –60  cm  Hg    mempunyai  hambatan  termal  lebih  besar artinya  perbedaan  temperatur free  stream (aliran  bebas)  dengan  temperatur  dinding  tabung    dalam  sisi luar adalah  kecil  sehingga waktu yang dibutuhkan untuk es mencair lebih lama.

## DAFTAR PUSTAKA
Brodkey, Robert S, (1998), Transport Phenomena , McGraw-Hill, Singapore. Chang S.W. dkk, (1999), An Eksperiment Study Of Heat Transfer In Reciprocating Square Duct Fitted with Ribs Skewed to the Flow , Journal of Heat Transfer, vol. 121 pp. 232-236. Campo A, Zamora B, (2000), Enhanced Natural Convection in a Vertical Rectangular Cavity on Account of the Mixing of Two Pure Gases , ASME. De Witt,  Incropera, (1996), Fundamentals of Heat  and Mass Transfer , John Willy & Sons, Singapore. Darling R B, (2002), Vacuum Systems , Journal Home page. Fox w, R. Mc Donald, Alan T, (1994), Introduction to Fluid Mechanics , John Willey & Sons, Singapore. Holman J P, (1988), Perpindahan kalor , Terjemahan Djasifi E, Erlangga, Jakarta. Laccarino  g,  Ooi  A,  (1998), Heat  Transfer  Predictions  in  Cavities ,  Center  for  Turbulent  Research, University of New South Wales, Australia.
Gambar 4. Validasi Eksperimen

## SIMPULAN
Hasil  penelitian  yang  telah  dilakukan  dapat disimpulkan  bahwa  suatu  rongga  yang  diberi perlakuan variasi derajat kevakuman akan mampu memberikan kontribusi sebagai isolator. Derajat kevakuman –60 cm Hg mempunyai laju perpindahan panas konveksi alamiah lebih rendah dibandingkan dengan derajat kevakuman –20 cm Hg dan –40 cm Hg, artinya derajat kevakuman –60 cm Hg mempunyai harga hambatan termal lebih besar (cocok sebagai isolator) khusunya dalam proses penyimpanan atau tingkat penyimpanannya lebih optimal. Aspek ratio rongga ( L d ) = 5.90 mempunyai harga bilangan Nusselt yang lebih kecil dibandingkan dengan aspek ratio rongga ( L d ) = 11.8, artinya aspek ratio rongga artinya aspek ratio rongga
5.90  mempunyai  kecepatan  pergerakan  molekul membawa energi lebih rendah. Rongga yang kecil (aspek ratio yang besar = 11.8), kecepatan aliran udara di dalam rongga relatif lebih besar karena bilangan  Nusselt  berbanding  langsung  dengan bilangan Reynolds. Hasil uji validasi eksperimen, dapat disampaikan bahwa semakin besar derajat kevakuman  (-60  cm  Hg)  maka  waktu  mencair es lebih lama, bila dibandingkan dengan derajat kevakuman yang lebih rendah (–20 cm Hg atau –40 cm Hg). Tekanan kevakuman –60 cm Hg mempunyai hambatan termal lebih besar artinya perbedaan temperatur  free  stream  (aliran  bebas)  dengan temperatur dinding tabung dalam sisi luar adalah kecil sehingga waktu yang dibutuhkan untuk es mencair lebih lama.

## DAFTAR PUSTAKA
Brodkey, and Robert S, 1998. Transport Phenomena,
McGraw-Hill, Singapore.
Chang S.W. dkk, 1999. An Eksperiment Study of Heat Transfer In Reciprocating Square Duct Fitted with Ribs Skewed to the Flow, Journal of Heat Transfer, vol. 121 pp. 232–236.
Campo  A,  and  Zamora  B,  2000.  Enhanced  Natural Convection  in  a  Vertical  Rectangular  Cavity  on Account of the Mixing of Two Pure Gases, ASME.
De Witt, and Incropera, 1996. Fundamentals of Heat and Mass Transfer, John Willy & Sons, Singapore.
Darling RB, 2002. Vacuum Systems, Journal Home page. Fox w, R. Mc Donald, and Alan T, 1994. Introduction to Fluid Mechanics, John Willey & Sons, Singapore. Holman J P, 1988. Perpindahan kalor, Terjemahan Djasifi E, Erlangga, Jakarta.
Laccarino g, Ooi A, 1998. Heat Transfer Predictions in Cavities, Center for Turbulent Research, University of New South Wales, Australia.
Oronzio Manca, dkk., 2002. Effect on Natural Convection of  the  distance  Between  an  Inclined  Discretely Heated Plate and a Parallel Shroud Below, Journal Heat Transfer, ASME.
Ousthuizen  H  Patrick,  and  David  Naylor,  1999.  An Introduction to Convective Heat Transfer Analysis, McGraw-Hill, Singapore.
Roth  A,  1989.  Vacuum  Technologi,  North  Holland, Amsterdam.
Setterfield Barry, (2002), Exploring the Vacuum, Journal of Theoritis, Journal home Page.
Wang Q.W. dkk., 2000. An Experiment Investigation of Natural convection in a Cubic Inclined Enclosure with Multiple Isolated Plates, Journal Heat Transfer, ASME.
Xundan Shi, dkk, 2003. Laminar Natural Convection Heat Transfer in a differentially heated square cavity due to a thin fin on the hot wall, Journal Heat Transfer, ASME.
Yarwood J, 1955. High Vacuum Technique, Chapman & Hall, London.