MATAKULIAH
: PERPINDAHAN PANAS II
MODEL
:
ETS-01
TUGAS/QUIZ MATERI
: KONDUKSI TRANSIEN
AWAL TUGAS
: 06 April 2017
AKHIR TUGAS
: 08 April 2017
NAMA
: Muhammad fiqih arif rochman
NBI/KLAS
: 1421404581 / B
email
:
[email protected] [email protected]
Judul : KONDUKSI TRANSIEN . PENGANTAR .
Pernahkah anda merasakan panas di sendok anda saat anda sedang mengaduk kopi yang panas yang baru saja mendidih?. Jika pernah maka coba anda teliti apa yang terjadi, kenapa panas yang harusnya ada di air kopi tersebut itu bisa sampai menuju tangan anda?. Semua itu terjadi karena adanya perpindahan panas dari air panas lewat ke sendok alumunium kita hingga hingga ke kulit kita. Perpindahan panas tersebut ada tiga jenis yaitu konduksi ,konveksi, dan radiasi, dan perpindahan yang tadi ta di terjadi adalah perpindahan panas secara konduksi dimana panas mengalir dari benda padat tanpa perpindahan molekulnya.Didalam perpindahan panas konduksi ada beberapa faktor yang mempengaruhi mulai dari bahan, suhu, lama waktu, serta bentuk. Berikut ini akan dibahas faktor bahan yang mempengaruhi konduksi, konduksi, beserta contoh soal dan pengertiannya. Sebuah perpindahan panas yang terjadi ada dua jenis yaitu dalam kondisi stedy state dan tidak. Dimana kondisi stedy state adalah kondisi dimana factor waktu diabaikan. Tapi pada kenyataannya factor waktu sangat;ah penting didalam sebuah perpindahan panas. Karena sebagai contoh sebuah batang besi yang dipanaskan selama 10 menit akan berbeda suhu, energy panas yang terjadi bila dibandingkan dengan batang besi yang di dipanaskan selama 30 menit. Banyak kegunaan factor waktu yang digunakan dalam kehidupan sehari-hari perihal perpindahan panas yaitu contohnya sebuah oven, dalam mengoperasikan sebuah s ebuah oven maka kita akan mengatur berapa lama kita akan memanggang roti atau makanan kita di dalam oven karena jika lama waktu dan suhu yang kita gunakan tidak pas maka makanan yang kita masak akan rusak. KOMPETENSI .
Setelah menyelesaikan bab ini saya memahami topic t opic yang berkaitan dengan KONDUKSI yaitu : 1. Persamaan umum konduksi sesuai factor bentuk 2. Persamaan umum konduksi transient 3. Persamaan umum konduksi transient benda semi infinit Persamaan umum…… 1. Q = K. S . 2. 3.
Δ
− = − − ,− =erf − √ .
TEORI Panas
Telahdiketahuidapatberpindahdaritempatdengantemperaturlebihtinggiketempatdengantempeaturleb ihrendah.Hokum Percampuranpanasjugaterjadikarenapanasituberpindah,sedangkanpadakalorimeter,perindahanpanas dapatterjadidalambentukpertukaranpanasdenganluarsistem.Jadi Pemberianataupenguranganpanastidaksajamengubahtemperaturataufasazatsuatubendasecaralokal,m elainkanpanasitumerambatkeataudaribagianlainbendaatautempatlain.Peristiwa Inidisebutperindahanpanas.Menurutpenyelidikan,perpindahantenagapanasdapatdibagidalambeberap agolong-ancaraperpindahan.Panas Itudapatmerambatdarisuatubagiankebagianlainmelaluizatataubendayangdiam.Panas Jugadapatdibawaolehpartikel-partikelzatyangmengalir.Pada Radiasipanas,tenagapanasberpindahmelaluipancaranyangmerupakanjugasatucaraperindahanpanas. Setelah ini kita bahas konduksi dan selengkapnya Faktor Bentuk
Faktor bentuk konduksi didefinisikan sebagai besaran yang diagunakan untuk mengoreksi perpindahan kalor konduksi dengan bentuk geometri tertentu. Perhitungan faktor bentuk akan menghasilkan nilai (S) yang selanjutnya dimasukkan untuk menghitung perpindahan panas konduksi dalam keadaan stedi state dimana faktor waktu diabaikan sehingga rumusnya ialah
Q = K . S . (T 2 – T1) Keterangan: Q = Laju panas (Watt) K = Konduktivitas termal bahan (W/m 0K)
S = Nilai Faktor bentuk (m) T = Suhu (0K)
Konduksi Transient Analisi temperatur suatu titik yang berubah setiap waktu adalah tanggung jawab transient analysis
Pada bab-bab sebelumnya kita hanya fokus pada pembahasan dimana waktu tidak mempengaruhi temperatur, artinya temperatur hanya fungsi dari posisi, biasa dituliskan T ( x ,y ,z ) . Kasus-kasus seperti inilah yang disebut steady. Jika faktor waktu diperhitungkan sehingga temperatur merupakan fungsi dari posisi dan waktu, T ( x ,y ,z ,t ) maka akan disebut transient. Misalnya, temperatur sebuah logam yang baru dikeluarkan dari cetakan dan diletakkan di luar (udara lingkungan), maka temperaturnya akan bervariasi tergantung setelah berapa lama logam tersebut diletakkan. Jika dibandingkan dengan kondisi steady, analysis transient memerlukan usaha ekstra. Metode-metode menganalisis perpindahan panas konduksi transient akan menjadi tanggung jawab bab ini. Pada bagian awal akan diperkenalkan kenapa kita membutuhkan analyisis transient. Kemudian teknik analysis yang biasa digunakan seperti analysis seragam ( lumped analysis system ) dan metode grafik akan dibahas. Selamat membaca.
I.
Analisis Seragam
Ada kalanya perubahan temperatur dalam suatu materi dapat diasumsikan seragam. Misalnya sebuah bola besi yang berukuran kecil ( ball bearing ) awalnya mempunyai temperatur 200 oC saat dikeluarkan dari cetakan dan diletakkan di suatu ruangan. Setelah beberapa menit temperaturnya akan turun, misalnya menjadi 100 oC. Dapat dipastikan temperatur didalam besi tidak semuanya 100 oC, dipermukaan bisa saja 98 oC, di bawah permukaan bisa saja 100oC, dan di dalam intinya malah masih 102 oC. Dalam analysis temperatur bola ini bisa dianggap seragam dan diwakili oleh satu angka, yaitu 100 oC. Sebaran temperatur ini sangat tergantung pada konduktivitas bahan bola tersebut. Jika konduktivitasnya sangat besar, maka temperatur di permukaan bola akan segera sama dengan temperatur di permukaan. Sebaliknya jika konduktivitas nya rendah sebarannya akan semakin lebar. Bayangkan jika bola tersebut terbuat dari batu bata. Temperatur di permukaan batu bisa saja sudah 100 oC, tetapi temperatur di dalam intinya masih bertahan 200 oC. Analysis yang membuat asumsi bahwa temperatur di dalam benda yang dibahas adalah seragam disebut lumpedsystem analysis (saya menerjemahkannya sebagai analysisseragam). Sejujurnya tidak ada benda yang bersifat seperti ini, karena untuk bisa segera seragam maka dibutuhkan tahanan termal nol atau konduktivitas tak terhingga. Tetapi perbedaan temperatur yang dianggap kecil bisa diasumsikan seragam. Asumsi ini hampir mirip dengan analisis transient. Untuk menghilangkan variabel waktu dari analysis, maka diasumsikan steady. Di sini untuk menghilangkan variabel posisi ( x, y, z ) dari persamaan, maka asumsi seragam dipergunakan. Apakah hal ini tidak menyalahi aturan. Jawabnya selama batasan masalah dipenuhi maka nilai pendekatan bisa dipergunakan. Misalnya untuk benda-benda yang kecil dan mempunyai Sustainable Energy Research Group
konduktivitas tinggi. Koefisien konveksi di permukaan juga turut berpengaruh. Sejauh mana asumsi seragam ini dapat diterima akan dibahas lagi nantinya. Misalkan sebuah benda yang mempunyai temperatur T diletakkan di ruangan yang temperaturnya T , untuk memudahkan analisis anggap T T , maka temperatur benda akan turun seiring dengan bertambahnya waktu atau benda mengalami pendinginan. Temperatur benda, sifat fisik dan konveksi dengan lingkungan ditampilkan pada Gambar 1. Pada gambar As bukanlah luas penampang, tetapi luas permuakaan benda yang berbatasan dengan udara lingkungan.
As
T
Q
h
Gambar 1 Proses pendinginan sebuah benda di udara lingkangan Energi panas akan berpindah dari bagian luar benda yang bersentuhan dengan udara luar secara konveksi dengan koefisien h . Pada analisis ini karena koefisien konduksi benda tinggi dianggap temperatur benda akan turun dari keadaan awal sampai temperatur temperatur tertentu yang besarnya seragam untuk seluruh benda. Jika asumsi ini dipenuhi maka dapat disebutkan bahwa:
47
Mechanical Engineering, USU
Perpindahan Panas I
Prepared by: Himsar AMBARITA
pada satuan waktu tertentu, panas yang berpindah dari permukaan benda ke lingkungan akan sama dengan
penurunan energi di dalam benda. Fakta ini dapat dirumuskan dengan persamaan berikut: T − hA s(T − T ∞) = mc (1) p t Tanda minus menyatakan benda mengalami penurunan energi. Jika limit perbedaan waktunya mendekati dibuat mendekati nol, maka persamaan (1) dapat ditulis menjadi: −hA s(T − T ∞)dt = mc (2) pdT Untuk memudahkan analisis diajukan parameter berikut: θ =T −T ∞, akan membuat d θ =dT (3) Substitusi persamaan (3) ke dalam persamaan (2) akan memberikan: (4) − hAs dt = d θ mc p θ t i = 0 sampai =
t
−
T i hAs
−
mc p
0
t
waktu
dimana
temperatur
T ∞dan temperatur saat ituθ = T (t ) − T ∞. θ d θ
dt =
∫
θ i
Pada gambar ditampilkan grafik penurunan suhu benda yang ditempatkan pada suhu lingkungan. Benda pertama mempuyai massa yang lebih kecil dari benda kedua dan
dapat diintegralkan untuk waktu dari
Persamaan (4) θi
Gambar 2 Grafik penurunan suhu
awal
ketiga, atau m1
(5)
1
Jika semua parameter beriku h, As, m dan c p adalah konstant atau tidak berubah meskipun waktu dan
maka
hasil pengintegralan
hAs
θ − t = ln θi mc p Atau dapat dimodifikasi menjadi: T−T T − T
θ
=
θ i
∞
= exp −
∞
i
proses pendinginaninidapatdihitung menggunakan persamaan: t
Q =∫ qdt =
(6)
0
ρVc p =
hAs
hAs
=e −t τ
(10)
0
Q = −hAs(T i− T ∞)∫ e
mc
−tτ
t
(11)
0
kurung
Jika diintegralkan akan didapat persamaan berikut: Q = mc p(T i − T ∞)1 − e−tτ
(9) T i− T ∞ Untuk menjelaskan persamaan ini akan dibuat grafik perubahan suhu terhadap waktu dari tiga buah benda yang mempunyai τ yang berbeda.
(12)
I.1. Keabsahan Analisis Seragam
(8)
Dengan menggunakan konstanta ini, maka persamaan (7) dapat ditulis menjadi: T − T ∞
∫ − hAs(T − T ∞)dt
t p
mc p
t
s
kurawal) bias disederhakan lagi dengan mengajukan satu parameter yang biasa disebut konstanta waktu ( time constant ): =
dengan
Gunakan persamaan (9) untuk mengganti (T −T ∞) pada persamaan ini.
hA
Bagian yang didalam kurung (masih didalam
τ
3
konstanta waktu maka semakin lama waktu yang dibutuhkan untuk mendinginkannya. Sementara total energi yang dipindahkan selama
θ
temperatur benda berubah, persamaan ( 5) adalah:
2
Perhatikan Gambar 3 berikut ini, yang menunjukkan proses perpindahan panas konduksi di dalam sebuah benda dan perpindahan panas konveksi kelingkungan.
Gambar 3 Perpindahan panas dari solid ke lingkungan Pada bagian (a), perpindahan panas konveksi ke lingkungan lebih kecil daripada konduksi di dalam benda. Panas dari dalam benda akan tertahan dipermukaan, akibatnya perbedaan temperatur di dalam benda padat akan kecil. Analisis temperatur seragam akan dapat diterima di Sustainable Energy Research Group
48
Mechanical Engineering, USU
sini. Pada bagian (b)perpindahan panas konveksi di permukaan jauh lebih besar daripada konduksi di dalam benda, akibatnya akan terjadi gradient tempertur yang besar. Maka analisis seragam tidak berlaku untuk kasus ini. Parameter yang digunakan untuk membandingkan besar konduksi dan konveksi pada kasus ini disebut bilangan
3. Hitung waktu pendin ginan dengan menggunakan persamaan (9): T − T T − T t ∞ ∞= − =e −tτ ⇒ ln − ∞ T i − T ∞ i τ 50 − 30 t = −3180 ×ln
Biot (Biot number ), yaitu perbandingan tahanan kondukdi dan tahanan konveksi: (Lc kA ) hLc cond (13) Bi = = = conv k (1 hA)
Bi ≤0,1
ini
diperlukan untuk
menentukan panjang dan kecepatan konveyor yang dibutuhkan. Jika masih terlalu lama, m aka anda bias mengusulkan kecepatan angin yang sesuai untuk mencapai tujuan yang diinginkan.
Lc = V ol Asdisebut panjang karakteristik. Jika
Dimana
= −3180(−2,6) = 8008s
300 − 30 Atau sekitar 2 jam 14 menit Catatan: Analisis waktu
(14)
bilangan Biot kecil maka konduksi akan lebih kuat dari konveksi. Tetapi sebaliknya jika bilangan Biot besar konduksi lebih lemah dari konveksi. Analisis seragam disarankan untuk digunakan jika memenuhi persamaan berikut: Sebelum menggunakan persamaan (9) untuk menghitung temperature sebuah benda, maka persamaan (14) harus dipastikan dulu terpenuhi.
Contoh 1 Sebagai seorang konsultan, anda diminta oleh sebuah industri yang memproduksi bola-bola besi bantalan luncur untuk menentukan lama pendinginan yang dibutuhkan bola tersebut setelah keluar dari cetakan. Bola yang diproduksi adalah terbuat dari stainless steel dengan diameter 1cm, 3
dengan sifat fisik sebagai berikut ρ= 7854 kg/m , c p=434
II. Persamaan Umum Analisis Seragam Bagian pertama dari bab ini, yang menjelaskan analisis seragam adalah merupakan fase pemanasan dari masalah yang sebenarnya di bidang engineering. Alasannya adalah: pada bagian tersebut hanya pengaruh konveksi yang diperhitungkan sebagai “pelarian” panas dari dalam benda padat. Maka dengan hanya konveksi yang diperhitungkan terbentuklah persamaan (9). Pada kasus yang sebenarnya ada tiga faktor yang mungkin mengubah persamaan (9), yaitu: (1) adanya sumber panas di dalam benda, (2) aliran flux panas melalui permukaan benda dan (3) efek radiasi yang menjadi pelarian panas dari permukaan benda ke lingkungannya. Jika semua faktor ini diperhitungakan, maka mekanisme yang terjadi pada Gambar 1 akan berubah menjadi seperti pada Gambar 4 berikut ini:
J/kgK , k =60 W/mK . Temperatur bola saat dikeluarkan dari cetakan adalah 300 oC. Begitu bola keluar proses pendinginan dimulai dengan menghembuskan udara luar 30oC dengan kecepatan 5m/s. Berdasarkan buku perpindahan panas konveksi kecepatan ini akan menghasilkan koefisien konveksi 17,9 W/m 2K. Karena industri tersebut ingin membeli konveyor yang sesuai, maka anda diminta menghitung waktu yang dibutuhkan bola-bola tersebut mencapai temperatur 50 oC.
q′′
3
p
k =60 W/mK , h =17,9 W/m2K . T
C ,T = 30 o C ,
dan D= 1× 10−2 m.
300o
=
=
434 J/kgK , ∞
i
Ditanya: t = ? untuk mencapai temperatur akhir T = 50o C Langkah yang dilakukan adalah: 1. Hitung Lc untuk menghitung bilangan Bi 4
πr 3 1 L = Vol = 3 = r =1,67×10−2m c 4πr As 3 17,9 × 1,67 × 10 −2 hL =
2.
c
=
−
= 4,96 × 10
60 k Karena lebih kecil dari 0,1 maka analisis seragam bias digunakan. Hitung konstanta waktu dengan persamaan (8). Persamaan ini bias diubah dengan menggunakan Lc .
ρVc ρL c p = c p =3180 τ =
T
rad
sur
Eg
h
qconv
T
Gambar 4 Aliran panas pada sebuah benda
Penyelesaian: Asumsi yang digunakan sifat fisik konstan. Data yang diberikan soal adalah: ρ= 7854 kg/m , c
q
A s
Selisih energi yang masuk dan yang dibangkitkan dengan energi yang keluar dari benda tersebut akan digunakan untuk menaikkan temperatur benda tersebut. Kalimat ini dapat dituliskan dalam bentuk persamaan berikut: dT q′′ A f + E g − q conv − q rad = ρVc p 15 dt Jika perpindahan panas konveksi dan radiasi dijabarkan dengan menggunakan rumus masing-masing maka persamaan ini menjadi: 4 4 dT (16) − T ) = ρVc − hA (T − T ) − εσA (T q′′ A + E p
dt Perhatikan bahwa luas permukaan untuk flux dan luas permukaan utuk radiasi dan konveksi dibedakan. Persamaan (16) dikategorikan sebagai persamaan differensial biasa orde pertama yang bersifat non-linier dan tidak homogen. Masalahnya sekarang adalah persamaan ini f
g
s
∞
s
s
shAs h Sustainable Energy Research Group
49
Mechanical Engineering, US
Profil temperatur sama
belum bisa diselesaikan secara analitik. Untuk saat ini hanya ada satu cara menyelesaiakannya yaitu dengan cara numerik yang akan dijelaskan di Bab 7. Seandainya dilakukan penghilangan beberapa komponen dari persamaan (16), maka persamaan ini akan mungkin diselesaiakan secara analitik. Misalnya tidak ada flux, tidak ada panas yang dibangkitkan, dan tidak ada konveksi (atau sangat kecil dibandingkan radiasi), maka persamaan (16) berubah menjadi: dT 4 4 − εσA(T − T ) = ρVc ( ) dt sur s p Jika persamaan ini diintegralakan dan dilakukan modifikasi akan menjadi: =
ρVc p
T + T
4εAσT
4
n
sur
T
− 2 tan −
T sur − T
T ( x , t ) T ( x , t )
T ( x , t )
T ( x , t )
y
−
z x
sur
ssur
mana
= n
T i + tan − T i
sur +
sur
−1
T i sur
Penyederhanaan lain yang mungkin
dari persamaan
Gambar 5 Distribusi temperatur transient 1 dimensi
temperatur benda yang dianalisis cukup rendah. Maka persamaan (16) dapat disederhanakan menjadi:
dθ + aθ − b = 0 dt θ = T − T ∞ Dimana
(19) ,
′′ b =(q A f
a =hAs ρVc p . Persamaan (19)
+ E g)
ρVc p
,
dan
ini dapat diselesaikan
dengan teknik transformasi, hasilnya adalah: T − T ∞ = exp(−at ) + b a 1 − exp(−at )
x L
x =0
(16) adalah jika perpindahan panas radiasi yang diabaikan. Kasus penyederhanaan seperti bisa dilakukan jika
(20)
T i − T ∞ T i− T ∞ Perhatikan, jika seandainya tidak ada panas yang dibangkitkan di dalam benda dan juga flux dianggap nol, maka parameter b akan sama dengan 0. Maka persamaan (20) akan kembali menjadi persamaan (7). Contoh 2 Jika pada soal yang ada pada contoh 1 efek flux panas akibat persentuhan permukaan bola dengan lantai tidak diabaikan dan dianggap 1% permukaan bola melepas flux 1W/m2 ke lantai maka tentukalah waktu yang dibutuhkan bola sekarang untuk mencapai suhu 50 oC.
III. Persamaan dengan pengaruh ruang (Spatial effects)
Pada pembahasan masalah transient di bagian sebelumnya kita telah mengeliminasi pengaruh ruang dengan memanfaatkan bilangan Biot. Adakalanya kasus pengaruh ruang tidak dapat diabaikan lagi. Pada bagian ini kita akan memasukkan pengaruh ruang ke dalam analisis. Seperti yang sudah diketahui pengaruh ruang mempunyai tiga dimensi, dalam kordinat kartesian disimbolkan x , y , dan z . Tetapi hal ini tidak dilakukan sekaligus, tetapi hanya memasukkan satu dimensi saja, misalnya hanya x . Kasus seperti ini mungkin saja terjadi, misalnya untuk kasus plat yang sangat panjang dan sangat lebar, seperti yang ditampilkan pada Gambar 5.
Pada gambar plat dianggap sangat panjang (arah- y ), akibatnya pada waktu kapanpun profil temperaturnya akan diasumsikan sama untuk setiap nilai y . Demikian juga jika plat sangat lebar sehingga profil temperatur akan sama pada nilai z berapapun. Oleh karena itu, temperatur benda di setiap posisi pada waktu tertentu hanya merupakan fungsi dari x dan waktu t atau T ( x , t ) . Dengan asumsi ini maka persamaan energi (yang diturunkan pada Bab II) dapat disederhakan menjadi: 1 ∂T ∂ T = (21) 2 ∂ x α ∂t Nilai awal dan nilai batas (pada pertengahan konduksi nol dan pada permukaan konduksi sama dengan konveksi ke lingkungan) untuk persamaan ini dapat dirumuskan:
T ( x ,0)= T i ∂T
= 0 dan
∂ x
(22) −k ∂T
=
hT (L,t ) − T
∂ x x = L
x =0
∞
(23)
Bilangan Tanpa Dimensi Untuk mempermudah langkah-langkah analysis, biasanya persamaan differensial diubah ke dalam bentuk tanpa dimensi. Pada tulisan ini bilangan tanpa dimensi yang digunakan adalah: T − T ∞ x * , dan t = , x *= L T − T
θ *=
i
∞
αt L
= Fo
(24)
Semua bilangan ini tidak lagi memiliki dimensi tetapi hanya konstanta. Khusus bagian terakhir waktu tanpa dimensi ini disebut bilangan Fourier. Jika persamaan (24) didifferensialkan dan disubstitusi ke persamaan (21), akan didapat persamaan: ∂
*
2 *
θ
∂ x
∂θ 2
*
=
∂Fo
(25)
Nilai awal dan kondisi batas pada persamaan (22) dan persamaan (23) menjadi:
Sustainable Energy Research Group
50
Mechanical Engineering, USU
Perpindahan Panas I
Prepared by: Himsar AMBARITA
* * θ ( x ,0)=1 * ∂θ *
(26)
= 0dan
∂ x
*
x =0
* ∂θ * ∂ x
= − Bi θ* (1,t *) *
x
(27)
Penyelesaian eksak dengan menyelesaiakan persamaan (25) untuk bentuk-bentuk geometri yang sederhana dapat dijumpai pada buku-buku perpindahan panas. Pada bagian ini akan ditulis kembali penyelesaian utuk bidang seperti yang ditampilkan pada Gambar 6. Sebuah plat dengan temperatur awal T i ditempatkan di daerah yang mengalami konveksi h dengan temperatur lingkungan T ∞ . Untuk kasus ini penyelesian eksak persamaan (25) adalah: ∞
(28)
n=1
Dimana koefisien C n dihitung dengan persamaan berikut:
C n =
4 sin ξ n
(29)
2ξ n+ sin(2ξ n )
T
T
T ( x , t ) h
h
x
x *0
x *1
Q max
1
bagian
0
1
DimanaC 1
1.0025 1.005 1.0099 1.0148 1.0197 1.0246 1.0483 1.0712 1.0931 1.1143 1.1345 1.1539 1.1724 1.1902 1.2071 1.3384 1.4191 1.4698 1.5029 1.5253 1.5411 1.5526 1.5611 1.5677 1.5919 1.5973 1.5993 1.6002 1.6015 1.6021
0.173 0.2445 0.345 0.4217 0.486 0.5423 0.7593 0.9208 1.0528 1.1656 1.2644 1.3525 1.432 1.5044 1.5708 2.0288 2.2889 2.4556 2.5704 2.6537 2.7165 2.7654 2.8044 2.8363 2.9857 3.0372 3.0632 3.0788 3.1102 3.1416
C 1 1.003 1.006 1.012 1.0179 1.0239 1.0298 1.0592 1.088 1.1164 1.1441 1.1713 1.1978 1.2236 1.2488 1.2732 1.4793 1.6227 1.7202 1.787 1.8338 1.8673 1.892 1.9106 1.9249 1.9781 1.9898 1.9942 1.9962 1.999 2.0
tengah plat datar, yang
(33)
ξ 1 dan V adalah volume
Contoh 3 Pada proses manufaktur plat-plat baja, setebal 4cm yang baru diproses didinginkan dengan meniupkan udara lingkungan 30 oC hingga koefisien konveksi yang terjadi pada masing-masing sisinya 100 W/m 2K. Temperatur awal plat adalah 300 oC. Anggap panjang dan lebar plat cukup besar hingga dapat dianalisis sebagai kasus transient 1 dimensi. Sifat fisik yang digunakan k =55 W/mK dan α =18×10−6m2/s. Setelah 20 menit Tentukanlah: (a)temperatur di pusat plat, (b) temperatur permukaan, dan (c) 3 energi yang dikeluarkan dari plat jika volume tiap plat 2m .
(32)
1
dapat
ξ 1
1 − θ * sin ξ 1
Dimana Qmax= ρc pV (T i − T ∞), benda.
0
dirumuskan dengan: θ =C exp(−ξ Fo)
=
0
Penyelesaian yang dinyatakan dalam bentuk deret, seperti pada persamaan (28) masih jauh dari praktis untuk digunakan. Oleh karena itu diperlukan pendekatan. Pada kasus khusus dengan bilangan Forier Fo>0,2 , persamaan (28) dapat diwakili dengan suku pertama deretnya, dengan kesalahan kurang dari 2% . θ = θ cos(ξ x ) (31) 0
C 1
Perpindahan Energi Besar energi panas yang telah dipindahkan selama waktu t dapat dihitung dengan persamaan berikut:
Dan parameter ξ n (eigenvalue) adalah akar-akar positif dari persamaan karakteristiknya: ξ ntan ξ n=Bi (30)
adalah temperatur
ξ 1
Bola
*
Gambar 6 Konveksi transien 1D pada plat datar
θ
C 1
0.0998 1.0017 0.1412 0.141 1.0033 0.1995 0.1987 1.0066 0.2814 0.2425 1.0098 0.3438 0.2791 1.013 0.396 0.3111 1.0161 0.4417 0.4328 1.0311 0.617 0.5218 1.045 0.7465 0.5932 1.058 0.8516 0.6533 1.0701 0.9408 0.7051 1.0814 1.0184 0.7506 1.0918 1.0873 0.791 1.1016 1.149 0.8274 1.1107 1.2048 0.8603 1.1191 1.2558 1.0769 1.1785 1.5995 1.1925 1.2102 1.7887 1.2646 1.2287 1.9081 1.3138 1.2403 1.9898 1.3496 1.2479 2.049 1.3766 1.2532 2.0937 1.3978 1.257 2.1286 1.4149 1.2598 2.1566 1.4289 1.262 2.1795 1.4961 1.2699 2.288 1.5202 1.2717 2.3261 1.5325 1.2723 2.3455 1.54 1.2727 2.3572 1.5552 1.2731 2.3809 ∞ 1.5708 1.2732 2.4048 Catt: Satuan ξ adalah radian
III.1 Plat datar dengan konveksi
θ =∑C nexp(−ξ n2Fo) cos(ξ n x *)
ξ 1
Silinder
0.01 0.02 0.04 0.06 0.08 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 20 30 40 50 100
=L
Persamaan (25) dengan nilai awal dan nilai batas pada persamaan (26) dan persamaan (27) akan diselesaiakan secara eksak pada bagian berikutnya.
*
Plat Datar
Bi
dihitung
Penyelesaian: Asumsi konduksi transient 1D dengan sifat fisik ko nstan. Data soal L=4 2 = 2 × 10−2 m , waktu pendinginant
denganmenggunakan
persamaan (29). Nilai C 1 dan akar pertama persamaan karakteristik ξ 1 adalah fungsi dari bilangan Biot, seperti yang ditampilkan pada Tabel 1.
Sustainable Energy Research Group
=1200s.
51
Mechanical Engineering, USU
Perpindahan Panas I
Prepared by: Himsar AMBARITA
Gunakan persamaan (32) untuk menghitung temperatur pertengahan plat. Caranya: 1. Hitung Bilngan Fo dan Bi
Fo = α t = 54, karena
> 0,2 persamaan pendekatan
2
Persamaan pendekatan yang dapat digunakan asalkan memenuhu syarat Fo=αt r 02 * * * θ = θ 0 J 0(ξ 1r )
(34) silinder, yang dirumuskan
*
θ 0 adalah temperatur
dapat digunakan. Bi = hL = 0,0364 k 2. Hitung ξ dan C 1 dengan menginterpolasi Tabel 1 Interpolasi antara 0,02 dan 0,04
Bi
1
0.02 0.0364 0.04
1
0.141 1.0033 0.188314 1.006006 0.1987 1.0066
Catatan persamaan (30) dapat juga digunakan untuk menghitung ξ . Tetapi teknik trial and error harus digunakan. 3. Hitung θ dengan menggunakan persamaan (32)
pusatdengan:
θ
C exp(−ξ Fo)
(35)
=
0
1
1
Dimana C 1 dan ξ 1 merupakan fungsi bilangan Biot yang
ditampilkan pada Tabel 1. Dan J 0 adalah fungsi Bessel enis pertama. Sementara besar perpindahan panas setelah waktu tertentu dapat dihitung dengan persamaan: 2θ * Q =1−
0
J (ξ 1
1
(36)
)
ξ 1
max
Dimana J 1 adalah fungsi Bessel jenis pertama.
0 =
θ *
1,006 exp(−0,1883
×
54) = 0,148
0
4. Hitung temperatur di pertengahan plat T−
III. 3 Bola (Sphere)
T −30 θ * = T ⇒ 0,148 = 0 T i − T 300 − 30 o Maka jawaban (a) T =70 C. 5.
T ( r , t ) r 0
Hitung temperatur permukaan dengan menghitung θ *
h *
r
θ * = θ 0*cos(ξ 1 x *) = 0,148 × cos(0,188 ×1) = 0,146 T
r
0
T
−
θ *= T T i − T
6.
r
0,146 = T −30 300 − 30 Maka jawaban (b) T =69,3oC. Catatan: perbedaan temperatur di permukaan dan ditengah plat sangat kecil. Hal ini bisa terjadi karena memang konduktivitas plat baja sangat baik. Hitung energi yang dikeluarkan dari plat dengan menggunakan persamaan (33). 9 (T − T ) = k V (T − T )=1,65×10 Q J = ρc V pi i max α Q = 1 − θ * sin ξ 1 ⇒ Q =1,41×109J
Q
⇒
0
θ
T ( r , t )
1
adalah temperatur pusat bola, yang
dirumuskan
dengan: θ =C exp(−ξ Fo) 1
(38)
1
max
=
1−
3θ 0
ξ 1
*
sin(ξ 1 ) −ξ 1 cos(ξ 1 )
(39)
3
Catatan Penyelesaian secara eksak untuk ketiga bentuk geometri sederhana ini (plat datar, silinder panjang, dan bola) hanya untuk kondisi batas konveksi di permukaan dengan udara luar yang mempunyai temperatur T .
r 0 h r
r
ξ 1r
0
Q
Penyelesaian pendekatan juga dapat dirumuskan untuk silinder dengan panjang tak terhingga seperti pada Gambar 7.
*
0
*
Sementara besar perpindahan panas setelah waktu tertentu dapat dihitung dengan persamaan:
III.2 Silinder
r
Untuk kasus bola seperti yang ditampilkan pada Gambar 7, persamaan menghitung temperaturnya sebagai fungsi waktu ditampilkan pada persamaan berikut: * θ *= θ * 1 sin( ξ r ) (37)
0
ξ 1
max
Gambar 8 Konveksi transient 1D pada Bola
T
0
Gambar 7 Konveksi transient 1D pada Silinder
Jawaban dari penyelesaian eksak adalah berupa deret yang akan konvergen pada nilai akhir. Jawaban seperti ini adalah tidak praktis, maka diajukan jawaban pendekatan dengan hanya menggunakan suku pertama dari jawaban yang berbentuk deret tersebut. Hal ini bisa diterima selama bilangan Fourier Fo≥0,2 , kesalahan yang didapat biasanya hanya kurang dari 2%. Diluar dari syarat ini, tidak
Sustainable Energy Research Group
52
Mechanical Engineering, USU
Perpindahan Panas I
Prepared by: Himsar AMBARITA
θ 0* = T (0,t )−T ∞ T (0,t )=365 C T i − T ∞
disarankan menggunakan penyelesaian pendekatan ini. Harus dicari cara lain untuk menyelesaiakannya. Masalah kondisi batas di permukaan juga masih terbatas pada kasus dengan konveksi di permukaan. Masih ada dua kemungkinan kasus dengan kondisi batas lain yang tentunya akan memberikan penyelesaiaan yang berbeda. Kedua kasus ini adalah: (a) kasus dengan kondisi batas tempertur permukaan konstan dan (b) kasus dengan kondisi batas flux panas konstan. Untuk kedua kasus ini tidak akan dibahas pada buku ini. Tetapi bagi mahasiswa yang serius ingin mempelajarinya disarankan untuk mengacu pada Referensi I.
5.
2
dengan koefisien konveksi h= 80 W/m K. Tentukan temperatur pusat poros setelah 45 menit dan panas yang dipindahkan per satuan panjang poros.
Penyelesaian: Asumsi yang digunakan poros mengalami konduksi transient 1 dimensi dan sifat fisik konstan. =14,9 Sifat material yang digunakan : W/mK,
−6 2 ρ =7900 kg/m3,c p=477J/kgK, danα =3,95×10 m /s.Data 2 soal: t = 45 menit (2700s), r = 0,1 m, h= 80 W/m K, T i =600oC, dan T ∞=200oC. L =1m (persatuan panjang). Gunakan persamaan (35) untuk menghitung temperatur pusat silinder. Tahap penyelesaian soal ini adalah: 1. Hitung bilangan Fo dan Bi
αt
Fo =
r
Maka temperatur inti akan turun dari 600 oC menjadi 365oC setelah 45menit. Jawaban ini berbeda 5oC dengan Referensi II, karena pada buku tersebut soal ini diselesaikan dengan menggunakan grafik. Panas yang dipindahkan persatuan panjang dapat dihitung dengan mengunakan persamaan (36) 2 2 V =π r L =3,14×10− m3
Q max= ρc pV (T i − T ∞) = 43,378 MJ
Q
Karena penyelesaian eksak dan pendekatan masalahmasalah konduksi transient ini sangat terbatas, maka saya lebih suka menggunakan pendekatan numerik untuk menyelesaiakan masalah konduksi transient ini. Penyelesaian secara numerik ini akan dibahas pada bab berikutnya dari buku ini.
Contoh 4 (Diambil dari Referensi II sebagai pembanding)Sebuah poros berbentuk silinder yang panjang dengan diameter 20cm terbuat dari stainless steel 304 o dikeluarkan dari oven produksi dengan temperatur 600 C. o Poros didinginkan di ruangan dengan temperatur 200 C
o
⇒
=
* 1 − 2θ 0 J (ξ ) ⇒Q= 30 MJ 1
1
ξ 1 Dimana J 1(0,969)=0,43adalah fungsi Bessel jenis pertama orde pertama. Dapat dilihat pada tabel atau dengan menggunakan fungsi Bessel pada Excell. max
IV. Benda Padat Semi infinite Di daerah yang beriklim tropis, seperti Indonesia, temperatur tanah (bumi) umumnya lebih redah dari temperatur udara lingkungan. Jika di pagi hari permukaan bumi mendapat penyinaran dari radiasi matahari, maka temperatur permukaan tanah akan naik. Hal ini akan mengakibatkan terjadinya perpindahan panas konduksi dari permukaan tanah ke dalam tanah. Proses perpindahan panas ini dapat dianggap hanya 1 dimensi dan sisi lainnya jauh tak terhingga ke dalam tanah. Kasus semacam ini disebut konduksi transient pada benda padat semi infinite. Artinya pada bagian permukaan bendanya tertentu, sementara pada bagian dalam dimensinya tak terhingga. Diagram perpindahan panasnya ditampilkan pada Gambar 9. Temperatur di dalam benda akan merupakan fungsi jaraknya dari permukaan dan waktu. Tetapi karena bendanya semi infinite, temperatur pada jarak yang menuju tak terhingga akan sama dengan temperatur a wal ( initial ).
= 1,07
T
2
hr
Bi =
=0,537 k
h
Nilai Fo sebenarnya tidak melampaui 0,2 sebagai syarat penggunaan pendekatan. Tetapi kita tidak punya alternatif lain dan memberi asumsi tambahan berlaku pendekatan agar persamaan (35) tetap dapat digunakan.
ξ 1
0.50
0.941 0.537
T ( x ,t )
0.60
1.018
x T ( x → ∞,t )= T i
C 1 1.114
0.969
Gambar 9 Konduksi transient pada benda semi-infinite
1.122 1.135
Seperti yang ditunjukkan pada gambar ada tiga kemungkinan kondisi yang mungkin terjadi, yaitu (a) temperatur permukaan konstan T s, (b) terjadi konveksi dengan udara lingkungan, dan (c) flux panas konstan. Penyelesaian untuk kasus ini akan disesuaikan dengan kondisi batas ini.
*
3. Hitung θ dengan menggunakan persamaan (35) 0
= C exp(−ξ Fo)=0,412
θ 0
1
q ′′
Ts
2. Hitung ξ dan C 1 dengan menginterpolasi Tabel 1 (untuk silinder) Interpolasi antara 0,5 dan 0,6
Bi
1
4. Hitung temperatu di pusat dengan menjabarkan θ 0
Sustainable Energy Research Group
53
Mechanical Engineering, USU
Prepared by:
Perpindahan Panas I
Himsar AMBARITA
a.
Temperatur permukaan konstan s Penyelesaian konduksi transient pada benda semiinfinite ini dapat dilakukan dengan menggunakan konsep similarity variable. suatu Dimisalkan variabel
q′′ = k (T s − T i )
παt b.
Permukaan semi infinite dengan konveksi
η = x (4α t) 1 2yang harus dimasukkan ke persamaan (21). 2
Untuk kasus dengan kondisi batas konveksi
1 ∂T
∂ T 2
permukaan, persamaan temperatur dapat dihitung dengan
∂ x = α ∂t
(21)
Dengan menggunakan chain rule akan didapat persamaan persamaan berikut: ∂T dT ∂η =
∂ x
(43)
dη ∂ x
1
2
4αt dη
=
∂T =
dan
∂ x
dT ∂η
=
dT 2
2t 4αt dη d η∂t Substitusi persamaan-persaman ini ke persamaan (21)
akan memberikan: 2
d T = −2η dT d η2 d η Nilai batas untuk persamaan ini adalah: T (η=0)= T s dan T (η→ ∞) =T i
h α t
(44)
2
+
x
h
(45)
αt
αt k Pada persamaan ini digunakan erfc yang artinya fungsi error tambahan (complement error function ). Artinya jika diketahui suatu bilangan sembarang u , maka erfc u= 1 − erf u.
(40)
c. (41)
Persamaan (40) dapat diubah bentuknya menjadi bentuk yang dapat diintegralkan sekali, yaitu: d dT dη
d η Jika diintegralkan, menjadi:
+
k
k
α
Dimana EC = erfc
4αt dη
∂t
dT
− EC exp
1 d T
x
=−
= erfc
i
T ∞− T i
2
∂ T
dT
cara yang sama dan hasilnya adalah: T − T x hx
Flux panas konstan
Untuk kasus dengan flux panas konstan, persamaan temperaturnya adalah: 2q′′ αt π T − T =
s k
−
exp 4αt
x
q′′ x
2
− x
s k
(46)
erfc 2 αt
= − 2ηdη
lndT dη= −η +C Pada dT 2 = C 1exp ( − η Atau ) d η Integralkan sekali lagi akan memberikan: (agar berbeda, variabelnya akan diganti dengan sebuah variabel dummy u) T = C 1 ∫∞ exp(−u2 )du+C 2 0
Dengan menggunakan kondisi batas yang pertama di persamaan (41), maka akan di dapat C 2 = T s. Masukkan kondisi batas kedua akan didapat: T i = C 1∫∞exp(−u2)du + T s 0
Jika diarrange, bentuknya bisa ditulis menjadi persamaan berikut: (T i −T s ) 2(T i −T s ) ∞ 2 C 1= ∫ exp(−u )du = π
Contoh 5 (Diambil dari Referensi II) daerah-daerah yang mempunyai musim dingin temperatur udara bisa mencapai temperatur di bawah 0 o C hal ini bisa menimbulkan masalah pada pipa aliran air, dimana air akan membeku di dalam pipa yang akan menghentikan aliran air dan dapat juga menyebabkan pipa pecah. Untuk menghindari hal ini terjadinya biasanya pipa ditanam di dalam tanah yang lebih hangat. Masalahnya adalah temperatur permukaan tanah akan berangsur-angsur turun dan jika kedalaman pipa tidak cukup dalam tetap akan terjadi pembekuan. Misalkan suatu pada suatu daerah permukaan tanah tertutup salju selama 3 bulan (lamanya musim dingin) hingga permukaannya sekitar -10 oC. Temperatur awal tanah saat memasuki musim dingin 15oC. adalah tanah Jika koefisien konduksi 6 − 2 α =0,15×10 k =0,4 W/mK dan m /s tentukan kedalaman minimum pipa air harus ditanam di daerah tersebut untuk menghindari pembekuan air.
0
Dengan menggunakan konstanta ini, maka jawaban untuk persamaan (40) adalah: T − T s
=
2
∫
exp(−u2 )du= erf η
(42)
η
0
T i − T s π Dimana bagian yang di sebelah kanan disebut sebagai fungsi error Gauss. Sementara heat flux perpindahan panas pada permukaand apatdihitun gdenganmenggun akan persamaan: = −k (T − T q′′ = − ∂T ∂η ) d (erf η) s
i
∂ x = 0
s
Jika diselesaiakan akan memberikan:
d η∂ x
η =0
Penyelesaian: Diagram contoh soal ini ditampilkan pada gambar berikut:
Sustainable Energy Research Group
54
Mechanical Engineering, USU
Perpindahan Panas I
Prepared by: Himsar AMBARITA
o T = −10 C
Teknik Superposisi
s
x Pipa air
r 0
Asumsi yang digunakan: sifat fisik konstan, tempertur permukaan tanah konstant pada -15 oC, dan 1 dimensi. Agar air di dalam pipa tidak membeku tempertur paling tidak di atas titik beku, T ( x,t ) = 0oC. Persamaan yang sesuai untuk soal ini adalah persamaan (39). 1. Hitung lamanya musim dingin dalam detik t 6 =3×30×24×60×60=7,776×10 s 2. Gunakan persamaan (39) T − T x s
Syarat utama menggunakan analisis transient 1 dimensi adalah dimensi benda yang dibahas tak terhingga ke di kedua sisi arah dimensi yang dihilangkan. Misalnya kalau distribusi temperaturnya adalah T ( x , t ) , karena
y dan z dihilangkan maka panjang benda ke arah y dan z adalah tak terhingga. Atau kalau kasusnya semi infinite,panjang benda tak terhingga hanya ke satu sisi dari dimensi yang dihilangkan. Sementara jika ukuran benda yang dibahas sudah tertentu, seperti silinder dengan panjang L dan
= erf
T i − T s x erf
4α t T − T =
s
0 − (−10) =
= 0,4
T i − T s 15 − (−10) 4α t Nilai error function yang sesuai adalah jika parameternya 0,371. Nilai ini bisa dilihat di tabel atau dengan menggunakan fungsi Erf (x) pada excell. 3.
Gambar 10 Teknik superposisi membentuk silinder sepanjang L
jari-jari r 0 seperti pada Gambar 10(b), kasusnya sudah
merupakan multi dimensi, tepatnya T ( x , r , t ) . Pada Gambar 10(a) ditunjukkan bahwa silinder dengan panjang
L dan jari-jarir 0merupakan irisan dari plat datar yangpanjang dengan tebal W ( W =L ) dengan silinder sangat panjang dengan yang berjari-jari r 0 .
Hitung x x = 0,371 4α t
x =0,371×
Dengan menggunakan analogi irisan ini, maka distribusi temperatur pada silinder sepanjang L dapat dirumuskan dengan persamaan:
4 × 0,15 × 10−6× 7,776 × 106= 0,8 m
θ *( x ,r ,t )
Maka jika pipa ditanam pada kedalaman lebih dari 0,8 m tidak akan membeku meskipun temperatur di permukaan
Sil , L
= θ *( x ,t )
Plat
− 10 Cselama 3 bulan karena tertutup salju.
θ * x ,r ,t
Sil , L
V. Transient Multi Dimensi
sedangkan hanya kasus 1 dimensi T ( x , t ) kita sudah hampir menyerah dan minta tolong pada jawaban pendekatan, apakah mungkin menyelesaikan T ( x , y , z , t ) lagi? Memang
θ * x ,t plat
tanpa dimensi yang
= (T ( x , r , t ) −T ∞ ) (T i − T ∞) (T ( x , t ) −T ∞ )
=
(T i −T ∞ ) (T (r , t ) −T )
*
θ r , t Sil =
(T i −T ) Bagaimana menghitung parameter-parameter yang disebelah kanan telah dibahas pada bagian sebelumnya. Analogi irisan juga dapat digunakan untuk menghitung besar perpindahan panas pada kasus-kasus multi dimensi. Persamaan untuk kasus dua dimensi adalah:
Q
benar, untuk saat ini hampir tidak mungkin menyentuh jawaban T ( x , y , z , t ) secara eksak. Yang akan kita lakukan
=
max
Q
+ max
2D
55
Q
1− max
1
Dan untuk tiga dimensi:
di bagian ini adalah teknik super posisi. Artinya kita akan menggabung jawaban-jawaban dari kasus-kasus 1 dimensi untuk mendapatkan jawaban multi dimensi. Teknik inilah yang disebut teknik superposisi, yang diperkenalkan oleh Langston 1982.
Sustainable Energy Research Group
(47)
Sil
Dimana θ * adalah temperatur dirumuskan sebagai berikut:
o
Pada bagian awal kita telah merumuskan konduksi transient tanpa pengaruh posisi dan bagian kedua konduksi transient dengan pengaruh posisi hanya 1 dimensi. Pada bagian ini kita akan membahas kasus transient dengan multi dimensi. Dengan kata lain kita ingin menyelesaian distribusi temperatur T ( x , y , z , t ) . Mungkin anda bertanya:
× θ *(r ,t )
Q (48) max
2
1
Mechanical Engineering, USU
Perpindahan Panas I
Prepared by: Himsar AMBARITA
Q
Q
=
Q
1−
+
max
max
max
3D
+
Bi =
1− max
(49)
1
Q
Q
1−
max 1
Penyelesaian: Soal ini tidak dapat lagi diselesaikan dengan mengunakan asumsi silinder sangat panjang (kasus 1D), karena panjang poros sudah tertentu. Soal ini harus diselesaikan dengan menggunakan teknik superposisi antara poros sangat panjang dan plat datar sangat panjang. Asumsi yang digunakan di sini adalah sifat fisik konstan. Data yang diberikan soal: Poros silinder r =2,5 cm, o o L =15cm, T =500 C, T =30 C, dan h =60W/m K. Sifat
∞
digunakan :k =14,9W/mK,
materialyang 3
ρ =7900 kg/m ,c p=477J/kgK, dan α =3,95×10
−6
2
m /s.
Cara menyelesaiakan soal ini: 1. Hitung bilangan Fo dan Bi untuk silinder panjang
αt
Fo =
r
Hitung ξ dan C 1 dengan menginterpolasi Tabel 1 Interpolasi antara 0,2 dan 0,3
ξ 1
Bi
= 11,4
2
hr
Bi
ξ 1
0.1
0.4417 0.101
0.2
C 1 0.443
0.617
1.0246 1.025 1.0483
Hitung θ 0* dengan menggunakan persamaan (32)
7.
θ 0*=C 1exp(−ξ 12Fo)=0,604 Gabungkan kedua parameter ini untuk mendapatkan θ 0*untuk silinder dengan menggunakan persamaan(47):
θ *(0,0,t )
3. Hitung θ 0 dengan menggunakan persamaan (35) θ 0*= C 1exp(−ξ 12 Fo)=0,11
Sil , L
Atau θ 0* Sil ,L 8.
=
*
θ (0,t )
=
Plat
θ *= T − T ∞ 0
× θ *(0, t )
Sil
=
0,604 × 0,11
0,0664
Hitung temperatur menjabarkan θ 0* .
di
pusat
0,0664 = T −30
⇒
silinder
⇒
T
dengan o
61,2 C
=
− T ∞
500 − 30 Angka ini merupakan jawaban untuk soal (a) i
9.
Panas maksimum yang berpindah dapat dihitung dengan menggunakan persamaan (48)
Q
=
max
Q
+
Q
max
2D
1− max
1
2
10. Hitung Q Qmax1dengan persamaan(33): sin ξ Q * 1 =0,415 = 1 − θ 0 max
=1−
max
Q max
menggunakan
1
1
11. Hitung Q Qmax2dengan persamaan(369: 2θ * Q 0
ξ 1
2
menggunakan
J 1 (ξ 1 ) =0,893
12. Gabungkan kedua parameter ini: Q
= 0,415 + 0,893(1 − 0,415)=0,9374
max
sil ,L
13. Hitung perpindahan panas maksimum − V = πr 2× L =1,96×10 4m3
= ρc pV (T i − T ∞)=347,58kJ max
*
1.0033 1.006006 1.0066
6.
Bi =
2.
=0,101 k Hitung ξ dan C 1 dengan menginterpolasi Tabel 1 (untuk silinder) Interpolasi antara 0,1 dan 0,2
C 1
0.02 0.141 0.0364 0.188314 0.04 0.1987
1
Contoh 6 Misalkan kelak, anda berminat menjadi pengusaha dengan mendirikan pabrik untuk memproduksi poros yang akan anda pasarkan ke pabrik-pabrik kelapa sawit di Indonesia. Salah satu poros produksi anda berbentuk silinder mempunyai diameter 5cm dan panjang 10cm terbuat dari stainless steel 304. Setelah mengalami proses produksi dengan temperatur awal 500 oC, poros-poros tersebut anda dinginkan dengan meniupkan udara lingkungan 30 oC hingga terjadi konveksi dengan koefisien 60 W/m 2K. Setelah 30 menit, tentukanlah temperatur inti poros dan total panas yang sudah berpindah dari poros
i
5.
max
3
= 0,201
k
max
2
1
Q
hLc
Q
14. Hitung perpindahan panas yang telah terjadi: Q
Q=
×Q
ax
= 325,8kJ
max
4.
Hitung Fo dan Bi untuk kasus plat datar panjang Di sini tebal plat sama dengan panjang silinder (perhatikan Gambar 10). Maka Lc =5 cm
Fo =
α t Ls
2
= 2,84
VI. Beberapa bentuk lain multi dimensi Beberapa bentuk multi dimensi yang merupakan irisan dari bentuk 1 dimensi dapat dilihat pada sub bagian berikut: VII. Soal –soal Latihan 1. Soal Latihan Menyusul
TUGAS : 1. Sebuah pipa horizontal yang memiliki suhu permukaan ( 67 C + 2 NBI AKHIR ) dan diameter ( 25 cm + 2 NBI AKHIR ) dikubur pada kedalaman 1,2 m dalam tanah yang memiliki konduktivitas termal ks = 1,8 W/m.C. Suhu permukaan tanah adalah 15 C. Hitunglah kehilangan panas per satuan panjang.
Penyelesaian: Diketahui:
T1 = 67 + 81 = 148 0C
T2 = 15 0C
D = 25 + 81 = 106cm = 1,06m Z = 1,2 m Ks = 1,8 W/m . 0C Karena nilai z < 3D/2 maka sesuai batasan didalam table factor ukuran L harus lebih besar daripada D. dan bila sesuai dengan persatuan panjang maka ukuran L = 1m Ditanya : … Q ? Jawab:
Q = K. S . (T1-T2) Nilai S (faktor bentuk)
. . = − .,.., = − , = , − , = , , = , , = ,
Maka Q = K. S . (T1-T2) Q = 1,8 .12,56 . (148 – 15) Q = 3006,864 W
2. Sebuah as berbentuk batangan silinder panjang berdiameter 20 cm yang terbuat dari baja stainless 304 keluar dari sebuah oven pada suhu seragam ( 600 C + 2 NBI AKHIR ). Batangan ini kemudian dibiarkan mendingin secara lambat di dalam suatu ruangan bersuhu 200 C dengan koefisien perpindahan panas konveksi rata rata h = 80 W/m2.C. Tentukan suhu di pusat poros setelah 30 menit. Juga tentukan perpindahan panas per satuan panjang poros selama itu ( 30 menit ).
Penyelesaian: Diketahui:
D = 20cm = 0,2 m
Ti = 600 + 81 = 681 0C
Cp = 477 J/Kg 0C
T∞ = 200 0C
K = 14,9 W/m 0C
ρ =
h = 80 W/ m 20C
α = 3,95 x 10 -6 m2/s
t = 30 menit = 1800s
L = 1m (persatuan panjang)
7900 Kg/m 3
Ditanya : a. Pada t = 30 menit, berapakah suhu bagian tengan silinder b.Perpindahan panas persatuan panjang poros selama 30 menit?
Jawab: a. Hitung bilangan Bi dan Fo
ℎ . = ., =0,537 ,
Bi =
− 1800 0,00711 . 3, 9 5 10 = = 0,1 = 0,01 =0,711
Dengan Bi = 0,537 maka system lumped tidak bisa dipakai, Fo > 0,2, maka kita bisa menghitung suhu permukaan dengan pendekatan satu suku. Dalam hal ini kita melihat table di bagian silinder. b. Hitung λ1 dan A1Dari Table Bi dan Fo didapat:
Bi
λ1
A1
0,5
0,9408
1,1143
0,537
0,9695 1,1217
0,6
1,0184
1,1345
∞ = − ,= ∞
c. Hitung
dengan menggunakan persamaan berikut:
,=1,1217 −, ,
= 0,276
d. Hitung temperature pussat dengan menjabarkan
0,1800 = ∞+ 1∞ 0,1800 = 200+ 0,276 681200
= 332,75 0C
Maka temperature pada T (0,1800)= 332,75 0C e. Panas yang dipindahkan ppersatuan panjang dapat dihitung dengan
menggunakan persamaan
= . . = ∞ = 7900 .477 .3,14 10−681200 = 56,91 =1 2 = 1 2 276 0,9695 = 56,91 1 20,.0,9695 = 3,14 x 10 -2
NIlai J1 didapat dari hasil interpolasi berikut: Bi
J1
0,9
0,4059
0,9695 0,4295 1,0
0,4400
52 .0,43) = 56,91 (1 0,0,95695 = ,
3. Sebuah lempeng besi yag mulanya berada pada suhu seragam ( 550 F + 2 NBI AKHIR ) tiba tiba diturunkan suhu permukaannya menjadi 100 F. Hitunglah waktu yang diperlukan untuk mencapai suhu 200 F pada suatu kedalaman 5 cm. Hitung pula total panas yang dikeluaarkan dari besi per satuan luas selama waktu itu
Penyelesaian: Diketahui:
Ti = 550 + 81 = 631 0F T0= 100 0F T (x,t) = 200 0F X = 5 cm = 5 x 0,032 = 0,16 ft α = 23,1 x 10 -6 m2/s = 23,1 x 10 -6 x 38750,1 = 0,89 Ft 2 / jam K = 80,2 W/m K = 80,2 x 0,58 = 46,5 Btu/Jam. ft. F
Ditanya:
a. berapa waktu yang diperlukan (t) … ? b. Total panas yang dikeluarkan persatuan luas selama t … ?
Jawab :
, =erf 2√ . , = 200100 = 0,188 = erf 631100 2√ . Lihat table nilai fungsi error Gauss pada angka 0,188, maka didapat angka.
2√ .
2√ .
0,17901
0,16
0,188
0,1681
0,20094
0,18
√ . 0, 1 6 = [2 0,191]0,89 / = 0,10,450256.0,89 = 0,129 = 7,7 0,1681 =
Masukkan harga α dan x, maka t yang diperlukan dapat dicari :
Total panas persatuan luas
Δ = 2. . ⁄ Δ = 2.46,5 100622 3,14 . 0,89⁄0,129 Δ = 48546 =10440 ⁄ 4,65
