PERPINDAHAN PANAS (HEAT TRANSFER)

Download Konveksi, Radiasi). 2. Aplikasi perpindahan panas dalam Industri. Dasar-dasar mempelajari perpindahan panas: • Persamaan differensial biasa...

0 downloads 710 Views 535KB Size
PERPINDAHAN PANAS (HEAT TRANSFER) Luqman Buchori, ST, MT Jurusan Teknik Kimia Fakultas Teknik UNDIP Semarang

REFERENSI 1. 2. 3. 4. 5. 6.

Kern, D.Q., “Process Heat Transfer”, International Student Edition, McGraw Hill Kogakusha, Ltd., New York. Holman, J.P., “Heat Transfer”, sixth edition, McGraw Hill, Ltd., New York, 1986. Mikheyev, M., “Fundamentals of Heat Transfer”, John Willey & Sons Inc., New York, 1986. Incopera De Witt, “Fundamentals of Heat Transfer”, John Willey & Sons Inc., New York, 1981. Ozisik, “Heat Transfer, a basic approach”, 1984. McAdams, W.H., “Heat Transmision”, 3rd edition, McGraw Hill Book Company, Inc., New York.

MATERI KULIAH 1. Dasar-dasar perpindahan panas (Konduksi, Konveksi, Radiasi). 2. Aplikasi perpindahan panas dalam Industri Dasar-dasar mempelajari perpindahan panas: • Persamaan differensial biasa/parsial • Mekanika fluida • Konsep neraca energi thermodinamika

Definisi : Ilmu yang mempelajari tentang laju perpindahan panas diantara material/benda karena adanya perbedaan suhu (panas dan dingin)

Panas akan mengalir dari tempat yang suhunya tinggi ke tempat yang suhunya lebih rendah

KEGUNAAN ILMU PERPINDAHAN PANAS Z Untuk merencanakan alat-alat penukar panas (heat exchanger). Z Untuk menghitung kebutuhan media pemanas/ pendingin pada suatu reboiler atau kondensor dalam kolom destilasi. radiasi Z Untuk perhitungan furnace/dapur. Z Untuk perancangan ketel uap/boiler. Z Untuk perancangan alat-alat penguap (evaporator). Z Untuk perancangan reaktor kimia – –

Eksotermis Endotermis

butuh pendingin butuh pemanas

MEKANISME PERPINDAHAN PANAS 1. Konduksi (hantaran) 2. Konveksi 3. Radiasi (sinaran)

1. KONDUKSI Adalah proses perpindahan panas jika panas mengalir dari tempat yang suhunya tinggi ke tempat yang suhunya lebih rendah, dengan media penghantar panas tetap. Dasar : Hukum Fourier ⎛ dT ⎞ q k = k A ⎜⎜ − ⎟⎟ ⎝ dx ⎠

atau

qk A

=k

⎛ dT ⎞ ⎟ ⎜− ⎜ dx ⎟ ⎠ ⎝

Contoh perpindahan panas konduksi

Perpindahan panas konduksi pada bahan dengan ketebalan berbeda, mana yang lebih lama naik suhunya ?

Perpindahan panas konduksi pada bahan dengan panjang berbeda, mana yang lebih lama panasnya ?

Perpindahan panas konduksi pada bahan dengan ∆ suhu berbeda, mana yang lebih cepat konduksinya ?

2. KONVEKSI Yaitu perpindahan panas yang terjadi antara permukaan padat dengan fluida yang mengalir di sekitarnya, dengan menggunakan media penghantar berupa fluida (cairan/gas) Dasar : Hukum Newton

qc = hc A⎛⎜Tw −Ts ⎞⎟ ⎝ ⎠

atau

qc A

= hc ⎛⎜ Tw − Ts ⎞⎟ ⎝



Contoh peristiwa perpindahan secara konveksi

Pergerakan udara pada peristiwa perpindahan konveksi dengan sumber panas pada salah satu sudutnya

Macam-macam Konveksi : 1. Konveksi bebas/konveksi alamiah (free convection/natural convection) perpindahan panas yang disebabkan oleh beda suhu dan beda rapat saja dan tidak ada tenaga dari luar yang mendorongnya. Contoh : plat panas dibiarkan berada di udara sekitar tanpa ada sumber gerakan dari luar

2. Konveksi paksaan (forced convection) perpindahan panas aliran gas atau cairan yang disebabkan adanya tenaga dari luar Contoh : plat panas dihembus udara dengan kipas/blower

3. RADIASI Adalah perpindahan panas yang terjadi karena pancaran/sinaran/radiasi gelombang elektromagnetik, tanpa memerlukan media perantara Dasar : Hukum Stefan-Boltzman

qr

4 = εσ AT

PERPINDAHAN PANAS KONDUKSI, KONVEKSI, RADIASI Panas yang dipancarkan dan dipantulkan

Panas radiasi dari matahari

Perpindahan panas konveksi alami dan/atau konveksi paksaan

Perpindahan panas konduksi ke tanah melalui blok beton

PERPINDAHAN PANAS KONDUKSI PERPINDAHAN PANAS KONDUKSI, STEADY STATE (TUNAK), KOORDINAT SATU DIMENSI z Meliputi : - bidang datar (x, y, z) - silinder (r, z, θ) - bola (r, θ, φ) Hukum Fourier untuk perpindahan panas konduksi :

q = −k A dT dx

Koordinat Cartesian ¾ arah x :

qx = −k A dT dx

¾ arah y:

¾ arah z :

q y = −k A dT dy

qz = −k AdT dz

Koordinat Silinder ¾ arah r :

q r = −k A dT dr

¾ arah θ:

k q = − A dT θ r dθ

¾ arah z :

q z = −k A dT dz

Koordinat Bola ¾ arah r :

q r = −k A dT dr

¾ arah θ:

k q = − A dT θ r dθ

¾ arah φ :

k q =− A dT φ r sin θ dφ

Konduktivitas Thermal (Daya Hantar Panas)

Adalah sifat bahan yang menunjukkan seberapa cepat bahan itu dapat menghantarkan panas konduksi Pada umumnya nilai k dianggap tetap, namun sebenarnya nilai k dipengaruhi oleh suhu (T). Konduktor → bahan yang mempunyai konduktivitas yang baik Contoh : logam Isolator → bahan yang mempunyai konduktivitas yang jelek Contoh : asbes

PERPINDAHAN PANAS KONDUKSI PADA BIDANG DATAR 1. Perpindahan Panas Konduksi Pada Satu Bidang Datar (Slab) q profil suhu

∆T

q

∆x

∆T Hk. Fourier : q = −k A dT = −kA dx ∆x

q=−

∆T ∆x kA

Laju perpindahan panas, q → aliran Temperatur → potensial konduktivitas thermal, k tebal bahan, ∆x luas permukaan, A

tahanan

Analogi listrik (Hk. Ohm) →

I= V R



Aliran = potensial tahanan

q=−

∆T ∆x kA

Bila aliran panas dinyatakan dengan analogi listrik menjadi : →q

T1

T2 R

⎛T − T ⎞ ⎜ ⎟ ∆T 1 ⎝ 2 ⎠ q=− =− ∆x R kA

∆T T1 − T2 q= = R ∆x kA

Contoh Soal : Salah satu permukaan sebuah plat tembaga yang tebalnya 3 cm mempunyai suhu tetap 400oC, sedangkan suhu permukaan yang sebelah lagi dijaga tetap 100oC. Berapa panas yang berpindah melintas lempeng itu?

2. Perpindahan Panas Konduksi Pada Satu Seri Bahan z Aliran panas dilewatkan pada bidang datar yang disusun berlapis-lapis dengan bahan yang berbeda-beda. z Aliran panas masuk dengan suhu T1 dan keluar dengan suhu T4. Suhu antar muka masing-masingnya adalah T2 dan T3. z Contoh : pada konstruksi furnace, boiler, dll.

A

B

C

T1 T2 q

q

kA T3 kB kC

∆xA

∆xB

T4

∆xC

Analogi listrik bahan yang disusun secara seri : q T1

T2

RA

T4

T3

RB

RC

Persamaan aliran panas untuk seluruh bidang datar adalah :

∆T q = menyeluruh ∑R th Rth adalah jumlah tahanan thermal. Untuk bahan yang disusun seri : Rth = RA + RB + RC + … Persamaan aliran panas untuk bidang yang disusun seri adalah :

∆T ∆T menyeluruh q= = RA +RB +RC ∑R th

T1 − T4 q= ∆x A ∆ x B ∆x C + + k A A k BA k CA

Pada keadaan steady state, panas yang masuk pada sisi muka sebelah kiri harus sama dengan panas yang meninggalkan sisi muka sebelah kanan, qinput = qoutput sehingga,

q = qA = qB = qC ∆T ∆TA ∆TB ∆TC q= = = = RA RB RC ∑R th

T1 − T2 qA = ∆x A k AA

T2 − T3 qB = ∆x B k BA

T3 − T4 qC = ∆x C k CA

Contoh Soal: Dinding furnace dilapisi oleh 3 lapisan : firebrick dengan ketebalan 6 in (k=0.95 Btu/h.ft.oF), insulating brick (k=0.4 Btu/h.ft.oF) dan common brick (k=0.8 Btu/h.ft.oF). Suhu masuk firebrick, T1 = 1800oF, suhu maksimum insulating brick, T2 = 1720oF dan suhu T3 = 280oF . z Hitunglah ketebalan lapisan insulating brick ! z Jika common brick tebalnya 9 in, hitunglah suhu keluar !

3. Perpindahan Panas Konduksi Melalui Bahan yang Disusun Seri dan Paralel Dinding yang terdiri atas beberapa macam bahan yang dihubungkan seri dan paralel dialiri panas. Perpindahan panas konduksi dianggap berlangsung hanya satu arah (arah x). T0

T1

T2

T3

T4 4a

2a 4b q

1 2b

∆x1

q

3

∆x2

4c

∆x3

∆x4

Analogi listrik untuk susunan seri dan paralel : R2a T0

Rk2

Rk1 T2

T1 R1

T3

R4a R4b

T4

R3 R2b

R4c

Untuk menyelesaikan susunan di atas, maka tahanan yang disusun paralel harus diselesaikan lebih dahulu sehingga pada akhirnya akan terbentuk susunan seri.

1 1 1 1 Untuk susunan paralel : R = R + R + R + ..... 1 2 3 Persamaan aliran panas untuk susunan di atas adalah :

q=

∆T ∆T = R1 + R k1 + R 3 + R k 2 ∑R th

∆x1 R1 = k1A1 ∆x 3 R3 = k 3A 3

∆x 2 R k1 = k 2a A 2a + k 2b A 2b ∆x 4 R k2 = k 4a A 4a + k 4b A 4b + k 4c A 4c

Penyelesaian persamaan aliran panas untuk susunan seri dan paralel adalah :

T0 − T4 q= ∆x 3 ∆x1 ∆x 2 ∆x 4 + + + k1A1 k 2a A 2a + k 2b A 2b k 3A3 k 4a A 4a + k 4b A 4b + k 4c A 4c

PERPINDAHAN PANAS KONDUKSI PADA SILINDER 1. Perpindahan Panas Konduksi pada Silinder Berongga Suatu silinder panjang berongga dengan jari-jari dalam ri, jari-jari luar ro dan panjang L dialiri panas sebesar q. Suhu permukaan dalam Ti dan suhu permukaan luar To. L To

ro ri Ti

Analogi listrik :

→q Ti

To

R

Aliran panas hanya berlangsung ke arah radial (arah r) saja. Luas bidang aliran panas dalam system silinder ini adalah : Ar = 2πrL Sehingga hukum Fourier menjadi : ⎛



q = kAr ⎜⎜ − dT ⎟⎟ = −k 2πrL dT dr ⎝ dr ⎠ Kondisi batas (Boundary Condition, BC) : T = Ti (i) r = ri (ii) r = ro T = To Dengan kondisi batas di atas, persamaan aliran panas untuk koordinat silinder adalah :

2πkL⎛⎜ T − To ⎞⎟ ⎝ i ⎠ q= ln⎛⎜ ro r ⎞⎟ i⎠ ⎝

atau

2πkL⎛⎜ T − To ⎞⎟ ⎠ ⎝ i q= 2,3 log ⎛⎜ ro r ⎞⎟ i⎠ ⎝

T − To ∆T = i q= ⎛ ⎞ R ⎜⎜ r ⎟⎟ ln r th o i ⎝ ⎠ 2πkL

ln⎛⎜ ro r ⎞⎟ Dalam hal ini tahanan thermalnya adalah : R = ⎝ i⎠ th 2πkL r D Jika D adalah diameter silinder maka : o = o r D i i Persamaan aliran panas dapat ditulis, q=

2πkL⎛⎜ T − To ⎞⎟ ⎝ i ⎠ ln⎛⎜ Do D ⎞⎟ i⎠ ⎝

atau

2πkL⎛⎜ T − To ⎞⎟ ⎝ i ⎠ q= 2,3 log ⎛⎜ Do D ⎞⎟ i⎠ ⎝

Jika diameter dalam silinder (Di) > 0,75 diameter luar (Do), aliran panas bisa dicari dengan : T − To i q= ⎛ ⎞ ⎜ Do − D ⎟ 2 i⎠ ⎝

πkL⎛⎜ D + Do ⎞⎟ 2 ⎝ i ⎠

2. Perpindahan Panas Konduksi pada Dinding Lapis Rangkap Berbentuk Silinder Sebuah silinder yang suhu permukaannya relatif tinggi dapat diisolasi dengan beberapa macam bahan yang disusun seri. L

kC kB kA

T1 r1

r2

A

T2 r3

T3

B r4

C

T4

q

Analogi listrik :

T1

T2

RA

T4

T3

RB

RC

Persamaan aliran panas untuk dinding lapis rangkap berbentuk silinder adalah :

∆T ∆T menyeluruh q= = RA +RB +RC ∑R th

RA =

ln(r3 r2 ) RB = 2πk BL

ln(r2 r1) 2πk AL

RC =

ln(r4 r3 ) 2πk CL

sehingga,

q=

(

ln r2 r1

T1 − T4 ln r3 r2

)+ (

2πk A L

2πk B L

) + ln(r4 r3 ) 2πk C L

atau

q=

(

ln r2 r1 kA

2πL⎛⎜ T1 − T4 ⎞⎟ ⎝



) + ln(r3 r2 ) + ln(r4 r3 ) kB

kC

qinput = qoutput sehingga,

∆T ∆TA ∆TB ∆TC q= = = = RA RB RC ∑R th

q=

T1 − T4 ∑R

th

T3 − T4 T2 − T3 T1 − T2 = = = ln r2 r1 ln r3 r2 ln r4 r3

(

)

2πk A L

(

2πk B L

)

(

2πk C L

)

Contoh soal : Sebuah pipa uap panas mempunyai suhu dalam 250oC. Diameter dalam pipa adalah 8 cm, tebalnya 5,5 mm. Pipa itu dilapisi dengan lapisan isolasi yang mempunyak k = 0,5 W/m.oC setebal 9 cm, diikuti dengan lapisan lain dengan k = 0,25 W/m.oC setebal 4 cm. Suhu luar isolasi adalah 20oC. Hitunglah kehilangan kalor per satuan panjang andaikan k = 47 W/m.oC untuk pipa !

PERPINDAHAN PANAS KONDUKSI PADA BOLA 1. Perpindahan Panas Konduksi pada Bola Berongga Suatu bola berongga dengan jari-jari dinding dalam ri, jari-jari dinding luar ro dan panjang L dialiri panas sebesar q. Suhu permukaan dalam Ti dan suhu permukaan luar To. To ro

ri

Ti

→q

Analogi listrik :

Ti

To

R

Aliran panas hanya berlangsung ke arah radial (arah r) saja. Luas bidang aliran panas adalah : Ar = 4πr2 Sehingga hukum Fourier menjadi : ⎛ dT ⎞ q = kAr ⎜⎜ − ⎟⎟ = −k 4πr 2 dT dr ⎝ dr ⎠

Kondisi batas (Boundary Condition, BC) : T = Ti (i) r = ri (ii) r = ro T = To Dengan kondisi batas di atas, persamaan aliran panas untuk koordinat bola adalah : T − To 4πk⎛⎜ T − To ⎞⎟ ∆ T ⎝ i ⎠ = i q= q= 1 −1 1− 1 R th r ro r ro i i 4πk Dalam hal ini tahanan thermalnya adalah :

1 −1 ro − r r ro i i R = = th 4πk 4πk r ro i

2. Perpindahan Panas Konduksi pada Dinding Lapis Rangkap Berbentuk Bola T4 r4

T3

Sebuah bola yang suhu permukaannya relatif tinggi dapat diisolasi dengan beberapa macam bahan.

r3 r2 r1

T2

T1

k1

Analogi listrik :

k2 k3

q T1

T2

R1

T4

T3

R2

R3

Persamaan aliran panas untuk dinding lapis rangkap berbentuk bola adalah : ∆T ∆T menyeluruh q= = R1 + R 2 + R 3 ∑R th sehingga,

q=

1 −1 r1 r2 4πk1

+

T1 − T4 1 −1 r2 r3 4πk 2

+

1 −1 r3 r4

atau

q=

4πk 3

4π⎛⎜ T1 − T4 ⎞⎟ ⎠



1 −1 r1 r2 k1

+

1 −1 r2 r3 k2

qinput = qoutput

∆T ∆T1 ∆T2 ∆T3 q= = = = R1 R 2 R 3 ∑R th q=

T1 − T4 ∑R

th

T3 − T4 T2 − T3 T1 − T2 = = = 1 −1 1 −1 1 − 1 r3 r4 r1 r2 r2 r3 4πk1

4πk 2

4πk 3

+

1 −1 r3 r4 k3

Contoh Soal :

Sebuah bola lowong terbuat dari alumunium (k = 202 W/m.oC) dengan diameter dalam 4 cm dan diameter luar 8 cm. Suhu bagian dalam adalah 100oC dan suhu luar 50oC. Hitunglah perpindahan kalornya !

PERPINDAHAN PANAS KONDUKSI DAN KONVEKSI SECARA SIMULTAN

KOEFISIEN PERPINDAHAN PANAS MENYELURUH (OVERALL HEAT TRANSFER COEFFICIENT, U)

Adalah merupakan aliran panas menyeluruh sebagai hasil gabungan proses konduksi dan konveksi. Koefisien perpindahan panas menyeluruh dinyatakan dengan W/m2.oC (Btu/h.ft2.oF)

1. KOEFISIEN PERPINDAHAN PANAS MENYELURUH PADA BIDANG BATAR Suatu bidang datar, salah satu sisinya terdapat fluida panas A dan sisi lainnya terdapat fluida B yang lebih dingin. TA

T1

Fluida A

Fluida B k

h2

q h1

T2

TB

q

Analogi listrik :

TA

T1

RA

TB

T2

R12

RB

Perpindahan panas menyeluruh dinyatakan dengan :

A⎛⎜ TA − TB ⎞⎟ TA − TB ⎝ ⎠ q= = 1 1 + ∆x + 1 + ∆x + 1 h1A kA h 2A h1 k h2 Selain itu

q = UA ∆Tmenyeluruh

sehingga koefisien perpindahan panas menyeluruh dapat dinyatakan dengan :

1 U= 1 + ∆x + 1 h1 k h2

Untuk bidang datar yang disusun seri,

q=

TA − TB

1

h 1A

+ ∑ ⎛⎜ ∆x ⎝

⎞ + kA ⎟⎠

1

= h 2A

A⎛⎜ TA − TB ⎞⎟ ⎝

1

h1

+ ∑ ⎛⎜ ∆x ⎝



⎞ + k ⎟⎠

1

h2

sehingga koefisien perpindahan panas menyeluruh dapat dinyatakan dengan :

U=

U=

1

h1

1 + ∑ ⎛⎜ ∆x ⎝

⎞ + k ⎟⎠

1

h2

1 ⎛ A⎜⎜ R C 1 ⎝

⎞ + ∑ R k +R C ⎟⎟ 2 ⎠

2. KOEFISIEN PERPINDAHAN PANAS MENYELURUH PADA SILINDER Suatu silinder berongga terkena lingkungan konveksi di permukaan bagian dalam dan luar oleh fluida A dan fluida B. Suhu kedua fluida, TA dan TB. Zat alir mengalir melalui pipa pada suhu TA. Perpindahan panas dari zat alir ke pipa secara konveksi diteruskan lewat pipa secara konduksi dan selanjutnya ke zat alir yang ada di luar pipa pada suhu TB secara konveksi. L

r1

r2

Analogi listrik :

TA

q

T1

TA

T1

TB

T2

T2

T

TB r

RC1

Rk

RC2

Perpindahan panas menyeluruh dari zat alir di dalam pipa ke zat alir di luar pipa adalah

q=

TA − TB ln ⎛⎜ r2 r1 ⎞⎟ 1 1 ⎠ + ⎝ + h1A1 2πkL h 2A 2

Luas permukaan untuk perpindahan panas zat alir : ‰ di dalam pipa, A1 = 2πr1L ‰ di luar pipa, A2 = 2πr2L sehingga,

q=

TA − TB ln ⎛⎜ r2 r1 ⎞⎟ 1 1 ⎠ + ⎝ + h1 2π r1L 2πkL h 2 2π r2 L

=

2πL⎛⎜ TA − TB ⎞⎟ 1 + h1r1

⎝ ln ⎛⎜ r2 ⎝

k



r1 ⎞⎟ ⎠

+

1 h 2 r2

Koefisien perpindahan panas menyeluruh dapat didasarkan atas bidang dalam atau bidang luar tabung. ‰ Bidang dalam,

A1 (TA − TB )

q=

=

⎛ ⎞ A1 1 A1 ln ⎜⎝ r2 r1 ⎟⎠ + + h1 2πkL h 2A 2 1 U1 = ⎛r r ⎞ r ln r1 1 1 ⎜⎝ 2 1 ⎟⎠ + + h1 k h 2 r2 ‰ Bidang luar, A 2 (TA − TB ) q= = ⎛ ⎞ A 2 A 2 ln ⎜⎝ r2 r1 ⎟⎠ 1 + + h1A1 2πkL h2

U2 =

r2 h1r1

+

1 r2 ln ⎛⎜ r2 r1 ⎞⎟ ⎝

k



+

2πr1L⎛⎜ TA − TB ⎞⎟

1 h2

1 + h1

⎝ r1 ln ⎛⎜ r2 ⎝

k



r1 ⎞⎟ ⎠

+

r1 h 2 r2

2π r2 L⎛⎜ TA − TB ⎞⎟ ⎝

r2 h1r1

+



r2 ln ⎛⎜ r2 r1 ⎞⎟ ⎝

k



+

1 h2

3. KOEFISIEN PERPINDAHAN PANAS MENYELURUH PADA BOLA Analogi listrik : T1 r1 TA

r2

q T2

TA

TB

T1

RA

TB

T2

R12

RB

Perpindahan panas menyeluruh dari zat alir di dalam pipa ke zat alir di luar pipa adalah

q= 1 + h1A1

TA − TB 1 −1 r1 r2 4πk

+

1 h 2A2

Koefisien perpindahan panas menyeluruh, ‰ Bidang dalam, 4π r12 ⎛⎜ TA − TB ⎞⎟ A1 (TA − TB ) ⎝ ⎠ q= = ⎛ ⎞ 2 ⎛⎜ 1 − 1 ⎞⎟ A1 ⎜⎜ 1r − 1r ⎟⎟ r 2 r ⎟ A r ⎜ r 1 1 1 2⎠ + 1 2⎠ + 1 + ⎝ 1 + ⎝ 1 h1 4πk h 2 A 2 h1 k h 2 r2 2 1 U1 = ⎛ ⎞ r12 ⎜⎜ 1r − 1r ⎟⎟ r12 1 2 ⎝ 1 ⎠ + + h1 k h 2 r2 2 ‰ Bidang luar, 4π r2 2 ⎛⎜ TA − TB ⎞⎟ A 2 (TA − TB ) ⎝ ⎠ q= = ⎛ ⎞ 2 ⎛⎜ 1 − 1 ⎞⎟ ⎜ 1 − 1 ⎟ A r 2 r ⎟ r ⎟ A2 r2 2 ⎜⎝ r1 2 ⎜⎝ r1 1 2 2⎠ + 1 ⎠ + + + h1A1 4πk h 2 h1r12 k h2

1

U2 = r2 2 h1r12



+



r2 2 ⎜⎜

1 − 1 ⎟ r r ⎟ 2⎠ + ⎝ 1

k

1 h2

Contoh soal : ¾ Sebuah bola lowong terbuat dari alumunium (k = 202

W/m.oC) dengan diameter dalam 4 cm dan diameter luar 8 cm. Suhu bagian dalam adalah 100oC dan suhu luar 50oC. Hitunglah perpindahan kalornya! ¾ Jika bola di atas dilapisi dengan bahan isolasi yang mempunyai k = 50 mW/m.oC setebal 1 cm. Bagian luar isolasi ini bersentuhan dengan lingkungan yang mempunyai h = 20 W/m2.oC dan Ts = 10oC. Bagian dalam bola tetap mempunyai suhu 100oC, hitunglah perpindahan kalor dalam kondisi ini!

TEBAL ISOLASI KRITIS 1. SILINDER TERISOLASI Sebuah pipa bundar dipasang selapis isolasi di sekelilingnya. Suhu dinding dalam isolasi adalah Ti sedang suhu luarnya terkena konveksi sebesar Ts.

h, Ts

ri Ti

rc

T

Analogi listrik untuk pipa terisolasi adalah

ln⎛⎜ rc r ⎞⎟ Rk = ⎝ i ⎠ 2πkL

q Ti

T

Rk

Ts

Rh =

Rh

1 2π rcLh

Persamaan perpindahan panas untuk pipa terisolasi adalah :

q=

∆Tmenyeluruh ∑R

=

T − Ts i

ln ⎛⎜ rc r

⎞ ⎟ i⎠



th

2πkL

q=

+

1 2π rc Lh

2πL⎛⎜ T − Ts ⎞⎟ ⎝ i

⎞ ⎟ i⎠

ln ⎛⎜ rc r ⎝

k



+

1 rc h

Untuk menentukan jari-jari kritis isolasi (rc) agar perpindahan panasnya maksimum dapat dicari dengan 2 cara, yaitu

dq = 0 drc

atau

Jari-jari kritis diperoleh :

dR = 0 drc

rc = k h

Artinya, perpindahan panas maksimum dari pipa terjadi ketika jarijari kritis sama dengan ratio konduktivitas thermal isolasi dengan koefisien perpindahan panas permukaan.

k Jika rc < h rc > k h

perpindahan panas meningkat dengan penambahan tebal isolasi. perpindahan panas menurun dengan penambahan tebal isolasi.

2. BOLA TERISOLASI Sebuah bola dipasang selapis isolasi di sekelilingnya. Suhu dinding dalam isolasi adalah Ti sedang suhu luarnya terkena konveksi sebesar Ts.

h, Ts

Analogi listrik untuk bola terisolasi adalah q

ri

rc Ti

Ti

T

Ts

T

Rk

1 −1 r rc i Rk = 4πk

Rh

Rh =

1 4π rc2h

Persamaan perpindahan panas untuk bola terisolasi adalah : q=

∆Tmenyeluruh ∑R

th

T − Ts i = 1 −1 r rc 1 i + 4πk 4π rc 2 h

4π⎛⎜ T − Ts ⎞⎟ ⎝ i ⎠ q= 1 −1 r rc 1 i + k rc 2h

Untuk menentukan jari-jari kritis isolasi (rc) agar perpindahan panasnya maksimum dapat dicari dengan 2 cara, yaitu

dq = 0 drc Jari-jari kritis diperoleh :

atau

rc = 2k h

dR = 0 drc

Contoh soal : Sebuah benda berbentuk pipa berdiameter 5 cm dan bersuhu 200oC diisolasi dengan menggunakan asbes (k = 0,17 W/m.oC). Benda tersebut terkena udara kamar yang suhunya 20oC dengan h = 3,0 W/m2.oC. ‰ Turunkan persamaan untuk jari-jari kritis isolasi tersebut ! ‰ Hitunglah jari-jari kritis isolasi asbes ! ‰ Hitung panas yang hilang pada jari-jari kritis ! ‰ Hitung panas yang hilang jika tanpa isolasi !

PERPINDAHAN PANAS KONVEKSI Cara-cara meramalkan nilai koefisien perpindahan kalor konveksi, h

KONVEKSI PAKSA (FORCED CONVECTION FLOW SYSTEM) Z ALIRAN DI ATAS PLAT RATA Daerah transisi

Daerah laminar

Daerah turbulen

U∞ U∞

U

U

Berbagai daerah aliran lapisan batas di atas plat rata

Pengelompokan aliran yang mengalir di atas plat diketahui dari bilangan Reynolds

Re =

U ∞ .x ρ.U ∞ .x = µ υ

U∞ = kecepatan aliran bebas x = jarak dari tepi depan υ = µ/ρ = viskositas kinematik Transisi dari aliran laminar menjadi turbulen terjadi bila Re > 5.105 Untuk aliran sepanjang plat rata, lapisan batas selalu turbulen untuk Re ≥ 4. 106 dimana :

Z ALIRAN DALAM TABUNG

Aliran berkembang penuh

Untuk aliran turbulen biasanya

Re d =

U m .d U m .d.ρ = > 2300 υ µ

Z LAPISAN BATAS PADA PLAT RATA Lapisan Batas Termal

Daerah dimana terdapat gradien suhu dalam aliran akibat proses pertukaran kalor antara fluida dan dinding Lapisan Batas Hidrodinamik

Daerah aliran dimana gaya-gaya viscous dirasakan T∞

δt

Tw = suhu dinding T∞ = suhu fluida di luar lapisan batas termal δt = tebal lapisan termal

Tw

qw dT = −k dy A

w

Angka Prandtl

Parameter yang menghubungkan ketebalan relatif antara lapisan batas hidrodinamik dan lapisan batas termal υ µ ρ Cp.µ = = k α k ρCp h .x Nu x = x k Pr =

Angka Nusselt :

Untuk plat yang dipanaskan pada keseluruhan panjangnya :

Nu x = 0,332 Pr Re x berlaku untuk fluida yang mempunyai angka Prandtl antara 0,6 – 50. 13

12

12 12 Untuk angka Prandtl yang rendah : Nu x = 0,530 Pr Re x

Untuk Angka Prandtl yang tinggi :

Nu x =

0,3387 Re x

12

Pr 1 3

⎡ ⎛ 0,0468 ⎞ 2 3 ⎤ ⎟ ⎥ ⎢1 + ⎜ Pr ⎠ ⎥⎦ ⎢⎣ ⎝

14

Koefisien perpindahan kalor rata-rata dan angka Nusselt bisa diperoleh dengan :

h = 2hx

Nu L = 2 Nu x = 0,664 Re L

12

Pr 1 3

dimana

ρ.U ∞ .L Re L = µ

Analisa di atas didasarkan atas pengandaian bahwa sifat-sifat fluida konstan di seluruh aliran. Jika terdapat perbedaan menyolok antara kondisi dinding dan kondisi aliran bebas, sifat-sifat tersebut dievaluasi pada suhu film, Tf yaitu rata-rata aritmatik antara suhu dinding dan suhu aliran bebas.

Tf =

Tw + T∞ 2

Beda suhu rata-rata sepanjang plat dapat dihitung dengan :

Tw − T∞ =

qw L k 0,6795 Re L

12

Pr1 3

Z ALIRAN TURBULEN DALAM TABUNG Untuk aliran turbulen yang sudah jadi atau berkembang penuh : Bilangan Reynolds :

Re d =

Bilangan Nusselt : Nu d =

ρ Um d µ

hd k

Nu d = 0,023 Re d 0,8 Pr n Nilai n

: n = 0,4 n = 0,3

untuk pemanasan untuk pendinginan

Perpindahan kalor per satuan panjang :

q = h π d (Tw − Tb ) L

Contoh Soal : Udara pada 27oC dan 1 atm mengalir di atas sebuah plat rata dengan kecepatan 2 m/s. Jika plat dipanaskan keseluruhan panjangnya hingga mencapai suhu 60oC, hitunglah panas yang dipindahkan pada (a) 20 cm pertama plat, dan (b) 40 cm pertama plat.

KONVEKSI BEBAS (NATURAL CONVECTION) Konveksi yang terjadi karena proses pemanasan yang menyebabkan fluida berubah densitasnya (kerapatannya) dan bergerak naik Gerakan fluida dalam konveksi bebas terjadi karena gaya bouyancy (apung) yang dialaminya apabila kerapatan fluida di dekat permukaan perpindahan kalor berkurang sebagai akibat proses pemanasan.

[ PLAT/SILINDER VERTIKAL g.β(Tw − T∞ )L3 Bilangan Grashoff : GrL = υ2 dimana :

g = percepatan gravitasi ϑ = viskositas kinematik β = 1/T = koefisien ekspansi volume (K-1)

Koefisien perpindahan kalor dievaluasi dari :

q w = h A (Tw − T∞ ) Koefisien perpindahan kalor konveksi bebas rata-rata untuk berbagai situasi dinyatakan dalam bentuk : hL Nu f = C (Grf Prf )m = k f menunjukkan bahwa sifat-sifat untuk gugus tak berdimensi dievaluasi pada suhu film : T + T∞ Tf = w 2

Gr.Pr = Ra (Bilangan Rayleigh) Harga C dan m dapat dilihat pada tabel : Jenis Aliran Laminar

Gr.Pr (Ra)

C

M

104 – 109 109 – 1013

0,59 0,10

¼ 1/3

Korelasi yang lebih rumit diberikan oleh Churchill dan Chu :

Nu = 0,68 +

Nu

12

0,670 Ra 1 4

[1 + (0,492 / Pr ) ]

= 0,825 +

9 16 4 9

0,387 Ra 1 6

[1 + (0,492 / Pr ) ]

9 16 8 27

untuk 10-1 < RaL < 109

untuk 10-1 < RaL < 1012

[ PLAT HORISONTAL Plat horisontal dengan permukaan panas menghadap ke atas :

Nu L = 0,13 (GrL Pr )1 3

untuk GrL.Pr < 2 x 108

Nu L = 0,16 (GrL Pr )

untuk 2 x 108 < GrL.Pr < 1011

13

Plat horisontal dengan permukaan panas menghadap ke bawah :

Nu L = 0,58 (GrL Pr )

untuk 106 < GrL.Pr < 1011

15

Jangan lupa bahwa :

Nu L =

hL k

q = h A (Tw − T∞ )

[ SILINDER HORISONTAL 14 g β (Tw − T∞ )d3 ( ) Nu = 0 , 53 Gr Pr Grd = d d 2 υ q = h π d (Tw − T∞ ) L

h=

k Nu d d

[ KONVEKSI BEBAS DARI BOLA Nilai Nusselt rata-rata untuk bola isotermal ke udara :

Nu f =

hd = 2 + 0,392 Grf 1 4 kf

untuk 1 < Grf < 105

Dengan memasukkan angka Prandtl diperoleh :

Nu f = 2 + 0,43 (Grf Prf )

14

Untuk rentang yang lebih tinggi : Nu f = 2 + 0,50 (Grf Prf )1 4

untuk 3 x 105 < Gr Pr < 8 x 108

PERPINDAHAN PANAS RADIASI

Radiasi ≅ pancaran ≅ sinaran ≅ ilian Radiasi thermal → radiasi elektromagnetik yang dipancarkan oleh suatu benda

karena suhunya. Radiasi selalu merambat dengan kecepatan cahaya, 3 x 1010 cm/s. Kecepatan ini sama dengan hasil perkalian panjang gelombang dengan frekuensi radiasi :

c = λν dimana : c = kecepatan cahaya λ = panjang gelombang ( = 10-8 cm) ν = frekuensi Perambatan radiasi thermal berlangsung dalam bentuk kuantum dan setiap kuantum mengandung energi sebesar

E = hν h = konstanta Planck, 6,625 x 10-34 J.s Setiap kuantum dianggap sebagai suatu partikel yang mempunyai energi, massa dan momentum seperti molekul gas → photon Sehingga, pd hakekatnya radiasi merupakan pancaran yg disebabkan oleh gas photon yang mengalir dari satu tempat ke tempat lain.

Dengan teori relatifitas dan thermodinamika statistik maka akan diperoleh suatu rumus yang disebut Hukum Stefan-Boltzmann dimana energi total yang dipancarkan oleh suatu benda sebanding dengan pangkat empat suhu absolut :

Eb = σ T4 Dilihat dari daya emisinya, benda terbagi ke dalam 3 macam : 1. Benda putih sempurna (absolutely white) → menyerap sinar, tanpa mengemisikan kembali. Emisivitas (ε) = 0 2. Benda abu-abu (gray body) 0<ε<1 3. Benda hitam (blackbody) → menyerap 100%, mengemisikan 100%. Emisivitas (ε) = 1

SIFAT-SIFAT RADIASI Sifat-sifat benda yang menerima energi radiasi : radiasi datang

dipantulkan/refleksi (ρ)

diserap/absorpsi (α)

diteruskan/transmisi (τ)

ρ= faktor refleksi (refleksivitas) α = faktor absorpsi (absorpsivitas) τ = faktor transmisi (transmisivitas)

ρ + α + τ =1 Kebanyakan benda padat tidak meneruskan radiasi thermal, τ = 0, sehingga

ρ + α =1

Sifat-sifat radiasi benda, 1. Benda yang sifatnya dapat menyerap energi yang datang seluruhnya (100%) disebut benda hitam (blackbody) α=1 ; ρ=0 Emisi benda hitam, ε = 1 → ε = α = 1 2. Benda yang dapat memantulkan energi yang datang 100% disebut benda putih sempurna (absolutely white) ρ=1 ; α=0 3. Benda yang diantara black body dan white body disebut benda abu-abu (grey body) 0<ε<1

IDENTITAS KIRCHHOFF Emisivitas (ε) suatu benda sama dengan absorpsivitas (α)-nya pada suhu yang sama Emisivitas suatu benda (ε) → perbandingan antara energi yang dapat dipancarkan oleh benda itu pada suhu T dibandingkan dengan ε= E E energi yang dipancarkan oleh b benda hitam pada suhu yang sama

Energi yang dipancarkan oleh suatu benda selalu lebih kecil dari energi yang dipancarkan oleh benda hitam sehingga harga ε ≤ 1.

FAKTOR PANDANGAN (Fm-n)

‰ Faktor bentuk (shape factor) ‰ Faktor pandang (view factor) ‰ Faktor sudut (angle factor) ‰ Faktor konfigurasi (configuration factor) ‰ Faktor geometris (geometry factor)

Eb1

Eb2

T1 A1

T2 A2

Pertukaran energi antara dua permukaan yang mempunyai suhu yang berlainan

Permukaan 1 dan permukaan 2 saling meradiasi energi di permukaan 1 bisa sampai di permukaan 2 dan sebaliknya. F1-2 = fraksi energi yang meninggalkan permukaan 1 dan diterima oleh permukaan 2. F2-1 = fraksi energi yang meninggalkan permukaan 2 dan diterima oleh permukaan 1 Fm-n = fraksi energi yang meninggalkan permukaan m dan diterima oleh permukaan n

Energi yang meninggalkan permukaan 1 dan sampai di permukaan 2 adalah : Eb1A1F12 Energi yang meninggalkan permukaan 2 dan sampai di permukaan 1 adalah : Eb2A2F21 Pertukaran energi nettonya adalah : q1-2 = Eb1A1F12 - Eb2A2F21 Pada 2 permukaan m dan n berlaku hubungan resiprositas AmFmn = AnFnm Sehingga pertukaran kalor nettonya menjadi : q1-2 = A1F12(Eb1-Eb2) = A2F21(Eb1-Eb2)

HUBUNGAN BERBAGAI FAKTOR BENTUK Benda-benda tidak bisa memandang dirinya sendiri : F11 = F22 = F33 = … = 0 Jika Fij adalah fraksi energi total yang meninggalkan permukaan i dan sampai di permukaan j maka : n

∑ Fij = 1

j=1

Untuk lengkung tiga permukaan dapat kita tuliskan : F11 + F12 + F13 = 1 F11 = 0

F13 = 1 – F12 F21 + F22 + F23 = 1

F22 = 0

F23 = 1 – F21

Dari hubungan resiprositas : A1F12 = A2F21

PERTUKARAN KALOR ANTARA BENDA TAK HITAM Pada perpindahan kalor radiasi antara permukaan hitam, semua energi radiasi yang menimpa permukaan itu diserap. Pada benda tak hitam, tidak seluruh energi yang jatuh di permukaan diserap; sebagian dipantulkan kembali ke permukaan lain dalam system dan sebagian mungkin dipantulkan keluar system. Diandaikan semua permukaan bersifat difus (baur, menyebar) dan mempunyai suhu seragam, emisivitas dan refleksivitas konstan di seluruh permukaan. Didefinisikan : G = iradiasi

panas radiasi total yang menimpa suatu permukaan sebuah benda per satuan waktu per satuan luas J = radiositas

panas radiasi total yang meninggalkan suatu permukaan sebuah benda per satuan waktu per satuan luas Dianggap seluruh permukaan mempunyai G dan J yang sama.

Radiositas → jumlah energi yang dipancarkan (emisi) dan energi yang

dipantulkan (refleksi) apabila tidak ada energi yang diteruskan (transmisi, τ = 0) α+ρ=1 ρ=1-α=1-ε

sehingga

J = εEb + ρG = εEb + (1 - ε)G

J − εEb G= 1− ε Energi netto yang meninggalkan permukaan adalah :

q = J −G A = εEb + (1− ε)G − G = εEb − εG

Masukkan persamaan G, akan diperoleh :

q = εA ⎛⎜⎝ E b − J ⎞⎟⎠ 1− ε Dari persamaan di atas diperoleh

− J ⎞⎟⎠ 1− ε εA

⎛⎜ E q=⎝ b



Arus =

beda potensial tahanan permukaan

Jaringan permukaan : →q Eb

J

1− ε εA

Pertukaran energi radiasi antara permukaan A1 dan A2 A1

A2 J1

J2 F12

F21

Energi yang meninggalkan permukaan 1 dan mencapai permukaan 2 adalah : J1A1F12 Energi yang meninggalkan permukaan 2 dan mencapai permukaan 1 adalah : J2A2F21 Pertukaran kalor netto antara kedua permukaan adalah q12 = J1A1F12 – J2A2F21

Dari hubungan resiprositas : A1F12 = A2F21 Sehingga : q12 = A1F12(J1 – J2) = A2F21(J1 – J2)

q=

(J1 − J 2 ) 1 A1F12

Jaringan ruang



Arus = beda potensial tahanan ruang

→q J1

J2

1 A1F12

Jaringan radiasi merupakan gabungan antara jaringan permukaan dan jaringan ruang. Kedua unsur jaringan itu merupakan pokokpokok metode jaringan radiasi (radiation network method).

PERPINDAHAN PANAS RADIASI ANTARA DUA PERMUKAAN Perpindahan panas antara dua permukaan dan tidak ada permukaan lain di lingkungannya q Eb1

J1

1− ε1 ε1A1

Eb2

J2

1 A1F12

1− ε 2 ε2A2

Pertukaran panas nettonya adalah : qnet =

Eb1 − E b2 Eb1 − E b2 = 1− ε1 1 − ε2 ∑R + 1 + ε1A1 A1F12 ε2A2

σ⎛⎜ T14 − T2 4 ⎞⎟ ⎝ ⎠ qnet = 1− ε1 1 − ε2 + 1 + ε1A1 A1F12 ε2A2

Contoh Soal : Dua buah piring sejajar berdiameter 60 cm, terpisah pada jarak 15 cm. Suhu pada permukaan bagian atas adalah 250 K dan suhu pada permukaan bagian bawah adalah 300 K. Andaikan semua permukaan hitam, berapakah laju perpindahan kalornya ?

PERPINDAHAN PANAS RADIASI ANTARA TIGA PERMUKAAN q Eb1

1− ε1 ε1A1

Eb2

J2

J1

1− ε 2 ε2A2

1 A1F12 1 A1F13

1 A 2F23 J3

1− ε 3 ε 3A 3 Eb3

Untuk menghitung perpindahan panas antara tiga benda ini dapat diselesaikan dengan menerapkan hukum arus Kirchhoff : Jumlah semua arus yang memasuki suatu node ialah nol. Node I :

Eb1 − J1 J2 − J1 J3 − J1 + + =0 1− ε1 1 1 A1F12 A1F13 ε1A1

Node II :

J1 − J 2 Eb2 − J 2 J3 − J2 + + =0 1 1 − ε2 1 A1F12 A2F23 ε A 2 2

Node III:

J1 − J3 J2 − J3 Eb3 − J3 + + =0 1 1− ε3 1 A1F13 A2F23 εA 3 3

PERPINDAHAN PANAS RADIASI ANTARA DUA BIDANG DATAR YANG DIHUBUNGKAN DENGAN BIDANG YANG TIDAK DAPAT MENGHANTARKAN PANAS TETAPI DAPAT MEMANTULKAN SEMUA PANAS YANG DITERIMA q Eb1

1− ε1 ε1A1

Eb2

J2

J1

1− ε 2 ε2A2

1 A1F12 1 A1F13

1 A 2F23

J3= Eb3

J3 tidak dihubungkan dengan tahanan permukaan radiasi karena permukaan 3 tidak bertukaran energi, sehingga J3 = Eb3 = σ T34

Contoh : Dua buah plat yang berada dalam ruangan yang besar. Karena luas ruang A3 sangat besar maka tahanan ruang

1− ε3 = 0 sehingga Eb3 = J3 ε3A3 Untuk menghitung aliran panas pada masing-masing permukaan, kita cari radiositas J1 dan J2 dengan menggunakan hukum arus Kirchhoff. Eb1 − J1 J2 − J1 J3 − J1 + + =0 Node J1 : 1− ε1 1 1 A1F12 A ⎛⎜1− F ⎞⎟ ε1A1 1⎝ 12 ⎠

Node J2 :

E −J J1 − J2 Eb2 − J2 + + b3 2 = 0 1 1 − ε2 1 A1F12 A ⎛⎜1− F ⎞⎟ εA 2 2

2⎝

21 ⎠

Eb1 − J1 = q Panas total yang dilepas plat 1 : 1 1− ε1 ε1A1 E −J Panas total yang dilepas plat 2 : q = b2 2 2 1 − ε2 ε 2A 2 Panas yang diterima dinding kamar :

q 3 = q1 + q 2 atau

q3 =

J1 − J 3 J 2 − J 3 J − E b3 J − E b3 + = 1 + 2 1 1 1 1 A 1 F13 A 2 F23 A 1 (1 − F12 ) A 2 (1 − F21 )

Contoh Soal : Dua buah plat sejajar, ukuran 0,5 x 1,0 m berjarak 0,5 m satu sama lain. Plat yang satu dipelihara pada suhu 1000oC dan yang satu lagi pada 500oC. Emisivitas plat itu masing-masing 0,2 dan 0,5. Kedua plat itu terletak di dalam sebuah ruang yang sangat besar yang dinding-dindingnya dipelihara pada suhu 27oC. Kedua plat itu saling bertukaan kalor satu sama lain. Tentukan perpindahan netto ke setiap plat dan ke ruang !