SIMULASI PENENTUAN POSISI DAN DAYA KIPAS UNTUK KESERAGAMAN POLA ALIRAN UDARA PANAS PADA PENGERING ERK TIPE RAK DENGAN METODE COMPUTATIONAL FLUID DYNAMICS (CFD)
ALFREDO
DEPARTEMEN TEKNIK MESIN DAN BIOSISTEM FAKULTAS TEKNOLOGI PERTANIAN INSTITUT PERTANIAN BOGOR BOGOR 2014
PERNYATAAN MENGENAI SKRIPSI DAN SUMBER INFORMASI SERTA PELIMPAHAN HAK CIPTA Dengan ini saya menyatakan bahwa skripsi berjudul Simulasi Penentuan Posisi dan Daya Kipas untuk Keseragaman Pola Aliran Udara Panas pada Pengering ERK Tipe Rak dengan Metode Computational Fluid Dynamics (CFD) adalah benar karya saya dengan arahan dari dosen pembimbing dan belum diajukan dalam bentuk apa pun kepada perguruan tinggi mana pun. Sumber informasi yang berasal atau dikutip dari karya yang diterbitkan maupun tidak diterbitkan dari penulis lain telah disebutkan dalam teks dan dicantumkan dalam Daftar Pustaka di bagian akhir skripsi ini. Dengan ini saya melimpahkan hak cipta dari karya tulis saya kepada Institut Pertanian Bogor. Bogor, Februari 2014 Alfredo NIM F14090140
ABSTRAK ALFREDO. Simulasi Penentuan Posisi dan Daya Kipas untuk Keseragaman Pola Aliran Udara Panas pada Pengering ERK Tipe Rak dengan Metode Computational Fluid Dynamics (CFD). Dibimbing oleh DYAH WULANDANI. Pengering Efek Rumah Kaca (ERK) adalah bangunan pengering berdinding transparan, di dalamnya terdapat plat absorber sebagai pengumpul panas dan wadah produk (rak atau bak) serta kipas untuk mengeluarkan uap air hasil pengeringan atau sebagai penyeragaman suhu di dalam pengering. Sumber energi pada pengering ERK diperoleh dari surya dan pembakaran biomassa. Masalah ketidakmerataan aliran udara panas di dalam ruang pengering, masih menjadi salah satu yang selalu terjadi. Pemecahan masalah ketidakmerataan aliran udara panas dapat dilakukan dengan metode Computational Fluid Dynamics (CFD). Penelitian ini bertujuan untuk mensimulasikan aliran udara panas dalam ruang pengering, lalu memvalidasikan hasil tersebut dan memodifikasi posisi kipas untuk sebaran udara panas yang lebih seragam. Metode simulasi dengan CFD dapat digunakan untuk mengetahui pola aliran udara panas dalam pengering. Validasi suhu dilakukan pada pukul 08.00, 12.00 dan 16.00 dengan kondisi tanpa beban, menunjukkan hasil simulasi yang sesuai dengan kondisi nyata. Modifikasi posisi penempatan kipas menggunakan 5 skenario. Pada simulasi modifikasi dengan kondisi memakai beban pada pukul 12.00 hari pertama pengeringan. Menunjukkan bahwa modifikasi 5, yaitu posisi kipas pada outlet atas sebelah kiri dan kanan pengering memberikan hasil distribusi suhu udara pengering yang paling seragam. Kata kunci: Aliran udara panas, CFD, ERK
ABSTRACT ALFREDO. Simulation of Determination Fan Position for Uniformity Pattern Heat Air Flow at Greenhouse Effect (GHE) Solar Dryer Rack Type with Computational Fluid Dynamics (CFD). Supervised by DYAH WULANDANI. Greenhouse Effect (GHE) Solar Dryer is a transparent wall structure, consists of an absorber plate as solar collector, product holders (rack, tray or batch) and fans to discharge vapor evaporated from product or to uniform heat air flow in dryer. The GHE solar dryer uses thermal energy source from the sun and biomass combustion. Un-uniformity of heat air flow in dryer chamber is still to be a big problem. Un-uniformity of heat air flow problem can be solved by using Computational Fluid Dynamics (CFD) method. The objective of the research is to simulate heat air flow in drying chamber then validate these results and modify position of fan for heat air distribution more uniform. Simulation of Computational Fluid Dynamics can be used to determine the heat air flow pattern in dryer. Validation of temperature performed at eight, twelve and sixteen o’clock with no load condition, have shown similar trend between simulation result and experiment of drying temperature distribution. Modification of the position of the fan placement using 5 scenario. Simulation performed with load condition at
twelve o’clock on the first day of drying. Simulation results show that modification 5, the position of the fan on the outlet on the left and the right of drying wall showed the best uniform of drying temperature distribution. Keywords: CFD, GHE, heat air flow
SIMULASI PENENTUAN POSISI DAN DAYA KIPAS UNTUK KESERAGAMAN POLA ALIRAN UDARA PANAS PADA PENGERING ERK TIPE RAK DENGAN METODE COMPUTATIONAL FLUID DYNAMICS (CFD)
ALFREDO
Skripsi sebagai salah satu syarat untuk memperoleh gelar Sarjana Teknologi Pertanian pada Departemen Teknik Mesin dan Biosistem
DEPARTEMEN TEKNIK MESIN DAN BIOSISTEM FAKULTAS TEKNOLOGI PERTANIAN INSTITUT PERTANIAN BOGOR BOGOR 2014
Judul Skripsi : Simulasi Penentuan Posisi dan Daya Kipas untuk Keseragaman Pola Aliran Udara Panas pada Pengering ERK Tipe Rak dengan Metode Computational Fluid Dynamics (CFD) Nama : Alfredo NIM : F14090140
Disetujui oleh
Dr Ir Dyah Wulandani, M.Si Pembimbing Akademik
Diketahui oleh
Dr Ir Desrial, M. Eng Ketua Departemen
Tanggal Lulus:
PRAKATA Puji dan syukur penulis panjatkan kepada Tuhan Yesus atas segala karuniaNya sehingga karya ilmiah ini berhasil diselesaikan. Tema yang dipilih dalam penelitian yang dilaksanakan sejak bulan Februari 2013 ini ialah pengering ERK, dengan judul Simulasi Penentuan Posisi dan Daya Kipas untuk Keseragaman Pola Aliran Udara Panas pada Pengering ERK Tipe Rak dengan Metode Computational Fluid Dynamics (CFD) Terima kasih penulis ucapkan kepada Ibu Dr Ir Dyah Wulandani, M.Si selaku pembimbing, yang telah memberikan bimbingan, arahan, motivasi kepada penulis, Dr Ir Emmy Darmawati, M.Si dan Dr Muhamad Yulianto ST, MT sebagai dosen penguji yang telah memberikan saran untuk penelitian ini. Ungkapan terima kasih juga disampaikan kepada mama, abang, dan adik atas segala dukungan dan doanya. Terima kasih juga penulis ucapkan kepada teman-teman sebimbingan (Elsamila, Stephani, Adit, dan Sueritah), TEP Orion 46, Komlit 46 (Eka, Victor, Jorex, Hera, dll), Komlit 47, 48, 49 atas semangat, doa dan dukungan yang telah diberikan. Kepada pak Darma, pak Harto dan teknisi departemen lainnya atas bantuannya. Semoga karya ilmiah ini bermanfaat.
Bogor, Februari 2014 Alfredo
DAFTAR ISI DAFTAR TABEL
xi
DAFTAR GAMBAR
xi
DAFTAR LAMPIRAN
xii
DAFTAR SIMBOL
xiii
PENDAHULUAN
1
Latar Belakang
1
Perumusan Masalah
2
Tujuan Penelitian
3
Manfaat Penelitian
3
Ruang Lingkup Penelitian
3
TINJAUAN PUSTAKA
3
Pengeringan
3
Aplikasi CFD pada Pengeringan Efek Rumah Kaca
4
CFD
5
Computational Fluid Dynamics (CFD)
5
Proses Simulasi CFD
6
METODE
6
Waktu dan Tempat
6
Bahan
7
Alat
7
Prosedur Penelitian
8
Persiapan Alat
8
Pengukuran dan Perhitungan
8
Pembuatan Model Geometri dan Pendefinisian Kondisi Batas
9
Percobaan Pengeringan Tanpa Beban (Percobaan 1 dan 2)
9
Simulasi dan Validasi Tanpa Beban
9
Simulasi dan Validasi Dengan Beban
9
Simulasi Modifikasi Penempatan Posisi Kipas dan Uji Keseragaman
10
Validasi Model
10
HASIL DAN PEMBAHASAN Tahapan Awal Simulasi
11 11
Tahapan Pra-pemrosesan (Preprocessing)
11
Asumsi, Kondisi Batas dan Perhitungan Parameter Kondisi Batas
13
Asumsi
13
Kondisi Batas
13
Perhitungan Parameter Kondisi Batas
14
Penyebaran Suhu Rak Hasil Pengukuran
16
Profil Sebaran Suhu Rak Bernomor Ganjil (Pengukuran 1)
16
Iradiasi Surya (pengukuran 1)
19
Profil Sebaran Suhu Rak Bernomor Genap (Pengukuran 2)
20
Iradiasi Surya (pengukuran 2)
22
Kontur Suhu Ruang Pengering ERK Hasil Simulasi CFD pada Kondisi Tanpa Beban 23 Suhu Ruang Pengering pada Pukul 08.00
23
Suhu Ruang Pengering pada Pukul 12.00
24
Suhu Ruang Pengering pada Pukul 16.00
25
Kontur Suhu Ruang Pengering ERK Hasil Simulasi CFD pada Kondisi dengan Beban 27 Validasi Model
28
Modifikasi Posisi Kipas Menggunakan Simulasi CFD
30
Modifikasi 1 Terdapat Dua Kipas di Tiap Ruang Pengering
30
Modifikasi 2 Terdapat Empat Kipas di Tiap Ruang Pengering
32
Modifikasi 3 Terdapat Satu Kipas di Tiap Ruang Pengering
33
Modifikasi 4 Kipas di Outlet Tiap Ruang Pengering
36
Modifikasi 5 Kipas di Outlet Atas Tiap Ruang Pengering
36
SIMPULAN DAN SARAN
37
Simpulan
37
Saran
38
DAFTAR PUSTAKA
38
LAMPIRAN
40
RIWAYAT HIDUP
75
DAFTAR TABEL 1 2 3 4 5 6 7 8 9 10 11 12
Tingkat keseragaman suhu rak bernomor ganjil Tingkat keseragaman suhu rak bernomor genap Distribusi suhu rak pada simulasi tanpa beban pukul 08.00 Distribusi suhu rak pada simulasi tanpa beban pukul 12.00 Distribusi suhu rak pada simulasi tanpa beban pukul 16.00 Distribusi suhu rak pada simulasi dengan beban Distribusi suhu rak hasil simulasi CFD pada modifikasi 1 Distribusi suhu rak hasil simulasi CFD pada modifikasi 2 Distribusi suhu rak pada modifikasi 3 ΔP 4.905, 50 dan 100 Koefisien relasi antara pressure jump (ΔP) dengan CV Distribusi suhu rak hasil simulasi CFD pada modifikasi 4 Distribusi suhu rak hasil simulasi CFD pada modifikasi 5
19 22 24 25 26 28 31 33 35 35 36 37
DAFTAR GAMBAR 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27
Hasil simulasi CFD pada pengering ERK (Wulandani, 2005) Alat pengering Efek Rumah Kaca (ERK) hybrid Diagram alir penelitian Tampilan isometri geometri 3D pengering ERK Tampilan depan geometri 3D pengering ERK Tampilan hasil mesh geometri pengering ERK Diagram alir simulasi Posisi dan penomoran rak Grafik suhu rak I D Grafik suhu rak II B Grafik suhu rak III D Grafik suhu rak IV B Profil iradiasi surya pengukuran 1 Grafik suhu rak I B Grafik suhu rak II D Grafik suhu rak III B Grafik suhu rak IV D Profil iradiasi surya pengukuran 2 Kontur suhu simulasi tanpa beban pukul 08.00 (Z = 0.5325 m) Kontur suhu simulasi tanpa beban pukul 08.00 pada sumbu X Kontur suhu simulasi tanpa beban pukul 12.00 (Z = 0.5325 m) Kontur suhu simulasi tanpa beban pukul 12.00 pada sumbu X Kontur suhu simulasi tanpa beban pukul 16.00 (Z = 0.5325 m) Kontur suhu simulasi tanpa beban pukul 16.00 pada sumbu X Kontur suhu simulasi dengan beban (Z = 0.5325 m) Kontur suhu simulasi dengan beban pada sumbu X Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi (pukul 08.00) 28 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi (pukul 12.00)
4 7 11 12 12 13 16 17 17 18 18 18 20 20 21 21 21 23 24 24 25 25 26 26 27 27 28 28
29 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi (pukul 16.00) 30 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi menggunakan beban 31 Kontur suhu pada modifikasi 1 v = 0.5 m/s 32 Kontur suhu pada modifikasi 1 v = 1.5 m/s 33 Kontur suhu pada modifikasi 2 v = 0.5 m/s 34 Kontur suhu pada modifikasi 1 v = 1.5 m/s 35 Kontur suhu pada modifikasi 3 porous jump 50 36 Kontur suhu pada modifikasi 3 porous jump 100 37 Kontur suhu pada modifikasi 3 porous jump 4.905 38 Grafik hubungan pressure jump dengan CV 39 Kontur suhu pada modifikasi 4 40 Kontur suhu pada modifikasi 5
29 29 31 31 32 32 34 34 34 35 36 37
DAFTAR LAMPIRAN 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
Perhitungan data input CFD pada pukul 08.00 tanpa beban Tabel validasi suhu pada pukul 08.00 kondisi tanpa beban Tabel validasi suhu pada pukul 12.00 kondisi tanpa beban Tabel validasi suhu pada pukul 16.00 kondisi tanpa beban Tabel validasi suhu pada pukul 12.00 kondisi dengan beban Data intensitas cahaya matahari Data input CFD (Boundary condition) pukul 08.00 kondisi tanpa beban Data input CFD (Boundary condition) pukul 12.00 kondisi tanpa beban Data input CFD (Boundary condition) pukul 16.00 kondisi tanpa beban Data input CFD (Boundary condition) pukul 12.00 kondisi dengan beban Tabel ΔP, efisiensi, daya kipas dan keseragaman suhu modifikasi 3 Data pengukuran suhu (oC) rak pengukuran 1 Data pengukuran suhu (oC) rak pengukuran 1 (lanjutan) Distribusi suhu rak pada pengukuran 1 Data pengukuran suhu (oC) rak pengukuran 2 Data pengukuran suhu (oC) rak pengukuran 2 (lanjutan) Distribusi suhu rak pada pengukuran 2 Data pengukuran suhu (oC) atap, dinding, HE dan lantai pengukuran 1 Data pengukuran suhu (oC) atap, dinding, HE dan lantai pengukuran 2 Tingkat keseragaman suhu rak bernomor ganjil selama pengukuran Tingkat keseragaman suhu rak bernomor genap selama pengukuran Tahapan Penggunaan ANSYS Fluent untuk Simulasi
DAFTAR SIMBOL Re ρ v D μ
Bilangan Reynold Massa jenis (kg/m3) Kecepatan (m/dtk) Diameter spesifik (m) Viskositas dinamik (kg/ms)
40 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 65 66 66 67
h Nu k l q/A I τα Td T∞ α Dp ε C2 k ΔP D ΔTln Gr g β T x ν Ra Pr η P ṁ
Koefisien pindah panas konveksi (W/m2K) Bilangan Nusselt Konduktivitas panas (W/mK) Ketebalan dinding (m) Fluks panas (W/m2) Iradiasi surya (W/m2) Hasil kali koefisien tembus cahaya penutup transparan Suhu dinding (oC) Suhu lingkungan (oC) Permeabilitas permukaan (m2) Diameter produk (m) Porositas tumpukan produk Koefisien porous jump (1/m) Koefisien kehilangan pada penukar panas Penurunan tekanan (Pa) Diameter pipa HE (m) Bilangan Grashoff Gaya gravitasi (m/s2) Koefisien (1/oC) Suhu (K) Panjang pipa HE (m) Viskositas kinematik (m2/s) Bilangan Rayleigh Bilangan Prandtl Efisiensi kipas Daya kipas (Watt) Laju aliran massa (kg/dtk)
PENDAHULUAN Latar Belakang Pengeringan hasil pertanian secara khusus adalah pengeluaran air dari suatu hasil pertanian sampai tercapai kadar air keseimbangan dengan udara lingkungan atau sampai tingkat kadar air dimana jamur, serangga dan enzim yang bersifat merusak tidak dapat aktif lagi. Pada umumnya keadaan ini dicapai pada tingkat kadar air 12-14% basis basah (Henderson dan Perry, 1976). Bahan pangan yang dihasilkan dari produk-produk pertanian pada umumnya mengandung kadar air. Kadar air jika tidak dihilangkan dapat mempengaruhi kondisi fisik bahan pangan. Proses pengeringan akan mengakibatkan produk yang dikeringkan mengalami perubahan warna, tekstur, rasa dan aroma. Faktor-faktor yang mempengaruhi pengeringan terdiri dari faktor udara pengering dan sifat bahan. Faktor yang berhubungan dengan udara pengering adalah suhu, kecepatan aliran udara pengering, dan kelembaban udara, sedangkan faktor yang berhubungan dengan sifat bahan yaitu ukuran bahan, kadar air awal, dan tekanan parsial dalam bahan (Fellow, 2001). Pada umumnya, pengeringan produk pertanian masih dilakukan secara tradisional yaitu pengeringan dengan cara penjemuran di bawah sinar matahari secara langsung. Usaha untuk melakukan proses pengeringan dapat pula dilakukan dengan menggunakan pengering efek rumah kaca (ERK), yaitu bangunan berbentuk segi empat, silinder atau kerucut terpancung dan berdinding transparan. Dinding plastik transparan ini berfungsi sebagai perangkap gelombang panjang yang dipancarkan oleh sinar matahari, sehingga terjadi akumulasi panas yang dimanfaatkan untuk pengeringan produk pertanian. Salah satu keuntungan dengan menggunakan pengering ERK adalah produk lebih bersih dan tidak terkontaminasi dengan bahan yang tidak dikehendaki. Dalam proses pengeringan masalah distribusi aliran panas yang belum merata menyebabkan menurunnya kinerja alat pengering, nampan-nampan pada pengering tipe rak dapat menyebabkan terjadinya perbedaan suhu antara bagian atas dan bagian bawah rak pengering. Wulandani dan Nelwan (2009) melaporkan bahwa, tingkat keseragaman suhu udara pengering pada alat pengering tipe rak yang sama digunakan dalam penelitian ini saat kondisi tanpa beban yang dinyatakan dengan nilai standar deviasi, pada posisi vertikal antara rak atas dan rak bawah adalah sebesar 2.4oC, 1.2oC dan 2.7oC pada masing-masing kolom rak. Sedangkan pada posisi horisontal, diperoleh standar deviasi sebesar 2.3oC dan 3.1oC, masing-masing pada rak atas dan rak bawah. Hartini (2010) dalam penelitiannya menunjukkan bahwa terjadi perbedaan suhu antara rak bagian atas dengan rak bagian bawah sebesar 4oC pada kondisi tanpa beban. Sedangkan Aritesty (2013) dalam penelitiannya menunjukkan perbedaan suhu antara rak bagian atas dengan rak bagian bawah pada pengering ERK saat pukul 12.00 dengan kondisi tanpa beban adalah sekitar 4oC pada hari pertama dan sekitar 22.7oC pada pengamatan hari kedua. Sitanggang (2010) melakukan pengujian dan simulasi suhu, kecepatan udara pada mesin pengering produk pertanian energi surya tipe rak bertingkat menggunakan program CFD. Kriteria mesh yang digunakan berupa tipe pave dengan elemen quad dengan ukuran yang disesuaikan
2 kebutuhan. Hasilnya menunjukkan bahwa model rak mesin pengering yang disusun secara zig-zag pada analisa simulasi modifikasi jauh lebih baik dibandingkan dengan model rak yang tersusun bertingkat yang diuji secara langsung di lapangan. Yunus dan Kayiem (2013) mensimulasikan pengering ERK tipe rak bertingkat dengan CFD. Kriteria mesh yang digunakan berupa tipe pave dengan elemen tri dengan ukuran interval sebesar 1.5 – 3. Model yang digunakan berupa P-1 radiation model, aliran steady dan kompresibel. Hasil simulasi riset tersebut menunjukkan error suhu sebesar 11.3% - 31.6%. Masalah ketidakseragaman kadar air produk hasil pengeringan juga merupakan akibat dari ketidakmerataan aliran panas di dalam ruang pengering, khususnya tipe rak. Untuk memecahkan masalah tersebut perlu dilakukan suatu analisis sifat dan pola aliran serta distribusi udara di dalam ruang pengering. Distribusi dan pola aliran udara diduga ditentukan oleh geometri ruang pengering, lokasi penempatan inlet dan outlet, penempatan dan kapasitas daya kipas, susunan rak, dan lokasi sumber panas di dalam bangunan pengering. Dengan mensimulasikan distribusi dan pola aliran udara ini, akan dapat ditentukan desain bangunan beserta penempatan parameter-parameter penentu diatas secara tepat sehingga tujuan keseragaman mutu produk terpenuhi. Pemecahan analisis aliran udara dapat dilakukan dengan metode CFD (Computational Fluid Dynamics). Metode CFD menggunakan analisis numerik yaitu kontrol volume sebagai elemen dari integrasi persamaan-persamaan yang terdiri dari persamaan keseimbangan massa, momentum dan energi (Versteeg dan Malalasekera, 1995). Dengan demikian penyelesaian persamaan untuk benda dua dimensi atau tiga dimensi lebih cepat dan dapat dilakukan secara simultan. Melalui CFD dapat ditentukan desain dengan penempatan kipas, sistem pemanas, sistem penyaluran udara dan wadah produk dan penghematan energi akan tercapai dan efisiensi meningkat.
Perumusan Masalah Performansi suatu pengering dapat dipengaruhi oleh dua hal. Pertama suhu udara pengering di dalam pengering ERK sangat berfluktuasi karena sangat dipengaruhi oleh keberadaan surya. Iradiasi surya sifatnya selalu berubah dan besar iradiasinya sangat dipengaruhi oleh waktu, lokasi dan musim. Kedua adalah distribusi aliran panas dalam ruang pengering yang belum merata, khususnya pada pengering tipe rak. Nampan-nampan pada tipe rak ini dapat menyebabkan distribusi udara yang kurang baik dan menurunkan kinerja pengeringan (Wulandani, 2005). Penelitian mengenai aliran udara panas di dalam pengering ERK dengan metode CFD digunakan untuk mendapatkan pola aliran udara serta distribusi suhu. Distribusi dan pola aliran udara diduga ditentukan oleh geometri ruang pengering, lokasi penempatan kipas. Dengan mensimulasikan distribusi dan pola aliran udara ini, akan dapat ditentukan posisi kipas sehingga tujuan keseragaman udara panas dapat terpenuhi.
3 Tujuan Penelitian Tujuan penelitian ini adalah untuk mensimulasikan pola aliran udara panas yang terjadi dalam ruang pengering. Memvalidasikan antara hasil pengukuran dengan hasil simulasi menggunakan CFD dan melakukan simulasi penempatan posisi kipas untuk mendapatkan aliran udara panas yang lebih seragam dalam ruang pengering.
Manfaat Penelitian Hasil dari penelitian ini diharapkan dapat meningkatkan performa pengeringan pada pengering ERK tipe rak terutama untuk kegiatan pengeringan produk-produk pertanian; kedua, hasil penelitian ini bermanfaat bagi pemakai alat pengering ERK, terutama bagi para petani maupun industri dan ketiga adalah hasil disain dari simulasi dinamika fluida dalam pengering ERK ini dapat dimanfaatkan sebagai dasar kerangka acuan dalam menentukan disain optimal suatu pengering produk-produk pertanian, untuk skala kecil maupun besar serta untuk pemanfaatan lainnya.
Ruang Lingkup Penelitian Ruang lingkup dari penelitian ini adalah mensimulasikan pola aliran udara panas dalam ruang pengering ERK tipe rak pada kondisi tanpa beban saat pukul 08.00, 12.00 dan pukul 16.00. Dasar dari penentuan pukul tersebut adalah ingin diketahuinya kondisi dan hasil simulasi CFD pada posisi dan intensitas matahari dari kondisi yang berbeda. Simulasi pada kondisi dengan beban saat pukul 12.00 hari pertama pengeringan. Pemilihan pukul 12.00 hari pertama pengeringan pada kondisi menggunakan beban adalah karena pada kondisi tersebut, bahan yang dikeringkan masih berada pada tahap dimana kadar air cukup tinggi dan intensitas matahari tinggi. Memvalidasikan pola udara panas hasil simulasi pada kondisi tanpa beban dan kondisi menggunakan beban. Mensimulasikan penempatan posisi dan daya kipas dengan kondisi memakai beban pada pukul 12.00 pada hari pertama pengeringan.
TINJAUAN PUSTAKA Pengeringan Pengeringan hasil pertanian secara khusus adalah pengeluaran air dari suatu hasil pertanian sampai tercapai kadar air keseimbangan dengan udara lingkungan atau sampai tingkat kadar air dimana jamur, serangga dan enzim yang bersifat merusak tidak dapat aktif lagi. Pada umumnya keadaan ini dicapai pada tingkat kadar air 12-14% basis basah (Henderson dan Perry, 1976).
4 Hall (1957) dan Broker (1974) menyatakan, pengeringan adalah suatu proses pengurangan kadar air sampai pada batas tertentu sehingga memperlambat laju kerusakan bahan pertanian akibat aktivitas biologis dan kimia sebelum bahan diolah atau dimanfaatkan. Proses pengeringan merupakan proses pengeluaran air dari bahan pertanian menuju kadar air keseimbangan dengan udara di sekelilingnya atau pada tingkat dimana mutu bahan pertanian dapat dijaga dari serangan kapang, aktivitas serangga, dan enzim (Henderson dan Perry, 1976). Menurut Brooker et al., (1974), salah satu parameter yang mempengaruhi waktu pengeringan yang dibutuhkan adalah suhu udara pengering. Laju penguapan air bahan dalam pengering sangat ditentukan oleh kenaikan suhu. Bila suhu pengering dinaikkan maka panas yang dibutuhkan untuk penguapan air bahan menjadi berkurang. Suhu udara pengering berpengaruh terhadap lama pengeringan dan kualitas bahan hasil pengeringan. Makin tinggi suhu udara pengering maka proses pengeringan makin singkat. Biaya pengeringan dapat ditekan pada kapasitas yang besar jika digunakan pada suhu tinggi, selama suhu tersebut tidak sampai merusak bahan.
Aplikasi CFD Pada Pengeringan Efek Rumah Kaca Wulandani, et al (2003) menggunakan metode CFD dalam menganalisis distribusi suhu dan kecepatan aliran udara dalam ruang pengering ERK. Dalam riset ini dilakukan simulasi posisi inlet dan outlet pada bangunan pengering dan perubahan ukuran model menjadi prototipe bangunan pengering dengan memakai dua skenario simulasi. Hasil riset tersebut menunjukkan bahwa model simulasi dapat menggambarkan kondisi nyata sebesar 89%. Wulandani (2005) telah menggunakan teknik CFD untuk mensimulasi udara pengering pada pengeringan ERK tipe rak. Dalam riset ini dilakukan analisis distribusi udara yang mencakup kecepatan, suhu dan RH serta dilanjutkan dengan melakukan validasi model tersebut terhadap hasil percobaan. Analisis ini penting untuk mengoptimasikan bentuk saluran udara yang harus didisain untuk menyeragamkan aliran udara pada pengering, sehingga diperoleh keseragaman kadar air yang berarti juga keseragaman kualitas biji. Gambar 1 menunjukkan salah satu hasil simulasi CFD berupa sebaran suhu pada pengering tipe rak yang dilakukan Wulandani (2005).
Kipas tengah Kipas atas Inlet Penukar panas Kipas bawah
Outlet
Rak
Gambar 1 Hasil simulasi CFD pada pengering ERK (Wulandani, 2005)
5 Gambar 1 terdiri dari dua buah inlet masing-masing berukuran 0.1 m x 1 m pada ketinggian 1.4 m. Dua buah outlet masing-masing berukuran 0.2 m x 0.8 m pada ketinggian 0.8 m pada dinding yang berseberangan dengan inlet. Tiga buah kipas dengan diameter masing-masing 0.2 m. Kipas bawah dengan daya 100 W, kipas tengah berdaya 40 W dan kipas atas dengan daya 100 W. Hasil simulasi CFD pada Gambar 1 menghasilkan nilai rata-rata suhu pada seluruh rak sebesar 45.4oC dengan nilai ragam sebesar 1.6 oC. Teknik CFD dalam analisis distribusi aliran udara juga digunakan oleh Nugraha (2005) dalam mensimulasikan pola aliran udara, suhu dalam ruang pengering ERK tipe rak pada proses pengeringan mahkota dewa. Hasil riset menunjukkan error rata-rata untuk validasi suhu sebesar 5.10%. Widodo (2009) dalam risetnya juga menggunakan metode CFD dalam mengkaji pola sebaran aliran udara panas pada pengering ERK hibrid tipe rak berputar. Widodo menggunakan enam skenario alat pengering berdasarkan bentuk geometri ruang pengering dan posisi inlet outlet untuk mendapatkan model pengering yang optimum.
CFD Computational Fluid Dynamics (CFD) Ditinjau dari istilah, Computational Fluid Dynamics (CFD) terdiri dari dua kata : • Computational : Segala sesuatu yang berhubungan dengan matematika dan metode numerik atau komputasi. • Fluid Dynamics : Dinamika dari segala sesuatu yang mengalir Ditinjau dari istilah diatas, CFD berarti suatu teknologi komputasi yang memungkinkan untuk mempelajari dinamika dari benda-benda atau zat-zat yang mengalir. Secara definisi, CFD adalah ilmu yang mempelajari cara memprediksi aliran fluida, perpindahan panas, reaksi kimia, dan fenomena lainnya dengan menyelesaikan persamaan-persamaan matematika (model matematika) (Tuakia, 2008). Terdapat 3 keuntungan penggunaan CFD, antara lain : 1. Pemahaman mendalam (Insight) Analisis CFD mampu mendesain sebuah sistem atau alat yang sulit untuk dibuat prototype-nya atau sulit untuk dilakukan pengujian. Analisis ini mampu masuk secara virtual ke dalam alat atau sistem yang dirancang. 2. Prediksi menyeluruh (Foresight) Computational Fluid Dynamics (CFD) adalah alat untuk memprediksi yang akan terjadi pada alat atau sistem yang didesain dengan satu atau lebih kondisi batas dan dapat segera menentukan desain optimal. 3. Efisiensi waktu dan biaya (Efficiency)
6 Prediksi menyeluruh (Foresight) yang diperoleh CFD mampu membantu untuk mendesain lebih cepat dan hemat uang. Analisis CFD akan lebih efisien waktu riset dan desain sehingga akan mencapai sasarannya. Proses Simulasi CFD Proses simulasi CFD terdapat tiga tahapan, yaitu : 1. Pra-pemprosesan (Preprocessing) Pra-pemprosesan adalah membuat model dalam paket CAD (Computer aided design), membuat grid yang sesuai dan membuat kondisi batas dan sifatsifat fluidanya. 2. Pencarian solusi (solving) Tahap solving adalah tahap dilakukan perhitungan mengenai kondisikondisi yang telah diterapkan pada tahap pra-pemprosesan guna mendapatkan solusi. Solusi teknik numerik untuk mencari solusi di dalam CFD terdiri atas beberapa metode, yaitu difference, finite element dan spectral method (Tuakia, 2008). 3. Pasca pemprosesan (Post processing) Tahap pasca pemrosesan adalah tahap penyajian hasil simulasi CFD dengan visualisasi warna untuk memudahkan dalam menganalisis. Tampilan hasil yang didapat meliputi : 1. Hasil geometri dan grid yang terbentuk. 2. Plot berdasarkan vektor. 3. Plot berdasarkan kontur. 4. Plot berdasarkan permukaan (2D atau 3D). Prosedur yang terdapat pada semua pendekatan program CFD, yaitu : 1. Pembuatan geometri dari model/problem. 2. Bidang atau volume yang diisi oleh fluida dibagi menjadi sel-sel kecil (meshing) 3. Pendefinisian model fisik 4. Pendefinisian kondisi-kondisi batas 5. Persamaan-persamaan matematika yang membangun CFD diselesaikan secara iteratif, bisa dalam kondisi tunak (steady state) atau transien 6. Analisis dan visualisasi dari solusi CFD (Tuakia, 2008)
METODE Waktu dan Tempat Kegiatan penelitian dilaksanakan mulai bulan Februari 2013 sampai Oktober 2013 dan bertempat di Laboratorium Lapang Siswadhi Soepardjo
7 Leuwikopo dan Laboratorium Energi, Departemen Teknik Mesin dan Biosistem, Fakultas Teknologi Pertanian, IPB, Bogor.
Bahan Bahan yang diperlukan adalah minyak tanah, kayu sebagai pemanas tambahan yang berasal dari biomassa. Data-data fisik dari temulawak dan polikarbonat sebagai data input simulasi.
Alat 1. Alat pengering yang digunakan adalah alat pengering Efek Rumah Kaca (ERK) Hybrid yang dilengkapi dengan tungku sebagai pemanas tambahan hasil rancangan Wulandani dan Nelwan (2009). Alat pengering ini terdiri dari beberapa bagian utama yaitu bangunan rumah kaca dengan dimensi : p = 3.1 m, l = 1.065 m dan t = 2.315 m. Dinding dan atap terbuat dari bahan tembus cahaya (polycarbonate) dengan tebal 0.001 m. Plat penutup pada bagian bawah terbuat dari plat besi (esser) dicat hitam berfungsi sebagai absorber . Produk yang akan dikeringkan ditempatkan pada rak bersusun vertikal dengan dimensi p = 0. 5 m , l = 0.6 m, berjumlah 144 buah terbagi kedalam delapan kolom susunan rak. Masing-masing kolom terdiri dari 18 level tersebar ke dalam dua segmen kiri dan kanan. Tungku biomassa ditempatkan ditengah bangunan yang bertujuan untuk menghasilkan energi panas yang dapat menjangkau kedua segmen ruang pengering. Udara panas yang berasal dari tungku ditransfer melalui heat exchanger (HE) yang terletak diatas tungku. Udara panas tersebut disirkulasikan menggunakan kipas. Gambar alat pengering dapat dilihat pada Gambar 2.
Atap
Dinding transparan Heat exchanger Rak
Lantai absorber Tungku
Gambar 2 Alat pengering Efek Rumah Kaca (ERK) Hybrid 2. – Perangkat lunak a. Proses pembuatan gambar dan meshing : Gambit 2.4.6 b. Proses simulasi dan analisis CFD : ANSYS 12.1 : Fluent 6.3 c. Penunjang analisis data hasil simulasi : Ms. Office & Ms. Excel 2007 - Perangkat keras
8 a. Komputer Dell b. Laptop Axioo 3. Alat pengukuran Termokopel tipe CC yang dapat mengukur suhu antara -200oC – 300oC, hybrid recorder Yokogawa tipe 3058, anemometer model Kanomax tipe 6011, pyranometer model EKO tipe MS-401, multimeter, penggaris, meteran, termometer raksa.
Prosedur Penelitian Pada penelitian ini dilakukan dengan beberapa tahapan yang dapat dilihat pada diagram alir penelitian pada Gambar 3. Persiapan Alat Persiapan alat dilakukan untuk membersihkan alat pengering, memperbaiki beberapa kerusakan yang ada, mengecek kesiapan instrument pengukuran, menentukan titik-titik pengukuran pada alat pengering, melakukan pengukuran dimensi alat serta mengindetifikasi bahan konstruksi alat pengering. Pengukuran dan Perhitungan Pengukuran dilakukan sebanyak 3 kali. Pengukuran pertama dilakukan untuk mendapatkan data-data yang diperlukan pada proses simulasi CFD dan validasi kondisi tanpa menggunakan beban. Pengukuran dilakukan pada pukul 08.00, 12.00 dan 16.00 yang meliputi pengukuran kecepatan angin pada kipas inlet dan ruangan pengering, iradiasi matahari, suhu rak, dinding, atap, lantai serta pipa penukar panas. Pengukuran kedua dan ketiga merupakan pengukuran untuk mengetahui profil sebaran suhu rak pada ruang pengering. Perhitungan digunakan untuk mendapatkan nilai-nilai yang dimasukkan pada simulasi. Perhitungan-perhitungan yang diperlukan diperoleh dari persamaan dibawah ini : 1. Perhitungan aliran laminar ρvD Re = μ
2. Perhitungan parameter kondisi batas pada dinding 2 a. Koefisien pindah panas konveksi (h) 𝑁𝑁𝑁𝑁 𝑥𝑥 𝑘𝑘 h= 𝑙𝑙 2 b. Fluks panas 𝑞𝑞 = (I x τα) – [h x (Td – T∞ )] 𝐴𝐴 3. Perhitungan parameter kondisi batas pada rak. Rak berupa plat berlubang dianggap sebagai porous jump. Nilai α dan C2 dinyatakan dalam persamaan yang terdapat pada Wulandani (2005). 3 a. Permeabilitas permukaan (α) α=
D 2p
ε2
150 (1−ε)2
3 b. Koefisien porous jump (C2)
9 C2 =
3.5 (1−ε) Dp ε 3
4. Perhitungan parameter kondisi batas pada heat exchanger (HE). Nilai k dinyatakan dalam persamaan yang terdapat pada Wulandani (2005). 4 a. Koefisien kehilangan (k) 2 x ΔP k= 2 ρ x v kipas
4 b. Koefisien pindah panas konveksi (h) h = (Nu x k)/D 4 c. Fluks panas (q/a) q/a = h x ΔTln
Pembuatan Model Geometri dan Pendefinisian Kondisi Batas Alat pengering ERK yang disimulasikan terlebih dahulu dibuatkan model geometri 3D dengan menggunakan perangkat lunak Gambit 2.4.6 yang meliputi : atap, lantai, dinding, pipa-pipa heat exchanger, inlet dan outlet beserta rak-rak pengering. Setelah model geometri dibuat, dilakukan pendefinisian bidang batas setiap bagian geometri, seperti atap, dinding dan lantai didefinisikan sebagai wall, pipa-pipa heat exchanger sebagai radiator, inlet sebagai velocity inlet, outlet sebagai outflow dan rak-rak pengering sebagai porous jump. Setelah pendefinisian selesai dilakukan proses meshing atau proses pembuatan grid-grid pada model geometri. Percobaan Pengeringan Tanpa Beban (Percobaan 1 dan 2) Percobaan pengeringan tanpa beban 1 dan 2 adalah percobaan pengeringan yang dilakukan untuk mengetahui profil sebaran suhu rak dari pukul 08.00 – pukul 17.00. Percobaan 1 dilakukan pengukuran suhu rak nomor urut ganjil pada kolom rak ID, IIB, IIID dan IV B. Percobaan II pengukuran suhu rak nomor urut genap pada kolom rak IB, IID, IIB dan IVD. Untuk skema dari posisi rak tersebut dapat dilihat pada Gambar 8. Simulasi dan Validasi Tanpa Beban Tahapan simulasi kondisi tanpa beban adalah mensimulasikan hasil pengukuran dari percobaan pengeringan tanpa beban (percobaan pendahuluan) pada pukul 08.00, 12.00 dan 16.00. Dasar dari penentuan jam tersebut adalah ingin diketahuinya kondisi dan hasil simulasi CFD pada posisi dan intensitas matahari dari kondisi yang berbeda. Hasil simulasi tersebut yang berupa parameter suhu lalu divalidasikan dengan hasil pengukuran pada titik-titik yang telah ditentukan. Simulasi dan Validasi Dengan Beban Tahapan simulasi kondisi menggunakan beban adalah mensimulasikan hasil pengukuran dari percobaan pengeringan menggunakan temulawak sebagai beban berdasarkan hasil penelitian Aritesty (2013) pada pukul 12.00 hari pertama pengeringan. Pemilihan pukul 12.00 adalah karena pada kondisi tersebut, bahan yang dikeringkan masih berada pada tahap dimana kadar air cukup tinggi dan intensitas matahari tinggi sehingga kondisi tersebut memungkinkan terjadinya ketidakseragaman suhu yang tinggi pada antar rak. Hasil simulasi berupa parameter suhu yang telah didapatkan lalu divalidasikan, dan apabila telah valid
10 hasil simulasi tersebut digunakan untuk tahapan simulasi modifikasi posisi penempatan kipas. Simulasi Modifikasi Penempatan Posisi Kipas dan Uji Keseragaman Modifikasi penempatan posisi kipas pada penelitian ini menggunakan lima skenario. Pertama, terdapat dua kipas di tiap ruang pengering kanan dan kiri. Kedua, terdapat empat kipas di tiap ruang pengering. Ketiga, terdapat satu kipas di tiap ruang pengering. Keempat, posisi kipas berada di outlet samping tiap ruang pengering. Kelima, posisi kipas berada di outlet atas tiap ruang pengering. Hasil simulasi dari tiap skenario tersebut diuji nilai keseragamannya dan dilihat pola aliran udara panasnya. Keseragaman yang terbesar atau nilai koefisien variasi terkecil dan sebaran yang paling seragam itulah yang dijadikan rekomendasi. Validasi Model Validasi model dilakukan dengan membandingkan suhu rak hasil simulasi dengan suhu rak hasil pengukuran di lapang. Validasi pada kondisi tanpa beban menggunakan data pada hasil pengukuran pertama yang sudah terlampir pada tabel validasi (Lampiran 2 – Lampiran 4 ) dan tabel data input CFD (Lampiran 7 – Lampiran 10). Sedangkan hasil pengukuran lapang untuk validasi dan data input kondisi batas CFD dengan kondisi memakai beban menggunakan hasil dari penelitian Aritesty (2013) pada pukul 12.00 percobaan hari pertama pengeringan. Hasil tersebut telah diringkas dan terlampir pada Lampiran 5 untuk validasi model dan Lampiran 10 untuk data input CFD. Keakuratan hasil simulasi dengan hasil pengukuran dinyatakan dalam persentase error. Error dinyatakan dengan persamaan: Error =
suhu simulasi (℃)−suhu pengukuran (℃) suhu pengukuran (℃)
𝑥𝑥 100%
Menurut Karmeli dalam Nuryawati (2011), distribusi suhu hasil simulasi juga divalidasi dengan menghitung koefisien keseragaman dan koefisien variasi. Apabila nilai koefisien keseragaman (CU) lebih dari 75% dan nilai koefisien variasi (CV) kurang dari 15%, maka dikatakan bahwa hasil simulasi telah baik. Persamaannya dinyatakan sebagai berikut: CV =
σ |μ |
CU = 100 �1 −
∑𝑛𝑛𝑖𝑖=1|𝑋𝑋 𝑖𝑖 − 𝜇𝜇 | ∑𝑛𝑛𝑖𝑖=1 𝑋𝑋 𝑖𝑖
�
dimana CV adalah koefisien variasi (coefficient of variation), CU adalah koefisien keseragaman (coefficient of uniformity) (%), σ adalah standar deviasi, μ adalah rata-rata suhu hasil simulasi (oC), Xi adalah suhu hasil pengukuran (oC), dan n adalah jumlah data.
11 Mulai
Persiapan alat
Pengukuran dan perhitungan
Data yang diperlukan untuk tahapan simulasi
Percobaan pengeringan tanpa beban (percobaan pendahuluan) Pembuatan model geometri bangunan pengering dan pendefinisian kondisi batas Simulasi tanpa beban
Valid (Error suhu < 10%)
Tidak
Ya Data profil sebaran suhu rak di kedua ruang pengering Data riset Aritesty (2013)
Percobaan pengeringan tanpa beban (percobaan 1 dan 2) Percobaan pengeringan dengan beban menggunakan temulawak Simulasi dengan beban
Valid (Error suhu < 10%)
Tidak
Pendefinisian kondisi batas
Ya Simulasi modifikasi daya dan posisi penempatan kipas dan uji keseragaman CV dan standar deviasi terkecil serta pola aliran udara yang paling seragam
Standar keseragaman -Nilai CV -Nilai Standar deviasi
Rekomendasi
Selesai
Gambar 3 Diagram alir penelitian
HASIL DAN PEMBAHASAN Tahapan Awal Simulasi Tahapan Pra-pemrosesan (Preprocessing) 1. Pembuatan model geometri 3 dimensi bangunan pengering ERK dan rak disesuaikan dengan dimensi aslinya.
12
Gambar 4 Tampilan isometri geometri 3D pengering ERK
Gambar 5 Tampilan depan geometri 3D pengering ERK 2. Pembuatan grid (mesh) model, jenis mesh harus dapat memenuhi kriteria baik sehingga konvergensi proses simulasi lebih cepat tercapai. Kriteria yang digunakan adalah nilai worst element ≤ 0.85. Pengering ERK yang dibuatkan ke mesh model menggunakan elemen : Tet/Hybrid, tipe : Tgrid, ukuran interval : 0.04 cm – 0.05 cm. Pemilihan jenis mesh didasarkan pada pemilihan yang dilakukan secara otomatis oleh software Gambit berdasarkan skema dan geometri yang tersusun serta ukuran mesh yang didapatkan berdasarkan trial and error untuk mendapatkan kualitas / nilai worst element ≤ 0.85.
13
Gambar 6 Tampilan hasil mesh geometri pengering ERK Asumsi, Kondisi Batas dan Perhitungan Parameter Kondisi Batas Asumsi 1. Udara tidak termampatkan (incompressible), ρ konstan 2. Bilangan Prandtl udara konstan (panas jenis, konduktivitas dan viskositas udara konstan) 3. Udara bergerak dalam keadaan steady 4. Udara lingkungan dianggap konstan selama simulasi 5. Kecepatan angin dianggap konstan 6. Aliran udara dianggap laminar. Aliran laminar jika Re < 2000 (Holman, 1997). Nilai dan perhitungan aliran laminar pada ruang pengering ditunjukkan pada Lampiran 1. Kondisi Batas Pengering ERK dibatasi oleh dinding yang berbentuk bangunan segiempat dengan atap miring serta outlet yang terletak di kanan kiri bangunan berbentuk trapesium. Contoh parameter penentu kondisi batas pada analisis ini adalah diambil dari kondisi batas untuk data validasi pada pukul 08.00 : 1. Inlet sekaligus kipas dianggap sebagai velocity inlet : 3 m/s 2. Outlet dianggap sebagai outflow dengan ratio bukaan 0.1 3. Dinding pada geometri pengering ERK terdapat dua bagian, yaitu dinding ruang pengering sebelah kanan dan kiri. Parameter kondisi batas untuk dinding sebelah kanan adalah sebagai berikut : - Fluks panas = 259 W/m2 - Suhu dinding = 34.6 oC - Ketebalan dinding = 0.001 m - Koefisien pindah panas konveksi (h) = 2.98 W/m2 oC - Emisivitas = 0.61 (Hartini, 2010) - Laju pembentukan panas = 0 W/m2 4. Rak dianggap sebagai porous jump. Kondisi batas porous jump adalah kondisi batas yang digunakan untuk mewakili kondisi suatu lapisan tipis yang memiliki karakteristik
14 kecepatan (penurunan tekanan), yang biasanya digunakan untuk pemodelan sekat, layar, penyaring dan radiator yang tidak terkait dengan pindah panas (berdasarkan ANSYS FLUENT User Guide). Rak dengan kondisi tanpa beban memliki parameter sebagai berikut : - Permeabilitas permukaan (α) = 0 m2 - Tebal rak = 0 m - Koefisien porous jump (C2) = 0 (1/m) 5. Pipa penukar panas (heat exchanger) yang terdapat di tengah bangunan pengering sebagai radiator terdiri dari dua bagian, yaitu bagian kanan dan kiri bangunan pengering. Parameter kondisi batas untuk HE (heat exchanger) bagian kanan adalah sebagai berikut : - Koefisien kehilangan (k) = 42.02 - Koefisien pindah panas konveksi (h) = 48.62 W/m2 oC - Suhu penukar panas = 31.1 oC - Fluks panas = 213.08 W/m2 Untuk data input kondisi batas CFD secara lengkap pada kondisi tanpa beban dan dengan beban dapat dilihat pada Lampiran 7 – Lampiran 10. Perhitungan Parameter Kondisi Batas 1. Perhitungan parameter kondisi batas pada dinding menggunakan persamaan sebagai berikut : a. Koefisien pindah panas konveksi (h) Koefisien pindah panas konveksi pada dinding dianggap sebagai konveksi bebas. Nilai h ini dinyatakan dalam persamaan yang diberikan Churchill dan Chu dalam Holman (1997) : gβTx 3
Gr = 2 ν Ra = Gr x Pr Nu = �0.825 + h=
𝑁𝑁𝑁𝑁 𝑥𝑥 𝑘𝑘
0.387 𝑅𝑅𝑅𝑅 1/6
[1+(0.492/𝑃𝑃𝑃𝑃 )9/16 ]
2
8/27 �
𝑙𝑙
Iradiasi surya didefinisikan sebagai fluks dari atap dan dinding – dinding bangunan pengering. Fluks panas dapat dihitung melalui persamaan sebagai berikut: 𝑞𝑞 = (I x τα) – h ( Td – T∞ ) 𝐴𝐴
2. Perhitungan parameter kondisi batas pada penukar panas. a. Koefisien kehilangan (k) Pη ΔP = v inlet x A inlet
k=
2 x ΔP
ρ x v 2kipas
b. Koefisien pindah panas konveksi (h)
15 Koefisien pindah panas konveksi pada penukar panas dihitung sebagai aliran menyilang rangkunan tabung (tube bank) yang dinyatakan di persamaan dalam Cengel (2003) : ρf =
p
RT
SD = �𝑆𝑆𝑆𝑆2 + (𝑆𝑆𝑆𝑆/2)2 AD = (SD – D) L AT = (ST – D) L 2 AD > AT v maks = Re =
𝑆𝑆𝑆𝑆
𝑆𝑆𝑆𝑆−𝐷𝐷
ρf v maks D
𝑣𝑣
μ
Nu = 0.35 (ST/SL)0.2 (Re)0.6 (Pr)0.36 (Pr/Prs)0.25 Faktor koreksi = 0.76 NuD = Nu x faktor koreksi h = (NuD x k)/(D) c. Fluks panas (q/a) A=NxπxDxL ṁ = ρ∞ x v x n x SL x L Te = THE – [(THE – Tinlet) x (Exp (–A x h))/(ṁ x Cp)] ΔTln = [(THE – Te) – (THE – Tinlet)]/[Ln (THE – Te)/(THE – Tinlet)] q/a = h x ΔTln Rincian perhitungan parameter kondisi batas pada kondisi pukul 08.00 tanpa beban dapat dilihat pada Lampiran 1. Secara keseluruhan tahapan dalam simulasi dapat dilihat pada diagram alir di Gambar 7. Sedangkan untuk tahapan penggunaan ANSYS Fluent untuk simulasi dapat dilihat pada Lampiran 20.
16 Mulai
Pembuatan geometri dan meshing
Pendefinisian bidang batas pada geometri
Pengecekan mesh
Kualitas mesh baik ? (worst element <0.85)
Tidak
worst element > 0.85 dan mendekati nilai 1
Ya Data sifat fisik Penentuan kondisi batas
Proses numerik
Ya
Iterasi error ? Tidak
Plot distribusi temperatur
Selesai
Gambar 7 Diagram alir simulasi
Penyebaran Suhu Rak Hasil Pengukuran Profil Sebaran Suhu Rak Bernomor Ganjil (Pengukuran 1) Hasil pengukuran suhu rak pengering ERK dikelompokkan berdasarkan ruang pengering, posisi rak dan urutan posisi rak dari posisi paling atas. Pengambilan data dilakukan pada hari yang berbeda. Gambar 8 memperlihatkan skema dari posisi rak.
17
1
IV
17
III II
B I
D
Gambar 8 Posisi dan penomoran rak Keterangan : Angka romawi I – IV, menunjukkan baris ruang pengering Huruf D dan B, menunjukkan posisi rak. D = Depan, B = Belakang Angka 1- 17 menunjukkan posisi rak dimulai dari urutan paling atas. Grafik perubahan suhu rak bernomor ganjil dapat dilihat pada Gambar 9 – Gambar 12. Data suhu rak pada pengukuran 1 dapat dilihat pada Lampiran 12 dan 13.
Gambar 9 Grafik suhu rak I D
18
Gambar 10 Grafik suhu rak II B
Gambar 11 Grafik suhu rak III D
Gambar 12 Grafik suhu rak IV B Gambar 9 – Gambar 12 menunjukkan profil suhu rak yang memiliki pola yang cenderung seragam. Profil memperlihatkan perubahan suhu, dimana terjadi penurunan suhu mulai dari rak yang paling atas sampai yang terbawah. Tabel 1 memperlihatkan tingkat keseragaman suhu tiap rak selama pengukuran yang dinyatakan dengan nilai standar deviasi. Rak-rak yang berada di ruang pengering I,II cenderung lebih seragam dibandingkan rak-rak yang berada
19 di ruang pengering III,IV. Hal itu terlihat dari nilai keseragaman yang lebih kecil dan terjadi akibat beberapa sebaran data suhu seperti atap, dinding dan lantai pada ruang pengering I,II selama pengukuran lebih kecil dari ruang pengering III,IV. Keseragaman suhu rak pada waktu yang sama selama pengukuran dapat dilihat pada Lampiran 20. Data suhu lantai, atap, dinding dan HE selama pengukuran dapat dilihat pada Lampiran 18. Tabel 1 Tingkat keseragaman suhu rak bernomor ganjil ID Standar deviasi Rak
IIB Standar deviasi Rak
IIID Standar deviasi Rak
IVB Standar deviasi Rak
1 3 5 7 9 11 13 15 17
1 3 5 7 9 11 13 15 17
1 3 5 7 9 11 13 15 17
1 3 5 7 9 11 13 15 17
(oC) 8.81 6.77 7.17 7.06 4.76 5.13 5.00 4.34 4.11
Ratarata
5.91
Ratarata
(oC)
(oC) 8.40 4.66 3.46 3.27 3.48 2.94 2.68 2.69 2.31 3.77
Ratarata
(oC)
(oC) 3.90 4.35 3.84 4.45 3.88 4.64 3.11 3.50 2.74 3.91
Ratarata
(oC)
(oC) 8.82 5.73 4.76 4.58 4.27 3.92 3.22 3.43 2.86 4.62
(oC)
Dari keseluruhan data yang terukur, rak nomor 17 memiliki nilai suhu yang paling kecil serta standar deviasi yang paling kecil. Hal ini dapat terjadi, karena rak 17 terdapat pada posisi paling bawah, tertutupi oleh rak-rak yang ada diatasnya, sehingga penerimaan sinar matahari menjadi kecil dan tidak tidak terlalu terpengaruh. Keseragaman paling kecil terjadi dari pukul 12.00 – 13.30, dimana puncak intensitas matahari terjadi (Gambar 13) . Profil suhu yang terjadi bersifat fluktuatif. Hal ini sangat dipengaruhi terhadap besar kecilnya penyinaran matahari. Adapun penyinaran matahari tersebut sangat dipengaruhi oleh faktor penutupan awan. Untuk mengetahui distribusi suhu rak selama pengukuran 1 dapat dilihat pada Lampiran 14. Iradiasi Surya (Pengukuran 1) Radiasi surya yang terukur selama proses percobaan terlihat seperti pada Gambar 13.
20
Gambar 13 Profil iradiasi surya pengukuran 1 Selama pengambilan data iradiasi surya, besarnya iradiasi yang terukur bersifat fluktuatif. Iradiasi terbesar ada pada pukul 13.30 sebesar 914.29 W/m2 sedangkan yang terendah pada pukul 17.00 sebesar 28.57 W/m2. Iradiasi rata-rata yang diperoleh sebesar 490.98 (W/m2). Faktor utama terjadinya iradiasi surya yang bersifat fluktuatif adalah faktor cuaca terutama penutupan awan. Iradiasi pada pukul 13.00 mengalami penurunan yang drastis, tetapi mengalami kenaikan kembali pada pukul 13.30, hal ini terjadi karena adanya penutupan oleh awan gelap yang sesaat. Profil Sebaran Suhu Rak Bernomor Genap (Pengukuran 2) Gambar 14 – Gambar 17 menunjukkan grafik suhu rak bernomor genap.
Gambar 14 Grafik suhu rak I B
21
Gambar 15 Grafik suhu rak II D
Gambar 16 Grafik suhu rak III B
Gambar 17 Grafik suhu rak IV D Gambar 14 – Gambar 17 menunjukkan profil suhu rak yang memiliki pola yang cenderung seragam. Profil memperlihatkan perubahan suhu, dimana terjadi penurunan suhu mulai dari rak yang paling atas sampai yang terbawah dan cenderung memiliki bentuk pola sebaran yang sama. Keseragaman suhu rak dapat dilihat pada Tabel 2. Keseragaman suhu rak pada waktu yang sama selama pengukuran dapat dilihat pada Lampiran 21.
22 Tabel 2 Tingkat keseragaman suhu rak bernomor genap IB Standar deviasi Rak
2 4 6 8 10 12 14 16
(oC) 6.09 2.55 4.13 3.80 3.58 3.80 3.66 3.27
Ratarata
3.86
o
( C)
IID Standar Rak deviasi
2 4 6 8 10 12 14 16
(oC) 6.34 5.74 4.15 3.23 3.32 3.10 2.62 2.39
Ratarata
3.86
o
( C)
IIIB Standar Rak deviasi
2 4 6 8 10 12 14 16
(oC) 7.31 6.54 6.35 6.01 5.30 5.45 5.15 3.15
Ratarata
5.66
o
( C)
IVD Standar Rak deviasi
2 4 6 8 10 12 14 16
(oC) 2.74 6.33 6.48 6.32 5.28 5.44 4.71 3.10
Ratarata
5.05
o
( C)
Tingkat keseragaman pada rak genap memperlihatkan, hasil yang sama pada pengukuran rak ganjil atau pengukuran 1. Rak pada bagian I,II memiliki tingkat keseragaman yang lebih kecil dibandingkan bagian III,IV. Hal ini juga disebabkan karena beberapa data distribusi suhu seperti atap, dinding dan lantai selama pengukuran lebih kecil dibandingkan ruang pengering III,IV. Keseragaman suhu rak terkecil terjadi pada pukul 12.30 – 13.00, dimana puncak intensitas matahari terjadi sama seperti pada pengukuran 1 (Gambar 18). Data pengukuran suhu atap, lantai, HE dan dinding pada ruang pengering selama pengukuran dapat dilihat pada Lampiran 19. Iradiasi Surya (Pengukuran 2) Selama pengambilan data iradiasi surya, besarnya iradiasi yang terukur bersifat fluktuatif. Iradiasi terbesar ada pada pukul 13.00 sebesar 985.71 W/m2 sedangkan yang terendah pada pukul 17.00 sebesar 85.71 W/m2. Iradiasi rata-rata yang diperoleh sebesar 532.33 W/m2. Faktor utama terjadinya iradiasi surya yang bersifat fluktuatif adalah faktor cuaca terutama penutupan awan. Iradiasi pada pukul 13.30 mengalami penurunan yang drastis, tetapi mengalami kenaikan kembali pada pukul 14.00, hal ini terjadi karena adanya penutupan oleh awan gelap yang hanya sesaat. Data pengukuran iradiasi surya pada pengamatan pertama dan kedua dapat dilihat pada Lampiran 6. Grafik iradiasi surya pengamatan kedua dapat dilihat pada Gambar 18.
23
Gambar 18 Profil iradiasi surya pengukuran 2 Berdasarkan hasil pengukuran 1 dan 2 didapatkan bahwa tingkat keseragaman suhu rak pada ruang pengering I,II lebih baik dibandingkan ruang pengering III,IV. Hal itu terlihat dari nilai standar deviasi pada ruang pengering I,II lebih kecil dibandingkan ruang pengering III,IV.
Kontur Suhu Ruang Pengering ERK Hasil Simulasi CFD Pada Kondisi Tanpa Beban Simulasi CFD pengering ERK bertujuan untuk mengetahui pola dan sebaran suhu yang terjadi pada ruang pengering. Simulasi dengan kondisi tanpa beban dilakukan dengan tiga waktu yang berbeda yaitu pada pukul 08.00, 12.00 dan 16.00. Data input kondisi batas simulasi tanpa beban dapat dilihat pada Lampiran 7, 8 dan 9. Suhu Ruang Pengering Pada Pukul 08.00 Kontur suhu pada potongan sumbu z 0.5325 m (Gambar 19) atau pada pertengahan bangunan pengering memperlihatkan perubahan suhu yang terjadi pada atas ruang pengering menuju ke bawah dan perbedaan kontur suhu pada ruang pengering I,II dan III,IV. Ruang pengering I,II memiliki kontur suhu yang lebih beragam dibandingkan ruang pengering III,IV. Pada ruang pengering I,II suhu mengalami penurunan mulai dari atas menuju pertengahan bangunan, lalu mengalami kenaikan hingga ke bagian lantai. Hal ini terjadi karena pada bagian atas suhu terpengaruh oleh radiasi surya sedangkan pada bagian bawah dipengaruhi oleh suhu lantai. Kontur suhu yang berbeda antara kanan dan kiri disebabkan karena matahari pada pukul 08.00 berada di arah ruang I,II ruang pengering, sehingga penerimaan sinar matahari lebih dominan dibandingkan ruang pengering III,IV.
24
Gambar 19 Kontur suhu simulasi tanpa beban pukul 08.00 (Z = 0.5325 m)
Gambar 20 Kontur suhu simulasi tanpa beban pukul 08.00 pada sumbu X Distribusi suhu rak di kedua ruang pengering pada simulasi pukul 08.00 diperlihatkan pada Tabel 3. Tabel 3 Distribusi suhu rak pada simulasi tanpa beban pukul 08.00 Rata-rata (oC) Ruang pengering III,IV Ruang pengering I,II
Pukul 08.00 Standar deviasi (oC)
Maksimum (oC)
Minimum CV (%) (oC)
31.26
0.54
33.81
30.72
1.73
35.08
0.96
38.3
33.17
2.73
Suhu Ruang Pengering Pada Pukul 12.00 Simulasi pada pukul 12.00 memperlihatkan pola sebaran suhu (Gambar 21) yang cenderung sama dengan simulasi pukul 08.00. Gambar 21 memperlihatkan pola sebaran pada kedua ruang pengering cenderung sama, walaupun pada ruang pengering I,II sedikit lebih beragam. Hal ini terlihat dari distribusi suhu rak yang terjadi pada kedua ruang pengering. Distribusi suhu rak pada kedua pengering dapat dilihat pada Tabel 4.
25
Gambar 21 Kontur suhu simulasi tanpa beban pukul 12.00 (Z = 0.5325 m)
Gambar 22 Kontur suhu simulasi tanpa beban pukul 12.00 pada sumbu X Tabel 4 Distribusi suhu rak pada simulasi tanpa beban pukul 12.00 Rata-rata (oC) Ruang pengering III,IV Ruang pengering I,II
Pukul 12.00 Standar deviasi (oC)
Maksimum (oC)
Minimum CV (%) (oC)
42.66
1.4
49.7
41.56
3.28
43.46
1.7
51.94
42.26
3.9
Suhu Ruang Pengering Pada Pukul 16.00 Hasil simulasi pada pukul 16.00 menunjukkan hasil yang cukup berbeda dibandingkan hasil simulasi pada pukul 08.00 dan pukul 12.00. Hal ini bisa dilihat dari pola sebaran suhu pada Gambar 23 dan Gambar 24. Sebaran suhu pada ruang pengering III,IV lebih beragam dibanding ruang pengering I,II yang sangat beragam. Ruang pengering I,II memiliki kontur suhu yang sama, terlihat dari warna yang dihasilkan hanya berupa satu warna dan berada pada kisaran suhu 39 o C. Hal ini berbeda jauh dari hasil pengukuran dimana sebaran suhu yang
26 dihasilkan cukup beragam. Sedangkan pada ruang pengering III,IV kontur suhu yang dihasilkan cukup sama seperti pada hasil simulasi pukul 08.00 dan 12.00. Hal ini kemungkinan dapat disebabkan karena adanya beberapa pengaturan kondisi batas dan pengaturan lainnya pada tahapan simulasi menggunakan Ansys Fluent yang kurang sesuai dengan kondisi saat pengukuran. Untuk mengetahui distribusi suhu rak pada kedua ruang pengering yang dihasilkan dari simulasi, dapat dilihat pada Tabel 5. Tabel 5 Distribusi suhu rak pada simulasi tanpa beban pukul 16.00 Rata-rata (oC) Ruang pengering III,IV Ruang pengering I,II
Pukul 16.00 Standar deviasi (oC)
Maksimum (oC)
Minimum CV (%) (oC)
34.36
1.14
38.46
33.16
3.32
39.46
0.089
39.55
39.13
0.23
Gambar 23 Kontur suhu simulasi tanpa beban pukul 16.00 (Z = 0.5325 m)
Gambar 24 Kontur suhu simulasi tanpa beban pukul 16.00 pada sumbu X
27 Kontur Suhu Ruang Pengering ERK Hasil Simulasi CFD Pada Kondisi Dengan Beban Simulasi CFD pada kondisi menggunakan beban yaitu mensimulasikan pola sebaran suhu di ruang pengering ERK pada saat pukul 12.00 hari pertama pengeringan. Beban atau komiditi yang digunakan adalah temulawak. Data input kondisi batas dan data hasil pengukuran untuk keperluan simulasi CFD dan validasi diperoleh dari hasil riset yang dilakukan Aritesty (2013). Kontur suhu pada Gambar 25 dan Gambar 26 memperlihatkan suhu pada bagian atas ruang pengering memiliki nilai yang paling tinggi, lalu menurun menuju bawah ruang pengering. Pola sebaran suhu pada simulasi ini memiliki kecerendungan yang sama terhadap pola sebaran pada hasil simulasi tanpa beban pukul 08.00 dan pukul 12.00. Hasil simulasi ini juga memiliki kecerendungan yang sama dengan kondisi yang sebenarnya, dimana suhu pada rak atas memiliki suhu yang paling tinggi, lalu menurun ke rak bawah.
Gambar 25 Kontur suhu simulasi dengan beban pukul 12.00 (Z = 0.5325 m)
Gambar 26 Kontur suhu simulasi dengan beban pada sumbu X Simulasi CFD menggunakan beban juga menghasilkan distribusi suhu rak di kedua ruang pengering. Distribusi suhu rak di kedua ruang pengering dapat dilihat pada Tabel 6. Ruang pengering I,II lebih seragam dibandingkan ruang pengering III,IV.
28 Tabel 6 Distribusi suhu rak pada simulasi dengan beban Rata-rata (oC) Ruang pengering III,IV Ruang pengering I,II
Pukul 12.00 Standar deviasi (oC)
Maksimum (oC)
Minimum CV (%) (oC)
41.29
1.48
47.57
38.66
3.57
40.15
0.94
44.6
39.16
2.35
Validasi Model Validasi diperlukan untuk mengetahui tingkat error yang terjadi pada hasil simulasi terhadap hasil pengukuran. Simulasi pada pukul 08.00 dengan kondisi tanpa beban memiliki tingkat rata-rata kesalahan secara keseluruhan sebesar 8.33% dengan standar deviasi error sebesar 9.52%. Adapun tingkat error terbesar ada pada rak 1 D II sebesar 31.50 % sedangkan yang terkecil pada rak 17 B I sebesar 0.45 %.
Gambar 27 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi (pukul 08.00) Tingkat keseragaman pada hasil simulasi menunjukkan nilai CV sebesar 5.77%. Tingkat error rata-rata secara keseluruhan pada simulasi pukul 12.00 menunjukkan nilai sebesar 7.09 % dengan standar deviasi sebesar 7.76. Nilai error terbesar ada pada rak 1 B I sebesar 22.07 % dan error terkecil ada pada rak 9 B III sebesar 0.77 %.
Gambar 28 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi (pukul 12.00)
29 Tingkat keseragaman pada hasil simulasi menunjukkan nilai CV sebesar 4.78%. Untuk simulasi pada pukul 16.00 dengan kondisi tanpa beban didapatkan nilai error rata-rata keseluruhan sebesar 12.54% dengan standar deviasi sebesar 10.03, nilai terbesar ada pada rak 1 D III sebesar 30.60 % dan yang terkecil ada pada rak 10 D I.
Gambar 29 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi (pukul 16.00) Tingkat keseragaman pada hasil simulasi menunjukkan nilai CV sebesar 7.03%. Simulasi dengan kondisi menggunakan temulawak sebagai beban, didapatkan nilai rata-rata error sebesar 5.34%, standar deviasi sebesar 2.62. Nilai error terbesar ada pada rak 1 D III sebesar 11.75% dan error terkecil ada di rak 8 B III.
Gambar 30 Perbandingan suhu hasil pengukuran dengan suhu hasil simulasi menggunakan beban Tingkat keseragaman pada hasil simulasi menunjukkan nilai CV sebesar 4.08%. Berdasarkan nilai keseragaman dan tingkat error pada seluruh validasi model, didapatkan bahwa simulasi pada pukul 12.00 dengan tanpa beban dan pada pukul 12.00 pada kondisi menggunakan beban merupakan model simulasi yang baik. Dikarenakan nilai CV dan CU sudah sesuai dengan standar validasi yang ada. Sehingga untuk model modifikasi, dapat dilakukan dengan menggunakan model simulasi pada pukul 12 dengan kondisi menggunakan beban. Model validasi pada
30 pukul 08.00 dan 16.00 dapat dikatakan cukup baik, karena nilai error yang terjadi masih berkisar di angka 10% dan nilai CV dibawah 15%. Error yang terjadi antara hasil simulasi dengan hasil pengukuran pada simulasi pada pukul 08.00 dan 16.00 disebabkan pada kondisi batas dinding, fluks panas yang dihitung berdasarkan iradiasi matahari. Sedangkan pada pukul 08.00 dan 16.00 iradiasi matahari tidak terlalu besar pengaruhnya terhadap pemanasan di dalam pengering serta posisi matahari yang tidak sejajar. Error terjadi juga disebabkan karena data masukkan atau boundary condition yang diperlukan untuk simulasi ke dalam Fluent berupa data yang cukup sederhana, sedangkan yang terjadi di alam atau kondisi sesungguhnya adalah kondisi yang kompleks dan rumit. Misalkan, untuk data boundary condition berupa dinding atau wall di dalam Fluent hanya dimasukkan satu buah data suhu dan bersifat konstan. Sedangkan kondisi suhu sesungguhnya pada dinding ruang pengering sangatlah kompleks dan berubah-ubah. Menurut Widodo (2009), error yang terjadi pada sebuah simulasi CFD dapat diakibatkan karena kondisi lingkungan yang berubah seperti iradiasi matahari, kecepatan udara masuk, suhu lingkungan dan kelembaban udara lingkungan. Sementara itu, nilai masukkan pada CFD sudah diatur pada kondisi suhu, kecepatan udara masuk dalam ruang pengering pada kondisi tertentu. Kondisi tersebut tidak berubah selama pengukuran CFD. Sehingga terdapat selisih pengukuran antara pengukuran lapang dan pengukuran CFD. Untuk data rekapitulasi error hasil simulasi dapat dilihat pada Lampiran 2 – Lampiran 5.
Modifikasi Posisi Kipas Menggunakan Simulasi CFD Berdasarkan hasil sebaran suhu yang diperoleh dari pengukuran dan simulasi CFD serta validasi simulasi CFD yang dihasilkan, maka dilakukan beberapa modifikasi penempatan posisi kipas yang bertujuan agar keseragaman aliran udara panas dapat tercapai. Modifikasi 1. Terdapat Dua Kipas di Tiap Ruang Pengering Tiap ruang pengering memiliki dua kipas yang dipasang pada bagian dinding atas. Kipas pada ruang pengering III dan IV ditempatkan pada posisi X = 65 cm, Y = 187 cm, Z = 0 cm dan Z = 106.5 cm. Kipas pada ruang pengering bagian I dan II ditempatkan pada posisi X = 245 cm, Y = 187 cm, Z = 0 cm dan Z = 106.5 cm. Kipas yang digunakan memiliki diameter 8 cm dan berdaya 30 watt serta kecepatan kipas yang disimulasikan adalah 0.5 m/s dan 1.5 m/s berdasarkan kecepatan minimum dan maksimum yang dihasilkan kipas tersebut. Gambar 31 dan Gambar 32 menunjukkan kontur suhu pada kecepatan 0.5 dan 1.5 m/s di potongan sumbu Z pada titik Z = 0.5325.
31
Posisi kipas modifikasi 1
Posisi kipas modifikasi 1
Gambar 31 Kontur suhu pada modifikasi 1 v = 0.5 m/s
Gambar 32 Kontur suhu pada modifikasi 1 v = 1.5 m/s Gambar 31 dan Gambar 32 menunjukkan pola sebaran suhu panas yang dihasilkan tidak jauh berbeda. Sebaran udara panas pada kedua ruang pengering cenderung sama, walaupun pada ruang pengering III,IV sebarannya lebih rapat dibandingkan pada ruang pengering I,II. Sebaran suhu ruang pengering pada modifikasi 1 berada di kisaran antara 36oC-55.60oC dan hanya menjangkau daerah pertengahan ruang pengering. Untuk mengetahui distribusi suhu rak yang dihasilkan dapat dilihat pada Tabel 7. Tabel 7 Distribusi suhu rak hasil simulasi CFD pada modifikasi 1
v (m/s) 0.5 1.5
Ruang pengering III,IV Rata- Standar CV rata deviasi (%) o o ( C) ( C) 41.65 1.23 2.95 41.80 1.29 3.09
Ruang pengering I,II Rata- Standar CV rata deviasi (%) o o ( C) ( C) 40.89 1.15 2.81 40.86 1.03 2.52
Ratarata (oC) 41.27 41.33
Keseluruhan Standar deviasi (oC) 1.24 1.26
CV (%) 3.00 3.05
Tabel 7 menunjukkan bahwa distribusi suhu rak pada kedua ruang pengering tidak jauh berbeda dan cenderung hampir sama. Rata-rata suhu yang terjadi pada ruang pengering I,II sedikit lebih rendah dibandingkan ruang
32 pengering III,IV. Hal ini disebabkan karena data kondisi batas ruang pengering I,II yang berupa suhu atap dan dinding yang disimulasikan pada simulasi kondisi pengeringan dengan beban lebih rendah dibandingkan ruang pengering III,IV. Data input kondisi batas tersebut dapat dilihat pada Lampiran 10. Modifikasi 2. Terdapat Empat Kipas di Tiap Ruang Pengering Tiap ruang pengering memiliki empat kipas dengan dua kipas yang dipasang pada bagian dinding atas depan dan belakang. Kipas pada ruang pengering bagian kiri ditempatkan pada posisi X = 33 cm, Y = 187 cm, Z = 0 cm dan Z = 106.5 cm serta X = 98 cm, Y = 187 cm, Z = 0 cm dan Z = 106.5 cm. Kipas pada ruang pengering bagian kanan ditempatkan pada posisi X = 213 cm, Y = 187 cm, Z = 0 cm dan Z = 106.5 cm serta X = 278 cm, Y = 187 cm, Z = 0 cm dan Z = 106.5 cm. Kipas yang digunakan memiliki diameter 8 cm, berdaya 30 watt dan kecepatan kipas yang disimulasikan adalah 0.5 m/s dan 1.5 m/s berdasarkan kecepatan minimum dan maksimum yang dihasilkan kipas tersebut.
Posisi kipas modifikasi 2
Posisi kipas modifikasi 2
Gambar 33 Kontur suhu pada modifikasi 2 v = 0.5 m/s
Gambar 34 Kontur suhu pada modifikasi 2 v = 1.5 m/s Gambar 33 dan Gambar 34 menunjukkan pola sebaran suhu di kedua ruang pengering. Sebaran suhu panas yang terjadi tidak cukup seragam, hal ini terlihat udara yang lebih panas atau kontur suhu yang berwarna hijau tidak dapat
33 menjangkau sampai ke seluruh ruang pengering dan sebaran yang terjadi lebih rendah dibandingkan sebaran suhu panas pada modifikasi 1. Hal ini kemungkinan disebabkan pengaruh jumlah kipas. Penggunaan empat kipas di tiap ruang pengering, mengakibatkan udara panas yang berada di atas ruang pengering lebih cenderung berakibat mengurangi suhu dibandingkan menyebarkannya. Untuk mengetahui distribusi suhu rak yang dihasilkan dapat dilihat pada Tabel 8. Tabel 8 Distribusi suhu rak hasil simulasi CFD pada modifikasi 2
v (m/s) 0.5 1.5
Ruang pengering III,IV Rata- Standar CV rata deviasi (%) o o ( C) ( C) 41.49 0.91 2.19 41.64 0.93 2.23
Ruang pengering I,II Rata- Standar CV rata deviasi (%) o o ( C) ( C) 40.73 0.96 2.36 40.64 0.84 2.07
Ratarata (oC) 41.12 41.14
Keseluruhan Standar CV deviasi (%) o ( C) 0.99 2.41 1.00 2.43
Tabel 8 menunjukkan distribusi suhu rak di kedua ruang pengering tidak jauh berbeda dan cenderung hampir sama serta distribusi suhu rak pada ruang pengering I,II lebih rendah dibandingkan ruang pengering III,IV. Dibandingkan dengan distribusi suhu modifikasi 1, distribusi suhu modifikasi 2 pada kedua ruang pengering maupun secara keseluruhan cenderung lebih rendah. Terlihat dari nilai rata-rata, standar deviasi dan CV yang lebih rendah. Modifikasi 3. Terdapat Satu Kipas di Tiap Ruang Pengering Posisi kipas pada modifikasi ini terletak diatas tiap ruang pengering. Letak kipas ini terdapat pada X = 73 cm, Y = 209 cm, Z = 53 cm dan X = 237 cm, Y= 209 cm, Z = 53 cm. Kipas pada modifikasi ini di dalam simulasi CFD didefinisikan sebagai fan dengan kondisi batas berupa kecepatan minimum dan maksimum serta penurunan tekanan (pressure jump). Kipas yang disimulasikan memiliki diameter 8 cm, dengan kecepatan minimum 0.5 m/s dan kecepatan maksimum 1.5 m/s. Adapun nilai pressure jump yang disimulasikan adalah sebesar 50, 100, dan 4.905. Nilai 50 dan 100 digunakan berdasarkan trial dan error, sedangkan nilai 4.905 didapatkan dari nilai pengukuran langsung di ruang pengering tersebut. Hasil simulasi menunjukkan, pada kondisi batas pressure jump 50 dan 100 kontur suhu yang terjadi sangat seragam bahkan cenderung tidak normal. Hal ini disebabkan karena kondisi pada pressure jump 50 dan 100 terlalu besar dari kondisi kenyataannya. Simulasi pada kondisi pressure jump 4.905 menunjukkan hasil yang lebih baik. Kontur suhu modifikasi 3 pada potongan sumbu Z = 0.5325 m dapat dilihat pada Gambar 35 – Gambar 37.
34
Posisi kipas modifikasi 3
Posisi kipas modifikasi 3
Gambar 35 Kontur suhu pada modifikasi 3 pressure jump 50
Gambar 36 Kontur suhu pada modifikasi 3 pressure jump 100
Gambar 37 Kontur suhu pada modifikasi 3 pressure jump 4.905
35 Tabel 9 Distribusi suhu rak pada modifikasi 3 ΔP 4.905, 50 dan 100 ΔP
4.905 50 100
Ruang pengering III,IV RataStandar CV rata deviasi (%) o o ( C) ( C) 42.33 0.56 1.33 42.65 0.08 0.18 42.64 0.06 0.14
Ruang pengering I,II RataStandar CV rata deviasi (%) o o ( C) ( C) 41.61 0.31 0.75 41.1 0.05 0.13 41.07 0.04 0.11
Ratarata (oC) 41.97 41.88 41.86
Keseluruhan Standar deviasi (oC) 0.58 0.78 0.79
CV (%) 1.38 1.86 1.89
CV
Tabel 9 menunjukkan distribusi suhu rak yang dihasilkan dari tiga kondisi penurunan tekanan (ΔP) yang berbeda. Rata-rata suhu rak yang terjadi hampir sama untuk semua kondisi. Tetapi untuk nilai standar deviasi dan CV pada kondisi ΔP 50 dan 100 memiliki nilai yang jauh lebih rendah dibandingkan ΔP 4.905. Hal ini menunjukkan bahwa sebaran suhu yang terjadi pada ΔP 50 dan 100 sangat seragam. Untuk melihat pengaruh antara nilai pressure jump, efisiensi kipas dan daya kipas dengan hasil sebaran yang terjadi, dilakukan simulasi dengan menggunakan pressure jump 5 – 20. 1,3 1,2 1,1 1 0,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0
CV ruang pengering kiri CV ruang pengering kanan
5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Pressure jump Gambar 38 Grafik hubungan pressure jump dengan CV Tabel 10 Koefisien korelasi antara pressure jump (ΔP) dengan CV ΔP CV ruang pengering kiri ΔP CV ruang pengering kanan
ΔP 1 -0.833 ΔP 1 -0.870
CV Kiri 1 CV Kanan 1
Gambar 38 dan Tabel 10 memperlihatkan pengaruh antara ΔP dengan keseragaman. Grafik menunjukkan bahwa semakin besar nilai ΔP yang disimulasikan, maka semakin besar keseragaman yang terjadi. Hal itu ditunjukkan dengan menurunnya nilai CV seiring bertambahnya ΔP. Hal tersebut juga
36 dibuktikan dengan nilai koefisien korelasi yang bernilai negatif, yang berarti hubungan antara ΔP dengan CV berbanding terbalik. Hubungan antara ΔP dengan distribusi suhu rak yang dihasilkan pada nilai ΔP 4.905, 5 sampai 20, 50 dan ΔP 100 dapat dilihat pada Lampiran 11. Modifikasi 4. Kipas di Outlet Tiap Ruang Pengering Posisi penempatan kipas pada modifikasi 4 berada di kedua outlet atas dan bawah di sisi kiri dan kanan pengering. Kipas yang disimulasikan berdiameter 22 cm berdaya 80 watt, dengan kondisi batas berupa kecepatan sebesar 3 m/s. Kontur suhu pada potongan Z = 0.5325 m dapat dilihat pada Gambar 39.
Posisi kipas modifikasi 4
Posisi kipas modifikasi 4
Gambar 39 Kontur suhu pada modifikasi 4 Gambar 39 memperlihatkan sebaran suhu yang terjadi di kedua ruang pengering. Udara panas yang berada di atas ruang pengering mengalir kebawah ke seluruh ruang pengering. Akan tetapi akibat adanya aliran udara di bawah yang disebabkan oleh kipas di outlet bawah, menyebabkan udara panas tersebut mengalami hambatan, sehingga tidak dapat menjangkau daerah bawah ruang pengering. Distribusi suhu rak yang terjadi pada modifikasi 4 dapat dilihat pada Tabel 11. Tabel 11 Distribusi suhu rak hasil simulasi CFD pada modifikasi 4 Ruang pengering kiri Rata- Standar rata deviasi CV (oC) (oC) 42.34 0.86 2.03
Ruang pengering kanan RataStandar rata deviasi CV (oC) (oC) 42.32 0.79 1.87
Keseluruhan Rata- Standar rata deviasi CV (oC) (oC) 42.33 0.84 1.98
Modifikasi 5. Kipas di Outlet Atas Tiap Ruang Pengering Modifikasi ini posisi kipas berada di tiap outlet atas di sisi kanan dan kiri pengering. Kipas yang digunakan berdiameter 22 cm berdaya 80 watt dengan kecepatan yang disimulasikan sebesar 3 m/s. Gambar 40 memperlihatkan sebaran suhu yang terjadi bila kipas ditempatkan pada outlet atas sisi kanan dan kiri pengering. Sebaran suhu sudah cukup seragam antara ruang pengering kanan dan
37 kiri, walaupun ruang pengering kanan cenderung lebih merata ke seluruh ruangan. Distribusi suhu rak pada modifikasi 5 dapat dilihat pada Tabel 12.
Posisi kipas modifikasi 5
Posisi kipas modifikasi 5
Gambar 40 Kontur suhu pada modifikasi 5 Tabel 12 Distribusi suhu rak hasil simulasi CFD pada modifikasi 5 Ruang pengering kiri RataStandar CV rata deviasi (%) o o ( C) ( C) 43.92 0.43 0.98
Ruang pengering kanan RataStandar CV rata deviasi (%) o o ( C) ( C) 44.28 0.33 0.75
Ratarata (oC) 44.10
Keseluruhan Standar CV deviasi (%) o ( C) 0.42 0.95
SIMPULAN DAN SARAN Simpulan Berdasarkan hasil dari penelitian ini dapat disimpulkan bahwa simulasi pada pukul 12.00 baik dengan kondisi tanpa beban maupun dengan beban memiliki sebaran suhu yang seragam di ruang pengering sebelah kiri dan kanan serta memiliki nilai error yang terkecil sebesar 7.09% dan 5.34%. Sedangkan pada pukul 08.00 dan 16.00 terjadi perbedaan sebaran suhu antara ruang pengering sebelah kanan dan kiri. Hal ini disebabkan karena dipengaruhi arah datang matahari. Penyinaran matahari mempengaruhi tingkat keseragaman suhu pada rak pengering. Hasil validasi menunjukkan bahwa simulasi CFD pada kondisi menggunakan beban pada pukul 12.00 memiliki nilai error terkecil dari keseluruhan simulasi yaitu sebesar 5.34%. Sedangkan pada kondisi tanpa beban pukul 12.00 memiliki error terkecil sebesar 7.09% dan pukul 16.00 memiliki error yang terbesar yaitu sebesar 12.54%. Error pada simulasi CFD dapat terjadi karena data masukkan yang diperlukan untuk keperluan simulasi berupa data yang cukup sederhana dan bersifat konstan, sedangkan kondisi yang terjadi di lapang adalah kompleks dan bersifat berubah-ubah.
38 Modifikasi 5 yaitu posisi kipas berada di tiap outlet atas di sisi kanan dan kiri pengering merupakan posisi penempatan kipas yang ideal dalam mendistribusikan udara panas secara merata yang terlihat dari nilai CV yang dihasilkan sebesar 0.95%, standar deviasi sebesar 0.42oC dan memiliki pola distribusi suhu yang paling seragam. Suhu yang terjadi pada modifikasi 5 berkisar antara 38.9 oC – 55.52oC dengan rata-rata 44.10oC, dimana kisaran suhu tersebut baik digunakan untuk pengeringan temulawak.
Saran 1. Hasil simulasi pada kondisi tanpa beban pukul 16.00 belum menunjukkan kontur suhu yang baik. Perlu dilakukan pengaturan kembali beberapa parameter kondisi yang terdapat pada program Fluent dan simulasi tambahan agar mengurangi tingkat error yang terjadi dan kontur suhu yang lebih mendekati kenyataan di lapang. 2. Dilakukannya simulasi tambahan (berupa penambahan iterasi) dan pengaturan ulang beberapa parameter kondisi batas yang terdapat pada Fluent untuk proses 5 modifikasi posisi kipas, terutama pada modifikasi 3 untuk mendapatkan hasil simulasi yang lebih baik. Penambahan skenario posisi penempatan kipas beserta kriteria kipasnya agar mendapatkan hasil keseragaman suhu yang lebih baik dan optimisasi efisiensi serta penghematan energi. 3. Penggunaan kipas saat proses pengeringan disarankan digunakan pada siang hari, karena pada saat itu keseragaman yang terjadi paling kecil. 4. Perlu dilakukan uji lapang terhadap hasil simulasi modifikasi 5 agar diketahui tingkat error yang terjadi dan keseragaman yang terjadi.
DAFTAR PUSTAKA Aritesty E. 2013. Uji Performansi Alat Pengering Efek Rumah Kaca (ERK) Hybrid Tipe Rak untuk Pengeringan Temulawak (Curcuma xanthorizza Roxb.) [skripsi]. Bogor (ID): Institut Pertanian Bogor. Brooker DB. Bakker-Arkema, C.W.Hall. 1974. Drying Cereal Grain. Westport, Connecticut [USA]. The AVI Publishing Co. Inc., Cengel YA. 2003. Heat Transfer, Published by Mc-Graw Hill Companies, Inc., 1221 Avenue of the Americas, New York, NY 10020. Fellow PJ. 2001. Food Processing Technology, Principles and Practices. CRC Press, Boca Raton, Boston, New York, Washington. Hall, C.W. 1957. Processing Equipment for Agricultural Product. Edward Brothers, Inc. Arbor, Michigan. Hartini S. 2010. Uji Performansi Alat Pengering Efek Rumah Kaca (ERK) Hybrid dan Tungku Biomassa sebagai Sistem Pemanas Tambahan untuk Pengeringan Biji Pala (Myristica sp) [skripsi]. Bogor (ID): Institut Pertanian Bogor.
39 Henderson MS, Perry ME. 1976. Agricultural Process Engineering. Third Edition. The AVI PublishingCompany, Inc., Wesport, Connecticut, USA. Holman JP. 1997. Perpindahan Kalor. Edisi keenam. Alih Bahasa: Jasjfi, E. Penerbit Erlangga. Jakarta. Mark JE. 2007. Physical Properties of Polymers Handbook Second Edition [Internet]. [diunduh 2013 Mei 1]. Tersedia pada : http://bib.convdocs.org/docs/24/23048/conv_1/file1.pdf Nugraha IB. 2005. Simulasi Pola Aliran Udara, RH dan Suhu Ruang Pengering dengan Teknik Computational Fluid Dynamics (CFD) Pada Proses Pengeringan Mahkotadewa [skripsi]. Bogor (ID): Institut Pertanian Bogor. Nuryawati T. 2011. Simulasi Distribusi Suhu dan Pola Aliran Udara dalam Rumah Tanaman Tipe Modified Arch Menggunakan Computational Fluid Dynamics [tesis]. Bogor (ID): Institut Pertanian Bogor. Sitanggang H. 2010. Pengujian dan Simulasi Mesin Pengering Produk Pertanian Tenaga Surya dengan Kapasitas 25 Kg/Siklus Menggunakan Program CFD [skripsi]. Medan (ID): Universitas Sumatera Utara. Tuakia F. 2008. Dasar-Dasar Computational Fluid Dynamics Menggunakan Fluent, Informatika, Bandung. Widodo P. 2009. Kajian Pola Sebaran Aliran Udara Panas pada Model Pengering Efek Rumah Kaca Hibrid Tipe Rak Berputar Menggunakan Computational Fluid Dynamics [tesis]. Bogor (ID): Institut Pertanian Bogor. Wulandani D, Nelwan LO, Kamaruddin A, Indra SA. 2003. Analisis Distribusi Suhu dan Kecepatan Aliran Udara dalam Ruang Pengering Berenergi Surya Menggunakan CFD. Buletin Keteknikan Pertanian. 17(1):68-76. Wulandani D. 2005. Kajian Distribusi Suhu, RH dan Aliran Udara Pengering untuk Optimasi Disain Pengering Efek Rumah Kaca [disertasi]. Bogor (ID): Institut Pertanian Bogor. Wulandani D, Nelwan LO. 2009. Rancang Bangun Kolektor Surya Tipe Plat Datar dan Konsentrator Surya Penghasil Panas pada Pengering ProdukProduk Pertanian. Bogor (ID): Prosiding Seminar Hasil-Hasil Penelitian IPB. Versteeg, H.K. dan W. Malalasekera. 1995. An introduction to computational fluid dynamics. The finite volume method. Longman Sc. & Techical. Malaysia. Yunus YM, Al-Kayiem HH. 2013. Simulation of Hybrid Solar Dryer. Proceeding of International Conference on Energy and Environment [Internet]; [diunduh 2014 Jan 5]; 16 (2013):012143.doi:10.1088/17551315/16/1/012143. Tersedia pada : iopscience.iop.org/17551315/16/1/012143/pdf/1755-1315_16_1_012143.pdf
40 Lampiran 1 Perhitungan data input CFD pada pukul 08.00 tanpa beban 1) Aliran udara ruang pengering sebelah kanan dianggap laminar : Re = 1261.81 Sifat termofisik udara : suhu inlet = 24.9 oC suhu outlet = 33 oC 24.9 + 33 suhu fluida operasi = = 28.95 oC = 301.95 K 2
Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl
Nilai 1.1704 kg/m3 1.0058 kj/kg oC 0.0264 w/ oC 1.855 x10-5 kg/ms 15.89 x10-6 m2/s 0.7076
Jika p dan w masing-masing adalah panjang dan jarak antar rak, dan kecepatan udara yang lewat di antara rak mempunyai kecepatan 0.12 m/detik - Diameter spesifik (D) = 4A/P 4 pw = = Re = =
2(p+w) 4 (0.5 x 0.1)
2 (0.5 +0.1)
= 0.17
ρvD μ
1.1704 x 0.12 x 0.17 1.855x10 −5
= 1261.81
2) Aliran udara ruang pengering sebelah kiri dianggap laminar : Re = 1070.38 Sifat termofisik udara : suhu inlet = 24.9 oC suhu outlet = 27 oC 24.9 + 27 suhu fluida operasi = = 25.95 oC = 298.95 K 2
Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl
Nilai 1.1823 kg/m3 1.00569 kj/kg oC 0.0262 w/ oC 1.841 x 10-5 kg/ms 15.60 x 10-6 m2/s 0.7083
Jika p dan w masing-masing adalah panjang dan jarak antar rak, dan kecepatan udara yang lewat di antara rak mempunyai kecepatan 0.12 m/detik Diameter spesifik (D) = 4A/P 4 pw = 2(p+w)
41 = Re = =
4 (0.5 x 0.1)
2 (0.5 +0.1)
= 0.17
ρvD μ
1.1823 x 0.12 x 0.17 1.841x10 −5
= 1070.38
3) Dinding sebelah kanan Suhu dinding = 34.6 oC Suhu udara lingkungan = 25 oC 34.6 + 25 Suhu fluida operasi = = 29.8 oC = 302.8 K 2
Sifat Nilai Massa jenis 1.1674 kg/m3 Panas jenis 1.0059 kj/kg oC Konduktivitas panas 0.0265 w/ oC Viskositas dinamik 1.859 x 10-5 kg/ms Viskositas kinematik 15.79 x 10-6 m2/s Bilangan Prandtl 0.7074
Koefisien pindah panas konveksi pada dinding dianggap sebagai konveksi bebas. Nilai h ini dinyatakan dalam persamaan yang diberikan Churchill dan Chu dalam Holman (1997) a) Mencari nilai koefisien pindah konveksi (h) Gr = =
gβTx 3 ν2
1
9.8 ( 302 .8) (34.6−25)(1.87)3 (15.97x10 −6 )2
= 7.97x109
Ra = Gr x Pr = 7.97x109 x 0.7074 = 5.64 x109 Nu = [0.825 + = [0.825 +
0.387 𝑅𝑅𝑅𝑅 1/6
]1/2
[1+(0.492/𝑃𝑃𝑃𝑃 )9/16 ]8/27 0.387 (5.64𝑥𝑥10 9 )1/6 1/2
[1+(0.492/0.7074)9/16 ]8/27 16.328
= [0.825 + ]1/2 1.1932 = 210.479
]
42 h=
𝑁𝑁𝑁𝑁 𝑥𝑥 𝑘𝑘
𝑙𝑙 210.479 𝑥𝑥 0.0265
= 1.87 = 2.98 W/ m2 oC b) Fluks panas 𝑞𝑞
𝐴𝐴
= (I x τα) – h ( Td – T∞ ) = (462 x 0.61) – 2.98 (34.6 – 25) = 253.2 W/m2
4) Dinding sebelah kiri Suhu dinding = 30.8 oC Suhu udara lingkungan = 25 oC 30.8 + 25 Suhu fluida operasi = = 27.9 oC = 300.9 K 2
Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl
Nilai 1.1742 kg/m3 1.0058 kj/kg oC 0.0263 w/ oC 1.850 x 10-5 kg/ms 15.78 x 10-6 m2/s 0.7078
Koefisien pindah panas konveksi pada dinding dianggap sebagai konveksi bebas. Nilai h ini dinyatakan dalam persamaan yang diberikan Churchill dan Chu dalam Holman (1997) a) Mencari nilai koefisien pindah konveksi (h) Gr = =
gβTx 3
ν2 1 9.8 ( 300 .9) (30.8 – 25)(1.87)3 (15.78x10 −6 )2
= 4.96x109
Ra = Gr x Pr = 4.96x109 x 0.7078 = 3.51x109 Nu = [0.825 + = [0.825 + h=
𝑁𝑁𝑁𝑁 𝑥𝑥 𝑘𝑘
= 181.488
𝑙𝑙 181.488 𝑥𝑥 0.0263
= 1.87 = 2.55 W/ m2 oC
0.387 𝑅𝑅𝑅𝑅 1/6
[1+(0.492/𝑃𝑃𝑃𝑃 )9/16 ]8/27 0.387 (3.51𝑥𝑥10 9 )1/6
]1/2 1/2
] [1+(0.492/0.7078)9/16 ]8/27
43
b) Fluks panas 𝑞𝑞
𝐴𝐴
= (I x τα) – h ( Td – T∞ ) = (462 x 0.61) – 2.55 (30.8 – 25) = 267 W/m2
5) HE ruang pengering sebelah kanan Suhu he = 31.1 oC Suhu inlet = 24.9 oC 31.1 + 24.9 Suhu fluida operasi = = 28 oC = 301 K 2 Suhu fluida HE = 31.1 + 273 = 304.1 K Suhu fluida ruang pengering = 24.9 + 273 = 297.9 K Sifat Massa jenis Massa jenis (ρ1) Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl Bilangan Prandtl (Prs)
Nilai 1.1738 kg/m3 1.1873 kg/m3 1.0058 kj/kg oC 0.02632 w/ oC 1.851x10-5 kg/ms 15.79 x10-6 m2/s 0.708 0.707
a) Koefisien kehilangan (k) Penurunan tekanan (ΔP) dihitung dari daya 2 kipas, atas dan bawah (P = 80 W) Ainlet = 0.29 m x 0.29 m = 0.0841 m2 ΔP untuk 1 buah inlet = =
k= =
2 x ΔP
Pη
v inlet x A inlet 80 x 0.7
3 x 0.0841
= 221.96
ρ x v 2kipas
2 x 221.96
1.1738 x 32
= 42.02
b) Koefisien pindah panas konveksi (h) 31.1 + 24.9 Suhu fluida operasi = = 28 oC = 301 K ρf =
p
RT
2
44 101325
= 287 x 301 = 1.173 kg/m3 * ST = 11.4063 cm , L = 187 cm * SL = 6.7501 cm , D = 3.25 cm * SD = �𝑆𝑆𝑆𝑆2 + (𝑆𝑆𝑆𝑆/2)2 = �6.75012 + (11.4063/2)2 = 8.837 cm AD = (SD – D) L = (8.837 – 3.25) x 187 = 1044.739 cm2 AT = (ST – D) L = (6.7501 – 3.25) x 187 = 1525.228 cm2 2 AD > AT 𝑆𝑆𝑆𝑆 v maks = 𝑣𝑣
𝑆𝑆𝑆𝑆−𝐷𝐷 11.4063
= 𝑥𝑥 3 11.4063−0.0325 = 4.195 m/s
Re =
ρf x v maks x D μ 1.173 x 4.195 x 0.0325
= 1.851 x 10 −5 = 8641.1
1000 < Re > 2x105 Nu = 0.35 (ST/SL)0.2 (Re)0.6 (Pr)0.36 (Pr/Prs)0.25 = 0.35 (11.4063/6.7501)0.2 (8641.1)0.6 (0.708)0.36 (0.708/0.707)0.25 = 79 Karena, ada 2 baris. Maka nilai Nu dikalikan faktor F = 0.76 NuD = F x Nu = 0.76 x 79 = 60.04 h=
𝑁𝑁𝑁𝑁 𝐷𝐷 𝑥𝑥 𝑘𝑘
𝐷𝐷 60.04 𝑥𝑥 0.02632
= 0.0325 = 48.62 w/m2 oC c) Fluks panas
45 Total dari HE adalah 17 buah. Untuk satu unit panjang tabung (L = 1.87 m), area permukaan pindah panas dan laju aliran udara adalah A=NπDL = 17 π (0.0325)(1.87) = 3.25 m2 ṁ = ρ1 V (n SL L) = 1.1873 x 3 x (2 x 0.067501 x 1.87) = 0.90 Te = Ts – (Ts – Ti) exp �–
AS h
ṁ Cp
�
= 31.1 – (31.1 – 24.9) exp �– = 25.89
ΔTln = =
3.25 x 48.62
0.90 x 1005.8
�
�𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑒𝑒 � − �𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑖𝑖 �
ln ��𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑒𝑒 �/�𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑖𝑖 �� �31.1 – 25.89� − �31.1 – 24.9�
ln ��31.1 – 25.89�/�31.1 – 24.9��
= 5.69 oC Q/A = h ΔTln = 48.62 x 5.69 = 276.60 W/m2
5) HE ruang pengering sebelah kiri Suhu he = 30.8 oC Suhu inlet = 24.9 oC 30.8 + 24.9 Suhu fluida operasi = = 27.85 oC = 300.85 K 2 Suhu fluida HE = 30.8 + 273 = 303.8 K Suhu fluida ruang pengering = 24.9 + 273 = 297.9 K Sifat Massa jenis Massa jenis (ρ1) Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl Bilangan Prandtl (Prs)
Nilai 1.1744 kg/m3 1.1873 kg/m3 1.0058 kj/kg oC 0.02630 w/ oC 1.850x10-5 kg/ms 15.78 x10-6 m2/s 0.708 0.707
a) Koefisien kehilangan (k) Penurunan tekanan (ΔP) dihitung dari daya 2 kipas, atas dan bawah (P = 80 W) Ainlet = 0.29 m x 0.29 m = 0.0841 m2 ΔP untuk 1 buah inlet =
Pη
v inlet x A inlet
46 =
2 x ΔP
k=
80 x 0.7
3 x 0.0841
= 221.96
ρ x v 2kipas
2 x 221.96
=
1.1744 x 32
= 42
b) Koefisien pindah panas konveksi (h) 30.8 + 24.9 Suhu fluida operasi = = 27.85 oC = 300.85 K ρf =
p
RT 101325
2
= 287 x 300.85 = 1.174 kg/m3 * ST = 11.4063 cm , L = 187 cm * SL = 6.7501 cm , D = 3.25 cm * SD = �𝑆𝑆𝑆𝑆2 + (𝑆𝑆𝑆𝑆/2)2 = �6.75012 + (11.4063/2)2 = 8.837 cm AD = (SD – D) L = (8.837 – 3.25) x 187 = 1044.739 cm2 AT = (ST – D) L = (6.7501 – 3.25) x 187 = 1525.228 cm2 2 AD > AT 𝑆𝑆𝑆𝑆 v maks = 𝑣𝑣
𝑆𝑆𝑆𝑆−𝐷𝐷 11.4063
= 𝑥𝑥 3 11.4063−0.0325 = 4.195 m/s
Re =
ρf x v maks x D
μ 1.174 x 4.195 x 0.0325
= 1.850 x 10 −5 = 8648.65
1000 < Re > 2x105 Nu = 0.35 (ST/SL)0.2 (Re)0.6 (Pr)0.36 (Pr/Prs)0.25 = 0.35 (11.4063/6.7501)0.2 (8648.65)0.6 (0.708)0.36 (0.708/0.707)0.25
47 = 79.04 Karena, ada 2 baris. Maka nilai Nu dikalikan faktor F = 0.76 NuD = F x Nu = 0.76 x 79.04 = 60.07 h=
𝑁𝑁𝑁𝑁 𝐷𝐷 𝑥𝑥 𝑘𝑘
𝐷𝐷 60.07𝑥𝑥 0.02630
= 0.0325 = 48.62 w/m2 oC c) Fluks panas Total dari HE adalah 17 buah. Untuk satu unit panjang tabung (L = 1.87 m), area permukaan pindah panas dan laju aliran udara adalah A=NπDL = 17 π (0.0325)(1.87) = 3.25 m2 ṁ = ρ1 V (n SL L) = 1.1873 x 3 x (2 x 0.067501 x 1.87) = 0.90 Te = Ts – (Ts – Ti) exp �–
AS h
ṁ Cp
�
= 30.8 – (30.8 – 24.9) exp �– = 25.84
ΔTln = =
3.25 x 48.62
0.90 x 1005.8
�
�𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑒𝑒 � − �𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑖𝑖 � ln ��𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑒𝑒 �/�𝑇𝑇𝑆𝑆 – 𝑇𝑇𝑖𝑖 �� �30.8 – 25.84� − �31.1 – 24.9�
ln ��30.8 – 25.84�/�31.1 – 24.9��
= 5.41 oC Q/A = h ΔTln = 48.62 x 5.41 = 263.23 W/m2
6) Perhitungan parameter kondisi batas pada rak dianggap sebagai porous jump. Rak dengan beban a) Permeabilitas permukaan α= =
𝐷𝐷𝑝𝑝2
𝜀𝜀 2
150 (1−𝜀𝜀)2 0.0298 2 150
0.838 2
(1−0.838)2
48 = 1.326 x 10-4 b) Koefisien porous jump 3.5 (1−𝜀𝜀) C2 = 3 𝐷𝐷𝐷𝐷 𝜀𝜀
=
3.5 (1−0.838)
0.0298 𝑥𝑥 (0.838)3
= 32.352
49 Lampiran 2 Tabel validasi suhu pada pukul 08.00 kondisi tanpa beban Rak 1BI 9BI 10 D I 17 B I
Suhu (Simulasi) o C 33.61 34.34 34.37 37.07
1 D II 9 D II 10 B II 17 D II
33.77 34.15 34.20 37.30
1 D III 9 B III 10 D III 17 B III
32.53 31.91 31.89 30.97
1 B IV 8 B IV 9 D IV 17 D IV
32.92 31.65 31.63 30.94
Suhu (Pengukuran) o C 44.2 36.8 35.7 36.9 Error rata-rata Error maksimum Error minimum Error standar deviasi 49.3 30.4 31.2 30.4 Error rata-rata Error maksimum Error minimum Error standar deviasi 31.5 32.8 30 30.2 Error rata-rata Error maksimum Error minimum Error standar deviasi 31.9 32.7 31.8 31.1 Error rata-rata Error maksimum Error minimum Error standar deviasi
Error rata-rata secara keseluruhan Error maksimum secara keseluruhan Error minimum secara keseluruhan Error standar deviasi secara keseluruhan
Error 23.95 6.69 3.72 0.45 8.70 23.95 0.45 10.48 31.50 12.34 9.63 22.69 19.04 31.50 9.63 10.03 3.29 2.72 6.29 2.54 3.71 6.29 2.54 1.75 3.19 3.22 0.54 0.52 1.87 3.22 0.52 1.54 8.33 31.50 0.45 9.52
50 Lampiran 3 Tabel validasi suhu pada pukul 12.00 kondisi tanpa beban Rak 1BI 9BI 10 D I 17 B I
1 D II 9 D II 10 B II 17 D II
1 D III 9 B III 10 D III 17 B III
1 B IV 8 B IV 9 D IV 17 D IV
Suhu (Simulasi) o C 48.31 42.59 42.51 42.73
Suhu (Pengukuran) o C 62 43.1 43.9 44.4 Error rata-rata Error maksimum Error minimum Error standar deviasi 46.89 57.2 42.74 41.5 42.64 42.2 42.79 41.3 Error rata-rata Error maksimum Error minimum Error standar deviasi 45.44 58.2 42.02 41.7 41.95 41 42.42 41.5 Error rata-rata Error maksimum Error minimum Error standar deviasi 46.69 55.5 41.99 46.2 41.88 43.8 42.33 41.9 Error rata-rata Error maksimum Error minimum Error standar deviasi Error rata-rata secara keseluruhan Error maksimum secara keseluruhan Error minimum secara keseluruhan Error standar deviasi secara keseluruhan
Error 22.07 1.18 3.16 3.77 7.54 22.07 1.18 9.75 18.02 2.99 1.05 3.61 6.42 18.02 1.05 7.81 21.93 0.77 2.32 2.21 6.81 21.93 0.77 10.11 15.87 9.12 4.38 1.02 7.60 15.87 1.02 6.44 7.09 22.07 0.77 7.76
51 Lampiran 4 Tabel validasi suhu pada pukul 16.00 kondisi tanpa beban Rak 1BI 9BI 10 D I 17 B I
1 D II 9 D II 10 B II 17 D II
1 D III 9 B III 10 D III 17 B III
1 B IV 8 B IV 9 D IV 17 D IV
Suhu (Simulasi) o C 39.40 39.52 39.52 39.26
Suhu (Pengukuran) o C 45.9 42 39.8 38.8 Error rata-rata Error maksimum Error minimum Error standar deviasi 39.44 45 39.51 41 39.50 42.4 39.22 38.3 Error rata-rata Error maksimum Error minimum Error standar deviasi 34.63 49.9 33.72 41.9 33.81 41.6 37.35 39.9 Error rata-rata Error maksimum Error minimum Error standar deviasi 35.04 49.6 33.50 43.8 33.52 42.2 36.50 38.3 Error rata-rata Error maksimum Error minimum Error standar deviasi Error rata-rata secara keseluruhan Error maksimum secara keseluruhan Error minimum secara keseluruhan Error standar deviasi secara keseluruhan
Error 14.15 5.90 0.71 1.18 5.49 14.15 0.71 6.23 12.37 3.64 6.84 2.41 6.31 12.37 2.41 4.45 30.60 19.52 18.73 6.39 18.81 30.60 6.39 9.90 29.35 23.52 20.57 4.71 19.54 29.35 4.71 10.54 12.54 30.60 0.71 10.03
52 Lampiran 5 Tabel validasi suhu pada pukul 12.00 kondisi dengan beban Rak 1DI 8DI 17 D I
1 D II 8 D II 17 D II
1 D III 8 B III 16 BIII
1 D IV 8 D IV 16 D IV
Suhu (Simulasi) o C 42.57 39.49 40.90
Suhu (Pengukuran) o C 45.9 42 39.8 Error rata-rata Error maksimum Error minimum Error standar deviasi 41.78 45 39.63 41 41.11 42.4 Error rata-rata Error maksimum Error minimum Error standar deviasi 44.04 49.9 41.22 41.9 40.42 41.6 Error rata-rata Error maksimum Error minimum Error standar deviasi 45.09 49.6 41.14 43.8 40.27 42.2 Error rata-rata Error maksimum Error minimum Error standar deviasi Error rata-rata secara keseluruhan Error maksimum secara keseluruhan Error minimum secara keseluruhan Error standar deviasi secara keseluruhan
Error 7.25 5.98 2.77 5.33 7.25 2.77 2.31 7.15 3.33 3.03 4.50 7.15 3.03 2.30 11.75 1.62 2.85 5.41 11.75 1.62 5.53 9.09 6.06 4.57 6.57 9.09 4.57 2.30 5.34 11.75 1.62 2.62
53 Lampiran 6 Data intensitas cahaya matahari Waktu (jam) 8:00 8:30 9:00 9:30 10:00 10:30 11:00 11:30 12:00 12:30 13:00 13:30 14:00 14:30 15:00 15:30 16:00 16:30 17:00
P1 (W/m2) 142.86 328.57 257.14 385.71 671.43 714.29 642.86 600.00 800.00 414.29 314.29 914.29 542.86 800.00 485.71 585.71 614.29 85.71 28.57
P2 (W/m2) 200.00 185.71 228.57 242.86 357.14 328.57 214.29 614.29 814.29 800.00 985.71 485.71 928.57 942.86 928.57 942.86 585.71 242.86 85.71
54 Lampiran 7 Data input CFD (Boundary condition) pukul 08.00 kondisi tanpa beban Dinding (Wall) Kanan
Kiri Satuan Parameter W/m2 Heat flux m Wall thickness o C Temperature W/m3 Heat generation rate W/moC Heat transfer coefficient o C Free stream temperature Atap (Wall) Kanan Kiri Parameter Nilai Satuan Parameter 264.37 W/m2 Heat flux Heat flux 0.001 m Wall thickness Wall thickness o 31.1 C Temperature Temperature 0 W/m3 Heat generation rate Heat generation rate 2.86 W/moC Heat transfer coefficient Heat transfer coefficient o 25 C Free stream temperature Free stream temperature Lantai (Wall) Kanan Kiri Parameter Nilai Satuan Parameter 49.3 W/m2 Heat flux Heat flux 0.001 m Wall thickness Wall thickness o 38.9 C Temperature Temperature 0 W/m3 Heat generation rate Heat generation rate 3.55 W/moC Heat transfer coefficient Heat transfer coefficient o 25 C Free stream temperature Free stream temperature Heat exchanger (Radiator) Kanan Kiri Parameter Nilai Satuan Parameter 276.60 W/m2 Heat flux Heat flux o 31.1 C Temperature Temperature 42.02 Loss coefficient Loss coefficient 48.62 W/moC Heat transfer coefficient Heat transfer coefficient Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature
Nilai 253.2 0.001 34.6 0 2.98 25
Outlet (Outflow) Parameter Kanan 0.1 Flow rate weighting Rak (Porous jump) Parameter Nilai 0 Face permeability 0 Porous medium thickness 0 Pressure jump coefficient
kiri 0.1 Satuan m2 m (C2) (1/m)
Nilai 267 0.001 30.8 0 2.55 25
Satuan W/m2 m o C W/m3 W/moC o C
Nilai 240.89 0.001 37.1 0 3.38 25
Satuan W/m2 m o C W/m3 W/moC o C
Nilai 15.55 0.001 30.5 0 2.83 25
Satuan W/m2 m o C W/m3 W/moC o C
Nilai 263.23 30.8 42 48.62
Satuan W/m2 o C W/moC
Inlet (Velocity inlet) Parameter Nilai Satuan 3 m/s Velocity o 24.9 C Temperature
55 Lampiran 8 Data input CFD (Boundary condition) pukul 12.00 kondisi tanpa beban Dinding (Wall) Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Kanan Parameter Heat flux Temperature Loss coefficient Heat transfer coefficient
Kiri Satuan Parameter Nilai 2 W/m 454.2 Heat flux m 0.001 Wall thickness o C 41.4 Temperature W/m3 Heat generation rate 0 W/moC Heat transfer coefficient 1.48 o C 42.5 Free stream temperature Atap (Wall) Kiri Nilai Satuan Parameter Nilai 2 372.95 W/m 391.78 Heat flux 0.001 m 0.001 Wall thickness o 63.5 C 59.4 Temperature 0 W/m3 Heat generation rate 0 3.79 W/moC Heat transfer coefficient 3.6 o 42.5 C 42.5 Free stream temperature Lantai (Wall) Kiri Nilai Satuan Parameter Nilai 2 0.58 W/m 0.24 Heat flux 0.0005 m 0.0005 Wall thickness o 42.9 C 42.7 Temperature 0 W/m3 Heat generation rate 0 1.44 W/moC Heat transfer coefficient 1.21 o 42.5 C 42.5 Free stream temperature Heat exchanger (Radiator) Kiri Nilai Satuan Parameter Nilai 2 411.70 W/m 352.13 Heat flux o 40.6 C 39.6 Temperature 16.75 16.72 Loss coefficient 64.21 W/moC Heat transfer coefficient 64.23 Nilai 453.1 0.001 42.1 0 1.09 42.5
Outlet (Outflow) Parameter Kanan kiri 0.1 0.1 Flow rate weighting Rak (Porous jump) Parameter Nilai Satuan 0 m2 Face permeability 0 m Porous medium thickness 0 (C2) (1/m) Pressure jump coefficient
Inlet (Velocity inlet) Parameter Nilai Satuan 4.82 m/s Velocity o C Temperature 33.7
Satuan W/m2 m o C W/m3 W/moC o C
Satuan W/m2 m o C W/m3 W/moC o C
Satuan W/m2 m o C W/m3 W/moC o C
Satuan W/m2 o C W/moC
56 Lampiran 9 Data input CFD (Boundary condition) pukul 16.00 kondisi tanpa beban Dinding (Wall) Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Kanan Parameter Heat flux Temperature Loss coefficient Heat transfer coefficient
Kiri Satuan Parameter Nilai 2 W/m 103.4 Heat flux m 0.001 Wall thickness o C 39.3 Temperature W/m3 Heat generation rate 0 W/moC Heat transfer coefficient 2.57 o C 33 Free stream temperature Atap (Wall) Kiri Nilai Satuan Parameter Nilai 2 102.98 W/m 101.21 Heat flux 0.001 m 0.001 Wall thickness o 38.9 C 39.4 Temperature 0 W/m3 Heat generation rate 0 2.81 W/moC Heat transfer coefficient 2.87 o 33 C 33 Free stream temperature Lantai (Wall) Kiri Nilai Satuan Parameter Nilai 2 16.81 W/m 18.61 Heat flux 0.0005 m 0.0005 Wall thickness o 38.9 C 39.4 Temperature 0 W/m3 Heat generation rate 0 2.85 W/moC Heat transfer coefficient 2.91 o 33 C 33 Free stream temperature Heat exchanger (Radiator) Kiri Nilai Satuan Parameter Nilai 2 459.25 W/m 500.42 Heat flux o 40.6 C 41.4 Temperature 26.58 26.61 Loss coefficient 55.89 W/moC Heat transfer coefficient 55.87 Nilai 103 0.001 39.4 0 2.58 33
Outlet (Outflow) Parameter Kanan 0.1 Flow rate weighting Rak (Porous jump) Parameter Nilai 0 Face permeability 0 Porous medium thickness 0 Pressure jump coefficient
kiri 0.1 Satuan m2 m (C2) (1/m)
Inlet (Velocity inlet) Parameter Nilai Satuan 3.82 m/s Velocity o 31.7 C Temperature
Satuan W/m2 m o C W/m3 W/moC o C
Satuan W/m2 m o C W/m3 W/moC o C
Satuan W/m2 m o C W/m3 W/moC o C
Satuan W/m2 o C W/moC
57
Lampiran 10 Data input CFD (Boundary condition) pukul 12.00 kondisi dengan beban Dinding (Wall) Kanan Parameter
Kiri Nilai 263.3 0.001 39 0 1.45 38
Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature
Satuan W/m2 m o C W/m3 W/moC o C
Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature Atap (Wall)
Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature
Nilai 220.63 0.001 51 0 3.39 38
Satuan Parameter 2 W/m Heat flux m Wall thickness o C Temperature W/m3 Heat generation rate W/moC Heat transfer coefficient o C Free stream temperature Lantai (Wall)
Kanan Parameter Heat flux Wall thickness Temperature Heat generation rate Heat transfer coefficient Free stream temperature
Nilai 10.16 0.0005 42 0 2.54 38
Satuan Parameter 2 W/m Heat flux m Wall thickness o C Temperature W/m3 Heat generation rate W/moC Heat transfer coefficient o C Free stream temperature Heat exchanger (Radiator)
Kanan Parameter Heat flux Temperature Loss coefficient Heat transfer coefficient
Nilai 282.63 42.4 43.51 48.24
Satuan W/m2 o C W/moC
Outlet (Outflow) Kanan 0.1 Flow rate weighting Rak (Porous jump) Parameter Nilai 1.326 x 10-4 Face permeability 0.00122 Porous medium thickness 32.352 Pressure jump coefficient
Parameter Heat flux Temperature Loss coefficient Heat transfer coefficient
kiri 0.1 Satuan m2 m (C2) (1/m)
Nilai 259.4 0.001 40.7 0 1.96 38
Satuan W/m2 m o C W/m3 W/moC o C
Kiri Nilai 200.47 0.001 55.6 0 3.65 38
Satuan W/m2 m o C W/m3 W/moC o C
Kiri Nilai 3.49 0.0005 39.7 0 2.05 38
Satuan W/m2 m o C W/m3 W/moC o C
Kiri Nilai 106.10 38.4 43.24 48.3
Satuan W/m2 o C W/moC
Inlet (Velocity inlet) Parameter Nilai Satuan 3 m/s Velocity o C Temperature 36
58
Lampiran 11 Tabel ΔP, efisiensi, daya kipas dan keseragaman suhu modifikasi 3 ΔP
4.905 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 50 100
Efisiensi Daya kipas
0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7 0.3 0.7
0.0428 0.0183 0.0436 0.0187 0.0524 0.0224 0.0611 0.0262 0.0698 0.0299 0.0785 0.0337 0.0873 0.0374 0.0960 0.0411 0.1047 0.0449 0.1134 0.0486 0.1222 0.0524 0.1309 0.0561 0.1396 0.0598 0.1484 0.0636 0.1571 0.0673 0.1658 0.0711 0.1745 0.0748 0.4363 0.1870 0.8727 0.3740
Ruang pengering kiri Rata- Standar CV rata deviasi
Ruang pengering kanan RataStandar CV rata deviasi
Keseluruhan RataStandar CV rata deviasi
42.33
0.56
1.33
41.61
0.31
0.75
41.97
0.58
1.38
42.27
0.5
1.18
41.55
0.35
0.83
41.91
0.56
1.33
42.25
0.47
1.1
41.54
0.27
0.66
41.9
0.52
1.24
42.13
0.42
1.01
41.59
0.24
0.57
41.86
0.44
1.04
42.04
0.47
1.13
41.61
0.22
0.53
41.83
0.43
1.02
42.1
0.46
1.1
41.62
0.18
0.44
41.86
0.43
1.02
42.03
0.37
0.89
41.66
0.17
0.41
41.84
0.35
0.83
42.04
0.27
0.63
41.66
0.16
0.4
41.85
0.29
0.69
42.02
0.27
0.63
41.66
0.15
0.36
41.84
0.28
0.67
42.01
0.25
0.58
41.69
0.15
0.36
41.85
0.26
0.62
42.36
0.14
0.32
41.58
0.13
0.31
41.97
0.41
0.98
41.95
0.27
0.63
41.69
0.19
0.45
41.82
0.26
0.63
41.97
0.27
0.64
41.71
0.13
0.31
41.84
0.25
0.59
41.96
0.26
0.61
41.7
0.13
0.31
41.83
0.24
0.57
41.94
0.23
0.54
41.73
0.14
0.33
41.84
0.22
0.51
41.98
0.19
0.44
41.74
0.12
0.29
41.86
0.2
0.47
41.95
0.26
0.62
41.75
0.09
0.23
41.85
0.22
0.52
42.65
0.08
0.18
41.1
0.05
0.13
41.88
0.78
1.86
42.64
0.06
0.14
41.07
0.04
0.11
41.86
0.79
1.89
59
60
61 Lampiran 14 Distribusi suhu rak pada pengukuran 1 Rak Standar deviasi (oC) Rata-rata (oC) Maksimum (oC) Minimum (oC) ID 3 4.66 43.01 52.20 35.00 5 3.46 40.80 47.00 34.80 7 3.27 40.63 46.10 35.40 9 3.48 38.45 43.60 31.20 11 2.94 37.68 43.20 33.20 13 2.68 36.34 41.80 31.20 15 2.69 35.07 40.70 30.00 17 2.31 33.59 37.70 29.50 IIB 3 4.35 42.93 49.40 34.70 5 3.84 40.62 47.30 32.10 7 4.45 40.16 47.80 31.70 9 3.88 38.74 45.60 30.80 11 4.64 38.34 46.70 30.00 13 3.11 34.53 43.70 29.90 15 3.50 36.61 42.80 30.80 17 2.74 33.62 40.10 28.10 IIID 1 8.82 49.77 64.40 36.20 3 5.73 44.28 54.80 34.10 5 4.76 41.39 50.00 33.20 7 4.58 39.32 46.70 30.30 9 4.27 37.90 44.90 29.60 11 3.92 36.42 42.30 29.00 13 3.22 35.72 41.50 30.20 15 3.43 35.00 41.90 29.20 17 2.86 35.17 39.80 30.10 IVB 1 8.81 52.65 72.20 36.60 3 6.77 45.42 56.30 32.00 5 7.17 45.39 56.90 31.80 7 7.06 42.64 54.10 31.20 9 4.76 41.97 50.30 35.00 11 5.13 40.60 48.60 30.60 13 5.00 40.06 47.80 30.20 15 4.34 38.78 46.30 30.50 17 4.11 38.14 44.90 31.00
62
63
64 Lampiran 17 Distribusi suhu rak pada pengukuran 2 Rak Standar deviasi (oC) Rata-rata (oC) Maksimum (oC) Minimum (oC) IB 2 6.09 44.47 52.80 35.90 4 2.55 32.95 37.40 29.80 6 4.13 42.45 48.90 35.40 8 3.80 41.44 46.60 34.30 10 3.58 41.17 46.20 34.80 12 3.80 40.80 45.80 33.40 14 3.66 39.73 44.10 32.60 16 3.27 39.52 45.10 31.80 IID 2 6.34 42.79 53.60 34.10 4 5.74 40.92 49.40 32.10 6 4.15 38.37 44.40 32.10 8 3.23 38.58 43.70 33.90 10 3.32 38.47 43.30 32.70 12 3.10 38.18 43.00 33.60 14 2.62 37.91 41.80 33.80 16 2.39 38.06 44.60 33.70 IIIB 2 7.31 44.25 55.50 33.10 4 6.54 42.22 52.10 32.50 6 6.35 41.39 50.80 30.30 8 6.01 40.56 48.30 29.20 10 5.30 39.45 45.80 29.00 12 5.45 39.77 46.60 30.80 14 5.15 39.16 45.40 30.10 16 3.15 33.02 38.30 27.20 IVD 4 6.33 43.65 52.20 34.60 6 6.48 40.96 49.80 32.00 8 6.32 39.88 48.10 30.10 10 5.28 39.51 47.30 32.20 12 5.44 38.98 46.90 30.80 14 4.71 37.07 44.60 29.60 16 3.10 33.94 38.90 30.10
65
66
67 Lampiran 22 Tahapan penggunaan ANSYS Fluent untuk simulasi Menjalankan ANSYS Fluent Menjalankan program Fluent, dengan cara memilih ikon Fluent pada komputer. Memilih 3D Dimension, lalu mencentang seluruh pilihan yang ada di Display Option dan memilih Serial Processing.
Tampilan awal Fluent (Fluent Launcher) Pengaturan Umum (General Setting) Proses pengaturan diawali dengan membuka file yang ingin disimulasikan. File yang dibuka adalah model geometri yang telah di meshing pada tahap prapemrosesan. Melakukan pengaturan satuan unit yang diinginkan. Model yang telah dibuka, dilakukan pengecekan model apakah model geometri tersebut memiliki kesalahan. Salah satu indikator model telah benar adalah total volume bernilai positif. Pengecekan model telah selesai, dilakukan tahapan pengaturan skala model ke dalam satuan meter.
Tampilan membuka file model geometri (read mesh)
68
Tampilan pengaturan unit
Tampilan pengecekan grid/mesh
Tampilan pengaturan skala Pemilihan Model (Models) Pemilihan model memiliki beberapa jenis, tergantung pada proses atau aliran. Model energy dipilih, karena model yang akan disimulasikan memerlukan analisis tentang temperatur dan perpindahan panas. Model viscous yang dipilih adalah laminar, karena aliran fuida yang terjadi dalam ruang pengering ERK adalah aliran laminar. Aliran fluida dalam pengering ERK dianggap laminar jika RE < 2000 (Wulandani, 2005).
69
Tampilan pemilihan model energy
Tampilan pemilihan model viskos Penambahan Material (Materials) Penambahan material yang diperlukan adalah material dinding pada pengering ERK yaitu polikarbonat. Hal ini dikarenakan dalam database Fluent tidak terdapat karakteristik material polikarbonat. Properti fisik yang diperlukan dalam menambahkan jenis material polikarbonat adalah kerapatan, panas jenis dan konduktivitas. Properti fisik polikarbonat Sifat fisik Densitas Panas jenis Konduktivitas
Nilai 1210 kg/m3 1250 J/kg.oC 0.228 W/m.oC
Sumber (Mark, 2007) (Mark, 2007) Pengukuran
70
Tampilan penambahan material baru Pengaturan Kondisi Batas (Boundary Conditions) Kondisi batas/boundary conditions menentukan suatu aliran dan variable panas pada batas-batas model fisik/geometri yang telah dibuat. Kondisi batas yang terdapat pada model pengering ERK ini adalah velocity inlet, outflow, porous jump, radiator, wall.
Tampilan pengaturan kondisi batas velocity inlet
Tampilan pengaturan kondisi batas outflow
71
Tampilan pengaturan kondisi batas porous jump.
Tampilan pengaturan kondisi batas radiator.
Tampilan pengaturan kondisi batas wall. Pengaturan Metode Solusi (Solution Methods) dan Kontrol Solusi (Solution Controls) Pengaturan metode solusi merupakan proses interpolasi perhitungan pada titik-titik simpul grid untuk mendapatkan nilai yang kontinu. Skema interpolasi yang digunakan adalah SIMPLE. Pengaturan parameter kontrol solusi adalah faktor under-relaxation untuk menstabilkan proses iterasi.
72
Tampilan pengaturan Solution Methods
Tampilan pengaturan Solution Controls. Pengaturan Inisialisasi (Initialization) Proses inisialisasi (tebakan awal) dihitung dari semua kondisi batas (all zones) dengan nilai awal sesuai dengan nilai yang dimasukkan pada kondisi batas tersebut.
Tampilan pengaturan Solution Initialization
73 Pengaturan Perhitungan Simulasi (Calculation Activities) Pengaturan perhitungan meliputi kriteria konvergensi dan banyaknya iterasi yang akan dilakukan. Kriteria konvergensi menggunakan nilai 0.001 untuk semua persamaan. Nilai 0.001 m/detik untuk persamaan kecepatan, 0.001 Pa untuk persamaan tekanan dan 10-6 pada persamaan energi. Kriteria konvergensi adalah perbedaan antara tebakan awal dan hasil akhir iterasi. Pengaturan perhitungan juga dapat digunakan untuk mengatur jumlah iterasi yang akan disimpan dan banyaknya iterasi yang akan disimpan secara otomatis. Hal ini berguna untuk mencegah data iterasi yang hilang apabila komputer mati atau rusak secara tiba-tiba.
Tampilan pengaturan Calculation Activities Proses Perhitungan (Run Calculation) Tahapan dimana perhitungan semua persamaan, batas-batas kondisi dari model yang sudah ditentukan dan berhenti sesuai jumlah iterasi yang diinginkan atau setelah mencapai kondisi konvergen.
Tampilan Run Calculation
74 Hasil (Results) Tahapan ini adalah tahapan dimana hasil dari perhitungan ditampilkan dalam bentuk grafik dan animasi, plot data dan berupa laporan/reports. Pengaturan hasil berupa grafik dan animasi meliputi, tampilan geometri dalam bentuk grid/mesh yang terbentuk, tampilan hasil simulasi ke dalam bentuk kontur.
Tampilan pengaturan Results Panduan lengkap tentang penggunaan program Fluent dapat dilihat pada menu help, pilih konten PDF lalu pilih file user guide atau tutorial guide.
75
RIWAYAT HIDUP Penulis dilahirkan di Jakarta pada tanggal 15 April 1991 pasangan Amelar Siregar (Almarhum) dan Rosmian Purba. Penulis adalah putra ketiga dari empat bersaudara. Tahun 2009 penulis lulus dari SMAN 58 Ciracas, Jakarta Timur dan pada tahun yang sama penulis lulus seleksi masuk Institut Pertanian Bogor melalui jalur SNMPTN dan diterima di Departemen Teknik Pertanian, Fakultas Teknologi Pertanian. Selama mengikuti perkualiahan, penulis aktif sebagai anggota Unit Kegiatan Mahasiswa Persekutuan Mahasiswa Kristen (PMK). Penulis juga menjadi koordinator di Komisi Literatur Persekutuan Mahasiswa Kristen IPB tahun 2012, dan menjadi anggota Litbang PMK IPB tahun 2013. Penulis juga pernah mengikuti seminar EBTKE CONEX pada tahun 2013. Penulis melaksanakan praktik lapangan pada bulan Juli – Agustus 2012 di PT. Mekar Unggul Sari, Bogor dengan judul Mempelajari Aspek Keteknikan Pertanian dan Pemanfaatan Energi pada Budidaya Melon dengan Sistem Hidroponik di PT. Mekar Unggul Sari, Bogor.