Pengujian Karakteristik Perpindahan Panas Dan Faktor Gesekan Pada Penukar Kalor Pipa Konsentrik Dengan Sisipan Pita Terpilin Berlubang
Pengujian Karakteristik Perpindahan Panas Dan Faktor Gesekan Pada Penukar Kalor Pipa Konsentrik Dengan Sisipan Pita Terpilin Berlubang
Pengujian Karakteristik Perpindahan Panas Dan Faktor Gesekan Pada Penukar Kalor Pipa Konsentrik Dengan Sisipan Pita Terpilin Berlubang
Keywords : Abstract :
Nusselt number This research was conducted to examine the characteristics of heat transfer
Reynolds number and friction factor in the annular channel concentric tube heat exchanger
friction factor with a classic twisted tape insert and perforated twisted tape inserts. Test
perforated twisted tape insert section was the single pass concentric tube heat exchanger with inner tube
and outer tube made of aluminum. Dimensions of outer tube; outer diameter
of 21.87 mm and inner diameter of 20.67 mm, and dimensions of inner tube,
outer diameter of 15.84 mm and inner diameter of 14.34 mm. The length of
heat exchanger was 2,000 mm and the length of pressure drop measurement
in the inner tube was 2,010 mm. Flows in the inner tube and in annulus was
counter flow. Working fluid in the inner tube was hot water which its inlet
temperature was maintained at 60°C, whereas in the annulus was cold water
at 28o C. Both of classic twisted tape insert and perforated twisted tape insert
with a twist ratio of 4.0 were made of aluminum strips with a thickness of
0.76 mm and width of 12.61 mm which it was twisted so forming twist with
the length of pitch was 50.35 mm. Perforated twisted tape inserts were varied
with holes diameter of 4 mm and 6.5 mm, respectively, where the distance
between the center holes of 4 cm. Twisted tape insert installed in the inner
tube of the concentric tube heat exchanger. The results showed that at the
same Reynolds number, the addition of classic twisted tape insert and
perforated twisted tape insert with holes diameter of 4 mm and 6,5 mm in the
inner tube increasing the average Nusselt numbers were 79,20%, 68,7%, and
57,8% than the inner tube without a twisted tape insert (plain tube),
respectively. At the same pumping power, the average Nusselt number in the
inner tube with the addition of perforated twisted tape inserts with holes
diameter of 6.5 decreased 1.62%, whereas the addition of perforated twisted
tape inserts with holes diameter of 4 mm classic twisted tape inserts
increased the average Nusselt number were 3.2% and 5,84% compared with
plain tube, respectively. At the same Reynolds number, the addition of classic
twisted tape inserts and perforated twisted tape insert with holes diameter of
4 mm and 6,5 mm in the inner tube increased the average friction factor
339.25%, 298.49% and 269.6% than plain tube, respectively. At the same
pumping power, the addition of classic twisted tape inserts and perforated
twisted tape insert with holes diameter of 4 mm and 6.5 mm in inner tube
increased the average friction factor 485.39%, 416.48 % and 362.72% than
plain tube, respectively.
isian material yang dapat membuat aliran di dalam Istanto et al (2011) menyelidiki karakteristik
pipa penukar kalor menjadi lebih turbulen. Di sisi perpindahan panas dan faktor gesekan pada penukar
lain, semakin turbulen sebuah aliran fluida maka kalor pipa konsentrik dengan penambahan sisipan
penurunan tekanan (pressure drop) yang terjadi pita terpilin panjang setengah penuh dan panjang
antara sisi masuk dan sisi keluar dari aliran fluida penuh. Istanto et al (2011) menyelidiki karakteristik
tersebut semakin besar. perpindahan panas dan faktor gesekan pada penukar
Penurunan tekanan berpengaruh pada besarnya kalor pipa konsentrik saluran persegi dengan
energi yang harus diberikan pompa atau fan untuk penambahan sisipan pita terpilin klasik dan
mengalirkan fluida dalam penukar kalor atau disebut berlubang. Istanto et al (2011) menyelidiki pengaruh
daya pemompaan (pumping power). Semakin besar rasio pilinan (twist ratio) terhadap karakteristik
penurunan tekanan, maka semakin besar daya perpindahan panas dan faktor gesekan pada penukar
pemompaan yang diperlukan. Penurunan tekanan kalor pipa konsentrik dengan penambahan sisipan
fluida mempunyai hubungan langsung dengan pita terpilin.
perpindahan panas dalam penukar kalor, operasi, Penelitian ini dilakukan untuk menguji
ukuran, dan faktor – faktor lain, termasuk karakteristik perpindahan panas dan faktor gesekan
pertimbangan ekonomi. Oleh sebab itu peningkatan pada penukar kalor pipa konsentrik dengan geometri
koefisien perpindahan panas konveksi dengan pipa dalam (inner tube) dan pipa luar (outer tube)
meningkatkan turbulensi aliran fluida dalam pipa adalah annular dengan penambahan sisipan pita
penukar kalor harus dikaitkan dengan besarnya terpilin klasik dan sisipan pita terpilin berlubang di
penurunan tekanan yang dihasilkan. Efektivitas pipa dalam.
optimum dari modifikasi ini adalah nilai koefisien
perpindahan panas konveksi yang tinggi diikuti METODE PENELITIAN
dengan kenaikan penurunan tekanan yang minimum. Skema alat pengujian karakteristik perpindahan
Salah satu jenis sisipan yang banyak digunakan panas dan faktor gesekan pada penukar kalor pipa
untuk meningkatkan perpindahan panas dalam konsentrik dengan penambahan sisipan pita terpilin
sebuah penukar kalor adalah sisipan pita terpilin di pipa dalam dapat dilihat pada gambar 1. Peralatan
(twisted tape insert). Pipa penukar kalor dengan penelitian terdiri dari 3 sistem, yakni sistem
penambahan sisipan pita terpilin telah digunakan pengukuran, sistem lintasan pipa dalam, dan sistem
secara luas sebagai alat pemutar aliran secara lintasan annulus. Lintasan aliran air panas di pipa
kontinyu untuk meningkatkan laju perpindahan dalam adalah sebuah lintasan tertutup. Air panas
panas. Penambahan sisipan pita terpilin dalam yang berada di tangki air panas digerakkan oleh
sebuah pipa penukar kalor merupakan metode pasif, pompa air, mengalir melewati pipa dalam dan
dimana tidak memerlukan tambahan energi dari luar. kembali ke tangki air panas. Lintasan aliran air
Sisipan pita terpilin sering digunakan dalam penukar dingin di annulus adalah lintasan terbuka. Aliran air
kalor karena harganya murah, perawatannya mudah, dingin menggunakan metode gravitasi, yaitu aliran
dan ringkas. air dingin berasal dari tangki air dingin yang
Banyak peneliti yang menyelidiki tentang diletakkan di pada ketinggian tertentu. Air dingin
peningkatan perpindahan panas pada sebuah pipa yang keluar dari annulus langsung dibuang.
atau pada penukar kalor menggunakan sisipan pita Skema seksi uji dapat dilihat pada gambar 2.
terpilin. Sarma et al (2002) menyelidiki pendekatan Penukar kalor berupa pipa konsentrik satu laluan
baru untuk memprediksi koefisien perpindahan dengan arah aliran-aliran fluida melalui pipa dalam
panas konveksi dalam sebuah pipa dengan dan annulus berlawanan arah (counter flow). Pipa
penambahan sisipan pita terpilin untuk berbagai dalam dan pipa luar terbuat dari aluminium, dengan
perbandingan nilai pitch dan diameter. Naphon panjang berturut-turut 2.300 mm dan 1.940 mm.
(2006) menyelidiki karakteristik perpindahan panas Dimensi pipa luar; diameter luar 21,87 mm dan
dan penurunan tekanan pada penukar kalor pipa diameter dalam 20,67 mm, dan dimensi pipa dalam;
ganda dengan dan tanpa sisipan pita terpilin. diameter luar 15,84 mm dan diameter dalam 14,34
Noothong et al (2006) menyelidiki pengaruh mm. Panjang penukar kalor 2.000 mm dan jarak
penambahan sisipan pita terpilin terhadap pengukuran beda tekanan di pipa dalam 2.010 mm.
karakteristik perpindahan panas dan faktor gesekan Ukuran celah annulus 2,42 mm dengan diameter
dalam sebuah penukar kalor pipa konsentrik. Rahimi hidrolik 4,83 mm. Panjang pengukuran beda tekanan
et al (2009) menyelidiki karakteristik perpindahan di pipa dalam 2.002 mm. Pengujian dilakukan
panas dan faktor gesekan pada penukar kalor dengan dengan arah penukar kalor mendatar. Temperatur
penambahan classic, perforated, notched dan jagged masukan air panas melalui pipa dalam dijaga konstan
twisted tape insert. Ahamed et al (2007) menyelidiki 60oC, sedangkan temperatur masukan air dingin di
perpindahan panas aliran turbulen pada sebuah pipa annulus adalah ± 28oC.
dengan penambahan sisipan pita terpilin berlubang.
MEKANIKA 9
Volume 10 Nomor 1, September 2011
Gambar 1. Skema alat uji penukar kalor pipa konsentrik dengan sisipan pita terpilin
DAFTAR NOTASI
Gambar 11. Hubungan f dengan Re pada daya Ai = luas permukaan dalam pipa dalam (m2)
pemompaan yang sama Ao = luas permukaan luar pipa dalam (m2)
e. Karakteristik unjuk kerja termal (η η) Cp,c = panas jenis air dingin di annulus
Karakteristik unjuk kerja termal untuk pipa (kJ/kg.oC)
dalam dengan penambahan sisipan pita terpilin dapat Cp,h = panas jenis air panas di pipa dalam
dilihat pada gambar 12. Dari gambar 12, dapat (kJ/kg.oC)
dilihat bahwa penambahan sisipan pita terpilin klasik di = diameter dalam pipa dalam (m)
di pipa dalam menghasilkan η yang lebih baik do = diameter luar pipa dalam (m)
dibandingkan penambahan sisipan pita terpilin f = faktor gesekan
berlubang. Nilai η rata-rata pipa dalam dengan hi = koefisien perpindahan panas konveksi
penambahan sisipan pita terpilin klasik adalah 1,05 , rata-rata di pipa dalam (W/m2.oC)
sedangkan penambahan sisipan pita terpilin ho = koefisien perpindahan panas konveksi
berlubang diameter lubang 4 mm, η rata-rata sebesar rata-rata di annulus (W/m2.oC)
1,03. Sedangkan untuk penambahan sisipan pita hp = koefisien perpindahan panas konveksi
terpilin diameter lubang 6,5 mm, η rata-rata sebesar rata-rata di pipa dalam tanpa sisipan pita
0,98. Penambahan sisipan pita terpilin klasik di pipa terpilin (W/m2.oC)
dalam menghasilkan η terbesar. Fenomena ini serupa hs = koefisien perpindahan panas konveksi
dengan penelitian Rahimi dkk (2009). Hal ini rata-rata di pipa dalam dengan sisipan pita
menunjukkan bahwa sisipan pita terpilin klasik lebih terpilin (W/m2.oC)
baik sebagai turbulator untuk pipa dalam dari ki = konduktivitas termal rata-rata air panas di
penukar kalor pipa konsentrik dalam meningkatkan pipa dalam (W/m.o C)
perpindahan panas pada daya pemompaan yang Lt = panjang pengukuran penurunan tekanan
sama. di pipa dalam (m)
= laju aliran massa air dingin di annulus
(kg/s)
= laju aliran massa air panas di pipa dalam
(kg/s)
Nui = bilangan Nusselt rata-rata di pipa dalam
Nuo = bilangan Nusselt rata-rata di annulus
Qc = laju perpindahan panas di annulus (W)
Qh = laju perpindahan panas di pipa dalam (W)
Q loss = kehilangan panas konveksi di pipa dalam
(W)
Gambar 12. Hubungan η dengan Re
Re = bilangan Reynolds
Tb,o = temperatur bulk rata-rata di annulus (oC)
KESIMPULAN
Tc,in = temperatur air dingin masuk annulus (oC)
Penambahan sisipan pita terpilin klasik di pipa
Tc,out = temperatur air dingin keluar annulus (oC)
dalam dari penukar kalor pipa konsentrik saluran
Th,in = temperatur air panas masuk pipa dalam
annular menghasilkan laju perpindahan kalor,
(oC)
penurunan tekanan, faktor gesekan dan unjuk kerja
Th,out = temperatur air panas keluar pipa dalam
termal yang lebih besar dibandingkan dengan
(oC)
,
penambahan sisipan pita terpilin berlubang untuk
keseluruhan nilai bilangan Reynolds. Penambahan = temperatur rata-rata dinding luar pipa
sisipan pita terpilin berlubang dapat menghasilkan dalam (oC)
penurunan tekanan yang lebih kecil dibandingkan Ui = koefisien perpindahan panas overall
penambahan sisipan pita terpilin klasik. Nilai unjuk berdasarkan permukaan dalam pipa
kerja termal untuk penambahan sisipan pita terpilin dalam (W/m2.oC)
Uo = koefisien perpindahan panas overall
MEKANIKA 14
Volume 10 Nomor 1, September 2011
berdasarkan permukaan luar pipa dalam tape insert”, Prosiding Seminar Nasional
(W/m2.o C) Teknoin 2011, Fakultas Teknologi Industri
V = kecepatan rata-rata air panas di pipa UII, 19 November 2011
dalam (m/s) Manglik, R. M. and Bergles, A. E., 1993, Heat
A = laju aliran volumetrik air panas di pipa transfer and pressure drop correlations for
dalam (m3/s) twisted tape insert in isothermal tubes. Part 1:
<=>?= = daya pemompaan (W) laminar flows. Trans. ASME, J. Heat
Transfer, ,Vol. 116, pp. 881–889.
ρ = densitas air panas di pipa dalam (kg/m3)
Manglik, R. K. and Bergles, A. E., 1993, Heat
η = unjuk kerja termal transfer and pressure drop correlations for
µ = viskositas dinamik air panas di pipa twisted-tape inserts in isothermal tubes: Part
dalam (kg/m.s) II: Transition and turbulent flows. Trans.
∆P = penurunan tekanan (Pa) ASME, J. Heat Transfer, Vol. 115, pp.890–
∆TLMTD = beda temperatur rata-rata logaritmik (oC) 896.
Murugesan, P, Mayilsamy, K., Suresh, S.,
DAFTAR PUSTAKA Srinivasan, P.S.S., 2009, Heat transfer and
Ahamed, J.U., Rashid Sarkar, M.A., Khan, M.H., pressure drop characteristics of turbulent flow
Wased, M.A., 2007, Heat transfer in turbulent in a tube fitted with trapezoidal-cut twisted
flow through tube with perforated twisted tape insert, International Journal of Academic
tape insert, Proceedings of the International Research, 1, 123-128
Conference on Mechanical Engineering Naphon, P., 2006, Heat transfer and pressure drop in
(ICME2007) 29-31 Desember 2007, Dhaka, the horizontal double pipes with and without
Bangladesh. twisted tape insert, International
Incropera, Frank, P., David P. De Witt, 2007, Communications in Heat and Mass Transfer,
Fundamentals of Heat and Mass Transfer, 6th Vol. 33, pp. 166–175
Edition, John Wiley & Sons Noothong, W., Eiamsa-ard, S., and Promvonge, P.,
Istanto, T., Juwana, W.E., Yaningsih, I., (2011), 2006, Effect of Twisted-tape Inserts on Heat
Pengujian karakteristik perpindahan panas Transfer in a Tube, The 2nd Joint
dan faktor gesekan pada penukar kalor pipa International Conference on “Sustainable
konsentrik dengan half length dan full length Energy and Environment (SEE 2006), 21-23
twisted tape insert, Prosiding Seminar November 2006, Bangkok, Thailand
Nasional Teknik Mesin 6, Universitas Kristen Rahimi, Masoud, Sayed Reza Shabaniana, Ammar
Petra Surabaya, 16 Juni 2011. Abdulaziz Alsairafib, 2009, Experimental and
Istanto, T., Juwana, W.E., Yaningsih, I., (2011), CFD studies on heat transfer and friction
Pengujian karakteristik perpindahan panas factor characteristics of a tube equipped with
dan faktor gesekan pada penukar kalor pipa modified twisted tape inserts, Chemical
konsentrik dengan classic twisted tape insert Engineering and Processing: Process
dan perforated twisted tape insert, Prosiding Intensification, 48, 762–770.
Seminar Nasional Thermofluid 2011, Jurusan Sarma, P.K., Subramanyam, T., Kishore, P.S.,
Teknik Mesin dan Industri Fakultas Teknik Dharma Rao, V., Kakac, S, 2002, A new
UGM Yogyakarta, 4-5 Oktober 2011 method to predict convective heat transfer in
Istanto, T., Juwana, W.E., Yaningsih, I., (2011), a tube with twisted tape inserts for turbulent
Pengaruh twist ratio terhadap karakteristik flow, International Journal of Thermal
perpindahan panas dan faktor gesekan pada Sciences, 41, 955–960.
penukar kalor pipa konsentrik dengan twisted