PENDEKATAN CFD UNTUK OPTIMASI KESERAGAMAN ALIRAN UDARA PADA PENGERING GABAH TIPE BAK
AHMAD FANSURI
DEPARTEMEN TEKNIK MESIN DAN BIOSISTEM FAKULTAS TEKNOLOGI PERTANIAN INSTITUT PERTANIAN BOGOR BOGOR 2013
PERNYATAAN MENGENAI SKRIPSI DAN SUMBER INFORMASI SERTA PELIMPAHAN HAK CIPTA Dengan ini saya menyatakan bahwa skripsi berjudul Pendekatan CFD untuk Optimasi Keseragaman Aliran Udara pada Pengering Gabah Tipe Bak adalah benar karya saya dengan arahan dari 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, Desember 2013 Ahmad Fansuri NIM F14090127
ABSTRAK AHMAD FANSURI. Pendekatan CFD untuk Optimasi Keseragaman Aliran Udara pada Pengering Gabah Tipe Bak. Dibimbing oleh LEOPOLD OSCAR NELWAN. Tujuan utama dari penelitian ini adalah mencari desain pengering yang memiliki keseragaman aliran udara yang optimum. Pendekatan CFD digunakan untuk menganalisis pola aliran udara pada tiga bentuk penampang bak yaitu bujur sangkar, lingkaran, dan persegi panjang. Porositas tumpukan pada CFD ditentukan secara trial and error. Penentuan porositas tersebut berdasarkan nilai koefisien korelasi yang terbaik pada kurva pressure drop hasil simulasi dengan persamaan Shedd sehingga nilai porositas yang diperoleh ialah 30%. Kapasitas bak yang dirancang dapat menampung gabah sebanyak 150 kg pada tingkat ketebalan 50 cm. Berdasarkan simulasi CFD diperoleh kecepatan udara lebih merata pada bentuk penampang lingkaran. Koefisien korelasi pressure drop antara simulasi CFD dengan pengukuran ialah sebesar 0.9957. Koefisien korelasi kecepatan udara antara pengukuran dan simulasi CFD ialah 0.87 sementara nilai error ialah 14.9%. Uji kinerja pengeringan gabah sebanyak 128.5 kg dengan bulk density 555 kg/m3 dimana kadar air awal 28.3% b.b. dilakukan selama 12 jam hingga mencapai kadar air keseimbangan rata-rata 14% b.b. Kata kunci: CFD, gabah, pengering
ABSTRACT AHMAD FANSURI. CFD Approach for Optimization Airflow Uniformity on Flat-Bed Dryer for Paddy. Supervised by LEOPOLD OSCAR NELWAN. The main purpose of this research was to find a dryer design that has the optimum airflow uniformity. CFD approach was used in analyzing the model of airflow on three kinds of chamber shape of flat-bed dryer involving square, circle, and rectangle. The porosity used in CFD was changed by using trial and error method. It’s value was determined based on the best coefficient of determination (COD) between pressure drop resulted from simulation and Shedd equation where the result was 30%. The capacity of flat-bed dryer was designed to accomodate 150 kg paddy with a thickness of 50 cm. The CFD simulation result showed that the circle-shape chamber provided the most uniform airflow pattern. The COD between experimental and simulated pressure drop value was 0.9957. The COD between experimental and simulated air velocity value was 0.87 while the error value was 14.9%. The drying process of 128.5 kg paddy with a bulk density of 555 kg/m3 from the moisture content of 28.3% w.b. to 14% w.b. was carried out for 12 hours. Keywords: CFD, dryer, rice
PENDEKATAN CFD UNTUK OPTIMASI KESERAGAMAN ALIRAN UDARA PADA PENGERING GABAH TIPE BAK
AHMAD FANSURI
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 2013
Judul Skripsi : Pendekatan CFD untuk Optimasi Keseragaman Aliran Udara pada Pengering Gabah Tipe Bak Nama : Ahmad Fansuri NIM : F14090127
Disetujui oleh
Dr Leopold Oscar Nelwan, STP, MSi Pembimbing I
Diketahui oleh
Dr Ir Desrial, MEng Ketua Departemen
Tanggal Lulus:
PRAKATA Puji dan syukur penulis panjatkan kepada Allah subhanahu wa ta’ala atas segala karunia-Nya sehingga karya ilmiah ini berhasil diselesaikan. Tema yang dipilih dalam penelitian yang dilaksanakan sejak bulan Maret 2013 ini ialah merancang dengan pendekatan computational fluid dynamics, dengan judul Pendekatan CFD untuk Optimasi Keseragaman Aliran Udara pada Pengering Gabah Tipe Bak. Terima kasih penulis ucapkan kepada Bapak Dr Leopold Oscar Nelwan, STP, MSi selaku dosen pembimbing, yang telah memberikan bimbingan, arahan, dan motivasi kepada penulis, Dr Ir Usman Ahmad, MAgr dan Dr Muhamad Yulianto, ST, MT sebagai dosen penguji yang telah memberikan arahan dan saran kepada penulis. Penelitian ini adalah bagian dari penelitian Kerja sama Kemitraan Penelitian dan Pengembangan Pertanian Nasional (KKP3N) No. 715/LB.620/I.1/2/2013 berjudul Pengembangan Sistem Pengeringan Gabah Mandiri Energi Menggunakan Sistem Heat Pump Absorpsi dan Pengoperasian Terkendali Berenergi Gasifikasi Sekam yang dipimpin oleh Dr Leopold Oscar Nelwan, STP, Msi. Ungkapan terima kasih juga disampaikan kepada ayah (H. Rojalih Hasan), ibu (Rohimah), kakak (Lina dan Rini), adik (Eli dan Sari) atas segala doa, dukungan, dan kasih sayangnya. Terima kasih juga penulis ucapkan kepada teman-teman satu bimbingan (Dziyad, Kala, dan Angela), teman-teman yang telah membantu selama penelitian (Pijar, Desi, Dian, Tika, Adi SN, dll), para teknisi Departemen TMB (Pak Parma, Pak Darma, Pak Harto dan Mas Firman) serta segenap teman-teman TEP Orion 46 atas semangat dan dukungan yang telah diberikan. Khusus untuk seseorang yang terkasih, terima kasih untuk semua semangatnya. Semoga karya ilmiah ini bermanfaat.
Bogor, Desember 2013 Ahmad Fansuri
DAFTAR ISI DAFTAR TABEL
vii
DAFTAR GAMBAR
vii
DAFTAR LAMPIRAN
viii
PENDAHULUAN
1
Latar Belakang
1
Perumusan Masalah
2
Tujuan Penelitian
2
Manfaat Penelitian
2
Ruang Lingkup Penelitian
2
TINJAUAN PUSTAKA
2
Pengeringan
2
CFD
4
METODE
5
Waktu dan Tempat
5
Bahan
5
Alat
5
Prosedur Penelitian
5
HASIL DAN PEMBAHASAN
16
Simulasi Kecepatan Udara
16
Pemilihan Desain
23
Validasi Pressure Drop
25
Validasi Kecepatan Aliran Udara
28
Hasil Uji Kinerja Pengeringan Gabah
29
SIMPULAN DAN SARAN
32
Simpulan
32
Saran
32
DAFTAR PUSTAKA
32
LAMPIRAN
33
RIWAYAT HIDUP
53
DAFTAR TABEL 1
Data pressure drop CFD, pengukuran dan kurva Shedd
26
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 28 29 30 31 32 33 34
Gambar skematik simulasi 1 Gambar skematik simulasi 2 Gambar skematik simulasi 3 Tahapan pelaksanaan penelitian Domain dan mesh pengering penampang bujur sangkar Domain dan mesh pengering penampang lingkaran Domain dan mesh pengering penampang persegi panjang Skematik lokasi pengukuran kadar air Gambar skematik manometer pipa U Kontur kecepatan udara pada simulasi 1 Vektor distribusi kecepatan udara pada simulasi 1 Vektor distribusi kecepatan udara sisi bawah plenum pada simulasi 1 Vektor kecepatan udara arah sumbu-y posisi outlet simulasi 1 Kontur kecepatan udara pada simulasi 2 Vektor distribusi kecepatan udara pada simulasi 2 Vektor distribusi kecepatan udara sisi bawah plenum pada simulasi 2 Vektor kecepatan udara arah sumbu-y posisi outlet simulasi 2 Kontur kecepatan udara pada simulasi 3 Vektor distribusi kecepatan aliran udara pada simulasi 3 Vektor distribusi kecepatan udara sisi bawah plenum pada simulasi 3 Vektor kecepatan udara arah sumbu-y posisi outlet simulasi 3 Profil kecepatan aliran udara simulasi 1 Profil kecepatan aliran udara simulasi 2 Profil kecepatan aliran udara pada 3 ketinggian simulasi 3 Perbandingan standar deviasi kecepatan udara dalam lapisan tumpukan pada simulasi CFD Standar deviasi rata-rata pada ketiga simulasi CFD Validasi pressure drop CFD dengan kurva Shedd (1953) Validasi pressure drop CFD dengan pressure drop Pengukuran Kontur tekanan pada simulasi 2 Validasi kecepatan aliran udara hasil simulasi CFD dengan hasil pengukuran Penurunan kadar air rata-rata dalam tumpukan gabah Penurunan kadar air rata-rata dalam tiga lapisan tumpukan Suhu udara dalam tiga lapisan tumpukan gabah terhadap waktu RH udara pada tiga lokasi pengukuran terhadap waktu
6 7 8 9 10 10 10 14 15 17 17 18 18 19 19 20 20 21 21 22 22 23 24 24 25 25 27 27 28 28 29 30 31 31
DAFTAR LAMPIRAN 1 2 3 4
Hasil simulasi 1 CFD Hasil simulasi 2 CFD Hasil simulasi 3 CFD Perbandingan keragaman kecepatan aliran udara dan tekanan dalam pengering pada ketiga simulasi CFD 5 Data validasi kecepatan aliran udara hasil pengukuran dan simulasi CFD serta nilai error dan standar deviasinya 6 Hasil pengukuran kadar air bijian selama 16 jam 7 Data kecepatan aliran udara hasil pengukuran dan simulasi CFD pada beberapa variasi kecepatan udara 8 Data pengukuran suhu uji kinerja pengeringan gabah 9 Gambar teknik pengering gabah tipe bak 10 Sistem pengering tipe bak lingkaran yang diujicoba dalam penelitian 11 Alat-alat ukur yang digunakan dalam penelitian
33 35 37 39 40 41 44 45 47 50 52
PENDAHULUAN Latar Belakang Gabah merupakan komoditas pertanian yang sangat penting. Gabah dimanfaatkan sebagai bahan pangan yang diolah menjadi beras. Oleh karena itu, gabah menjadi komoditas perdagangan yang sangat penting sehingga perlu diperhatikan salah satu aspek penting dalam proses pasca panen yakni pengeringan. Pengeringan dilakukan untuk mencegah perkecambahan biji, untuk mempertahankan kualitas bijian, dan untuk mencapai level kadar air dimana tidak memungkinkan bateri dan jamur berkembang (Hall 1970). Proses pengeringan dapat dilakukan menggunakan mesin pengering. Berbagai macam bentuk mesin pengering beredar di masyarakat. Mesin pengering tipe bak adalah salah satu tipe pengering yang sering digunakan dalam proses pengeringan gabah. Salah satu parameter penting dari optimalnya kinerja mesin pengering ialah bentuk bak pengering. Bentuk bak pengering sangat menentukan sebaran aliran udara sehingga berpengaruh pada keseragaman kadar air pada tumpukan bahan. Salah satu permasalahan dalam pengeringan tumpukan ialah keseragaman kadar air yang sulit tercapai. Keseragaman kadar air pada pengeringan tumpukan tebal sulit tercapai antara lapisan bawah, tengah dan atas. Pada penelitian ini akan dirancang bentuk bak pengering yang dapat memberikan sebaran kecepatan udara dan suhu yang seragam pada posisi tengah dan pinggiran bak berdasarkan simulasi CFD (Computational Fluid Dynamics). CFD memprediksi aliran berdasarkan model matematika melalui persamaan diferensial parsial, metode numerik dan tools perangkat lunak (solvers, tools predan postprocessing). CFD adalah alat untuk memprediksi apa yang akan terjadi pada alat atau sistem dengan satu atau lebih kondisi batas. CFD sebenarnya mengganti persamaan-persamaan diferensial parsial dari kontinuitas, momentum, dan energi dengan persamaan-persamaan aljabar. CFD merupakan pendekatan dari persoalan yang asalnya kontinum (memiliki jumlah sel tak terhingga) menjadi model yang diskrit (jumlah sel terhingga) (Tuakia 2008). CFD dijadikan tahap desain skala laboratorium menggunakan software analisis CFD. Desain yang terpilih ialah yang memiliki sebaran kecepatan udara yang lebih merata pada model tumpukan gabah. Pada penelitian ini dibuat suatu model pengeringan gabah tipe bak. Model pengering yang dianalisis diantaranya pengering tipe bak berpenampang bujur sangkar, lingkaran dan persegi panjang. Kapasitas pengering yang dirancang dapat menampung gabah sebanyak 150 kg. Pola keseragaman kecepatan udara menjadi hal yang ingin dicapai pada desain pengering gabah. Melalui analisis CFD dapat diketahui pendekatan kondisi yang sebenarnya jika desain tersebut digunakan untuk pengeringan gabah. Validasi perlu dilakukan untuk mengetahui ketepatan data simulasi dengan data pengujian lapang. Diharapkan dengan simulasi dan rancangan yang dibuat dapat menduga pola sebaran kecepatan udara dalam tumpukan gabah pada pengering sehingga dapat ditentukan desain yang paling baik.
2 Perumusan Masalah Salah satu tahapan yang memerlukan banyak waktu dan biaya dalam pembuatan pengering buatan ialah menentukan bentuk dan analisis perhitungannya. Tujuan yang ingin dicapai dalam desain pengering ialah keseragaman aliran udara pada tumpukan bahan. Namun metode perhitungan biasa tidak dapat menduga bagaimana pola aliran yang terjadi dalam desain pengering. CFD digunakan untuk menganalisis pola dan sebaran aliran udara pada berbagai bentuk penampang bak pengering rancangan. Proses analisis CFD dimulai dengan pembuatan geometri untuk menentukan domain selanjutnya dilakukan komputasi. Hasil yang diperoleh dari analisis ini diantaranya kecepatan udara di titik yang ingin diketahui nilainya serta pola aliran yang terjadi di dalam domain rancangan.
Tujuan Penelitian Penelitian yang dilakukan bertujuan untuk menentukan desain bak pengering gabah tipe bak yang memiliki keseragaman aliran udara yang optimum berdasarkan simulasi CFD.
Manfaat Penelitian Manfaat penelitian ini ialah diperoleh desain dan model pengering gabah tipe bak yang memiliki pola dan sebaran kecepatan aliran udara yang seragam. Kecepatan aliran udara yang seragam berakibat pada penyebaran suhu udara pengering yang merata di seluruh tumpukan.
Ruang Lingkup Penelitian Ruang lingkup dari penelitian ini adalah untuk mendapatkan bentuk rancangan pengering gabah yang memiliki pola aliran udara yang seragam dari tiga bentuk penampang yang disimulasikan. Validasi dilakukan terhadap bentuk penampang yang dipilih untuk membandingkan hasil yang diperoleh dari simulasi dengan hasil pengukuran.
TINJAUAN PUSTAKA Pengeringan Periode pengeringan terdiri dari dua periode utama yaitu periode pengeringan konstan dan periode pengeringan menurun. Dalam periode pengeringan konstan, pengeringan terjadi pada permukaan bijian dan sama dengan penguapan dari kandungan permukaan air bebas. Titik yang menandai akhir dari periode konstan terjadi ketika difusi kelembaban dalam produk menurun di bawah yang diperlukan untuk mengembalikan kelembaban di permukaan. Besarnya
3 tingkat pengeringan periode ini tergantung pada area yang terkena, perbedaan kelembaban permukaan, koefisien transfer massa dan kecepatan udara pengeringan. periode pengeringan menurun masuk setelah periode pengeringan konstan. Kadar air kritis terjadi antara periode konstan dan periode menurun. Kadar air kritis adalah kadar air minimum dari biji-bijian yang mempertahankan aliran air bebas ke permukaan biji-bijian yang sama dengan laju maksimum penguapan uap air dalam biji-bijian di bawah kondisi pengeringan (Hall 1957). Pada pengeringan lapisan tebal penambahan kadar air akan dihilangkan dari lapisan kering hingga tercapai keseimbangan kadar air. Sedikit kadar air akan dihilangkan dan beberapa dapat bertambah ke zona basah sampai zona pengeringan tercapai (Hall 1957). Parameter-parameter pengeringan yang berpengaruh terhadap waktu yang dibutuhkan untuk mengurangi kadar air biji-bijian hingga kadar air yang diinginkan adalah (Hall 1957) : a. Suhu udara pengering Suhu udara pengering yang diukur terdiri dari suhu bola kering dan suhu bola basah. Suhu bola kering adalah suhu udara atau produk yang ditunjukkan oleh termometer yang tidak terpengaruh oleh kandungan uap air dari udara. Suhu bola basah adalah suhu yang ditunjukkan oleh termometer dengan bagian sensor yang ditutupi dengan lapisan tipis air dan udara bergerak melaluinya sampai diperoleh suhu yang stabil selama penguapan. Perbedaan antara suhu bola kering dan bola basah disebut depresi bola basah (Hall 1957). b. Kelembaban relatif udara pengering Kelembaban relatif udara adalah rasio tekanan parsial uap air terhadap tekanan uap jenuh pada suhu bola kering tertentu, biasanya dinyatakan sebagai persentase (Hall 1957). Jika suhu ditingkatkan saat kadar air terjaga konstan maka kelembaban relatif akan menurun. Udara yang memiliki kelembaban relatif rendah paling efektif digunakan untuk pengeringan (Hall 1970). c. Porositas Porositas didefinisikan sebagai fraksi volume dari ruang kosong udara dan disajikan sebagai perbandingan dari volume kosong udara terhadap volume total (Champ 1996). d. Kadar air bahan Jumlah uap air di produk ditetapkan berdasarkan berat air dan biasanya dinyatakan dalam persen. Ada dua metode yang menunjuk kadar air bahan yaitu basis basah dan basis kering. Kandungan uap air pada basis basah diperoleh dengan membagi berat air yang terkandung dalam bahan dengan berat keseluruhan materi. Untuk menghitung kadar air basis basah digunakan persamaan 1 (Hall 1957). Kadar Air (basis basah) =
w w
keterangan : Ww : berat dari air Wd : berat dari bahan kering
d
x 100% .......................................... (1)
4 Persen kelembaban pada basis kering ditentukan dengan membagi berat air dengan berat bahan kering. Untuk menghitung kadar air basis kering digunakan persamaan 2 (Hall 1957). w
Kadar Air (basis kering) = x 100% ............................................... (2) d Basah basis digunakan untuk penunjukan komersial dan juga digunakan oleh Federal Grain Standard. Kadar air basis kering digunakan terutama pada penelitian dan dalam persamaan yang berhubungan dengan variasi kadar air. Oleh karena itu, metode basis kering untuk mengungkapkan kadar air yang digunakan dalam persamaan pengeringan. Kadar air pada basis kering selalu lebih besar daripada basis basah (Hall 1957).
CFD Computational Fluid Dynamics atau CFD adalah analisis sistem yang melibatkan aliran fluida, perpindahan panas dan fenomena terkait seperti reaksi kimia dengan cara simulasi berbasis komputer (Versteeg dan Malalasekera 1995). CFD sebenarnya mengganti persamaan-persamaan diferensial parsial dari kontinuitas, momentum, dan energi dengan persamaan-persamaan aljabar (Tuakia 2008). Pada umumnya terdapat tiga tahapan yang harus dilakukan dalam simulasi CFD menurut Tuakia (2008), yaitu : 1. Preprocessing Preprocessing merupakan langkah pertama dalam membangun dan menganalisis sebuah model CFD. Teknisnya adalah membuat model dalam paket CAD (Computer Aided Design), membuat mesh yang sesuai, kemudian menerapkan kondisi batas dan sifat-sifat fluidanya. 2. Solving Solvers (program inti pencari solusi) CFD menghitung kondisi-kondisi yang diterapkan pada saat preprocessing. 3. Postprocessing Postprocessing adalah langkah terakhir dalam analisis CFD. Hal yang dilakukan pada langkah ini adalah mengorganisasi dan menginterpretasi data hasil simulasi CFD yang bisa berupa gambar, kurva, dan animasi. Prosedur yang terdapat pada semua pendekatan program CFD (Tuakia 2008), 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.
5
METODE Waktu dan Tempat Kegiatan penelitian dilaksanakan mulai bulan Maret 2013 sampai dengan Agustus 2013. Pengambilan data uji kinerja pengering dilakukan di Laboratorium Teknik Energi Terbarukan, Departemen Teknik Mesin dan Biosistem, Fakultas Teknologi Pertanian, IPB.
Bahan Bahan yang digunakan dalam penelitian ini adalah gabah varietas IR 64 dengan kadar air rata-rata 28.33%b.b. sebanyak 128.5 kg.
Alat Peralatan yang digunakan dalam penelitian ini adalah notebook ASUS K42J (Operating System dan Microsoft Office), software design (Solidworks 2011 dan Gambit 2.4.6), software analisis CFD (Ansys Fluent v12.1), hybrid recorder Yokogawa, termokopel tipe CC (Copper Constanta), timbangan digital EK-1000, anemometer Kanomax tipe 6011, blower 1 phase, satu unit pengering hasil rancangan, kawat kasa, sampel picker, oven drying.
Prosedur Penelitian Pengering yang akan didesain memiliki komponen yakni bak, plenum, dan lantai pengering. Kapasitas bak pengering yang dirancang dapat menampung 150 kg gabah dengan ketebalan gabah 50 cm. Bulk density gabah ialah 577 kg/m3 (Hall 1957). Tinggi plenum yang dirancang ialah 18 cm, diameter lubang pemasukan udara ialah 10 cm. m ss V= = ≈ 0.25 m3 l m vo ume . m L= = = 0.5 m2 t . m Keterangan : V : Volume bak (m3) L : Luas alas (m2) Bentuk dasar penampang pengering yang akan dirancang yaitu : 1. Bujur sangkar Skenario simulasi 1 dilakukan pada bentuk penampang bujur sangkar. Gambar skematik rancangan pengering simulasi 1 disajikan pada Gambar 1. Dimensi penampang ditentukan melalui perhitungan dibawah ini : s = √ = √ . m ≈ 0.7 m = 70 cm
6
Keterangan : s : Panjang sisi bujur sangkar (cm)
Gambar 1 Gambar skematik simulasi 1 2. Lingkaran Skenario simulasi 2 dilakukan pada bentuk penampang lingkaran. Gambar skematik rancangan pengering simulasi 2 disajikan pada Gambar 2. Dimensi penampang ditentukan melalui perhitungan dibawah ini : =√
x
=√
Keterangan : : Diameter lingkaran (cm)
= . 9 m ≈ 80 cm
7
Gambar 2 Gambar skematik simulasi 2 3. Persegi panjang Skenario simulasi 3 dilakukan pada bentuk penampang persegi panjang. Asumsi yang digunakan ialah sisi panjangnya berukuran 2 kali sisi lebarnya. Gambar skematik rancangan pengering simulasi 3 disajikan pada Gambar 3. Dimensi penampang ditentukan melalui perhitungan dibawah ini : p=2xƖ L = p x Ɩ = 2 x Ɩ x Ɩ = 2Ɩ2 Keterangan : . m p : panjang sisi (cm) Ɩ =√ =√ = 0.5 m = 50 cm Ɩ : lebar sisi (cm) p = 2 x Ɩ = 2 x 50 cm = 100 cm
8
Gambar 3 Gambar skematik simulasi 3 1.
2.
3. 4.
Tahapan yang dilakukan dalam penelitian ini meliputi 4 tahapan yaitu : Menggambar geometri model menggunakan perangkat lunak Gambit. Tahap selanjutnya ialah melakukan simulasi aliran udara dengan menggunakan perangkat lunak Ansys Fluent. Pembuatan rancang bangun (pabrikasi) pengering. Desain yang dipabrikasi dipilih dari hasil simulasi CFD yang memiliki pola aliran udara yang seragam. Dalam penelitian ini tidak melakukan simulasi CFD mengenai perubahan kadar air yang terjadi pada bahan. Pengujian kinerja pengering meliputi kadar air bahan, kecepatan aliran udara, pressure drop, suhu dan lama pengeringan. Analisis data kecepatan udara dilakukan untuk melihat hasil yang diperoleh dari pengujian lapang dan simulasi CFD. Validasi data bertujuan untuk melihat keakuratan hasil uji lapang dan simulasi CFD. Adapun tahapan perancangan hingga validasi data ditunjukkan oleh Gambar 4.
9
Gambar 4 Tahapan pelaksanaan penelitian 1. Simulasi CFD Parameter yang diukur dalam simulasi CFD ialah kecepatan aliran udara. Simulasi difokuskan pada sebaran pola aliran udara pada tumpukan sehingga asumsi tumpukan berpori gabah dalam simulasi tidak mengandung kadar air. Oleh sebab itu suhu tidak menjadi parameter yang diukur dalam simulasi CFD. Penentuan rancangan pengering menggunakan CFD melalui tahap-tahap sebagai berikut : a. Pembuatan Geometri, Meshing dan Boundary Conditions Geometri yang telah dibuat akan dilakukan proses pembuatan Grid/Mesh. Pembuatan mesh dilakukan menggunakan perangkat lunak Gambit 2.4.6. Mesh/Grid menggunakan element yang tidak terstruktur dengan model tet/hybrid dan tipe Tgrid, interval size : 2. Mesh model tet/hybrid dengan tipe Tgrid sebagian besar terdiri dari elemen tetrahedral, heksahedral, piramida dan wedge. Sebagian besar volume dapat langsung
10 di-mesh dengan tipe ini tanpa harus dipisah menjadi beberapa volume tertentu. Sehingga tahap pemberian mesh dapat menjadi lebih mudah. Selanjutnya grid tersebut dilengkapi dengan kondisi batas. Bentuk domain dan mesh setiap pengering disajikan pada Gambar 5, Gambar 6 dan Gambar 7.
Gambar 5 Domain dan mesh pengering penampang bujur sangkar
Gambar 6 Domain dan mesh pengering penampang lingkaran
Gambar 7 Domain dan mesh pengering penampang persegi panjang Proses selanjutnya penentuan kondisi batas simulasi. Penentuan kondisi batas simulasi bertujuan untuk membatasi bagian yang akan dianalisis oleh perangkat lunak. Kondisi batas yang diterapkan ialah sebagai berikut : 1) Wall atau dinding yakni kondisi batas dinding berfungsi untuk memisahkan antara regional fluida dan solid.
11 2) 3) 4) 5)
Velocity inlet pada lubang pemasukan udara dari blower. Pressure outlet yaitu posisi pengeluaran udara diatas tumpukan. Porous jump pada lantai berlubang. Pembagian dua cell zone yakni cell zone gabah sebagai asumsi tumpukan gabah dan cell zone udara. Pengaturan tentang cell zone dijelaskan pada tahap selanjutnya.
b. Penyelesaian dengan Ansys Fluent Proses pembuatan simulasi dilakukan menggunakan perangkat lunak Ansys Fluent v12.1. Proses analisis terdiri dari tahap-tahap sebagai berikut : 1) General setting Pengaturan diawali dengan memasukkan mesh model lalu penskalaan mesh model dalam cm, mengatur solver type pressure-based, velocity formulation absolute, time steady, nilai gravitasi pada sumbu y = -9.81 m/s2. Skala unit dalam satuan cm. 2) Pemilihan model Model simulasi yang diaktifkan ialah model persamaan energi dan model viskos. Model viskos yang dipilih yaitu model viscous k-epsilon (2 eqn) standard dan near-wall treatment standard wall functions. Model kepsilon dikembangkan oleh Launder & Spalding. Alasan pemilihan model ini karena dinilai cukup stabil, cukup ekonomis, memiliki akurasi yang cukup memadai untuk berbagai jenis aliran turbulen membuat model ini sering digunakan pada simulasi aliran fluida dan perpindahan panas (Tuakia 2008). Parameter model constants pada k-epsilon dipilih default. 3) Penambahan material Material yang ditambahkan ialah gabah yang dimasukkan secara manual berdasarkan rujukan pustaka. Nilai yang dimasukkan ialah : Material type : solid Material name : gabah Density : 577 kg/m3 (Hall 1957) 4) Pengaturan cell zone conditions Zona dalam domain dibagi menjadi dua bagian, yaitu zona udara dan zona gabah (porous zone). Zona gabah merupakan cell fluida yang diasumsikan memiliki porositas tertentu dengan mengaktifkan porous zone pada panel Fluid. Porous formulation dipilih physical velocity agar diperoleh hasil simulasi yang lebih akurat. Input pada porous media ialah : Memilih material fluida yakni udara yang melewati media berpori. Memilih relative velocity resistance formulation. Memasukkan nilai viscous resistance coefficients ( α) d inertial resistance coefficients (C2). Memasukkan nilai porositas dari media berpori. Memilih material media berpori yakni gabah. Penentuan parameter input cell zone gabah ialah sebagai berikut : Viscous resistance coefficients (1/ α) de pers m (F UENT ver.12.1) :
12 ( - . )
= 119,642,868.5 1/m2 p . . Inertial resistance coefficients (C2) dengan persamaan (FLUENT ver.12.1) :
1/α =
=
- . . ( - . ) . C2 = p = . = 19,023.2 1/m . Keterangan : : porositas material Dp : diameter rata-rata gabah = 0.477 cm = 0.00477 m (Margana 2010) Nilai porositas gabah ini ditentukan secara trial and error dengan bertolok ukur pada nilai korelasi kurva pressure drop antara simulasi CFD dengan persamaan pressure drop kurva Shedd (persamaan 8). 5) Pengaturan boundary conditions Zona pada boundary conditions meliputi inlet, outlet, porous jump, interior dan wall. Simulasi menggunakan velocity inlet dan pressure outlet. Inlet dianggap sebagai velocity inlet dengan kecepatan udara masuk 13.2 m/detik. Metode spesifikasi kecepatan yang dipilih ialah components dimana memasukkan nilai kecepatan fluida di sisi masuk dan arah sumbu koordinat (sumbu negatif z). Selanjutnya data thermal dengan nilai temperatur udara masuk sebesar 40oC. Outlet dianggap sebagai pressure outlet. Data yang dimasukkan pada outlet ialah data thermal dengan nilai Backflow Total Temperatur sebesar 30oC, sedangkan nilai tekanan gauge adalah 0 pascal. Lantai berlubang pengering dianggap sebagai porous jump. Pengaturan Porous Jump dengan memasukkan data face permeability (m2), porous medium thickness (cm), dan pressure-jump coefficient (C2) (1/m). Persen lubang pada lantai diasumsikan sebesar 70%. Perhitungan parameter input porous jump ialah sebagai berikut : Permeabilitas permukaan (α) dihitung dengan persamaan (FLUENT ver.12.1) :
α =
p
=
.
.
= 2.289x10-7 m2
( - ) ( - . ) Koefisien porous jump (C2) dihitung dengan persamaan (FLUENT ver.12.1) :
C2 =
. ( - . )
=
.
- .
= 1,020.4 1/m p . . keterangan : : Asumsi persen lubang pada plat lantai = 0.7 Dp : Asumsi diameter lubang pada plat = 3 mm = 0.003 m
13
6)
7)
8)
9)
Udara lingkungan dianggap konstan selama simulasi, yaitu pada 30oC. Suhu inlet : 40oC Suhu Outlet : 30oC Suhu fluida operasi : (40+30)/2 = 35oC Tekanan udara : 1 atm = 101.325 kPa Pengaturan solution methods dan solution controls Pengaturan metode solusi digunakan untuk memilih pola interpolasi yang dilakukan pada node mesh. Skema pressure-velocity coupling ialah SIMPLE dengan pengaturan default pada spatial discretization. Dalam skema ini dipilih first-order upwind scheme dimana merupakan skema interpolasi yang paling ringan dan cepat mencapai konvergen. Pengaturan solution controls menggunakan nilai default. Pengaturan solution monitors Pengaturan monitors pada residuals meliputi penentuan kriteria konvergensi. Kriteria konvergensi menggunakan nilai 0.001 untuk semua persamaan residual dan 10-6 pada persamaan energi. Kriteria konvergensi adalah perbedaan antara tebakan awal dan hasil akhir hasil iterasi. Pengaturan solution initialization Proses inisialisasi dilakukan sebagai langkah awal dengan menghitung semua nilai dari kondisi batas yang telah dimasukkan. Langkah yang dilakukan ialah memilih compute from-inlet. Pengaturan solution run calculation Proses perhitungan dari persamaan dan model yang dipilih dan berhenti hingga tercapai konvergen atau sesuai waktu iterasi yang telah ditentukan. Iterasi dilakukan maksimum sebanyak 1000 kali.
c. Tahapan Penyajian Data Tampilan hasil yang didapat berasal dari Ansys CFD-Post v12.1, meliputi : 1) Hasil geometri yang terbentuk 2) Plot kontur kecepatan dan tekanan 3) Plot vektor kecepatan dan tekanan 4) Nilai titik sampel kecepatan dan tekanan 2. Pembuatan Pengering Pembuatan pengering merupakan proses rancang bangun (pabrikasi). Desain yang dibuat merupakan desain yang memiliki sebaran aliran udara yang paling merata dari ketiga simulasi CFD ditandai dengan nilai standar deviasi rata-rata yang paling kecil. Pembuatan gambar teknik pengering menggunakan software SolidWorks 2011. Gambar teknik pengering disajikan pada Lampiran 9. 3. Uji Kinerja Pengering a. Kadar air bahan Kadar air bahan diukur meliputi kadar air awal dan penurunannya tiap dua jam. Alat yang digunakan untuk mengambil sampel gabah ialah sampel picker. Metode yang digunakan ialah metode langsung dengan oven drying
14 pada suhu konstan 105oC. Berat rata-rata sampel gabah yang diambil ialah 5-6 gram. Massa padatan gabah diperoleh dengan mengeringkan sampel gabah selama 24 jam. Lokasi pengukuran kadar air bahan ialah pada titik tumpukan bawah, tengah dan atas dengan masing-masing 5 pengukuran yang ditunjukkan pada Gambar 8.
Gambar 8 Skematik lokasi pengukuran kadar air, suhu, RH dan pressure gauge Kadar air produk selama pengeringan dihitung berdasarkan : Kadar air (%bb) =
m r m r mp d t
x
Kadar air (%bk) =
m r mp d t
............................................... (4)
x
b. Kecepatan aliran udara Titik pengukuran aliran udara yaitu : 1) Lubang pemasukan udara (inlet),
................................... (3)
15 2) Tumpukan bagian bawah dan tengah sebanyak 4 titik yaitu dengan jarak 15 cm dari tepi bak, serta lapisan atas sebanyak 5 titik pengukuran termasuk diantaranya posisi tengah. c. Suhu Suhu udara masuk diatur pada kisaran 39-45oC. Bagian-bagian yang diukur suhunya ialah udara lingkungan sekitar pengering, inlet, plenum dan di dalam tumpukan gabah. Alat yang digunakan untuk mengukur suhu ialah termokopel. Hybrid recorder digunakan untuk menampilkan data pengukuran termokopel. Penempatan termokopel yaitu : 1) Lubang pemasukan udara (inlet). 2) Tumpukan gabah bagian bawah, tengah, dan atas masing-masing 5 titik pengukuran. Koordinat pengukuran ialah sama seperti lokasi pengambilan sampel kadar air. 3) Udara lingkungan. d. RH Lokasi pengukuran relative humidity (RH) ialah : 1) Plenum. 2) Di atas tumpukan gabah. 3) Udara lingkungan sekitar pengering. Alat yang digunakan untuk mengukur RH ialah termokopel bola basah dan bola kering. Termokopel bola basah dibuat dengan menyelubungi ujung sensor termokopel menggunakan kain yang dijaga agar tetap basah. e. Pressure drop Lokasi pengukuran tekanan udara ialah pada plenum. Alat yang digunakan untuk mengukur tekanan udara ialah manometer pipa U. Nilai tekanan udara pada plenum diukur berdasarkan perbedaan ketinggian permukaan air pada kolom pipa U seperti ditunjukkan pada Gambar 9.
Gambar 9 Gambar skematik manometer pipa U
16 Perbedaan ketinggian ( h dikonversi ke nilai tekanan menggunakan persamaan : p= h................................................................................................. (5) keterangan : p : tekanan gauge (Pascal) : massa jenis fluida (kg/m3) : percepatan gravitasi (m/s2) h : beda ketinggian permukaan fluida di dalam pipa (m) 4. Analisis data Data kecepatan udara dan tekanan yang diperoleh dari hasil pengujian kinerja dibandingkan dengan data hasil simulasi CFD. Validasi dilakukan pada desain yang terpilih dari simulasi CFD. Tujuan dari validasi ini ialah untuk melihat keakuratan data simulasi terhadap data pengujian. Data-data tersebut dibuat korelasi antara kecepatan udara ukur dan kecepatan udara CFD, tekanan udara ukur dan tekanan udara CFD. Perhitungan terhadap nilai kesalahan (error) dilakukan untuk membandingkan pengukuran dan simulasi CFD. Persamaan yang akan digunakan untuk mengukur error sebagai berikut : u ur F Error = x 100% .................................................................... (6) u ur Standar Deviasi : S = √
∑
( -̅)
............................................................ (7)
Keterangan : v ukur : kecepatan udara pengukuran (m/s) v CFD : kecepatan udara simulasi CFD (m/s) : kecepatan udara ke-i (m/s) ̅ : kecepatan udara rata-rata (m/s) n : jumlah data
HASIL DAN PEMBAHASAN Simulasi Kecepatan Udara 1. Aliran kecepatan udara simulasi 1 Distribusi kecepatan udara dalam pengering pada simulasi ditunjukkan pada kontur kecepatan (Gambar 10) dan vektor kecepatan (Gambar 11). Nilai hasil simulasi selengkapnya dapat dilihat pada Lampiran 1.
17
Gambar 10 Kontur kecepatan udara pada simulasi 1
Gambar 11 Vektor distribusi kecepatan udara pada simulasi 1 Udara dengan suhu 40oC dihembuskan melalui lubang inlet menuju plenum kemudian melewati lantai berlubang disebarkan ke seluruh bagian bak tumpukan. Kecepatan aliran udara tertinggi terdapat di bagian lubang inlet, yaitu pada kisaran 13.2 m/s yang ditunjukkan oleh warna merah. Aliran udara ini mengalami turbulensi di dalam plenum sebelum melewati lantai pengering yang berlubang. Kecepatan udara di dalam plenum mengalami penurunan yaitu pada kisaran 10.57 m/s yang diambil pada titik tengah plenum. Saat melewati lantai berpori kecepatan udara mengalami penurunan berkisar pada 0.669 m/s pada lapisan bawah, 0.687 m/s pada lapisan tengah dan 0.185 m/s pada lapisan di atas tumpukan.
18
Gambar 12 Vektor distribusi kecepatan udara sisi bawah plenum pada simulasi 1 Berdasarkan Gambar 12, vektor kecepatan aliran udara pada plenum terjadi turbulensi yang tinggi saat aliran masuk dari inlet membentur dinding plenum sisi belakang. Arah aliran kemudian berputar ke arah sisi kanan dan kiri plenum mengikuti bentuk penampang plenum berbalik ke arah datangnya aliran. Kecepatan aliran udara setelah melewati lantai berlubang mengalami penurunan. Hal itu ditandai dengan perubahan warna yang menandai penurunan kecepatan. Arah vektor saat melewati lantai berlubang bergerak miring ke samping. Arah vektor aliran kecepatan udara pada tumpukan gabah tampak lurus ke atas. Arah vektor aliran kecepatan udara di atas tumpukan gabah sedikit tidak beraturan. Hal tersebut terlihat dari adanya aliran yang mengarah ke atas, ke samping dan sedikit ke arah bawah. Pada posisi outlet sisi belakang terlihat arah vektor aliran udara yang lebih tinggi dibanding sisi lainnya. Pada posisi outlet tengah terjadi cekungan karena arah vektor pada lokasi ini terdistribusi ke arah belakang dan sedikit ke arah bawah. Pada posisi tengah seperti ada tekanan aliran dari atas yang bertemu dengan aliran dari bawah sehingga aliran terdesak ke sisi belakang. Hal tersebut dapat dilihat pada Gambar 13.
Gambar 13 Vektor kecepatan udara arah sumbu-y posisi outlet simulasi 1
19 2. Aliran kecepatan udara simulasi 2 Distribusi kecepatan udara dalam pengering pada simulasi ditunjukkan pada kontur kecepatan (Gambar 14) dan vektor kecepatan (Gambar 15). Nilai hasil simulasi selengkapnya dapat dilihat pada Lampiran 2.
Gambar 14 Kontur kecepatan udara pada simulasi 2
Gambar 15 Vektor distribusi kecepatan udara pada simulasi 2 Udara dengan suhu 40oC dihembuskan melalui lubang inlet menuju plenum kemudian melewati lantai berlubang disebarkan ke seluruh bagian bak tumpukan. Kecepatan aliran udara tertinggi terdapat di bagian lubang inlet, yaitu pada kisaran 13.5 m/s yang ditunjukkan oleh warna merah. Aliran udara ini mengalami turbulensi di dalam plenum sebelum melewati lantai pengering yang berlubang. Kecepatan udara di dalam plenum mengalami penurunan yaitu pada kisaran 10 m/s yang diambil pada titik tengah plenum. Saat melewati lantai berpori kecepatan udara mengalami penurunan berkisar pada 0.610 m/s
20 pada lapisan bawah, 0.685 m/s pada lapisan tengah dan 0.243 m/s pada lapisan di atas tumpukan.
Gambar 16 Vektor distribusi kecepatan udara sisi bawah plenum pada simulasi 2 Berdasarkan Gambar 16, vektor kecepatan aliran udara pada plenum terjadi turbulensi yang tinggi saat aliran masuk dari inlet membentur dinding plenum sisi belakang. Arah aliran kemudian berputar ke arah sisi kanan dan kiri plenum mengikuti bentuk plenum dimana kecepatan aliran pada dinding plenum yang paling tinggi. Kecepatan aliran udara setelah melewati lantai berlubang mengalami penurunan. Hal itu ditandai dengan perubahan warna yang menandai penurunan kecepatan. Arah vektor saat melewati lantai berlubang bergerak miring ke samping cenderung tidak beraturan. Arah vektor aliran kecepatan udara pada tumpukan gabah tampak lurus ke atas. Arah vektor aliran kecepatan udara di atas tumpukan gabah sedikit tidak beraturan. Hal tersebut terlihat dari adanya aliran yang mengarah ke atas, ke samping dan sedikit ke arah bawah. Pada posisi outlet belakang terlihat arah vektor aliran udara yang lebih tinggi dibanding sisi lainnya. Pada posisi outlet tengah terjadi cekungan karena arah vektor pada lokasi ini terdistribusi ke arah belakang dan sedikit ke arah bawah. Pada posisi tengah seperti ada tekanan aliran dari atas yang bertemu dengan aliran dari bawah sehingga aliran terdesak ke sisi belakang. Hal tersebut dapat dilihat pada Gambar 17.
Gambar 17 Vektor kecepatan udara arah sumbu-y posisi outlet simulasi 2
21 3. Aliran kecepatan udara simulasi 3 Distribusi kecepatan udara dalam pengering pada simulasi ditunjukkan pada kontur kecepatan (Gambar 18) dan vektor kecepatan (Gambar 19). Nilai hasil simulasi selengkapnya dapat dilihat pada Lampiran 3.
Gambar 18 Kontur kecepatan udara pada simulasi 3
Gambar 19 Vektor distribusi kecepatan aliran udara pada simulasi 3 Udara dengan suhu 40oC dihembuskan melalui lubang inlet menuju plenum kemudian melewati lantai berlubang disebarkan ke seluruh bagian bak tumpukan. Kecepatan aliran udara tertinggi terdapat di bagian lubang inlet, yaitu pada kisaran 13.2 m/s yang ditunjukkan oleh warna merah. Aliran udara ini mengalami turbulensi di dalam plenum sebelum melewati lantai pengering yang berlubang. Kecepatan udara di dalam plenum mengalami penurunan yaitu pada kisaran 6.78 m/s yang diambil pada titik tengah plenum. Saat melewati lantai berpori kecepatan udara mengalami penurunan berkisar pada 0.706 m/s
22 pada lapisan bawah, 0.673 m/s pada lapisan tengah dan 0.183 m/s pada lapisan di atas tumpukan.
Gambar 20 Vektor distribusi kecepatan udara sisi bawah plenum pada simulasi 3 Berdasarkan Gambar 20, vektor kecepatan aliran udara pada plenum terjadi turbulensi yang tinggi saat aliran masuk dari inlet. Namun aliran tidak sampai membentur sisi ujung plenum. Arah aliran kemudian menyebar ke arah sisi dinding kanan dan kiri plenum mengikuti bentuk plenum dimana kecepatan aliran pada dinding plenum yang paling tinggi. Kecepatan aliran udara setelah melewati lantai berlubang mengalami penurunan. Kecepatan aliran pada posisi dekat dengan inlet lebih tinggi dibanding di ujung. Arah vektor saat melewati lantai berlubang bergerak miring ke samping cenderung tidak beraturan. Arah vektor aliran kecepatan udara pada tumpukan gabah tampak lurus ke atas. Arah vektor aliran kecepatan udara di atas tumpukan gabah sedikit tidak beraturan. Hal tersebut terlihat dari adanya aliran yang mengarah ke atas, ke samping dan sedikit ke arah bawah. Pada posisi outlet belakang terlihat arah vektor aliran udara yang lebih tinggi dibanding sisi lainnya. Pada posisi outlet tengah terjadi cekungan karena arah vektor pada lokasi ini terdistribusi ke arah belakang dan sedikit ke arah bawah. Pada posisi tengah seperti ada tekanan aliran dari atas yang bertemu dengan aliran dari bawah sehingga aliran terdesak ke sisi belakang. Hal tersebut dapat dilihat pada Gambar 21.
Gambar 21 Vektor kecepatan udara arah sumbu-y posisi outlet simulasi 3
23 Pemilihan Desain Tingkat keragaman kecepatan aliran udara dan penurunan tekanan pada ketiga simulasi didapatkan dari nilai rata-rata hasil simulasi. Nilai keragaman dan masing-masing standar deviasi tersebut digunakan untuk mengevaluasi pengaruh bentuk geometri terhadap keseragaman kecepatan aliran udara dan penurunan tekanan yang terjadi pada kondisi operasi yang sama. Berdasarkan nilai rata-rata dan standar deviasi hasil simulasi maka dipilih desain yang akan dipabrikasi. Kontur kecepatan aliran udara pada tiga lapisan ketebalan secara vertikal dan horisontal dari hasil simulasi untuk simulasi 1, 2 dan 3 dapat dilihat pada Gambar 22, 23 dan 24. Pada Simulasi 1 kecepatan aliran udara terkecil terdapat pada ketebalan 20 cm dengan nilai 0.678 m/s, sementara kecepatan terbesar terdapat pada ketebalan 65 cm dengan nilai 0.688 m/s, sedangkan rata-rata kecepatan aliran udara adalah sebesar 0.565 m/s. Standar deviasi rata-rata kecepatan aliran udara pada simulasi 1 sebesar 0.0139 m/s.
Gambar 22 Profil kecepatan aliran udara simulasi 1 Untuk simulasi 2, kecepatan aliran udara terkecil berada pada ketebalan 20 cm sebesar 0.661 m/s, sedangkan kecepatan terbesar berada pada ketebalan 65 cm sebesar 0.686 m/s. Kecepatan aliran udara rata-rata pada simulasi 2 adalah sebesar 0.560 m/s. Standar deviasi rata-rata kecepatan aliran udara pada simulasi 2 sebesar 0.0131 m/s.
24
Gambar 23 Profil kecepatan aliran udara simulasi 2 Untuk simulasi 3, kecepatan aliran udara terkecil berada pada ketebalan 20 cm sebesar 0.669 m/s, sedangkan kecepatan terbesar berada pada ketebalan 45 cm sebesar 0.674 m/s. Kecepatan aliran udara rata-rata pada simulasi 3 adalah sebesar 0.556 m/s. Standar deviasi rata-rata kecepatan aliran udara pada simulasi 3 sebesar 0.016 m/s.
Gambar 24 Profil kecepatan aliran udara pada 3 ketinggian simulasi 3 Perbandingan standar deviasi kecepatan aliran udara dalam lapisan tumpukan pada tiga simulasi disajikan secara grafis pada Gambar 25, sementara penyajian datanya dapat dilihat pada Lampiran 4.
Stdandar deviasi kecepatan udara (m/s)
25 0.04 0.035 0.03 0.025 0.02 0.015 0.01 0.005
0 20 Simulasi 1
45 65 Ketinggian (cm) Simulasi 2
70
Simulasi 3
Gambar 25 Perbandingan standar deviasi kecepatan udara dalam lapisan tumpukan pada simulasi CFD
Standar deviasi rata-rata
0.0160 0.0150 0.0140
0.0130 0.0120 0.0110 simulasi 1
simulasi 2
simulasi 3
Simulasi CFD Gambar 26 Standar deviasi rata-rata pada ketiga simulasi CFD Pada Gambar 26 disajikan perbandingan nilai standar deviasi kecepatan udara pada ketiga simulasi. Berdasarkan Gambar 26, simulasi 1 memiliki nilai standar deviasi rata-rata sebesar 0.0139, simulasi 2 sebesar 0.031 dan simulasi 3 sebesar 0.0155. Simulasi 2 memiliki standar deviasi rata-rata yang paling kecil. Hal tersebut menunjukkan bahwa simulasi 2 memiliki tingkat keseragaman kecepatan udara yang lebih baik dibandingkan simulasi 1 dan 3. Validasi Pressure Drop Validasi pressure drop dilakukan dua tahap yaitu pada awal penentuan nilai porositas model tumpukan gabah dan pada pengukuran saat uji kinerja pengering gabah. Validasi yang dilakukan diawal yakni bertolok ukur pada nilai korelasi
26 kurva pressure drop antara simulasi dengan persamaan airflow resistance (ASAE D272.3 Mar 1996). Persamaan airflow resistance yang digunakan ialah :
=
................................................................................... (8) o e( ) nilai a dan b merupakan konstanta dalam airflow resistance. Nilai a ialah 2.57x104 Pa.s2/m3 dan nilai b ialah 13.2 m2.s/m3 (Shedd 1953). Validasi data pressure drop CFD dan data pressure drop berdasarkan persamaan airflow resistance dilihat pada Gambar 27. Berdasarkan trial and error yang dilakukan dalam simulasi ditetapkan nilai porositas pada model simulasi tumpukan gabah yaitu 30%. Adapun nilai porositas dari padi dan beras berdasarkan data penelitian ASAE Wratten et al (1969) dalam Champ (1996) ialah antara 46 sampai 60% tergantung pada varietas dan kadar air. Nilai porositas yang diperoleh pada pendekatan CFD cukup berbeda dengan data literatur. Hal ini dapat disebabkan oleh varietas padi dan kadar air yang berbeda serta pendefinisian bentuk porous pada cell zone gabah yang tidak sama persis pada bentuk porous tumpukan gabah yang sebenarnya. Beberapa peneliti menemukan bahwa laju aliran udara, permukaan dan karakteristik bentuk dari bijian, ukuran dan konfigurasi kekosongan tumpukan, kadar air, variasi ukuran partikel, ketebalan bahan, tingkat pengemasan, dan jumlah dari material asing mempengaruhi pressure drop dari laju aliran udara yang melewati biji-bijian (Champ 1995). Pada ketebalan bahan dan kadar air yang sama, terlihat bahwa laju aliran udara mempengaruhi pressure drop yang terjadi pada tumpukan gabah. Laju aliran yang tinggi menunjukkan pressure drop yang terjadi juga tinggi. Sedangkan semakin rendah laju aliran udara menunjukkan pressure drop yang terjadi juga semakin rendah. Hal tersebut dapat dilihat pada Tabel 1. Tabel 1 Data pressure drop CFD, pengukuran dan kurva Shedd airflow (m3/s m2)
pressure drop CFD (Pa/m)
pressure drop ukur (Pa/m)
pressure drop kurva Shedd (Pa/m)
0.21 0.2 0.1 0.05 0.03
899.2 782.1 292.4 131.8 72.4
882.9 804.4 333.5 157.0 58.9
830.3 795.7 305.4 126.8 69.3
pressure drop Shedd, Pa/m
27
1000.0 y = 0.9544x + 9.776 R² = 0.9936
800.0 600.0 400.0 200.0 0.0 0.0
200.0 400.0 600.0 800.0 1000.0 pressure drop CFD, Pa/m
Gambar 27 Validasi pressure drop CFD dengan kurva Shedd (1953)
Pressure drop plenum ukur (Pa/m)
Berdasarkan grafik yang ditunjukkan oleh Gambar 27, nilai korelasi yang dimiliki oleh validasi tersebut ialah sebesar 0.9936 dan persamaannya mendekati 1. Sehingga asumsi pada simulasi CFD telah cukup baik digunakan sebagai prediksi kondisi yang sebenarnya. Nilai pressure drop hasil simulasi CFD selanjutnya divalidasi dengan pressure drop yang diperoleh dari hasil pengukuran. Validasi pressure drop hasil simulasi CFD dengan hasil pengukuran dilihat pada Gambar 28. 1000.0
y = 0.9855x + 18.08 R² = 0.9957
800.0 600.0 400.0 200.0 0.0 0.0
200.0 400.0 600.0 800.0 1000.0 Pressure drop CFD (Pa/m)
Gambar 28 Validasi pressure drop CFD dengan pressure drop pengukuran Berdasarkan grafik yang ditunjukkan oleh Gambar 28, nilai korelasi yang dimiliki oleh validasi tersebut ialah sebesar 0.9957 dan persamaannya mendekati 1. Sehingga hasil simulasi CFD memiliki nilai yang mendekati hasil pengukuran. Profil kontur distribusi tekanan pada simulasi secara vertikal disajikan pada Gambar 29. Pada Gambar 29 terlihat bahwa tekanan udara yang paling tinggi berada pada plenum. Tekanan terlihat semakin menurun setelah melewati plat berlubang dan terus menurun sepanjang tumpukan.
28
Gambar 29 Kontur tekanan pada simulasi 2
Validasi Kecepatan Aliran Udara Hasil validasi kecepatan aliran udara CFD berupa kontur distribusi kecepatan aliran udara, menunjukkan besar dan arah aliran udara sebagai media pengering yang digunakan dalam pengering gabah. Validasi model dilakukan dengan membandingkan data pengukuran dan data simulasi CFD pada 13 titik pengukuran. Nilai hasil simulasi dan hasil pengukuran kecepatan udara pada bidang-xz dengan ketinggian-y dapat dilihat pada Lampiran 5.
Gambar 30 Validasi kecepatan aliran udara hasil simulasi CFD dengan hasil pengukuran
29 Berdasarkan Gambar 30 terlihat bahwa kecepatan aliran udara hasil simulasi hampir mengikuti data pengukuran, walaupun secara garis besar memiliki nilai yang berbeda cukup signifikan. Perbedaan ini salah satunya disebabkan oleh penentuan asumsi media berpori tumpukan pada simulasi CFD. Pendefinisian bentuk pori oleh CFD berbeda dengan kondisi yang sebenarnya sehingga menyebabkan perbedaan hasil yang diperoleh. Gambar 30 menunjukkan penyajian grafis hasil validasi kecepatan udara dengan membandingkan kecepatan udara hasil pengukuran dan hasil simulasi pada tiga ketebalan dalam tumpukan. Berdasarkan Gambar 33 terlihat bahwa kecepatan udara pada ketinggian 70 cm atau tepat di atas tumpukan gabah sangat rendah dibanding kedua lapisan lainnya. Pada lapisan 20 dan 45 cm atau di dalam tumpukan gabah, laju aliran udara yang tinggi dikarenakan adanya pori dalam tumpukan gabah. Pori tumpukan gabah tersebut berupa cell zone gabah yang menggunakan porous formulation physical velocity sehingga terjadi peningkatan kecepatan aliran udara pada seluruh wilayah berpori. Pada lapisan 70 cm atau di luar cell zone gabah, kecepatan lebih rendah karena secara default terpilih opsi superficial velocity. Pada superficial velocity, kecepatan aliran udara tidak mengalami peningkatan seperti kecepatan udara yang terjadi di dalam pori. Perbedaan antara data kecepatan udara hasil pengukuran dan kecepatan udara hasil simulasi CFD diberikan dalam nilai standar deviasi rata-rata sebesar 0.051 m/s. Penyimpangan dari hasil validasi dinyatakan dalam rata-rata error sebesar 14.9%.
Hasil Uji Kinerja Pengeringan Gabah
Penurunan Kadar Air (%bb)
1. Perubahan Kadar Air Gabah Data pengukuran perubahan kadar air bijian dalam pengering selama 16 jam pengeringan disajikan pada Lampiran 6. Secara grafis dapat dilihat pada Gambar 31. Hasil pengukuran kadar air pada pengering di lapangan menunjukkan terjadinya penurunan yang stabil dengan rata-rata penurunan 2%b.b. tiap dua jam pada awal pengeringan. Kadar air keseimbangan rata-rata tumpukan gabah 14% tercapai setelah pengeringan selama 12 jam. 30 25 20 15 10 5 0 0
5 10 Waktu pengeringan (jam) Perubahan Kadar Air
15
Kadar Air 14%
Gambar 31 Penurunan kadar air rata-rata dalam tumpukan gabah
20
30
Penurunan Kadar Air (%bb)
Penurunan kadar air tiap lapisan disajikan secara grafis pada Gambar 32. Terlihat bahwa penurunan kadar air pada lapisan bawah terjadi paling cepat dibandingkan lapisan tengah dan atas. Kondisi ini dikarenakan lapisan bawah berada pada posisi paling dekat udara panas dihembuskan. Pada lapisan tengah penurunan kadar air terjadi diantara lapisan bawah dan atas. Lain halnya penurunan kadar air pada lapisan atas yang cenderung lambat pada awal pengeringan lalu menurun saat kadar air lapisan bawah telah melewati garis keseimbangan kadar air 14%. Penurunan kadar air pada lapisan atas tumpukan terutama dipengaruhi oleh kelembaban udara yang meningkat. Peningkatan kelembaban udara ini terjadi karena udara dari lapisan dibawahnya membawa uap air hasil pengeringan gabah. 30 25 20 15 10 5 0 0 lapisan bawah
5
10 15 Waktu pengeringan (jam) lapisan tengah
lapisan atas
20 K.A. 14%
Gambar 32 Penurunan kadar air rata-rata dalam tiga lapisan tumpukan 2. Suhu udara pengukuran dalam tumpukan gabah Data pengukuran suhu pada pengering selama 16 jam pengeringan disajikan pada Lampiran 8. Secara grafis dapat dilihat pada Gambar 33. Hasil pengukuran suhu pada tumpukan dengan tiga lapisan pengukuran memperlihatkan bahwa terjadi kenaikan suhu pada tumpukan seiring lamanya waktu pengeringan. Hal ini disebabkan oleh penurunan kelembaban di dalam tumpukan gabah. Lapisan bawah memiliki suhu udara yang paling tinggi karena posisinya paling dekat dengan plenum. Lapisan atas memiliki suhu udara yang paling rendah akibat kelembaban yang meningkat dari lapisan dibawahnya. Sedangkan lapisan tengah memiliki suhu udara diantara kedua lapisan. Hall (1957) mengungkapkan bahwa ketika udara panas bergerak pada produk, panas ditransfer ke produk dan terjadi penguapan air sehingga meningkatkan kelembaban relatif udara. Proses ini dikenal sebagai panas simultan dan transfer massa. Efek pengeringan dari udara ini terkait dengan rasio kadar air udara terhadap kadar air jenuh pada temperatur yang sama.
Suhu (oC)
31 45 40 35 30 25 20 15 10 05 00 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Waktu pengeringan (Jam) Bawah Tengah Atas Gambar 33 Suhu udara dalam tiga lapisan tumpukan gabah terhadap waktu 120 RH Udara (%)
100 80 60 40
20 0 0
1
2
3
4
5 6 7 8 9 10 11 12 13 14 15 16 Waktu pengeringan (Jam) RH Udara Lingkungan RH Udara Plenum RH Udara Keluar Tumpukan
Gambar 34 RH udara pada tiga lokasi pengukuran terhadap waktu Data perubahan RH dapat dilihat pada Lampiran 8. Secara grafis dapat dilihat pada Gambar 34. Udara lingkungan dengan RH tinggi dihembuskan ke dalam plenum setelah terlebih dahulu dipanaskan sehingga nilai RH-nya menurun. Pada proses tersebut udara mengalami proses pengeringan. Udara plenum dengan RH yang rendah tersebut dilewatkan di dalam tumpukan gabah. Di dalam tumpukan gabah, udara menerima uap air yang dikandung gabah sehingga RH udara mengalami peningkatan setelah keluar dari tumpukan gabah. Seiring penurunan kadar air gabah menuju periode konstan, terjadi penurunan RH udara yang keluar dari tumpukan gabah. Hal tersebut menunjukkan bahwa uap air dari gabah yang diterima udara pengering mengalami penurunan dan menuju konstan. Kondisi tersebut terjadi setelah pengeringan berjalan selama 11 jam.
32
SIMPULAN DAN SARAN Simpulan Berdasarkan hasil penelitian ini dapat disimpulkan bahwa simulasi pola aliran udara pada simulasi geometri penampang lingkaran memiliki hasil yang lebih seragam. Nilai error rata-rata validasi kecepatan udara ialah 14.9%. Hal tersebut dikarenakan pendekatan model pori oleh simulasi tidak sama dengan kondisi sebenarnya. Koefisien korelasi pada pressure drop menunjukkan hasil yang cukup baik yaitu 0.9957. Pengujian kinerja pengering gabah sebanyak 128.5 kg pada kadar air awal 28.3% b.b. dilakukan selama 12 jam hingga mencapai kadar air keseimbangan rata-rata 14% b.b.
Saran Saran yang dapat diberikan adalah sebaiknya dilakukan penelitian lanjutan mengenai alternatif bentuk lantai pengering dan pendekatan asumsi tumpukan mengandung kadar air sehingga dapat disimulasikan perubahan kadar air.
DAFTAR PUSTAKA Champ B.R., Hightley E., Johnson G.I., ed. 1995. Grain Drying in Asia. Proceedings of an International conference held at the FAO Regional Office for Asia and the Pacific, Bangkok, Thailand, 17-20 October 1995. ACIAR Proceedings No. 71, 410p. Hall C.W. 1957. Drying Farm Crops. Michigan: Agricultural Consulting Associates, Inc. Hall C.W. 1970. Handling and Storage Food Grain in Tropical and Subtropical Areas. Roma: FAO. Hall C.W. 1980. Drying And Storage Of Agricultural Crops. AVI Publishing Company, INC. Westport, Connecticut. Margana Cahyawan C.E., Sukmawaty. 2010. Karakteristik dan Simulasi Mesin Pengering Gabah Tipe Vertikal Kontinyu dengan Aliran Udara Panas Berlawanan. Prosiding Seminar Nasional PERTETA 2010. Purwakarta. Nurba D. 2008. Analisis Distribusi Suhu, Aliran Udara, RH, dan Kadar Air dalam In-Store Dryer (ISD) untuk Biji Jagung. Thesis. Program Pascasarjana. Institut Pertanian Bogor. Shedd C.K. 1953. Resistance of Grains and Seeds to Air Flow. Agricultural Engineering 34(9):616-619. Tuakia F. 2008. Dasar-Dasar CFD Menggunakan Fluent. Bandung(ID) : Bandung Informatika.
33 Lampiran 1 Hasil simulasi 1 CFD Poin 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 28 29 30 31 32 33 34 35 36 37 38 39 40
Koordinat X (cm) Y (cm) Z (cm) 0 20 0 15 20 0 30 20 0 -15 20 0 -30 20 0 -15 20 15 0 20 15 15 20 15 -30 20 30 0 20 30 30 20 30 -15 20 -15 0 20 -15 15 20 -15 -30 20 -30 0 20 -30 30 20 -30 0 45 0 15 45 0 30 45 0 -15 45 0 -30 45 0 -15 45 15 0 45 15 15 45 15 -30 45 30 0 45 30 30 45 30 -15 45 -15 0 45 -15 15 45 -15 -30 45 -30 0 45 -30 30 45 -30 0 65 0 15 65 0 30 65 0 -15 65 0 -30 65 0 -15 65 15
v (m/s)
P (Pa)
0.669 0.680 0.665 0.638 0.664 0.618 0.701 0.647 0.701 0.675 0.690 0.670 0.716 0.668 0.704 0.759 0.661 0.687 0.687 0.687 0.687 0.687 0.686 0.686 0.686 0.685 0.685 0.685 0.689 0.689 0.689 0.690 0.690 0.690 0.690 0.690 0.687 0.689 0.686 0.688
448.227 449.493 449.346 448.817 449.955 449.523 447.400 447.703 450.115 446.723 450.180 455.624 454.864 456.603 460.423 458.891 459.634 215.744 217.343 216.549 216.798 218.475 217.752 219.44 215.941 215.779 215.338 215.975 218.062 217.333 215.263 218.827 217.83 218.412 29.1279 29.4777 29.7088 29.4306 30.1509 29.783
Rata-rata & St. Deviasi
Rata V St. Dev V
0.678 0.032
Rata P St. Dev P
451.972 4.620
Rata V St. Dev V
0.687 0.002
Rata P St. Dev P
217.109 1.302
34 Poin
Koordinat X (cm) Y (cm) Z (cm) 0 65 15 15 65 15 -30 65 30 0 65 30 30 65 30 -15 65 -15 0 65 -15 15 65 -15 -30 65 -30 0 65 -30 30 65 -30 0 70 0 15 70 0 30 70 0 -15 70 0 -30 70 0 -15 70 15 0 70 15 15 70 15 -30 70 30 0 70 30 30 70 30 -15 70 -15 0 70 -15 15 70 -15 -30 70 -30 0 70 -30 30 70 -30
41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 Rata-rata Rata standar deviasi
v (m/s)
P (Pa)
0.688 0.687 0.686 0.685 0.686 0.690 0.693 0.692 0.686 0.685 0.686 0.186 0.187 0.206 0.186 0.209 0.188 0.184 0.189 0.215 0.214 0.216 0.216 0.211 0.239 0.233 0.235 0.230 0.565 0.0139
28.023 28.2133 30.5131 29.9495 30.3844 29.7668 30.0242 30.6294 31.1284 31.5082 30.1662 0.716 0.609 0.403 0.633 0.423 0.486 0.591 0.481 0.224 0.311 0.229 0.479 0.539 0.653 0.845 1.159 0.874 174.883 1.7645
Rata-rata & St. Deviasi
Rata V St. Dev V
0.688 0.002
Rata P St. Dev P
29.881 0.895
Rata V St. Dev V
0.208 0.019
Rata P St. Dev P
0.568 0.240
35 Lampiran 2 Hasil simulasi 2 CFD Poin 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 28 29 30 31 32 33 34 35 36 37 38 39 40
Koordinat X (cm) Y (cm) Z (cm) 0 20 0 15 20 0 30 20 0 -15 20 0 -30 20 0 -15 20 15 0 20 15 15 20 15 -25 20 25 0 20 30 25 20 25 -15 20 -15 0 20 -15 15 20 -15 -25 20 -25 0 20 -30 25 20 -25 0 45 0 15 45 0 30 45 0 -15 45 0 -30 45 0 -15 45 15 0 45 15 15 45 15 -25 45 25 0 45 30 25 45 25 -15 45 -15 0 45 -15 15 45 -15 -25 45 -25 0 45 -30 25 45 -25 0 65 0 15 65 0 30 65 0 -15 65 0 -30 65 0 -15 65 15
v (m/s)
P (Pa)
0.598 0.631 0.659 0.627 0.660 0.640 0.642 0.617 0.662 0.696 0.662 0.661 0.701 0.672 0.698 0.722 0.698 0.685 0.685 0.685 0.685 0.685 0.684 0.684 0.684 0.684 0.683 0.684 0.687 0.687 0.687 0.688 0.689 0.688 0.688 0.685 0.685 0.686 0.686 0.685
429.228 430.370 432.287 430.275 431.392 428.921 428.173 430.258 429.887 430.737 428.816 434.283 437.170 435.327 436.628 441.655 436.513 214.862 213.666 214.777 214.794 215.196 215.444 214.379 214.968 213.240 213.952 213.750 215.497 215.123 215.772 214.818 216.991 214.517 28.617 28.611 29.212 27.471 27.774 28.512
Rata-rata & St. Deviasi
Rata v St. Dev v
0.661 0.034
Rata P St. Dev P
432.466 3.803
Rata v St. Dev v
0.686 0.002
Rata P St. Dev P
214.809 0.892
36 Poin
Koordinat X (cm) Y (cm) Z (cm) 0 65 15 15 65 15 -25 65 25 0 65 30 25 65 25 -15 65 -15 0 65 -15 15 65 -15 -25 65 -25 0 65 -30 25 65 -25 0 70 0 15 70 0 30 70 0 -15 70 0 -30 70 0 -15 70 15 0 70 15 15 70 15 -25 70 25 0 70 30 25 70 25 -15 70 -15 0 70 -15 15 70 -15 -25 70 -25 0 70 -30 25 70 -25
41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 Rata-rata Rata Standar Deviasi
v (m/s)
P (Pa)
0.685 0.685 0.684 0.684 0.684 0.688 0.689 0.687 0.684 0.685 0.685 0.239 0.185 0.201 0.186 0.200 0.186 0.208 0.199 0.205 0.201 0.205 0.203 0.201 0.205 0.217 0.235 0.223 0.560 0.0131
26.272 27.285 27.137 27.521 27.525 28.228 27.485 27.294 28.853 28.101 29.117 0.867 0.472 0.274 0.480 0.276 0.332 0.624 0.445 0.180 0.182 0.193 0.392 0.457 0.361 0.567 0.523 0.615 168.911 1.4200
Rata-rata & St. Deviasi
Rata V St. Dev V
0.686 0.002
Rata P St. Dev P
27.942 0.801
Rata v St. Dev v
0.206 0.015
Rata P St. Dev P
0.426 0.184
37 Lampiran 3 Hasil simulasi 3 CFD Koordinat
Poin 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 28 29 30 31 32 33 34 35 36 37 38 39
X (cm) 0 -20 -10 10 20 -15 0 15 0 -20 0 20 -15 0 15 0 -20 0 20 0 -20 -10 10 20 -15 0 15 0 -20 0 20 -15 0 15 0 -20 0 20 0
Y (cm) 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 45 65
Z (cm) 0 0 0 0 0 15 15 15 30 45 45 45 -15 -15 -15 -30 -45 -45 -45 0 0 0 0 0 15 15 15 30 45 45 45 -15 -15 -15 -30 -45 -45 -45 0
v (m/s)
P (Pa)
0.707 0.665 0.644 0.680 0.665 0.631 0.634 0.606 0.623 0.679 0.677 0.675 0.642 0.699 0.698 0.699 0.694 0.695 0.694 0.674 0.673 0.673 0.673 0.673 0.668 0.668 0.668 0.665 0.664 0.664 0.664 0.679 0.679 0.679 0.683 0.684 0.684 0.684 0.675
438.358 436.528 435.479 436.746 436.202 427.845 426.523 426.865 422.729 428.046 426.384 428.758 444.228 443.982 444.228 446.197 448.707 447.085 449.288 208.534 207.792 208.250 208.476 208.230 207.392 206.788 208.261 205.151 204.098 205.811 204.920 209.150 210.958 209.878 212.734 213.472 214.058 213.663 27.642
Rata-rata & St. Deviasi
Rata v St. Dev v
0.669 0.030
Rata P St. Dev P
436.536 8.785
Rata v St. Dev v
0.674 0.007
Rata P St. Dev P
208.822 3.001
38 Poin
Koordinat X (cm)
40 -20 41 -10 42 10 43 20 44 -15 45 0 46 15 47 0 48 -20 49 0 50 20 51 -15 52 0 53 15 54 0 55 -20 56 0 57 20 58 0 59 -20 60 -10 61 10 62 20 63 -15 64 0 65 15 66 0 67 -20 68 0 69 20 70 -15 71 0 72 15 73 0 74 -20 75 0 76 20 Rata-rata Rata standar deviasi
Y (cm)
Z (cm)
65 65 65 65 65 65 65 65 65 65 65 65 65 65 65 65 65 65 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70 70
0 0 0 0 15 15 15 30 45 45 45 -15 -15 -15 -30 -45 -45 -45 0 0 0 0 0 15 15 15 30 45 45 45 -15 -15 -15 -30 -45 -45 -45
v (m/s)
P (Pa)
0.673 0.673 0.674 0.673 0.670 0.669 0.670 0.668 0.666 0.666 0.666 0.678 0.677 0.678 0.683 0.680 0.677 0.680 0.183 0.201 0.222 0.188 0.200 0.196 0.182 0.195 0.199 0.208 0.202 0.209 0.259 0.185 0.196 0.215 0.228 0.224 0.226 0.556 0.016
25.995 26.845 26.064 27.466 26.161 25.397 27.196 25.138 26.047 26.580 26.988 27.046 27.186 26.843 26.988 28.584 27.744 29.193 0.640 0.401 0.833 0.553 0.428 0.370 0.558 0.387 0.533 0.186 0.208 0.187 1.059 0.653 0.587 0.591 1.218 1.559 1.214 168.224 3.292
Rata-rata & St. Deviasi
Rata v St. Dev v
0.673 0.005
Rata P St. Dev P
26.900 1.004
Rata v St. Dev v
0.206 0.019
Rata P St. Dev P
0.640 0.379
39 Lampiran 4 Perbandingan keragaman kecepatan aliran udara dan tekanan dalam pengering pada ketiga simulasi CFD Keragaman kecepatan aliran udara Ketinggian (cm) Rata-rata (m/dtk) Standar Deviasi (m/s) Rata-rata (m/dtk) Standar Deviasi (m/s) Rata-rata (m/dtk) Standar Deviasi (m/s)
20 45 Simulasi 1 0.6781 0.6873 0.0324 0.0018 Simulasi 2 0.6615 0.6855 0.0337 0.0018 Simulasi 3 0.6688 0.6737 0.0303 0.0072
65
70 Rata-rata
0.6878 0.0023
0.2084 0.0190
0.5654 0.0139
0.6856 0.0015
0.2058 0.0154
0.5596 0.0131
0.6734 0.0053
0.2063 0.0193
0.5555 0.0155
Keragaman tekanan udara Ketinggian (cm) Rata-rata (Pa) St. Deviasi (Pa) Rata-rata (Pa) St. Deviasi (Pa) Rata-rata (Pa) St. Deviasi (Pa)
20
45 Simulasi 1 451.9718 217.1095 4.6204 1.3022 Simulasi 2 432.4659 214.8086 3.8034 0.8922 Simulasi 3 436.5357 208.8219 8.7855 3.0008
65
70
Rata-rata
29.8815 0.5679 174.8827 0.8949 0.2405 1.7645 27.9419 0.4259 168.9106 0.8007 0.1837 1.4200 26.9001 0.6402 168.2245 1.0040 0.3793 3.2924
40 Lampiran 5 Data validasi kecepatan aliran udara hasil pengukuran dan simulasi CFD serta nilai error dan standar deviasi Koordinat Poin
Simbol
1 B1 2 B2 3 B3 4 B4 5 T1 6 T2 7 T3 8 T4 9 A1 10 A2 11 A3 12 A4 13 A5 Maksimum Minimum Rata-rata Korelasi (R)
x (cm)
y (cm)
z (cm)
0 -25 0 25 0 -25 0 25 0 -25 0 25 0
20 20 20 20 45 45 45 45 70 70 70 70 70
25 0 -25 0 25 0 -25 0 25 0 -25 0 0
Kecepatan udara (13.2 m/s) vv-ukur SD Error CFD 0.63 0.70 0.014 11.1% 0.82 0.65 0.021 20.7% 0.62 0.73 0.014 17.7% 0.69 0.63 0.007 8.7% 0.61 0.68 0.007 11.5% 0.56 0.68 0.014 21.4% 0.61 0.68 0.042 11.5% 0.53 0.68 0.085 28.3% 0.18 0.20 0.007 11.1% 0.21 0.20 0.049 4.8% 0.21 0.23 0.000 9.5% 0.26 0.25 0.007 3.8% 0.18 0.24 0.000 33.3% 0.82 0.73 0.120 33.3% 0.18 0.20 0.007 3.8% 0.47 0.50 0.051 14.9% 0.87
41 Lampiran 6 Hasil pengukuran kadar air bijian selama 16 jam
Jam ke
Lapisan
Poin
0(Awal)
Bawah
2
Tengah
Atas
Bawah
4
Tengah
Atas
Bawah 6 Tengah
B-De B-Be B-Ka B-Ki B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki B-Te T-De T-Be
%b.k
39.53 24.04 24.14 24.51 23.94 25.34 37.22 35.85 37.41 39.00 37.76 38.38 37.26 36.64 37.73 37.35 20.18 19.69 20.47 20.17 21.18 29.65 29.30 34.73 31.75 33.58 34.68 36.86 35.61 37.12 36.76 18.27 17.78 18.03 18.78 18.64 25.06 25.95
Ratarata (%b.k) 39.53
24.39
37.45
37.47
20.34
31.80
36.21
18.30
26.68
Kadar Air %b.b Ratarata (%b.b) 28.33 28.33 19.38 19.44 19.68 19.61 19.32 20.21 27.13 26.39 27.22 27.24 28.06 27.41 27.73 27.15 26.82 27.26 27.39 27.19 16.79 16.45 16.99 16.90 16.79 17.48 22.87 22.66 25.78 24.11 24.10 25.14 25.75 26.93 26.26 26.58 27.07 26.88 15.45 15.09 15.28 15.47 15.81 15.71 20.04 21.05 20.61
Ratarata (%b.b) 28.33
24.70
22.53
20.55
42
Jam ke
Lapisan
Tengah 6 Atas
Bawah
8
Tengah
Atas
Bawah
10
Tengah
Atas
12
Bawah
Poin T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki
%b.k
26.52 27.74 28.13 31.82 33.09 32.92 33.41 36.67 15.40 14.91 16.07 15.96 16.09 22.02 22.73 23.81 24.64 25.29 24.88 27.61 27.21 27.95 28.37 14.46 13.42 14.51 15.04 14.00 19.49 19.48 21.44 21.35 20.85 24.27 23.74 23.31 24.12 26.42 13.08 12.43 13.65 13.43
Ratarata (%b.k) 26.68
33.58
15.68
23.70
27.20
14.28
20.52
24.37
13.12
Kadar Air %b.b Ratarata (%b.b) 20.96 21.05 21.72 21.95 24.14 24.86 24.77 25.13 25.04 26.83 13.34 12.97 13.85 13.56 13.76 13.86 18.05 18.52 19.23 19.15 19.77 20.18 19.92 21.64 21.39 21.38 21.84 22.10 12.63 11.83 12.67 12.50 13.07 12.28 16.31 16.30 17.66 17.02 17.59 17.25 19.53 19.19 18.91 19.59 19.43 20.90 11.56 11.06 11.60 12.01 11.84
Ratarata (%b.b)
20.55
18.03
16.37
14.86
43
Jam ke
Lapisan
Poin
Bawah
B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te B-De B-Be B-Ka B-Ki B-Te T-De T-Be T-Ka T-Ki T-Te A-De A-Be A-Ka A-Ki A-Te
Tengah 12
Atas
Bawah
14
Tengah
Atas
Bawah
16
Tengah
Atas
%b.k
13.02 17.65 16.43 19.07 18.63 19.10 20.53 21.06 21.34 21.75 22.29 11.83 12.50 11.75 12.08 11.85 16.23 15.37 17.20 16.51 17.57 20.46 17.82 18.24 18.43 19.94 11.62 11.26 11.64 11.26 11.29 14.53 13.79 15.38 14.99 14.82 16.71 17.26 17.15 16.91 17.87
Ratarata (%b.k) 13.12
18.18
21.39
12.00
16.58
18.98
11.41
14.70
17.18
Kadar Air %b.b Ratarata (%b.b) 11.52 11.60 15.00 14.11 16.02 15.37 15.71 16.04 17.03 17.40 17.58 17.62 17.86 18.22 10.58 11.11 10.51 10.71 10.78 10.60 13.97 13.33 14.68 14.22 14.17 14.94 16.98 15.12 15.43 15.94 15.57 16.62 10.41 10.12 10.42 10.24 10.12 10.14 12.69 12.12 13.33 12.82 13.04 12.91 14.31 14.72 14.64 14.66 14.47 15.16
Ratarata (%b.b)
14.86
13.63
12.57
44
Lampiran 7 Data kecepatan aliran udara hasil pengukuran dan simulasi CFD pada beberapa variasi kecepatan udara
No. 1 2 3 4 5 6 7 8 9 10 11 12 13
13.2
Posisi 25,20,0 -25,20,0 0,20,25 0,20,-25 25,45,0 -25,45,0 0,45,25 0,45,-25 0,70,0 25,70,0 -25,70,0 0,70,25 0,70,-25 rata2 error
Velocity (m/s) 6.4
12.8
3.2
1.92
Ukur
CFD
error
Ukur
CFD
error
Ukur
CFD
error
Ukur
CFD
error
0.69 0.82 0.63 0.62 0.53 0.56 0.61 0.61 0.18 0.26 0.21 0.18 0.21
0.63 0.65 0.70 0.73 0.68 0.68 0.68 0.68 0.24 0.25 0.20 0.20 0.23
8.0% 20.2% 12.0% 17.7% 28.3% 22.5% 11.5% 11.5% 37.1% 2.0% 2.4% 11.1% 12.2% 15.1%
0.61 0.99 0.63 0.52 0.57 0.62 0.60 0.50 0.23 0.19 0.21 0.25 0.26
0.60 0.63 0.67 0.69 0.65 0.65 0.65 0.65 0.23 0.24 0.18 0.18 0.22
1.0% 36.4% 7.2% 34.0% 14.0% 5.7% 8.9% 30.0% 0.0% 26.3% 14.3% 26.5% 13.7% 16.8%
0.42 0.50 0.36 0.34 0.35 0.40 0.38 0.35 0.20 0.11 0.11 0.17 0.17
0.31 0.32 0.34 0.34 0.32 0.33 0.32 0.33 0.16 0.16 0.15 0.14 0.16
26.2% 35.4% 4.2% 1.5% 8.6% 16.5% 14.7% 4.3% 17.9% 52.4% 42.9% 17.6% 5.9% 19.1%
0.23 0.35 0.22 0.21 0.25 0.22 0.24 0.20 0.06 0.05 0.05 0.06 0.04
0.15 0.19 0.18 0.21 0.17 0.16 0.18 0.16 0.07 0.10 0.09 0.09 0.10
34.8% 44.9% 20.5% 2.4% 32.7% 29.5% 23.4% 20.5% 16.7% 122.2% 70.0% 54.5% 150.0% 47.9%
Uku r 0.07 0.10 0.06 0.08 0.08 0.09 0.10 0.13 0.05 0.06 0.05 0.06 0.04
CFD
error
0.09 0.10 0.10 0.10 0.10 0.10 0.10 0.10 0.05 0.05 0.04 0.03 0.06
38.5% 2.1% 66.7% 25.0% 33.3% 17.6% 0.0% 20.0% 11.1% 16.7% 11.1% 50.0% 71.4% 28.0%
Lampiran 8 Data suhu selama pengeringan gabah suhu (oC) jam keNo. Posisi Awal 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 8:30 9:30 10:30 11:30 12:30 13:30 14:30 15:30 16:30 17:30 18:30 19:30 20:30 21:30 22:30 23:30 0:30 7 0,0,0 34.0 39.7 39.4 39.6 39.2 39.3 38.4 39.1 40.2 40.4 40.6 40.5 42.5 41.5 41.3 40.6 40.3 8 20,0,0 32.3 39.2 39.1 39.5 39.2 39.0 38.2 39.0 40.1 40.3 42.4 40.3 42.2 41.3 41.0 40.4 40.1 9 -20,0,0 34.2 39.9 39.6 39.4 39.1 39.0 38.2 38.9 40.1 40.2 40.4 40.5 42.4 41.5 41.3 40.5 40.3 10 0,0,20 34.7 39.9 39.6 39.9 39.5 39.4 38.5 39.2 40.3 40.3 40.5 40.5 42.5 41.6 41.3 40.5 40.2 11 0,0,-20 36.1 40.3 39.9 40.2 40.1 39.5 38.6 39.3 40.5 40.5 40.6 40.8 42.9 42.0 41.7 40.9 40.6 rata-rata 34.3 39.8 39.5 39.7 39.4 39.2 38.4 39.1 40.2 40.3 40.9 40.5 42.5 41.6 41.3 40.6 40.3 12 0,25,0 28.3 27.9 28.2 30.1 32.3 32.9 33.5 34.1 36.6 36.5 36.8 37.1 38.4 39.1 38.9 39.2 39.2 13 20,25,0 28.2 27.9 28.3 30.2 32.2 32.8 33.4 34.3 36.5 36.4 36.5 37.0 38.3 39.0 38.8 39.1 39.0 14 -20,25,0 28.1 27.9 28.5 30.8 33.3 33.5 33.8 34.2 37.2 37.0 37.2 37.7 38.9 39.5 39.2 39.3 39.4 15 0,25,20 28.3 27.9 28.9 29.8 32.4 33.2 33.6 33.8 36.9 36.8 37.0 37.3 38.5 39.2 39.0 38.9 39.0 16 0,25,-20 28.4 28.3 30.4 30.7 33.1 34.0 34.3 34.9 37.8 37.6 37.8 38.1 39.5 39.9 39.6 39.5 39.6 rata-rata 28.3 28.0 28.9 30.3 32.7 33.3 33.7 34.3 37.0 36.9 37.1 37.4 38.7 39.3 39.1 39.2 39.2 17 0,50,0 28.5 28.2 27.7 27.8 28.0 27.6 28.2 29.6 31.6 33.0 33.2 33.4 34.7 35.8 35.9 37.1 37.3 18 20,50,0 28.2 27.8 27.4 27.0 27.1 28.1 30.2 31.4 34.6 34.7 34.8 35.4 36.7 37.9 37.7 37.9 38.1 19 -20,50,0 28.0 27.5 27.1 27.6 27.7 27.3 28.4 30.5 34.2 34.7 34.9 35.6 36.8 37.7 37.7 37.9 38.1 20 0,50,20 28.0 27.6 27.1 27.1 26.8 26.7 27.0 27.9 31.8 33.1 33.4 33.8 35.0 35.7 36.7 36.6 36.8 21 0,50,-20 28.1 27.9 27.5 28.1 27.9 27.4 27.5 29.5 33.8 34.3 24.7 34.8 36.2 37.2 37.4 37.7 37.8 rata-rata 28.2 27.8 27.4 27.5 27.5 27.4 28.3 29.8 33.2 34.0 32.2 34.6 35.9 36.9 37.1 37.4 37.6 Inlet 22 44.6 45.0 44.9 39.1 41.1 39.3 41.3 42.7 43.6 42.5 41.1 43.0 43.8 43.1 42.9 41.8 41.5 RH Plenum (%) 37.63 35.81 22.54 41.35 41.63 42.75 32.94 33.7 41.81 39.34 37.72 39.21 37.07 34.06 36.33 35.85 35.12 23 Tbb 29.9 29.6 25.2 27.7 34.7 31.4 32.6 33.7 31.5 33.8 33.8 33.7 34.1 33.1 33.2 31.7 31.3 24 Tbk 43.5 43.8 43.6 39.4 42.8 39.2 40.5 42.3 39.5 42.1 42.4 42.0 42.9 42.2 41.9 40.9 40.5 RH keluar (%) 98.51 97.02 95.53 96.25 96.23 92.48 88.72 76.72 70.14 69.35 65.11 63.48 62.55 43.77 58.9 56.37 65.73 45
46
No.
Posisi
25 Tbb Tbk 26 RH Lingkungan (%) Tbk 27 Tbb 28
suhu (oC) jam ke1 2 3 4 5 6 7 8 9 Awal 8:30 9:30 10:30 11:30 12:30 13:30 14:30 15:30 16:30 17:30 28.2 27.9 27.5 27.3 27 26.2 25.2 24.9 27.8 28.3 28.4 28.3 28.1 27.8 27.5 27.2 26.7 28.2 32.5 33.2 72.15 68.55 71.64 62.9 55.35 45.91 40.79 37.76 50.68 53.25 30.6 31.7 32.0 34.4 35.3 37.4 37.8 38.0 37.3 35.1 26.4 26.8 27.6 28.2 27.5 27.2 26.3 25.7 28.2 26.9
10 18:30 27.9 33.6 77.52 33.2 29.7
11 19:30 27.6 33.6 56.73 31.6 24.6
12 20:30 28.4 34.7 73.5 34.8 30.5
13 21:30 25.9 36.4 77.25 27.9 24.7
14 22:30 29 36.2 72.42 28.7 24.7
15 23:30 29 36.8 75.4 28.2 24.7
16 0:30 29 34.7 75.73 27.5 24.1
47
Lampiran 9 Gambar Teknik Pengering
Bak
Lubang sampel
Lantai berlubang Inlet Plenum
Lubang manometer plenum
Lubang manometer inlet Roda
Isometric View
DEPARTEMEN TEKNIK MESIN DAN BIOSISTEM FAKULTAS TEKNOLOGI PERTANIAN INSTITUT PERTANIAN BOGOR Nama Mesin : Pengering Gabah Tipe Bak Dirancang
Ahmad Fansuri
Digambar
Ahmad Fansuri
Skala : 1 : 10
Diperiksa
Leopold O. Nelwan
Satuan : cm
Disetujui
Leopold O. Nelwan Nama
Nama Bagian Gambar Piktorial Pengering Gabah Tipe Bak
12/09/2013 Material
A4 Tanggal
Paraf
No. Gambar :
48
1 3 4 8
NO.
2
5
NAMA BAGIAN
6
BAHAN
JUMLAH
UKURAN
1 Bak Pengering
Besi Plat Esser
1
tebal 2 mm
2 Lubang sampel bahan
Besi Sok Drat Dalam + Tutup
8
1 inch
3 Lubang Inlet
Besi Plat Esser
1
tebal 2 mm
4 Lantai Pengering
Besi Plat berlubang
1
tebal 1.2 mm, lubang 2 mm
5 Plenum
Besi Plat Esser
1
tebal 2 mm
6 Lubang manometer plenum
Besi Sok Drat Dalam + Tutup , Keni Besi Sok Drat Dalam + Tutup
1
3/4 inch
7 Lubang manometer inlet
7
KETERANGAN
3/4 inch
8 Roda
4
DEPARTEMEN TEKNIK MESIN DAN BIOSISTEM FAKULTAS TEKNOLOGI PERTANIAN INSTITUT PERTANIAN BOGOR Nama Mesin : Pengering Gabah Tipe Bak Dirancang
Ahmad Fansuri
Digambar
Ahmad Fansuri
Skala : 1 : 15
Diperiksa
Leopold O. Nelwan
Satuan : cm
Disetujui
Leopold O. Nelwan Nama
Nama Bagian Gambar Komposisi Assembly Pengering Gabah Tipe Bak
12/09/2013
Material A4 Tanggal
Paraf
No. Gambar :
22 47 9
18
90
49
40 80 Front View
5 10 Right View
DEPARTEMEN TEKNIK MESIN DAN BIOSISTEM FAKULTAS TEKNOLOGI PERTANIAN INSTITUT PERTANIAN BOGOR Nama Mesin : Pengering Gabah Tipe Bak Dirancang
Ahmad Fansuri
Digambar
Ahmad Fansuri
Skala : 1 : 12
Diperiksa
Leopold O. Nelwan
Satuan : cm
Disetujui
Leopold O. Nelwan Nama
Nama Bagian
Gambar Ortogonal Pengering Gabah Tipe Bak
12/09/2013 Material
A4 Tanggal
Paraf
No. Gambar :
50 Lampiran 10 Sistem Pengering Tipe Bak Lingkaran yang diujicoba dalam penelitian
a. Pengering Tipe Bak Lingkaran yang terintegrasi pemanas dan blower
b. Bagian dalam bak pengering dan peletakkan rangka besi tempat termokopel
51
c. Pengukuran kecepatan udara di dalam tumpukan gabah dengan anemometer
d. Pengukuran kecepatan udara di atas tumpukan gabah dengan anemometer
52 Lampiran 11 Alat-alat ukur yang digunakan dalam penelitian
a. Hybrid recorder
a. Timbangan digital
c. Sampel picker
b. Anemometer
b. Oven drying
53
RIWAYAT HIDUP Penulis dilahirkan dengan nama lengkap Ahmad Fansuri pada 20 November 1989 di Jakarta. Penulis merupakan anak ketiga dari lima bersaudara, putera dari H.Rojalih Hasan (Bapak) dan Rohimah (Ibu). Penulis menamatkan sekolah menengah atas di SMA Negeri 78 Jakarta tahun 2009 dan melanjutkan studi di Departemen Teknik Pertanian, Fakultas Teknologi Pertanian, IPB pada tahun yang sama melalui jalur ujian Seleksi Nasional Masuk Perguruan Tinggi Negeri (SNMPTN). Penulis pernah melakukan praktik lapangan di PT PG Rajawali II Unit PSA Palimanan dengan judul Aspek Keteknikan dalam Pengolahan Etanol di PSA Palimanan, Cirebon. Selama menjadi mahasiswa penulis aktif sebagai anggota Engineering Design Club 2011/2012 dan Koordinator Multimedia Engineering Design Club 2012/2013. Penulis juga pernah menjadi asisten praktikum mata kuliah Gambar Teknik tahun ajaran 2011/2012 dan 2012/2013.
