PENENTUAN PUTARAN PADA PENGERING SURYA TIPE RAK BERPUTAR DAN PEMODELAN PENGERING DENGAN METODE COMPUTATIONAL FLUID DYNAMICS SUERITAH HENNY WANTI

PENENTUAN PUTARAN PADA PENGERING SURYA TIPE RAK BERPUTAR DAN PEMODELAN PENGERING DENGAN METODE COMPUTATIONAL FLUID DYNAMICS

SUERITAH HENNY WANTI

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 Penentuan Putaran Pada Pengering Surya Tipe Rak Berputar dan Pemodelan Pengering dengan metode Computational Fluid Dynamics 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, April 2014 Sueritah Henny Wanti NIM F14090026

ABSTRAK SUERITAH HENNY WANTI. Penentuan Putaran pada Pengering Surya Tipe Rak Berputar dan Pemodelan Pengering dengan Metode Computational Fluid Dynamics. Dibimbing oleh DYAH WULANDANI. Pengering surya 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. Sumber energi pada pengering ERK diperoleh dari surya dan pembakaran biomasa. Masalah ketidakmerataan suhu rak di dalam ruang pengering, masih selalu terjadi. Pemecahan masalah ketidakmerataan aliran udara panas dapat dilakukan dengan memutar rak pada tipe rak berputar. Tujuan penelitian ini adalah untuk menentukan putaran rak yang tepat dan membuat pemodelan simulasi CFD untuk menggambarkan sebaran suhu pada rak pengering. Penentuan putaran rak pada penelitian ini menggunakan 6 skenario. Hasil terbaik diperoleh pada skenario 5 pada putaran 135o dengan putaran satu jam sekali. Pengering pada skenario 5 ini memberikan hasil distribusi suhu rak yang seragam yang dinyatakan oleh nilai standar deviasi suhu terendah. Validasi suhu pengeringan hasil simulasi CFD yang dilakukan pada pukul delapan dan dua belas menunjukkan kecenderungan yang sama dengan suhu pengeringan hasil percobaan dengan eror 10%. Kata kunci: Pengering ERK, Rak berputar, Suhu rak, CFD.

ABSTRACT SUERITAH HENNY WANTI. Determination of rotation of the rotating rack Solar Dryer and Dryer Modelling by using Computational Fluid Dynamics Method. 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) that can be rotated vertically, and fans to discharge vapor evaporated from product The GHE solar dryer uses thermal energy source from the sun and biomass combustion. Un-uniformity of rack temperature problem can be solved by rotate the rack within the dryer. The aim of the study is to determine the properly of rack rotation and the temperature distribution by using Computational Fluid Dynamics (CFD) method. There are six scenarios of rack rotation condition to be simulated. The best result is scenario 5 at 135o with rotation per hour shown by the best uniformity of drying temperature indicated by the lowest temperature standar deviation. Validation of CFD result drying temperature whichis performed at eight and twelve o’clock without product, have shown similar trend to the experiment drying temperature distribution indicated by error value of 10%. Keywords: CFD, GHE, Rack Temperature

PENENTUAN PUTARAN PADA PENGERING SURYA TIPE RAK BERPUTAR DAN PEMODELAN PENGERING DENGAN METODE COMPUTATIONAL FLUID DYNAMICS

SUERITAH HENNY WANTI

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 : Penentuan Putaran pada Pengering Surya Tipe Rak Berputar dan Pemodelan Pengering dengan Metode Computational Fluid Dynamics Nama : Sueritah Henny Wanti NIM : F14090026

Disetujui oleh

Dr Ir Dyah Wulandani, MSi Pembimbing Akademik

Diketahui oleh

Dr Ir Desrial, MEng Ketua Departemen

Tanggal Lulus:

PRAKATA Puji dan syukur penulis panjatkan kepada Tuhan Yang Maha Kuasa atas segala karunia-Nya sehingga karya ilmiah ini berhasil diselesaikan. Tema yang dipilih dalam penelitian yang dilaksanakan sejak bulan Februari 2013 ini ialah pengering ERK, dengan judul Penentuan Putaran pada Pengering Surya Tipe Rak Berputar dan Pemodelan Pengering dengan Metode Computational Fluid Dynamics Dengan telah selesainya karya ilmiah ini, penulis ingin menyampaikan ucapan terima kasih yang sebesar-besarnya kepada: 1. Dr. Ir. Dyah Wulandani, M.Si selaku pembimbing yang memberikan bimbingan, arahan, motivasi kepada penulis. 2. Prof. Dr. Ir. Bambang Pramudya, M.Eng dan Dr. Muhamad Yuliantono, ST, MT selaku penguji terima kasih atas saran, bimbingan dan arahan, motivasi kepada penulis. 3. Pak Harto, Pak Darma, Mas Firman terima kasih atas bantuannya selama penelitian ini berlangsung. 4. Mama, Papa, Kakak (Sartika Mega S.Pd) adik (Robert Kenedy), abang (BRIPKA. Syah Rudi). 5. Terima kasih untuk hasian Benny Purba S.P. Untuk dukungan semangat dan saran mengerjakan skripsi, dan waktu untuk mendengarkan keluh kesah selama ini. 6. Teman-teman satu bimbingan (Elsamila, Stephani, Adit, dan Alfredo). 7. Teman-teman Vina Rondang, Gina Lupita, Gina Anisa, Tiara Etika, Riris, Raisa, Denny, Romi, Sumiharjon Simbolon, Heraldy Risva, Arnold, Aldha Hermiyanti, Ni’made Citta, Ni’putu dian, Eti, Yeti Aryani. 8. Teman-teman seperjuangan Teknik Mesin dan Biosistem angkatan 46 terima kasih atas kebersamaannya, bantuan dan semangatnya bagi penulis. Akhirnya penulis berharap semoga tulisan ini bermanfaat dan memberikan kontribusi nyata terhadap ilmu pengetahuan.

Bogor, April 2014 Sueritah Henny Wanti NIM F14090026

DAFTAR ISI DAFTAR TABEL

viii

DAFTAR GAMBAR

viii

DAFTAR LAMPIRAN

ix

PENDAHULUAN

1

Latar Belakang

1

Perumusan Masalah

1

Tujuan Penelitian

2

Manfaat Penelitian

2

Ruang Lingkup Penelitian

2

TINJAUAN PUSTAKA

2

Pengeringan

2

Pengering Efek Rumah Kaca Hybrid Tipe Rak Berputar

3

Computational Fluid Dynamics (CFD)

4

Proses Simulasi CFD

4

METODOLOGI PENELITIAN

5

Waktu dan Tempat Penelitian

5

Peralatan

5

Prosedur Penelitian

6

Persiapan Alat

6

Percobaan Pengeringan

7

Parameter yang Diukur pada Percobaan Pengeringan

7

Parameter yang diukur adalah sebagai berikut:

7

Pemodelan CFD

7

Perhitungan Parameter Kondisi Batas simulasi CFD

8

Pembuatan Model Geometri dan Pendefinisian Kondisi Batas Simulasi CFD 8 Simulasi CFD Udara Panas

9

Analisis Sebaran Suhu Hasil Simulasi CFD

9

Validasi Model Simulasi CFD

9

Posisi Titik Pengukuran ERK Tipe Rak Berputar HASIL DAN PEMBAHASAN

10 11

Suhu Rak Pengering

11

Hasil Pengukuran Skenario 1

11


13


14


16


17


18

Penentuan Putaran Rak Terbaik

20

Kelembaban Udara dan Iradiasi Surya pada Ruang Pengering

20

Pemodelan Simulasi CFD Suhu Rak Pengering

21

Bentuk Model Pengering dengan Gambit

21

Hasil Simulasi Pukul 08.00

21

Hasil Simulasi Pukul 12.00

25

Validasi Simulasi CFD

27

SIMPULAN DAN SARAN

29

Simpulan

29

Saran

29

DAFTAR PUSTAKA

29

RIWAYAT HIDUP

55

DAFTAR TABEL Tingkat keragaman suhu rak Skenario 1 Tingkat keragaman suhu rak Skenario 2 Tingkat keragaman suhu rak Skenario 3 Tingkat keragaman suhu rak Skenario 4 Tingkat keragaman suhu rak Skenario 5 Tingkat keragaman suhu rak Skenario 6 Distribusi suhu rak hasil simulasi CFD pukul 08.00 8 Distribusi suhu rak hasil simulasi CFD pukul 12.00 1 2 3 4 5 6 7

12 14 15 17 18 19 28 28

DAFTAR GAMBAR 1 2 3

Alat pengering Efek Rumah Kaca (ERK)-Hybrid tipe rak berputar Titik Pengukuran ERK Tipe Rak Berputar Diagram alir prosedur penelitian

3 10 11

Profil suhu rak Skenario 1 Profil suhu rak Skenario 2 Profil suhu rak Skenario 3 Profil suhu rak Skenario 4 Profil suhu rak Skenario 5 Profil suhu rak Skenario 6 Rata-rata standar deviasi masing-masing skenario Simulasi Skenario 1 Simulasi Skenario 2 Simulasi Skenario 3 Simulasi Skenario 4 Simulasi Skenario 5 Simulasi Skenario 6 Simulasi Skenario 1 Simulasi Skenario 2 Simulasi Skenario 3 Simulasi Skenario 4 Simulasi Skenario 5 Simulasi Skenario 6

4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

13 14 16 17 18 19 20 22 23 23 23 24 24 25 25 26 26 26 27

DAFTAR LAMPIRAN 1 2 3 4 5

Contoh analisis perhitungan parameter yang digunakan untuk simulasi CFD 31 Tabel distribusi suhu rak hasil simulasi CFD pada pukul 08.00 38 Tabel distribusi suhu rak hasil simulasi CFD pada pukul 12.00 44 Tabel kondisi batas (boundary conditions) pengukuran pukul 08.00 dan pukul 12.00 50 Tabel nilai kecepatan angin lingkungan, RH, Kecepatan angin outflow, Iradiasi matahari dan Suhu pengering pada masing-masing skenario. 52

DAFTAR SIMBOL Re ρ v D μ h Nu k l q/A I

Bilangan Reynold Massa jenis (kg/m3) Kecepatan (m/dtk) Diameter spesifik (m) Viskositas dinamik (kg/ms) Koefisien pindah panas konveksi (W/m2K) Bilangan Nusselt Konduktivitas panas (W/mK) Ketebalan dinding (m) Fluks panas (W/m2) Iradiasi surya (W/m2)

τα Td T∞ α Dp ε C2 k ΔP Gr g β T ν Ra Pr ṁ

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) Bilangan Grashoff Gaya gravitasi (m/s2) Koefisien (1/oC) Suhu (K) Viskositas kinematik (m2/s) Bilangan Rayleigh Bilangan Prandtl Laju aliran massa (kg/dtk)

PENDAHULUAN Latar Belakang Pengeringan merupakan operasi rumit yang memerlukan keseimbangan antara ketiga parameter yakni suhu, kecepatan aliran dan RH udara pengering. Dasar proses pengeringan adalah terjadinya proses penguapan air bahan ke udara karena perbedaan kandungan uap air antara udara dengan bahan yang dikeringkan. Penelitian tentang Efek Rumah Kaca (ERK) perlu mengacu pada dasar-dasar pindah panas, massa dan momentum, dikaitkan dengan pengetahuan tentang sifat bahan dan mutu (Brooker et al.1974). Seorang peneliti, selain mengandalkan pengetahuanya, maka percobaan di laboratorium dan berbagai percontohan yang dirangkai dengan pengalaman lapang sangat penting untuk mengembangkan dan penerapan suatu pengering untuk mencapai spesifikasi mutu yang lebih baik, laju produksi yang lebih tinggi, biaya energi yang lebih rendah dan kondisi lingkungan yang berubah. Usaha untuk melakukan proses pengeringan dapat pula dilakukan dengan menggunakan pengering Efek Rumah Kaca (ERK), salah satu keuntungan dengan menggunakan pengering ERK adalah produk lebih bersih dan tidak terkontaminasi dengan bahan yang tidak dikehendaki. Masalah ketidakseragaman kadar air produk hasil pengeringan juga merupakan akibat dari ketidakmerataan suhu rak di dalam ruang pengering, khususnya tipe rak bahwa perbedaan suhu rak atas lebih besar dibandingkan suhu rak bagian bawah. Perbedaan suhu rak ERK tipe rak berputar, suhu rak pada arah vertikal lebih besar dibandingkan suhu rak pada arah horizontal (Kamarudin et al.1994). Untuk memecahkan masalah tersebut perlu dilakukan suatu analisis penentuan sudut putaran pada sudut-sudut tertentu dalam mengatasi sebaran suhu rak di dalam ruang pengering. Sebaran suhu udara di dalam ruang pengering sangat penting untuk mengetahui pola sebaran aliran udara panas serta tingkat keragamannya. Pengukuran untuk melihat pola sebaran udara pengering secara bersamaan sulit untuk dilakukan karena keterbatasan alat ukur. Melalui pemodelan simulasi CFD (Computational Fluid Dyanamics) pada alat pengering Efek Rumah Kaca (ERK) tipe rak berputar dengan sistem hybrid (surya), sebaran suhu udara di setiap posisi rak dalam ruang pengering atau produk yang dikeringkan dapat diduga. Metode CFD menggunakan analisis numeric yaitu kontrol volume sebagai elemen dari intergrasi persamaan yang terdiri dari persamaan keseimbangan massa, momentum dan energi (Versteeg dan Malalasekera 1995). CFD pada penelitian ini menggunakan software Gambit dan software Fluent yang mampu melakukan simulasi sebaran udara panas dan kecepatan aliran udara kondisi nyata ke dalam virtual model atau prototipe dari sebuah sistem. Perumusan Masalah Suhu rak pengering di dalam pengering ERK tipe rak berputar sangat berfluktuasi karena perubahan iradiasi surya yang sangat dipengaruhi oleh

2 waktu, lokasi dan musim. Hal ini menyebabkan perbedaan suhu rak antara rak bagian atas berbeda dengan rak bagian bawah. Sehingga diperlukan penentuan putaran rak dengan sudut-sudut tertentu agar suhu pada setiap rak memiliki suhu yang seragam. Sebaran suhu pada pengering perlu diketahui untuk melihat keragaman suhu yang terjadi, hal ini sulit diketahui karena terbatasnya peralatan ukur. Oleh karena itu, melalui metode CFD, sebaran suhu ini dapat diprediksi. Tujuan Penelitian Tujuan penelitian ini adalah untuk menentukan putaran rak yang tepat, membuat pemodelan simulasi CFD untuk menggambarkan sebaran suhu pada rak pengering. Serta melakukan validasi antara hasil pengukuran dan perhitungan dengan CFD. Manfaat Penelitian Hasil dari penelitian ini diharapkan mampu memberikan rekomendasi pada operasi pengering dalam memperoleh hasil pengeringan yang merata dan efisien pada ERK tipe rak berputar. Hasil penelitian ini bermanfaat bagi pemakai alat pengering ERK tipe rak berputar, terutama bagi para petani maupun industri. Ruang Lingkup Penelitian Penelitian ini meliputi metode penentuan putaran rak pada pengering ERK selama proses pengeringan untuk memperoleh hasil pengeringan yang merata dan efisien. Enam skenario percobaan putaran rak dilakukan untuk menentukan putaran rak terbaik yang dinyatakan oleh keragaman suhu rak terkecil. Penentuan sebaran suhu pada setiap lokasi di dalam ruang pengering terutama pada bagian rak perlu diketahui. Metode CFD dengan software Gambit dan software Fluent dapat digunakan untuk menduga sebaran udara panas dan kecepatan aliran udara kondisi nyata ke dalam virtual model atau prototipe dari sebuah sistem. Validasi perhitungan CFD dinyatakan valid jika error yang dihasilkan lebih kecil dari 10% ketika dibandingkan dengan hasil pengukuran.

TINJAUAN PUSTAKA Pengeringan Hall (1957) dan Brooker (1974) menyatakan pengeringan merupakan proses penurunan kadar air bahan sampai mencapai kadar air tertentu sehingga dapat memperlambat laju kerusakan produk akibat aktivitas biologi dan kimia. Ada tiga hal yang mempengaruhi proses pengeringan yaitu kecepatan udara, suhu udara, dan kelembaban udara. Dasar

3 pengeringan yaitu kecepatan udara, suhu udara, dan kelembaban udara. Dasar proses pengeringan adalah terjadinya proses penguapan air bahan ke udara karena perbedaan kandungan uap air antar udara dengan bahan yang dikeringkan. Selama proses pengeringan terjadi dua proses yaitu proses pindah panas dan pindah massa air yang terjadi secara simultan. Panas dibutuhkan untuk menguapkan air bahan yang akan dikeringkan (Fellow 2001). Penguapan terjadi karena suhu bahan lebih rendah dari pada suhu udara di sekelilingnya. Proses pindah panas diperlukan untuk memindahkan massa uap air dari permukaan ke udara. Pindah panas terjadi karena tekanan uap air di dalam bahan lebih tinggi dari pada udara. Mekanise pengeringan diterangkan melalui teori tekanan uap, air yang diuapkan terdiri dari air bebas dan air terikat. Air bebas berada di permukaan telah habis, maka terjadi migrasi air karena perbedaan tekanan pada bagian dalam dan bagian luar (Henderson dan Perry 1976). Pengering Efek Rumah Kaca Hybrid Tipe Rak Berputar Pengering efek rumah kaca di Institut Pertanian Bogor pertama kali dikembangkan oleh Abdullah 1999 dari Departemen Teknik Mesin dan Biosistem. Dalam perkembangannya, pengering ini didisain dengan berbagai tipe sesuai dengan kebutuhan produk yang dikeringkan, diantaranya adalah pengering tipe rak untuk panili (Mursalim 1995), cengkeh (Wulandani 2005), kakao (Nelwan 1997), tipe bak untuk kakao (Manalu 1999), kopi (Wulandani 1997). Dalam rangka untuk menyeragamkan aliran udara di dalam ruang pengering terutama pada tipe rak, (Wulandani 2009) merancang pengering ERK tipe rak berputar secara vertikal. Pengering rumah kaca ini terdiri dari tiga bagian utama, yaitu bangunan rumah kaca, silinder dengan rak pengering dan bagian pemanas tambahan. Rumah kaca tersebut berukuran panjang x lebar x tinggi (2.15 m x 1.75 m x 1.9 m). Tungku untuk tempat pembakaran biomassa dan tanki untuk memanaskan air sebagai supply panas. Pompa air digunakan untuk sirkulasi sedangkan radiator untuk pembangkit panasnya. Bangunan rumah kaca berfungsi sebagai pengumpul panas. Pengering Efek Rumah Kaca (ERK) hybrid tipe rak berputar dapat dilihat pada Gambar 1.

Gambar 1 Alat pengering Efek Rumah Kaca (ERK) hybrid tipe rak berputar

4 Computational Fluid Dynamics (CFD) Computational Fluid Dynamics (CFD) adalah ilmu yang mempelajari cara memperdiksi aliran fluida, perpindahan panas, reaksi kimia, dan fenomena lainya dengan menyelesaikan persamaan-persamaan matematika atau model matematika (Tuakia 2008). Menurut (Versteeg dan Malalasekera 1995) CFD adalah sistem ananlisis yang meliputi aliran fluida, perpindahan panas dan fenomena seperti reakasi kimia yang berdasarkan simulasi komputer. Pada saat sekarang aplikasi CFD sudah banyak diterapkan pada disain mesin, ruang pembakaran, gas turbin dan tungku pembakaran. Ditinjau dari istilah CFD bisa berarti suatu teknologi komputasi yang memungkinkan untuk mempelajari dinamika dari benda-benda atau zat-sat yang mengalir. Perangkat lunak (software) CFD mampu untuk melakukan simulasikan aliran fluida, perpindahan panas, perpindahan massa, bendabenda bergerak, aliran multifasa, reaksi kimia, interaksi fluida dengan struktur dan sistem akustik hanya dengan pemodelan di komputer. Software CFD ini dapat membuat virtual prototipe dari sebuah sistem atau alat yang ingin dianalisis dengan menerapkan kondisi nyata di lapang. Software CFD akan memberikan data-data, gambar-gambar atau kurva-kurva yang menunjukkan prediksi dari performansi keandalan sistem desain. Hasil analisis CFD dapat berupa prediksi kualitatif maupun prediksi kuantitatif tergantung dari persoalan dan data input. Berikut adalah tiga keuntungan penggunaan CFD. 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) 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 Seorang engineer dalam merancang mesin memiliki berbagai alasan untuk mendapatkan pengembangan persamaan-persamaan, namun masalah utamanya dalam mendesain sistem termal adalah: (1) menyediakan proses simulasi termal (2) mengembangkan peryataan secara matematis untuk optimasi. Simulasi secara nyata dan problem-problem optimasi harus dilaksanakan dengan menggunakan komputer. Berikut adalah tahapan proses simulasi CFD. 1. Prapemprosesan (Preprocessing) Prapemprosesan adalah membuat model dalam paket CAD (Computer Aided Design), membuat grid yang sesuai dan membuat kondisi batas dan sifat-sifat fluidanya.

5 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. Pascapemprosesan (Postprocessing) Tahap pascapemprosesan adalah tahap penyajian hasil simulasi CFD dengan visualisasi warna untuk memudahkan dalam menganalisis. Berikut adalah tampilan hasil CFD. a. Hasil geometri dan grid yang terbentuk. b. Plot berdasarkan vektor. c. Plot berdasarkan kontur. d. Plot berdasarkan permukaan (2D atau 3D). Semua program pendekatan CFD dilakukan melalui prosedur sebagai berikut: 1. Sofware Gambit melakukan: a. Pembuatan geometri model b. Bidang atau volume yang diisi oleh fluida dibagi menjadi sel-sel kecil (meshing) 2. Sofware Fluent melakukan: a. Pendefinisian model fisik, misalnya : persamaan-persamaan gerak, entalpi, konversi spesies untuk zat-zat yang didefinisikan. b. Pendefinisian kondisi-kondisi batas, termasuk di dalamnya sifat-sifat dan prilaku dari batas-batas model. Untuk kasus transient, kondisi awal juga didefinisikan. c. Persamaan-persamaan matematika yang membangun CFD diselesaikan secara iterasi, biasa dalam kondisi tunak (steady state) atau transient. d. Analisis dan visualisasi dari solusi CFD.

METODOLOGI PENELITIAN Waktu dan Tempat Penelitian Kegiatan penelitian dilaksanakan mulai Februari 2013 sampai Agustus 2013 dan berlokasi di Laboratorium Lapang Siswadhi Soepardjo, Leuwikopo dan Laboratorium Energi, Departemen Teknik Mesin dan Biosistem, Fakultas Teknologi Pertanian, Institut Pertanian Bogor. Peralatan Alat pengering yang digunakan pada penelitian ini adalah alat pengering Efek Rumah Kaca (ERK) Hybrid rancangan Wulandani, et al. 2009 (Gambar 1). Deskripsi ERK tipe rak berputar adalah sebagai berikut: 1. Bangunan pengering ERK berdimensi 2.15 m x 1.75 m x 1.9 m merupakan bangunan pengering ERK hybrid yang di dalamnya terdapat

6

2. 3. 4.

5.

1. 2. 3. 4. 5. 6. 7.

rak-rak yang digantungkan pada silinder berputar sehingga dengan adanya putaran as menyebabkan rak-rak pengering dapat bergerak mengikuti rotasi gerakan silindernya. Pada bangunan tersebut di dalamnya terdapat 3 bagian silinder, dalam setiap ruang pengering terdapat 8 gantungan rak pengering dengan posisi rak atas dan bawah sehingga rak-rak pengering berjumlah 48 rak. Rangka bangunan, rangka silinder dan rangka rak terbuat dari besi siku dan besi hollow bercat hitam. Dinding terbuat dari plastik polikarbonat. Rak berupa aluminium berlubang, diberi bingkai dari kayu dengan dimensi 0.3m x 0.6m x 0.03 m, jumlah rak di setiap ruang adalah 16 buah, sehingga total rak adalah 48 buah. Penukar panas, berupa radiator dengan sistem silang kompak, berdimensi 0.3 m x 0.35 m x 0.02 m. Fluida kerja yang dipakai pada penukar panas adalah air. Kipas, sebuah kipas dengan daya 80 W terdapat di bawah radiator, berdiameter 0.3 m. Kipas outlet berupa exhaust fan berdiameter 0.12 m dengan daya 30 W. Pada pengering ini digunakan 3 buah kipas radiator, 3 buah kipas outlet dan 3 buah kipas pengaduk. Tungku untuk tempat pembakaran biomassa dan tangki untuk memanaskan air agar dapat digunakan untuk sirkulasi sedangkan radiator untuk pembangkit panasnya. Peralatan ukur yang digunakan dalam penelitian meliputi: Thermocouple tipe CC. Termometer bola basah dan bola kering. Chino Recorder Yokogawa tipe 3058 dengan ketelitian 0.1oC Pyranometer model EKO tipe MS-401 Multimeter digital model YEW tipe 2506 A dengan ketelitian 0.01 mV Hot wire Anemometer model Lutron tipe AM-4204HA dengan ketelitian 0.01 m/dt. Perangkat lunak yang di pakai yaitu software Gambit 2.4.6 diperlukan untuk proses pembuatan gambar dan meshing, ANSYS 12.1: Fluent 6.3 untuk melakukan proses simulasi dan analisis menggunakan Ms.Office dan Ms. Excel 2007. Prosedur Penelitian

Pada penelitian ini dilakukan dengan beberapa tahapan, dapat dilihat pada diagram alir penelitian 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 titik pengukuran lantai, lantai, dinding depan (pintu), dinding kanan, dinding kiri, dinding belakang, atap, ruang pengering. Dan melakukan pengukuran dimensi alat serta mengindetifikasi bahan konstruksi alat pengering.

7 Percobaan Pengeringan Percobaan pengeringan dilakukan sebanyak 6 skenario. 1 Skenario 1 : Rak dalam keadaan diam (tidak diputar). 2 Skenario 2 : Rak diputar 45o selama 1 jam sekali. 3 Skenario 3 : Rak diputar 45o selama 2 jam sekali. 4 Skenario 4 : Rak diputar 90o selama 1 jam sekali. 5 Skenario 5 : Rak diputar 135o selama 1 jam sekali. 6 Skenario 6 : Rak diputar 180o selama 1 jam sekali. Penamaan rak diberi label A, B, C, D, E, F, G dan H. Posisi penamaan rak sama untuk semua skenario. Penamaan posisi rak dalam posisi yang tetap walaupun rak dilakukan pemutaran sudut sesuai dengan skenario. Pada saat pengambilan data rak dalam posisi diam (tidak diputar) dan posisi rak akan berpindah dengan penamaan yang sudah diberikan. Percobaan dilakukan tanpa beban selama 1 hari dari pukul 08.00 WIB sampai 16.00 WIB dengan interval pengukuran 0.5 jam. Pemutaran rak dilakukan dengan interval 45o hal ini dikarenakan pengunci poros pemutar rak didesain sebesar 45o. Perlakuan putaran rak dilakukan sejauh 180o hal ini dikarenakan bentuk rotasi pergerakan silinder berbentuk lingkaran sehingga dilakukan pemutaran rak sejauh 180o untuk mewakili pemutaran rak pada bagian atas, tengah, dan bawah. Pemutaran rak dilakukan selama 1 jam sekali agar ergonomis saat pemutaran rak, dan lama pengambilan data saat pengukuran. Parameter yang Diukur pada Percobaan Pengeringan Parameter yang diukur adalah sebagai berikut: 1. Suhu yaitu meliputi suhu udara pada rak atas, rak tengah, rak bawah, lantai, inlet, outlet, suhu ruang pengering, suhu dinding, suhu atap, suhu lantai pengering, dan suhu lingkungan. 2. Kecepatan meliputi kecepatan udara pada inlet, outlet, depan kipas kiri, kipas tengah dan kipas kanan, udara pada rak atas, rak tengah, rak bawah, RH ruang pengering, dan RH udara lingkungan. 3. Radiasi surya. Titik pengukuran suhu rak pada rak atas 3 titik pengukuran dan rak bawah 2 titik pengukuran. Pengukuran kondisi batas pada lantai, atap, dinding kanan, dinding kiri, pintu (dinding depan), dan suhu ruang pengering. Posisi titik pengukuran dapat dilihat pada Gambar 2. Pemodelan CFD Pengukuran suhu-suhu tersebut digunakan pula untuk menvalidasikan simulasi CFD.Simulasi CFD dilakukan pada kondisi tanpa menggunakan beban. Pemodelan CFD menggunakan asumsi: 1. Udara tidak termampatkan (incompressible), ρ konstan. 2. Bilangan Prandtl udara konstan (panas jenis, konduktivitas, dan vikositas udara konstan). 3. Udara bergerak dalam keadaan steady. 4. Udara lingkungan dianggap konstan selama simulasi. 5. Kecepatan angin dianggap konstan.

8 6. Aliran udara dianggap laminar. Aliran laminar jika Re<2000 (Holman 1997). Hasil pengukuran di atas selanjutnya menjadi input untuk menyusun pemodelan CFD diantaranya sebagai input kondisi batas. Nilai dan perhitungan aliran laminer pada ruang pengering ditunjukkan pada Lampiran 1. Perhitungan Parameter Kondisi Batas simulasi CFD Kondisi batas untuk pemodelan CFD meliputi beberapa parameter. Perhitungan parameter kondisi batas pada dinding menggunakan persamaan sebagai berikut : a. Perhitungan aliran laminer. (Holman1997). Re = ..............................................................................................(1) b.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) : Gr = ............................................................................................ (2) Ra = Gr x Pr...................................................................................... (3) Nu =

.............................................. (4)

h= .............................................................................................. (5) Iradiasi surya didefinisikan sebagai fluks dari atap dan dinding bangunan pengering. Fluks panas dapat dihitung melalui persamaan sebagai berikut: = (I x τα) – h ( Td – T∞ )................................................................. (6) c. Perhitungan parameter kondisi batas pada rak. Rak berupa plat berlubang dianggap sebagai porous jump. Nilai α dan C2 dinyatakan dalam persamaan yang terdapat pada Wulandani (2009). a. Permeabilitas permukaan (α) α=

....................................................................................(7)

b. Koefisien porous jump (C2) C2 = ......................................................................................(8) Rincian perhitungan parameter kondisi batas pada kondisi tanpa beban dapat dilihat pada Lampiran 1. Pembuatan Model Geometri dan Pendefinisian Kondisi Batas Simulasi CFD Alat pengering ERK yang akan disimulasikan terlebih dahulu dibuatkan model geometri 3D dengan menggunakan perangkat lunak Gambit 2.4.6 yakni atap, lantai, dinding, atap, inlet dan outlet beserta rakrak pengering. Setelah model geometri dibuat, dilakukan pendefinisian bidang batas setiap bagian geometri, seperti atap, dinding dan lantai didefinisikan sebagai wall, inlet sebagai velocity inlet, outlet sebagai outflow

9 dan rak-rak pengering sebagai porous jump. Setelah pendefinisian selesai dilakukan proses meshing. Simulasi CFD Udara Panas Program fluent 6.3 melakukan beberapa proses sebagai berikut: 1. Mendefinisikan: a. Model, dimana didalamnya ditentukan solver (2D atau 3D), pemakaian energi, viskositas model (laminer atau turbulen) b. Menentukan jenis fluida dan material penyusun bagunan pengering yang digunakan serta sifat termofisiknya. c. Menentukan kondisi operasi (Operating conditions) d. Memasukkan nilai-nilai kondisi batas (Boundary conditions) terhadap domain yang sudah dibuat dengan program Gambit. 2. Melakukan proses inisialisasi. 3. Melakukan proses iterasi. 4. Melihat penampilan hasil simulasi dalam bentuk grid, kontur (suhu kecepatan angin),vektor (suhu, kecepatan). 5. Mendapatkan informasi data yang terkait hasil simulasi untuk keperluan validasi Plot (XY Plot, Residual) Analisis Sebaran Suhu Hasil Simulasi CFD Dari informasi sebaran suhu dengan beberapa skenario yang dilakukan maka dapat dilihat dari pola sebaran suhu udara panas yang dihasilkan pada penampilan hasil simulasi CFD. Validasi Model Simulasi CFD Validasi model dilakukan dengan membandingkan suhu rak hasil simulasi dengan suhu rak hasil pengukuran di lapang. Keakuratan hasil simulasi dengan hasil pengukuran dinyatakan dalam persentase error. Error dinyatakan dengan persamaan: Error=

......................... (9)

Menurut Karmeli (2011) dalam Nuryawati (2011), distribusi suhu hasil simulasi juga divalidasi dengan menghitung koefisien keragaman dan koefisien variasi. Apabila Nilai koefisien variasi (CV) kurang dari 15%, maka dikatakan bahwa hasil simulasi telah baik. Persamaan dinyatakan sebagai berikut: CV =

............................................................................................ (10)

Dimana CV adalah koefisien variasi (coefficient of variation), σ adalah standar deviasi, μ adalah rata-rata suhu hasil simulasi (oC).

10 Posisi Titik Pengukuran ERK Tipe Rak Berputar

H

A

G Arah kanan ruang pengering

B

Arah kiri ruang pengering F

C D

E

Gambar 2 Titik Pengukuran ERK Tipe Rak Berputar Keterangan: Pada saat pengukuran di lapangan rak diberi tanda yakni A, B, C, D, E ,F, G dan H. Pemberian tanda posisi rak sesuai tingkatan posisi rak pada bagian atas, tengah dan bawah. Rak A pada posisi rak teratas sebelah kanan dan rak H pada posisi rak teratas sebelah kiri, rak B pada posisi dibawah rak A sebelah kanan dan rak G pada posisi bawah rak H sebelah kiri, rak F pada posisi bagian tengah sebelah kiri, rak C pada posisi tengah sebelah kanan, rak D pada posisi rak terbawah sebelah kanan dan rak E pada posisi rak terbawah sebelah kiri. Pada pada saat pemutaran posisi rak akan berubah sesuai skenario yang dilakukan dengan penamaan yang sama. Terdapat tiga ruang rak pengering, dalam satu ruang pengering terdapat 8 jumlah gantungan rak pengering dalam posisi rak atas dan bawah sehingga terdapat 48 rak. Titik pengukuran suhu rak pada rak atas 3 titik pengukuran dan rak bawah 2 titik pengukuran. Pengukuran kondisi batas pada lantai, atap, dinding kanan, dinding kiri, pintu (dinding depan), dan suhu ruang pengering.

11

Mulai

Simulasi CFD

Experimen pengering ERK

Tidak

Valid (Error suhu < 10%)

Ya Skenario terbaik

Selesai

Gambar 3 Diagram alir prosedur penelitian

HASIL DAN PEMBAHASAN Suhu Rak Pengering Hasil Pengukuran Skenario 1 Pada skenario 1 rak tidak dilakukan pemutaran (diam), sehingga posisi rak dalam keadaan tetap. Tabel 1 memperlihatkan tingkat keragaman suhu rak selama pengukuran yang dinyatakan dengan nilai standar deviasi. Nilai standar deviasi rata-rata skenario 1 adalah 0.52oC pukul 08.00 WIB.

12 Gambar 4 memperlihatkan pengukuran suhu rak pada pukul 08.00 WIB pada rak A, adalah suhu rak tertinggi. Pada pukul 12.00 WIB suhu rak tertinggi masih terjadi pada rak A karena posisi rak A berada di bagian paling atas, karena tidak ada perubahan posisi. Nilai standar deviasi pengukuran pukul 12.00 WIB menunjukkan keragaman tertinggi, yaitu 2.48oC. Rata-rata nilai standar deviasi selama 1 hari pengukuran suhu pada skenario 1 sebesar 1.70oC. Nilai iradiasi matahari pada pukul 08.00WIB adalah 285.71W/m2dan pada pukul 12.00 WIB adalah 314.285 W/m2. Pada pukul 12.00 WIB, nilai iradiasi matahari yang dihasilkan rendah, hal ini dikarenakan pada pukul 12.00 WIB mendung (matahari tertutup oleh awan). Adapun suhu yang terukur pengukuran pukul 12.00 WIB suhu yang dihasilkan besar, hal ini dikarenakan pada saat pengukuran pukul 12.00 WIB dilakukan tambahan pemanas yakni biomassa pada saat kondisi cuaca mendung, tambahan pemanas dilakukan pada saat pengukuran pukul 11.00 WIB -13.00 WIB. Berikut adalah profil sebaran suhu rata-rata setiap rak pada pengukuran tiga titik bagian rak atas dan dua titik bagian rak bawah. Suhu maksimum 44.6oC, suhu minimum 42.3oC dan suhu rata-rata 43.09oC pengukuran pukul 08.00 WIB. Pada suhu maksimum 69.8oC, suhu minimum 62oC dan suhu rata-rata 67.07oC pengukuran pukul 12.00 WIB. Tabel 1 Tingkat keragaman suhu rak Skenario 1 Waktu (Pukul)

Suhu rata-rata (oC)

41.0 52.4 57.5 57.4 63.1 64.2 58.3 54.8 48.4 40.8

Ratarata (oC) 41.4 49.6 56.0 59.3 67.0 67.4 59.3 55.6 49.2 43.2

0.52 2.42 1.11 4.14 2.49 1.92 0.86 1.68 0.72 1.73

54.9

56.3

1.78

Rak

A

B

C

D

E

F

G

H

08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 17.00 Ratarata (oC)

42.3 50.9 56.6 59.6 69.3 70.3 61.1 56.7 50.3 46.2

41.6 49.5 56.1 63.8 68.8 69.1 59.9 55.9 50.0 45.2

42.7 49.5 56.2 61.6 69.2 68.3 59.2 55.4 49.4 43.2

41.4 49.1 56.0 59.3 69.2 68.2 59.3 55.2 49.4 43.2

41.9 49.4 54.5 58.7 67.8 67.6 59.4 55.5 49.3 42.6

40.2 48.6 55.0 57.4 64.8 66.1 59.1 55.6 48.7 42.0

40.5 47.8 56.4 57.0 64.1 65.8 58.7 55.4 48.2 42.0

58.3

57.6

57.1

56.6

56.1

55.5

55.2

SD (oC)

13

Waktu (Pukul)

Gambar 4 Profil suhu rak Skenario 1 Hasil Pengukuran Skenario 2 Skenario 2 dilakukan dengan memutar rak sebesar 45o setiap jam sekali. Tabel 2 memperlihatkan tingkat keragaman suhu rak selama pengukuran dinyatakan dengan nilai standar deviasi. Nilai standar deviasi pengukuran pukul 11.00 WIB lebih tnggi dibanding pengukuran pukul 12.00 WIB hal ini dikarenakan nilai iradisi pengukuran pada pukul 10.00 WIB sangat rendah, kemudian meningkat pada saat pengukuran pukul 11.00 WIB sehingga suhu yang diterima pada saat pengukuran pukul 11.00 WIB signifikan berubah dan menyebabkan rak tidak menerima suhu yang merata. Nilai radiasi matahari pengukuran pukul 10.00 WIB sebesar 590.00 W/m2 dan nilai iradiasi pengukuran pukul 11.00 WIB sebesar 782.00 W/m2 dan nilai standar deviasi pengukuran pukul 12.00 WIB suhu yang diterima rak dengan peningkatan nilai iradiasi tidak terlalu berubah mengakibatkan suhu rak lebih seragam. Pengukuran pukul 11.00 WIB suhu rak tertinggi pada rak G dan H dan pukul 12.00 WIB pada rak D dan E. Pemutaran rak memiliki kelebihan yaitu rak dapat menempati posisi panas matahari yang bergantian di dalam ruang pengering, pemutaran rak mengakibatkan posisi rak tengah, bawah dan atas menempati posisi bergantian. Nilai standar deviasi dipengaruhi jumlah iradiasi pada saat pengukuran, nilai iradiasi surya ratarata pukul 08.00 WIB adalah 571.43W/m2 dan nilai iradiasi pada pukul 12.00 WIB adalah 890.42W/m2. Gambar 5 menunjukkan profil suhu ratarata pada setiap rak, pengukuran pukul 08.00 suhu maksimum 40.7oC, suhu minimum 35.12oC dan suhu rata-rata 38.14oC dan pengukuran pukul 12.00 WIB suhu maksimum 66.7oC, suhu minimum 65oC dan suhu rata-rata 62.75oC.

14 Tabel 2 Tingkat keragaman suhu rak Skenario 2 Waktu (Pukul)

Suhu rata-rata (oC)

Rak

A

B

C

D

E

F

G

H

08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 Ratarata (oC)

40.6 46.5 46.3 63.9 66.7 72.2 64.2 52.3 38.7

40.6 45.6 45.4 62.5 66.4 70.0 61.8 50.4 39.6

40.3 44.4 43.4 60.5 69.1 65.4 62.8 51.9 42.2

39.7 43.4 42.3 58.4 72.6 68.6 65.4 52.6 44.0

38.6 42.6 43.6 60.1 72.5 72.0 67.3 53.3 45.3

38.6 43.5 45.4 59.9 66.1 73.2 68.1 56.7 44.1

39.7 44.5 46.6 65.2 66.9 74.4 67.6 54.4 42.6

40.6 46.5 47.4 65.6 67.3 73.4 66.5 54.4 41.7

Ratarata (oC) 39.8 44.6 45.0 62.0 68.4 71.1 65.4 53.2 42.2

55.72

54.6

53.5

54.2

55.2

56.1

56.9

56.4

55.3

SD (oC) 1.75 1.39 1.71 2.89 1.47 2.88 2.24 1.80 2.27 2.04

Waktu (Pukul)

Gambar 5 Profil suhu rak Skenario 2 Hasil Pengukuran Skenario 3 Pengukuran skenario 3 dengan melakukan pemutaran rak sebesar 45o melakukan pemutaran rak dua jam sekali. Hal yang membedakan pada skenario 2 yakni lama pemutaran rak. Tabel 3 memperlihatkan tingkat keragaman dinyatakan dengan standar deviasi. Nilai suhu rak pada skenario 3 kurang seragam dari pada pengukuran skenario 2. Nilai rata-rata standar deviasi suhu rak skenario 3 pengukuran pukul 08.00 WIB sebesar 1.75oC dan pukul 12.00 WIB sebesar 1.60oC, nilai rata-rata standar deviasi pengukuran skenario 3 adalah 2.12oC. Nilai standar deviasi yang dihasilkan lebih besar dibanding hasil pengukuran standar deviasi skenario 2, hal ini dikarenakan perbedaan lama waktu pemutaran yakni pemutaran rak pada skenario 2 rak diputar 1 jam sekali sedangkan pemutaran rak skenario 3 dilakukan selama 2 jam sekali. Suhu maksimum pada rak A dan B pukul 08.00 WIB dan pukul 12.00 WIB posisi rak G dan H terletak pada posisi

15 atas pengering, pengukuran pukul 13.00 WIB posisi rak menerima suhu maksimum pada G dan H, posisi rak pada pengukuran pukul 13.00 WIB dan pengukuran pukul 12.00 WIB menempati posisi suhu rak tertinggi hal ini dikarenakan pemutaran sudut sebesar 45o pada pemutaran rak selama 2 pukul sekali. Pengukuran pada skenario 2 posisi rak pengukuran pada pukul 12.00 WIB dan pengukuran pukul 13.00 WIB, suhu rak maskimum pengukuran pukul 12.00 WIB pada rak D dan E sedangkan pengukuran pukul 13.00 WIB pada rak G dan H. Salah satu kelebihan pemutaran rak adalah posisi rak menempati tempat yang berpindah-pindah sehingga masing-masing rak menerima panas matahari bergantian. Lama pemutaran rak menentukan besarnya suhu standar deviasi yang dihasilkan. Gambar 6 memperlihatkan profil suhu rata-rata rak, suhu rata-rata rak dipengaruhi iradiasi surya rata-rata saat dilakukan pengukuran. Pengukuran pukul 08.00 WIB iradiasi sebesar 471.43W/m2, nilai iradiasi pada pukul 12.00 WIB adalah 760.42W/m2. Pengukuran pukul 08.00 WIB suhu maksimum 40.7oC, suhu minimum 35.2oC dan suhu rata-rata 35.14 oC dan pengukuran pukul 12.00 WIB suhu maksimum 67.52oC, suhu minimum 65oC dan suhu rata-rata 65oC. Tabel 3 Tingkat keragaman suhu rak Skenario 3 Waktu (Pukul) Rak 08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 Ratarata (oC)

Suhu rata-rata (oC)

40.5 46.5 46.2 63.9 67.6 69.0 64.2 52.3 38.7

39.9 45.5 45.4 62.5 66.5 70.0 61.7 50.4 39.6

39.4 44.4 43.4 60.5 66.3 65.4 62.8 51.8 42.2

37.6 43.4 42.3 58.4 66.2 68.5 65.1 52.6 44.1

37.0 42.6 43.6 60.1 67.4 70.0 67.2 53.3 45.3

36.4 43.5 45.3 59.8 68.7 73.2 68.1 56.7 44.1

38.7 44.5 46.5 65.1 69.8 74.4 67.5 54.3 42.6

35.6 46.4 47.6 65.5 69.6 73.4 66.5 53.2 40.4

Ratarata (oC) 38.1 44.6 45.0 62.0 67.8 70.5 65.4 53.1 42.1

57.2

56.0

54.5

55.1

56.0

57.0

58.2

57.8

56.52

A

B

C

D

E

F

G

H

SD ( (oC) 1.75 1.39 1.71 2.89 1.60 2.88 2.24 1.80 2.27 2.12

16

Waktu (Pukul)

Gambar 6 Profil suhu rak Skenario 3 Hasil Pengukuran Skenario 4 Skenario 4 rak diputar selama 1 jam dengan sudut putaran 90o, tabel 5 memperlihatkan nilai standar deviasi yang dihasilkan pengukuran pukul 08.00 WIB 0.68oC dan pengukuran pukul 12.00 WIB adalah 2.95oC. Nilai standar deviasi pengukuran pukul 08.00 WIB pada skenario 4 lebih rendah dibanding pengukuran skenario 1, 2 dan 3. Hal ini dikarenakan nilai iradiasi saat pengukuran lebih rendah dibanding nilai iradiasi skenario 1, 2 dan 3. Rata-rata nilai standar deviasi pengukuran skenario pukul 09.00 WIB-16.00 WIB sebesar 2.19oC. Nilai suhu maksimum pada rak H dan A pengukuran pukul 12.00 WIB, posisi rak H dan B berada pada posisi bagian paling atas rak pengering dan pengukuran pukul 08.00 WIB pada rak A dan B berada pada posisi paling atas rak pengering. Kesergaman suhu rata-rata rak terlihat pada Gambar 7 yakni suhu rak kurang merata. Iradiasi surya ratarata pukul 08.00 WIB 430.01W/m2 dan nilai iradiasi pengukuran pukul 12.00 WIB adalah 970.42W/m2. Pengukuran pukul 08.00 WIB suhu maksimum 41.3oC, suhu minimum 37oC suhu rata-rata 39.9oC dan pengukuran pukul 12.00 WIB suhu maksimum 70.7oC, suhu minimum 66.2oC , suhu rata-rata 69.41oC.

17 Tabel 4 Tingkat keragaman suhu rak Skenario 4 Waktu (Pukul)

Suhu rata-rata (oC) SD

39.2 58.5 62.5 73.3 70.7 77.6 62.9 54.3 54.2

Ratarata (oC) 39.9 58.8 62.7 72.2 68.7 73.5 63.9 56.4 54.1

62.6

62.0

2.1

Rak

A

B

C

D

E

F

G

H

08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 Ratarata (oC)

41.25 57.4 57.7 69.5 70.6 73.4 62.3 57.6 54.5

40.4 56.5 61.5 72.2 70.2 73.1 62.8 57.8 57.5

40.1 55.1 63.1 74.4 69.8 70.1 63.7 58.5 54.4

39.5 56.1 63.7 72.6 68.5 72.5 65.5 58.4 53.7

39.5 58.6 64.7 73.3 67.4 73.2 65.3 57.5 51.3

39.6 58.6 64.7 73.3 66.3 70.3 65.5 54.2 53.3

39.4 69.5 63.7 68.9 66.2 77.9 63.5 53.8 53.8

62.2

62.0

61.9

61.6

61.9

62.3

61.6

(oC) 0.6 1.4 3.1 1.9 2.9 2.7 1.2 1.9 2.0

Waktu (Pukul)

Gambar 7 Profil suhu rak Skenario 4 Hasil Pengukuran Skenario 5 Pengukuran pada skenario 5 rak diputar dengan sudut 135o selama 1 jam sekali. Tabel 5 memperlihatkan nilai standar deviasi dihasilkan pukul 08.00 WIB sebesar 1.64oC dan pengukuran pukul 12.00 WIB 1.85oC. Ratarata nilai standar deviasi pengukuran skenario 5 sebesar 1.67oC nilai standar deviasi pengukuran skenario 5 merupakan nilai standar deviasi terkecil dibanding skenario yang lain (skenario 1 sampai skenario 6). Hal ini dikarenakan pemutaran rak sejauh 135o rak berpindah jauh sehingga rak pada posisi tengah dan bawah bergeser ke atas bagian rak pengering, pemutaran rak sejauh 135o mempengaruhi suhu rak pada bagian tengah dan bawah berganti posisi menempati posisi rak atas dalam ruang pengering. Pada Gambar 8 memperlihatkan suhu rak maksimum pada rak A dan E pengukuran pukul 08.00 WIB dan rak E dan F pengukuran pukul 12.00 WIB. Masing-masing skenario menghasilkan posisi suhu rak maksimum

(

18 yang berbeda-beda sesuai dengan pergeseran sudut yang dilakukan saat pemutaran rak. Adapun kesulitan saat pengukuran yakni saat melakukan pemutaran rak, termokopel tergulung dan lepas dari posisi pengukuran. Pengukuran pukul 08.00 WIB iradiasi surya rata-rata 450.05W/m2 dan nilai iradiasi pengukuran pukul 12.00 WIB 970.42 W/m2. Pengukuran pukul 12.00 WIB suhu maksimum 49.75oC, suhu minimum 49.7oC, suhu rata-rata 49.75oC dan pengukuran pukul 08.00 WIB dan suhu maksimum 77.3oC, suhu minimum 73oC dan suhu rata-rata 74.68oC . Tabel 5 Tingkat keragaman suhu rak Skenario 5 Suhu rata-rata (oC)

Waktu (Pukul) Rak

A

B

C

D

E

F

G

H

08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 Ratarata (oC)

49.7 49.5 56.4 65.7 70.5 75.1 68.5 58.7 47.5

49.5 47.1 55.0 66.9 73.6 76.0 69.1 57.6 45.1

48.8 48.2 58.4 65.5 70.1 75.6 71.4 57.3 47.4

47.4 49.7 58.2 61.9 73.6 74.9 72.3 57.4 48.6

47.5 49.7 57.2 63.8 77.1 73.6 71.3 55.7 52.3

46.4 50.1 55.0 65.6 77.0 76.1 69.5 56.3 54.5

46.0 50.2 53.2 66.9 73.1 77.2 68.3 57.2 51.9

45.1 49.7 54.5 68.4 73.1 74.8 67.3 56.6 49.7

Ratarata (oC) 47.5 49.3 56.0 65.6 73.5 74.4 69.7 57.1 49.6

60.2

59.5

60.5

60.4

60.5

61.0

60.7

60.1

60.7

SD (oC) 1.6 1.0 1.8 1.9 1.8 1.0 1.7 0.9 3.0 1.6

Waktu (Pukul)

Gambar 8 Profil suhu rak Skenario 5 Hasil Pengukuran Skenario 6 Pengukuran pada skenario 6 rak diputar sebesar 180o diputar selama satu jam sekali. Hasil pengukuran pada skenario 6 suhu maksimum, pada rak A dan rak B pukul 12.00 WIB dan pukul 08.00 WIB. Keragaman rak

19 diyatakan pada Gambar 9 nilai suhu rak atas, tengah dan bawah kurang seragam. Suhu yang dihasilkan pada setiap pukul memiliki perbedaan yang besar. Rak bagian atas dan bawah akan menempati posisi yang bergantian, sedangkan rak bagian posisi tengah saat dilakukan pemutaran berada di posisi bagian tengah. Nilai standar deviasi pengukuran pukul 08.00 WIB 1.09oC dan pengukuran pukul 12.00 WIB 1.8oC. Nilai rata-rata standar deviasi rata-rata pengukuran pemutaran rak pukul 09.00 WIB - 16.00 WIB sebesar 2.09oC, nilai standar deviasi lebih tinggi dibanding pengukuran skenario 5. Kesulitan pada saat pengukuran yakni saat pemutaran termokopel tergulung dan putus dari rak hal ini dikarenakan rak diputar terlalu jauh. Iradiasi surya rata-rata pukul 08.00 WIB adalah 400.05W/m2 dan nilai iradiasi pada pukul 12.00 WIB adalah 856.00W/m2. pengukuran pukul 08.00 WIB suhu maksimum 42.1oC, suhu minimum 38.8oC, suhu rata-rata 40.1 oC, dan pengukuran pukul 12.00 WIB suhu maksimum 71.9oC, suhu minimum 68.7oC dan suhu rata-rata 70.56oC. Tabel 6 Tingkat keragaman suhu rak Skenario 6 Suhu rata-rata (oC)

Rak Waktu (Puku) 08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 Ratarata (oC)

A

B

42.1 50.1 57.6 52.2 71.6 60.8 61.4 58.7 47.1

41.3 53.2 57.6 55.2 71.6 61.6 61.7 59.5 46.7

55.4

56.5

C

39.2 52.0 57.4 56.6 71.5 63.6 61.6 59.0 45.5

Ratarata (oC) 40.1 53.9 57.1 56.1 70.5 64.1 58.6 60.1 45.8

(oC) 1.0 1.9 0.6 2.4 1.8 2.9 3.3 1.1 2.3

56.7

56.7

2.0

D

E

F

G

H

40.4 55.2 57.3 55.4 71.2 62.1 61.4 59.6 49.6

40.1 54.7 56.7 56.8 69.6 64.4 55.6 61.3 42.2

39.6 55.6 56.0 60.4 68.8 69.5 57.4 61.7 42.7

39.1 55.3 56.6 56.7 69.6 67.8 57.6 61.5 45.6

38.9 55.2 57.5 55.4 70.5 63.3 52.4 59.5 46.5

57.2

55.9

55.9

57.3

55.8

Waktu (Pukul)

Gambar 9 Profil suhu rak Skenario 6

SD

(

20 Penentuan Putaran Rak Terbaik Nilai standar deviasi rata-rata pengukuran masing-masing skenario terlihat pada Gambar 10. Nilai standar deviasi terendah pada pengukuran skenario 5 sebesar 1.67oC dengan pemutaran rak sejauh 135o selama 1 jam sekali, nilai standar deviasi terendah kedua yakni pengukuran pada skenario 1 nilai standar deviasi sebesar 1.70oC, pengukuran suhu rak pada skenario 2, 3 dan 4 menghasilkan nilai standar deviasi terbesar. Penentuan putaran terbaik dihasilkan dari nilai standar terendah saat pengukuran, hal ini menunjukkan tingkat keragaman suhu masing-masing rak pada posisi rak yang berbeda. Sehingga semakin besar standar deviasi semakin besar pula tidak seragamnya. Pemutaran rak memiliki kelebihan yakni posisi rak akan menempati letak yang bergantian sehingga masing-masing rak akan memperoleh panas yang lebih merata.

Gambar 10 Rata-rata standar deviasi masing-masing skenario Kelembaban Udara dan Iradiasi Surya pada Ruang Pengering Kelembaban udara ruang pengering dipengaruhi oleh suhu dan tekanan uap udara.Suhu yang meningkat di dalam ruang pengering menyebabkan tekanan uap udaranya juga meningkat sehingga menyebabkan tekanan uap yang tinggi bergerak menuju tekanan yang lebih rendah yaitu kearah outletnya. Kelembaban udara ruang pengering selama proses pengeringan berlangsung secara terus menerus akan mengurangi uap air yang ada di dalam ruang pengering, dengan demikian kelembaban udaranya dari semula tinggi menjadi rendah karena pergerakan udara dari tekanan tinggi menuju tekanan rendah. Kelembaban udara rata-rata ruang pengering yang diperoleh pada skenario 1adalah 52.82%, skenario 2 adalah 65.03%, skenario 3 adalah 50.78%, skenario 4 adalah 55.89%, skenario 5 adalah 45.37% dan skenario 6 adalah 48.52%. Faktor yang menyebabkan kelembaban udara ruang pengering berubah yaitu: iradiasi matahari yang berpengaruh terhadap perubahan suhu ruang pengering dan udara lingkungan, peningkatan suhu dalam ruang pengering oleh pengaruh penyebaran panas secara konveksi . Suhu dan kelembaban dalam ruang pengering mempengaruhi laju pengeringan. Pada penelitian ini suhu ruang pengering tanpa bahan memiliki suhu rata-rata 50oC-73oC, sedangkan suhu rata-rata pengering

21 mengunakan bahan produk pada pengeringan sawut ubi jalar menghasilkan suhu pengering rata-rata 40oC-60oC (Utari 2013). RH rata-rata pengering tanpa mengunakan bahan sebesar 53.06% sedangkan RH rata-rata pengering pada sawut ubi jalar sebesar 46.49%. Hal ini menunjukkan hubungan semakin besar suhu pada saat pengeringan maka kelembaban pada ruang pengering akan semakin rendah. Hubungan antara waktu pengamatan terhadap radiasi matahari mempengaruhi laju pengeringan. Nilai iradiasi matahari rata-rata masing-masing skenario 700W/M2. Nilai radiasi matahari, RH dan kecepatan angin terlampir pada Lampiran 4. Pemodelan Simulasi CFD Suhu Rak Pengering Bentuk Model Pengering dengan Gambit Bentuk domain 3D ruang pengering diperoleh dari proses pembentukkan geometri ruang pengering menggunakan software Gambit 2.2.30. Domain adalah betuk dasar bangunan pengering dengan kondisi batas pada dinding, atap, lantai, inlet outlet fan dan rak-rak pengering. Pembentukkan grid ruang pengering dengan nilai mesh 60 mm, karena ukuran worst mesh dihasilkan lebih kecil dari 0.85, nilai meshing yang sudah dibuat sudah baik (Tuakia F. 2008). Pembentukkan grid ini menggunakan Tet/ Hybrid yang berarti aliran fluida akan melalui jaringan mengikuti bentuk ruang pengeringnya dan jaringan tersebut dapat berbentuk segiempat, segilima, atau segienam sesuai bentuk ruang pengeringnya. Bentuk simulasi domain utama model pengering serta grid ruang pengering sama untuk 6 skenario. Grid domain memberikan jalur lintasan yang dilalui oleh rambatan atau aliran udara panas maupun aliran udara yang bergerak dari inlet menuju outletnya melalui pencampuran dan initial condition untuk membentuk distribusi dan kecepatan aliran udara. Hasil Simulasi Pukul 08.00 Simulasi CFD pengering ERK bertujuan untuk mengetahui pola dan sebaran suhu yang terjadi pada ruang pengering. Bentuk Simulasi dengan kondisi tanpa beban dilakukan dengan waktu yang berbeda yaitu pada pukul 08.00 WIB dan pukul 12.00 WIB. Data input kondisi batas simulasi yakni suhu dinding kanan, dinding kiri, dinding depan (pintu), lantai dan atap masing-masing skenario terlampir pada Lampiran 4. Nilai suhu rata-rata pengukuran dengan suhu simulasi hampir sama, Perbedaan nilai simulasi dan nilai suhu pengukuran dinyatakan dalam error. Suhu rata-rata pengukuran terbesar pada skenario 5 sebesar 47.27oC standar deviasi pengukuran sebesar 1.64oC dan nilai suhu simulasi sebesar 45.62oC dengan standar deviasi simulasi 2.21oC. Besarnya nilai suhu simulasi dan suhu pengukuran dipengaruhi oleh besarnya iradiasi saat pengukuran dilapangan. Besar standar deviasi simulasi mempengaruhi warna tampilan dihasilkan saat simulasi. Hasil simulasi menghasilkan nilai standar deviasi yang lebih besar dibandingkan dengan nilai standar deviasi hasil pengukuran. Hasil pengukuran didasarkan pada 6 skenario pengukuran hasil standar deviasi terkecil, sedangkan pada simulasi didasarkan pada kualitas mesh pembentukan rak pengering sebesar 60 mm, kualitas mesh yang diperoleh

22 sebesar 0.70659. Nilai kualitas mesh dihasilkan lebih kecil dari 0.85, ini berarti meshing yang sudah dibuat sudah baik ( Tuakia F 2008). Gambar 11-16 dibawah ini memperlihatkan perbedaan distribusi suhu rak yang dinyatakan dalam tingkatan warna yang dihasilkan saat simulasi untuk masing-masing skenario dengan posisi titik pengukuran pada masing-masing skenario. Sebaran suhu udara di rak pengering dipengaruhi oleh kecepatan udara yang masuk ke ruang pengering (inlet), kecepatan udara yang keluar dari ruang pengering (outlet). Aliran udara di dalam ruang pengering memiliki besaran yang digambarkan dari nilai kecepatan aliran udara dan juga memiliki arah aliran yang digambarkan dengan degradasi warna pada ruang pengering. Hasil simulasi CFD dapat diyatakan dalam tampilan kontur dan vektor, vektor diyatakan untuk mengetahui pergerakkan udara yang berada di dalam ruang pengering yakni pergerakan udara yang masuk dari inlet dan keluar menuju otlet. Kecepatan angin dan besarya iradiasi matahari terlampir pada Lampiran 5. Kontur pada CFD ditampilkan untuk melihat gradasi warna pada rak pengering agar mengetahui besarnya panas yang diterima rak pada saat pengukuran. Keuntungan dari pada kontur pada penelitian ini yakni mengetahui besarnya suhu titik-titik pengukuran selain pada titik pengukuran di lapangan, sehingga dapat memprediksi besarnya suhu rak pengering dalam kondisi yang berbeda dengan waktu pengukuran yang berbeda untuk masing-masing skenario. Nilai boundary conditions untuk masing-masing skenario berbeda sesuai dengan nilai hasil pengukuran dilapangan dengan perhitungan yakni nilai heat flux, poros jump, heat transfer, velocity inlet dan outflow. Nilai poros jump pada semua skenario sama hal ini dikaranakan rak yang pergunakan sama untuk semua skenario, sedangkan nilai heat flux, heat transfer, velocityinlet dan outflow berbeda sesuai dengan perhitungan analisis data pengukuran di lapang. Berikut adalah gambar kontur hasil simulasi pada pengukuran pukul 08.00 WIB.

Gambar 11 Simulasi Skenario 1

23




24


Gambar 16 Simulasi Skenario 6 Hasil simulasi pada pukul 08.00 WIB menunjukkan hasil warna yang cukup berbeda, skenario 1 warna simulasi yang dihasilkan berwarna orange mendominasi seluruh bidang rak. Skenario 2 dan 3 warna simulasi berwarna kuning, mendominasi seluruh bidang rak. Pada skenario 5 suhu simulasi yang dihasilkan berwarna orange kurang merata dikarenakan nilai standar deviasi yang dihasilkan pada saat simulasi lebih besar dibanding nilai standar deviasi skenario yang lain. Skenario 4 warna simulasi yang dihasilkan berwarna biru. Warna biru memiliki suhu yang lebih rendah dibanding simulasi berwarna kuning dan orange. Pada skenario 6 warna simulasi yang dihasilkan berwarna hijau, warna hijau mendominasi bagian rak. Warna simulasi yang dihasilkan menampilkan besarnya suhu yang dihasilkan dari visualisasi tampilan kontur temperatur yang dimiliki Fluent 6.3 dengan nilai inputan kondisi batas pada saat melakukan iterasi. Nilai suhu kondisi batas pada Lampiran 4 pengukuran pukul 08:00 WIB.

25 Hasil Simulasi Pukul 12.00 Suhu rata-rata pengukuran sebesar 75.68oC dan suhu simulasi sebesar 72.5oC. Perbedaan suhu yang terukur dan suhu simulasi yang dihasilkan diyatakan pada nilai error yang dihasilkan. Nilai error yang dihasilkan tidak melebihi 10% maka simulasi ini dapat dikatakan baik. Gambar 17-22 dibawah ini menunjukkan perbedaan suhu simulasi masingmasing skenario.



26




27

Gambar 22 Simulasi Skenario 6 Hasil simulasi pengukuran pukul12.00 WIB menampilkan skenario 1 hasil simulasi berwarna orange, warna yang dihasilkan tidak merata hal ini menunjukkan sebaran suhu rak pengering kurang merata, dibanding hasil simulasi pada pukul 08.00 WIB dengan perlakuan rak dalam keadaan diam (tidak diputar), dengan standar deviasi simulasi sebesar 2.21oC merupakan nilai standar deviasi terbesar dari semua skenario. Skenario 2 dan 3 memiliki suhu sebaran rak yang merata yakni berwarna orange, dengan perlakuan rak diputar 45o. Skenario 4 sebaran suhu yang dihasilkan kurang merata, warna sebaran rak yang dihasilkan berwarna kuning dan orange kurang merata, suhu orange lebih tinggi dibanding suhu yang berwarna kuning dengan pemutaran rak 90o. Skenario 5 sebaran suhu yang dihasilkan berwarna kuning merata rak diputar 135o. Skenario 6 sebaran suhu yang dihasilkan berwarna hijau merata dengan perlakuan rak diputar 180o. Nilai warna pada pengukuran pukul 08.00 WIB dan pukul12.00 WIB berbeda hal ini sesuai dengan tingkatan suhu yang dihasilkan pada tingkat visualisai pada Fluent 6.3. Nilai kondisi batas saat melakukan simulasi terlampir pada Lampiran 4 untuk pengukuran pukul 12.00 WIB. Validasi Simulasi CFD Validasi diperlukan untuk mengetahui tingkat error yang terjadi pada hasil simulasi terhadap hasil pengukuran. Validasi pada kondisi tanpa beban menggunakan data pada hasil pengukuran pertama yang sudah terlampir pada tabel validasi terlampir pada Lampiran 2 pengukuran pukul 08.00 WIB dan Lampiran 3 pengukuran pukul 12.00 WIB. Tabel 7 dan Tabel 8 menunjukkan distribusi rak hasil simulasi CFD pada pukul 08.00 WIB dan pukul 12.00 WIB dengan nilai error rata-rata skenario.

28 Tabel 7. Distribusi suhu rak hasil simulasi CFD pukul 08.00 Perlakuan Waktu (Pukul) Skenario 1 Skenario 2 Skenario 3 Skenario 4 Skenario 5 Skenario 6

Suhu pengukuran (oC) 08.00

Suhu simulasi (oC)

Error (%)

CV (%)

08.00

08.00

08.00

42.78 39.85 38.15 39.90 47.27 40.10

49.40 46.40 46.40 31.30 45.60 42.30

15.47 16.43 21.56 21.55 3.53 5.48

9.67 12.11 12.82 9.03 4.72 4.39

Tabel 8. Distribusi suhu rak hasil simulasi CFD pukul 12.00 Perlakuan Waktu (Pukul) Skenario 1 Skenario 2 Skenario 3 Skenario 4 Skenario 5 Skenario 6

Suhu pengukuran (oC) 12.00

Suhu simulasi (oC)

Error (%)

CV (%)

12.00

12.00

12.00

67.08 68.48 67.80 69.41 74.68 70.56

68.3 72.9 72.9 75.6 67.0 65.9

1.81 6.40 6.27 7.55 6.7 7.55

2.52 8.11 7.74 9.73 9.65 7.78

Nilai error rata-rata yang dihasilkan pada masing-masing skenario tidak lebih dari 10% maka simulasi pada masing-masing skenario baik. Menurut Widodo (2009), error yang terjadi pada sebuah simulasi CFD dapat diakibatkan karena kondisi lingkungan yang berubah seperti iradiasi matahari. Nilai error pada skenario 1,2, dan 3 pengukuran pukul 08.00 WIB lebih dari 10%, hal ini dikarenakan kondisi cuaca pada saat pengukuran pada sinar matahari yang mempengaruhi nilai heat flux pada saat melakukan simulasi. Nilai radiasi matahari pada masing-masing dinding pengering yang beragam sehingga mengakibatkan adanya dinding yang tidak terkena sinar matahari pada pagi hari. Pada pagi hari nilai iradiasi matahari mempengaruhi suhu yang diserap oleh dinding pengering, dinding pengering yang tidak terkena matahari akan memperoleh panas yang lebih rendah. Dinding depan, atap memperoleh panas yang lebih tinggi dibanding dinding kiri dan dinding kanan, sedangkan lantai akan menyerap panas yang lebih cepat dikarenakan lantai dicat berwarna hitam dan bersifat menyerap panas (absorber). Nilai error pengukuran pukul 12.00 WIB tidak melebihi dari 10% hal ini dikarenakan penerimaan sinar matahari. Tingkat keragaman pada hasil simulasi menunjukkan nilai koefisien variasi (CV) maka kesergaman rak dapat dikatakan baik. Keragaman yang terbesar atau nilai

29 koefisien variasi terkecil dan sebaran yang paling seragam menjadi rekomendasi hasil penelitian yakni pada skenario 5, dihasilkan dari grafik sebaran suhu rata-rata rak dan standar deviasi yang dihasilkan

SIMPULAN DAN SARAN Simpulan Berdasarkan data yang diperoleh dari penelitian ini dapat diambil kesimpulan yakniputaran rak yang tepat berdasarkan nilai standar deviasi terendah, pengukuran dihasilkan pada pengukuran skenario 5 dengan perlakuan pemutaran rak sejauh 135o lama pemutaran rak selama 1 jam sekali menghasilkan suhu yang lebih seragam terlihat pada grafik dan nilai rata-rata standar deviasi yang dihasilkan masing-masing skenario. Nilai standar deviasi terendah merupakan skenario terbaik pada skenario 5 sebesar 1.67oC. Simulasi dengan teknik CFD pada model pengering ERK dilakukan dengan kondisi steady untuk mengamati pola distribusi suhu rak di dalam ruang pengering. Simulasi CFD memberikan gambaran suhu rak seperti yang diharapkan yakni menunjukkan suhu rak saat pengukuran cenderung sama dengan hasil simulasi yang ditampikan dengan kontur simulasi Fluent 6.3. Validasi nilai hasil pengukuran simulasi dengan nilai hasil pengukuran dilapangan menghasilkan nilai error lebih kecil 10% pada pengukuran pukul 12.00 WIB maka dapat dikatakan baik. Saran Diperlukan penelitian untuk membandingkan sebaran suhu rak yang dihasilkan untuk ERK (Efek Rumah Kaca) tipe rak berputar dengan supply panas berasal dari radiator dengan menggunakan CFD dengan melakukan perhitungan pada masing-masing suhu boundary conditions pada kondisi keadaan lingkungan sekitar.

DAFTAR PUSTAKA Abdullah K. 1999. Recent development of GHE solar drying in Indonesia Grassroot Project. Di dalam: ADC’99, editor. Proceedings of the First Asian-Australia Drying Conference: Bali, Indonesia. Bali (ID): [No halaman tidak diketahui]. Brooker DB, Bakker-arkema FW, Hall CW. 1974. Drying Cereal Grain. Westport (USA). The AVI Publishing Co. Inc. Cengel YA. 2003. Heat Transfer. New York (USA): Mc-Graw Hill Companies.

30 Fellow PJ. 2001. Food Processing Technology Principles and Practices. New York (USA): CRC Press. Hall CW. 1957. Processing Equipment for Agricultural Product. Michigan (USA): Edward Brothers Inc. Arbor. Henderson MS, Perry ME. 1976. Agricultural Process Engineering. Third Edition. Wesport (USA): The AVI Publishing Company. Holman JP. 1997. Perpindahan Kalor. Penerjemah: Jasjfi, E, editor. Jakarta (ID): Penerbit Erlangga. Terjemahan dari: Heat Transfer. Edisi ke-6. Kamaruddin A, Thamrin, F. Wenur. D. Wulandani. 1994. Optimasi dalam perencanaan alat pengering hasil pertanian dengan energi surya. Laporan Akhir Penelitian Hibah Bersaing I. Bogor (ID): Ditjen Dikti, Departemen Pendidikan dan Kebudayaan IPB. Manalu LP, Kamaruddin A. 2001. Model simulasi proses pengeringan kakao memakai pengering surya tipe efek rumah kaca. Buletin Keteknikan Pertanian: 15 (Nomor 3). Mark JE. 2007. Physical properties of polymers. [Internet]. Second edition; [diunduh 2013 Mei 1]. Tersedia pada: http://bib.convdocs.org/docs/24/23048/conv_1/file1.pdf. Mursalim. 1995. Uji performansi sistem pengering energi surya dan tungku batubara dengan bangunan tembus cahaya sebagai pembangkit panas untuk pengering Vanilla planifolia [Tesis]. Bogor (ID): Institut Pertanian Bogor. Nelwan LO. 1997. Pengeringan kakao dengan energi surya menggunakan rak pengering dengan kolektor tipe efek rumah kaca [Tesis]. Bogor (ID): Institut Pertanian Bogor. 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 (CFD) [Tesis]. Bogor (ID): Institut Pertanian Bogor. Tuakia F. 2008. Dasar-Dasar Computational Fluid Dynamics Menggunakan Fluent. Bandung (ID): Informatika. Versteeg HK dan W. Malalasekera. 1995. An Introduction to Computational Fluid Dynamics. Malaysia (MY): Longman Sc. Widodo P. 2009. Kajian pola sebaran aliran udara panas pada model pengering efek rumah kaca hibrid tipe rak berputar menggunakan computational fluid dynamics (CFD) [Tesis]. Bogor (ID): Institut Pertanian Bogor. Wulandani, D. 1997. Analisis pengeringan pada alat pengering Coffe sp efek rumah kaca berenergi surya [Tesis]. Bogor (ID): Institut Pertanian Bogor. Wulandani, D. 2005. Kajian distribusi suhu, RH dan aliran udara pengering untuk optimasi disain pengering efek rumah kaca [Disertasi]. Bogor (ID): Institut Pertanian Bogor.

31 Wulandani, D 2009. Pengembangan alat pengering efek rumah kaca (ERK) hybrid tipe rak berputar untuk penyeragaman aliran udara. Di dalam Widodo, editor. Prosiding Seminar Hasil-hasil Penelitian IPB. Bogor (ID): Institut Pertanian Bogor.

32 Lampiran 1 Contoh analisis perhitungan parameter yang digunakan untuk simulasi CFD 1. Aliran udara ruang pengering: Suhu inlet = 33.3oC Suhu outlet = 28 oC Suhu fluida operasi = =30.65oC Dp= Re=

= =

=1.875 =41741.934 (laminar) Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

Nilai 1.1423 kg/m3 1.0058 kj/kg oC 0.0263w/ oC 1.898 x10-5 kg/ms 15.89 x10-6 m2/s 0.7076

Sumber: Sifat-sifat udara pada (Table A-15, Cengel) 2. Plat absorber Heat transfer coefficient (h) t =300C (free stream temperature) V = 0.37 m/s tp=47 0C t= = =38.5 oC Re=

=

=62057.586

Nu=0.664Re0.5 Pr0.33= 0.664(62057.586)0.5(0.72771)0.33=118.719 Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

Nilai 1.1324 kg/m3 1.0058 kj/kg oC 0.0265w/ oC 1.6879x10-5 kg/ms 14.89 x10-6 m2/s 0.72771

Sumber: Sifat-sifat udara pada (Table A-15, Cengel)

3. Atap Heat transfer coefficient (h) U=0.367m/s hw= 4.5+2.9u hw= 4.5+2.9(0.37)=5.573W/m2 Nilai Heat Flux = ((1000/7)*3 = 428.571 W/m2

33 Lampiran 1 (Lanjutan) 4. Dinding kanan Heat transfer coefficient (h) Pindah panas konveksi pada dinding tegak diyatakan dalam (Chapman, 1974 dalam Wulandani, 2005) : h= Nu Nu= [0.68-0.67Ra1/4[1+( )(9/16)](-4/9)] Untuk 0109 Ra=Gr.Pr Gr= L3gβ Suhu lingkungan =30 oC Suhu dinding = 43.2 oC Suhu fluida operasi = 36.6 oC tm= =36.6+273=309.6K β= =

=3.2299x10-3

Gr=(0,865)3 (9.8)( 3.229910x-3)

=9.71x108

Ra=Gr.Pr=9.71x10-8 x0.727216=7.06x108 Nu= (0.825-0.387Ra1/6[1+( )(9/16)](-8/27]2 = (0.825-0.387(7.061x08)1/6[1+( =110.7260 h= Nu =110.7260 =1.1678 fluks panas =k

)(9/16)](-8/27]2

=0.0263684 =348.06288 W/m2 Nilai Heat Flux = ((1000/7)*3 = 428.571 W/m2 Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

Nilai 1007 kg/m3 1.0058 kj/kg oC 0.0265w/ oC 1.6879x10-5 kg/ms 1.6700 x10-5 m2/s 0.7272


34 Lampiran 1 (Lanjutan) 5. Dinding kiri Heat transfer coefficient (h) Pindah panas konveksi pada dinding tegak diyatakan dalam (Chapman, 1974 dalam Wulandani, 2005) : h= Nu Nu= [0.68-0.67Ra1/4[1+( )(9/16)](-4/9)] Untuk 0109 Ra=Gr.Pr Gr= L3gβ Suhu lingkungan =30 0C Suhu dinding = 49.8 0C Suhu fluida operasi = 39.9 0C tm= =39.9+273=312.9K β= =

=3.1959x10-3

Gr=(0,865)3 (9.8)( 3.1959x10-3)

=1.01x109

Ra=Gr.Pr=1.01x109x0.728074=1.01x109 Nu= (0.825-0.387Ra1/6[1+( )(9/16)](-8/27]2 = (0.825-0.387(1.01x109)1/6[1+( =123.6966 h= Nu =123.6966 =1.3167 fluks panas =k

)(9/16)](-8/27]2

=0.0266126 =526.92948 W/m2 Nilai Heat Flux = ((1000/7)*3 = 428.571 W/m2 Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

Nilai 1007 kg/m3 1.242 kj/kg oC 0.02661w/ oC 1.7010x10-5 kg/ms 12.89 x10-6 m2/s 0.728074


35 Lampiran 1 (Lanjutan) 6. Dinding depan Heat transfer coefficient (h) Pindah panas konveksi pada dinding tegak diyatakan dalam (Chapman, 1974 dalam Wulandani, 2005) : h= Nu Nu= [0.68-0.67Ra1/4[1+( )(9/16)](-4/9)] Untuk 0109 Ra=Gr.Pr Gr= L3gβ Suhu lingkungan =30 oC Suhu dinding = 43.8 oC Suhu fluida operasi = 36.9 oC tm= =36.9+273=309.9K β= =

=3.1565x10-3

Gr=(0,865)3 (9.8)( 3.1565x10-3)

=1.01x109

Ra=Gr.Pr=1.01x109x0.728074=2.95x108 Nu= (0.825-0.387Ra1/6[1+( )(9/16)](-8/27]2 = (0.825-0.387(2.95x108)1/6[1+( =84.8274 h= Nu

)(9/16)](-8/27]2

=84.8274 =0.884661976 fluks panas β= =0.0266126 =135.57648 W/m2 Nilai Heat Flux = ((1000/7)*3 = 428.571 W/m2 Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

Nilai 1.1423 kg/m3 1.0058 kj/kg oC 0.0266w/ oC 1.1273x10-5 kg/ms 13.09 x10-6 m2/s 0.7280


36 Lampiran 1 (Lanjutan) 7. Dinding belakang Heat transfer coefficient (h) Pindah panas konveksi pada dinding tegak diyatakan dalam (Chapman, 1974 dalam Wulandani, 2005) : Nu= [0.68-0.67Ra1/4[1+( )(9/16)](-4/9)] Untuk 0109 Ra=Gr.Pr Gr= L3gβ Suhu lingkungan =30 oC Suhu dinding = 35.2 oC Suhu fluida operasi = 32.6 oC tm= =32.6+273=305.6K β= =

=3.27225x10-3

Gr=(0,865)3 (9.8)( 3.27225x10-3)

=4.05x108

Ra=Gr.Pr=4.05x108x0.727472=2.95x108 Nu= (0.825-0.387Ra1/6[1+( )(9/16)](-8/27]2 = (0.825-0.387(2.95x108)1/6[1+( = 84.827439 h= Nu =84.827439 =0.884661 fluks panas β= =0.0260724 =135.57648W/m2 Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

)(9/16)](-8/27]2

Nilai 1.15412 kg/m3 1.0058 kj/kg oC 0.0260w/ oC 1.1541 x10-5 kg/ms 13.45 x10-6 m2/s 0.7274


37 Lampiran 1 (lanjutan) 8. Perhitungan parameter kondisi batas pada rak dianggap sebagai porous jump a. Permeabilitas permukaan (α) α= = = 3.9x10-8 b. Koefisien porous jump C2 = = =27021.58 9. Properti udara dalam ruang pengering Sifat Massa jenis Panas jenis Konduktivitas panas Viskositas dinamik Viskositas kinematik Bilangan Prandtl

Nilai 1.14804 kg/m3 1007 kj/kg oC 0.0261w/ oC 1.89132 x10-5 kg/ms 1.6474 x10-6 m2/s 0.7270

Sumber: Sifat-sifat udara pada (Table A-15, Cengel) 10. Properti dinding dengan bahan polikarbonat Kerapatan =1210 kg/m3 Panas jenis = 1250 j/kg0C Konduktivitas panas = 0.21 W/Mk

38 Lampiran 2 Tabel distribusi suhu rak hasil simulasi CFD pada pukul 08.00 Skenario 1 A1 A2 A3 A4 A5 B1 B2 B3 B4 B5 C1 C2 C3 C4 C5 D1 D2 D3 D4 D5 E1 E2 E3 E4 E5 F1 F2 F3 F4 F5 G1 G2 G3 G4 G5 H1 H2 H3 H4 H5 Rata-rata (oC)

Suhu Terukur (oC) 42.6 43.2 42.3 42.3 43.4 43 43.6 43 43.3 42.5 4.3 43.2 42.7 43.2 43.3 43.5 43 42.3 43.2 43.3 42.8 42.9 43.2 43.2 44.4 43.2 43.2 42.6 43.2 44.4 42.9 43.8 44.6 43 42.8 42.4 42.9 43 42.4 42.9 41.47

Suhu Simulasi (oC)

τ (oC)

µ(oC)

CV (%)

49.3 49.2 49.4 49.3 49.3 49.3 49.2 49.0 49.2 49.2 49.2 47.4 47.2 48.2 48.2 48.1 49.2 49.2 49.1 49.2 49.2 49.3 49.3 49.4 49.2 49.4 49.1 49.2 49.3 49.2 49.1 49.2 49.3 49.2 49.4 49.3 49.3 49.3 49.3 49.3 49.4

4.73 4.24 5.02 4.94 4.17 4.45 3.95 4.24 4.17 4.73 4.87 2.96 3.18 3.53 3.46 3.25 4.38 4.87 4.17 4.17 4.52 4.52 4.31 4.38 3.39 4.38 4.17 4.66 4.31 3.81 4.38 4.66 4.87 4.38 4.66 4.87 4.52 4.45 4.87 4.52 4.23

0.43 0.43 0.42 0.42 0.44 0.43 0.44 0.43 0.44 0.43 0.42 0.45 0.45 0.44 0.44 0.45 0.43 0.42 0.43 0.44 0.43 0.43 0.43 0.43 0.45 0.43 0.43 0.43 0.43 0.44 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43 0.43

10.96 9.66 11.72 11.53 9.47 10.21 8.93 9.66 9.48 10.96 11.34 6.51 7.03 7.89 7.88 7.19 9.48 11.34 9.48 9.48 10.40 10.40 9.84 9.84 7.52 10.02 9.48 7.52 10.03 8.57 10.03 10.77 11.38 10.03 10.77 11.34 10.40 10.21 11.34 10.40 9.67

39 Lampiran 2 (Lanjutan) Skenario 2

A1 A2 A3 A4 A5 B1 B2 B4 B5 C1 C2 C3 C4 D1 D2 D3 D4 D5 E1 E2 E3 E4 E5 F1 F2 F3 F4 F5 G1 G2 G3 G4 G5 H1 H2 H3 H4 H5 Rata-rata (oC)

Suhu Terukur (oC) 40.5 40.7 40.3 40.3 40.2 39.8 39.8 40.2 40.1 38.8 38.4 38.2 38.9 37.8 37.6 38.2 38.4 38.9 38.6 38.7 38.9 39.8 38.9 39.3 39.7 38.7 38.6 38.7 38.7 38.9 38.9 38.9 38.1 38.7 38.7 38.6 38.7 38.7 43.0

Suhu Simulasi (oC) 46.4 46.3 46.7 45.3 46.2 46.3 46.2 46.8 46.9 46.2 46.3 46.1 46.2 46.2 46.1 46.2 46.2 46.1 46.1 46.9 46.2 46.2 46.1 46.2 46.8 46.3 46.2 46.8 46.2 46.3 46.3 46.3 46.2 46.3 46.2 46.1 46.2 46.1 49.0

τ o ( C)

µ(oC)

CV (%)

4.17 3.96 4.53 3.54 4.24 4.60 4.53 4.67 4.81 5.23 5.59 5.59 5.16 5.94 6.01 5.66 5.52 5.09 5.30 5.80 5.16 4.53 5.09 4.88 5.02 5.37 5.37 5.73 5.30 5.23 5.23 5.23 5.73 5.37 5.30 5.30 5.30 5.23 4.23

0.44 0.44 0.43 0.44 0.44 0.43 0.43 0.43 0.43 0.42 0.41 0.41 0.42 0.41 0.41 0.41 0.42 0.42 0.42 0.41 0.42 0.43 0.42 0.43 0.42 0.42 0.42 0.41 0.42 0.42 0.42 0.42 0.41 0.42 0.42 0.42 0.42 0.42 0.42

9.56 9.01 10.49 7.95 9.75 10.69 10.51 10.87 11.25 12.46 13.47 13.48 12.26 14.52 14.74 13.68 13.27 12.07 12.67 14.05 12.26 10.51 12.07 11.47 11.84 12.86 12.86 13.85 12.66 12.46 2.46 12.46 13.89 12.86 12.66 12.67 12.66 12.47 12.11

40 Lampiran 2 (Lanjutan) τ (oC)

µ(oC)

CV (%)

40.5 40.7 40.3 40.2 40.5 39.8 39.8 40.2 40.2 40 38.8 40.7 38.7 39.4 38.7 38.3 39.7 38.9 39.8 39.8 38.7 38.7 36.7 38.3 37.05 38.8 38.79 38.1 38.7 38.6 38.9 38.7 38.6 38.4 38.65 38 38.7 38.8 38.7 38.7

Suhu Simulasi (oC) 46.3 46.8 46.7 46.9 46.9 46.8 46.4 46.2 46.2 46.89 46.4 46.8 46.9 49.6 46.9 46.7 46.2 46.1 46.7 46.4 46.1 46.8 46.3 46.1 46.5 46.78 46.31 46.9 46.3 46.1 46.2 46.8 46.7 46.78 46.8 46.8 46.7 46.9 46.8 46.8

4.10 4.31 4.53 4.74 4.53 4.95 4.67 4.24 4.24 4.87 5.37 4.31 5.80 7.21 5.80 5.94 4.60 5.09 4.88 4.67 5.23 5.73 6.79 5.52 6.68 5.64 5.32 6.22 5.37 5.30 5.16 5.73 5.73 5.93 5.76 6.22 5.66 5.73 5.73 5.73

0.44 0.43 0.43 0.43 0.43 0.43 0.43 0.44 0.44 0.43 0.42 0.43 0.41 0.40 0.41 0.41 0.43 0.42 0.43 0.43 0.42 0.41 0.40 0.42 0.40 0.41 0.42 0.41 0.42 0.42 0.42 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41

9.38 9.92 10.49 11.05 10.48 11.64 10.88 9.75 9.75 11.42 12.85 9.92 14.05 18.16 14.05 14.48 10.70 12.07 11.45 10.88 12.47 13.85 17.13 13.28 16.77 13.61 12.70 15.32 12.86 12.67 12.26 13.85 13.83 14.44 13.96 15.33 13.65 13.85 13.85 13.85

38.15

46.4

4.24

0.42

12.82

Skenario 3

SuhuTeruku r(oC)

A1 A2 A3 A4 A5 B1 B2 B3 B4 B5 C1 C2 C3 C4 C5 D1 D2 D3 D4 D5 E1 E2 E3 E4 E5 F1 F2 F3 F4 F5 G1 G2 G3 G4 G5 H1 H2 H3 H4 H5 Rata-rata (oC)

41 Lampiran 2 (Lanjutan) Skenario 4 A1 A2 A3 A4 A5 B1 B2 B3 B4 B5 C1 C2 C3 C4 C5 D1 D2 D3 D4 D5 E1 E2 E3 E4 E5 F1 F2 F3 F4 F5 G1 G2 G3 G4 G5 H1 H2 H3 H4 H5 Rata-rata (oC)

Suhu Terukur (oC) 41.30 41.20 41.20 41.25 40.40 40.50 41.90 40.93 40.00 40.20 40.10 39.70 39.40 39.80 39.63 39.30 39.80 39.80 39.63 39.60 39.70 39.70 39.67 39.10 39.80 39.04 39.31 39.00 39.40 39.70 39.37 38.90 38.90 38.90 38.10 38.70 38.70 38.60 38.70 38.70

Suhu Simulasi (oC) 31.30 31.20 31.90 31.24 31.24 31.40 31.50 31.80 31.90 31.89 31.70 31.67 31.78 31.67 31.67 31.89 31.76 31.78 31.67 31.67 31.67 31.67 31.98 31.89 31.78 31.78 31.89 31.67 31.89 31.78 31.78 31.76 31.56 31.67 31.67 31.54 31.40 31.50 31.60 31.58

τ (oC)

µ(oC)

CV (%)

7.07 7.07 6.58 7.08 6.48 6.43 7.35 6.46 5.73 5.88 5.94 5.68 5.39 5.75 5.63 5.24 5.69 5.67 5.63 5.61 5.68 5.68 5.44 5.10 5.63 5.14 5.22 5.11 5.33 5.66 5.36 5.07 5.13 5.15 4.52 5.05 5.12 5.05 5.01 5.04

0.66 0.66 0.65 0.66 0.65 0.64 0.67 0.64 0.63 0.63 0.63 0.63 0.62 0.63 0.63 0.62 0.63 0.63 0.63 0.63 0.63 0.63 0.62 0.61 0.63 0.61 0.62 0.62 0.62 0.62 0.62 0.61 0.62 0.61 0.60 0.61 0.62 0.61 0.61 0.61

10.72 10.71 10.18 10.72 10.02 9.98 11.06 10.03 9.14 9.32 9.39 9.06 8.69 9.15 9.00 8.50 9.07 9.06 9.00 8.97 9.06 9.06 8.76 8.32 9.06 8.36 8.52 8.42 8.60 8.97 8.66 8.24 8.42 8.32 7.56 8.25 8.38 8.19 8.20 8.22

43.08

49.0

4.26

0.63

9.03


Suhu Terukur (oC) 49.8 49.7 49.8 49.8 49.5 49.6 49.7 49.7 49 48.6 48.8 48.6 47.6 47.3 48.3 48.6 47.8 47.3 47.3 47.7 46.5 46.3 46.7 46.7 47.3 46 46.1 46.3 46.1 45 45.2 45.2 45.8 38.4 38.65 38 38.7 38.8 38.7 38.7

Suhu Simulasi (oC) 45.6 45.4 45.3 45.9 45.3 45.3 45.4 45.2 45.2 45.89 45.4 45.8 45.9 45.6 45.9 45.7 45.2 45.1 45.7 45.4 45.1 45.8 45.3 45.1 45.5 45.78 45.31 45.9 45.3 45.1 45.2 45.8 45.7 45.78 45.8 45.8 45.7 46.9 46.8 46.8

τ (oC)

µ(oC)

CV (%)

2.97 3.04 3.18 2.76 2.97 3.04 3.04 3.18 2.69 1.92 2.40 1.98 1.20 1.20 1.70 2.05 1.84 1.56 1.13 1.63 0.99 0.35 0.99 1.13 1.27 0.16 0.56 0.28 0.57 0.07 0.00 0.42 0.07 5.22 5.06 5.52 4.95 5.73 5.73 5.73

0.55 0.55 0.55 0.54 0.55 0.55 0.55 0.55 0.54 0.53 0.54 0.53 0.52 0.52 0.53 0.53 0.53 0.52 0.52 0.53 0.52 0.51 0.52 0.52 0.52 0.50 0.51 0.50 0.51 0.50 0.50 0.49 0.50 0.42 0.42 0.41 0.42 0.41 0.41 0.41

5.44 5.55 5.79 5.08 5.44 5.55 5.55 5.79 4.96 3.62 4.47 3.73 2.32 2.32 3.23 3.86 3.48 2.97 2.19 3.10 1.92 0.70 1.92 2.19 2.45 0.31 1.10 0.56 1.11 0.14 0.00 0.86 0.14 12.44 11.98 13.30 11.69 13.85 13.85 13.85

38.15

46.4

4.24

0.51

4.72


Suhu Terukur (oC) 42.1 42.1 42.1 42.3 42.4 41.3 41.3 41.3 41.4 42.5 40.4 40.4 40.4 40.5 40.1 40.2 40.9 40.7 40.8 39.8 39.4 39.7 39.7 39.6 39.2 39 39.6 39.7 39.7 39.1 38.8 38.9 38.7 38.4 39.3 39.2 39.7 39.7 39.8 39.60

Suhu Simulasi (oC) 42.40 42.98 42.79 42.67 42.67 42.67 42.67 42.89 42.87 42.78 42.45 42.60 42.89 42.50 42.50 42.67 42.78 42.67 42.67 42.68 42.67 42.68 42.56 42.68 43.50 43.78 43.78 43.68 43.60 44.50 44.67 44.12 44.12 43.23 43.23 43.56 43.67 43.56 43.56 42.50

τ (oC)

µ(oC)

CV (%)

0.21 0.62 0.49 0.26 0.19 0.97 0.97 1.12 1.04 0.20 1.45 1.56 1.76 1.41 1.70 1.75 1.33 1.39 1.32 2.04 2.31 2.11 2.02 2.18 3.04 3.38 2.96 2.81 2.76 3.82 4.15 3.69 3.83 3.42 2.78 3.08 2.81 2.73 2.66 2.05

0.50 0.49 0.49 0.50 0.50 0.48 0.48 0.48 0.48 0.50 0.48 0.47 0.47 0.48 0.47 0.47 0.48 0.48 0.48 0.47 0.46 0.47 0.47 0.46 0.45 0.45 0.45 0.45 0.46 0.44 0.43 0.44 0.44 0.44 0.45 0.45 0.45 0.46 0.46 0.47

0.43 1.27 0.99 0.53 0.38 2.00 2.00 2.34 2.15 0.40 3.05 3.28 3.74 2.97 3.60 3.71 2.78 2.92 2.77 4.37 5.01 4.53 4.34 4.69 6.75 7.59 6.54 6.19 6 8.69 9.56 8.37 8.74 7.69 6.11 5 6.18 5.99 5.62 4.40

42.1

42.40

0.21

0.47

4.39

44 Lampiran 3 Tabel distribusi suhu rak hasil simulasi CFD pada pukul 12.00 Skenario 1 A1 A2 A3 A4 A5 B1 B2 B3 B4 B5 C1 C2 C3 C4 C5 D1 D2 D3 D4 D5 E1 E2 E3 E4 E5 F1 F2 F3 F4 F5 G1 G2 G3 G4 G5 H1 H2 H3 H4 H5 Rata-rata (oC)

Suhu Terukur (oC) 69.3 69.2 69.2 69.5 69.6 68.9 68.5 68.6 69.3 69.1 69.6 69.8 69 69.2 69 69 69.7 69 68.9 68.6 68.3 68.2 68 67.3 67.2 65.4 65.3 64.3 65.3 64 64.2 64 64.3 64.3 64.1 64.2 63.2 63.1 63.1 62

Suhu Simulasi (oC) 68.90 68.70 68.70 68.00 68.70 68.32 68.23 68.90 68.43 68.72 68.32 68.70 68.30 68.32 68.20 68.32 63.20 68.70 67.23 68.12 68.12 68.10 66 66.23 66.24 67.23 65.30 65.78 67.23 65.23 67.24 67.14 67.13 67.89 67.23 67.87 67.34 67.13 67.13 68.90

τ (oC)

µ(oC)

CV (%)

0.28 0.35 0.35 1.06 0.64 0.41 0.19 0.21 0.62 0.27 0.91 0.78 0.49 0.62 0.57 0.48 4.60 0.21 1.18 0.34 0.13 0.07 1.19 0.76 0.68 1.29 0.00 1.05 1.36 0.87 2.15 2 2.00 2.54 2.21 2.60 2.93 2 4.85 4.88

0.50 0.50 0.50 0.51 0.51 0.50 0.50 0.50 0.51 0.50 0.51 0.51 0.51 0.51 0.51 0.50 0.55 0.50 0.51 0.50 0.50 0.50 0.51 0.51 0.51 0.49 0.50 0.49 0.49 0.49 0.48 0.48 0.48 0.47 0.48 0.47 0.47 0.47 0.47 3.45

0.56 0.70 0.70 2.08 1.26 0.81 0.38 0.43 1.21 0.53 1.78 1.53 0.98 1.23 1.12 0.95 8.34 0.42 2.30 0.67 0.25 0.14 2.32 1.49 1.34 2.66 0.00 2.14 2.81 1.77 4.50 4.66 4.18 5.36 4.64 5.49 6.24 6.06 6.06 10.84

67.08

68.3

1.23

0.52

2.52


Suhu Terukur (oC) 67.8 67.5 67.4 67.4 66.7 66.5 66.4 66.2 65.2 65.3 65 65.23 66.3 66.2 66.3 66.12 67.4 67.8 67 66.12 68.9 68.9 68.3 66.12 69.8 69.8 70 66.24 68.9 68.3 68.8 68.4 68.3 68.2 68.12 68.2 68.12 68.76 67.2 67.3

Suhu Simulasi (oC) 72.90 72.30 72.13 72.30 72.34 72.30 72.19 72.90 72.13 71.23 72.30 71.34 71.34 71.34 71.23 72.34 71.23 71.24 72.34 72.13 72.13 72.39 72.19 73.12 73.12 73.12 78.90 76.87 75.90 74.30 73.40 72.34 72.39 72.34 72.90 72.89 73.20 72.90 72.90 71.00

τ (oC)

µ(oC)

CV (%)

3.61 3.39 3.34 3.46 3.99 4.10 4.09 4.74 4.90 4.19 5.16 4.32 3.56 3.63 3.49 4.40 2.71 2.43 3.78 4.25 2.28 2.47 2.75 4.95 2.35 2.35 6.29 7.52 4.95 4.24 3.25 2.79 2.89 2.93 3.38 3.32 3.59 2.93 4.03 2.62

0.47 0.47 0.47 0.47 0.46 0.46 0.46 0.45 0.45 0.46 0.45 0.46 0.46 0.46 0.47 0.46 0.47 0.48 0.46 0.46 0.48 0.48 0.47 0.45 0.48 0.48 0.44 0.43 0.45 0.46 0.47 0.47 0.47 0.47 0.47 0.47 0.47 0.47 0.46 0.47

7.76 7.27 7.16 7.43 8.65 8.92 8.90 10.43 10.84 9.15 11.48 9.45 7.67 7.83 7.49 9.62 5.72 5.11 8.15 9.27 4.78 5.19 5.81 10.95 4.92 4.92 14.19 17.45 10.91 9.23 6.94 5.89 6.13 6.21 7.23 7.09 7.72 6.21 8.74 5.52

68.48

72.90

3.74

0.46

8.11


Suhu Terukur (oC) 67.80 67.50 67.40 67.24 67.57 66.70 66.50 66.40 67.30 66.73 65.20 65.30 65.00 67.30 65.70 66.30 66.20 66.30 67.90 66.68 67.40 67.80 67.00 67.89 67.52 68.90 68.90 68.30 68.23 68.58 69.80 69.80 70.00 69.87 69.87 68.90 68.30 68.80 67.2 67.3

Suhu Simulasi (oC) 72.89 72.13 72.13 72.14 72.18 72.30 72.30 72.67 72.54 72.35 72.34 72.35 72.67 72.34 72.13 72.13 72.13 72.35 72.34 72.13 72.13 72.39 72.19 73.12 73.12 73.12 78.90 76.87 75.90 74.30 73.40 72.34 72.39 72.34 72.90 72.89 72.35 72.78 72.89 72.17

τ (oC)

µ(oC)

CV (%)

3.60 3.27 3.34 3.46 3.26 3.96 4.10 4.43 3.71 3.98 5.05 4.99 5.42 3.56 4.55 4.12 4.19 4.28 3.14 3.86 3.34 3.25 3 3.70 3.96 2.98 7.07 6.06 5.42 4.04 2.55 1.80 1.69 1.75 2.14 2.82 2.86 2.81 4.02 3.44

0.47 0.47 0.47 0.47 0.47 0.46 0.46 0.46 0.46 0.46 0.45 0.45 0.45 0.47 0.46 0.46 0.46 0.46 0.47 0.46 0.47 0.47 0.46 0.46 0.46 0.47 0.44 0.44 0.45 0.46 0.48 0.48 0.48 0.48 0.48 0.47 0.47 0.47 0.46 0.47

7.74 7.00 7.16 7.43 6.97 8.58 8.92 9.70 7.99 8.63 11.20 11.05 12.13 7.66 9.98 8.97 9.1 34 6.6 35 7.16 6.93 7.91 7.97 8.57 6.33 16.19 13.64 12.07 8.76 5.35 3.72 3.50 3.62 4.47 5.97 6.07 5.95 8.73 7.39

67.80

72.89

3.60

0.47

7.74


Suhu Terukur (oC) 75.20 75.90 75.67 75.55 70.40 70.80 70.12 70.44 69.90 69.70 69.80 68.40 68.60 68.97 68.50 67.60 67.30 67.89 67.60 66.40 66.30 67.89 66.86 66.20 66.30 66.45 66.32 70.60 70.90 68.57 69.80 69.80 70.00 69.87 69.87 68.90 68.30 68.80 67.20 67.30

Suhu Simulasi (oC) 75.60 75.56 75.67 75.67 75.89 75.70 75.30 75.89 75.89 75.78 75.78 75.32 75.34 75.12 75.13 75.39 75.78 75.89 75.89 75.32 75.24 75.65 75.67 75.9 75.34 75.38 75.30 75.80 75.69 75.34 75.76 75.13 75.14 75.19 75.19 75.36 75.65 75.43 72.89 76.12

τ (oC)

µ(oC)

CV (%)

0.28 0.24 0.00 0.08 3.88 3.46 3.66 3.85 4.24 4.30 4.23 4.89 4.77 4.35 4.69 5.51 6.00 5.66 5.86 6.31 6.32 5.49 6.23 6.86 6.39 6.31 6.35 3.68 3.39 4.79 4.21 3.77 3.63 3.76 3.76 4.57 5.20 4.69 4.02 6.24

0.50 0.50 0.50 0.50 0.46 0.47 0.47 0.46 0.46 0.46 0.46 0.45 0.46 0.46 0.46 0.45 0.44 0.45 0.45 0.44 0.44 0.45 0.44 0.44 0.44 0.44 0.44 0.47 0.47 0.46 0.46 0.46 0.47 0.46 0.46 0.46 0.45 0.46 0.46 0.44

0.57 0.48 0.00 0.17 8.37 7.41 7.87 8.30 9.20 9.35 9.18 10.78 10.47 9.47 10.28 12.29 13.50 12.65 13.17 14.31 14.35 12.23 14.09 15.73 14.53 14.33 14.43 7.90 7.20 10.53 9.15 8.11 7.80 8.10 8.10 9.99 11.51 10.28 8.73 14.11

69.41

75.6

4.40

0.46

9.73


Suhu Terukur (oC) 74.60 74.20 73.00 74.30 74.12 77.20 76.00 77.20 77.12 77.12 73.10 73.10 73.12 73.13 74.10 73.20 73.12 73.12 73.24 77.00 77.13 77.13 77 .12 77.13 77.13 77.14 73.10 73.20 78.90 78.78 73.30 73.00 73.14 73.14 73.78 73.16 73.19 73.17 73.76 74.60

Suhu Simulasi (oC) 67.00 67.54 67.45 67.80 67.34 67.10 67.80 67.87 67.45 67.34 67.23 67.89 67.23 67.23 67.13 67.89 67.65 67.89 67.23 67.12 67.13 67.14 67 67.89 67.89 67.89 67.23 67.13 67.89 67.45 67.13 67.90 67.23 67.13 67.54 67.23 67.14 67.13 67.13 67.89 67.00

τ (oC)

µ(oC)

CV (%)

5.37 4.71 3.92 4.60 4.79 7.14 5.80 6.60 6.84 6.92 4.15 3.68 4.16 4.17 4.93 3.75 3.87 3.70 4.25 6.99 7.07 7.06 6 6.53 6.53 6.53 7.01 4.22 3.75 8.10 8.24 3.82 4.08 4.25 3.96 4.63 4.26 4.29 4.27 4.15 5.37

0.56 0.55 0.54 0.55 0.55 0.58 0.56 0.57 0.57 0.57 0.54 0.54 0.54 0.54 0.55 0.54 0.54 0.54 0.54 0.57 0.57 0.57 0.57 0.57 0.57 0.57 0.57 0.54 0.54 0.58 0.59 0.54 0.54 0.54 0.54 0.55 0.54 0.55 0.5 0.54 0.56

9.65 8.57 7.25 8.39 8.71 12.41 10.35 11.60 11.96 12.08 7.63 6.84 7.66 7.67 8.93 6.96 7.16 6.87 7.80 12.18 12.31 12.30 12.16 11.49 11.50 11.50 12.21 7.75 6.96 13.84 14.04 7.07 7.52 7.80 7.31 8.44 7.81 7.54 7.84 7.64 9.65


Suhu Terukur (oC) 71.9 71.4 71.2 71.3 71.7 71.5 71.4 71.4 71.4 71.5 71 71.2 71.4 69.7 69.5 69.7 69.8 69.12 68.7 68.9 68.1 68.2 68.7 69.7 69.5 69.8 69.8 69.8 70.6 70.4 67.1 67.8 67.9 71.6 71.4 71.2 71.8 71.9 73.76 73.13

Suhu Simulasi (oC) 65.90 65.89 65.89 65.23 65.89 65.12 65.90 65.32 65.89 65.79 65.13 65.89 65.89 65.90 65.13 65.89 65.45 65.87 65.89 65.78 65.98 65.40 65.13 65.14 65.24 65.90 65.80 65.80 65.80 65.80 65.90 65.89 65.90 65.13 65.13 65.13 65.89 65.90 65.80 65.43

τ (oC)

µ(oC)

CV (%)

4.24 3.90 3.75 4.29 4.11 4.51 3.89 4.30 3.90 4.04 4.15 3.75 3.90 2.69 3.09 2.69 3.08 2.30 1.99 2.21 1.50 1.98 2.52 3.22 3.01 2.76 2.83 2.83 3.39 3.25 0.85 1.35 1.41 4.57 4.43 4.29 4.18 4.24 5.63 5.44

0.55 0.54 0.54 0.55 0.54 0.55 0.54 0.55 0.54 0.54 0.55 0.54 0.54 0.53 0.53 0.53 0.53 0.52 0.52 0.52 0.52 0.52 0.53 0.54 0.53 0.53 0.53 0.53 0.54 0.53 0.51 0.51 0.52 0.55 0.55 0.55 0.54 0.55 0.56 0.56

7.78 7.19 6.95 7.85 7.55 8.22 7.18 7.87 7.19 7.43 7.62 6.95 7.19 5.08 5.79 5.09 5.77 4.38 3.81 4.21 2.90 3.80 4.79 6.03 5.66 5.21 5.33 5.33 6.33 6.08 1.67 2.63 2.75 8.32 8.09 7.85 7 7.78 10.04 9.74

70.56

65.90

4.24

0.55

7.78

50 Lampiran 4 Tabel kondisi batas (boundary conditions) pengukuran pukul 08.00 dan pukul 12.00 Skenario 1 Posisi titik pengukuran Waktu (Pukul) Atap Dinding kanan Dinding kiri Dinding depan Dinding belakang Lantai Ruang pengering

Suhu terukur (oC) 08.00 47.5 43.2 43.2 45.5 42.6 46.5 41.1

12.00 67.2 63.2 62.3 63.2 60.2 63.1 63

Skenario 2 Posisi titik pengukuran Waktu (Pukul) Atap Dinding kanan Dinding kiri Dinding depan Dinding belakang Lantai Ruang pengering

Suhu terukur (oC) 08.00 12.00 47.5 67.2 43.2 63.2 43.2 62.3 45.5 63.2 42.6 60.2 46.5 63.1 41.1 63

Posisi titik pengukuran Waktu (Pukul) Atap Dinding kanan Dinding kiri Dinding depan Dinding belakang Lantai Ruang pengering

Suhu terukur (oC) 08.00 12.00 40.9 59.6 40.1 52.8 40.2 57.3 40.2 58.3 39 52.8 40.1 60 39 58

Skenario 3

51 Lampiran 4 (Lanjutan) Skenario 4 Posisi titik pengukuran Waktu (Pukul) Atap Dinding kanan Dinding kiri Dinding depan Dinding belakang Lantai Ruang pengering

Suhu terukur (oC) 08.00 12.00 40.6 60 35.7 50.7 35.9 53.6 36.9 50.3 36.9 50.1 40.5 56.1 35 57


Suhu terukur (oC) 08.00 45.5 41.3 45.5 42.9 42.9 43.2 43

12.00 60 55.7 57.6 55.9 50.3 55.1 55


Suhu terukur (oC) 08.00 40.9 40.8 40.4 40.6 40.8 40.5 40

12.00 60 59 59 59 52 60 58

52 Lampiran 5 Tabelnilai kecepatan angin lingkungan, RH, Kecepatan angin outflow, Iradiasi matahari dan Suhu pengering pada masingmasing skenario. Skenario 1 Waktu (Pukul)

RH (%)

Iradiasi matahari (W/M2)

Kecepatan angin lingkungan (m/s2)

Kecepatan angin outflow kiri (m/s2)

Kecepatan angin outflow tengah (m/s2)


08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 17.00

79.44 58.49 53.43 48.31 51.77 54.18 41.52 44.66 51.38 67

285.71 356.90 346.67 312.15 314.28 316.56 327.3 456.89 457.9 509.78

0.37 0.34 0.15 0.3 0.3 0.6 0.2 0.3 0.5 0.6

2.15 2.44 3.92 2.69 1.52 1.59 2.05 2.05 0.6 2.03

1.16 1.5 3.94 0.48 0.72 0.51 2.17 2.47 1.16 1.3

2.08 2.06 3.97 1.33 1.54 1.02 3.32 2.26 0.6 1.3

Waktu (Pukul)

RH (%)






08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00

86.03 74.46 57.48 83.39 51.2 60.41 51.2 51.99 42.3

571.43 568.00 590.00 782.00 890.42 679.00 768.00 879.00 658.99

0.02 0.1 0.02 0.14 0.13 0.02 0.02 0.02 0.03

0.02 0.2 0.01 0.02 0.02 0.01 0.01 0.2 0.34

0.1 0.05 0.01 0.03 0.03 0.01 0.01 0.21 0.72

0.03 0.1 0.01 0.04 0.02 0.01 0.04 0.23 0.56

Skenario 2

53 Lampiran 5 (Lanjutan) Skenario 3 Waktu (Pukul)

RH (%)




08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00

86.03 63.96 57.43 40.49 44.44 57.48 40.49 40.89 39.40

471.43 514.28 657.14 700.00 760.42 942.85 857.14 714.98 450.20

0.02 0.1 0.02 0.02 0.02 0.14 0.12 0.3 0.12

0.12 0.03 0.1 0.03 0.01 0.03 0.03 0.04 0.04

Kecepatan angin outflow tengah (m/s2) 0.02 0.02 0.2 0.02 0.02 0.03 0.02 0.06 0.03


Skenario 4 Waktu (Pukul)

RH (%)




08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00

64.33 66.89 53.04 51.28 55.90 66.89 84.77 56.24 51.68

450.05 500.00 571.00 785.70 970.42 714.30 624.00 571.40 428.60

0.03 0.10 0.14 0.90 0.10 0.12 0.13 0.24 0.13

0.01 0.04 0.15 0.20 0.14 0.17 0.14 0.30 0.30



54 Lampiran 5 (lanjutan) Skenario 5 Waktu (Pukul)

RH (%)




08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00

74.99 62.09 50.78 55.59 44.44 31.83 31.01 29.86 55.99

400.05 571.42 771.42 771.42 856.00 842.85 771.42 61.28 614.28

0.03 0.1 0.02 0.02 0.04 0.20 0.10 0.02 0.30

0.02 0.20 0.10 0.20 0.02 0.01 0.10 0.31 0.20





Skenario 6 Waktu (Pukul)

RH (%)




08.00 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00

47.4 62.09 34.18 40.89 34.18 31.83 68.33 60.41 68.33

450.05 571.42 642.85 771.42 970.42 771.42 714.28 614.28 571.42

0.01 0.1 0.02 0.02 0.12 0.14 0.20 0.13 0.30

0.01 0.20 0.01 0.20 0.03 0.02 0.12 0.01 0.20

55

RIWAYAT HIDUP Sueritah Henny Wanti, dilahirkan di Tanah Raja Sumatera Utara tanggal 24 Desember 1991 dari ayah H. Sianipar dan ibu T. Simanjuntak, sebagai putri ketiga dari empat bersaudara. Pendidikan formal mulai ditempuh di TK Tanah Raja pada tahun 1995-1997. Selanjutnya, penulis menempuh pendidikan di SD Negeri III Sei rampah pada tahun 1997-2003. Pada tahun 2003-2006 melanjutkan di SMP Negeri 1 Sei Rampah. Pada tahun 2006-2009 melanjutkan di SMA Negeri II Tebing Tinggi. Pada tahun 2009 diterima di IPB melalui jalur Undangan Seleksi Masuk IPB di Departemen Teknik Pertanian yang sekarang berganti nama menjadi Teknik Mesin dan Biosistem, Fakultas Teknologi Pertanian. Selama mengikuti perkuliahan, penulis aktif dalam kepengurusan organisasi intra-kampus diantaranya sebagai Komisi Diaspora PMK IPB, GAMASINTAN dan HIMATETA. Penulis juga pernah menjadi asisten praktikum mata kuliah Teknik Irigasi dan Drainase (TMID), Praktimum Terpadu Mata Kuliah Bahan Petanian (PTMBT), Teknik Green House (TGH) pada tahun 2013. Pada tahun 2012, penulis melaksanakan praktik lapang di PT. Perkebunan Nusantara IV Unit Usaha Adolina-Sumatera Utara dengan judul “Mempelajari Aspek Keteknikan Pertanian pada Proses Pemanfaatan Limbah Kelapa Sawit”. Kemudian penulis melakukan penelitian dan penulisan tugas akhir dengan judul “Penentuan Putaran pada Pengering Surya Tipe Rak Berputar dan Pemodelan Pengering dengan Metode Computational Fluid Dynamics”.

PENENTUAN PUTARAN PADA PENGERING SURYA TIPE RAK BERPUTAR DAN PEMODELAN PENGERING DENGAN METODE COMPUTATIONAL FLUID DYNAMICS SUERITAH HENNY WANTI

Recommend Documents