PERHITUNGAN AWAL DESAIN TERMAL PENUKAR PANAS SISTEM PENDINGIN RRI-50 Sukmanto Dibyo, Gregorius Bambang Heru, Pusat Teknologi Keselamatan Reaktor Nuklir
[email protected]
ABSTRAK PERHITUNGAN AWAL DESAIN TERMAL PENUKAR PANAS SISTEM PENDINGIN RRI-50. Desain proses aliran sistem pendingin RRI-50 telah dibuat. Makalah ini melengkapi dokumen untuk desain komponen unit penukar panas sistem pendingin reaktor. Dua unit penukar panas dipilih jenis shell-tube dan didesain dengan kapasitas masing-masing 25,1MWt. Desain dikerjakan dengan metoda penentuan koefisien transfer panas pada penukar panas. Dalam desain ini, rugi tekanan aliran merupakan salah satu faktor yang juga harus dipenuhi. Data masukan yang digunakan adalah temperatur dan laju aliran dari desain proses sistem pendingin RRI-50.Hasil desain telah diperoleh sebagai berikut: tipe shell-tube (1-2 pass) aliran lawan arah dilengkapi 6 baffle penyekat dengan koefisien transfer panas 1377,1W/m2oC, luas permukaan transfer panas 1336,4m2, beda temperatur air panas dan air pendingin 13,6oC (∆TLmtd), rugi tekanan desain 95,74kPa (sisi tube) dan 67,37kPa (sisi shell). Dari hasil perhitungan desain termal maka dapat digunakan penukar panas yang dimensinya mendekati data desain standar yang ada untuk tipe shell-tube. Kata Kunci : desain termal, penukar panas, RRI-50
ABSTRACT THERMAL DESIGN PRELIMINARY CALCULATION FOR HEAT EXCHANGER OF RRI-50 COOLING SYSTEM. The design of flow cooling system process RRI-50 has been made. This paper to complete the documents for the heat exchanger unit components of reactor coolant system. Two units of heat exchangers shell-tube type are selected and designed in which the capacity of each 25,1MWt. Design is done by the method of determining the heat transfer coefficient in the heat exchanger. The design, flow pressure loss is also one of the factor that must be accepted.Thermal design of heat exchanger shell-tube type is selected that includes data operating parameters and specifications.The input data used are temperature and flow rate of the RRI-50 taken from the design of process cooling system. Design result obtained are as follows : shell-tube exchanger type (1-2 pass) counter-current flow provided by six baffle plates with the heat transfer coefficient of 1377,1W/m.2oC, heat transfer area of 1336,4m2, temperature difference between hot water and coolant water is 13,6oC (∆TLmtd), the pressure drop design are 95,74kPa (tube side) and 67,37kPa (shellside). From the result of thermal design calculation then the heat exchanger can be selected in which the dimension in closed with the existing data of shell-tube type. Keywords : thermal design, heat exchanger, RRI-50
1
PENDAHULUAN Desain proses sistem pendingin merupakan bagian desain pada konseptual desain Reaktor Riset Inovatif daya 50 MWt (RRI-50).Desain penukar panas merupakan salah satu desain pendukung dalam pelaksanaan desain reaktor RRI-50. Sistem proses pendinginan reaktor riset diantaranya meliputi komponen utama yaitu kolam reaktor di dalam bejana reaktor, pompa sirkulasi sistem pendingin primer dan sekunder, penukar panas, tangki tunda N-16 dan menara pendingin (cooling tower). Pembuatan diagram aliran proses (flow-sheet) merupakan langkah desain pada tahap awal untuk menetapkan kondisi terhadap parameter proses pendingin (neraca laju aliran pendingin dan neraca panas)[1]. Pengambilan panas dari reaktor didukung oleh sistem pendingin primer dan sekunder melalui media antarmuka penukar panas. Di dalam sistem proses pendinginan, unit penukar panas untuk mentransfer panas dari aliran pendingin primer ke aliran pendingin sekunder merupakan hal yang sangat penting. Sampai saat ini penukar panas yang banyak digunakan dalam industri besar dan yang memiliki banyak data dukung referensinya adalah jenis shelltube. Penukar panas jenis ini memiliki keunggulan diantaranya mudah perawatan, luas permukaan besar, turbulensi aliran yang cukup tinggi, dan rugi tekanan aliran yang rendah[2]. Setiap penukar panas selalu di desain untuk mampu mentransfer panas secara optimal, permukaan transfer panas yang minimum dan kondisi operasi yang efektif. Dari uraian tersebut maka tujuan makalah ini adalah menghitung desain termal unit penukar panas jenis shell-tube pada sistem pendingin RRI-50 untuk memenuhi
2
ketentuan parameter operasi berdasarkan desain diagram proses pendingin RRI-50. Di dalam desain penukar panas, parameter yang harus dipenuhi adalah luas permukaan transfer panas yang diperlukan dan rugi tekanan aliran pendingin [3]. Disamping hal di atas, kriteria lain yang perlu dipertimbangkan terhadap persyaratan desain adalah kondisi operasi dari desain sistem proses pendingin yang harus dipenuhi. Dengan menggunakan metoda analitik yang tidak terlalu rumit melalui persamaan konstitutif dan persamaan empirik yang sudah baku maka hasil desain termal penukar panas dapat diperoleh. Disamping itu penggunaan fluida pendingin air pada kedua sisi penukar panas sistem pendingin reaktor menjadikan perhitungan awal desain termal ini mudah dikerjakan.
TEORI Diskripsi Penukar Panas Penukar panas shell-tube adalah salah satu tipe penukar panas, yang mana penukar panas ini merupakan jenis yang paling umum dan cocok untuk aplikasi tekanan tinggi. Penukar panas terdiri dari shell (bejana silindris) dengan bundle tube di dalamnya. Air pendingin dilewatkan melalui tube, dan air panas dari reaktor mengalir melalui shell pada arah yang berlawanan atau bersilangan untuk mentransfer panas antara dua fluida tersebut. Untuk meningkatkan efisiensi transfer panas, maka penukar panas ini dapat dilengkapi dengan baffle [4]. Air pada sisi shell mengalir menyilang (crossflow) terhadap bundel tube. Adapun pada sisi tube, aliran masuk dan keluar melalui channel head yang dilengkapi dengan partisi aliran. Secara umum sketsa ilustrasi penukar panas ditunjukkan pada Gambar 1.
Gambar 1.Sketsa penukar panas tipeshelltube
METODOLOGI Perhitungan Desain Langkah pertama yang harus dikerjakan dalam proses desain termal adalah mengidentifikasi input data operasi (laju aliran fluida, temperatur dan rugi tekanan aliran). Selanjutnya menentukan jenis atau tipe penukar panas yang dikehendaki dan kemudian melakukan perhitungan desain berdasarkan kondisi parameter operasi tersebut. Prosedur perhitungan desain penukar panas secara lengkap ditampilkan pada diagram Gambar 2. Langkah I diperoleh dari hasil desain proses sistem pendingin menggunakan ChemCad.6.4.1 [1]. Langkah II merupakan proses desain termal penukar panas yang pada dasarnya dalam perhitungan ini adalah menghitung koefisien transfer panas dari fluida panas ke fluida pendingin dan menghitung luasan transfer panas yang diperlukan. Diagram proses perhitungan ini dapat dengan mudah dikerjakan dengan pemrograman komputer. Meskipun penentuan hasil akhir desain penukar panas diputuskan pada langkah III (setelah garis
kotak) namun harus dilakukan perhitungan pada langkah II dengan benar. Catatan bahwa faktor lain di luar perhitungan desain yang juga dipertimbangkan adalah perioda cleaning-up, ruang (space) yang tersedia, syarat standar ASME safety code harus dipenuhi. Disamping itu unit komponen penukar panas senantiasa memiliki spesifikasi dan ukuran standar tertentu.
Gambar 2. Diagram perhitungan desain
3
Perssamaan-persaamaan pentinng dalam dessain penuukar panas disampaikann pada uraaian beriikut [3],
Q = M C (tin − toutt )............................. . ..........(1) ( A=
A=
Lintasaan pada aliraan sisi shell lebih koomplek, olehh karena itu digunakan d korelasi daari Kern yaang ditentukaan oleh diam meter ekkivalen, keceepatan crossfflow, faktor friksi, f juumlah aliran menyilang (crosses), ( dannfluid prroperties:
Q ..........................................(2) U o ∆TLMTD Q .......... .......... .......... .........((3) U o ∆TLMTD 1
Uo =
........( 4) ⎛ Ai ⎞ (1 / hshell ) + (1 / htube )⎜⎜ + Rd ⎝ Ao ⎠ Keteerangan Q = beban panas air, Joule/s A = luasan transffer panas, m2 M = laju aliran massa m pendinggin, kg/s C = kapasitas pan nas spes.air, Joule/ J kg.oC o t= teemperatur airr, C o Uo= koefisien traansfer panas, W/m W 2 C o ∆TLmtd =beda temp peratur air paanas-dingin, C L h = koefisien traansfer panas bagian sisi tuube o dan shell, W/m2 C o Rd = koefisien fouling, m2 C/W d termal penukar pannas, Dalam desain pengggunaan korrelasi rugi tekanan aliran merrupakan salah h satu faktorr penentu. Rugi R tekaanan dapat disebabkan oleh beberaapa benttuk hilangny ya energi kaarena friksi dan d peruubahan luas penampang aliran. Rugi R tekaanan pada sisi s tube diihitung denggan mennggunakan persamaan p (55) seperti paada aliraan pipa beriku ut [5] : ∆P = (4 f
L ρV 2 ) .......... .......... ......( 5 ) D 2 gc
Keteerangan : f= faktor fa friksi ρ= densiitas air,kg/m3 L= panj.kanal, p m D= diam meter kanal, m V=kkecepatan.alirran,m/s..gc=koonst.percepataan gravvitasi, m/s2
4
Keterangan : K G = kecepatan massa air, kgg/m2.s N = jumlah bafffle D = diameter shell, s m Deq = diameterr ekivalen linttasan aliran, m H HASIL DESA AIN DAN PE EMBAHASAN N Data speesifikasi desaain termal pennukar paanas pada siistem pendinngin RRI-50 telah diiperoleh yaituudata temperratur, tekanann dan laaju aliran masssa pendinginn, kapasitas beban b paanas dan dim mensi penukaar panas. Addapun daataspesifikasii desain termal terrsebut diisampaikan pada Tabell 1, pada tabel teersebut tamppak bahwa penukar panas p diidesain untuuk kapasitas 25,1MWt yang diikutip dari data hasil dessain proses sistem peendingin RRI-50 R mennggunakan paket C ChemCad6.4.1 1 [1]. Dalam m desain terssebut, teerdapat 2 unit penukaar panas untuk u m mengatasi beban panas 50 MWt dari reeaktor R RRI-50. Penukarr panas ini tippe shell-tubee (1-2 paass) yang dikoneksikann dengan sistem peendingin prim mer dan sisttem sekunderr, air peendingin sekkunder dari menara m penddingin m mengalir padaa sisi tube dan d air penddingin prrimer tanpa mineral m (dem mineralized water) w diilewatkan padda sisi shell. Dengan D penetapan inni maka alirann air yang mempunyai m tinngkat kaandungan imppurities besarr dilewatkan pada siisi tube, hal h ini dim maksudkan untuk u m memudahkan proses peerawatan beerkala clleaning-up. Selanjutnya, S arah aliran yang diigunakan anttara sisi sheell dan sisi tube addalah berlaw wanan (countter current flow). f
Tabel 1. Spesifikasi Desain Penukar Panas Parameter
shell
tipe fluida sirkulasi
Tube Air pendingin 33,05 42,00 665,6
Air demineral
beban panas (MJoule/s) temperatur air panas masuk (oC) temperatur air panas keluar (oC) temperatur pendingin masuk (oC) temperatur pendingin keluar (oC) laju aliran air panas (kg/s) laju aliran pendingin (kg/s) ∆TLmtd (oC) Faktorfouling (m2 oC/W) angka Reynolds (-) Uoverall (W/m2C) shell- tube pass (-)
25,03 57,85 44,60 451,0 -
13,6 0,00034 37077
[5] 38668
1377,1 1
2
diameter (mm)
1495,6
20 ID / 22 OD
Pitchtube / layout jumlah tube (-) jumlah baffle (-) Jarak baffle (mm) luas transfer panas (m2) panjang tube(m) ∆P (kPa) * = maximum allowable
1,25OD / square pitch 1342 6 1200,0 1336,4 7,20 95,74 67,37
∆TLMTD yang dihitung dari data desain proses menunjukkan bahwa angka ini cukup marjin untuk desain awal penukar panas dalam mengatasi degradasi alat selama dioperasikan. Semakin kecil angka ∆TLMTD desain maka semakin baik kemampuan alat penukar panas karena perbedaan yang kecil antara temperatur air panas dan air pendinginnya. Gambar 3 menunjukkan kurva temperatur yang dipakai untuk menghitung temperatur ∆TLMTD.
Gambar 3 Kurva temperatur air panas dan air pendingin Pada bagian sisi shell penukar panas terdapat 6 baffles penyekat, hal ini
5
Catatan
shell-tube
ID,OD:in,outdia. Ref std=1524 mm [6] Ref std = 1624 [6]
103 *[6]
dimaksudkan untuk memperoleh efektivitas yang tinggi dalam mengambil panas dari sisi shell ke sisi tube. Dari tabel 1 terlihat pula bahwa Jarak baffle (baffle spacing) 1200mm dan diameter Shell 1495,6mm, ini berarti bahwa aliran pada sisi shell cukup baik karena mengalir secara menyilang (across) tegak lurus terhadap bundel tube. Penggunaan baffle tidak menghasilkan ∆P yang cukup tinggi, hal ini tampak pada angka ∆P sisi shell hasil perhitungan 67,37kPa yang masih memenuhi syarat karena kurang dari 103kPa [6]. Penukar panas ini ditentukan dengan bentuk konfigurasi tube square pitch, oleh karena itu aliran di sisi shell memiliki ∆P lebih rendah dibanding bentuk triangular pitch. Perlu disampaikan di sini bahwa triangularpitch dapat memberikan turbulensi aliran yang tinggi namun berakibat ∆P yang tinggi pula. Sementara itu hasil perhitungan koefisien transfer panas menunjukkan 1377,1W/m2.oC, angka ini cukup memenuhi kriteria karena berada pada rentang besaran
nilai dalam desain penukar panas yang menggunakan media air fasa tunggal pada sisi shell dan sisi tube (800 – 1500 W/m2 oC) [7]. Hasil perhitungan desain untuk diameter shell = 1495,6mm dan jumlah tube =1342, sehingga dari referensi dimensi standar penukar panas shell-tubedapat digunakan data dimensi yang mendekati angka tersebut yaitu untuk diameter shell adalah 1524mm dan jumlah tube 1624[6].
2.
3. 4.
KESIMPULAN Desain awal penukar panas untuk sistem pendingin RRI-50 telah dapat diselesaikan, penggunaan baffle penyekat masih dapat mengatasi rugi tekanan karena masih dibawah kriteria yang disyaratkan. Begitu pula transfer panasnya dari perhitungan telah sesuai dengan kriteria rentang koefisien transfer panas untuk pendingin air. Dari hasil perhitungan maka dapat digunakan penukar panas yang dimensinya mendekati data desain standar yang ada untuk tipe shell-tube.
5. 6.
7.
Sekunder RRI”, LAPTEK 2013, PTRKN Batan. BELL K.J, “Heat Exchanger Design Hand-Book (HEDH)”, Bab 3, 4 dan 5, HemispherePublish Corp, Washington DC, 1983. KERN DQ, “Process Heat Transfer”, International Student Edition, Mc. Graw Hill Book Co, New York, 1965. SUKMANTO D, Studi Leakage Baffle Pada Penukar Panas Shell-Tube, Presentasi Ilmiah Hasil Studi DoktorMaster BATAN, Jakarta, 9-11/12/1997. COULSON JM RICHARDSON JF “Chemical Engineering Design”, p.548553, Pergamon Press, New York. M. KONTOPOULOU,” Heat Transfer Analysis and Heat Exchanger Design Part 3”, Department of Chemical Engineering, Queen's University Kingston, Ontario Canada. Site Up-date: January 3, 2014, http://www.chemeng.queensu.ca/courses/ CHEE331+ 332+333 EngineeringPage,”Typical Overall Heat Transfer CoefficientU-Values”, 09/ 2014http://www.engineeringpage.com/te chnology/thermal/transfer.html
DAFTAR PUSTAKA 1. SUKMANTO D. dkk,” Desain Parameter Proses Operasi Sistem Pendingin
6
***
