PENGUJIAN KARAKTERISTIK PERPINDAHAN PANAS DAN PENURUNAN TEKANAN DARI INLINE CYLINDRICAL PIN FIN ARRAY DALAM RECTANGULAR CHANNEL

PENGUJIAN KARAKTERISTIK PERPINDAHAN PANAS DAN PENURUNAN TEKANAN DARI INLINE CYLINDRICAL PIN FIN ARRAY DALAM RECTANGULAR CHANNEL

SKRIPSI Diajukan sebagai salah satu syarat untuk memperoleh gelar Sarjana Teknik

Oleh: MAHMUDIN HUDA NIM: I0404048

JURUSAN TEKNIK MESIN FAKULTAS TEKNIK UNIVERSITAS SEBELAS MARET SURAKARTA 2009

PENGUJIAN KARAKTERISTIK PERPINDAHAN PANAS DAN PENURUNAN TEKANAN DARI INLINE CYLINDRICAL PIN FIN ARRAY DALAM RECTANGULAR CHANNEL

Disusun oleh:

Mahmudin Huda NIM. I0404048 Dosen Pembimbing I

Dosen Pembimbing II

Tri Istanto, S.T., M.T. NIP. 19730820 200012 1 001

Wibawa Endra J., S.T., M.T. NIP. 19700911 200003 1 001

Telah dipertahankan di hadapan Tim Dosen Penguji pada hari Selasa tanggal 28 Juli 2009. 1. Ir. Santoso, M.Eng.Sc. NIP. 19450824 198012 1 001

.............................................

2. Syamsul Hadi, S.T., M.T. NIP. 19710615 199802 1 002

.............................................

3. Eko Prasetya Budiana, S.T., M.T. NIP. 19710926 199903 1 002

.............................................

Mengetahui: Ketua Jurusan Teknik Mesin

Koordinator Tugas Akhir

Dody Ariawan, S.T., M.T. NIP. 19730804 199903 1 003

Syamsul Hadi, S.T., M.T. NIP. 19710615 199802 1 002

PERSEMBAHAN Kepada mereka yang telah berjasa, kepada mereka pula aku persembahkan hasil jerih payahku selama menempuh jenjang S1 ini yaitu sebuah skripsi yang akan menjadi karya terbesarku selama ini sehingga aku lulus dari Universitas Sebelas Maret ini dengan gelar Sarjana Teknik. Mereka adalah: 1. Dengan nama-Nya yang telah menciptakan alam semesta dan seisinya. Segala puji bagi Allah, tidak ada daya dan upaya kecuali dengan-Nya. Allahlah pemilik segala keagungan, kemuliaan, kekuatan dan keperkasaan. 2. Bapak Sumadi, Ibu Siti Maryam, karena beliaulah penulis terlahir di dunia dengan kelebihan dan kekurangannya. 3. Adikku Khoirun Nisa dan kakakku Nur Qomarudin (kalian adalah harta yang tak ternilai). 4. Mr. 3G and Prof. Bawa, yang telah mensupport material (makan-makan), spiritual, dan membimbing tanpa rasa letih dan selalu ceria. 5. Semua orang yang dekat dan kenal dengan penulis (mereka yang pernah bersama memberi pengalaman yang berarti, memberikan nasehat serta dukungan dalam kehidupan penulis).

MOTTO

“Demi masa. Sesungguhnya manusia itu benar-benar dalam kerugian, Kecuali orang-orang yang beriman dan mengerjakan amal saleh dan nasehat menasehati supaya mentaati kebenaran dan nasehat menasehati supaya menetapi kesabaran.” (Q.S. Al-‘Ashr: 1-3) “Maka sesungguhnya bersama kesulitan pasti ada kemudahan, maka bersama kesulitan pasti ada kemudahan” (Q.S. Al-Insyirah: 5-6) “Kamu adalah umat terbaik yang dilahirkan untuk manusia, menyuruh yang ma’ruf dan mencegah dari yang munkar dan beriman kepada Allah.” (Ali Imran: 110) “Mukmin yang kuat lebih baik dan lebih dicintai Allah daripada mukmin yang lemah; pada keduanya ada kebajikan.” (HR. Muslim) “Cogito, ergo sum. (I think, therefore I am).” (Renè Descartes) “If i die tomorrow I'd be allright Because i believe That after we're gone The spirit carries on” (Dream Theater) “Kalau bisa dikerjakan sekarang, kenapa harus besok?” “Aku hidup di masa sekarang untuk menuju masa depan, bukan masa lalu.” “Jika dengan mengenang kegagalan bisa membuat hari esok lebih baik, kenanglah. Namun jika tidak, buanglah.” “Selalu berpikir ke depan dan optimis.” “Doa, usaha, dan tawakal.” (From the deepest of my heart)

Investigation on Heat Transfer and Friction Factor Characteristics of Inline Cylindrical Pin Fin Array in Rectangular Channel Mahmudin Huda Mechanical Engineering Departement Sebelas Maret University Surakarta, Indonesia E-mail: [email protected] Abstract Surface heat transfer extension using fins often used in heat exchanger equipment to enhance heat transfer. This research was conducted to investigate the characteristics of heat transfer and pressure drop as well as thermal performance of inline cylindrical pin fin assembly in rectangular channel. It was used cylindrical pin fins with 12.7 mm of diameter and 75 mm of height. The parameters of this research were Reynolds number (3,000 – 37,500), depending on the averaged inlet air velocity and hydraulic diameter, the interfin pitch distance in the spanwise direction (Sx/D) which was kept constant at 2.95 and the interfin pitch distance in the streamwise direction (Sy/D = 1.97 – 3.94). The experimental result shown that increasing Reynolds number (Re) and decreasing the interfin pitch distance in the streamwise direction (Sy/D) increased the average convection heat transfer coefficient (h) and Nusselt number (Nu). The pressure drop (∆P) and friction factor (f) decreased while the interfin pitch distance in the streamwise direction (Sy/D) increased. Increasing Reynolds number would decrease thermal performance (η). At Sy/D = 1.97 for Re > 24,900, and Sy/D > 1.97 for Re > 12,400, the values of η were less than 1, i.e. the used of pin fin assembly would cause an energy loss rather than gain. The net energy gain (η was greater than 1) was only at Sy/D = 1.97 for Re < 24,900, and at Sy/D > 1.97 for Re < 12,400. A net energy gain up to 28 % was achieved at Sy/D = 1.97 for Re = 3,076. Key words: pin fin, rectangular channel, Reynolds number, friction factor, thermal performance.

vi

Pengujian Karakteristik Perpindahan Panas dan Penurunan Tekanan dari Inline Cylindrical Pin Fin Array dalam Rectangular Channel Mahmudin Huda Jurusan Teknik Mesin Fakultas Teknik Universitas Sebelas Maret Surakarta, Indonesia E-mail: [email protected] Abstrak Perluasan permukaan perpindahan panas menggunakan sirip-sirip sering digunakan dalam peralatan penukar panas untuk meningkatkan perpindahan panas. Penelitian ini dilakukan untuk menguji karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal dari cylindrical pin fin assembly yang disusun secara segaris dalam saluran segiempat. Digunakan sirip-sirip pin silinder dengan diameter 12,7 mm dan tinggi 75 mm. Parameter-parameter dalam penelitian ini adalah bilangan Reynolds (3.000 – 37.500) berdasarkan kecepatan udara masuk rata-rata dan diameter hidrolik, jarak antar titik pusat sirip dalam arah melintang aliran udara (Sx/D) yang dibuat konstan sebesar 2,95 dan jarak antar titik pusat sirip dalam arah aliran udara (Sy/D = 1,97 – 3,94). Hasil pengujian menunjukkan bahwa peningkatan bilangan Reynolds (Re) dan semakin kecilnya jarak antar titik pusat sirip arah aliran udara (Sy/D) akan meningkatkan koefisien perpindahan panas konveksi rata-rata (h) dan bilangan Nusselt (Nu). Nilai penurunan tekanan (∆P) dan faktor gesekan (f) semakin menurun dengan meningkatnya nilai Sy/D. Meningkatnya bilangan Reynolds akan menurunkan unjuk kerja termal (η).Untuk Sy/D = 1,97 pada Re > 24.900, dan Sy/D > 1,97 pada Re > 12.400, nilai-nilai η lebih kecil dari 1, yang berarti bahwa pemakaian pin fin assembly akan menyebabkan kehilangan energi daripada perolehan energi. Perolehan energi netto (nilai η lebih besar dari 1) hanya untuk Sy/D = 1,97 pada Re < 24.900, dan Sy/D > 1,97 pada Re < 12.400. Perolehan energi netto dapat dicapai hingga 28 % untuk nilai Sy/D = 1,97 pada Re = 3.076. Kata kunci: sirip pin, saluran segiempat, bilangan Reynolds, faktor gesekan, unjuk kerja termal.

vii

KATA PENGANTAR Puji dan syukur penulis panjatkan ke hadirat Allah SWT atas segala limpahan rahmat dan Karunia-Nya sehingga penulis dapat melaksanakan dan menyelesaikan skripsi ”Pengujian Karakteristik Perpindahan Panas dan Penurunan Tekanan dari Inline Cylindrical Pin Fin Array dalam Rectangular Channel” ini dengan baik. Skripsi ini disusun guna memenuhi persyaratan untuk memperoleh gelar Sarjana Teknik di Jurusan Teknik Mesin Universitas Sebelas Maret Surakarta. Dalam penyelesaian skripsi ini tidaklah mungkin dapat terselesaikan tanpa bantuan dari berbagai pihak, baik secara langsung maupun tidak langsung. Oleh karena itu pada kesempatan ini penulis ingin menyampaikan rasa terimakasih yang sebesar-besarnya kepada semua pihak yang telah membantu dalam menyelesaikan skripsi ini, terutama kepada: 1. Bapak Dody Ariawan, S.T., M.T, selaku Ketua Jurusan Teknik Mesin UNS Surakarta. 2. Bapak Tri Istanto, S.T., M.T. selaku Pembimbing I atas bimbingannya hingga penulis dapat menyelesaikan skripsi ini. 3. Bapak Wibawa Endra J, S.T., M.T. selaku Pembimbing II yang telah turut serta memberikan bimbingan yang berharga bagi penulis. 4. Bapak Teguh Triyono, S.T., M.T. selaku Pembimbing Akademis yang telah berperan sebagai orang tua penulis dalam menyelesaikan studi di Universitas Sebelas Maret ini. 5. Bapak Syamsul Hadi, S.T,. M.T. selaku koordinator Tugas Akhir. 6. Bapak Budi Kristiawan, S.T., M.T. dan Ibu Eliza yang telah banyak membantu meminjami alat-alat pendukung penelitian. 7. Seluruh dosen serta staf di Jurusan Teknik Mesin UNS, yang telah turut mendidik penulis hingga menyelesaikan studi S1. 8. Ayah, Ibu, kakak, dan adikku, atas doa restu, motivasi, dan dukungan material maupun spiritual selama penyelesaian Tugas Akhir. 9. Rekan Skripsi: Apras, Adit, Doddy, Aji’, Teddy, Wisnu, mas Fendy, dan mas Giyono yang telah bersama-sama mengerjakan penelitian ini, terima kasih atas bantuan kalian semua.

viii

10. Rekan Asisten Lab. Konversi Energi: Ridho, Syafiq-jembe, Rian, paijo Bisyri, Ervan, Jihad-jae gepok, Tendy; dan juga rekan asisten lab 2005: Zaki, Yusno, duet Indri-Topan, Ahmad yang telah menemani, merepotkan dan direpotkan dalam pembuatan alat penelitian ini. Terima kasih yang tak terkira atas bantuan kalian semua. 11. Rekan-rekan Teknik Mesin semua, khususnya angkatan 2004 terima kasih atas kerja samanya selama ini. 12. Facebook dan Jamaah Facebookiyah yang selalu menemaniku dan menyemangatiku selama mengambil data dan menyelesaikan laporan ini. 13. Semua pihak yang tidak dapat saya sebutkan satu per satu yang telah membantu pelaksanaan dan penyusunan laporan Tugas Akhir ini. Penulis menyadari bahwa dalam penyusunan skripsi ini masih jauh dari sempurna, maka kritik dan saran yang membangun sangat penulis harapkan untuk kesempurnaan skripsi ini. Semoga skripsi ini dapat berguna bagi ilmu pengetahuan dan kita semua. Amin. Surakarta, Juli 2009 Penulis

ix

DAFTAR ISI LEMBAR PENGESAHAN ................................................................................ iii PERSEMBAHAN ............................................................................................... iv MOTTO .............................................................................................................

v

ABSTRAK .......................................................................................................... vi KATA PENGANTAR ........................................................................................ viii DAFTAR ISI.......................................................................................................

x

DAFTAR TABEL............................................................................................... xii DAFTAR GAMBAR .......................................................................................... xiii DAFTAR NOTASI ............................................................................................. xv DAFTAR LAMPIRAN....................................................................................... xvii BAB I

BAB II

PENDAHULUAN ..............................................................................

1

1.1 Latar Belakang ............................................................................

1

1.2 Perumusan Masalah ....................................................................

2

1.3 Batasan Masalah .........................................................................

2

1.4 Tujuan dan Manfaat ....................................................................

3

1.5 Sistematika Penulisan .................................................................

4

DASAR TEORI ..................................................................................

5

2.1 Tinjauan Pustaka .........................................................................

5

2.2 Dasar Teori..................................................................................

6

2.2.1 Sirip ...................................................................................

6

2.2.2 Sirip Pin............................................................................. 11 2.2.3 Macam-Macam Bentuk Sirip pin ...................................... 12 2.2.3.1 Silinder ................................................................ 12 2.2.3.2 Kubus .................................................................. 12 2.2.3.3 Oblong................................................................. 13 2.2.3.4 Ellips ................................................................... 13 2.2.4 Aplikasi Sirip pin .............................................................. 14 2.2.5 Perpindahan Panas ............................................................ 15 2.2.6 Parameter Tanpa Dimensi ................................................. 16 2.2.7 Perhitungan Perpindahan Panas dan Faktor Gesekan pada Pin-Fin Assembly...................................................... 17

x

2.2.7.1 Perhitungan Perpindahan Panas (Heat Transfer) 17 2.2.7.2 Perhitungan Faktor Gesekan (Friction Factor)... 22 2.2.7.3 Perhitungan Unjuk Kerja Termal pada Pin-Fin Assembly.............................................................. 22 BAB III METODOLOGI PENELITIAN.......................................................... 24 3.1 Tempat penelitian........................................................................ 24 3.2 Alat penelitian ............................................................................. 24 3.3 Spesimen ..................................................................................... 28 3.4 Pelaksanaan Penelitian ................................................................ 30 3.4.1 Tahap Persiapan ................................................................ 30 3.4.2 Tahap Pengujian................................................................ 30 3.5 Teknik Analisis Data................................................................... 31 3.6 Diagram Alir Penelitian .............................................................. 32 BAB IV DATA DAN ANALISIS..................................................................... 33 4.1 Data Hasil Pengujian................................................................... 33 4.2 Perhitungan Data......................................................................... 39 4.3 Analisis Data ............................................................................... 47 4.3.1 Pengaruh Bilangan Reynolds dan Jarak Antar Titik Pusat Sirip dalam Arah Streamwise terhadap Karakteristik Perpindahan Panas ............................................................ 47 4.3.2 Pengaruh Bilangan Reynolds dan Jarak Antar Titik Pusat Sirip dalam Arah Streamwise terhadap Karakteristik Penurunan Tekanan........................................................... 49 4.3.3 Pengaruh Bilangan Reynolds dan Jarak Antar Titik Pusat Sirip dalam Arah Streamwise terhadap Unjuk Kerja Termal ............................................................................... 51 BAB V PENUTUP........................................................................................... 53 5.1 Kesimpulan ................................................................................. 53 5.2 Saran............................................................................................ 53 DAFTAR PUSTAKA ......................................................................................... 54 LAMPIRAN........................................................................................................ 56

xi

DAFTAR TABEL Tabel 3.1 Spesifikasi spesimen penelitian ........................................................ 29 Tabel 4.1 Data hasil pengujian spesimen 1 (Sx/D = 2,95; Sy/D = 1,97) ............ 34 Tabel 4.2 Data hasil pengujian spesimen 2 (Sx/D = 2,95; Sy/D = 2,36) ............ 35 Tabel 4.3 Data hasil pengujian spesimen 3 (Sx/D = 2,95; Sy/D = 2,95) ............ 36 Tabel 4.4 Data hasil pengujian spesimen 4 (Sx/D = 2,95; Sy/D = 3,94) ............ 37 Tabel 4.5 Data hasil pengujian spesimen 5 (pelat tanpa sirip).......................... 38

xii

DAFTAR GAMBAR Gambar 2.1

Beberapa contoh jenis extended surface: (a) sirip longitudinal dengan profil segiempat (b) pipa silindris dengan sirip berprofil segiempat (c) sirip longitudinal dengan profil trapesium (d) sirip longitudinal dengan profil parabola (e) pipa silindris dengan sirip radial berprofil segiempat (f) pipa silindris

dengan

sirip

radial

berprofil

trapesium (g)

cylindrical spine (h) truncated conical spine (i) truncated parabolic spine............................................................................ Gambar 2.2

7

Beberapa contoh jenis permukaan penukar panas kompak: (a) pipa silinder (b) pipa silinder dengan sirip-sirip radial (c) flat tube dengan sirip-sirip kontinyu (d) sirip plat (plate fin) (e) offset plate fin (f) crossed rod matrix..........................................

8

Gambar 2.3

Perbedaan-perbedaan gradien temperatur dalam sirip ................

9

Gambar 2.4

Sebuah susunan sirip pin............................................................. 11

Gambar 2.5

Susunan sirip pin (a) inline (b) staggered ................................... 12

Gambar 2.6

Perbandingan sirip pin silinder lurus dengan sirip pin silinder berfillet ........................................................................................ 12

Gambar 2.7

Perbandingan antara konfigurasi susunan staggered sirip pin kubus dan sirip pin diamond ....................................................... 13

Gambar 2.8

Konfigurasi susunan staggered menggunakan sirip pin oblong . 13

Gambar 2.9

Ukuran relatif dari circular fin, SEF dan N fin. .......................... 13

Gambar 2.10 Potongan melintang sudu turbin dengan pendinginan dalam (internal cooling) ........................................................................ 14 Gambar 2.11 Susunan sirip pin dalam suatu saluran udara segiempat dengan clearance nol .................................................................. 19 Gambar 3.1

Skema alat penelitian .................................................................. 24

Gambar 3.2

Satu set alat uji inline cylindrical pin fin assembly..................... 24

Gambar 3.3

Pelurus aliran udara (flow straightener)...................................... 25

Gambar 3.4

Fan hisap..................................................................................... 25

Gambar 3.5

Rheostat....................................................................................... 25

Gambar 3.6

Anemometer................................................................................ 26

xiii

Gambar 3.7

Pemanas listrik (electric heater) ................................................. 26

Gambar 3.8

Regulator pengatur tegangan listrik yang masuk heater............. 26

Gambar 3.9

Multitester digital ........................................................................ 27

Gambar 3.10 Amperemeter............................................................................... 27 Gambar 3.11 Manometer U dan posisi titik pengukuran tekanan..................... 27 Gambar 3.12 Termokopel tipe T....................................................................... 28 Gambar 3.13 Thermocouple reader .................................................................. 28 Gambar 3.14 Dimensi dan tata nama spesimen ................................................ 28 Gambar 3.15 Model spesimen .......................................................................... 29 Gambar 3.16 Spesimen 1 .................................................................................. 29 Gambar 3.18 Spesimen 2 .................................................................................. 29 Gambar 3.18 Spesimen 3 .................................................................................. 30 Gambar 3.19 Spesimen 4 ................................................................................. 30 Gambar 4.1

Posisi titik pengukuran (a) temperatur udara masuk saluran (b) temperatur base plate dan (c) temperatur udara keluar saluran. ........................................................................................ 33

Gambar 4.2

Pengaruh bilangan Reynolds terhadap koefisien perpindahan panas konveksi rata-rata pada Sx/D = 2,95.................................. 47

Gambar 4.3

Pengaruh bilangan Reynolds terhadap bilangan Nusselt pada Sx/D = 2,95 .................................................................................. 48

Gambar 4.4

Pengaruh nilai Sy/D terhadap bilangan Nusselt pada Sx/D = 2,95.............................................................................................. 49

Gambar 4.5

Pengaruh bilangan Reynolds terhadap penurunan tekanan pada Sx/D = 2,95.......................................................................... 50

Gambar 4.6

Pengaruh bilangan Reynolds terhadap faktor gesekan pada Sx/D = 2,95 .................................................................................. 50

Gambar 4.7

Pengaruh bilangan Reynolds terhadap unjuk kerja termal pada Sx/D = 2,95 .................................................................................. 51

xiv

DAFTAR NOTASI A

= luas penampang saluran udara (m2)

Afront

= luas frontal dari sirip-sirip (m2)

As

= luas seluruh permukaan yang kontak dengan udara dari susunan sirip pin (m2)

Cp

= panas jenis udara (J/kg.K)

D, d

= diameter sirip pin (m)

Dh

= diameter hidrolik dari saluran udara (m)

f

= faktor gesekan

fs

= faktor gesekan untuk spesimen tanpa sirip

fa

= faktor gesekan untuk spesimen dengan sirip

h

= koefisien perpindahan panas konveksi rata-rata (W/m2.K)

ha

= koefisien perpindahan panas konveksi spesimen dengan sirip (W/m2.K)

hs

= koefisien perpindahan panas konveksi spesimen tanpa sirip (W/m2.K)

H

= tinggi saluran udara atau sirip pin (m)

k

= konduktivitas termal udara (W/m.K)

L

= panjang plat dasar untuk susunan sirip pin (m)

Lt

= panjang seksi uji (m)

m&

= laju aliran massa udara (kg/s)

Nu

= duct Nusselt number

NuD

= pin Nusselt number

Nf

= jumlah total sirip pin dalam susunan sirip pin

∆P

= perbedaan tekanan statik (N/m2)

∆ Ps

= penurunan tekanan tanpa halangan (N/m2)

∆ Pa

= penurunan tekanan dengan halangan (N/m2)

Q

= laju perpindahan panas (Watt)

Qelect

= laju aliran panas dari listrik (W)

Qconv

= laju perpindahan panas konveksi (W)

Qloss

= laju aliran panas yang hilang dari sistem (heat loss) (W)

Qrad

= laju perpindahan panas radiasi (W)

Qcond

= laju perpindahan panas konduksi (W)

ReD

= pin Reynolds number

xv

Re

= duct Reynolds number

Res

= bilangan Reynolds untuk spesimen tanpa sirip

Rea

= bilangan Reynolds untuk spesimen dengan sirip

Sy

= jarak antar titik pusat sirip sepanjang arah aliran (streamwise direction) (mm)

Sx

= jarak antar titik pusat sirip yang diukur tegak lurus terhadap arah aliran (spanwise direction) (mm)

Tin

= temperatur inlet dari aliran udara (K)

Tout

= temperatur outlet dari aliran udara (K)

Tb

= temperatur base plate (K)

T

= temperatur (K)

Tp

= temperatur intermediate (K)

Ts

= temperatur lingkungan (K)

Tf

= temperatur udara rata-rata (K)

V

= kecepatan rata-rata udara dalam saluran udara (m/s)

Vmax

= kecepatan udara maksimum yang melalui sirip pin (m/s)

V&a

= laju aliran volumetrik di atas plat dengan halangan (blocks) (m3/s)

V&s

= laju aliran volumetrik di atas plat tanpa halangan (blocks) (m3/s)

Wb

= lebar plat dasar untuk susunan sirip pin (m)

x

= ketebalan bahan (m)

γ

= sudut orientasi pada susunan sirip pin oblong (o)

ε

= emisivitas

η

= unjuk kerja termal (%)

µ

= viskositas dinamik udara (kg/m.s)

ν

= viskositas kinematik udara (m2/s)

ρ

= massa jenis udara (kg/m3)

σ

= konstanta Stefan-Boltzmann = 5,669 x 10-8 (W/m2.K4)

xvi

DAFTAR LAMPIRAN Lampiran A. Data hasil pengujian Lampiran B. Properti berbagai jenis material Lampiran C. Hasil perhitungan data pengujian

xvii

BAB I 1

PENDAHULUAN I

1.1 Latar Belakang Perluasan permukaan perpindahan panas menggunakan sirip-sirip (fins) sering digunakan dalam peralatan penukar panas yang bertujuan untuk meningkatkan perpindahan panas antara permukaan utama dan fluida di sekitarnya. Sirip-sirip tersebut menonjol keluar dari sebuah permukaan dasar segiempat atau silindris. Berbagai tipe sirip alat penukar panas, mulai dari bentuk yang sederhana, seperti sirip segiempat (rectangular), silindris, annular, tirus (tapered) atau pin, sampai kombinasi dari berbagai geometri yang berbeda telah digunakan. Tipe sirip yang digunakan tergantung dari proses permesinan dan ruang yang tersedia dalam peralatan pembangkit panas yang terlibat dalam proses pendinginan. Salah satu tipe sirip alat penukar panas yang biasa digunakan adalah sirip pin. Sirip pin adalah elemen berbentuk silinder atau bentuk lainnya yang dipasang secara tegak lurus terhadap dinding alat penukar panas, dengan fluida pendingin mengalir dalam arah aliran melintang (crossflow) terhadap elemen tersebut. Siripsirip pin yang menonjol dari sebuah permukaan yang dipanaskan dapat meningkatkan luasan permukaan disipasi panas, sehingga meningkatkan unjuk kerja disipasi panas yang berdampak pada meningkatnya ketahanan (reliability) dan umur peralatan. Terdapat berbagai parameter yang menggolongkan sirip pin, seperti bentuk, tinggi, diameter, perbandingan tinggi-diameter (H/D) dan sebagainya. Sirip-sirip pin dapat disusun secara segaris (inline) maupun selangseling (staggered) terhadap arah aliran. Karakteristik perpindahan panas dan penurunan tekanan dalam sistem sirip pin telah menjadi bahan penelitian secara luas karena perannya yang penting dalam berbagai aplikasi keteknikan. Sirip pin dengan perbandingan tinggidiameter (H/D) antara 0,5 – 4 dikategorikan sebagai sirip pin pendek (short pin

fin), sedangkan sirip pin panjang (long pin fin) memiliki perbandingan tinggi diameter > 4. Perbandingan tinggi-diameter yang besar merupakan bagian yang menarik dalam aplikasi alat penukar panas dalam hal pencapaian koefisien perpindahan panas yang sangat tinggi. Sirip pin banyak digunakan dalam berbagai 1

2 aplikasi industri, khususnya dalam alat penukar panas kompak (compact heat

exchanger), trailing edge sudu turbin gas, dan beberapa sistem elektronik modern. Laju perpindahan panas dari suatu rakitan sirip pin (pin fin assembly) ke lingkungan tergantung pada distribusi temperatur pada sirip pin dan plat dasar (base plate), geometri sirip pin, shroud clearance (jarak antara ujung sirip pin dengan permukaan atas saluran udara), sifat-sifat fluida dan sirip pin, laju aliran udara, jarak antar titik pusat sirip (interfin pitch), susunan sirip pin, dan orientasi dari penukar panas (terutama untuk laju aliran udara yang rendah). Untuk plat dasar dengan temperatur tertentu, laju perpindahan panas dapat ditingkatkan dengan menaikkan koefisien perpindahan panas rata-rata, menaikkan luas permukaan perpindahan panas atau kedua-duanya. Kenaikan perpindahan panas dapat dicapai dengan cara konveksi paksa (forced convection) atau mengubah konfigurasi geometri dari alat penukar panas. Dalam praktiknya, cara-cara ini dibatasi oleh penurunan tekanan maksimum yang diijinkan melalui susunan sirip pin tersebut karena kenaikan perpindahan panas akan disertai penurunan tekanan. Energi yang hilang karena penurunan tekanan dapat melebihi energi yang didapatkan dari usaha peningkatan perpindahan panas tersebut. Penelitian ini bertujuan untuk menguji karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal dari cylindrical pin fin assembly yang disusun secara segaris (inline) dalam saluran segiempat (rectangular channel). 1.2 Perumusan Masalah Bagaimanakah pengaruh bilangan Reynolds dan jarak antar titik pusat sirip dalam arah aliran udara (streamwise direction) terhadap karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal dari cylindrical pin fin

assembly yang disusun secara segaris (inline) dalam saluran segiempat (rectangular channel). 1.3 Batasan Masalah Pada penelitian ini masalah dibatasi sebagai berikut: 1. Material sirip pin dan base plate yang digunakan adalah duralumin. 2. Dimensi base plate yang digunakan adalah; panjang 200 mm, lebar 150 mm dan tebal 6,5 mm.

3 3. Dimensi pin fin yang digunakan adalah; tinggi 75 mm dan diameter 12,7 mm, atau H/D = 5,9. 4. Jarak antara ujung sirip pin dengan permukaan atas saluran udara (shroud clearance) adalah nol. 5. Penelitian menggunakan alat saluran udara segiempat yang terdiri dari: 9 Saluran udara segiempat berdimensi penampang 150 mm x 75 mm. 9 Fan hisap. 9 Pemanas (heater). 9 Pelurus aliran udara (flow straightener). 9 Manometer U. 6. Permukaan dalam saluran udara dilapisi dengan bahan melamin yang halus, sehingga faktor gesekan diabaikan. 7. Permukaan luar saluran udara dimana seksi uji diletakkan diisolasi dengan glasswool dan styrofoam sehingga perpindahan panas ke lingkungan diabaikan. 8. Parameter yang dibuat konstan yaitu temperatur permukaan base plate sebesar 60 oC, temperatur udara masuk, jarak antar titik pusat sirip dalam arah melintang (spanwise direction) sebesar 37,5 mm. 9. Parameter yang divariasi adalah kecepatan aliran udara masuk yaitu sebesar 0,5 m/s, 1 m/s, 2 m/s, 3 m/s, 4 m/s, 5 m/s, 5,5 m/s dan 6 m/s serta jarak antar titik pusat sirip dalam arah aliran udara (streamwise

direction) yaitu sebesar 25 mm, 30 mm, 37,5 mm dan 50 mm. 10. Pengujian karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal dari cylindrical pin fin assembly dilakukan pada kondisi tunak (steady state). 11. Penelitian dilakukan dalam keadaan diam (static experiment) dan pada temperatur kamar. 1.4 Tujuan dan Manfaat Penelitian ini bertujuan untuk: 1. Mengetahui

pengaruh

bilangan

Reynolds

terhadap

karakteristik

perpindahan panas dan penurunan tekanan dari cylindrical pin fin

assembly dengan susunan sirip inline dalam rectangular channel.

4 2. Mengetahui pengaruh jarak antar titik pusat sirip dalam arah aliran udara (streamwise direction) terhadap karakteristik perpindahan panas dan penurunan tekanan dari cylindrical pin fin assembly dengan susunan sirip inline dalam rectangular channel. 3. Mengetahui pengaruh variasi bilangan Reynolds dan jarak antar titik pusat sirip dalam arah aliran udara terhadap unjuk kerja termal dari

cylindrical pin fin assembly dengan susunan sirip inline dalam rectangular channel. Hasil penelitian yang diperoleh diharapkan dapat memberikan manfaat sebagai berikut: 1. Mampu memberikan pengetahuan baru yang berguna dalam ilmu perpindahan panas, khususnya mengenai karakteristik perpindahan panas dan unjuk kerja termal dari cylindrical pin fin assembly dengan susunan sirip inline dalam rectangular channel. 2. Dapat diterapkan dalam sistem pendinginan sudu-sudu turbin gas dan sistem elektronik modern. 1.5 Sistematika Penulisan Sistematika penulisan Tugas Akhir ini adalah sebagai berikut: BAB I

: Pendahuluan, menjelaskan tentang latar belakang masalah, perumusan masalah, batasan masalah, tujuan dan manfaat penelitian.

BAB II

: Dasar teori, berisi tinjauan pustaka yang berkaitan dengan pengujian susunan sirip pin, dasar teori tentang sirip pin dan teori perhitungan perpindahan panas, penurunan tekanan dan unjuk kerja termal dari susunan sirip pin dalam saluran segiempat.

BAB III : Metodologi

penelitian,

menjelaskan

tempat

dan

pelaksanaan

penelitian, peralatan yang digunakan, langkah-langkah penelitian dan pengambilan data. BAB IV : Data dan analisis, menjelaskan data hasil pengujian, perhitungan data hasil pengujian serta analisis hasil dari perhitungan. BAB V

: Penutup, berisi tentang kesimpulan dan saran.

BAB II 2

DASAR TEORI

I

2.1 Tinjauan Pustaka

Solid block dapat meningkatkan perpindahan panas dari sebuah plat secara signifikan karena meningkatnya luasan permukaan perpindahan panas, tetapi menimbulkan energi yang terbuang karena adanya penurunan tekanan (pressure

drop) yang lebih besar dalam aliran dan lebih sedikitnya aliran udara yang kontak dengan plat. Selain itu, meningkatnya bilangan Reynolds menyebabkan menurunnya unjuk kerja termal (Sara,1999). Tahat, M. et al. (2000) melakukan penelitian tentang perpindahan panas kondisi tunak pada alat penukar panas tipe plat bersirip pin yang tersusun secara

inline dan staggered untuk menentukan desain optimum alat penukar panas tersebut. Dalam penelitian tersebut menggunakan spesimen berupa plat datar berdimensi 250 mm x 300 mm yang diberi sirip pin berbentuk silinder pejal dengan diameter 8 mm dan panjang 90 mm dengan jarak antar titik pusat sirip,

Sx/D = 9,86 – 63,44 dan Sy/D = 1,09 – 83,92. Spesimen diletakkan pada saluran udara segiempat (rectangular channel) yang berdimensi 405 mm x 100 mm x 3.000 mm dengan variasi kecepatan aliran udara 6 m/s, 7 m/s dan 7,8 m/s. Dari penelitian tersebut diperoleh bahwa laju panas yang hilang meningkat seiring dengan meningkatnya bilangan Reynolds, namun menurun seiring dengan meningkatnya jarak antar titik pusat sirip pin untuk arah streamwise dan spanwise. Bilen, K. et al. (2001) melakukan penelitian tentang karakteristik perpindahan panas, faktor gesekan (friction factor) dan unjuk kerja termal pada suatu permukaan bersirip di dalam saluran udara segiempat. Pada penelitian tersebut menggunakan spesimen berupa plat datar berdimensi 180 mm x 300 mm x 2 mm yang diberi sirip pin berbentuk silinder pejal dengan diameter 17 mm dan panjang 100 mm yang disusun secara inline dan staggered dengan jarak antar titik pusat sirip, Sx/D = 2,2 dan Sy/D = 1,96 – 4,41. Spesimen diletakkan dalam saluran udara segiempat (rectangular channel) yang berdimensi 180 mm x 100 mm x 2.000 mm. Pengujian dilakukan dengan variasi bilangan Reynolds sebesar 3.700 – 30.000. Dari penelitian tersebut diperoleh hasil bahwa dengan meningkatnya bilangan

Reynolds

akan

meningkatkan 5

bilangan

Nusselt,

yang

berarti

6 meningkatkan perpindahan panas. Tetapi dengan meningkatnya bilangan Reynolds, peningkatan perpindahan panas (heat transfer enhancement) dan unjuk kerja termal semakin menurun, dimana perpindahan panas maksimum terjadi pada variasi Sy/D = 2,94. Penambahan sirip pin pada suatu permukaan meningkatkan perpindahan panas dari permukaan tersebut sebagai hasil dari kenaikan luasan permukaan perpindahan panas, tetapi disertai dengan adanya penurunan tekanan yang lebih besar dalam saluran. Sahin, B. dan Demir, A. (2008) melakukan penelitian tentang peningkatan perpindahan panas dan penurunan tekanan pada plat bersirip pin di dalam suatu saluran udara segiempat. Penelitian dilakukan dengan menggunakan spesimen berupa plat datar berdimensi 250 mm x 250 mm x 6 mm yang diberi sirip pin berbentuk silinder pejal dengan diameter 15 mm dan panjang 100 mm yang disusun secara inline dan staggered dengan jarak antar titik pusat sirip, Sx/D = 2,2 dan Sy/D = 1,208 – 3,417. Spesimen diletakkan pada saluran udara segiempat yang berdimensi 250 mm x 100 mm x 3.140 mm dengan clearance ratio sebesar 0, 0,33 dan 1. Pengujian dilakukan dengan variasi bilangan Reynolds sebesar 13.500 – 42.000. Dari penelitian tersebut diperoleh bahwa penggunaan sirip pin berpenampang lingkaran dapat meningkatkan perpindahan panas. Efisiensi meningkat seiring dengan menurunnya bilangan Reynolds. Faktor gesekan dan bilangan Nusselt meningkat seiring dengan meningkatnya jarak antar titik pusat sirip pin dan clearance ratio. Parameter terpenting yang mempengaruhi perpindahan panas adalah bilangan Reynolds, tinggi sirip pin dan jarak antar titik pusat sirip pin, dimana perpindahan panas dapat ditingkatkan dengan mengontrol parameter-parameter tersebut. 2.2 Dasar Teori 2.2.1 Sirip Permukaan perpindahan panas yang diperluas (extended surface heat

transfer) adalah studi tentang komponen-komponen perpindahan panas berunjuk kerja tinggi yang berkenaan dengan berat, volume, dan biaya yang semakin kecil, dan perilakunya dalam berbagai kondisi lingkungan termal. Komponen-komponen tertentu telah diterapkan dalam berbagai macam aplikasi seperti dalam prosesproses kimia, refrigerasi, dan kriogenika, dalam peralatan listrik dan elektronika,

7 dalam tungku konvensional dan turbin gas, dalam proses pembuangan panas pada boiler, dan dalam modul bahan bakar nuklir. Dalam desain dan konstruksi dari berbagai macam peralatan perpindahan panas, bentuk-bentuk sederhana seperti silinder, batang dan plat biasa diterapkan pada aliran panas antara sumber panas dan penyerap panas (heat source and heat

sink). Permukaan-permukaan penyerap panas maupun pembuang panas masingmasing dikenal sebagai permukaan utama (prime surface). Apabila permukaan utama diperluas dengan permukaan tambahan seperti dalam gambar 2.1, maka gabungan antara kedua permukaan tersebut dinamakan permukaan yang diperluas (extended surface). Elemen yang digunakan untuk memperluas permukaan utama dikenal sebagai sirip (fin). Jika elemen sirip tersebut berbentuk kerucut atau silinder, sirip tersebut dinamakan spines atau pegs.

(a)

(e)

(b)

(f)

(c)

(g)

(d)

(h)

(i)

Gambar 2.1 Beberapa contoh jenis extended surface: (a) sirip longitudinal dengan profil segiempat (b) pipa silindris dengan sirip berprofil segiempat (c) sirip longitudinal dengan profil trapesium (d) sirip longitudinal dengan profil parabola (e) pipa silindris dengan sirip radial berprofil segiempat (f) pipa silindris dengan sirip radial berprofil trapesium (g) cylindrical spine (h) truncated conical spine (i) truncated parabolic spine

Kebutuhan untuk perlengkapan turbin gas, pengkondisian udara, dan kriogenika telah mendapatkan perhatian khusus dalam hal keringkasan permukaan alat penukar panas, terutama pada permukaan yang mengalami gradien tekanan yang kecil dalam fluida yang bersirkulasi melaluinya. Beberapa diantaranya ditunjukkan dalam gambar 2.2. Keringkasan (compactness) mengacu pada perbandingan luas permukaan perpindahan panas per satuan volume alat penukar panas.

8

(a)

(d)

(b)

(e)

(c)

(f)

Gambar 2.2 Beberapa contoh jenis permukaan penukar panas kompak: (a) pipa silinder (b) pipa silinder dengan sirip-sirip radial (c) flat tube dengan sirip-sirip kontinyu (d) sirip plat (plate fin) (e) offset plate fin (f) crossed rod matrix

Definisi awal dari Kays dan London (1950) menetapkan bahwa elemen alat penukar panas kompak adalah alat penukar panas yang mempunyai kelebihan 245 m2 per meter kubik alat penukar panas. Elemen alat penukar panas kompak telah tersedia lebih dari 4.100 m2 per meter kubik dibandingkan dengan 65 – 130 m2 per meter kubik untuk alat penukar panas konvensional dengan pipa 5/8 – 1 in. Kebanyakan elemen alat penukar panas kompak terdiri dari plat-plat permukaan utama atau pipa-pipa yang dipisahkan oleh plat, batang atau spines, yang juga bertindak sebagai sirip. Seperti pada gambar 2.2(d), setiap sirip dapat diperlakukan sebagai sirip tunggal dengan tinggi sirip sama dengan setengah dari jarak plat pemisah dan dengan plat pemisah bertindak sebagai permukaan utama. Sehingga, alat penukar panas kompak dipandang sebagai bentuk lain dari permukaan yang diperluas (extended surface). Apabila sirip dan permukaan utamanya ditempatkan pada lingkungan termal yang seragam, efektivitas permukaan sirip lebih kecil daripada permukaan utama. Hal ini dapat dilihat pada plat dengan sirip memanjang (longitudinal) pada penampang melintang segiempat pada gambar 2.3.

9 Hot face of plate

TP T Source (T1)

Fin

TP

T

Ts (Surroundings)

Gambar 2.3 Perbedaan-perbedaan gradien temperatur dalam sirip

Permukaan plat bagian dalam membuang panas dari sumber panas dengan koefisien perpindahan panas seragam dan temperatur T1, sedangkan permukaan plat bagian luar dan sirip membuang panas tersebut ke lingkungan sekelilingnya yang lebih dingin dengan koefisien perpindahan panas seragam dan temperatur Ts. Permukaan plat yang lebih dingin berada pada temperatur intermediate Tp, dan panas dari sumber panas meninggalkan plat karena adanya gradien temperatur, Tp

– Ts. Dengan cara yang sama, permukaan sirip dengan temperatur tertentu T, dan panas meninggalkan sirip karena adanya gradien temperatur T – Ts. Panas memasuki sirip melalui dasarnya (base), dan bergerak berpindah secara kontinyu melalui sirip secara konduksi. Hampir dalam setiap kasus, temperatur dasar sirip akan sangat mendekati sama dengan Tp. Panas yang diserap oleh sirip melalui dasarnya dapat mengalir menuju ujung sirip hanya jika ada gradien temperatur dalam sirip, yakni Tp lebih besar dari T. Untuk kondisi ini, karena temperatur T bervariasi dari dasar hingga ujung sirip, gradien temperatur T – Ts akan lebih kecil daripada Tp – Ts, sehingga satu satuan luas permukaan sirip akan kurang efektif dibandingkan dengan satu satuan luas permukaan plat atau permukaan utama. Efisiensi sirip didefinisikan sebagai perbandingan antara panas aktual yang hilang dari sirip terhadap panas ideal yang hilang jika temperatur seluruh permukaan sirip sama dengan temperatur dasarnya. Sirip dengan ukuran, bentuk, dan material tertentu memiliki efisiensi sirip yang berbeda-beda, dan efisiensi sirip akan bervariasi terhadap konduktivitas termal dan cara perpindahan panasnya berkenaan dengan lingkungannya. Saat ini telah terdapat beberapa referensi yang dibuat mengenai permukaan yang diperluas yang berisikan beberapa tipe

10 permukaan utama dan beberapa tipe sirip. Banyak pengetahuan mengenai aliran panas, profil temperatur, efisiensi dan optimasi parameter-parameter sirip dapat diperoleh dari analisis tiga geometri dasar sirip yang ditunjukkan dalam gambar 2.1 yaitu sirip longitudinal, sirip radial dan spines. Sirip dengan berbagai macam geometri dan konduktivitas termal akan memberikan tanggapan yang berbeda terhadap sumber panas dan penyerap panas (source and sink) yang seragam dan identik. Sebaliknya, terdapat banyak hal mengapa temperatur dan koefisien perpindahan panas dari sumber panas dan penyerap panas bisa bervariasi. Hal penting dalam menganalisis geometri sirip adalah asumsi-asumsi yang diambil untuk menentukan dan membatasi masalah dan tentunya menyederhanakan penyelesaiannya. Analisis dari tiga geometri dasar sirip dapat dilakukan dengan mengacu pada asumsi Murray-Gardner, yaitu: 1. Aliran panas dalam sirip dan temperaturnya tetap konstan terhadap waktu. 2. Material sirip adalah homogen, konduktivitas termal sama di segala arah, dan tetap konstan. 3. Koefisien perpindahan panas konveksi di permukaan sirip konstan dan seragam di keseluruhan permukaan sirip. 4. Temperatur dari medium lingkungan sirip konstan. 5. Ketebalan sirip adalah kecil, dibandingkan dengan tinggi dan panjangnya, sehingga gradien temperatur melewati ketebalan sirip dan perpindahan panas dari tepi sirip dapat diabaikan. 6. Temperatur dasar sirip adalah seragam. 7. Tidak ada tahanan kontak dimana dasar sirip digabung dengan permukaan utama. 8. Tidak ada sumber panas di dalam sirip itu sendiri. 9. Panas yang ditransfer melalui ujung sirip diabaikan dibandingkan dengan panas yang meninggalkan permukaan lateralnya. 10. Perpindahan panas ke atau dari sirip sebanding terhadap kelebihan temperatur antara sirip dan medium sekitar.

11 2.2.2 Sirip Pin Sirip pin adalah elemen berbentuk silinder atau bentuk lainnya yang dipasang secara tegak lurus terhadap dinding alat penukar panas, dengan fluida pendingin mengalir dalam arah aliran melintang (crossflow) terhadap elemen tersebut. Terdapat berbagai parameter yang menggolongkan sirip pin, seperti bentuk, tinggi, diameter, perbandingan tinggi-diameter (H/D) dan sebagainya. Sirip pin dengan perbandingan tinggi-diameter (H/D) antara 0,5 – 4 dikategorikan sebagai sirip pin pendek (short pin fin), sedangkan sirip pin panjang (long pin fin) memiliki perbandingan tinggi-diameter > 4. Perbandingan tinggi-diameter yang besar merupakan bagian yang menarik dalam aplikasi alat penukar panas dalam hal pencapaian koefisien perpindahan panas yang sangat tinggi.

Gambar 2.4 Sebuah susunan sirip pin

Sirip pin dapat disusun dalam dua arah utama. Pada gambar 2.5 sirip-sirip pin ditunjukkan dalam susunan segaris (inline) dan selang-seling (staggered). Sx adalah jarak antar titik pusat sirip yang diukur normal/tegak lurus terhadap arah aliran (spanwise direction), sedangkan Sy adalah jarak antar titik pusat sirip sepanjang arah aliran (streamwise direction).

12

Sy Sy

Sx Sx

(a)

(b)

Gambar 2.5 Susunan sirip pin (a) inline (b) staggered

2.2.3 Macam-Macam Bentuk Sirip Pin 2.2.3.1 Silinder Silinder lurus merupakan geometri sirip pin yang paling umum. Geometri sirip pin yang lain adalah silinder lurus dengan ujung difillet (filleted pin fin) dan silinder tirus (tapered cylindrical pin fin). Hubungan antara geometri sirip pin silinder lurus dengan geometri sirip pin silinder berfillet ditunjukkan dalam gambar 2.6.

Gambar 2.6 Perbandingan sirip pin silinder lurus dengan sirip pin silinder berfillet

2.2.3.2 Kubus Sirip pin berbentuk kubus (cubic pin fin) dapat disusun secara segiempat maupun diamond berdasarkan arah aliran. Diamond merupakan susunan segiempat yang diputar 45o. Gambar 2.7 menunjukkan sketsa kedua tipe susunan.

13

Gambar 2.7 Perbandingan antara konfigurasi susunan staggered sirip pin kubus dan sirip pin diamond

2.2.3.3 Oblong Sirip pin oblong merupakan perpaduan antara bentuk silinder dan bentuk kubus. Sirip pin oblong tersusun pada sudut orientasi yang berbeda, γ, berdasarkan arah aliran. Gambar 2.8 menunjukkan tata nama yang digunakan dalam sirip pin oblong.

Gambar 2.8 Konfigurasi susunan staggered menggunakan sirip pin oblong

2.2.3.4 Ellips Sirip pin ellips adalah bentuk silinder yang direntangkan dalam satu arah garis diameternya. Gambar 2.9 menunjukkan sketsa geometri circular fin dan dua bentuk sirip pin ellips. Dalam menguji sirip pin ellips, sumbu utama (major axis) segaris dengan arah aliran.

Gambar 2.9 Ukuran relatif dari circular fin, SEF dan N fin.

14 Bentuk-bentuk sirip pin ellips dapat dijelaskan sebagai berikut: a. Standard Elliptical Fin (SEF) Sirip pin ini mempunyai standar penampang ellips dengan panjang sumbu minor sama dengan diameter circular fin dan panjang sumbu utama adalah 1,67 kali panjang sumbu minor. Luas permukaan dari SEF dihitung 1,35 kali luas permukaan circular fin, tetapi luas frontal efektif sama dengan circular fin karena panjang sumbu minor sama dengan diameter circular fin. b. N fin Bentuk sirip ini diturunkan dari seri airfoil simetris 4 digit NACA. Panjang sumbu minor sama dengan diameter circular fin dan panjang sumbu utama 2,5 kali panjang sumbu minor. Luas permukaan N fin dihitung 1,85 kali lebih besar daripada luas permukaan circular fin. Luas frontal efektif adalah sama dengan circular fin. 2.2.4 Aplikasi Sirip Pin Perpindahan panas dari susunan sirip pin merupakan subjek yang sangat penting dengan banyak aplikasi keteknikan. Aplikasi tersebut mulai dari alat penukar panas kompak, boiler untuk turbin uap dan pendinginan internal secara konveksi dari turbin gas. Seperti pada gambar 2.10, sirip pin biasanya dimasukkan dalam ruang pendinginan dalam (internal cooling) dekat trailing edge dari sudu turbin untuk meningkatkan perpindahan panas. Hal ini memungkinkan sudu beroperasi dalam temperatur tinggi tanpa mengalami kerusakan, sehingga meningkatkan efisiensi termal dan daya output.

Gambar 2.10 Potongan melintang sudu turbin dengan pendinginan dalam (internal cooling)

15 Oleh karena pertimbangan aerodinamis, seperti pada gambar 2.10, trailing

edge dari sudu menuntut profil yang semakin mengecil. Untuk itu, ruang pendinginan dalam daerah ini harus dengan bentuk penampang trapesium. Pendingin dari pangkal sudu (blade base) bergerak memutar ke samping kemudian dikeluarkan dari slot ujung sudu, atau melalui saluran sirip pin kemudian keluar dari slot sepanjang trailing edge sudu. Namun, kebanyakan penelitian yang dilakukan adalah untuk sirip pin yang menggunakan saluran segiempat (rectangular channel) dengan aliran keluar yang lurus. 2.2.5 Perpindahan Panas Perpindahan panas adalah ilmu yang mempelajari tentang cara untuk meramalkan perpindahan (distribusi) energi berupa panas yang terjadi karena adanya perbedaan temperatur di antara benda atau material. Perpindahan panas dapat dibagi menjadi tiga macam, yaitu perpindahan panas secara konduksi, konveksi dan radiasi. Perpindahan panas secara konduksi adalah distribusi energi berupa panas yang terjadi pada benda atau medium yang diam (padat) bertemperatur tinggi ke bagian benda yang bertemperatur rendah atau terdapat gradien temperatur pada benda tersebut. Untuk

menghitung perpindahan panas konduksi dapat

dipergunakan rumus: Q=

k A ∆T x

dimana: Q

(2.1) = laju perpindahan panas (Watt)

k

= konduktivitas termal (W/m.oC)

A

= luasan perpindahan panas arah normal Q (m2)

∆T

= beda temperatur (oC)

x

= ketebalan bahan (m)

Perpindahan panas konveksi adalah distribusi energi berupa panas yang terjadi karena terdapat aliran fluida. Perpindahan panas konveksi dapat dihitung dengan rumus sebagai berikut: Q = h A ∆T

dimana: h

(2.2) = koefisien perpindahan panas konveksi (W/m2 oC)

16 Perpindahan panas konveksi ada dua macam, yaitu: 1. Konveksi alami (natural convection) Adalah perpindahan panas konveksi yang terjadi karena berubahnya densitas fluida tersebut yang disebabkan adanya pemanasan. 2. Konveksi paksa (forced convection) Adalah perpindahan panas konveksi konveksi yang berlangsung dengan bantuan peralatan mekanis, misalnya udara yang dihembuskan di atas plat oleh kipas. Perpindahan panas radiasi adalah distribusi energi berupa panas yang terjadi melalui pancaran gelombang cahaya dari suatu zat ke zat yang lain tanpa zat perantara. Untuk menghitung besarnya panas yang dipancarkan digunakan rumus:

Q =εσ AT4 dimana: Q

(2.3) = panas yang dipancarkan (Watt)

ε

= emisivitas (0 s.d. 1)

A

= luas perpindahan panas (m2)

T

= temperatur (K)

σ

= konstanta Stefan-Boltzmann = 5,669 x 10-8 W/m2.K4

Khusus untuk benda hitam sempurna menurut hukum Stefan-Boltzmann:

Q =σ AT4

(2.4)

2.2.6 Parameter tanpa Dimensi Persamaan perpindahan panas konveksi berkaitan dengan variabel penting yang dinamakan parameter tanpa dimensi (dimensionless). Parameter tanpa dimensi dalam kaitannya dengan perpindahan panas konveksi adalah: a. Bilangan Reynolds (Reynolds Number) Bilangan Reynolds didefinisikan sebagai perbandingan antara gaya inersia dengan gaya kekentalan, di dalam lapis batas kecepatan. Untuk kontrol volume yang berbeda pada lapis batas ini, gaya inersia diasosiasikan dengan sebuah kenaikan momentum dari fluida yang bergerak melewati kontrol volume.

17 Gaya inersia dalam bentuk ∂[( ρ u) u ] / ∂x dapat didekati dengan persamaan: FI = ρV 2 L . Gaya kekentalan diwakili dengan gaya geser dalam bentuk ∂ τ yx ∂ y = ∂ [ µ (∂u ∂y )] ∂y , dapat didekati dengan persamaan: Fs = µV L2 .

Perbandingan kedua gaya tersebut dapat ditulis:

FI ρV 2 L ρVL = = = ReL Fs µV L2 µ

(2.5)

b. Bilangan Nusselt (Nusselt Number) Bilangan Nusselt adalah bilangan tanpa dimensi yang menyatakan perbandingan antara koefisien perpindahan panas konveksi terhadap konduktivitas termal fluida. Bilangan ini menyediakan sebuah perhitungan tentang perpindahan panas konveksi yang terjadi pada permukaan. Bilangan Nusselt dirumuskan:

Nu =

hd k

(2.6)

2.2.7 Perhitungan Perpindahan Panas dan Faktor Gesekan pada Pin Fin Assembly 2.2.7.1 Perhitungan Perpindahan Panas (Heat Transfer) Kesetimbangan energi kondisi tunak (steady state) untuk permukaan uji yang dipanaskan secara elektrik adalah sebagai berikut: Q elect = Q conv + Q loss

dimana: Qelect

(2.7)

= laju aliran panas dari listrik (W)

Qconv

= laju perpindahan panas konveksi (W)

Qloss

= laju aliran panas yang hilang dari sistem (W)

Input panas listrik dapat dihitung dari tegangan listrik dan arus listrik yang disuplai ke permukaan. Kehilangan panas (heat loss) dari sistem dapat disebabkan karena radiasi dari permukaan maupun konduksi melalui dinding-dinding saluran ke atmosfer. Sehingga persamaan (2.7) dapat ditulis menjadi: Qelect = Qconv + Q rad + Q cond

dimana: Qrad

= laju perpindahan panas radiasi (W)

Qcond = laju perpindahan panas konduksi (W)

(2.8)

18 Pada penelitian yang serupa, Naik et al (1987) dan Hwang dan Liou (1995) melaporkan bahwa total heat loss radiasi dari permukaan uji yang serupa adalah sekitar 0,5% dari total input panas listrik, sehingga radiative heat loss diabaikan.

Heat loss karena konduksi dari sisi dinding-dinding dapat diabaikan dibandingkan dari permukaan bawah dari seksi uji, karena luas total sisi plat yang dipanaskan jauh lebih kecil dari luas permukaan bawah. Pada penelitian ini, permukaan bawah dari plat uji tidak dipapar ke aliran, dan disolasi dengan kombinasi lapisan isolator dan lapisan kayu, sehingga heat loss konduksi dapat diabaikan. Analisis data akan memuaskan jika persentase total heat loss,

Qelect − Qconv kurang dari Qconv

10% (Naphon, P., 2007). Maka persamaan (2.8) menjadi: Qelect = Qconv

(2.9)

Panas yang dipindahkan dari permukaan bersirip dengan cara konveksi adalah:

⎡ ⎛ T + Tout Qconv = h . As . ⎢Tb − ⎜⎜ in 2 ⎝ ⎣ dimana: Qconv

h

⎞⎤ ⎟⎟⎥ ⎠⎦

(2.10)

= laju perpindahan panas konveksi (W) = koefisien perpindahan panas konveksi rata - rata (W/m2.K)

As

= luas seluruh permukaan yang kontak dengan udara dari susunan sirip pin (m2)

Tb

= temperatur base plate (K)

Tin


Tout


Dari persamaan (2.10), Qconv dapat juga dinyatakan dengan:

Qconv = m& . C p . (Tout − Tin ) dimana: m&

= laju aliran massa udara (kg/s)

Cp

= panas jenis udara (J/kg.K)

Tin


Tout


(2.11)

19 Koefisien perpindahan panas konveksi rata-rata (h) dapat dihitung menggunakan kombinasi persamaan (2.10) dan (2.11), sehingga didapatkan bahwa:

h=

m& . C p . (Tout − Tin )

(2.12)

As . [Tb − ((Tout + Tin ) 2)]

Dari persamaan (2.12), laju aliran massa udara, m& , dapat dihitung dengan persamaan: m& = ρ . A .V

dimana: ρ

(2.13) = massa jenis udara (kg/m3)

A


V


Gambar 2.11 Susunan sirip pin dalam suatu saluran udara segiempat dengan clearance nol

Untuk clearance nol seperti pada gambar 2.11, maka A dihitung dengan rumus: (2.14)

A = H .W b

As adalah luas seluruh permukaan yang kontak dengan udara dari susunan sirip pin atau luas permukaan total dari permukaan plat dasar dan sirip, dan dapat dihitung dengan menggunakan persamaan:

As = Wb . L + π . d . H . N f −

π .d 2 . N f

d⎞ ⎛ As = Wb . L + π . d . N f . ⎜ H − ⎟ 4⎠ ⎝

4

(2.15) (2.16)

20 dimana: Wb

= lebar plat dasar untuk susunan sirip pin (m)

L

= panjang plat dasar untuk susunan sirip pin (m)

d


Nf

= jumlah total sirip pin dalam susunan sirip pin

H

= tinggi saluran udara atau sirip pin (m)

Dari persamaan (2.12), nilai-nilai Tb , Tin dan Tout diukur dari penelitian yang dilakukan menggunakan termokopel. Sedangkan sifat termofisik dari udara Cp dihitung berdasarkan temperatur udara rata-rata, T f = (Tin + Tout ) 2 menggunakan persamaan sebagai berikut:

C p = [9,8185 + 7,7 x 10 −4 (Tin + Tout ) 2] x 10 −2 Persamaan (2.17) berlaku untuk udara pada 250 K ≤

J kg . K

(2.17)

Tin + Tout ≤ 400 K dan pada 2

tekanan atmosfer. Parameter tanpa dimensi yang digunakan dalam perhitungan perpindahan panas untuk permukaan bersirip dihitung sebagai berikut: a. Bilangan Reynolds (Re) Dua jenis bilangan Reynolds digunakan untuk menggolongkan kondisi aliran. Pertama adalah bilangan Reynolds berdasarkan kecepatan rata-rata (V) dalam saluran halus (smooth duct) dan diameter hidrolik dari saluran (Dh) dan dinyatakan dengan:

Re =

V . Dh ν

(2.18)

Re =

ρ .V . Dh µ

(2.19)

dimana: Re

= duct Reynolds number

V


Dh


ν

= viskositas kinematik udara (m2/s)

ρ


µ

= viskositas dinamik udara (kg/m.s)

21 Kedua adalah bilangan Reynolds berdasarkan kecepatan maksimum melalui sirip pin dan ketebalan dari sirip pin, yaitu: Re D =

ρ . Vmaks . d µ

(2.20)

Vmaks dihitung dengan menggunakan persamaan: ⎛ A V maks = ⎜ ⎜ A − A front ⎝

dimana: ReD

Vmaks

⎞ ⎟ .V ⎟ ⎠

(2.21)

= pin Reynolds number

= kecepatan udara maksimum yang melalui sirip pin (m/s)

d


A


Afront

= luas frontal dari sirip-sirip (m2)

ReD telah digunakan secara luas dalam banyak studi perpindahan panas pada sirip pin, dan ReD tergantung pada jumlah pin dalam arah spanwise dan tinggi dari sirip pin. Re di atas biasa disebut sebagai duct Reynolds number dan ReD disebut sebagai pin Reynolds number. b. Bilangan Nusselt (Nu) Seperti pada definisi bilangan Reynolds, bilangan Nusselt rata-rata juga dinyatakan dengan duct Nusselt number dan pin Nusselt number, dimana berturut-turut dinyatakan dengan persamaan:

Nu =

h . Dh k

(2.22)

h.d k

(2.23)

Nu D =

dimana: Nu

= duct Nusselt number

NuD

= pin Nusselt number

Dh


k

= konduktivitas termal udara (W/m.K)

d


h

= koefisien perpindahan panas konveksi rata-rata (W/m2.K)

22 Untuk saluran udara segiempat, diameter hidrolik dihitung dengan persamaan: Dh =

4 .H .Wb 4. A = 2 (H + Wb ) P

(2.24)

Sifat fisik udara, µ dan k dihitung berdasarkan temperatur udara rata-rata, T f = (Tin + Tout ) 2 menggunakan persamaan sebagai berikut:

µ = [4,9934 + 4,483 x 10 −2 (Tin + Tout ) 2] x 10 −6 kg/m.s

(2.25)

k = [3,7415 + 7,495 x 10 −2 (Tin + Tout ) 2] x 10 −3 W/m.K

(2.26)

Persamaan 2.25 dan 2.26 berlaku untuk udara pada 250 K ≤

Tin + Tout ≤ 400 K dan 2

pada tekanan atmosfer. 2.2.7.2 Perhitungan Faktor Gesekan (Friction Factor) Penelitian penurunan tekanan (pressure drop) sepanjang seksi uji dalam saluran bersirip diukur di bawah kondisi aliran panas. Pengukuran ini dikonversi ke faktor gesekan (friction factor), f. Faktor gesekan ditentukan dari nilai pengukuran penurunan tekanan, ∆P, sepanjang seksi uji menggunakan persamaan: f =

∆P ⎡⎛ Lt ⎢⎜⎜ ⎢⎣⎝ D h

⎞⎛ V 2 ⎟⎟⎜ ρ ⎜ ⎠⎝ 2

dimana: f

⎞⎤ ⎟⎥ ⎟ ⎠ ⎥⎦

(2.27)

= faktor gesekan

∆P

= perbedaan tekanan statik (N/m2)

Lt

= panjang seksi uji (m)

Dh

= diameter hidrolik (m)

ρ


V


2.2.7.3 Perhitungan Unjuk Kerja Termal pada Pin-Fin Assembly Peningkatan perpindahan panas dicapai dengan disertai penurunan tekanan. Dalam banyak aplikasi praktis hal tersebut diperbolehkan, sehingga perlu untuk menentukan keuntungan ekonomis karena peningkatan perpindahan panas dan pengaruh sirip pin dan susunannya terhadap unjuk kerja energi menyeluruh dari sistem perpindahan panas sekarang melalui sebuah analisis unjuk kerja termal.

23 Daya pemompaan (pumping power) adalah daya yang dibutuhkan untuk mengalirkan fluida pendingin ke susunan sirip pin, dalam hal ini adalah daya

blower. Untuk sebuah daya pemompaan yang konstan, efektivitas peningkatan perpindahan panas dari permukaan bersirip dibandingkan dengan permukaan halus dapat dihitung dengan menggunakan persamaan sebagai berikut:

V&s . ∆Ps = V&a . ∆Pa

(2.28)

dimana V&s dan V&a berturut-turut adalah laju aliran volumetrik di atas plat tanpa halangan (blocks) dan laju aliran volumetrik di atas plat dengan halangan, sedangkan ∆ Ps dan ∆ Pa berturut-turut adalah penurunan tekanan tanpa halangan dan penurunan tekanan dengan halangan. Menggunakan persamaan Darcy untuk penurunan tekanan dan bilangan Reynolds untuk masing-masing geometri, dari hubungan antara permukaan dengan sirip dan permukaan halus untuk daya pemompaan yang sama, persamaan (2.31) dapat ditulis ulang menjadi:

f s . Res3 = f a . Rea3

(2.29)

Unjuk kerja termal peningkatan perpindahan panas untuk suatu daya pemompaan yang konstan dapat dinyatakan sebagai berikut:

η = (ha hs ) p dimana: ha

(2.30) = koefisien perpindahan panas konveksi dengan sirip (W/m2.K)

hs

= koefisien perpindahan panas konveksi tanpa sirip (W/m2.K)

Jika nilai η ≥ 1, teknik yang dipakai untuk menaikkan laju perpindahan panas adalah menguntungkan dari sudut pandang energi. Jika η ≤ 1, energi yang digunakan untuk menaikkan laju perpindahan panas lebih besar daripada yang diperoleh.

BAB III 3

METODOLOGI PENELITIAN

I

3.1 Tempat Penelitian Penelitian

dilakukan

di

Laboratorium

Perpindahan

Panas

dan

Termodinamika, Jurusan Teknik Mesin Fakultas Teknik Universitas Sebelas Maret Surakarta. 3.2 Alat Penelitian

Gambar 3.1 Skema alat penelitian

Gambar 3.2 Satu set alat uji inline cylindrical pin fin assembly

Alat yang digunakan dalam penelitian ini adalah: a. Saluran udara segiempat (rectangular channel) Saluran udara segiempat terbuat dari kayu yang permukaan bagian dalam dan luarnya dilapisi melamin. Dimensi penampang bagian dalam dari saluran udara segiempat adalah 150 mm x 75 mm x 2.000 mm.

24

25 b. Pelurus aliran udara (flow straightener) Pelurus aliran udara tersusun dari sedotan plastik berdiameter 5 mm, panjang 200 mm sehingga dimensi keseluruhan dari pelurus aliran udara adalah 150 mm x 75 mm x 200 mm.

Gambar 3.3 Pelurus aliran udara (flow straightener)

c. Fan hisap

Fan hisap merupakan modifikasi dari blower, sedemikian rupa sehingga blower dapat menghisap udara.

Gambar 3.4 Fan hisap

d. Rheostat Rheostat digunakan untuk mengatur putaran fan hisap agar didapatkan kecepatan udara yang diinginkan.

Gambar 3.5 Rheostat

26 e. Anemometer Anemometer digunakan untuk mengukur kecepatan aliran udara yang masuk ke dalam saluran udara segiempat.

Gambar 3.6 Anemometer

f. Pemanas listrik (electric heater) Pemanas terbuat dari pita nikelin dengan panjang 4.000 mm, lebar 3 mm dan tebal 1,3 mm yang dililitkan pada kertas mika tahan panas dengan dimensi panjang 200 mm, lebar 150 mm dan tebal 1 mm.

Gambar 3.7 Pemanas listrik (electric heater)

g. Regulator Regulator digunakan untuk mengatur tegangan listrik yang dialirkan ke

heater sehingga temperatur base plate dapat dijaga konstan pada setiap variasi kecepatan dan variasi jarak antar titik pusat sirip pin arah

streamwise.

Gambar 3.8 Regulator pengatur tegangan listrik yang masuk heater

27 h. Multitester digital Multitester digital digunakan untuk mengukur besarnya tegangan listrik yang dibutuhkan pemanas listrik untuk mencapai temperatur base plate yang diinginkan.

Gambar 3.9 Multitester digital

i. Amperemeter Amperemeter digunakan untuk mengukur besarnya arus listrik yang dibutuhkan pemanas listrik untuk mencapai temperatur base plate yang diinginkan.

Gambar 3.10 Amperemeter

j. Manometer U Manometer terbuat dari selang plastik berdiameter 5 mm yang kedua ujungnya ditempatkan pada awal dan akhir dari seksi uji, sehingga dapat mengukur besarnya beda tekanan yang terjadi antara keduanya. Fluida yang digunakan dalam manometer ini adalah solar.

Gambar 3.11 Manometer U dan posisi titik pengukuran tekanan

28 k. Termokopel Menggunakan termokopel tipe T sebanyak 17 buah, dimana 3 buah termokopel digunakan untuk mengukur temperatur udara inlet, 5 buah termokopel untuk mengukur temperatur udara outlet dan 9 buah termokopel untuk mengukur temperatur base plate yang direkatkan dengan lem Araldite.

Gambar 3.12 Termokopel tipe T

l. Thermocouple reader Alat ini digunakan untuk menunjukkan temperatur yang terukur oleh sensor termokopel.

Gambar 3.13 Thermocouple reader

3.3 Spesimen Spesimen berupa pin fin assembly dengan dimensi base plate panjang 200 mm, lebar 150 mm dan tebal 6,5 mm, sedangkan profil sirip adalah silinder pejal dengan diameter sirip 12,7 mm dan tinggi sirip 75 mm. Bahan base plate dan sirip adalah duralumin.

Gambar 3.14 Dimensi dan tata nama spesimen

29

Gambar 3.15 Model spesimen

Spesifikasi spesimen adalah sebagai berikut: Tabel 3.1 Spesifikasi spesimen penelitian

Spesimen 1 2 3 4 5

Sx 37,5 mm 37,5 mm 37,5 mm 37,5 mm

Sy H Wb 25 mm 75 mm 150 mm 30 mm 75 mm 150 mm 37,5 mm 75 mm 150 mm 50 mm 75 mm 150 mm Plat tanpa sirip

L 200 mm 200 mm 200 mm 200 mm

Gambar 3.16 Spesimen 1




30 3.4 Pelaksanaan Penelitian Pengujian dilakukan pada keadaan diam (statis). Sistem beroperasi pada temperatur udara masuk yang konstan sebesar 26 oC dan temperatur base plate yang konstan sebesar 60 oC. Alat penelitian harus dinetralkan terlebih dahulu seperti pada kondisi pengujian sebelumnya (kondisi awal) untuk periode pengujian selanjutnya. Pengambilan data penelitian berupa temperatur dilakukan hingga semua temperatur udara keluar seksi uji mencapai kondisi steady state. Data-data temperatur dan beda tekanan pada keadaan steady state inilah yang akan digunakan untuk analisis data penelitian. Dalam setiap pencatatan data akan diperoleh 17 data temperatur (3 data temperatur udara masuk seksi uji, 5 data temperatur udara keluar seksi uji dan 9 data temperatur base plate). 3.4.1 Tahap Persiapan Memastikan bahwa seluruh alat yang digunakan dalam pengujian, seperti

fan hisap, saluran udara, seksi uji, heater, regulator, rheostat, manometer dan alat pendukung lainnya telah terpasang dengan benar dan berfungsi dengan baik. Memastikan termokopel pengukur temperatur udara masuk dan keluar telah terhubung ke thermocouple reader serta termokopel pengukur temperatur permukaan base plate spesimen telah terpasang dengan benar. 3.4.2 Tahap Pengujian 1. Memasang spesimen ke dalam saluran udara. 2. Menghubungkan semua termokopel pengukur temperatur base plate dengan thermocouple reader. 3. Menyalakan heater sebagai pemanasan awal (preheating) sebesar 70 oC. 4. Menghidupkan fan hisap. 5. Mengatur kecepatan udara sebesar 0,5 m/s dengan mengatur putaran fan menggunakan rheostat. 6. Mengatur temperatur base plate pada temperatur 60 oC. 7. Mencatat seluruh data temperatur dan tekanan setiap 10 menit sampai didapatkan temperatur steady. 8. Mencatat tegangan listrik dan arus listrik yang disuplai ke heater dan fan hisap. 9. Mencatat beda tekanan yang terjadi yang terukur pada manometer U.

31 10. Mematikan heater setelah selesai mengambil data. 11. Mematikan fan. 12. Mengulangi pengujian untuk variasi kecepatan udara yang lain (1 m/s, 2 m/s, 3 m/s, 4 m/s, 5 m/s, 5,5 m/s dan 6 m/s). 13. Mengulangi langkah (3) – (11). 14. Mengulangi pengujian untuk variasi jarak antar titik pusat sirip pin arah

streamwise, Sy yang lain (30 mm, 37,5 mm dan 50 mm). 15. Mengulangi langkah (1) – (11). 16. Mengulangi pengujian untuk spesimen tanpa sirip. 17. Mengulangi langkah (1) – (4). 18. Mengatur daya pemompaan. 19. Mengulangi langkah (6) – (11). 3.5 Proses Analisis Data Berdasarkan data hasil pengujian, yaitu berupa kecepatan aliran udara, temperatur rata-rata udara masuk seksi uji, temperatur rata-rata udara keluar seksi uji, temperatur rata-rata base plate, beda tekanan serta tegangan listrik dan arus listrik yang disuplai ke heater, maka dapat dilakukan perhitungan dan analisis mengenai karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal dari inline cylindrical pin fin assembly. Dari hasil perhitungan tersebut akan dibuat grafik-grafik yang menunjukkan pengaruh jarak antar titik pusat sirip pin dalam arah aliran (streamwise direction, Sy/D) terhadap koefisien perpindahan panas rata-rata (h) dengan bilangan Reynolds (Re), bilangan Nusselt (Nu) dengan bilangan Reynolds, faktor gesekan (f) dengan bilangan Reynolds, dan efisiensi (η) dengan bilangan Reynolds, serta grafik bilangan Nusselt dengan Sy/D. Kemudian dari grafik-grafik tersebut akan dilakukan analisis.

32 3.6 Diagram Alir Penelitian Mulai

Persiapan: Inline cylindrical pin fin assembly

• •

Variasi: Kecepatan aliran udara; 0,5 m/s, 1 m/s, 2 m/s, 3 m/s, 4 m/s, 5 m/s, 5,5 m/s dan 6 m/s. Jarak antar titik pusat sirip pin dalam arah aliran udara (streamwise); 25 mm, 30 mm, 37,5 mm dan 50 mm.

Pengumpulan data: • Temperatur udara masuk seksi uji (Tin), temperatur udara keluar seksi uji (Tout) dan temperatur permukaan plat • Beda tekanan udara masuk dan udara keluar seksi uji (∆P)

Analisis data: • Laju aliran panas dari listrik (Qelect) • Laju perpindahan panas konveksi (Qconv) • Koefisien perpindahan panas konveksi rata-rata (h) • Bilangan Reynolds (Re) • Bilangan Nusselt (Nu) • Faktor Gesekan (f) • Unjuk kerja termal dari Inline cylindrical pin fin assembly (η)

Hasil analisis data: Karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal

Kesimpulan

Selesai

BAB IV 4

DATA DAN ANALISIS

I

Pada bab ini akan dianalisis mengenai pengaruh bilangan Reynolds dan jarak antartitik pusat sirip dalam arah aliran udara (streamwise direction) terhadap karakteristik perpindahan panas dan penurunan tekanan serta unjuk kerja termal dari cylindrical pin fin assembly yang disusun secara segaris (inline) dalam saluran segiempat (rectangular channel). Pengujian dilakukan dengan variasi kecepatan aliran udara masuk antara 0,5 – 6 m/s, dan jarak antar titik pusat sirip dalam arah aliran udara (streamwise

direction, Sy/D) yaitu sebesar 25 mm, 30 mm, 37,5 mm dan 50 mm. Data yang diperoleh dalam pengujian ini, yaitu kecepatan udara masuk, temperatur udara masuk, temperatur udara keluar, temperatur base plate, penurunan tekanan serta tegangan dan arus listrik yang disuplai ke heater dan fan hisap. Sistem dijalankan sampai didapatkan temperatur pada kondisi steady pada tiap variasi pengujian. Proses pengambilan data adalah setiap 10 menit hingga tercapai kondisi steady. 4.1 Data Hasil Pengujian Pengujian

dilakukan

di

Laboratorium

Perpindahan

Panas

dan

Termodinamika Jurusan Teknik Mesin Fakultas Teknik Universitas Sebelas Maret Surakarta. Dari hasil pengamatan temperatur udara masuk seksi uji, temperatur udara keluar seksi uji, temperatur base plate, penurunan tekanan, kecepatan aliran udara masuk serta tegangan listrik dan arus listrik yang disuplai ke heater dan fan hisap saat pengujian pada kondisi steady, diperoleh data seperti pada tabel 4.1 – 4.4 di bawah ini.

(a)

(b)

Gambar 4.1 Posisi titik pengukuran (a) temperatur udara masuk saluran (b) temperatur base plate dan (c) temperatur udara keluar saluran.

33

(c)

34 1. Spesimen 1 Tabel 4.1Data hasil pengujian spesimen 1 (Sx/D = 2,95; Sy/D = 1,97) Kecepatan aliran udara (m/s) 0,5

1

2

3

4

5

5,5

6

Tegangan heater (V)

41

53

70

77

87

93

96

99

Arus heater

2,2

2,8

3,6

4,1

4,5

4,9

5

5,2

56

70

86

98

113

134

155

185

1,1

1,3

1,35

1,65

1,75

1,85

1,95

2,05

Beda tinggi fluida manometer (mm)

0,5

0,75

1,5

2

3,25

4,75

6

7

Tin,1 (oC)

26,2

26,3

26,3

26,4

26,3

26,2

26,3

26,2

26,2

26,2

26,3

26,4

26,3

26,2

26,2

26,2

26,2

26,2

26,2

26,3

26,2

26,2

26,2

26,1

26,2

26,23

26,27

26,37

26,27

26,2

26,23

26,13

60,8

58,9

58,8

56

55,6

53,5

53,3

52,8

61,8

60

59,3

57,4

57,5

58,8

58,2

57,5

58,8

56,9

56,8

56,3

56

55,2

54,6

55,2

60,4

57,1

58,5

59,8

58,7

58,4

58,7

58,2

56,7

57,4

57,7

61,4

61,8

62,2

62,5

62

57,7

61,5

58,7

60,6

59,2

60,6

59,8

59,3

63,3

64,6

63,6

63,4

63,6

63,4

63,8

65,2

62,3

64,7

61,7

63,2

64,4

66,8

67

66,5

60,1

60,3

60,8

62,8

63,2

64,8

64,9

64,4

60,21

60,16

59,54

60,1

60

60,41

60,31

60,3

36,8

35

33,7

32,2

32,2

31,5

31

30,8

36

34,2

33,8

32

31,3

30,8

30,7

30,6

36,3

34,3

33,9

32,3

31,4

31

30,8

30,7

37

35,2

35

33,8

33,1

33

33,2

33

36,6

34,9

33,5

32,1

32,1

31,4

30,9

30,6

36,54

34,8

33,6

32,48

32,02

31,54

31,32

31,14

(A)

Tegangan fan Arus fan

(V)

(A)

o

Tin,2 ( C) o

Tin,3 ( C) o

Tin,rata-rata

( C)

Tbase,1 (oC) o

Tbase,2 ( C) o

Tbase,3 ( C) o

Tbase,4 ( C) o

Tbase,5 ( C) o

Tbase,6 ( C) o

Tbase,7 ( C) o

Tbase,8 ( C) o

Tbase,9 ( C) o

Tbase, rata-rata ( C) Tout,1 (oC) o

Tout,2 ( C) o

Tout,3 ( C) o

Tout,4 ( C) o

Tout,5 ( C) o

Tout, rata-rata ( C)

35 2. Spesimen 2 Tabel 4.2 Data hasil pengujian spesimen 2 (Sx/D = 2,95; Sy/D = 2,36) Kecepatan aliran udara (m/s) 0,5

1

2

3

4

5

5,5

6

Tegangan heater (V)

36

48

61

70

75

80

85

88

Arus heater

2

2,6

3,2

3,6

4

4,2

4,5

4,6

55

69

86

98

113

134

153

186

1,1

1,3

1,35

1,65

1,75

1,85

1,9

2,05


0,3

0,6

1,25

1,75

2,5

3,5

4,25

5,25

Tin,1 (oC)

26,3

26,2

26,3

26,3

26,3

26,3

26

26,3

26,3

26,2

26,2

26,2

26,2

26,3

25,8

26,3

26,3

26,1

26,2

26,2

26,2

26,2

25,8

26,2

26,3

26,17

26,23

26,23

26,23

26,27

25,87

26,27

59,7

60,1

60

60

59,6

59,2

58,9

59

59,2

59,5

59,6

59,3

58,6

58,2

58

58

60,1

60,2

60,7

60,1

59,5

59,5

60

59,5

59,8

59,1

58,8

58,8

59,3

59,5

59,6

60

60,5

60,3

60,6

60,9

60,6

60,6

60,6

60,9

61,8

61,1

61,8

62,9

63,5

64,5

65

64,6

60,4

60,3

60,2

60,2

59

58,8

58,9

60

60,7

60,5

59,8

60

59,3

59,1

59,2

59,6

60,2

59,8

59,8

61,1

60,8

61,1

61,7

61,4

Tbase, rata-rata ( C)

60,43

60,2

59,93

60,43

59,7

59,67

59,93

60,33

Tout,1 (oC)

34,8

33,8

31,9

31,3

30,7

30,3

29,8

30,2

34,4

32,8

31,5

30,4

30

29,6

29,5

29,7

34,5

33

31,8

30,6

30,1

29,8

29,7

29,8

35

34

32,9

31,7

31,2

30,9

30,3

30,6

34,9

33,6

31,8

31,2

30,8

30,4

30

30,3

34,8

33,53

32,17

31,17

30,7

30,37

30

30,23

(A)


(V)

(A)

o

Tin,2 ( C) o

Tin,3 ( C) o

Tin,rata-rata

( C)

Tbase,1 (oC) o

Tbase,2 ( C) o

Tbase,3 ( C) o

Tbase,4 ( C) o

Tbase,5 ( C) o

Tbase,6 ( C) o

Tbase,7 ( C) o

Tbase,8 ( C) o

Tbase,9 ( C) o

o

Tout,2 ( C) o

Tout,3 ( C) o

Tout,4 ( C) o

Tout,5 ( C) o



1

2

3

4

5

5,5

6

Tegangan heater (V)

34

43

56

64

70

73

78

82

Arus heater

1,8

2,4

3

3,4

3,7

3,8

4,1

4,3

56

70

86

98

114

134

154

186

1,1

1,3

1,35

1,65

1,75

1,85

1,95

2,05


0,25

0,5

1

1,5

2,25

3,5

4,25

5,25

Tin,1 (oC)

26,3

26

26,3

26,3

26,2

26,3

26,3

26,3

26,2

26

26,2

26,3

26

26,3

26,2

26,2

26,2

26

26,2

26,2

26

26,2

26,2

26,1

26,23

26

26,23

26,27

26,07

26,27

26,23

26,2

58

59,7

59

59,1

57,1

57,2

58

57,8

61,4

62,6

62,2

62,4

61,4

62

63,3

60,2

57,3

54,3

58,1

58,2

56,6

56,7

58,3

58,5

61,3

62

62

61,5

59

58,5

59,2

60

59,7

60

60

59

57

56,5

57,2

58,7

60,3

61,2

60,6

59,9

60

61,3

62,7

63,6

59,5

60,3

59,6

63,7

64,2

60,6

60,9

60,2

59,8

59,6

58,1

58,8

60,4

56

57,9

57,6

61,8

62

60,1

63

64,6

65

66,9

66,4

(A)


(V)

(A)

o

Tin,2 ( C) o

Tin,3 ( C) o

Tin,rata-rata

( C)

Tbase,1 (oC) o

Tbase,2 ( C) o

Tbase,3 ( C) o

Tbase,4 ( C) o

Tbase,5 ( C) o

Tbase,6 ( C) o

Tbase,7 ( C) o

Tbase,8 ( C) o

Tbase,9 ( C) o


60,37 60,63

59,27

61,83

63,07

60,53

61,9

61,4

Tout,1 (oC)

33,2

32

31,2

30,2

29,6

29,3

29,2

29,3

33

31,4

30,5

29,8

29,4

29,1

29,3

29,2

33,3

31,7

30,7

30,2

29,4

29,1

29,4

29,4

33,6

32,6

31,2

31,5

30,8

30,5

30,7

30,5

33,4

32,3

31,4

30,5

29,7

29,4

29,4

29,4

33,43

32,2

31,1

30,73

29,97

29,67

29,83

29,77

o

Tout,2 ( C) o

Tout,3 ( C) o

Tout,4 ( C) o

Tout,5 ( C) o



1

2

3

4

5

5,5

6

Tegangan heater (V)

32

41

51

59

67

71

73

74

Arus heater

1,7

2,3

2,8

3,1

3,6

3,7

3,8

3,9

56

69

86

98

113

133

154

185

1,1

1,3

1,35

1,65

1,75

1,85

1,9

2,04


0,2

0,3

2

3

2,25

3,5

4

5

Tin,1 (oC)

26,5

26,6

26,2

26,1

26,4

26

26

25,8

26,5

26,6

26,2

26,1

26,3

26

26

25,8

26,5

26,6

26,2

26,1

26,4

26

26

25,8

26,5

26,6

26,2

26,1

26,37

26

26

25,8

58

58

59,5

61,3

59,4

59,3

59,2

59,2

61,1

61

60,8

59,4

60,4

60,6

60,4

60

59,2

58,5

56,9

58,6

57,4

57,8

57,8

57,9

59

61,3

61,9

61

61,8

61,6

61,3

61,1

60,5

59

58,8

58

58,7

58,6

58,4

58,3

59,2

62,4

61

61,2

62,8

64,3

63,4

63,4

61,5

61,3

62

60,8

61,4

61,2

60,8

60,6

61,8

61,2

62

61,2

62,2

62,3

62,2

62

60,5

57,2

57,2

57,1

58,5

58,1

58,2

58

61,27

59,9

60,4

59,7

60,7

60,53

60,4

60,2

33

32,5

30,4

29,4

29,4

28,6

28,5

28,2

32,2

31

29,6

29

29,2

28,4

28,4

28

32,4

31,2

29,8

29,2

29,4

28,7

28,6

28,2

33,2

32,7

32

31,2

31,4

30,4

30,3

30

32,9

32,4

30,2

29,5

29,6

28,9

28,7

28,4

32,83

32,1

30,67

29,97

30,13

29,33

29,2

28,87

(A)


(V)

(A)

o

Tin,2 ( C) o

Tin,3 ( C) o

Tin,rata-rata

( C)

Tbase,1 (oC) o

Tbase,2 ( C) o

Tbase,3 ( C) o

Tbase,4 ( C) o

Tbase,5 ( C) o

Tbase,6 ( C) o

Tbase,7 ( C) o

Tbase,8 ( C) o

Tbase,9 ( C) o

Tbase, rata-rata ( C) Tout,1 (oC) o

Tout,2 ( C) o

Tout,3 ( C) o

Tout,4 ( C) o

Tout,5 ( C) o


38 5. Spesimen 5 Tabel 4.5 Data hasil pengujian spesimen 5 (pelat tanpa sirip) Kecepatan aliran udara (m/s) 0,6

1,1

2,2

3,3

4,4

5,6

6,3

7

Tegangan heater (V)

19

25

34

40

47

52

55

57

Arus heater

1,1

1,5

1,9

2,1

2,5

2,8

3

3,1

56

70

86

98

113

134

154

186

1,1

1,3

1,35

1,65

1,75

1,85

1,95

2,05


0,1

0,15

0,3

0,5

0,75

1

1,2

1,3

Tin,1 (oC)

26,2

26,3

25,9

26,1

25,9

26,3

26

26

26,2

26,2

25,9

26,1

26

26,3

26

25,9

26,1

26,2

25,8

26

25,8

26,2

25,8

25,8

26,17 26,23

25,87

26,07

25,9

26,27

25,93

25,9

(A)


(V)

(A)

o

Tin,2 ( C) o

Tin,3 ( C) o

Tin,rata-rata

( C)

Tbase,1 (oC)

59

58,8

59,2

57,2

57,3

57,4

57,8

58,3

56,5

56,5

56,4

56,9

56,5

55,7

55,9

55,6

63,3

64,1

64,7

64,4

64,6

64,8

65

65,2

60,5

59,8

58,3

58,8

59,1

59,7

59,8

60

54,8

54,4

53,4

52,6

52,3

51,8

51,9

51,3

64,7

66,9

67,1

68,1

68,5

69,1

69,4

70,2

64,3

65

63,6

63,8

63,6

64,1

64

64,4

62

61,9

61,8

61,6

61,7

62

62

62,2

63,67

64,6

64,17

64,5

64,6

65,07

65,13

65,6


26,2

26,3

25,9

26,1

25,9

26,3

26

26,4

Tout,1 (oC)

27,3

27,3

27,1

27,1

27

27

27

27

26,5

26,4

26,4

26,3

26,3

26,3

26,3

26,2

27,4

27,3

27

27

26,8

26,8

26,8

26,7

32,6

32,3

31,1

31

31

31

30,8

30,6

27,4

27,5

27,3

27,2

27,2

27,1

27,1

27,1

28,24 28,16

27,78

27,72

27,66

27,64

27,6

27,52

o

Tbase,2 ( C) o

Tbase,3 ( C) o

Tbase,4 ( C) o

Tbase,5 ( C) o

Tbase,6 ( C) o

Tbase,7 ( C) o

Tbase,8 ( C) o

Tbase,9 ( C) o

o

Tout,2 ( C) o

Tout,3 ( C) o

Tout,4 ( C) o

Tout,5 ( C) o


39 4.2 Perhitungan Data Data spesimen dan seksi uji: Panjang seksi uji (Lt)

= 250 mm = 0,25 m

Tinggi sirip (H)

= 75 mm

Diameter sirip (d)

= 12,7 mm = 0,0127 m

= 0,075 m

Panjang spesimen (L) = 200 mm = 0,2 m Lebar spesimen (Wb)

= 150 mm = 0,15 m

Contoh perhitungan: 1. Spesimen 1 (Sx/D = 2,95; Sy/D = 1,97; Nf = 28) pada kecepatan udara 0,5 m/s Data hasil pengujian: Tegangan heater = Vh = 41 V

Tin, rata-rata = Tin = 26,2 oC = 299,2 K

Arus heater

= Ih = 2,2 A

Tout,, rata-rata = Tout = 36,5 oC = 309,5 K

Tegangan fan

= Vf = 56 V

Tbase, rata-rata = Tb = 60,2 oC = 333,2 K

Arus fan

= If = 1,1 A

Beda ketinggian fluida manometer = h

= 0,5 mm

• Pumping power Pfan = V f . I f . cosϕ

= 56 V x 1,1 A x 0,8 = 49,3 W • Temperatur film Tf =

(T

=

in

+ Tout ) 2

(299,2 + 309,5) K 2

= 304,4 K • Properti udara pada temperatur film

ρ = 1,1469 kg m 3 (tabel Incropera) C p = [9,8185 + 7,7 x 10 −4 x (Tin + Tout ) 2] x 10 2

= [9,8185 + 7,7 x 10 −4 x 304,4] x 10 2 = 1005,29 J kg.K

40 k = [3,7415 + 7,495 x 10 −2 x (Tin + Tout ) 2] x 10 −3 = [3,7415 + 7,495 x 10 −2 x 304,4] x10 −3

= 0,0266 W m.K

µ = [4,9934 + 4,483 x 10 −2 x (Tin + Tout ) 2] x 10 −6 = [4,9934 + 4,483 x10−2 x 304,4] x10−6

= 0,00001864 kg m.s

• Luas penampang melintang saluran udara A = H .Wb

= 0,075 m x 0,15 m = 0,01125 m 2 • Luas total permukaan perpindahan panas d⎞ ⎛ As = Wb . L + π . d . N f . ⎜ H − ⎟ 4⎠ ⎝

0,0127 m ⎞ ⎛ = 0,15 m x 0,2 m + 3,14 x 0,0127 m x 28 x ⎜ 0,075 m − ⎟ 4 ⎝ ⎠

= 0,11 m 2 • Diameter hidrolik saluran udara Dh =

4A P

=

4 .H .Wb 2 (H + Wb )

=

4 x 0,075 m x 0,15 m 2 x (0,075 m + 0,15 m )

= 0,1 m • Laju aliran panas dari heater Qelect = Vh . I h . cosϕ

= 41 V x 2,2 A x 0,8 = 72,2 W

41

• Laju aliran massa udara m& = ρ . A . V

= 1,1469 kg m 3 x 0,01125 m 2 x 0,5 m s = 0,0065 kg s

• Laju perpindahan panas konveksi

Qconv = m& . C p . (Tout − Tin ) = 0,0065 kg s x 1005,29 J kg.K x (309,5 − 299,2 ) K

= 67,1 W • Heat losses yang terjadi pada seksi uji Qloss = =

Qelect − Qconv x 100 % Qconv 72,2W − 67,1 W x 100 % 67,1 W

= 7,6 % • Koefisien perpindahan panas konveksi rata-rata

ha =

=


As . [Tb − ((Tout + Tin ) 2)]

0,0065 kg s x 1005,29 J kg.K x (309,5 − 299,2) K 0,11 m 2 x [333,2 K − ((309,5 + 299,2) K 2)]

= 21,1 W m 2 .K • Bilangan Reynolds Re =

=

ρ .V . Dh µ 1,1469 kg m 3 x 0,5 m s x 0,1 m 2 0,00001864 kg m.s

= 3076,6

42

• Bilangan Nusselt Nu =

=

h . Dh k

21,1 W m 2 .K x 0,1 m 2 0,0266 W m.K

= 79,4 • Penurunan tekanan ∆P = ρ . g . h = 800 kg m 3 x 9,81 m s 2 x 0,0005 m

= 3,92 Pa • Faktor gesekan f =

=

∆P ⎡⎛ Lt ⎢⎜⎜ ⎢⎣⎝ Dh

⎞⎛ V 2 ⎟⎟⎜ ρ ⎜ ⎠⎝ 2

⎞⎤ ⎟⎥ ⎟ ⎠ ⎥⎦ 3,92 Pa

⎡⎛ 0,25 m ⎞⎛ (0,5 m s )2 3 ⎜ ⎜ ⎟ 1 , 1469 kg m x ⎢⎜ ⎟⎜ 2 ⎢⎣⎝ 0,1 m ⎠⎝

= 10,9

⎞⎤ ⎟⎥ ⎟ ⎠⎥⎦

43 2. Spesimen tanpa sirip pada pumping power = 49,3 W Data hasil pengujian: Tegangan heater = Vh = 19 V

Tin, rata-rata = Tin = 26,2 oC = 299,2 K

Arus heater

= Ih = 1,1 A

Tout,, rata-rata = Tout = 28,2 oC = 301,2 K

Tegangan fan

= Vf = 56 V

Tbase, rata-rata = Tb = 60,2 oC = 333,2 K

Arus fan

= If = 1,1 A

Beda ketinggian fluida manometer = h= 0,1 mm

• Temperatur film Tf =

=

(T

in

+ Tout ) 2

(299,2 + 301,2) K 2

= 300,2 K • Properti udara pada temperatur film

ρ = 1,1607 kg m 3 (tabel Incropera) C p = [9,8185 + 7,7 x 10 −4 x (Tin + Tout ) 2] x 10 2

= [9,8185 + 7,7 x 10 −4 x 300,2] x 10 2 = 1004,97 J kg.K

k = [3,7415 + 7,495 x 10 −2 x (Tin + Tout ) 2] x 10 −3 = [3,7415 + 7,495 x 10 −2 x 300,2] x10 −3

= 0,0262 W m.K

µ = [4,9934 + 4,483 x 10 −2 x (Tin + Tout ) 2] x 10 −6 = [4,9934 + 4,483 x10 −2 x 300,2] x10 −6

= 0,00001845 kg m.s

• Luas penampang melintang saluran udara A = H .Wb

= 0,075 m x 0,15 m = 0,01125 m 2

44

• Luas total permukaan perpindahan panas As = L .Wb

= 0,2 m x 0,15 m = 0,03 m 2 • Diameter hidrolik saluran udara Dh =

4A P

=

4 .H .Wb 2 (H + Wb )

=

4 x 0,075 m x 0,15 m 2 x (0,075 m + 0,15 m )

= 0,1 m • Laju aliran panas dari heater Qelect = Vh . I h . cosϕ

= 19 V x 1,1 A x 0,8

= 16,7 W • Laju aliran massa udara m& = ρ . A .V

= 1,1607 kg m 3 x 0,01125 m 2 x 0,6 m s = 0,008 kg s

• Perpindahan panas konveksi

Qconv = m& . C p . (Tout − Tin ) = 0,008 kg s x 1004,97 J kg.K x (301,2 − 299,2 ) K

= 16,3 W

45

• Heat loss yang terjadi Qloss =

=

Qelect − Qconv x 100 % Qconv

16,7 W − 16,3 W x 100 % 16,3 W

= 2,4 % • Koefisien perpindahan panas

hs = =


As . [Tb − ((Tout + Tin ) 2)] 0,008 kg s x 1004,97 J kg.K x (301,2 − 299,2) K 0,03 m 2 x [333,2 K − ((301,2 + 299,2) K 2)]

= 16,5 W m 2 .K • Bilangan Reynolds Re =

=

ρ .V . Dh µ 1,1607 kg m 3 x 0,5 m s x 0,1 m 2 0,00001845 kg m.s

= 3774,4 • Bilangan Nusselt Nu = =

h . Dh k 16,5 W m 2 .K x 0,1 m 2 0,02624 W m.K

= 62,8 • Penurunan tekanan ∆P = ρ . g . h = 800 kg m 3 x 9,81 m s 2 x 0,0001 m

= 0,78 Pa

46 • Faktor gesekan f =

=

∆P ⎡⎛ Lt ⎢⎜⎜ ⎣⎢⎝ D h

⎞⎛ V 2 ⎟⎟⎜ ρ ⎜ ⎠⎝ 2

⎞⎤ ⎟⎥ ⎟ ⎠ ⎦⎥

0,78 Pa 2 ⎡⎛ 0,25 m ⎞⎛ kg (0,6 m s) ⎞⎟⎤ ⎟⎟⎜1,1607 3 x ⎢⎜⎜ ⎜ ⎟⎥ 2 m ⎠⎥⎦ ⎣⎢⎝ 0,1 m ⎠⎝

= 1,50 • Unjuk kerja termal pada pin-fin assembly

η = (ha hs ) p =

21,1 W m 2 .K 16,5 W m 2 .K

= 1,28 Selanjutnya data perhitungan untuk seluruh variasi pengujian dapat dilihat pada lampiran C.

47 4.3 Analisis Data 4.3.1 Pengaruh Bilangan Reynolds dan Jarak Antar Titik Pusat Sirip dalam Arah Streamwise terhadap Karakteristik Perpindahan Panas

Sirip-sirip dipasang secara vertikal pada permukaan uji sehingga memberikan nilai-nilai Sy/D sebesar 1,97, 2,36, 2,95, dan 3,94, sedangkan nilai Sx/D konstan sebesar 2,95. Pengaruh bilangan Reynolds terhadap karakteristik perpindahan panas pada fin pin assembly susunan segaris (inline) dapat dilihat pada gambar 4.2. Karakteristik perpindahan panas pada pin fin assembly dapat dilihat pada hubungan antara koefisien perpindahan panas konveksi (h) dan duct Reynolds number. Gambar 4.2 menunjukkan kelakuan koefisien perpindahan panas konveksi terhadap bilangan Reynolds pada jarak antar titik pusat sirip dalam arah streamwise yang berbeda-beda. 125

h (W/m2.K)

100

75

50

25

0 0

5

10

15

20

25

30

35

40

Re x 103

Gambar 4.2 Pengaruh bilangan Reynolds terhadap koefisien perpindahan panas konveksi rata-rata pada Sx/D = 2,95

Dari gambar 4.2 dapat dilihat bahwa dengan kenaikan bilangan Reynolds, maka nilai koefisien perpindahan panas konveksi semakin besar. Hal ini terjadi pada keseluruhan nilai Sy/D, dimana nilai koefisien perpindahan panas konveksi semakin tinggi untuk nilai Sy/D yang semakin kecil. Semakin besar nilai koefisien perpindahan panas konveksi, maka semakin besar laju perpindahan panas konveksi yang terjadi. Karakteristik perpindahan panas pada pin fin assembly juga dapat dilihat pada hubungan antara duct Nusselt number dan duct Reynolds number. Gambar 4.3 menunjukkan kelakuan bilangan Nusselt rata-rata terhadap bilangan Reynolds

48 pada jarak antar titik pusat sirip dalam arah streamwise yang berbeda-beda untuk susunan sirip segaris (inline). Gambar 4.4 menunjukkan kelakuan bilangan Nusselt rata-rata terhadap jarak antar titik pusat sirip dalam arah streamwise pada bilangan Reynolds yang berbeda-beda untuk susunan segaris. Dari gambar 4.3 dapat dilihat bahwa bilangan Nusselt meningkat dengan kenaikan bilangan Reynolds. Hal ini terjadi pada keseluruhan nilai Sy/D, dimana nilai bilangan Nusselt semakin tinggi untuk nilai Sy/D yang semakin kecil. Ini berarti bahwa dengan semakin kecil jarak antar titik pusat sirip dalam arah streamwise, (Sy), maka semakin besar nilai bilangan Nusselt yang terjadi. Fenomena ini juga terjadi pada penelitian yang dilakukan oleh Tanda (2001). Dengan semakin kecil nilai Sy/D maka jumlah sirip semakin banyak untuk luasan base plate yang sama. Faktor penambahan luasan permukaan perpindahan panas yang berasal dari luasan permukaan perpindahan panas sirip memberikan kontribusi nyata terhadap peningkatan laju perpindahan panas konveksi pada pin fin assembly susunan segaris tersebut. 450 375

Nu

300 225 150 75 0 0

5

10

15

20

25

30

35

40

3

Re x 10

Gambar 4.3 Pengaruh bilangan Reynolds terhadap bilangan Nusselt pada Sx/D = 2,95

Dari gambar 4.2 – 4.4 dapat dilihat bahwa bilangan Reynolds mempunyai pengaruh yang kuat terhadap laju perpindahan panas. Hal ini disebabkan dengan kenaikan laju aliran udara (kenaikan bilangan Reynolds), maka akan menurunkan ketebalan lapis batas (boundary layer) (Bilen, 2002).

49 450 375

Nu

300 225 150 75 0 1,75

2,25

2,75

3,25

3,75

4,25

4,75

Sy/D

Gambar 4.4 Pengaruh nilai Sy/D terhadap bilangan Nusselt pada Sx/D = 2,95

Dari data-data penelitian ini dapat diperoleh korelasi matematis untuk karakteristik perpindahan panas dari pin fin assembly susunan segaris (inline) dengan program SPSS 16. Dari hasil penelitian untuk plat dengan sirip-sirip pin, korelasi antara bilangan Nusselt (Nu) dengan bilangan Reynolds (Re), jarak antar titik pusat sirip (Sy) dan panjang spesimen uji (L) adalah sebagai sebagai berikut:

Nu = 0,214 Re 0,633 (S y L )

−0 , 427

(4.1)

Korelasi perpindahan panas pada persamaan (4.1) berlaku valid untuk range bilangan Reynolds 3.000 < Re < 37.500, L/Dh = 2 dan 1,97 < Sy/D < 3,94. 4.3.2 Pengaruh Bilangan Reynolds dan Jarak Antar Titik Pusat Sirip dalam Arah Streamwise terhadap Karakteristik Penurunan Tekanan

Pengaruh bilangan Reynolds dan jarak antar titik pusat sirip dalam arah streamwise terhadap penurunan tekanan (pressure drop) dan faktor gesekan pin fin assembly susunan segaris berturut-turut dapat dilihat pada gambar 4.5 dan 4.6. Kelakuan penurunan tekanan (∆P) terhadap bilangan Reynolds (Re) serupa dengan hasil penelitian yang dilakukan oleh Yang et al (2007). Dari gambar 4.5 dapat dilihat bahwa penambahan sirip-sirip dengan susunan segaris, menyebabkan penurunan tekanan (pressure drop) yang signifikan dibandingkan dengan permukaan tanpa sirip-sirip (smooth surface). Kelakuan faktor gesekan terhadap bilangan Reynolds pada gambar 4.6 serupa dengan hasil penelitian yang dilakukan oleh Kakac et al (1987). Nilai penurunan tekanan (∆P) dan faktor gesekan (f),

50 semakin menurun dengan kenaikan nilai Sy/D. Hal ini disebabkan dengan semakin besar nilai Sy/D, maka jumlah sirip-sirip pin akan semakin berkurang, sehingga tahanan terhadap aliran udara (resistance to flow) akan semakin berkurang (Bilen, 2002). 70 60

∆P (Pa)

50 40 30 20 10 0 0

5

10

15

20

25

30

35

40

45

3

Re x 10

Gambar 4.5 Pengaruh bilangan Reynolds terhadap penurunan tekanan pada Sx/D = 2,95

Dari gambar 4.6 dapat dilihat bahwa nilai Sy/D lebih berpengaruh dibandingkan bilangan Reynolds terhadap nilai faktor gesekan (f). Hal ini menunjukkan bahwa kenaikan faktor gesekan (f) seiring dengan berkurangnya nilai Sy/D pada dasarnya karena meningkatnya luas permukaan halangan dan efek halangan (blockage effect) akibat kenaikan jumlah sirip-sirip pin. 12 10

f

8 6 4 2 0 0

5

10

15

20

25

30

35

40

45

3

Re x 10

Gambar 4.6 Pengaruh bilangan Reynolds terhadap faktor gesekan pada Sx/D = 2,95

51 Dari data-data penelitian dapat dibuat korelasi matematis antara faktor gesekan (f) yang dihasilkan dengan bilangan Reynolds (Re), jarak antar titik pusat sirip (Sy) dan panjang spesimen uji (L) dengan program SPSS 16 sebagai berikut: f = 2597,024 Re −1,048 (S y L )

−1, 366

(4.2)

Korelasi faktor gesekan pada persamaan (4.2) berlaku valid untuk range bilangan Reynolds 3.000 < Re < 37.500, L/Dh = 2 dan 1,97 < Sy/D < 3,94. 4.3.3 Pengaruh Bilangan Reynolds dan Jarak Antar Titik Pusat Sirip dalam Arah Streamwise terhadap Unjuk Kerja Termal

Dari data penelitian dapat dianalisis mengenai pengaruh jarak antar titik pusat sirip dan susunan sirip terhadap unjuk kerja umum dari sistem dan dapat dievaluasi perolehan energi netto karena penambahan sirip-sirip. Peningkatan perpindahan panas disertai oleh kenaikan penurunan tekanan yang signifikan, dimana dapat mengeliminasi perolehan energi karena peningkatan laju perpindahan panas. Untuk tujuan aplikasi praktis, analisis unjuk kerja termal menjadi sebuah pemikiran yang berguna untuk menentukan perolehan energi netto karena adanya penambahan sirip-sirip. 1.50 1.25

η

1.00 0.75 0.50 0.25 0.00 0

5

10

15

20

25

30

35

40

3

Re x 10

Gambar 4.7 Pengaruh bilangan Reynolds terhadap unjuk kerja termal pada Sx/D = 2,95

Gambar 4.7 menunjukkan hubungan antara unjuk kerja termal (η) dengan dengan bilangan Reynolds (Re) pada jarak antar titik pusat sirip dalam arah streamwise yang berbeda-beda untuk susunan sirip segaris (inline). Perlu ditekankan lagi disini bahwa untuk perolehan energi netto yaitu untuk perpindahan panas yang efektif, nilai η harus lebih besar dari 1 (batas ambang

52 perolehan energi). Dari gambar 4.7 dapat dilihat bahwa nilai η menurun dengan kenaikan bilangan Reynolds (Re), dan nilai η bervariasi antara 0,589 dan 1,28 untuk seluruh Sy/D yang diteliti. Untuk Sy/D = 1,97 pada Re > 24.900, dan Sy/D > 1,97 pada Re > 12.400, nilai η lebih kecil dari 1 dan bervariasi antara 0,589 dan 0,98. Ini berarti bahwa pemakaian pin fin assembly dengan Sy/D = 1,97 pada Re > 24.900, dan Sy/D > 1,97 pada Re > 12.400 akan menyebabkan kehilangan energi daripada perolehan energi. Nilai η lebih besar dari 1 hanya untuk Sy/D = 1,97 pada Re < 24.900, dan Sy/D > 1,97 pada Re < 12.400. Sehingga direkomendasikan untuk memperbaiki efisiensi dari suatu sistem dengan menggunakan pin fin assembly susunan segaris dibatasi pada spesifikasi Sy/D = 1,97 pada Re < 24.900, dan Sy/D > 1,97 pada Re < 12.400. Perolehan energi netto dapat dicapai hingga 28 % untuk nilai Sy/D = 1,97 pada Re = 3.076.

BAB V 5

PENUTUP I

5.1 Kesimpulan

Berdasarkan analisis data dan pembahasan, dapat diambil kesimpulan mengenai pengujian karakteristik perpindahan panas dan penurunan tekanan pada inline cylindrical pin fin array dalam saluran segiempat sebagai berikut: 1. Sirip pin silinder susunan segaris (inline) meningkatkan perpindahan panas dari permukaan base plate sebagai hasil dari kenaikan luasan permukaan perpindahan panas, tetapi disertai dengan adanya penurunan tekanan (pressure drop) yang lebih besar dalam saluran segiempat. 2. Kenaikan bilangan Reynolds (Re) meningkatkan laju perpindahan panas, tetapi menurunkan unjuk kerja termal (η) untuk keseluruhan nilai Sy/D. 3. Faktor gesekan (f) meningkat seiring dengan berkurangnya nilai Sy/D. 4. Penurunan nilai Sy/D meningkatkan unjuk kerja termal (η). 5. Sirip pin silinder susunan segaris (inline) dapat mencapai perolehan energi netto hingga 28 % untuk nilai Sy/D = 1,97 pada Re = 3.076. 5.2 Saran

Berdasarkan pengalaman yang diperoleh dari penelitian tentang pengujian karakteristik perpindahan panas dan penurunan tekanan pada inline cylindrical pin fin array dalam saluran segiempat ini, direkomendasikan beberapa saran sebagai berikut: 1. Temperatur udara lingkungan perlu dijaga dengan lebih baik agar temperatur udara masuk saluran dapat lebih stabil, karena perbedaan temperatur udara lingkungan yang terlalu tinggi sangat berpengaruh terhadap hasil perhitungan. 2. Perlu adanya pengembangan penelitian mengenai pengaruh variasi geometri dan jarak antar titik pusat sirip pin terhadap perpindahan panas, penurunan tekanan, dan unjuk kerja termal dari pin fin assembly.

53

54 DAFTAR PUSTAKA

Bilen K., Akyol, U., and Yapici, S., 2002, Thermal Performance Analysis of A Tube Finned Surface, International Journal of Energy Research, Vol. 26, pp. 321-333. Holman, J.P., 1992, Perpindahan Kalor, Edisi 6, Erlangga, Jakarta. http://hypertextbook.com/physics/matter/density/, diakses 15 April 2009. Incropera, F.P., and DeWitt, D.P., 2007, Fundamentals of Heat and Mass Transfer, 6th Ed, John Willey and Sons, New York. Kakac, S., Shah, R.K., and Aung, W., 1987, Handbook of Single Phase Convective Heat Transfer, John Wiley and Sons, New York. Krauss, A.D., Aziz, A., and Welty, J., 2006, Extended Surface Heat Transfer, 5th Ed, John Wiley and Sons, Inc., England. Lyall, M.E., 2006, Heat Transfer from Low Aspect Ratio Pin Fins, M.S. Thesis, Departement of Mechanical Engineering, Virginia Polytechnic Institute and State University, Blacksburg, Virginia. Naphon P., and Sookkasem A., 2007, Investigation on Heat Transfer Characteristics of Tapered Cylinder Pin fin Heat Sinks, Energy Conversion and Management, Vol. 48 pp. 2671–2679. Sahin B., and Demir, A., 2008, Thermal Performance Analysis and Optimum Design Parameters of Heat Exchanger Having Perforated Pin Fins, Energy Conversion and Management, Vol. 49, pp. 1684-1695. Sahiti, N., Durst, F., and Dewan, A., 2005, Heat Transfer Enhancement by Pin Elements, International Journal of Heat and Mass Transfer, Vol. 48, pp. 4738-4747. Sara, O.N., 2003, Performance Analysis of Rectangular Ducts with Staggered Square Pin Fins, Energy Conversion and Management, Vol. 44, pp. 17871803. Tahat, M., Kodah Z.H., Jarrah, B.A., and Probert, S.D., 2000, Heat Transfers from Pin-Fin Arrays Experiencing Forced Convection, Applied Energy, Vol. 67, pp. 419-442.

55 Tanda, G., 2001, Heat Transfer and Pressure Drop in A Rectangular Channel with Diamond-Shaped Elements, International Jurnal of Heat and Mass Transfer, Vol. 44, pp. 3529-2541. Yang, K.S., Chu, W.H., Chen, I.Y., and Wang, C.C., 2007, A Comparative Study of the Airside Performance of Heat Sinks having Pin Fin Configurations, International Journal of Heat and Mass Transfer, Vol. 50, pp. 4661–4667.

56

LAMPIRAN

Lampiran A. Data hasil pengujian (ada di file Microsoft Excel)

57

77 Lampiran B. Properti berbagai jenis material Tabel B.1 Properti udara pada tekanan atmosfer (sumber: Fundamentals of Heat and Mass Transfer)

Tabel B.2 Massa jenis berbagai material (sumber: www.hypertextbook.com) material acetone acid, acetic (CH3COOH) acid, hydrochloric (HCl) acid, sulfuric (H2SO4) air, 100 K air, 200 K air, 293 K air, 300 K air, 500 K air, 1000 K alcohol, ethyl (grain) alcohol, isopropyl (rubbing) alcohol, methyl (wood) ammonia aluminum argon, gas, ~300 K argon, liquid, 87 K beer, pilsner, 4 °C benzene blood body fat bone butane butter carbon carbon dioxide, gas, +25 °C carbon dioxide, solid, −37 °C copper corn starch, loosely packed corn starch, tightly packed corn syrup diesel formaldehyde freon 12, liquid freon 12, vapor gasoline glycerine gold helium, gas, ~300 K helium, liquid, 4 K hydrogen (H2), gas, 300 K hydrogen (H2), liquid, 17 K honey iron iridium kerosene lard lead lithium lithium 6 deuteride

density (kg/m3) 790 1050 ???? 1390 3.556 1.746 1.207 1.161 0.696 0.340 789.2 785.4 791.3 771 2700 1.449 1430 1008 870 1035 918 1900 551 911 2250 1.799 1101 8960 540 630 1380 800 1130 1311 36.83 803 1260 19,300 0.164 147 0.082 71 1420 7870 22,400 810 919 11,350 534 820

material lungs mayonnaise, traditional mayonnaise, light methane, gas, +25 °C methane, liquid, -90 °C milk, cow, heavy cream milk, cow, light cream milk, cow, whole milk, cow, skim mercury monosodium glutamate nickel nitrogen (N2), gas, ~300 K nitrogen (N2), liquid, 74 K oil, vegetable, coconut oil, vegetable, corn oil, vegetable, olive oil, vegetable, palm oil, vegetable, peanut oil, vegetable, soya osmium oxygen (O2), gas, ~300 K oxygen (O2), liquid, 87 K perchlorethylene platinum plutonium, α salt (sodium chloride) silicon silicon dioxide (quartz) silicone silver skin sodium bicarbonate sugar, sucrose titanium tungsten uranium water, liquid, 100 °C water, liquid, 50 °C water, liquid, 30 °C water, liquid, 20 °C water, liquid, 10 °C water, liquid, 4 °C water, liquid, 0 °C water, ice, 0 °C water, ice, -50 °C water, ice, -100 °C water, sea zinc

density (kg/m3) 400 910 1000 0.656 162 994 1012 1030 1033 13,594 1620 8900 1.145 808 924 922 918 915 914 927 22,500 1.308 1155 1600 21,450 19,860 2165 2330 2600 993 10,490 1050 2200 1550 4500 19,300 19,050 958.40 988.03 995.65 998.21 999.70 999.98 999.84 916 922 927 1025 7140

Lampiran C. Hasil perhitungan data pengujian Tabel C.1 Perhitungan spesimen 1

No. 1 2 3 4 5 6 7 8

V (m/s) 0,5 1 2 3 4 5 5,5 6

Vh

Ih

(volt) 41 53 70 77 87 93 96 99

(A) 2,2 2,8 3,6 4,1 4,5 4,9 5 5,2

Vf

(volt) 56 70 86 98 113 134 155 205

If

(A) 1,1 1,3 1,35 1,65 1,75 1,85 1,95 2,15

Tin

Tout

Tb

(K) 299,2 299,2 299,3 299,4 299,3 299,2 299,2 299,2

(K) 309,5 307,8 306,6 305,5 305,0 304,5 304,3 304,1

(K) 333,2 333,2 332,5 333,1 333,0 333,4 333,3 333,3

Pfan

(W) 49,3 72,8 92,9 129,4 158,2 198,3 241,8 352,6

Tf

(K) 304,4 303,5 302,9 302,4 302,1 301,9 301,8 301,7

ρ

(kg/m3) 1,1469 1,1497 1,1516 1,1533 1,1543 1,1552 1,1555 1,1559

Cp

(J/kg.K) 1005,29 1005,22 1005,18 1005,14 1005,12 1005,09 1005,09 1005,08

k

µ

(W/m.K) 0,02655 0,02649 0,02645 0,02641 0,02639 0,02637 0,02636 0,02635

(kg/m.s) 0,00001864 0,00001860 0,00001857 0,00001855 0,00001854 0,00001853 0,00001852 0,00001852

Tabel C.2 Perhitungan spesimen 1 (lanjutan)

V (m/s) 0,5 1 2 3 4 5 5,5 6

A

(m2) 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125

As

(m2) 0,11 0,11 0,11 0,11 0,11 0,11 0,11 0,11

m&

(kg/s) 0,0065 0,0129 0,0259 0,0389 0,0519 0,0650 0,0715 0,0780

Qelect (W) 72,2 118,7 201,6 252,6 313,2 364,6 384,0 411,8

Qconv (W) 67,1 111,4 191,0 239,2 300,4 348,8 365,5 390,0

Qloss (%) 7,6 6,6 5,5 5,6 4,3 4,5 5,1 5,6

h

(W/m2.K) 21,1 34,1 58,5 70,7 88,3 100,3 105,1 111,8

Re

Nu

3076,6 6181,1 12400,6 18651,3 24905,3 31176,8 34311,4 37455,1

79,4 128,7 221,3 267,8 334,6 380,4 398,9 424,2

∆P (Pa) 3,92 5,89 11,77 15,70 25,51 37,28 47,09 54,94

f

η

10,93 4,10 2,04 1,21 1,10 1,03 1,08 1,06

1,28 1,22 1,18 1,09 0,98 0,92 0,83 0,76

79

Tabel C.3 Perhitungan spesimen 2

No. 1 2 3 4 5 6 7 8

V (m/s) 0,5 1 2 3 4 5 5,5 6

Vh

Ih

(volt) 36 48 61 70 75 80 85 88

(A) 2 2,6 3,2 3,6 4 4,2 4,5 4,6

Vf

(volt) 55 69 86 98 113 134 153 190

If

(A) 1,1 1,3 1,35 1,65 1,75 1,85 1,9 2,05

Tin

Tout

Tb

(K) 299,3 299,2 299,2 299,2 299,2 299,3 298,9 299,3

(K) 307,7 306,4 305,0 304,0 303,6 303,2 302,9 303,1

(K) 333,3 333,1 333,1 333,4 333,0 333,1 333,2 333,3

Pfan

(W) 48,4 71,8 92,9 129,4 158,2 198,3 232,6 311,6

Tf

(K) 303,5 302,8 302,1 301,6 301,4 301,2 300,9 301,2

ρ

(kg/m3) 1,1497 1,1521 1,1544 1,1560 1,1568 1,1573 1,1585 1,1574

Cp

(J/kg.K) 1005,22 1005,17 1005,11 1005,08 1005,06 1005,04 1005,02 1005,04

k

µ

(W/m.K) 0,02649 0,02644 0,02638 0,02635 0,02633 0,02632 0,02629 0,02632

(kg/m.s) 0,00001860 0,00001857 0,00001854 0,00001852 0,00001851 0,00001850 0,00001848 0,00001850


V (m/s) 0,5 1 2 3 4 5 5,5 6

A

(m2) 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125

As

(m2) 0,099 0,099 0,099 0,099 0,099 0,099 0,099 0,099

m&

(kg/s) 0,0065 0,0130 0,0260 0,0390 0,0521 0,0651 0,0717 0,0781

Qelect (W) 57,6 99,8 156,2 201,6 240,0 268,8 306,0 323,8

Qconv (W) 54,7 94,8 150,0 188,5 226,4 257,3 287,7 302,6

Qloss (%) 5,2 5,4 4,1 7,0 6,0 4,5 6,4 7,0

h

(W/m2.K) 18,6 31,7 48,9 60,1 72,5 81,9 90,0 95,3

Re

Nu

3090,7 6204,6 12455,1 18729,2 25004,1 31282,2 34477,9 37546,6

70,3 119,8 185,5 228,2 275,2 311,1 342,5 362,2

∆P (Pa) 2,35 4,71 9,81 13,73 19,62 27,47 33,35 41,20

f

η

6,55 3,27 1,70 1,06 0,85 0,76 0,76 0,79

1,13 1,14 0,98 0,92 0,81 0,75 0,71 0,64

80


No. 1 2 3 4 5 6 7 8

V (m/s) 0,5 1 2 3 4 5 5,5 6

Vh

Ih

(volt) 34 43 56 64 70 73 78 82

(A) 1,8 2,4 3 3,4 3,7 3,8 4,1 4,3

Vf

(volt) 56 70 86 98 114 134 154 186

If

(A) 1,1 1,3 1,35 1,65 1,75 1,85 1,95 2,05

Tin

Tout

Tb

(K) 299,2 299,0 299,2 299,3 299,1 299,3 299,2 299,2

(K) 306,3 305,0 304,0 303,4 302,8 302,5 302,6 302,6

(K) 332,9 333,2 333,0 333,6 333,0 332,3 333,5 333,3

Pfan

(W) 49,3 72,8 92,9 129,4 159,6 198,3 240,2 305,0

Tf

(K) 302,8 302,0 301,6 301,4 300,9 300,9 300,9 300,9

ρ

(kg/m3) 1,1522 1,1547 1,1560 1,1569 1,1583 1,1585 1,1583 1,1585

Cp

(J/kg.K) 1005,16 1005,10 1005,07 1005,05 1005,02 1005,02 1005,02 1005,02

k

µ

(W/m.K) 0,02643 0,02638 0,02635 0,02633 0,02630 0,02629 0,02630 0,02629

(kg/m.s) 0,00001857 0,00001853 0,00001851 0,00001850 0,00001848 0,00001848 0,00001848 0,00001848


V (m/s) 0,5 1 2 3 4 5 5,5 6

A

(m2) 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125

As

(m2) 0,087 0,087 0,087 0,087 0,087 0,087 0,087 0,087

m&

(kg/s) 0,0065 0,0130 0,0260 0,0390 0,0521 0,0652 0,0717 0,0782

Qelect (W) 49,0 82,6 134,4 174,1 207,2 221,9 255,8 282,1

Qconv (W) 46,0 78,3 124,6 163,8 194,5 210,4 242,5 264,1

Qloss (%) 6,4 5,4 7,9 6,3 6,5 5,5 5,5 6,8

h

(W/m2.K) 17,5 28,8 45,5 58,1 69,4 76,7 85,3 93,2

Re

Nu

3102,9 6231,1 12487,5 18757,4 25066,9 31341,9 34468,2 37608,9

66,2 109,1 172,8 220,8 263,9 291,6 324,3 354,4

∆P (Pa) 1,96 3,92 7,85 11,77 17,66 27,47 33,35 41,20

f

η

5,45 2,72 1,36 0,90 0,76 0,76 0,76 0,79

1,06 1,03 0,92 0,89 0,77 0,70 0,67 0,63

81


No. 1 2 3 4 5 6 7 8

V (m/s) 0,5 1 2 3 4 5 5,5 6

Vh

Ih

(volt) 32 41 51 59 67 71 73 74

(A) 1,7 2,3 2,8 3,1 3,6 3,7 3,8 3,9

Vf

(volt) 56 69 86 98 113 132 146 173

If

(A) 1,1 1,3 1,35 1,65 1,75 1,85 1,9 2,03

Tin

Tout

Tb

(K) 299,5 299,6 299,2 299,1 299,4 299,0 299,0 298,8

(K) 305,7 305,0 303,4 302,7 302,8 302,0 301,9 301,6

(K) 333,1 333,0 333,0 332,8 333,3 333,4 333,2 333,1

Pfan

(W) 49,3 71,8 92,9 129,4 158,2 195,4 221,9 281,0

Tf

(K) 302,6 302,3 301,3 300,9 301,1 300,5 300,5 300,2

ρ

(kg/m3) 1,1527 1,1538 1,1571 1,1585 1,1578 1,1597 1,1599 1,1608

Cp

(J/kg.K) 1005,15 1005,13 1005,05 1005,02 1005,03 1004,99 1004,98 1004,96

k

µ

(W/m.K) 0,02642 0,02640 0,02632 0,02629 0,02631 0,02626 0,02626 0,02624

(kg/m.s) 0,00001856 0,00001854 0,00001850 0,00001848 0,00001849 0,00001846 0,00001846 0,00001845


V (m/s) 0,5 1 2 3 4 5 5,5 6

A

(m2) 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125

As

(m2) 0,076 0,076 0,076 0,076 0,076 0,076 0,076 0,076

m&

(kg/s) 0,0065 0,0130 0,0260 0,0391 0,0521 0,0652 0,0718 0,0784

Qelect (W) 43,5 75,4 114,2 146,3 193,0 210,2 221,9 230,9

Qconv (W) 40,7 69,9 109,9 139,9 179,8 196,7 209,2 217,3

Qloss (%) 7,0 7,9 4,0 4,6 7,3 6,9 6,1 6,2

h

(W/m2.K) 17,6 30,0 45,7 57,7 73,6 78,8 84,2 87,2

Re

Nu

3105,3 6221,8 12508,4 18804,5 25045,6 31403,9 34553,5 37748,7

66,6 113,7 173,6 219,4 279,7 299,9 320,8 332,1

∆P (Pa) 1,57 2,35 5,89 11,77 17,66 27,47 31,39 39,24

f

η

4,36 1,63 1,02 0,90 0,76 0,76 0,72 0,75

1,0684 1,0781 0,9192 0,8867 0,8209 0,7219 0,6635 0,5892

82


No. 1 2 3 4 5 6 7 8

V

Vh

(m/s) 0,6 1,1 2,2 3,3 4,4 5,6 6,3 7

(volt) 19 25 34 40 47 52 55 60

Ih

(A) 1,1 1,5 1,9 2,1 2,5 2,8 3 3,3

Vf

(volt) 56 70 86 98 113 134 154 186

If

(A) 1,1 1,3 1,35 1,65 1,75 1,85 1,95 2,05

Tin

Tout

Tb

(K) 299,2 299,2 299,1 299,2 299,1 299,2 299,1 298,9

(K) 301,2 301,2 300,8 300,7 300,7 300,6 300,6 300,5

(K) 333,2 333,5 333,1 333,0 333,4 333,0 333,2 333,3

Pfan

(W) 49,3 72,8 92,9 129,4 158,2 198,3 240,2 305,0

Tf

(K) 300,2 300,2 299,9 300,0 299,9 299,9 299,8 299,7

ρ

(kg/m3) 1,1607 1,1607 1,1617 1,1615 1,1618 1,1617 1,1620 1,1624

Cp

(J/kg.K) 1004,97 1004,97 1004,94 1004,95 1004,94 1004,94 1004,94 1004,93

k

µ

(W/m.K) 0,02624 0,02624 0,02622 0,02622 0,02622 0,02622 0,02621 0,02620

(kg/m.s) 0,00001845 0,00001845 0,00001844 0,00001844 0,00001844 0,00001844 0,00001843 0,00001843


V (m/s) 0,6 1,1 2,2 3,3 4,4 5,6 6,3 7

A

(m2) 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125 0,01125

As

(m2) 0,03 0,03 0,03 0,03 0,03 0,03 0,03 0,03

m&

(kg/s) 0,0078 0,0144 0,0288 0,0431 0,0575 0,0732 0,0824 0,0915

Qelect (W) 16,7 30,0 51,7 67,2 94,0 116,5 132,0 158,4

Qconv (W) 16,3 27,8 49,5 64,4 90,2 108,4 126,9 149,0

Qloss (%) 2,4 7,9 4,4 4,3 4,3 7,5 4,0 6,3

h

(W/m2.K) 16,5 27,8 49,7 65,1 89,7 109,1 126,9 147,9

Re

Nu

3774,4 6920,0 13860,0 20784,1 27726,4 35283,8 39708,9 44149,9

62,8 106,1 189,6 248,1 342,0 416,1 484,2 564,5

∆P (Pa) 0,78 1,18 2,35 3,92 5,89 7,85 9,42 10,20

f 1,50 0,67 0,34 0,25 0,21 0,17 0,16 0,14

83

PENGUJIAN KARAKTERISTIK PERPINDAHAN PANAS DAN PENURUNAN TEKANAN DARI INLINE CYLINDRICAL PIN FIN ARRAY DALAM RECTANGULAR CHANNEL

Recommend Documents