Makalah Konveksi - Kelompok 9

MAKALAH PERPINDAHAN KALOR

Pemicu 2: Perpindahan Kalor Konveksi

Disusun oleh : Kelompok 9

Abu Bakar Ash Shiddiq

(1306449302) (1306449302)

Fachreza Maulana I.

(1406643085) (1406643085)

Farah Moulydia

(1306370650) (1306370650)

Pangiastika Putri Wulandari (1306370493) (1306370493) Yolla Miranda

(1306414841) (1306414841)

UNIVERSITAS INDONESIA DEPOK 2015

Makalah Perpindahan Kalor Konveksi KATA PENGANTAR

Puji syukur kehadirat Tuhan Yang Maha Esa yang telah melimpahkan rahmat dan hidayah-Nya sehingga kami dapat menyelesaikan makalah ini dengan baik. Penulisan makalah ini bertujuan untuk memenuhi tugas mata kuliah Perpindahan Kalor. Dalam makalah ini, kami ingin memaparkan jawaban Pemicu 2 mengenai Perpindahan Kalor Konveksi yang terdiri dari perpindahan kalor konveksi secara alami dan paksa. Kami mengucapkan terima kasih kepada Ibu Dianursanti dan Ibu Tania Surya Utami, sebagai dosen mata kuliah Perpindahan Kalor. Tidak lupa juga kami mengucapkan terima kasih kepada semua pihak yang telah membantu dan memberikan dukungan kepada kelompok kami sehingga dapat menyelesaikan makalah ini. Semoga makalah ini dapat bermanfaat bagi semua pihak. Kami menyadari dalam pembuatan karya tulis ini masih terdapat kekurangan. Oleh karena itu, kritik dan saran dari para pembaca yang membangun sangat kami harapkan.

Depok, April 2015

Penulis

Kelompok 9 – PK 02

i

Makalah Perpindahan Kalor Konveksi

DAFTAR ISI

Kata Pengantar ............................................. ................................................................... ............................................ ................................... .............

i

Daftar Isi ........................................... ................................................................. ............................................ ............................................ ........................ ..

ii

Teori Dasar............................ Dasar................................................... ............................................. ............................................ ................................... .............

1

Jawaban Pemicu ........................................... ................................................................. ............................................ ................................... .............

15

Kesimpulan .......................................... ................................................................ ............................................ ........................................... .....................

34

Daftar Pustaka .......................................... ................................................................ ............................................ ....................................... .................

35


ii


TEORI DASAR

KONVEKSI ALAMIAH 1. Prinsip-prinsip Konveksi 1.1

Lapisan Batas Laminar pada Plat Rata

Persamaan gerakan untuk lapisan batas dapat kita turunkan dengan membuat neraca gaya dan momentum pada unsure volume itu. Untuk menyederhanakan analisis pada gambar 5-4 dapat kita misalkan seperti berikut : a. Fluida tak mampu merapat dan aliran tunak b. Tidak terdapat perubahan tekanan di arah tegak lurus pada plat c. Viskositas tetap d. Gaya geser viskos diarah y dapat diabaikan Persamaan ini ialah persamaan momentum integral untuk lapisan batas. Jika tekanan

 =0=~ ~ …1  =  …2  = √ 4,64 …3

di seluruh aliran tetap, maka

Persamaan lapisan batas integral menjadi

Setelah disederhanakan didapatkan angka Reynold yaitu

1.2

Persamaan Energi Lapisan Batas

Gambar 1. unsur volume untuk analisis energi lapisan

Untuk menyederhanakan analisis dalam penyusunan neraca energi, buat asumsi sebagai berikut :  Aliran tunak tak mampu mampat  Viskositas, konduktivitas termal  Konduksi kalor pada arah x diabaikan


1

Makalah Perpindahan Kalor Konveksi Neraca energi yang dapat disusun adalah : Energi konveksi permukaan kiri + energi konveksi permukaan bawah + kalor konduksi permukaan bawah + kerja viskos netto pada unsur = energi konveksi pada muka kanan + energi konveksi permukaan atas + kalor dikonduksi pada permukaan atas

   …4   …5     …6

Formulasi gaya geser viskos adalah :

Jarak perpindahan per satuan waktu terhadap unsur volume kendali dx dy adalah :

Formulasi energi viskos netto yang dihasilkan adalah :

1.3

Lapisan Batas Turbulen

Pada fluida yang mengalir diatas plat datar, daerah turbulen adalah daerah yang berada jauh dari dinding plat. Apabila posisi fluida yang ditinjau semakin jauh dari dinding plat maka semakin turbulen lah daerah tersebut. Pada daerah turbulen variabel intensif fluida yang berperan sudah bersifat makroskopik, seperti viskositas Eddy dan konduktivitas termal Eddy. Dalam perpindahan kalor, batas bilangan Reynold aliran laminar adalah 5 x 10 5. Untuk nilai yang lebih besar dari batas tersebut, aliran dikategorikan turbulen. Seperti yang tampak pada gambar 9, untuk daerah yang dekat dengan permukaan plat, aliran memiliki karakteristik aliran laminar. Untuk daerah yang semakin jauh di atasnya, aliran mulai memasuki kategori daerah buffer / batas transisi. Semakin jauh aliran dari permukaan plat, aliran mulai bersifat turbulen sempurna, dan daerah ini ditunjukan oleh arsiran berwarna biru. Daerah ini juga disebut sebagai lapisan batas turbulen. Cara perhitungan ketebalan lapisan batas turbulen ini didasarkan pada 2 kondisi, yaitu: 1. Ketika lapisan batas berbentuk turbulen sempurna mulai dari ujung awal plat

... (7)


2

Makalah Perpindahan Kalor Konveksi 2. Ketika lapisan batas baru berbentuk turbulen mulai dari batas transisi

... (8) Cara perhitungan perpindahan panas pada aliran turbulen menggunakan analogi friksi-fluida. Untuk kasus plat isotermal dengan rentang 5 x 10 5 < ReL < 107 dan Recrit = 5 x 10 5 persamaan bilangan Nusseltnya adalah:

... (9)

Untuk rentang yang lebih tinggi yaitu 10 7 < ReL < 109 dan Recrit = 5 x 10 5 persamaannya menjadi:

... (10)

1.4

Lapisan Aliran Laminar dalam Tabung

Perpindahan kalor dalam kondisi aliran berkembang bila aliran itu tetap laminar. Suhu dinding adalah T w, jari-jari tabung adalah r o, dan kecepatan pada pusat tabung u0. Kita mengandaikan tekanan seragam pada setiap titik pada penampang

   1     1     = ∝  4   4   …11 ℎ= 1124 = 1148 …12  = ℎ =4,364 …13

tabung. Distribusi suhu yang diperoleh setekah dilakukan penurunan rumus adalah :

Suhu L imbak

Setelah melakukan penurunan rumus akan diperoleh persamaan berikut :

1.5

Aliran Turbulen dalam Tabung

Pada daerah dekat permukaan terdapat suatu sub-lapisan laminar, atau “film”, sedangkan inti tengah aliran bersifat turbulen. Untuk mengetahui distribusi suhu, kita harus menganalisis pengaruh pusaran turbulen dalam perpindahan kalor dan


3

Makalah Perpindahan Kalor Konveksi momentum. Untuk perhitungan, persamaan yang tepat untuk digunakan pada aliran

 =0,023 ,, …14

turbulen dalam tabung licin adalah sebagai berikut :

2. Sistem Konveksi Alamiah 2.1

Konveksi Bebas dari Bola

Rumus Empiris untuk perpindahan kalor konveksi bebas dari bola ke udara sebagai berikut

5 Untuk 1 < Gr f <10

 = ℎ =20,392  …15

Persamaan diatas dapat diubah dengan memasukan angka prandtl sehingga didapatkan

 =20,43 …16

Persamaan ini berlaku untuk perhitungan konveksi bebas pad gas. Akan tetapi dapat digunakan untuk zat cair apabila tidak ada informasi khusus untuk itu. Hasil dari perkalian angka garshof dan prandtl yang rendah, angka Nusselt mendekati 2. Nilai inilah yang didapatkan pada konduksi murni melalui fluida stagnan tak berhingga yang mengelilingi bola itu. Untuk rentang angka Rayleigh yang lebih tinggi, hasil eksperimen dari Amato dan Tien dengan air dapat ditunjukkan dengan persamaan berikut :

 =20,5 …17

Untuk 3 x 105 < Gr Pr < 8 x 10 8

KONVEKSI PAKSA 1. Rumus-rumus Empiris pada Konveksi Paksa 1.1

Rumus-rumus Empiris untuk Aliran dalam Pipa dan Tabung Gambar 2. Pengaruh pemanasan pada profil kecepatan aliran laminar dalam tabung (Sumber: Holman, J.P. 2010)

Untuk menghitung variasi sifat-sifat, Sieder dan Tate men yarankan rumus berikut:


4


 =0,027 , ,

... (18)

Semua sifat-sifat ditentukan pada suhu-limbak kecuali µ w ditentukan pada suhu dinding. Pada bagian inlet, dimana aliran belum terbentuk jenis alirannya (atau masih berkembang), Nusselt menyarankan rumus:

dengan L

,     ,     =0,036    adalah panjang



untuk 10 < < 400

tabung, dan d

adalah diameter

... (19) tabung. Sifat pada

persamaan ini ditentukan pada suhu rata-rata. Persamaan-persamaan

di

atas

memungkinkan

tetapi kesalahannya bisa mencapai

perhitungan

yang

sederhana,

25%. Petukhov mengembangkan persamaan

yang lebih teliti, tetapi lebih rumit. Untuk aliran yang sepenuhnya turbulen dalam tabung licin dirumuskan:

   = ,+,−  

... (20)

dengan n = 0,11 untuk T w> T b, n = 0,25 untuk T w< T b, dan n = 0 untuk fluks kalor tetap dan gas. Semua sifat ditentukan pada T f = (Tw + T b)/2, kecuali untuk µ w (viskositas pada dinding) dan µ b (viskositas pada suhu limbak). Faktor gesek

 =1,82log  1,64−  ,  (  =3,66 +, [())]        =1,86 Pr    

( friction factor ) untuk tabung licin persamaannya adalah:

... (21)

Hausen memberikan rumus empiris untuk aliran laminar yang berkembang penuh ( fully developed laminar flow) pada tabung dengan suhu tetap, yaitu: ... (22)

Persamaan empiris untuk perpindahan kalor aliran laminar dalam tabung yang dirumuskan oleh Sieder dan Tate adalah: (

1.2

...(23)

Aliran Menyilang Silinder dan Bola M ekani sme Per pin dahan Panas Konveksi Paksa

Fenomena pemisahan lapisan batas ditunjukkan dalam gambar 8 (Lampiran). Selama aliran fluida bergerak sepanjang bagian depan silinder ataupun bola, tekanan akan berkurang untuk kemudian meningkat lagi pada bagian belakang silinder. Hal ini menyebabkan bertambahnya kecepatan aliran bebas pada bagian depan silinder,


5

Makalah Perpindahan Kalor Konveksi dan berkurangnya kecepatan di bagian belakang silinder. Kenaikan tekanan dan penurunan kecepatan ini dihubungkan oleh persamaan Bernoulli, apabila ditulis  2  sepanjang garis aliran: dP  d  u   2 g    c 

Karena tekanan di seluruh lapisan-batas dianggap tetap, maka akan terlihat bahwa aliran balik akan bermula pada lapisan-batas dekat permukaan, artinya, momentum lapisan-lapisan fluida di dekat permukaan tidak cukup tinggi untuk bisa mengatasi peningkatan tekanan. Bila gradien kecepatan pada permukaan menjadi nol, maka aliran tersebut dikatakan mencapai titik pisah (gambar 9, Lampiran). Bila aliran bergerak terus melewati titik pisah, maka terjadilah fenomena aliran balik. Pada akhirnya, daerah aliran yang terpisah pada bagian belakang silinder menjadi turbulen dan bergerak secara acak. Proses terbentuknya aliran ini amat mempengaruhi perpindahan kalor dari silinder panas ataupun bola banas ke aliran fluida. Fenomena perpindahan kalor ini diselidiki oleh Giedt dan hasilnya dirangkum pada gambar 10 (Lampiran). Pada awal aliran, ketika bilangan reynolds masih rendah (70.800 dan 101.300), nilai koefisien perpindahan kalor konveksi h berada pada titik minimum. Kemudian ketika aliran mengalir ke bagian belakang silinder terjadi peningkatan nilai h karena adanya gerakan pusaran yang turbulen pada aliran yang kemudian terpisah. Perpindahan kalor meningkat dengan cepat ketika lapisan-batas menjadi turbulen, dan ketika terjadi peningkatan gerakan pusaran pada pemisahan aliran.

H ubungan Empiri s 

Silinder Menurut Knudsen dan Katz , nilai koefisien perpindahan kalor rata-rata h aliran gas dan zat cair dicari dengan angka Nusselt, yang dituliskan dalam bentuk persamaan berikut:

 u d  Nu   C     v  k f  f  h d

n 1/ 3

Pr

untuk 10-1 < Re f < 105

...(24)

dimana nilai konstanta C dan eksponen n digunakan untuk menggambarkan susunan berkas tabung, dan nilainya dapat dilihat dalam tabel 3 (Lampiran), dan sifat-sifat fluida dievaluasi pada suhu film. 

Silinder Tak Bundar


6

Makalah Perpindahan Kalor Konveksi Perpindahan kalor pada silinder tak bundar dapat dihitung dengan persamaan perpindahan kalor secara umum, dimana nilai h dicari dengan angka Nusselt pada persamaan 8. Konstanta yang digunakan pada persamaan ini dapat dilihat pada tabel 2 (Lampiran). 

Bola Dari berbagai persamaan yang ada, Whitaker merumuskan persamaan tunggal untuk gas dan zat cair yang mengalir melintasi bola: Nu = 2 + (0,4 Red 1/2 + 0,06 Red 2/3) Pr 0,4 ( μ∞/ μw)1/4

...(25)

yang berlaku untuk rentang 3,5 < Re d < 8x104 dan 0,7 < Pr <380.

1.3

Aliran Menyilang Rangkunan Tabung

Karakteristik Perpindahan kalor pada rangkunan tabung yang segaris atau selang seling dipelajari oleh Grinson dan atas dasar korelasi data dari berbagai peneliti, ia berhasil menyajikan datanya dalam bentuk persamaan (6-17) dan nilai

ℎ =~  …617

Konstanta C dan di daftar 6-4

Daftar 6-4 korelasi Grimson untuk perpindahan kalor dalam rangkunan tabung 10 baris atau lebih untuk digunakan 6-17


7


Angka Reynolds didasarkan pada kecepatan maksimum yang terjadi pada rangkunan tabung, yaitu kecepatan melalui bidang aliran yang minimum. Luas bidang ini bergantung dari susunan geometri tabung. Nomenklator pada daftar 6-4

~

ditunjukkan pada gambar 6-14. Jika jumlah tabung dalam baris lebih sedikit, maka perbandingan nilai h untuk baris tabung terhadap baris 10 tabung diberikan pada daftar 6-5.

Penurunan tekanan untuk aliran gas melintas rangkunan tabung dapat dihitung pada persamaan 6-31

 

Gmaks

  ,    2 ∆=    …631

: Kecepatan massa pada luas bidang aliran minimum (Kg/m2s)

: Densitas ditentukan pada kondisi aliran bebas ( Kg/m3) : Jumlah baris melintang : Viskositas aliran bebas rata-rata

Faktor gesek empiris f’ adalah

 =0,25  0,118,−, …632


8


 ={0,044 [ 0,08]/,+,/}−, …633

Untuk baris selang seling dan untuk baris segaris

Zukauskas memberikan informasi tambahan untuk berkas tabung dengan memperhitungkan rentang angka Reynolds yang luas dan perbedaan perbedaan sifat. Persamaan korelasinya mempunyai bentuk

1.4

/ ℎ   ,   =  =   

Perpindahan Kalor Logam-Cair

…634

Dalam beberapa tahun terakhir banyak perhatian dicurahkan kepada perpindahan kalor logam-cair karena tingginya laju perpindahan kalor yang dapat dicapai dengan media ini. Hal ini disebabkan oleh tingginya konduktivitas termal logam cair dibandingkan dengan fluida lain dan sebagai akibatnya, logam cair sangat sesuai untuk situasi di mana sejumlah besar energi harus harus dikeluarkan dari ruang yang relatif kecil. Selain itu, logam cair masih tetap berada dalam keadaan cair pada suhu yang lebih tinggi dibandingkan dengan fluida konvensional seperti air dan bahan bahan pendingin organik. Angka Prandtl untuk logam cair sangat rendah, sekitar 0,01, sehingga tebal lapisan batas termal jauh lebih besar dari lapisan batas hidrodinamik. Hal ini disebabkan oleh nilai konduktivitas termal yang tinggi pada logam cair. Oleh karena perbandingan t kecil, profil kecepatan berbentuk sangat tumpul pada sebagian besar lapisan-batas termal. Asumsi yang dapat digunakan yaitu :

/

=~

Pada keseluruhan lapisan-batas untuk menghitung suku transpor energi dalam persamaan energi internal adalah :

 [ ∫ ~ ]=  []  …26

2. Alat Penukar Kalor

Alat penukar panas merupakan suatu alat yang berfungsi untuk memindahkan panas antara dua fluida yang berbeda suhu melalui suatu sekat. Berdasarkan jenis aliran fluida yang digunakan alat penukar panas dibagi atas beberapa jenis, yaitu One through system, Closed recirculating system, dan Open recirculating system. Sedangkan jenis alat penukar


9

Makalah Perpindahan Kalor Konveksi panas berdasarkan bentuk alatnya adalah Tipe pipa rangkap, Tipe Selongsong-Tabung, dan Tipe plate and frame. 2.1

Faktor Pengotor

Setelah dipakai beberapa lama, permukaan perpindahan kalor pada HE dapat dilapisi oleh berbagai endapan atau permukaan itu mengalami korosi karena adanya interaksi antara fluida dengan bahan yang digunakan. Kedua hal tersebut dapat memberikan tahanan tambahan terhadap aliran kalor sehingga menurunkan kinerja HE. Pada shell-and-tube heat exchanger, fouling dapat terjadi baik pada bagian dalam maupun luar tube dan dapat terjadi pada bagian dalam shell . Pengaruh menyeluruh dari hal ini dapat dinyatakan secara matematis dengan fouling factor atau tahanan pengotoran, R f , yang harus diperhitungkan bersama tahanan termal lainnya dalam menghitung koefisien perpindahan kalor menyeluruh. Faktor pengotoran harus didapatkan dari percobaan, yaitu dengan menentukan U untuk

 = 1  1  …27

kondisi bersih dan kondisi kotor pada penukar kalor tersebut.

Nilai faktor pengotoran yang disarankan untuk berbagai fluida dapat dlihat pada buku Holman.

Gambar 3. Pipa (kanan) baru dipakai, (kiri) permukaan dalam pipa terkotori (Sumber : alat penukar panas.pdf)

2.2

Beda Suhu Rata-rata (LMTD)

Temperatur fluida panas maupun fluida dingin yang masuk heat exchanger biasanya selalu berubah-ubah. Untuk menentukan perbedaan temperatur tersebut digunakan perbedaan temperatur rata-rata atau Logarithmic Mean Temperature Difference (LMTD). LMTD

digunakan

dalam

perhitungan-perhitungan

heat

exchanger yang menunjukkan panas yang dipindahkan.


10


∆ = ln[ / ] …28 2.3

METODE NTU-EFEKTIVITAS

Pendekatan LMTD dalam analisis penukar kalor berguna bila suhu masuk dan suhu keluar diketahui atau dapat ditentukan dengan mudah, sehingga LMTD dapat dengan mudah dihitung dan aliran kalor, luas permukaan dan koefisien perpindahan kalor menyeluruh dapat ditentukan. Bila kita harus menentukan suhu masuk atau suhu keluar, analisis akan melibatkan prosedur iterasi karena LMTD merupakan fungsi logaritma. Efektivitas HE didefinnisikan sebagai :

     = ℎℎ    

Untuk perpindahan kalor aktual :

Q = mhch(Th1 – Th2) = mccc(Tc2-Tc1) (untuk HE aliran sejajar) ...(29) Q = mhch(Th1 – Th2) = mccc(Tc1-Tc2) (untuk HE aliran berlawanan) ...(30) Perpindahan kalor maksimum yang mungkin : Qmaks = (mc)min(Thmasuk - Tcmasuk ) 2.4

...(31)

Penukar Kalor Kompak

Hal terpenting dari penukar kalor kompak adalah memiliki luas permukaan yang sangat besar per satuan volume. Penukar kalor kompak sangat cocok untuk penerapan dalam aliran gas di mana nilai h rendah. Angka Stanton dan angka Reynold didasarkan atas kecepatan massa pada luas penampang aliran minimum dan diameter hidrolik.

=  …32 =   …33

Perbandingan luas aliran bebas terhadap luas frontal, yaitu :

2.5

Pertimbangan Perancangan Heat Exchanger

Secara umum, pertimbangan pada perancangan alat penukar kalor adalah :  Persyaratan Perpindahan kalor  Biaya  Ukuran fisis  Pressure Drop Characteristics


11

Makalah Perpindahan Kalor Konveksi Persyaratan perpindahan kalor harus selalu dipenuhi dalam setiap pemilihan atau perancangan penukar kalor. Pertimbangan yang cermat pada perancangan Heat Exchanger pada faktor-faktor di atas akan menghasilkan rancangan Heat Exchanger yang baik pula.


12

Makalah Perpindahan Kalor Konveksi JAWABAN PEMICU Bagian A Contoh Kasus

Dapatkah anda menggambarkan dan menjelaskan mekanisme perpindahan kalor yang terjadi pada peristiwa angin laut dan angin darat, serta persamaan-persamaan konveksi yang terlibat dalam penjelasan mekanisme tersebut?

Jawaban :

Proses terjadinya angin darat dan angin laut disebabkan oleh beda sifat fisis antara permukaan darat dan laut. Yaitu perbedaan sifat antara daratan dan lautan dalam menyerap dan melepaskan energi panas matahari. Daratan menyerap dan melepas energi panas lebih cepat daripada lautan. Periode angin darat dan angin laut adalah harian. 1. Angin Laut (the sea breeze)

Angin laut terjadi ketika pada pagi hingga menjelang sore hari, daratan menyerap energi panas lebih cepat dari lautan sehingga suhu udara di darat lebih panas daripada di laut. Akibatnya udara panas di daratan akan naik dan digantikan udara dingin dari lautan. Maka terjadilah aliran udara dari laut ke darat. 2. Angin Darat (the land breeze)


13

Makalah Perpindahan Kalor Konveksi Angin darat terjadi ketika pada malam hari energi panas yang diserap permukaan bumi sepanjang hari akan dilepaskan lebih cepat oleh daratan (udara dingin). Sementara itu di lautan energi panas sedang dalam proses dilepaskan ke udara. Gerakan konvektif tersebut menyebabkan udara dingin dari daratan bergerak menggantikan udara yang naik di lautan sehingga terjadi aliran udara dari darat ke laut.

Terjadinya angin darat dan angin laut juga melibatkan perpindahan kalor secara konduksi dan konveksi. Kalor jenis daratan (kalor jenis benda padat) lebih kecil daripada kalor jenis air laut, karenanya daratan lebih cepat panas ketika disinari matahari dan juga lebih cepat dingin ketika malam hari tiba. Daratan

yang

lebih

cepat

panas,

memanaskan udara yang berada di atas (kalor berpindah dari daratan ke udara secara konduksi). Suhu udara yang mendapat tambahan kalor meningkat dan udara memuai. Akibatnya massa jenis udara tersebut berkurang. Sebaliknya suhu air laut lebih dingin sehingga udara yang berada di atas permukaan air laut juga lebih dingin, dibandingkan udara yang berada di permukaan daratan. Udara di permukaan laut lebih dingin sehingga massa jenisnya lebih besar. Adanya perbedaan massa jenis udara menyebabkan udara yang berada di permukaan laut bergerak menuju daratan dan menekan udara di daratan ke atas. Semakin jauh dari permukaan bumi, jumlah udara semakin berkurang karena gaya gravitasi bumi semakin kecil. Karena jumlah udara semakin berkurang maka tekanan udara juga semakin kecil. Udara panas di daratan yang bergerak ke atas mengalami pendinginan karena semakin jauh dari permukaan bumi, tekanan udara semakin berkurang. Udara dingin tersebut kemudian bergerak lagi ke bawah, tidak menuju permukaan daratan tetapi menuju permukaan lautan yang mempunyai suhu yang lebih dingin. Proses ini terjadi secara terus menerus sehingga timbul aliran udara dari laut menuju darat. Singkatnya, udara di dekat permukaan laut bergerak ke daratan, udara di dekat permukaan daratan bergerak ke atas, udara yang berada di atas bergerak ke permukaan laut.


14

Makalah Perpindahan Kalor Konveksi Persamaan Konveksi Alami Pada Peristiwa Angin Darat dan Angin Laut

Asumsi : 

Permukaan darat dan laut adalah permukaan tak berhingga



Permukaan darat dan laut adalah permukaan bola (bumi)





Konveksi yang terjadi adalah konveksi alami vertikal dan horizontal dari permukaan ke udara Konveksi terjadi pada tekanan 1 atm

Persamaan yang dapat digunakan adalah :

 =20, 3 92

5 Gr f1/4 untuk 1 < Gr f < 10



Nuf =



Apabila subtitusi dengan angka Prandtl , didapat : 1/4 Nuf = 2 + 0,43 (Gr f Pr f)

Bagian B 2

1. Sebuah kolektor sinar matahari, berbentuk plat rata berukuran 1 m , terletak miring dengan sudut 20 terhadap horizontal. Permukaan panas berada pada suhu o

160 C dan tekanan 0,1 atm. Sejajar di atas permukaan panas tersebut, dipasang jendela transparan yang berfungsi melewatkan energi radiasi dari matahari. Jarak antara jendela transparan dengan permukaan panas adalah 8 cm. Suhu jendela transparan dipertahankan pada 40 C. Hitunglah perpindahan kalor konveksi alami yang terjadi antara permukaan panas dengan jendela transparan.

Jawaban:

Diketahui : T jendela (T∞) = 40oC T pada plat panas (Tw) = 160oC Tekanan = 0,1 atm A plat panas = 1 m 2 Jarak antara kedua plat = 8 cm = 0.08 m θ terhadap sumbu horizontal = 20o θ terhadap sumbu vertikal Ditanya

= 90o - 20o= 70o

:

Q antara permukaan panas dengan jendela transparan = ?


15

Makalah Perpindahan Kalor Konveksi Penyelesaian

:

 Fluida yang terperangkap adalah udarra  Perpindahan kalor konveksi alami terjadi pada ruang tertutup plat miring horizontal

 16040  = 2  =100=373.15

Melihat sifat-sifat udara pada daftar Daftar A.5, Holman hal.589 untuk suhu 375.15K didapatkan dari hasil interpolasi,

  2.172 10



k  0.0317 W  

1

T f



kg

5

m. s 0

C

m.

2.68 10



3

1



K

Pr  0.69

 

P . Mr R.T f

10100 Pa



287

J kg . K

Menghitung nilai

. 



0.0946kg / m3

.373.15 K

,

    .  .    .    .  =  .   − 9. 8  0. 0 946 1 0    .   . 2 . 6 8. 4 33. 1 5313. 1 5 .  0. 0 8 . 0 . 6 9 2 3  .  = − 1 0 2.172.  .. .=2.11.10 =20 =0.212

Sudut kemiringan plat terhadap sumbu horizontal

.

Menggunakan data 7.3 hal.320 didapatkan,


16


= 14  =  =0.21221100.cos20/ =2.516   =2.516     =2.516     =0,0797  0,0317

Besar perpindahan kalor konveksi alami yang terjadi antara permukaan panas dengan

     0, 0 797 1 6040 .  1             =  = 0.08 = 

jendela transparan:

2. Sebuah bola berdiameter 2,5 cm berada pada suhu 38°C, akan dibenamkan ke dalam suatu wadah yang berisi air dengan suhunya 15°C. a. Bagaimana anda menjelaskan pengaruh dimensi dan ukuran wadah tersebut terhadap mekanisme perpindahan kalor yang terjadi pada sistem di atas? 3

b. Jika wadah yang digunakan adalah suatu bejana berukuran 8x7x6 cm , bagaimana anda menentukan laju perpindahan kalornya? c. Apa yang menjadi pertimbangan anda dalam menentukan persamaan empiris yang akan digunakan untuk menyelesaikan problem di atas?

Jawaban :

a. Dimensi wadah pada perpindahan kalor dalam kasus ini tidak berpengaruh. Yang justru berpengaruh adalah jumlah air yang ada dalam bejana tersebut, dimana jumlah

 = =ℎ∆=∆′

air (massa) berhubungan dengan volume air, bukan volume bejana, karena:

Bagaimana pun bentuk wadahnya, yang menjadi pertimbangan adalah volume airnya, sebab dalam perpindahan kalor pada sistem ini dapat dinyatakan dengan:

b. Untuk menghitung laju perpindahan kalor, tahapan perhitungan yang dilakukan adalah:


17

Makalah Perpindahan Kalor Konveksi 1) Membuat sistem ( shell )

q- (dalam-luar)

2) Menghitung suhu film

3) Menentukan sifat air sebagai fluida yang didapat dari Daftar A.9 Buku Holman.

4) Menghitung GrPr

5) Dengan nilai GrPr demikian, berlaku:

c. Yang menjadi pertimbangan dalam menentukan rumus empiris seperti di atas adalah: 1) Sistem konveksi bebas pada bola 2) Berlaku untuk fluida air sesuai eksperimen Amato dan Tien


18

Makalah Perpindahan Kalor Konveksi 3) Rentang GrPr pada soal memenuhi syarat eksperimen Amato dan Tien untuk fluida air (3 x 10 5 < GrPr < 8 x 10 8)

0

3. Sebuah silinder vertical dengan tinggi 1,8 m; diameter 7,5 cm; dan suhu 93 C, 0

berada dalam lingkungan dengan suhu 30 C a. Hitunglah kalor yang dilepas melalui konveksi alami dari silinder ini b. Dapatkah silinder tersebut diperlakukan sebagai sebuah plat rata vertical? Berapakah diameter minimum yang harus dimiliki oleh silinder tersebut agar dapat diasumsikan sebagai sebuah plat rata vertical? c. Jika silinder tidak dapat dianalogikan dengan plat rata vertical, bagaimanakah menyelesaikan permasalahan diatas? Jawaban :

Diketahui: 

D= 0.075 m



Tw = 93 C



T



g = 9.8 m/s2

⁰ ∞ ⁰

= 30 C

Asumsi: 





Plat rata vertical bukan merupakan silinder vertical berdiri di atas tanah sehingga luas permukaan kontak langsung dengan fluida kondisi dalam keadaan tunak

a. Untuk menentukan laju perpindahan kalor, mula mula menentukan suhu filmnya

 =  2  = 93302 =61,5℃=334,5  −/ =19, 1 88  10 =0, 0 288 /℃ =  =2,99  10−− Pr= 0,689 terlebih dahulu yaitu:

Pada 334,5 K sifat sifat udara (pada Daftar A.5, Holman hal 589) dari hasil interpolasi didapatkan data:




19


          =  −     1 , 8 = 9,82,9919,10188 109330 − =2, 9 4  10   −   =0,5080,952Pr  −    =0,=154,508270,6890,9520,698 2,94  10   =  ℎ=  ℎ= 154,271,08,0288 ℎ=2,46832  ℃  =ℎ         =2, 4 6832 1 , 8 9 330 =503,83 



Maka nilai Grashof adalah



Maka nilai





Nilai

maka

Maka nilai perpindahan kalor adalah

b. Kriteria umum Silinder vertical dapat dianalogikan sebagai plat rata vertical apabila

 ≥ 35   = 7,51,108 − =0,042 35 = 2,94 35 10 =0,085   ≤ 

memenuhi syarat berikut:

Dengan menggunakan data yang sama kita dapat menghitung

Dapat dilihat bahwa

 

maka dari itu hal ini tidak memenuhi syarat plat rata

vertical. Sehingga harus menggunakan silinder vertical.


20

Makalah Perpindahan Kalor Konveksi Untuk memenuhi syarat plat rata vertical kita harus mencari D minimum yang dapat memenuhi persyaratan

 ≥ 35   =0,085 1,8 =0,085 =0,085  1,8 =0,153

Sehingga didapatkan diameter minumum untuk memenuhi persyaratan plat rata vertical sebesar 0,153 m.

c. Dengan menggunakan rumus untuk silinder vertical, mula mula menentukan terlebih

= 0,689 =2, 9 4  10 =2,03 10 / 0, 3 87 /  =0,825 [10,492/Pr ]/       0, 3 87 2 , 0 3  10   =0,825 [10,492/0,689 ]  =17,732 =314, 4 2 =  ℎ=   = 314,421,08,0288 ℎ=5,0307 / ℃ =ℎ 

dahulu nilai Ra

Nilai Ra yang didapat masuk kedalam range 10-1
Maka

adalah

Sehingga perpindahan kalor yang terjadi adalah


21


=5,03070,075  1,8  4 0,0759330 =135,74 

4. Suhu pada suatu permukaan dinding vertical 4 ft x 10 ft dipertahankan konstan 0

0

530 F sedangkan suhu udara sekeliling 70 F dan tekanan 1 atm a. Hitunglah kalor yang hilang dari permukaan dinding itu secara konveksi bebas ke udara b. Jika dinding itu disekat dengan bahan penyekat yang tebalnya 2 inc dan daya hantar panasnya (konduktivitas termal) = 0,121 BTU/jam ft

2 0

F. Hitunglah

kalor yang hilang secara konduksi dan konveksi bebas bila dianggap suhu pada 0

permukaan penyekat 250 F Jawaban :

Diketahui: 

4 ft x 10 ft



Tw = 530 F



T



g = 9.8 m/s2

⁰ ∞ ⁰

10ft

= 70 F

4ft

Asumsi: 

Plat Vertikal



Terjadi proses konduksi dan konveksi



Tinggi permukaan dinding vertical 10 ft



Lebar permukaan dinding vertical 4 ft



Bidang permukaan adalah plat vertical isothermal



Penyekat dipasang di satu sisi

a. Untuk menentukan laju perpindahan kalor, mula mula menentukan suhu filmnya

 =  2  = 53070 2 =300℉=421,85 

terlebih dahulu yaitu:

Pada 421,85 K sifat sifat udara (pada Daftar A.5, Holman hal 589) dari hasi l interpolasi didapatkan data:


22


−/ =2, 8 4  10 =0, 0 351 /℃ =  =2,37  10−− Pr= 0,686           =  −  = 9,82,37 102,845103070   3, 0 48 − =3,64  10 = = 3,64  10  0,686=2,49 10  ⁄   0, 3 87   ⁄   ̅ =0,825 [1  0,492⁄⁄]⁄ ̅⁄ =0,825 [1 0,3870,4920,2,⁄4968610⁄⁄]⁄ ̅ = ̅ (ℎ ) ̅ = ̅  …30 ℎ ℎ̅ = 680,833,00,048351W/m. ℃ =7,84 /℃    ̅   …31 =ℎ        =7,8453070=,  / 




Perhitungan Bilangan Rayleigh (Ra)



Maka nilai

Setelah mendapatkan bilangan Rayleigh (Ra), kita dapat mencari bilangan

Nusselt (Nu) untuk 10-1< RaL < 1012

680,83



Perhitungan Koefisien Perpindahan Kalor Konveksi Bebas Rata-rata



Perhitungan Kalor yang hilang secara konveksi alami dari dinding ke udara

b. K = 0,121 BTU/jam ft 2 0F X = 2 inch = 0,0508 m

2”


23

Makalah Perpindahan Kalor Konveksi 

Konduksi



Konveksi

= ∆∆∆  = 0,∆121530250  ==,0,0/ 508  

Untuk menentukan laju perpindahan kalor, mula mula menentukan suhu filmnya

 =  2  = 25070 2 =160℉=344,1  −/ =0,0295 /℃ =2, 0 16  10 =  =2,91  10−− Pr= 0,698          =  −  = 9,82,91 210,0162105070   3 , 0 48 − =3,576  10 = = 3,576  10  0,698=2,49 10 terlebih dahulu yaitu:

Pada 421,85 K sifat sifat udara (pada Daftar A.5, Holman hal 589) dari hasi l interpolasi didapatkan data:






Perhitungan Bilangan Rayleigh (Ra)



Perhitungan Bilangan Nusselt (Nu)

Setelah mendapatkan bilangan Rayleigh (Ra), kita dapat mencari bilangan Nusselt (Nu) untuk 10-1< RaL < 1012

̅⁄ =0,825 [1  0,03,48792⁄⁄⁄ ]⁄ ⁄     0, 3 87 2 , 4 9  10 ⁄   ̅ =0,825 [1  0,4920,⁄ 698⁄]⁄ ̅ = (ℎ̅) 703.52



Perhitungan Koefisien Perpindahan Kalor Konveksi Bebas Rata-rata


24




̅ = ̅  …30 ℎ  ℃  ℎ̅ = 703,523, 00,048295W/m. = 6, 8 1 / ℃   =ℎ̅   …31   = 6,8125070 =,  /   =       =, ,  =,  /

Perhitungan Kalor yang hilang secara konveksi alami dari dinding ke udara

5. Suatu alat pengukur panas dipakai untuk memanaskan sesuatu zat alir dari suhu 50°F dengan kecepatan W lb/jam yang menyebabkan aliran turbulen. Alat pengukur panas tersebut terdiri atas n buah pipa dengan diameter D ft dan panjang L ft. Jika kemudian dalam keadaan sama (W sama) pipa-pipa pada alat pengukur diganti dengan pipa berdiameter 0,5 D sedang jumlahnya tetap n pipa, hitung berapa % perubahan panjang pipa untuk mendapatkan pemanasan yang sama. Dalam hal ini dianggap h = U serta sifat zat alir tetap.

Jawaban :

Berdasarkan pernyataan dalam soal, beberapa hal yang dapat diasumsikan: 1) Karena koefisien kalor menyeluruh sama dengan koefisien konveksi, dapat diasumsikan ahwa perpindahakan kalor yang diperhitungkan hanyalah konveksi 2) Sistem terjadi pada aliran tabung 3) Aliran turbulen berkembang penuh

Sesuai daftar 6-8 buku Holman, untuk sistem dan asumsi di atas, dengan n pemanasan adalah 0,4, rumus yang berlaku yaitu:

 =0,023 0,8

Prasyarat dalam menggunakan rumus di atas haruslah memiliki: 0,6 < Pr < 100. Dengan melihat Lampiran A Daftar A-5, pada suhu 50°F = 283 K, didapat Pr = 0,713. Artinya, rumus di atas memang berlaku untuk soal ini.


25


 =0,023 0,80,4 0, 8 ℎ    1  =  =0,0230,8    0,4 ℎ1 =0,023  0,40,2 0, 8   ℎ2 =0,023    0,40,50,2 0, 8  0, 2 0, 4   0, 0 23   ℎℎ12 = 0,023 0,8 0,40,50,2 ℎ1 ℎ2 =0,87  = 1 2 ℎ10,1∆=ℎ 0, 5  ∆ 2 2 827=1, =0,7451 1,74 11 1 100%=74%

Untuk sistem pertama:

Untuk sistem kedua:

Maka perbandingan kedua sistem:

Karena laju perpindahan kalor di kedua sistem sama:

6. Dalam sebuah alat penukar kalor aliran silang, digunakan gas panas (Cp = 1,09 o

o

o

kJ/kg. C) untuk memanaskan 2,5 kg/detik air dari suhu 35 C menjadi 85 C. Gas o

o

masuk pada suhu 200 C dan keluar pada suhu 93 C. Koefisien perpindahan kalor 2 o

menyeluruh sebesar 180 W/m . C. Hitunglah luas penukar kalor dengan menggunakan : a. Pendekatan LMTD b. Metode NTU-Efektivitas

Jawaban :

Diketahui : Alat Penukar Kalor Aliran Silang (Cross Flow) Cpgas panas = 1,09 kJ/kg. oC


26


  ̇ 

= 200oC = 93oC = 85oC = 35oC = 2,5 kg/s = 180 W/m 2.oC

Ditanya : Luas alat penukar kalor dengan metode : a. LMTD b. NTU – efektivitas

Asumsi :

1. Jenis alat penukar kalor : Penukar kalor aliran silang sekali lintas, satu fluida campur dan satu fluida tak campur 2. Cpair = 4,18 kJ/kg. oC

Penyelesaian :

a. Luas Alat Penukar Kalor Menggunakan metode LMTD

̇ =  ∆    =2,5 4180 . ° 8535 ° =522.500 

 Menentukan perpindahan kalor total dari energi yang diserap oleh air :

 Menentukan beda suhu rata-rata (LMTD) :

 

= 200oC = 93oC = 85oC = 35oC


∆ = ln[ / ] ℃ 2 0085 ℃ ∆ = ln9[335 9335 ℃/20085 ℃] =83,27℃ 27

Makalah Perpindahan Kalor Konveksi  Mencari Faktor koreksi menggunakan grafik yaitu Gambar 10-11 pada buku

Heat Transfer Holman :

Dari Gambar 10-11, t 1 dan t2 menunjukkan fluida tak campur (air) dan T 1 dan T2 menunjukkan fluida campur (gas panas), sehingga : T1 = 200oC T2 = 93oC t1 = 35oC t2 = 85oC

=   8535 =0,303 = 20035 =   = 20093 8535 =2,14 =0,91

Maka faktor koreksinya adalah :


28


 Mencari luas yang diperlukan alat penukar kalor

=∆   = ∆  = 180 . ℃522.0,5009183,27℃ =38,31 

Dengan demikian, berdasarkan perhitungan menggunakan metode LMTD diperoleh luas alat penukar kalor sebesar 38,31 m2.

b. Luas Alat Penukar Kalor Menggunakan metode NTU – efektivitas 

 =   ̇ ̇   ∆ =  ∆       2,5 4180 . ° 8535 ° = ̇ 20093 ° ⁄ ̇̇ =4883, 1 8   =10450 ⁄   ⁄  ==̇ ̇=4883, 1 8   ⁄ =10450       =0,467 ⁄  = 4883, 1 8  10450 ⁄   ∆ 20093  = ∆ = 20035 =0,648

Menentukan Cmin dan Cmaks

Maka, fluida minimum adalah gas panas :



Menentukan nilai Jadi :



Menentukan nilai efektivitas


29




Menentukan nilai NTUmaks dari Gambar 10-14 pada Buku Holman

Fluida tak campur adalah air (C maks) dan fluida campur gas panas (C min), sehingga diperoleh :



 =1,35  =   .  1, 3=54883,  ⁄ 1 8     = 180 . ℃   =36,6 

Menentukan Luas Alat Penukar Kalor :

Berdasarkan perhitungan menggunakan metode NTU-efektivitas diperoleh luas alat penukar kalor sebesar 36,6 m2.


30

Makalah Perpindahan Kalor Konveksi 7. Sebuah sistem pemanas air menggunakan alat penukar kalor jenis selongsongo

tabung. Uap panas mengalir dalam satu lintasan selongsong pada suhu 120 C, o

sedangkan air masuk pada suhu 30 C dan melakukan empat lintasan tabung 2 o

dengan nilai U = 2000 W/m . C. a. Hitunglah luas penukar kalor, jika aliran air yang masuk sebesar 2,5 kg/detik o

dan air keluar pada suhu 100 C. b. Jika setelah beroperasi selama beberapa waktu alat penukar kalor tersebut 2 o

mengalami faktor pengotoran sebesar 0,0002 m . C/W, berapakah suhu air yang keluar pada kondisi tersebut?

Jawaban :

Diketahui :

  = ̇ 

Alat penukar kalor jenis selongsong-tabung =

120oC

= 100oC = 30oC

= 2,5 kg/s = 2000 W/m 2.oC

2o R f = 0,0002 m . C/W

Ditanya : a. Luas alat penukar kalor b. Suhu air yang keluar

Asumsi :

1. Jenis alat penukar kalor : Selongsong tabung (1 lintas selongsong dan 4 lintas tabung) 2. Cpair = Cpuap = 4,18 kJ/kg. oC 3.

̇ ̇ =

= 2,5 kg/s

4. Menggunakan metode LMTD

Penyelesaian : a. Luas Alat Penukar Kalor Menggunakan metode LMTD


31


̇ =  ∆    =2,5 4180 . ° 10030 ° =731.500 

 Menentukan perpindahan kalor total dari energi yang diserap oleh air :

 Menentukan beda suhu rata-rata (LMTD) :

 

= 120oC = 120oC = 100oC = 30oC

∆ = ln[ / ]    ℃  1 20100 ℃ ∆ = ln1[2030 12030 ℃/120100 ℃] =46,54℃

 Tidak dapat mencari faktor koreksi menggunakan grafik yaitu Gambar 10-4 pada

buku Heat Transfer Holman, karena suhu uap saat masuk dan keluar adalah sama, yaitu 120oC.

 Mencari luas yang diperlukan alat penukar kalor

=∆   = ∆  = 2000 731.. ℃500 46,54℃ =7,859 

Dengan demikian, berdasarkan perhitungan menggunakan metode LMTD diperoleh luas alat penukar kalor sebesar 7,859 m2.

b. Suhu air keluar jika faktor pengotor 0,0002 m2.oC/W  Mencari U (Koefisien Perpindahan Kalor) :

 = 1  1  0,0002= 1  20001

Faktor pengotoran didefinisikan sebagai berikut :

Sehingga :


32


 Mencari Cmin :

1 =0,0002 20001 1 = 100007  =1428,57     = ⁄ =̇⁄   =2,5 =10450   4180⁄  =     1 428, 5 7   7 , 8 59    = 10450 ⁄  =1,074 − =1 −, =1 =0,658 ∆  =0,658 = 12030 ∆ =59,22 ∆  =59,22   =3059,22 = 89,22

 Mencari nilai NTU :

NTU didefinisikan sebagai :

 Mencari nilai efektivitas alat penukar kalor :

Persamaan efektivitas alat penukar kalor :

Maka :

 Menentukan suhu air yang keluar :

Jadi, suhu air yang keluar pada alat penukar kalor adalah sebe sar 89,22oC.


33

Makalah Perpindahan Kalor Konveksi KESIMPULAN

1. Prinsip dasar dari konveksi paksa adalah adanya suatu alat yang memaksa kalor untuk berpindah (melalui suatu fluida), maka perpindahan kalor yang diinginkan dapat berlangsung lebih cepat dan efektif. 2. Analisis lapisan batas dibagi menjadi 3 kategori yaitu lapisan batas laminar, lapisan batas termal, dan lapisan batas turbulen. 3. Analisis lapisan batas adalah suatu metode untuk mencari nilai koefisien perpindahan panas konveksi melalui bilangan Nusselt dari sistem dengan keadaan yang berbeda beda (aliran laminar/turbulen, plat isotermal, fluks panas konstan, fluida dengan Pr > 100). 4. Perpindahan kalor konveksi secara alamiah dapat terjadi pada plat rata vertikal dan horizontal, silinder vertikal dan horizontal, permukaan miring, dan bola. 5. Prinsip kerja dari alat penukar kalor yaitu memindahkan panas dari 2 fluida pada temperatur berbeda di mana transfer panas dapat dilakukan secara langsung ataupun tidak langsung. 6. Kemampuan suatu alat penukar kalor (APK) dalam memindahkan panas dipengaruhi oleh beberapa faktor, antara lain: faktor pengotoran, koefisien perpindahan panas, luas permukaan permindahan panas, beda suhu rata-rata, jumlah lintasan, dan penurunan tekanan APK. 7. Fouling factor adalah akumulasi dari mineral ataupun karat dalam suatu alat penukar

kalor yang akan mengganggu kinerja alat. Selain fouling factor , penurunan tekanan juga akan mempengaruhi kinerja dengan cara melambatkan aliran. 8. Metode untuk menentukan laju kalor pada konveksi paksa yang terdapat pada Alat Penukar Kalor (APK) adalah metode beda suhu rata-rata (LMTD) dan metode NTUefektivitas.


34

Makalah Konveksi - Kelompok 9

Recommend Documents