ANALISIS DISTRIBUSI SUHU, ALIRAN UDARA, RH DAN KADAR AIR DALAM IN-STORE DRYER (ISD) UNTUK BIJI JAGUNG
DISWANDI NURBA
SEKOLAH PASCASARJANA INSTITUT PERTANIAN BOGOR BOGOR 2008
PERNYATAAN MENGENAI TESIS DAN SUMBER INFORMASI Dengan ini saya menyatakan bahwa tesis Analisis Distribusi Suhu, Aliran Udara, RH dan Kadar Air dalam In-Store Dryer (ISD) untuk Biji Jagung adalah karya saya dengan arahan dari komisi pembimbing dan belum diajukan dalam bentuk apapun kepada perguruan tinggi manapun. Sumber informasi yang berasal dan dikutip dari karya yang diterbitkan maupun tidak diterbitkan dari penulis lain telah disebutkan dalam teks dan dicantumkan dalam Daftar Pustaka di bagian akhir tesis ini.
Bogor,
Juli 2008
Diswandi Nurba NRP F151060061
ABSTRACT DISWANDI NURBA. Analysis of Temperature, Air Flow, RH and Water Content Distribution inside In-Store Dryer (ISD) for Shelled Corn. Academic advisor: DYAH WULANDANI, Y. ARIS PURWANTO and RAFFI PARAMAWATI. In-Store Dryer (ISD) is commonly used as second step in drying process of grains. ISD is usually utilized ambient temperature and consist of dryer and storage system. Temperature, air flow and RH are key parameters during drying process using ISD. The objective of this study were to analyze the distribution of temperature, air flow, RH and water content inside ISD and to analyze the quality of shelled corn during drying process. Computational Fluid Dynamic (CFD) was used to analyze the distribution of temperature, air flow and RH inside ISD. Deep bed drying process was used to analyze the distribution of water content. The capacity of ISD used in this study was 7500 kg of corn, with dimensions of 3.5 m in high and 2.5 m in diameter. ISD have 13 aeration pipes of air flow, consist of 9 input pipes and 4 output pipes. All walls of ISD were assumed to be in condition of adiabatic. Validation of water content that simulation with the measurement has been done at capacity of ISD were 1500 kg. The result showed that CFD simulation of temperature, air flow and RH have coefficient of correlation of 0.66, 0.73 and 0.66 respectively. Deep bed drying simulation of water content at Layer 10 and 40 have coefficient of correlation of 0.90 and 0.35 in rainy season, 0.88 and 0.84 in dry season. The results showed that only a minor changes in the quality of shelled corn dried and stored using ISD. Keyword: CFD, ISD, simulation, analysis of temperature, deep bed drying.
RINGKASAN DISWANDI NURBA. Analisis Distribusi Suhu, Aliran Udara, RH dan Kadar Air dalam In-Store Dryer (ISD) untuk Biji Jagung. Dibimbing oleh DYAH WULANDANI, Y. ARIS PURWANTO dan RAFFI PARAMAWATI. Pengeringan dan penyimpanan merupakan proses penting dalam penanganan pascapanen biji-bijian dan produk pertanian pada umumnya. Penanganan tersebut dimaksudkan untuk mengurangi susut dan turunnya kualitas hasil panen sehingga dapat bertahan lebih lama. Penelitian ini merupakan bagian dari penelitian “Rancang Bangun Alat Pengering Efek Rumah Kaca (ERK)-Hybrid dan In-Store Dryer (ISD) Terintegrasi untuk Biji-Bijian”. ISD sebagai satu unit dari sistem terintegrasi tersebut, memerlukan pengkajian secara spesifik untuk melihat performa operasinya secara lebih mendalam yang juga akan menjadi bahan evaluasi bagi pengembangan alat pengering dan penyimpan terintegrasi ini. Metode yang digunakan pada ISD adalah pemanfaatan udara lingkungan yang dihembuskan melalui tumpukan biji-bijian yang akan dikeringkan. Sebagai sebuah sistem pengeringan konvektif, yang mengandalkan aliran udara (gas) sebagai media utama untuk keberhasilan proses pengeringan maka penataan sistem saluran udara di dalam tumpukan biji-bijian dapat memberikan sebaran kondisi udara yang lebih seragam. Secara umum tujuan penelitian ini adalah untuk mengkaji distribusi aliran udara, suhu, RH dan kadar air di dalam ISD. Secara khusus penelitian ini terdiri dari beberapa tujuan antara lain : 1) melakukan simulasi dan validasi model distribusi aliran udara dan suhu di dalam ISD dengan menggunakan teknik Computational Fluid Dynamics (CFD), 2) mendapatkan nilai perhitungan dan validasi RH dengan Microsoft excel, 3) melakukan simulasi dan validasi kadar air jagung menggunakan model pengeringan tumpukan biji-bijian dengan pemrograman Visual Basic 6.0, 4) analisis mutu jagung hasil pengeringan dan peyimpanan di dalam ISD berdasarkan SNI. Simulasi sistem thermal 3D pada bagunan ISD dibuat dengan menggunakan software CFD yaitu Gambit 2.2.30 & Fluent 6.1.18. Simulasi dilakukan dengan dua kondisi terhadap pipa saluran udara yang ditempatkan di dalam ISD, yaitu: Simulasi 1 bangunan ISD dikondisikan sesuai dengan bangunan ISD di lapangan dengan 9 buah pipa input setengah berpori, sisi dari pipa yang berpori diposisikan menghadap ke dinding. Sementara 4 buah pipa output seluruhnya berpori. Simulasi 2 bangunan ISD dengan pipa input dan output yang seluruhnya berpori, simulasi ini merupakan modifikasi pada jenis pipa input. Pada simulasi ini seluruh pipa input dikondisikan memiliki pori keseluruhan pada permukaannya. Simulasi pengeringan tumpukan jagung dibuat dalam program komputer Visual Basic 6.0, untuk menyelesaikan persamaan-persamaan secara simultan meliputi: a) persamaan keseimbangan massa, b) persamaan laju pengeringan, c) persamaan keseimbangan panas dan d) persamaan laju perpindahan panas. Simulasi dilakukan dengan kondisi: ketebalan tumpukan 2.50 m, kadar air awal biji jagung 18%, suhu biji jagung 29.5 oC dan laju massa udara 12.7 kg/mnt-m2. Kondisi suhu udara dan RH input dibedakan untuk dua musim yaitu: pada musim
hujan dengan suhu udara masuk 31 oC dan RH 73%, sementara untuk musim kemarau dengan suhu udara masuk 33 oC dan RH 59.8%. Validasi dilakukan untuk membandingkan hasil pengukuran dan hasil simulasi pada titik-titik dan lokasi tertentu. Kriteria hasil validasi dianalisis dengan metode curve-fitting dan standar deviasi. Sementara besarnya error dalam validasi dihitung dengan persamaan mean absolute error (MAE), disamping itu juga ditentukan korelasi antara data simulasi dan hasil pengukuran dengan koefisien korelasi. Validasi data pada CFD meliputi suhu, aliran udara dan RH, sementara validasi data kadar air bahan dilakukan berdasarkan hasil simulasi model pengeringan tumpukan pada Visual Basic 6.0. Selanjutnya dilakukan pengujian mutu meliputi parameter-parameter yang menjadi persyaratan mutu jagung untuk perdagangan, untuk benih dan juga untuk dijadikan pakan ternak menurut SNI, meliputi; kadar air, butir rusak, butir warna lain, butir pecah, kotoran, kandungan nutrisi (kimia), dan kontaminasi aflatoxin. Disamping itu juga dilakukan uji viabilitas dengan metode perkecambahan. Hasil analisis distribusi suhu, aliran udara dan RH menunjukkan bahwa tingkat keseragaman sebaran suhu, aliran udara dan RH pada ISD dengan menggunakan jenis pipa input dengan pori seluruhnya lebih seragam dibandingkan jenis pipa input setengah berpori. Sementara validasi antara nilai simulasi CFD terhadap nilai pengukuran didapat korelasi untuk sebaran suhu, kecepatan aliran udara dan RH sebesar 0.66, 0.73 dan 0.66. Hasil analisis perubahan kadar air biji jagung dengan simulasi pengeringan tumpukan didapatkan: pada musim hujan membutuhkan waktu 150 jam untuk mengurangi kadar air dari 18% b.k mencapai kadar air keseimbangan 16.6% b.k. Sementara pada musim kemarau membutuhkan waktu 120 jam untuk mengurangi kadar air dari 18% b.k mencapai kadar air keseimbangan 13.5% b.k. Validasi antara hasil simulasi model pengeringan tumpukan terhadap nilai pengukuran pada layer 10 dan 40 didapatkan nilai korelasi sebesar 0.90 dan 0.35 pada musim hujan, pada musim kemarau didapatkan nilai korelasi sebesar 0.88 dan 0.84. Hasil pengujian mutu sampel jagung setelah proses pengeringan selama 40 jam dan penyimpanan selama 30 hari dalam ISD, kontaminasi aflatoxin pada jagung tidak jauh berubah, yaitu dari rata-rata 18.48 ppb sebelum proses menjadi 21.10 ppb. kandungan tersebut masih berada di bawah ambang batas toleransi yang ditetapkan SNI untuk pakan ternak sebesar 50 ppb. Sementara komposisi nutrisi yang diuji juga masih memenuhi standar SNI. Dari analisis mutu biji jagung menunjukkan bahwa proses pengeringan dan penyimpanan di dalam ISD mampu mempertahankan mutu biji jagung dengan baik sebagai bahan pakan ternak. Kata kunci : CFD, ISD, simulasi, analisis distribusi suhu, pengeringan tumpukan.
@ Hak Cipta milik IPB, tahun 2008 Hak Cipta dilindungi Undang-Undang 1. Dilarang mengutip sebagian atau seluruh karya tulis ini tanpa mencantumkan atau menyebutkan sumber a. Pengutipan hanya untuk kepentingan pendidikan, penelitian, penulisan karya ilmiah, penyusunan laporan, penulisan kritik, atau tinjauan suatu masalah b. pengutipan tersebut tidak merugikan kepentingan yang wajar IPB 2. Dilarang mengumumkan dan memperbanyak sebagian atau seluruh Karya tulis dalam bentuk apa pun tanpa izin IPB.
ANALISIS DISTRIBUSI SUHU, ALIRAN UDARA, RH DAN KADAR AIR DALAM IN-STORE DRYER (ISD) UNTUK BIJI JAGUNG
DISWANDI NURBA
Tesis sebagai salah satu syarat untuk memperoleh gelar Magister Sains pada Program Studi Ilmu Keteknikan Pertanian
SEKOLAH PASCASARJANA INSTITUT PERTANIAN BOGOR BOGOR 2008
Penguji Luar Komisi pada Ujian Tesis : Dr. Ir. Rokhani Hasbullah, M.Si
Judul Tesis Nama NRP
: Analisis Distribusi Suhu, Aliran Udara, RH dan Kadar Air dalam In-Store Dryer (ISD) untuk Biji Jagung : Diswandi Nurba, S.TP : F151060061
Disetujui Komisi Pembimbing
Dr. Ir. Dyah Wulandani, M.Si Ketua
Dr. Ir. Y. Aris Purwanto, M.Sc Anggota
Dr. Ir. Raffi Paramawati, M.Si Anggota
Diketahui Ketua Program Studi Ilmu Keteknikan Pertanian
Dekan Sekolah Pascasarjana
Prof. Dr. Ir. Armansyah H. Tambunan, M.Agr
Prof. Dr. Ir. Khairil A. Notodiputro, MS
Tanggal Ujian :
Tanggal Lulus :
PRAKATA Puji syukur senantiasa kita panjatkan ke hadirat Allah SWT sebagai pemilik segala kesempurnaan, dan shalawat serta salam kepada Rasulullah Muhammad SAW manusia mulia sebagai uswatun hasanah dalam kehidupan ini hingga akhir zaman kelak. Dengan mengucapkan Alhamdulillahirabbil ‘alamin dan mengharap ridho-Nya, penulis telah dapat menyelesaikan tesis dengan judul “Analisis Distribusi Suhu, Aliran Udara, RH dan Kadar Air Dalam In-Store Dryer (ISD) untuk Biji Jagung”. Pada kesempatan ini penulis mengucapkan rasa terimakasih yang mendalam kepada : 1.
Dr. Ir. Dyah Wulandani, M.Si selaku ketua komisi pembimbing, atas segala arahan dan bimbingannya yang sangat berharga bagi penulis selama pendidikan, penelitian dan penyelesaian tesis.
2.
Dr. Ir. Y. Aris Purwanto, M.Sc sebagai anggota komisi pembimbing, atas segala koreksian, bimbingan dan motivasinya.
3.
Dr. Ir. Raffi Paramawati, M.Si sebagai anggota komisi pembimbing, yang telah meluangkan waktunya, pemikiran dan masukan-masukan yang sangat membantu penulis dalam menyelesaikan penelitian dan penulisan tesis.
4.
Dr. Leopold Oscar Nelwan, STP.,M.Si selaku ketua peneliti pada proyek penelitian KKP3T atas kepercayaan terhadap penulis sebagai bagian dalam Tim Peneliti dan juga atas segala masukan dan arahannya.
5.
Prof. Dr. Ir. Armansyah H. Tambunan, M.Agr selaku ketua Program Studi Ilmu Keteknikan Pertanian.
6.
Dr. Ir. Rokhani Hasbullah, M.Si selaku dosen penguji luar komisi pembimbing pada ujian tesis, atas segala masukan dan saran bagi penulisan tesis ini.
7.
Depertemen Pendidikan Nasional RI, khususnya DIKTI melalui Program BPPS atas bantuan biaya pendidikan yang diberikan.
8.
Departemen Pertanian RI melalui Proyek Penelitian KKP3T Tahun 2007 yang telah membantu membiayai penelitian.
9.
Bapak dan Ibu Staf Pengajar Departemen Teknik Pertanian, dan khususnya kepada Staf Laboratorium Energi dan Elektrifikasi Pertanian FATETA IPB, Pak Harto, Mas Firman dan Mas Darma, terima kasih atas semua bantuan teknis selama perkuliahan, praktikum dan juga saat melaksanakan penelitian.
10. Teman-teman satu tim peneliti; Mas Lilik, Mas Deni, dan Kak Tamaria, terimakasih atas kerjasamanya. Seluruh teman-teman angkatan 2006; Mas Susanto, Mas Warji, Mas Surya, Mas Farry, dan Kak Riswanti, terimakasih atas kebersamaannya selama pendidikan. 11. Teman-teman Prodi TEP : Bang Hendri, Mas Nuruddin, Mas Bayu, Bang Yaziz, Bang Iqbal dan seluruh rekan-rekan yang tidak dapat saya sebutkan satu persatu, terimakasih semuanya. 12. Rekan-rekan IKAMAPA : Pak Samingan, Pak Ali, Bang Daud, Bang Faisal, Bang Safrizal dan seluruh anggota IKAMAPA. Ketulusan kasih sayang, pengorbanan serta do’a yang tiada henti dari Ayahanda dan Ibunda selama ini adalah penyemangat bagi penulis dalam menyelesaikan studi ini, ketulusan yang tidak mungkin akan terbalas. Do’a dan dukungan dari Kakak, Dinda dan Adik-adikku serta seluruh keluarga merupakan dorongan yang memberikan energi positif dalam menjalani pendidikan ini. Penulis telah berupaya optimal untuk menghasilkan karya ilmiah yang baik, namun demikian tentunya masih sangat banyak kekurangan dan jauh dari kesempurnaan, kiranya kritik, saran dan koreksian sangat kami harapkan demi perbaikan dan kontribusi terhadap ilmu pengetahuan di masa yang akan datang. Akhirnya dengan segala keterbatasan yang ada, seiring doa dan harapan penulis, semoga karya ilmiah ini dapat bermanfaat. Semoga kita semua selalu mendapat bimbingan-Nya dalam mengamalkan ilmu dengan baik dan menjadi hamba-Nya yang selalu bersyukur.
Bogor,
Juli 2008
Diswandi Nurba
RIWAYAT HIDUP Diswandi Nurba dilahirkan di Manjeng pada tanggal 28 April 1982, adalah putra kedua dari empat bersaudara dari Ayahanda Bangsawan dan Ibunda Nur Asiah. Penulis lulus dari SMU Negeri 1 Meulaboh pada tahun 2000 dan melanjutkan pendidikan Sarjana pada Jurusan Teknik Pertanian Fakultas Pertanian Universitas Syiah Kuala Banda Aceh. Selama pendidikan S1, Penulis menjadi Asisten pada Laboratorium Alat dan Mesin Pertanian tahun 2004 dan juga menjabat sebagai Sekjend Pemerintah Mahasiswa (PEMA) Unsyiah periode 2004-2005. Pada tahun 2005 Penulis menyelesaikan pendidikan S1 dan selanjutnya mengabdi sebagai Staf Pengajar pada Sekolah Tinggi Ilmu Pertanian (STIP) Teungku Dirundeng Meulaboh, disamping itu pada tahun 2006 Penulis juga mengisi waktu sebagai Staf Teknis Bidang Pemberdayaan Pertanian dan Perikanan Dewan Pengawas BRR NAD dan Nias. Pada medio Agustus tahun 2006, Penulis melanjutkan pendidikan Magister pada Program Studi Ilmu Keteknikan Pertanian Sekolah Pascasarjana Institut Pertanian Bogor, dengan sponsor BPPS Ditjen DIKTI. Selama pendidikan S2, Penulis menjabat sebagai Sekretaris Ikatan Keluarga Mahasiswa Pascasarjana Aceh (IKAMAPA) periode 2006-2007 dan Sekretaris Forum Mahasiswa Pascasarjana Ilmu Keteknikan Pertanian (FORMATETA) IPB periode 2007-2008.
xii
DAFTAR ISI Halaman DAFTAR TABEL ............................................................................................
xiii
DAFTAR GAMBAR .......................................................................................
xv
DAFTAR LAMPIRAN .................................................................................... xvii DAFTAR SIMBOL..........................................................................................
xx
I
PENDAHULUAN .................................................................................... 1.1 Latar Belakang ................................................................................... 1.2 Hipotesa ............................................................................................. 1.3 Tujuan ................................................................................................
1 1 5 5
II
TINJAUAN PUSTAKA ........................................................................... 2.1 Pengeringan........................................................................................ 2.1.1 Proses Pengeringan ................................................................... 2.1.2 Paramater Pengeringan ............................................................. 2.1.3 Aliran Udara Pengeringan ........................................................ 2.1.4 Karakteristik Pengeringan Jagung ............................................ 2.1.5 Sorpsi Isotermi .......................................................................... 2.1.6 Aktivitas Air ............................................................................. 2.2 Penyimpanan ...................................................................................... 2.2.1 Pengaruh Kadar Air terhadap Penyimpanan ............................ 2.2.2 Kelembaban dan Suhu Penyimpanan ....................................... 2.2.3 Hubungan Antara Penyimpanan dan Kerusakan Bahan Pakan ........................................................................................ 2.2.4 Persyaratan Mutu Jagung ......................................................... 2.3 Perkembangan Penelitian In-Store Dryer .......................................... 2.4 Teknik Simulasi Computational Fluid Dynamics (CFD) .................. 2.4.1 Pre-processor ........................................................................... 2.4.2 Solver ........................................................................................ 2.4.3 Post-processor .......................................................................... 2.5 Model Pengeringan Tumpukan (Deep Bed Drying) .......................... 2.5.1 Keseimbangan Massa ............................................................... 2.5.2 Laju Pengeringan ...................................................................... 2.5.3 Keseimbangan Panas ................................................................ 2.5.4 Laju Perpindahan Panas ...........................................................
6 6 7 8 10 11 12 13 14 16 17
III METODE PENELITIAN .......................................................................... 3.1 Waktu dan Tempat ............................................................................. 3.2 Bahan dan Alat ................................................................................... 3.2.1 Bahan ........................................................................................ 3.2.2 Alat ........................................................................................... 3.2.3 Bangunan In-Store Dryer ......................................................... 3.3 Prosedur Penelitian ............................................................................
31 31 31 31 31 31 34
18 18 20 22 22 23 25 25 27 28 28 29
xiii
3.3.1 3.3.2 3.3.3 3.3.4 3.3.5 3.3.6 3.3.7 3.3.8 3.3.9
Simulasi Sistem Thermal ISD .................................................. Simulasi Pengeringan Tumpukan (Lapis Tebal) Jagung .......... Percobaan Pengeringan dan Penyimpanan Jagung ................... Distribusi Udara........................................................................ Pengukuran Kecepatan Udara dan Suhu .................................. Pengukuran RH ........................................................................ Pengukuran Kadar Air Jagung.................................................. Validasi Model Simulasi .......................................................... Uji Mutu Jagung Hasil Pengeringan dan Penyimpanan di dalam In-Store Dryer ............................................................
35 35 36 36 37 38 39 39
IV HASIL DAN PEMBAHASAN................................................................. 4.1 Simulasi Aliran Udara pada ISD........................................................ 4.1.1 Pembentukan Grid Perhitungan ............................................... 4.1.2 Distribusi Suhu dan Kecepatan Udara Pengering Simulasi 1 ................................................................................. 4.1.3 Distribusi Suhu dan Kecepatan Udara Pengering Simulasi 2 ................................................................................. 4.1.4 Distribusi RH Udara Hasil Simulasi......................................... 4.1.5 Keragaman Kecepatan Udara, Suhu dan RH ........................... 4.2 Validasi Suhu, Kecepatan Aliran Udara dan RH ............................... 4.3 Perubahan Kadar Air Jagung pada Simulasi ...................................... 4.3.1 Kadar Air Jagung Simulasi 1 .................................................... 4.3.2 Kadar Air Jagung Simulasi 2 .................................................... 4.4 Perubahan Kadar Air Jagung Percobaan............................................ 4.4.1 Kadar Air Jagung Percobaan 1 ................................................. 4.4.1 Kadar Air Jagung Percobaan 2 ................................................. 4.5 Validasi Perubahan Kadar Air Jagung ............................................... 4.5.1 Validasi Kadar Air Jagung Percobaan 1 ................................... 4.5.2 Validasi Kadar Air Jagung Percobaan 2 ................................... 4.6 Analisis Mutu Jagung Hasil Pengeringan dan Penyimpanan ............ 4.6.1 Mutu Perdagangan .................................................................... 4.6.2 Mutu Benih ............................................................................... 4.6.3 Mutu Pakan Ternak ..................................................................
41 41 41
45 48 48 53 56 56 58 60 60 62 63 63 65 67 67 69 70
V SIMPULAN DAN SARAN ......................................................................... 5.1 Simpulan ............................................................................................ 5.2 Saran ..................................................................................................
73 73 73
DAFTAR PUSTAKA ......................................................................................
74
LAMPIRAN .....................................................................................................
77
40
42
xiv
DAFTAR TABEL Halaman 1
Suhu udara pengering beberapa jenis biji-bijian menurut tujuan penggunaannya .......................................................................................
10
Kebutuhan volume aliran udara pengering pada berbagai cara pengeringan. ...........................................................................................
11
3
Hasil perhitungan kadar air berdasarkan Persamaan (8) ........................
12
4
Aktivitas air (aw) minimum untuk pertumbuhan mikroba dan perkecambahan spora .............................................................................
14
5
Perubahan biologi dan kimia pada pakan konsentrat .............................
18
6
Persyaratan mutu jagung ........................................................................
19
7
Standar mutu jagung bahan baku pakan ternak ......................................
19
8
Spesifikasi persyaratan mutu benih jagung hibrida di laboratorium............................................................................................
19
Lokasi pipa-pipa penyalur udara dalam ISD menurut fungsinya pada bidang xz. .......................................................................................
33
Koordinat lokasi titik-titik pengukuran suhu, kecepatan udara dan RH dalam ISD dengan termokopel ..................................................
38
Koordinat lokasi titik-titik pengambilan sampel untuk pengukuran kadar air bahan ....................................................................
39
2
9 10 11
xv
DAFTAR GAMBAR Halaman 1
Proses pengeringan pada kurva psychrometric.......................................
7
2
Sorpsi isotermi yang menunjukkan hysterisis ........................................
13
3
Batas-batas suhu dan kadar air yang aman pada penyimpanan biji-bijian (Hall 1970) .............................................................................
17
4
Elemen pada bak (Bala 1997) .................................................................
26
5
Grid finite different untuk persamaan deep bed drying ..........................
27
6
Skema Bangunan ISD .............................................................................
32
7
Diagram Alir Penelitian ..........................................................................
34
8
Model ISD 3 Dimensi untuk simulasi CFD............................................
37
9
Pembentukan grid pada domain perhitungan .........................................
41
10
Distribusi suhu udara di dalam ISD Simulasi 1 ......................................
42
11
Pengaruh pipa input berpori sebagian terhadap sebaran suhu ISD pada Simulasi 1. ......................................................................................
43
12
Distribusi kecepatan udara didalam ISD pada Simulasi 1 ......................
44
13
Pengaruh pipa input berpori sebagian terhadap vektor aliran udara dalam ISD pada Simulasi 1...........................................................
44
14
Distribusi suhu udara di dalam ISD Simulasi 2 ......................................
45
15
Pengaruh pipa input berpori sebagian terhadap sebaran suhu ISD pada Simulasi 2 .......................................................................................
46
16
Distribusi kecepatan udara di dalam ISD pada Simulasi 2 .....................
47
17
Pengaruh pipa input berpori sebagian terhadap vektor aliran udara dalam ISD pada Simulasi 2...........................................................
48
18
Profil suhu pada 5 ketinggian Simulasi 1 ...............................................
49
19
Profil suhu pada 5 ketinggian Simulasi 2 ...............................................
49
20
Keragaman suhu pada kedua simulasi ....................................................
50
21
Profil kecepatan aliran udara pada 5 ketinggian Simulasi 1 ...................
51
22
Profil kecepatan aliran udara pada 5 ketinggian Simulasi 2 ...................
51
23
Keragaman kecepatan aliran udara pada kedua simulasi .......................
52
xvi
24
Keragaman RH udara pada kedua simulasi ............................................
53
25
Validasi suhu udara hasil simulasi terhadap suhu pengukuran ..............
54
26
Validasi kecepatan aliran udara hasil simulasi terhadap data pengukuran. ............................................................................................
54
27
Validasi RH hasil perhitungan terhadap RH hasil pengukuran ..............
55
28
Perubahan kadar air pada Simulasi 1 ......................................................
56
29
Perubahan kadar air pada Simulasi 2 ......................................................
58
30
Kadar air hasil pengukuran selama 50 jam pengeringan pada Percobaan 1.............................................................................................
60
Pengaruh fluktuasi RH terhadap perubahan kadar air bijian pada Percobaan 1.............................................................................................
61
Kadar air hasil pengukuran selama 40 jam pengeringan pada Percobaan 2.............................................................................................
62
Pengaruh fluktuasi RH terhadap perubahan kadar air bijian pada Percobaan 2.............................................................................................
63
Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 10 pada Percobaan 1 .................................................
64
Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 40 pada Percobaan 1 .................................................
64
Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 10 pada Percobaan 2 .................................................
66
Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 40 pada Percobaan 2 .................................................
66
Perbandingan parameter mutu hasil percobaan dengan SNI jagung untuk perdagangan ......................................................................
68
Hasil pengujian tingkat kemurnian benih dan viabilitas sebelum dan setelah proses dalam ISD dan perbandingan dengan SNI ...............
69
Perbandingan kandungan nutrisi (kimia) jagung antara hasil percobaan dan persyaratan SNI ..............................................................
71
31 32 33 34 35 36 37 38 39 40
xvii
DAFTAR LAMPIRAN Halaman 1.
Sistem Pengering ERK-Hybrid dan In-Store Dryer Terintegrasi ...........
78
2.
Profil aliran udara pada pipa setengah berpori (Brooker et al. 1992) .......................................................................................................
79
Arah aliran udara pada pipa input dan output (Brooker et al. 1992) .......................................................................................................
80
4.
Susunan pipa-pipa penyalur udara dalam ISD........................................
81
5.
Standar ASAE untuk ukuran dan kapasitas Silo ....................................
82
6.
Algoritma numerik volume hingga dengan metode SIMPLE (Versteeg & Malalasekera 1995) ............................................................
84
Asumsi, kondisi awal dan kondisi batas yang digunakan pada simulasi CFD ..........................................................................................
85
8.
Algoritma simulasi pengeringan tumpukan tebal ...................................
87
9.
Parameter yang digunakan dalam simulasi pengeringan tumpukan ................................................................................................
88
Kode program Visual Basic untuk simulasi pengeringan tumpukan ................................................................................................
90
11.
Interface program simulasi pengeringan tumpukan ...............................
95
12.
Perhitungan pressure drop dan tekanan statis kipas ...............................
96
13.
Lokasi titik pengukuran suhu dan kecepatan udara di dalam ruangan ISD ............................................................................................
100
14.
Lokasi titik pengambilan sampel pengukuran kadar air .........................
101
15.
Hasil Simulasi 1 CFD .............................................................................
102
16.
Hasil Simulasi 2 CFD .............................................................................
105
17.
Perbandingan keragaman suhu, kecepatan udara dan RH di dalam ISD pada kedua simulasi CFD. ....................................................
108
Data validasi suhu dan kecepatan udara hasil pengukuran dan simulasi CFD serta nilai error dan standar deviasinya...........................
109
3.
7.
10.
18.
xviii
19.
Data validasi RH udara hasil pengukuran dan perhitungan serta nilai error dan standar deviasinya. .........................................................
110
Perubahan kadar air setelah simulasi pengeringan selama 150 jam pada Simulasi 1. ...............................................................................
111
Perubahan kadar air setelah simulasi pengeringan selama 120 jam pada Simulasi 2. ...............................................................................
114
Hasil pengukuran nilai kadar air bijian selama 50 jam pada Percobaan 1.............................................................................................
117
Pengaruh fluktuasi RH terhadap perubahan kadar air pengukuran pada Percobaan 1 ................................................................
120
Hasil pengukuran nilai kadar air bijian selama 40 jam pada Percobaan 2.............................................................................................
121
Pengaruh fluktuasi RH terhadap perubahan kadar air pengukuran pada Percobaan 2 ................................................................
123
Validasi perubahan kadar air pengukuran dan simulasi pada Percobaan 1.............................................................................................
124
Validasi perubahan kadar air pengukuran dan simulasi pada Percobaan 2.............................................................................................
125
28.
Mutu jagung pada percobaan ..................................................................
126
29.
Kontaminasi aflatoxin pada jagung hasil uji laboratorium sebelum proses pengeringan dan penyimpanan ......................................
127
Kandungan abu, lemak dan serat pada jagung hasil uji laboratorium sebelum proses pengeringan dan penyimpanan ...............
128
Kandungan protein pada jagung hasil uji laboratorium sebelum proses pengeringan dan penyimpanan ....................................................
129
Kontaminasi aflatoxin pada jagung hasil uji laboratorium setelah proses pengeringan dan penyimpanan ....................................................
130
Kandungan abu, lemak dan serat pada jagung hasil uji laboratorium setelah proses pengeringan dan penyimpanan .................
131
Kandungan protein pada jagung hasil uji laboratorium setelah proses pengeringan dan penyimpanan ....................................................
132
Sistem Pengering Efek Rumah Kaca Hybrid (ERK-Hybrid) dan In-Store Dryer Terintegrasi yang di ujicoba dalam penelitian ...............
133
20. 21. 22. 23. 24. 25. 26. 27.
30. 31. 32. 33. 34. 35.
xix
36.
Alat-alat ukur yang digunakan dalam penelitian ....................................
134
xx
DAFTAR SIMBOL aw
: aktivitas air (desimal)
Cpa
: panas jenis udara kering (J/kgoK)
Cpg
: panas jenis bijian (J/kgoK)
Cpl
: panas jenis air pada bijian (J/kgoK)
Cpw
: panas jenis uap air (J/kgoK)
Cv
: kalor molekul (K)
C2
: koefisien porous jump plat (1/m)
D
: diameter spesifik (m)
Dp
: diameter lubang pada plat (m)
Ga
: laju aliran massa udara (kg/mnt m2)
H
: kelembaban mutlak (kg/kg)
hcv
: coefisien panas volumetric air (kJ/mnt-m3-K)
Hfg
: panas laten penguapan (kJ/kg)
i
: energi dalam (J)
i
: data ke-i
k
: konstanta pengeringan (dalam mnt-1)
k
: konduktivitas panas (W/m K)
La
: panas laten penguapan air (kJ/kg)
Lg
: panas laten penguapan dari bijian (kJ/kg)
M
: kadar air bijian basis kering (% b.k)
MAE
: mean absolute error
Me
: kadar air keseimbangan (% b.k)
Mw
: kadar air bijian basis basah (% b.b)
N
: jumlah data
p
: tekanan parsial air (Pa)
pw
: tekanan keseimbangan uap air (Pa)
Patm
: tekanan atmosfer (Pa)
Ps
: tekanan jenuh air (Pa)
Pv
: tekanan uap (Pa)
Pr
: bilangan Prandtl (desimal)
Qc
: nilai hasil simulasi
xxi
Qo
: nilai hasil pengukuran
R
: konstanta gas ideal (J/mol K)
Re
: bilangan Reynold (desimal)
RH
: kelembaban nisbi (%)
RHa
: kelembaban udara lingkungan (%)
RHr
: kelembaban udara pengering (%)
Si
: sumber gerakan energi dalam
SMx
: sumber gerakan momentum arah x
SMy
: sumber gerakan momentum arah y
SMz
: sumber gerakan momentum arah z
t
: waktu (mnt)
T
: suhu (oC)
Ta
: suhu udara (oC)
Tg
: Suhu bijian (oC)
u
: exponen
u
: kecepatan arah x (m/s)
v
: kecepatan arah y (m/s)
w
: kecepatan arah z (m/s)
Wd
: bobot bahan kering (kg)
Ww
: bobot bahan basah (kg)
x
: koordinat arah x (m)
y
: koordinat arah y (m)
z
: ketebalan tumpukan biji (m)
z
: koordinat arah z (m)
α
: permeabilitas permukaan plat (m2)
ε
: porositas bak (desimal)
ρd
: massa jenis bijian (kg/m3)
ρ
: densitas fluida (kg/m3)
μ
: viskositas dinamik (kg/m.s)
I PENDAHULUAN 1.1 Latar Belakang Salah satu permasalahan utama dalam pascapanen komoditi biji-bijian adalah susut panen dan turunnya kualitas, sehingga perlu diupayakan metode pengeringan dan penyimpanan yang baik untuk menjaga dan mempertahankan kuantitas dan kualitasnya. Pengeringan adalah proses pemindahan air dengan menggunakan panas atau aliran udara untuk menghambat pertumbuhan jamur dan bakteri
sehingga
tidak
dapat
berkembang
lagi
atau
memperlambat
perkembangannya (Hall 1980). Penyimpanan hasil pertanian berhubungan dengan waktu penggunaan, baik distribusi maupun konsumsi atau pengolahan lebih lanjut. Penyimpanan bertujuan agar bahan tidak mengalami kerusakan dan penyusutan selama masa simpannya. Pada umumnya, penyimpanan biji-bijian dilakukan setelah proses pengeringan hingga kadar air yang dianggap aman. Pengeringan dapat dilakukan menggunakan cara alamiah (penjemuran) ataupun cara buatan (artificial drying). Penjemuran merupakan cara pengeringan yang cukup murah akan tetapi ada faktor yang menjadi kendala penjemuran yang mencakup: kebutuhan lahan yang luas, kontaminasi bahan asing, tidak praktis pada daerah yang sering berubah cuacanya dan pada musim hujan praktis sulit dilakukan. Pengeringan buatan dengan energi konvensional pada umumnya dapat digunakan untuk mengatasi hal tersebut, akan tetapi kendala utamanya adalah biaya pengoperasian yang relatif tinggi, terutama ketika akhir-akhir ini biaya bahan bakar meningkat pesat. Selain itu pada biji-bijian terutama jagung pipilan, kesalahan penanganan pada pengeringan dan penyimpanan dapat menyebabkan tingginya kontaminasi mikotoksin terutama jenis aflatoxin yang berbahaya bagi kesehatan ternak dan manusia. Salah satu tipe pengering berenergi surya yang telah dikembangkan adalah pengering tipe Efek Rumah Kaca (ERK) (Abdullah 1993). Pengering ini merupakan bangunan dengan struktur terintegrasi antara kolektor surya dengan wadah produk yang dikeringkan. Pengembangan alat pengering surya tipe ERK pada skala penelitian dan lapangan telah dilakukan diantaranya untuk berbagai produk biji-bijian (gabah dan jagung), perkebunan (kopi, kakao dan cengkeh),
2
buah-buahan (pisang dan pepaya), benih (cabai dan mentimun) dan ikan (Abdullah 1995, 1998, 1999; Nelwan 1997, 2005; Wulandani 2005; Manalu 1999). Suhu udara pengering rata-rata berkisar antara 39-50 oC untuk berbagai lokasi, dengan waktu pengeringan berkisar antara 4-57 jam tergantung dari jenis produk yang dikeringkan. Untuk menjamin kontinuitas operasi, pengering ini juga dapat mengandalkan energi biomassa sebagai salah satu sumber energi termalnya. Alat pengering yang menggunakan dua sumber energi termal ini disebut sebagai pengering ERK-hybrid. Pengembangan ERK-Hybrid menjadi alternatif yang sangat baik dalam hal peralihan penggunaan sumber energi, disamping itu usaha konservasi energi juga penting dilakukan dalam proses pengeringan.
Untuk maksud tersebut, dapat
dilakukan pengeringan dua tahap, yaitu pengeringan dengan laju relatif tinggi kemudian diikuti dengan laju rendah. Menunda atau melakukan pengeringan pada laju rendah sesaat setelah panen merupakan hal yang cukup beresiko. Kadar air tinggi sangat rentan terhadap pertumbuhan mikroba yang dapat menurunkan kualitas biji, sehingga harus diturunkan pada tingkat kadar air tertentu untuk kemudian dapat dilakukan tahap kedua yaitu pengeringan dengan laju relatif lebih rendah. Pemindahan produk ke pengering tahap ke dua dapat menjadi solusi penghematan energi termal, dibandingkan apabila pengeringan dilakukan secara lengkap pada pengering ERK-hybrid. Energi termal yang dibutuhkan untuk pengeringan secara lengkap pada ERK-hybrid lebih besar untuk mendukung laju penurunan kadar air yang tinggi pada awal proses pengeringan, sehingga diperlukan bahan bakar biomassa agar proses dapat berjalan dengan baik. Sistem pengeringan tahap kedua pada umumnya dapat disebut sebagai pengering dalam penyimpan (In-Store Dryer/ISD). Pada kadar air sekitar 18%, biji-bijian termasuk jagung pipilan lebih aman untuk disimpan dalam jangka waktu yang relatif lebih lama pada suhu dan kelembaban umum di Indonesia. Apabila menggunakan asumsi suhu biji-bijian 27 oC umur simpan yang aman pada kadar air 18% dapat lebih dari 20 hari, sedangkan pada kadar air yang lebih tinggi (misalnya 20%) pada suhu yang sama umur aman simpan menjadi hanya kurang dari 10 hari (Brooker et al. 1992). Pada kondisi udara (suhu dan kelembaban) dan kadar air produk tertentu, udara lingkungan (tanpa pemanasan)
3
mempunyai potensi yang sangat besar untuk diterapkan sebagai media pengering berbagai produk bebijian termasuk jagung pipilan. Metode yang digunakan pada ISD umumnya menggunakan udara lingkungan yang dihembuskan melalui tumpukan biji-bijian yang akan dikeringkan. Dengan metode ini, penggunaan pemanas yang membutuhkan perawatan serta biaya operasi lebih tinggi dapat direduksi, selain itu juga dapat dilakukan penghematan energi secara signifikan karena rendahnya kebutuhan energi termal pada operasi, yang biasanya membutuhkan energi cukup tinggi untuk memanaskan udara. ISD sebagai sebuah sistem pengeringan konvektif, mengandalkan aliran udara (gas) yang merupakan kunci utama untuk keberhasilan proses pengeringan karena udara berfungsi sebagai pembawa panas dan uap air. Distribusi aliran udara yang kurang baik dapat menyebabkan ketidakseragaman kadar air dan menyebabkan pula tidak seragamannya kualitas produk, seperti diperlihatkan pada pengeringan rak untuk kakao (Nelwan 1997). Sistem pengeringan tumpukan (deep bed drying) pada ISD akan sangat rentan terhadap permasalahan ketidakseragaman kadar air seperti juga terjadi pada pengeringan tumpukan lainnya, sehingga perlu mekanisme yang baik untuk mengurangi masalah ini. Ketidakseragaman kadar air biasanya dapat diatasi dengan cara pengadukan, namun kebutuhan energi untuk proses pengadukan ini biasanya cukup besar. Hal ini dikemukakan oleh Manalu (1999), yang melakukan percobaan pengadukan dengan menggunakan motor pada pengeringan kakao. Dalam percobaan tersebut dilaporkan bahwa untuk menggerakkan tumpukan kakao sebesar 300-400 kg dibutuhkan motor dengan daya sebesar
1.5 hp.
Kenyataan ini membuat alternatif pengadukan menjadi tidak mungkin dilakukan pada ISD, karena konsep awalnya adalah penghematan energi, sehingga penataan saluran udara di dalam ISD menjadi alternatif yang sangat baik untuk memecahkan permasalahan tersebut. Penataan sistem saluran udara di dalam tumpukan biji-bijian dapat memberikan sebaran kondisi udara yang lebih seragam, karena secara prinsip penataan saluran udara dalam tumpukan biji adalah menyediakan rongga bebas untuk pergerakan udara, sehingga diharapkan udara menjadi bebas bergerak ke segala arah masuk dan keluar tumpukan biji. Sebenarnya konstruksi saluran udara
4
dapat dilakukan secara sederhana dan dengan bahan yang mudah diperoleh, akan tetapi banyaknya kombinasi saluran yang dapat dipilih membuat simulasi matematik menjadi penting untuk menghemat waktu dan biaya disain penataan saluran. Salah satu metode untuk mensimulasikan pola aliran udara, suhu dan tekanan dalam suatu ruang dapat dilakukan dengan metode CFD (Computational Fluid Dynamics). CFD adalah suatu analisis sistem yang meliputi aliran fluida, pindah panas dan fenomena lainnya seperti reaksi kimia yang menggunakan simulasi berbasis komputer. Sementara untuk melihat penyebaran kadar air pada tumpukan biji dalam ISD, dapat dilakukan simulasi dengan menggunakan model pengeringan tumpukan (Brooker et al. 1992), sehingga didapatkan gambaran pengeringan dan perubahan kadar air pada setiap lapisan di dalam tumpukan. Penyimpanan dalam silo besi telah umum digunakan untuk menyimpan produk biji-bijian seperti gabah dan jagung. Namun permasalahan penyimpanan menggunakan silo besi adalah mudahnya terjadi migrasi uap air dalam silo, sehingga kadar air pada bagian tertentu akan naik dan pada bagian lain akan menurun. Hal tersebut disebabkan oleh fluktuasi suhu udara dan RH lingkungan serta radiasi sinar surya. Kadar air yang tinggi pada penyimpanan dapat menurunkan
mutu
mikroorganisme.
karena
akan
Permasalahan
mudah tersebut
terjadinya perlu
perkembangbiakan
dipecahkan
dengan
mengkondisikan silo yang memiliki aerasi udara yang baik dan juga mengurangi efek pemanasan dinding silo akibat radiasi sinar surya. Sebagai sebuah silo, ISD yang dilengkapi dengan pipa-pipa aerasi udara dan dinding dengan insulator panas sehingga bersifat adiabatis, diharapkan mampu melakukan penyimpanan dan mempertahankan mutu produk dengan baik. Penelitian ini merupakan bagian dari penelitian “Rancang Bangun Alat Pengering Efek Rumah Kaca (ERK)-Hybrid dan In-Store Dryer (ISD) Terintegrasi untuk Biji-Bijian”. ISD sebagai satu unit dari sistem terintegrasi tersebut, memerlukan pengkajian secara spesifik untuk melihat performa operasinya secara lebih mendalam yang juga akan menjadi bahan evaluasi bagi pengembangan alat pengering dan penyimpan terintegrasi ini.
5
1.2 Hipotesa Distribusi aliran udara, suhu dan RH mempengaruhi keseragaman kadar air dan mutu biji-bijian di dalam ISD. Dengan simulasi CFD dapat diketahui sebaran aliran udara, suhu dan RH pada ISD, sementara distribusi kadar air dapat diketahui dengan simulasi model pengeringan tumpukan. Berdasarkan kedua hasil simulasi dan validasi terhadap data pengukuran maka akan dapat dijadikan rujukan bagi evaluasi dan pengembangan ISD. 1.3 Tujuan Secara umum tujuan penelitian ini adalah untuk mengkaji distribusi aliran udara, suhu, RH dan kadar air di dalam ISD. Secara khusus penelitian ini terdiri dari beberapa tujuan antara lain : 1. Melakukan simulasi dan validasi model distribusi aliran udara dan suhu di dalam ISD dengan menggunakan teknik CFD. 2. Mendapatkan nilai perhitungan dan validasi RH. 3. Melakukan simulasi dan validasi kadar air jagung menggunakan model pengeringan tumpukan biji-bijian dengan pemrograman Visual Basic 6.0. 4. Analisis mutu jagung hasil pengeringan dan penyimpanan di dalam ISD berdasarkan SNI.
II TINJAUAN PUSTAKA 2.1 Pengeringan Pengeringan adalah proses penurunan kadar air sampai tingkat kadar air tertentu. Secara spesifik pengeringan bertujuan untuk mengurangi kadar air bahan sampai batas dimana perkembangan mikroorganisme dan kegiatan enzim yang dapat menyebabkan kerusakan (fisika/kimia) terhambat atau terhenti, sehingga bahan dapat mempunyai waktu simpan yang lebih lama. Menurut Hall (1980), pada proses pengeringan komoditas pertanian terjadi dua proses dasar yaitu pindah panas untuk menguapkan cairan bahan dan pindah massa akibat adanya perbedaan tekanan uap. Beberapa hal yang harus diperhatikan dalam mengontrol perpindahan kadar air dalam bahan adalah: a) difusi antara cairan dan uap, b) gaya kapilaritas, c) gradien penyusutan dan tekanan uap, d) gravitasi, dan e) penguapan kadar air. Dalam suatu proses pengeringan, dikenal adanya suatu laju pengeringan yang dibedakan menjadi dua tahap utama, yaitu laju pengeringan konstan dan laju pengeringan menurun. Laju pengeringan konstan terjadi pada lapisan air bebas yang terdapat pada permukaan biji-bijian. Laju pengeringan ini terjadi sangat singkat selama proses pengeringan berlangsung, kecepatan penguapan air pada tahap ini dapat disamakan dengan kecepatan penguapan air bebas. Besarnya laju pengeringan ini tergantung dari: a) lapisan yang terbuka, b) perbedaan kelembaban antara aliran udara dan daerah basah, c) koefisien pindah massa, dan e) kecepatan aliran udara pengering. Laju pengeringan menurun terjadi setelah periode pengeringan konstan selesai. Pada tahap ini kecepatan aliran air bebas dari dalam biji ke permukaan lebih kecil dari kecepatan pengambilan uap air maksimum dari biji. Kadar air kritis (critical moisture content) menjadi batas antara laju pengeringan konstan dan laju pangeringan menurun (Hall 1980). Kadar air kritis adalah kadar air terendah pada saat kecepatan aliran air bebas dari dalam biji ke permukaan sama dengan kecepatan pengambilan uap air maksimum dari biji ( Henderson & Perry 1976).
7
2.1.1 Proses Pengeringan Proses pengeringan terjadi dengan cara penguapan air. Cara ini dilakukan dengan menurunkan kelembaban nisbi udara melalui aliran udara panas atau udara bertekanan sehingga tekanan uap air bahan lebih besar dari tekanan uap air udara. Perbedaan tekanan uap ini menyebabkan terjadinya aliran uap air dari bahan ke udara. Proses pengeringan biji-bijian dapat dianggap sebagai proses adiabatik. Selama proses pengeringan berlangsung, entalpi dan suhu bola basah udara pengering tetap, sedangkan suhu bola kering berkurang yang diikuti dengan kenaikan kelembaban mutlak, kelembaban nisbi, tekanan parsial uap air dan suhu pengembunan udara pengering. Terjadinya proses pengeringan dengan udara pengering yang dipanaskan pada kurva psikometrik dapat dilihat pada Gambar 1. Pengeringan dengan menggunakan udara alami berarti proses pemanasan udara (1)-(2) ditiadakan. Kenaikan suhu udara alami karena gesekan atau turbulensi udara dapat dianggap sebagai proses pemanasan udara sebelum masuk ruang pengering.
Gambar 1 Proses pengeringan pada kurva psychrometric Keterangan : (1)-(2) : Proses pemanasan udara (2)-(3) : Proses pengeringan i : udara masuk alat pengering p : udara pengering o : udara keluar dari alat pengering
8
Kelembaban relatif (RH) yang dinyatakan dalam persen merupakan perbandingan antara tekanan uap terhadap tekanan jenuh air pada suhu ruang pengering, yang dinyatakan dalam persamaan (Brooker et al. 1974): ................................................................................................ (1) sedangkan kelembaban mutlak (H) konstan, maka : .
....................................................................................... (2)
dimana 255.38 ≤ T ≤ 533.16 oK dan Pv < Patm, sehingga tekanan uap (Pv) juga konstan. Bila kelembaban udara lingkungan (RHa) dan kelembaban udara pengering (RHr), maka : ............................................................................................ (3) ................................................................ (4) dimana 273.16 ≤ T ≤ 533.16 oK (Keenan & Keyes 1936 dalam ASAE Standard 1994), dimana : R = 22105649.25 A = -27405.526 B = 97.5413 C = -0.146244
D = 0.12558 x 10-3 E = -0.48502 x 10-7 F = 4.34903 G = 0.39381 x 10 -2
2.1.2 Paramater Pengeringan Menurut Brooker et al. (1974), beberapa parameter yang mempengaruhi lama waktu yang dibutuhkan pada proses pengeringan antara lain: a. Suhu udara pengering Suhu udara pengering akan mempengaruhi laju penguapan air bahan dan mutu pengeringan. Semakin tinggi suhu maka panas yang digunakan untuk penguapan air akan meningkat sehingga waktu pengeringan akan menjadi lebih singkat. Agar bahan yang dikeringkan tidak sampai rusak, suhu harus dikontrol terus menerus. b. Kelembaban relatif (RH) udara pengering Kelembaban relatif menentukan kemampuan udara pengering untuk menampung kadar air bahan yang telah diuapkan. Jika RH semakin rendah maka
9
semakin banyak uap air yang diserap udara pengering, demikian juga sebaliknya. RH dan suhu pengering akan menentukan tekanan uap jenuh. Perbedaan tekanan uap air pada udara pengering dan permukaan bahan akan mempengaruhi laju pengeringan. Untuk proses pengeringan yang baik diperlukan RH yang rendah sesuai dengan kondisi bahan yang akan dikeringkan. c. Kecepatan aliran udara pengering Aliran udara pada proses pengeringan berfungsi membawa panas untuk menguapkan kadar air bahan serta mengeluarkan uap air hasil penguapan tersebut. Uap air hasil penguapan bahan dengan panas harus segera dikeluarkan agar tidak membuat jenuh udara pada permukaan bahan, yang akan mengganggu proses pengeringan. Semakin besar volume udara yang mengalir maka akan semakin besar kemampuannya dalam membawa dan menampung air dari permukaan bahan. d. Kadar air bahan Keragaman kadar air awal bahan sering dijumpai pada proses pengeringan dan hal ini juga menjadi suatu masalah. Beberapa hal yang dapat dilakukan untuk mengurangi masalah ini adalah dengan mengurangi ketebalan tumpukan bahan yang dikeringkan, mempercepat aliran udara pengering, menurunkan suhu udara pengering dan dilakukan pengadukan bahan. Kadar air akhir bahan merupakan tujuan akhir proses pengeringan, besarnya kadar air akhir akan menentukan lamanya proses pengeringan berlangsung. Menurut Brooker et al. (1974), Kadar air dapat dinyatakan dalam dua cara, yaitu kadar air basis basah (Mw) dan kadar air basis kering (M). Untuk dipasarkan biasanya kadar air biji-bijian ditentukan berdasarkan basis basah, sementara kadar air basis kering sering digunakan dalam perhitungan-perhitungan engineering. Untuk menghitung kadar air biji-bijian digunakan Persamaan (5) dan (6). .....................................................................................(5) ........................................................................................(6) Pada proses pengeringan sering dijumpai adanya variasi kadar air dari bijibijian yang dikeringkan. Variasi kadar air ini dipengaruhi oleh ketebalan tumpukan biji-bijian, kelembaban nisbi udara pengering, dan kadar air biji-bijian
10
itu sendiri. Brooker et al. (1974) mengemukakan bahwa variasi kadar air bijibijian yang dikeringkan dapat dikurangi dengan cara (1) menipiskan tumpukan biji-bijian, (2) menggunakan kecepatan aliran udara tinggi, (3) mempertahankan suhu udara pengering tetap rendah, dan (4) melakukan pengadukan. Kerusakan fisik dan kimia biji-bijian dapat terjadi akibat pengeringan pada suhu udara pengering yang melebihi batas suhu udara pengering yang diizinkan untuk setiap jenis biji-bijian seperti ditunjukkan pada Tabel 1. Tabel 1 Suhu udara pengering beberapa jenis biji-bijian menurut tujuan penggunaannya No
Jenis biji-bijian
1 2 3 4 5 6 7 8
Tongkol jagung Biji jagung Wheat Oats Barley Butir sorgum Kacang kedelai Padi
Suhu udara pengering maksimum (oC) Benih Dipasarkan Makanan ternak 43.3 54.4 82.2 43.3 54.4 82.2 43.3 60.0 82.2 43.3 60.0 82.2 40.6 40.6 82.2 43.3 60.0 82.2 43.3 48.9 43.3 43.3 -
Sumber : Hall (1970) 2.1.3 Aliran Udara Pengeringan Pada proses pengeringan, udara berfungsi sebagai pendistribusi panas untuk menguapkan kandungan air dari biji-bijian dan mengeluarkan uap air tersebut. Menurut Soemartono (1968), suhu udara dan kecepatan aliran udara pengering berpengaruh penting terhadap proses pengeringan. Air yang dikeluarkan dalam bentuk uap harus segera dipindahkan dan dijauhkan dari biji-bijian sehingga tidak menyebabkan udara jenuh pada permukaan biji-bijian yang dapat memperlambat pengeluaran air selanjutnya. Aliran udara yang cepat akan membawa uap air dari permukaan biji-bijian dan mencegah penjenuhan udara disekitar permukaan bijibijian. Volume udara yang lebih besar dapat menampung dan membawa uap air lebih banyak.
Semakin kering udara maka akan semakin cepat pula proses
pengeringan yang terjadi. Udara kering dapat menampung uap air lebih banyak dari pada udara lembab. Tekanan statik aliran udara pengering yang melalui tumpukan biji-bijian akan memiliki nilai yang berbeda pada saat udara pengering masuk dan keluar tumpukan biji-bijian. Perbedaan tekanan statik ini disebabkan
11
oleh adanya gesekan antara udara pengering dengan biji-bijian dan pengaruh turbulensi aliran udara pengering. Brooker et al. (1974) mengemukakan bahwa tekanan statik aliran udara pengering yang melalui tumpukan bebijian tergantung pada: (a) kecepatan aliran udara pengering, (b) karakteristik bentuk dan permukaan bebijian, (c) jumlah, ukuran dan konfigurasi ruang antar bebijian, (d) variasi ukuran bebijian dan (e) tebal tumpukan bebijian. Faktor lain yang dapat mempengaruhi tekanan statik aliran udara pengering adalah prosentase lubang lantai ruang pengering dan panjang pipa penyalur udara pengering (Hall & Davis 1979). Kebutuhan volume aliran udara pengering untuk biji-bijian menurut cara pengeringan dapat dilihat pada Tabel 2. Tabel 2 Kebutuhan volume aliran udara pengering pada berbagai cara pengeringan. Cara pengeringan Aerasi Tempering Udara pengering tanpa pemanasan Tumpukan tipis Udara pengering dengan pemanasan
Volume aliran udara (m3/m3 det) 2.67 x 10-4 0.0067 0.0267 0.0267 0.4005
Sumber : Brooker et al. (1974) 2.1.4 Karakteristik Pengeringan Jagung Brooker et al. (1992) mengemukakan suatu persamaan untuk konstanta pengeringan jagung yang diambil dari persamaan Pabis dan Henderson (1961) yaitu: 0.54
.
.....................................................................(7)
dimana k dalam (dtk-1) dan T dalam (oR). Kadar air keseimbangan (equilibrium moisture content) merupakan kadar air minimum yang dapat dicapai pada kondisi udara pengeringan yang tetap atau pada suhu dan kelembaban relatif yang tetap. Suatu bahan dalam keadaan setimbang apabila laju kehilangan air dari bahan ke udara sekelilingnya sama dengan laju penambahan air ke bahan dari udara di sekelilingnya. Kadar air pada keadaan setimbang disebut juga dengan kadar air keseimbangan atau keseimbangan higroskopis (Henderson & Perry 1979).
12
Salah satu persamaan kadar air keseimbangan pada jagung pipilan adalah persamaan Henderson termodifikasi (Brooker et al. 1992) yaitu: 1
exp
8.6541 10
49.810 100 .
.
.
.
.
.......... (8)
............................................. (9)
dimana Me adalah kadar air keseimbangan basis kering, T adalah suhu mutlak udara (oC) dan RH adalah kelembaban nisbi.
Apabila persamaan di atas
digunakan untuk kondisi udara alami yang umum di Indonesia, sebagai contoh pada suhu 30 oC dan RH 70%, maka nilai kadar air keseimbangan jagung pipilan yang diperoleh adalah 16.0% b.k. Tabel 3 menyajikan beberapa nilai yang diuji pada studi pendahuluan. Dari tabel tersebut, dapat dilihat bahwa udara alami tanpa pemanasan mempunyai potensi untuk menurunkan kadar air jagung pipilan sampai 13.8%b.k. Kadar air ini sudah memadai untuk penyimpanan dalam waktu yang cukup lama. Tabel 3 Hasil perhitungan kadar air berdasarkan Persamaan (8) T (oC) 25 30
RH (%) 60 70 80 60 70 80
Me (%b.k) 14.3 16.5 19.3 13.8 16.0 18.8
Penggunaan udara lingkungan tanpa pemanasan sebagai udara pengering telah diuji di Korea selama empat tahun (Kim et al. 1989).
Kondisi udara
lingkungan yang digunakan mempunyai suhu udara rata-rata 12.8-18.6 oC dengan RH rata-rata berkisar antara 63.3-72.0%. Dengan kondisi tersebut sebanyak 25003000 kg gabah dapat dikeringkan dari kadar air awal 17.2-21.9% sampai kadar air akhir 13.2-14.6%. 2.1.5 Sorpsi Isotermi Sorpsi isotermi adalah suatu plot kadar air keseimbangan terhadap kelembaban relatif pada suatu temperatur tertentu. Isotermi yang diperoleh dengan memaparkan padatan pada udara yang kelembabannya meningkat dikenal dengan isotermi adsorpsi, sedangkan isotermi yang diperoleh dengan memaparkan
13
p padatan padda udara yang kelembbabannya menurun m dikkenal dengan isotermi d desorpsi. Isootermi desorrpsi merupakkan perhatian n utama penngeringan kaarena kadar a padatan menurun secara air s proggresif. Kebaanyakan bahhan yang dikeringkan d m menunjukka an hysterisiss dimana keedua isoterm mi tersebut tidak samaa sebangun ( (Devahastin 2000). m isotermi ssorpsi tipikaal. Bentuk Gambaar 2 menunnjukkan benntuk umum t tersebut diciirikan oleh tiga wilayahh secara teggas, A, B daan C, yang merupakan m p pertanda meekanisme peengikatan aiir yang berbbeda pada tempat t terpiisah dalam m matrik padaatan. Pada wilayah w A, air terikat ku uat pada tem mpat tersebutt dan tidak d dapat digunaakan untuk reaksi. Padaa wilayah inni, terutama terdapat adssorpsi lapis t tunggal uap air dan tid dak tampak pperbedaan tegas t antara isotermi addsorpsi dan d desorpsi. Paada wilayah B, air terikkat lebih lonnggar. Penuurunan tekannan uap air h hingga di bawah b tekannan keseimbbangan uap air pada suuhu yang saama adalah k karena air teersebut terkuurung dalam kapiler yang g lebih kecill. Air dalam m wilayah C b bahkan terikkat lebih lo onggar dalaam kapiler yang lebihh besar, air ini dapat d digunakan u untuk reaksi dan sebagai pelarut (Devvahastin 20000).
Gambar 2 Sorpsi isoteermi yang meenunjukkan hysterisis 2 2.1.6 Aktivitas Air hastin (2000)) mengemukkakan bahwaa dalam penggeringan bahhan pangan Devah a atau pakan, ketersediaan n air untuk ppertumbuhann mikroorgannisme, perkeecambahan s spora dan koontribusi dallam beberapa reaksi kim mia, memerluukan perhatiaan penting. H ini dikkarena aktiivitas tersebbut akan mengakibatk Hal m kan kerusakkan bahan. A Aktivitas air didefinisikkan sebagai perbandinggan antara ttekanan parssial air (p)
14
pada sistem padatan basah terhadap tekanan keseimbangan uap air (pw) pada suhu yang sama, dalam persamaan dituliskan sebagai : ................................................................................................ (10) atau, .............................................................................. (11) Daftar nilai aw minimum terukur untuk pertumbuhan mikroba dan perkecambahan spora disajikan pada Tabel 4. Jika aw diturunkan dibawah nilai ini dengan cara pengeringan atau dengan menambahkan agen pengikat seperti gula, gliserol atau garam, maka pertumbuhan mikroba dapat dihambat. Akan tetapi seharusnya penambahan tersebut tidak mempengaruhi aroma, rasa atau kriteria mutu lainnya, sehingga proses pengeringan merupakan solusi yang baik untuk menurunkan aw pada bahan pangan dengan kadar air tinggi. Tabel
4 Aktivitas air (aw) minimum untuk pertumbuhan mikroba dan perkecambahan spora
Mikroorganisme Organisme penghasil lendir pada daging Spora Pseudomonas, Bacillus cereus Spora B.subtilis, C.botulinum C.Botulinum, Salmonela Bakteri pada umumnya Ragi pada umumnya Aspergillus niger Jamur pada umumnya Bakteri halofolik Jamur Xerofilik Ragi Osmifilik
Aktivitas air 0.98 0.97 0.95 0.93 0.91 0.88 0.85 0.80 0.75 0.65 0.62
Sumber : Brockmann 1973 dalam Devahastin 2000 2.2 Penyimpanan Penyimpanan adalah suatu cara pengamanan yang selalu berkaitan dengan waktu. Hasil pertanian terutama bebijian selama penyimpanan masih mengalami proses respirasi karena bahan tersebut masih hidup. Proses respirasi merupakan proses produksi energi yang digunakan oleh sel-sel tanaman, pada proses respirasi terjadi pemindahan energi dari ikatan kimia dalam bahan kepada ikatan kimia Adenosin Tri Phospat (ATP) yang berenergi tinggi dan langsung digunakan dalam proses kehidupan (Suseno 1974). Menurut Hall (1970), air dan panas yang dihasilkan dari proses respirasi akan menaikkan kadar air bahan dan suhu,
15
sehingga laju respirasi meningkat. Kadar air dan panas hasil respirasi membuat kondisi yang baik bagi pertumbuhan kapang. Wijandi (1988) mengemukakan bahwa penyimpanan dimaksudkan untuk menjaga dan mempertahankan mutu komoditi
yang
disimpan
dengan
jalan
menghindari,
mengurangi
atau
menghilangkan berbagai faktor yang dapat mengurangi mutu komoditi yang disimpan. Menurut Soesarsono (1977), penyimpanan dapat dibagi dalam berbagai tahapan/kelompok antara lain: berdasarkan perjalanan hasil panen, waktu, tempat, modifikasi udara dan berdasarkan teknologi. Dalam penyimpanan berdasarkan perjalanan hasil panen, dikenal penyimpanan tingkat panen, tingkat petani, tingkat pengumpul, tingkat penyalur, transit, tingkat pengecer dan tingkat konsumen. Berdasarkan waktu dilakukan penyimpanan jangka panjang, jangka menengah, jangka pendek, transit dan penyimpanan panjang. Penyimpanan berdasarkan tempat digolongkan menjadi penyimpanan di atas atmosfer, di dalam tanah, di udara dan di bawah permukaan air. Berdasarkan modifikasi udara dikenal penyimpanan alami, penyimpanan atmosfir yang dimodifikasi (modified atmosfer storage) dan penyimpanan atmosfr yang dikendalikan (control atmosfer storage). Sedangkan berdasarkan teknologi, penyimpanan dapat digolongkan menjadi penyimpanan tradisional dan penyimpanan modern. Cara penyimpanan modern merupakan pengembangan dari penyimpanan tradisional. Menurut Wiliam (1991), ada beberapa faktor yang berpengaruh pada penyimpanan biji-bijian antara lain: tipe dari bebijian, periode penyimpanan, metode penyimpanan, suhu lingkungan, kadar air bahan, kandungan bahan asing, proteksi fisik dan kelembaban relatif. Jagung dapat disimpan dalam beberapa cara seperti curah (pipilan), kemas (pipilan) dan gantung (dengan tongkol). Berdasarkan pengaruh udara lingkungan pada kondisi penyimpanan, penyimpanan dapat dibedakan menjadi penyimpanan udara bebas dan penyimpanan rapat udara (Thahir et al. 1988). Penyimpanan udara bebas adalah penyimpanan yang dilakukan pada kondisi udara bebas dengan suhu kamar, pada kondisi ini lingkungan berpengaruh langsung terhadap proses penyimpanan. Sistem penyimpanan udara bebas kurang menguntungkan bagi biji dengan kadar air awal rendah pada daerah dengan
16
kelembaban yang tinggi, karena kadar air biji akan naik menyesuaikan dengan kelembaban udara lingkungan. Kerusakan akan tetap terjadi meskipun telah diterapkan persyaratan penyimpanan yang cukup baik (Thahir et al. 1988). Penyimpanan rapat udara merupakan sistem penyimpanan dengan prinsip membatasi dampak negatif dari udara lingkungan sehingga laju kerusakan dapat dihambat. Penyimpanan ini juga sering disebut penyimpanan kedap udara. Tujuan dari sistem penyimpanan tersebut adalah untuk memperpanjang daya simpan. Kerusakan butir bijian terjadi karena kegiatan biologis hama, kapang dan bakteri. Kegiatan biologis berupa pernafasan dapat dihambat dengan cara kemasan diisi biji penuh, kadar air butiran rendah pada awal penyimpanan, digunakan wadah dengan sistem kedap udara. Tingkat pernafasan dapat dihambat dengan cara pemberian CO2, pengurangan O2. Keuntungan dari sistem penyimpanan ini memperpanjang daya simpan jagung dari tiga bulan menjadi paling sedikit 12 bulan, serta tidak membutuhkan insektisida dan fungisida, yang diperlukan adalah kadar air yang rendah (Thahir et al. 1988). 2.2.1 Pengaruh Kadar Air terhadap Penyimpanan Dalam mencegah kerusakan selama masa penyimpanan, pengendalian kadar air merupakan faktor terpenting. Pengendalian kadar air adalah faktor yang paling mudah
dan
murah
sebelum
dilakukan
penyimpanan
terhadap
bahan.
Perkembangan kapang dapat ditekan dengan adanya pengurangan kadar air selama penyimpanan (Wiliam 1991). Pengeringan yang berlanjut dengan menggunakan sinar matahari dapat menyebabkan biji-bijian retak dan kehilangan daya hidupnya (Covanic 1991 dalam Dharmaputra et al. 1997). Selama masa penyimpanan kadar air bahan pangan akan bergerak menuju kadar air keseimbangan. Henderson dan Perry (1976) mengemukakan bahwa kadar air keseimbangan terjadi pada saat biji-bijian tidak lagi menyerap atau melepaskan uap air. Pengeringan mekanis untuk menurunkan kadar air sampai 14% selama 2.5 hari efektif untuk mengontrol aflatoksin pada jagung yang diproduksi pada musim hujan. Untuk menghemat biaya, pengeringan mekanis dipadukan dengan metode pengeringan field drying, yaitu dengan membiarkan jagung tetap di pohon meskipun jagung tersebut telah siap panen. Perlakuan ini dimaksudkan untuk
17
mencapai titik keseimbangan dengan kelembaban relatif. Di Thailand selama musim hujan menggunakan field drying selama 1 sampai 4 minggu efektif mengeringkan jagung dari kadar air >26% menjadi 18 – 22%, juga menjaga dari kerusakan fisik serta mengontrol aflatoksin (Negler et al. 1986). 2.2.2 Kelembaban dan Suhu Penyimpanan Chikubu (1974) mengemukakan bahwa kelembaban dan suhu ruang merupakan faktor lingkungan yang penting dalam penyimpanan. Kelembaban lebih berperan dalam menentukan mutu bahan dan proses kerusakan selama penyimpanan. Kelembaban akan mempengaruhi kadar air bahan, dan kadar air bahan juga selalu dipengaruhi oleh kelembaban ruangan, sehingga terjadi suatu keseimbangan. Selain itu suhu ruangan juga sangat menentukan tingkat keseimbangan kelembaban dengan kadar air tersebut. Batas suhu dan kadar air aman pada penyimpanan biji-bijian dapat dilihat pada Gambar 3.
Gambar 3 Batas-batas suhu dan kadar air yang aman pada penyimpanan biji-bijian (Hall 1970) Kelembaban ruangan, suhu dan kadar air bahan selain mempengaruhi aktifitas di dalam bahan juga akan mempengaruhi kegiatan hidup organisme perusak. Setiap organisme perusak memerlukan syarat hidup tertentu sehingga dapat tumbuh dan berkembang biak dengan baik. Soesarsono (1977) mengemukakan bahwa pada umumnya aktivitas serangga tidak dipengaruhi oleh
18
kelembaban, tetapi dipengaruhi oleh perubahan suhu. Sedangkan jamur dipengaruhi oleh kadar air dan relatif tidak terpengaruh oleh suhu. Menurut Hall (1970), kondisi yang sesuai untuk mencegah kerusakan selama penyimpanan dan perdagangan adalah pada kelembaban sebesar 70%. Pada kelembaban 70% dan suhu 27 oC jagung memiliki kadar air keseimbangan 13.5%. Selanjutnya dijelaskan bahwa suhu yang tinggi (berkisar antara 21-43 oC) akan mempercepat kehidupan organisme, disamping itu reaksi kimia juga akan meningkat karena peningkatan suhu. Kenaikan suhu bahan juga disebabkan oleh kegiatan respirasi, aktifitas serangga, aktifitas kapang dan bakteri. 2.2.3 Hubungan Antara Penyimpanan dan Kerusakan Bahan Pakan Menurut Francis dan Wood (1982), kondisi lingkungan yang berpengaruh pada penyimpanan adalah suhu dan kelembaban relatif, dan hanya terpengaruh kecil oleh oksigen dan cahaya. Suhu dan kelembaban relatif tidak hanya berpengaruh terhadap laju perubahan kimia tapi juga berpengaruh pada perkembangan serangga dan kapang. Perubahan kimia berhubungan erat dengan aktivitas kapang dan serangga. Untuk lebih jelas dapat dilihat pada Tabel 5. Tabel 5 Perubahan biologi dan kimia pada pakan konsentrat Kadar air (%)
RH pada 20-30 oC (%)
Aktivitas Biologis
<8
30
Tidak nyata
8 – 14
30 – 70 >60
Serangan serangga Serangan serangga
Peningkatan asam urat Reaksi Maillard
12 – 20
70 – 90
Serangan serangga Pertumbuhan kapang
Produksi mikotoksin
20 – 25
90 – 95
Serangan serangga Pertumbuhan kapang
Peningkatan produksi mikotoksin
> 25
-
Pertumbuhan bakteri
Kehilangan fisik dan depolimerisasi pati dan protein
Aktivitas Kimia oksidasi lemak, peningkatan peroksida
Sumber : Francis dan Wood (1982) 2.2.4 Persyaratan Mutu Jagung Persyaratan mutu jagung untuk perdagangan menurut SNI dikelompokkan menjadi dua bagian yaitu persyaratan kualitatif dan persyaratan kuantitatif (Kristanto 2007). Persyaratan kualitatif jagung meliputi: 1. Produk harus terbebas dari hama dan penyakit.
19
2. Produk terbebas dari bau busuk maupun zat kimia lainnya (berupa asam). 3. Produk harus terbebas dari bahan dan sisa-sisa pupuk maupun pestisida. 4. Memiliki suhu normal. Sedangkan persyaratan kuantitatif jagung dapat dilihat pada Tabel 6. Tabel 6 Persyaratan mutu jagung No Komponen Utama 1 2 3 4 5
Kadar air Butir rusak Butir warna lain Butir pecah Kotoran
Persyaratan Mutu (% Maks) I II III IV 14 14 15 17 2 4 6 8 1 3 7 10 1 4 3 5 1 1 2 2
Sumber: SNI 01-03920-1995 dalam Kristanto, 2007 Standar Mutu jagung yang digunakan untuk bahan baku pakan meliputi zat makanan dan kandungan bahan berbahaya/racun serta kemurnian, standar tersebut dapat dilihat pada Tabel 7. Tabel 7 Standar mutu jagung bahan baku pakan ternak No 1 2 3 4 5 6 7 8 9 10
Komponen Kadar air (maksimum) % Kadar protein kasar (minimum) % Kadar serat kasar (maksimum) % Kadar abu (maksimum) % Kadar lemak (minimum) % Mikotoksin: a) Aflatoksin (maksimum) ppb b) Okratoksin (maksimum) ppb Butir pecah (maksimum) % Warna lain (maksimum) % Benda asing (maksimum) % Kepadatan minimum kg/m3
Persyaratan 14.0 7.5 3.0 2.0 3.0 50.0 5.0 5.0 5.0 2.0 700
Sumber: SNI 01-4483-1998 Untuk dijadikan benih, biji jagung hibrida harus memenuhi persyaratan sebagai berikut: Tabel 8 Spesifikasi persyaratan mutu benih jagung hibrida di laboratorium No 1 2 3 5
Komponen Kadar air (maksimum) Benih murni (minimum) Daya berkecambah (minimum) Kotoran Benih (maksimum)
Sumber: SNI 01-6944-2003
Persyaratan (%) 12.0 98.0 85.0 2.0
20
2.3 Perkembangan Penelitian In-Store Dryer Gagasan yang mutakhir mengenai proses penyimpanan yang disatukan dengan pengeringan telah banyak dilakukan dengan berbagai bentuk bangunan maupun metode pengeringan dan penyimpanannya. Di beberapa negara ASEAN yang beriklim tropis dan sub-tropis, telah berkembang penelitian serta percobaan untuk mengetahui sejauh mana sistem penyatuan proses pengeringan dan penyimpanan dapat mengurangi susut bahan pascapanen. Koto (1983) telah melakukan penelitian mengenai penyimpanan dalam silo besi kedap udara. Penelitian menggunakan gabah sebagai bahan uji ini bertujuan untuk melihat perubahan kadar air selama penyimpanan akibat pengaruh fluktuasi suhu udara dan radiasi sinar surya. Percobaan menggunakan gabah varietas bolon, yang disimpan dalam silo besi dengan diameter 150 cm dan tinggi 100 cm, yang diletakkan pada udara terbuka sehingga dinding silo dapat terkena sinar matahari langsung. Hasil pengamatan selama 100 hari penyimpanan menunjukkan adanya perbedaan kadar air antara hasil perhitungan sebesar 0.62% dan hasil pengamatan sebesar 0.58%. Pada lokasi pusat lapisan bawah silo dan di lokasi sepanjang 5 cm dari dinding terjadi penurunan kadar air sebesar 3%. Selama penyimpanan tidak terjadi perubahan warna beras, sementara peningkatan butir retak paling banyak terjadi di sekitar dinding dan paling sedikit pada lokasi pusat silo. Peningkatan populasi kapang sangat kecil dan pengaruhnya tidak nyata terhadap mutu beras. Beberapa percobaan modifikasi lumbung pengering telah dibuat oleh Universitas Gadjah Mada di Yogyakarta dan Balai Besar Pengembangan Mekanisasi Pertanian di Serpong. Soemangat et al. (1987) melakukan studi implementasi pengering tipe lahat dalam tanah untuk jagung. Karakteristik utama alat ini adalah; (a) berukuran 4.6 m x 2.1 m x 1.8 m, terdiri atas ruang piramida, tungku, plenum, cerobong dan atap, (b) kapasitas pengering 1 ton, (c) waktu pengeringan 13 jam pada suhu 68 oC untuk menurunkan kadar air dari 35% menjadi 17% basis basah, (d) beroperasi pada malam hari dan musim hujan pada suhu lingkungan 24 oC dan RH 96%. Pengering ini pertama kali dikembangkan oleh SUCA (Silliman University College of Agriculture) pada tahun 1984 di Afrika.
21
Lumbung pengering bahan bakar non-konvensional IRRI (Harlos et al. 1983) dikembangkan lebih lanjut oleh Jeon et al. (1983), dengan kapasitas pengeringan 8 ton/proses dapat mengeringkan gabah dari kadar air 20% menjadi 14% basis basah selama 6-12 jam pada suhu 39-42 oC dan laju hisapan udara sebesar 9.83 m3/ m3 /mnt. Komar (1988) meneliti sebuah alat penyimpan sekaligus pengering berupa sebuah sistem lumbung pengering gabah bahan bakar sekam. Penelitian ini menghasilkan suatu bangunan lumbung berukuran 3 m x 2 m x 3 m yang dapat menghasilkan antara lain: suhu udara panas 35-40 oC, RH 55-58.72%, laju aerasi 3.09 x 10-3 kg/dtk selama pengeringan dan suhu udara ruangan 30-33 oC, RH 56.89-60% dengan laju aerasi 1.12 x 10-4 kg/dtk selama penyimpanan. Lumbung dengan muatan 500 kg gabah ini dapat menurunkan kadar airnya dari 27.63% basis kering menjadi 15% b.k dalam jangka waktu 36 jam, dengan konsumsi bahan bakar selama pengeringan adalah 3.3 kg/jam, efisiensi panas tungku yang digunakan adalah 60%. Percobaan penyimpanan gabah dalam lumbung selama dua bulan, dengan memanfaatkan panas surya yang dipindahkan melalui atap seng gelombang ke ruang lumbung untuk menurunkan dan mempertahankan kadar air gabah. Penyimpanan ini menghasilkan indeks kerusakan antara 1-5 dan susut bahan kering antara 1-1.5, dari nilai indeks tersebut lumbung dapat digunakan untuk penyimpanan gabah jangka panjang. Widodo et al. (1994) di BBP MEKTAN Serpong melakukan analisis teknis dan ekonomis pada pengering padi dengan menggunakan Drying and Storage System (DS System). DS Sytem tersebut terdiri dari 8 kotak pengering, motor penggerak dan kipas penghembus. Pengering ini mampu mengeringkan gabah dari kadar air 27.75-28.97% b.b menjadi 11.4% b.b dalam waktu 16 jam, dengan laju pengeringan 0.98% per jam. Kelima model pengering dan penyimpan tersebut masih memanfaatkan mekanisme pindah panas konveksi alami dan aliran udara dengan sumber panas dari bahan bakar pada tungku maupun sinar matahari.
22
2.4 Teknik Simulasi Computational Fluid Dynamics (CFD) Computational Fluid Dynamics (CFD) adalah suatu analisis sistem yang meliputi aliran fluida, perpindahan panas dan fenomena lain seperti reaksi kimia yang menggunakan simulasi dengan bantuan software komputer. CFD telah dikenal sejak tahun 1960-an untuk mendisain mesin jet dan aircraft. Perkembangan selanjutnya metoda ini digunakan
untuk mendisain mesin
pembakaran internal, tabung pembakaran dalam turbin gas dan tungku, kendaraan bermotor dan aliran udara yang menyelimuti casing mobil. Metode CFD menggunakan analisa numerik yaitu kontrol volume sebagai elemen dari integrasi persamaan-persamaan yang terdiri dari persamaan keseimbangan massa, momentum dan energi (Versteeg & Malalasekera 1995). Wulandani (2005) telah menggunakan teknik CFD untuk mensimulasi udara pengering pada pengeringan ERK tipe rak. Dalam riset ini dilakukan analisis distribusi aliran udara yang mencakup kecepatan, suhu dan RH serta dilanjutkan dengan melakukan validasi model tersebut terhadap hasil percobaan. Analisis ini penting untuk mengoptimisasikan bentuk saluran udara yang harus didisain untuk menyeragamkan aliran udara pada pengering, sehingga diperoleh keseragaman kadar air yang berarti juga keseragaman kualitas biji. Dalam CFD, pola aliran udara digambarkan secara kuantitatif dalam besaran suhu dan kecepatan melalui persamaan diferensial berupa koordinat cartesian. Pemecahan secara matematik dalam CFD dilakukan melalui analisis numerik tiga dimensi dengan metode volume hingga melalui diskretisasi dan iterasi. Analisis distribusi dan simulasi suhu dan kecepatan udara pada ruangan ISD dalam CFD dapat dilakukan dengan menggunakan software gambit 2.2.30 (meshing dan boundary condition) dan fluent 6.1.18 (mendefinisikan model 3D, pemakaian energi, viscous model, jenis material dan sifat termofisik fluida, input nilai boundary condition, inisialisasi, iterasi dan visualisasi). Computational Fluid Dynamics (CFD) mengandung 3 komponen utama, yaitu : pre-processor, solver dan post-processror (Versteeg & Malalasekera 1995). 2.4.1 Pre-processor Komponen pre-processor merupakan komponen input dari permasalahan aliran ke dalam program CFD dengan menggunakan interface yang memudahkan
23
operator, berfungsi sebagai transformer input berikutnya ke dalam bentuk yang sesuai dengan pemecahan oleh solver. Pada tahapan pre-processor, dapat dilakukan hal-hal sebagai berikut: 1) mendefinisikan geometri daerah yang dikehendaki (perhitungan domain); 2) pembentukan grid (mesh) pada setiap domain; 3) pemilihan fenomena kimia dan fisik yang dibutuhkan; 4) menentukan sifat-sifat fluida (konduktivitas, viskositas, panas jenis, massa jenis dan sebagainya); 5) menentukan kondisi batas yang sesuai dengan keperluan. Ketepatan aliran dalam geometri yang dibentuk dalam CFD ditentukan oleh jumlah sel di dalam grid yang dibangun. Semakin besar jumlah sel, ketepatan atau ketelitian dari hasil pemecahan semakin baik. Mesh optimal tidak harus selalu seragam, dapat dilakukan dengan memperhalus mesh pada bagian yang memiliki variasi cukup besar dan semakin kasar untuk bagian yang relatif tidak banyak mengalami perubahan. 2.4.2 Solver Proses pada solver merupakan proses pemecahan secara matematika dalam CFD dengan software fluent 6.1.18 Metode yang digunakan adalah metode volume hingga (finite volume) yang dikembangkan dari metode beda hingga (finite difference) khusus. Proses pemecahan matematika pada solver digambarkan sebagai diagram alir metode SIMPLE (Semi-Implicit Method for Pressure-Linked Equation) (Lampiran 6). Proses pemecahan matematika pada solver memiliki 3 tahapan yaitu: 1) aproksimasi aliran yang tidak diketahui dilakukan dengan menggunakan fungsi sederhana; 2) diskretisasi dengan mensubstitusi hasil aproksimasi ke dalam persamaan aliran disertai dengan manipulasi matematis; 3) penyelesaian persamaan aljabar. Pada proses solver, terdapat 3 persamaan atur aliran fluida yang menyatakan hukum kekekalan fisika, yaitu : 1) massa fluida kekal; 2) laju perubahan momentum sama dengan resultansi gaya pada partikel fluida (Hukum II Newton); 3) laju perubahan energi sama dengan resultansi laju panas yang ditambahkan dan laju kerja yang diberikan pada partikel fluida (Hukum I Termodinamika).
24
Hukum Kekalan Massa 3 Dimensi Steady State Keseimbangan massa untuk elemen fluida dinyatakan sebagai berikut: laju kenaikan massa dalam elemen fluida = laju net aliran massa ke dalam elemen terbatas. Adapun bentuk matematis dapat ditulis sebagai berikut (Versteeg & Malalasekera 1995) : 0 ................................................................... (12) Persamaan (10) merupakan persamaan kontinyuitas untuk fluida. Ruas kiri menggambarkan laju netto massa keluar dari elemen melewati batas dan dinyatakan sebagai faktor konveksi. Persamaan Momentum 3 Dimensi Steady State Persamaan momentum dikembangkan dari persamaan Navier-Stokes dalam bentuk yang sesuai dengan metode finite volume (Versteeg & Malalasekera 1995) sebagai berikut : Momentum arah x: ................. (13) Momentum arah y : .................. (14) Momentum arah z: ............... (15)
Persamaan Energi 3 Dimensi Steady State Persamaan energi diturunkan dari Hukum I Termodinamika (Versteeg & Malalasekera 1995) yang menyatakan bahwa : laju perubahan energi partikel fluida = laju penambahan panas ke dalam partikel fluida ditambahkan dengan laju kerja yang diberikan pada partikel. Secara matematika dapat ditulis sebagai berikut :
25
....... (16) Persamaan state: Kecepatan fluida selalu mencari keseimbangan secara termodinamika, kecuali adanya gangguan. Jika digunakan variabel ρ dan p, maka persamaan state untuk p dan i (Versteeg & Malalasekera 1995) adalah sebagai berikut : , ,
.......................................................................................... (17) ............................................................................................ (18)
Untuk gas ideal, dimana :
dan
2.4.3 Post-processor Post-processor merupakan tahapan akhir dari simulasi CFD, tahap ini berupa penampilan hasil yang diperoleh dari proses sebelumnya dalam preprocessor dan solver. Tampilan tersebut dapat berupa: a) tampilan geometri domain dan grid; b) plot vektor; c) plot permukaan 2 dan 3 dimensi; d) pergerakan partikel; e) manipulasi pandangan; f) output warna. 2.5 Model Pengeringan Tumpukan (Deep Bed Drying) Brooker et al. (1992) mengemukakan bahwa suatu model pengeringan tumpukan diturunkan berdasarkan keseimbangan panas dan massa. Menurut Sharp (1982) dalam Napitupulu (1993), ada beberapa model pengeringan lapis tebal yang dapat digunakan untuk menggambarkan kondisi pengeringan tumpukan diantaranya adalah: model keseimbangan, model logaritmik dan model persamaan differensial parsial. Nugroho (1986) mengemukakan bahwa hasil simulasi yang didapatkan dengan model differensial parsial lebih mendekati hasil percobaan dibandingkan dengan model keseimbangan. Suatu
model
persamaan
differensial
parsial
merupakan
analisis
keseimbangan kalor dan massa melalui kontrol volume seperti yang ditunjukkan pada Gambar 4.
26
Gambar 4 Elemen pada bak (Bala 1997) Menurut Brooker et al. (1992), untuk mempermudah penyelesaian dari persamaan pengeringan yang ada, beberapa asumsi digunakan untuk menurunkan persamaan model tersebut yaitu: 1) selama proses pengeringan, penciutan volume bahan biji-bijian diabaikan, 2) gradien temperatur dalam setiap bahan diabaikan, 3) konduksi antar biji-bijian diabaikan, 4) aliran udara bersifat plug type dan konstan, 5)
/
dan
/
dapat diabaikan terhadap
/
dan
/
, 6)
dinding alat pengering bersifat adiabatis selama proses pengeringan dan kapasitas panasnya diabaikan, dan 7) kapasitas panas uap air dan biji-bijian dianggap tetap selama pengeringan. Bala (1997), juga mengemukakan beberapa asumsi untuk penurunan persamaan deep bed drying, antara lain: a) aliran udara satu dimensi, b) tidak ada kehilangan panas tegak lurus aliran udara, c) kehilangan panas konduksi di dalam bak diabaikan, d) panas spesifik bijian kering, moisture dan udara konstan, e) panas laten penguapan dipengaruhi oleh moisture content, f) penyusutan bijian di dalam bak dipengaruhi oleh moisture content, g) bulk density bijian dipengaruhi oleh penyusutan, h) kontribusi (dH/dt) dan (dTa/dt) diabaikan. Dari parameter keseimbangan entalpi dan massa pada elemen bak yang ditunjukkan pada Gambar 4, ada empat parameter yang belum ditentukan dalam model ini, yaitu: a) kadar air biji-bijian (M), b) kelembaban mutlak udara (H), c) suhu udara pengering (T) dan d) suhu bahan biji-bijian (Tg). Berdasarkan asumsi-asumsi di atas, Bala (1997) menurunkan persamaan model matematis diferensial parsial tipe bak, meliputi: a) keseimbangan massa, b) laju pengeringan, c) keseimbangan panas dan d) laju perpindahan panas.
27
2.5.1 Keseimbangan Massa Dalam unit waktu tertentu, aliran uap air yang masuk ke dalam elemen bak dituliskan sebagai
dan keluar sebagai
. Selisih antara
keduanya adalah uap air yang bertambah ke udara dari bijian. Sehingga dalam unit waktu tertentu dituliskan sebagai: ...........................................(19) Dengan menggunakan deret Taylor untuk H dan dengan mengabaikan semua batasan dz2, maka Persamaan (19) ditulis menjadi : ....................................................................(20) .............................................................................(21) Persamaan (21) jika ditulis dalam bentuk finite difference adalah: ∆
∆
∆
∆
...............................................................................(22)
dengan mengacu pada grid (Gambar 5) maka : ∆ ∆
∆
..................................................................(23)
Gambar 5 Grid finite different untuk persamaan deep bed drying
28
2.5.2 Laju Pengeringan Kadar air pada suatu lapis tipis biji-bijian sesuai dengan ekspresi persamaan lapis tipis yaitu: .......................................................................(24) integrasi persamaan ini dari step 1 ke 2 pada Gambar 5 menjadi: .......................................................................(25) sehingga bila t2-t1 = Δt, maka: ∆
∆
1
....................................................(26)
Persamaan (26) ditulis dalam page equation (Hall 1970; Van Rest & Isaacs 1968 dalam Bala 1997) maka: ∆
∆
...............................(27)
dimana : /
/
..........................................................(28)
2.5.3 Keseimbangan Panas Perubahan entalpi udara = perpindahan panas konveksi ke bijian – panas yang dibawa oleh uap ke udara. Panas yang masuk ke dalam elemen (z, z+dz) dalam suatu waktu tertentu adalah: ...................................(29) dan panas yang keluar adalah: ...........(30) dari Persamaan (29) dan Persamaan (30)
..............(31)
29
dengan menggunakan deret Taylor, Persamaan (31) ditulis sebagai :
...............(32) subsitusi (dH/dz)= -(ρd/Ga)(dM/dt):
...............................(33) Persamaan (33) menjadi : ..........................................................(34) /
jika
diasumsikan P dan Tg konstan untuk interval Δz (Gambar 5), Persamaan (34) menjadi: ......................................................................(35) integrasi step 2 ke 4 didapatkan: ............................................................................(36) jika z4-z2=Δz , maka: ∆
1
∆
...................................................(37)
2.5.4 Laju Perpindahan Panas Perubahan entalpi bijian = perpindahan panas konveksi ke bijian – panas yang di supply untuk menguapkan uap air ke udara. Pada awal waktu dt, panas bijian adalah : ................................................................(38) dan pada t+dt adalah: .............................................(39)
30
sehingga perubahan entalpi bijian:
..................................(40) dengan menggunakan deret Taylor, Persamaan (40) menjadi :
...............................(41) Persamaan (41) disusun ulang menjadi : .....................................(42) /
Jika ′
/ Diasumsikan P’ dan Q konstan terhadap waktu untuk interval dt, maka: ′
....................................................................(43)
Integrasi step 1 ke 2 (Gambar 5) didapatkan: ′
′
/ ′
′
/ ′ .............(44)
jika t2 - t1 = Δt, maka Persamaan (44) dapat dituliskan : ′∆
′
1
′∆
/ ′ ..................................(45)
Model pengeringan tumpukan bebijian pada prinsipnya disusun dengan membagi tumpukan (lapisan tebal) menjadi beberapa lapisan tipis, keluaran dari lapisan sebelumnya merupakan masukan untuk lapisan berikutnya dan pada setiap lapisan diasumsikan bahwa suhu dan kadar air bebijian dalam kondisi seragam. Model
persamaan
differensial
parsial
biasanya
digunakan
untuk
mempelajari pengeringan pada suhu tinggi. Tidak banyak model persamaan diferensial parsial yang memberikan hasil memuaskan untuk kondisi pengeringan dengan suhu rendah, kecepatan aliran udara rendah dan kedalaman bak cukup besar (Sharp 1982 dalam Napitupulu 1993).
III METODE PENELITIAN 3.1 Waktu dan Tempat Penelitian ini dilaksanakan pada bulan September 2007 sampai dengan Mei 2008, bertempat di Laboratorium Energi dan Elektrifikasi Pertanian Leuwikopo, Departemen Teknik Pertanian, Fakultas Teknologi Pertanian, Institut Pertanian Bogor, Bogor. 3.2 Bahan dan Alat 3.2.1 Bahan Bahan yang digunakan dalam penelitian ini adalah 3000 kg jagung pipil varietas hibrida dengan kadar air rata-rata 18% b.k (range 18-20% b.k), setelah melalui proses pengeringan dengan laju tinggi (suhu udara pengering 50-60 oC) pada pengering ERK-hybrid. 3.2.2 Alat Peralatan yang digunakan meliputi bangunan In-Store Dryer (ISD), termokopel (CA), hybrid recorder HR-2500E, chino recorder, anemometer merek Kanomax Model A541, moisture tester, oven drying tipe SS-204D, termometer (bola basah dan bola kering), timbangan digital Tipe EK-1200 A, kain kasa, kawat kasa, botol kaca 140 ml, jangka sorong, mistar ukur, note book dan personal computer (PC) dengan software Visual Basic 6.0, Gambit 2.2.30 & Fluent 6.1.18 3.2.3 Bangunan In-Store Dryer Bangunan ISD yang diuji adalah suatu bangunan silo yang berbentuk silinder dengan ukuran tinggi 3.50 m dan diameter 2.50 m, memiliki kapasitas 7.5 ton jagung. Seluruh dinding ISD terbuat dari plat esser yang dilapisi galvanis dengan ketebalan 0.002 m, yang diperkuat oleh rangka dari pipa-pipa besi. Dinding terdiri dari dua lapisan, yaitu bagian luar dan dalam. Diantara kedua lapisan dinding tersebut diisi dengan busa glasswool sebagai insulator agar pemanasan oleh radiasi matahari tidak mempengaruhi kondisi dalam bangunan ini, sehingga dinding dalam kondisi adiabatis. Pada bagian atas bangunan ini terdapat lubang sebagai outlet udara dan juga untuk lubang loading bahan dengan
32
diameter bukaan 0.60 m. Gambar 6 menunjukkan skema bangunan ISD. Untuk lebih jelas posisi ISD di dalam sistem pengering terintegrasi, dapat dilihat pada gambar Pengering Efek Rumah Kaca (ERK)-Hybrid dan In-Store Dryer (ISD) terintegrasi yang disajikan pada Lampiran 1. 3 8
Keterangan :
1
1 2 3 4 5 6 7 8
2
5
Pipa input Pipa output Outlet udara ISD Kipas ISD Katup penutup Lantai pengering berlubang Saluran outlet biji-bijian Pintu kontrol
6 7
4
Gambar 6 Skema Bangunan ISD Bagian dalam bangunan ISD ini dilengkapi dengan 13 buah pipa penyalur udara dengan posisi yang diatur untuk meratakan distribusi aliran udara di dalam ISD. Pipa-pipa tersebut terbuat dari plat esser berpori (diameter pori 0.004 m) yang digalvanis dengan ketebalan 0.002 m, pipa-pipa ini menurut fungsinya terdiri dari dua jenis; yaitu pipa input dan pipa output. Pipa input berjumlah 9 pipa dengan diameter 0.15 m dan tingginya 0.22 m dari lantai ISD. Pipa input ini berhubungan langsung dengan lantai pengering, ujung pipa bagian bawah merupakan bukaan sebagai input udara dari ruang plenum. Pipa input dikondisikan setengah berpori dengan bagian yang berpori menghadap ke dinding ISD. Pengkondisian pipa input setengah berpori didasarkan profil aliran udara seperti yang diberikan secara skematis oleh Brooker et al. (1992) yang dapat dilihat pada Lampiran 2. Sementara pipa output berjumlah 4 pipa dengan diameter 0.20 m dan tingginya 2.0 m, pipa output ini ditempatkan dalam posisi tidak berhubungan langsung dengan lantai ISD, sehingga pipa ini dianggap sebagai pipa melayang dengan jarak 0.30 m diatas lantai ISD. Penempatan pipa input dan output udara di dalam ISD adalah berdasarkan pada profil aerasi udara, sehingga
33
udara dapat menyebar dengan merata. Hal ini didasarkan pada arah aliran udara pada pipa input dan output melalui tumpukan biji seperti yang dikemukakan oleh Brooker et al. (1992) yang secara skematis dapat dilihat pada Lampiran 3. Lokasi dari pipa-pipa penyalur udara pada koordinat bidang xz dapat dilihat pada Tabel 9, sementara penyajian gambar susunan pipa-pipa penyalur udara dalam ISD dapat dilihat pada Lampiran 4. Tabel 9 Lokasi pipa-pipa penyalur udara dalam ISD menurut fungsinya pada bidang xz. x(m) z (m) Diameter (m) 0.8 0.8 0.15 0 - 0.8 0.15 -0.8 0.8 0.15 0.8 - 0.8 0.15 0.8 0.8 0.15 - 0.8 - 0.8 0.15 -0.8 0 0.15 Jumlah
Fungsi Unit x(m) z (m) Diameter (m) Input 1 0.8 0 0.15 Input 1 0 0 0.15 Input 1 0.4 -0.4 0.20 Input 1 0.4 0.4 0.20 Input 1 -0.4 0.4 0.20 Input 1 -0.4 -0.4 0.20 Input 1 7
Fungsi Input Input Output Output Output Output -
Unit 1 1 1 1 1 1 6
Lantai ISD berbentuk plenum yang dilengkapi dengan lubang unloading bahan. Lantai terbuat dari plat esser berpori yang digalvanis, dengan ukuran sesuai dengan ukuran bangunan ISD yaitu diameter atas 2.50 m, sementara tinggi plenum ini 0.20 m. Pada bagian bawah plenum terdapat dua buah lubang berbentuk persegi berukuran 0.32 m x 0.20 m yang berfungsi untuk unloading bahan setelah selesai proses pengeringan dan penyimpanan. Pada bagian bawah ISD, di bawah lantai terdapat kipas axial sebagai penghembus udara lingkungan kedalam sistem ISD untuk proses pengeringan dengan laju rendah. Ukuran dan kapasitas ISD yang digunakan dalam penelitian ini adalah untuk skala penelitian, sementara ukuran dan kapasitas silo yang banyak digunakan untuk skala lapangan dan skala komersial sesuai dengan standar ASAE dapat lihat pada Lampiran 5.
34
3.3 Prosedur Penelitian Prosedur atau langkah kerja pada penelitian ini ditunjukkan pada diagram alir penelitian (Gambar 7). Mulai
Studi Literatur
Simulasi CFD : Suhu, RH dan Aliran Udara Pengering
Validasi Tanpa Beban
Tidak
Valid Ya Pemrograman & Simulasi Model Pengeringan Tumpukan
Pengujian dan Validasi Model dengan beban (Jagung)
Tidak
Valid Ya
Suhu, RH, Kecepatan udara, Kadar air, Mutu pengeringan dan penyimpanan, hasil evaluasi akhir untuk pengembangan ISD
Selesai
Gambar 7 Diagram Alir Penelitian
In-Store Dryer
35
3.3.1 Simulasi Sistem Thermal ISD Simulasi sistem thermal pada bagunan ISD dibuat dengan menggunakan software CFD yaitu Gambit 2.2.30 & Fluent 6.1.18, dengan langkah-langkah sebagai berikut: a) pembuatan geometri bangunan ISD menggunakan Gambit 2.2.30, meliputi: bangunan, lantai, pipa-pipa penyalur udara, inlet/kipas (kecepatan udara masuk) dan outlet. b) pembuatan mesh volume pada geometri dengan menggunakan Gambit 2.2.30. c) pendefinisian variabel operasi (meliputi; sifat termal bahan, menentukan kondisi batas yang sesuai pada sel yang merupakan batas domain berdasarkan bentuk saluran dalam ruangan, model simulasi berupa model energi, dan satuan yang digunakan pada geometri) yang dibuat dengan menggunakan Fluent 6.1.18. d) inisiasi dilakukan pada Fluent 6.1.18 .e) penentuan zona permukaan hasil simulasi yang akan ditampilkan pada Fluent 6.1.18. Adapun asumsi, kondisi awal dan kondisi batas dalam simulasi aliran udara dengan CFD dapat dilihat pada Lampiran 7. Simulasi CFD dilakukan dengan dua kondisi terhadap pipa saluran udara yang ditempatkan di dalam ISD, yaitu: 1. Simulasi 1: Bangunan ISD dikondisikan sesuai dengan bangunan ISD di lapangan dengan 9 buah pipa input setengah berpori, sisi dari pipa yang berpori diposisikan menghadap ke dinding. Sementara 4 buah pipa output seluruhnya berpori. 2. Simulasi 2: Bangunan ISD dengan pipa input dan output yang seluruhnya berpori, Simulasi ini merupakan modifikasi pada jenis pipa input. Pada Simulasi ini seluruh pipa input dikondisikan memiliki pori keseluruhan. 3.3.2 Simulasi Pengeringan Tumpukan (Lapis Tebal) Jagung Simulasi pengeringan tumpukan (lapis tebal) jagung yang dibuat dalam program komputer Visual Basic 6.0, untuk menyelesaikan persamaan-persamaan secara simultan, yaitu; Persamaan 23, Persamaan 27, Persamaan 37 dan Persamaan 45. Algoritma dan parameter simulasi dapat dilihat pada Lampiran 8 dan 9, sementara kode program Visual Basic dan Interface dari program simulasi tersebut disajikan pada Lampiran 10 dan 11. Simulasi dilakukan dengan dua kondisi yaitu:
36
1. Simulasi 1: untuk kondisi proses pengeringan pada musim hujan, dilakukan dengan ketebalan tumpukan 2.50 m, dengan kadar air awal jagung 18% b.k, suhu jagung 29.5 oC, suhu udara masuk 31 oC dengan RH 73% dan laju massa udara 12.7 kg/mnt-m2. 2. Simulasi 2: untuk kondisi proses pengeringan pada musim kemarau, dilakukan dengan ketinggian tumpukan 2.50 m, dengan kadar air awal jagung 18% b.k, suhu jagung 29.5 oC, suhu udara masuk 33 oC dengan RH 59.8% dan laju massa udara 12.7 kg/mnt-m2. 3.3.3 Percobaan Pengeringan dan Penyimpanan Jagung Percobaan pengeringan dan penyimpanan jagung di dalam ISD dilakukan pada kapasitas 1500 kg dari kapasitas total ISD sebesar 7500 kg. Pengeringan dilakukan dalam 2 kali percobaan untuk waktu yang berbeda yaitu : Percobaan 1 : dilakukan pada musim hujan, yaitu pada bulan Nopember 2007 dengan massa jangung 1500 kg dan ketebalan tumpukan jagung 0.5 m. Percobaan 2 : dilakukan pada awal musim kemarau, yaitu pada bulan Mei 2008 dengan massa jagung 1500 kg dan ketebalan tumpukan jagung 0.5 m. 3.3.4 Distribusi Udara Sebagai instrumen pengering dan penyimpanan yang terintegrasi dengan alat pengering ERK-Hybrid, pengeringan dan penyimpanan dalam ISD dilakukan setelah mengalami tahapan pengeringan dengan laju tinggi pada pengering ERKHybrid. Setelah kadar air pada kisaran 18-20% b.k, jagung dari ERK-Hybrid dipindahkan ke dalam bangunan ISD untuk proses penyimpanan dan pengeringan dengan laju pengeringan rendah sampai mencapai kadar air pada kisaran 15-16% b.k (13-14% b.b), dan mempertahankannya dengan mengontrol hembusan laju udara lingkungan ke dalam ruang ISD. Pengontrolan ini dilakukan dengan cara menghidupkan dan mematikan (ON/OFF) kipas berdasarkan perbandingan suhu dan RH udara bagian dalam bangunan ISD dengan suhu dan RH lingkungan. Pengontrolan tersebut dilakukan dengan seperangkat sensor suhu dan RH yang terhubung dengan komputer, kipas dihidupkan pada kisaran suhu 25-34 oC dengan
37
RH antara 50-90%. Secara rata-rata dalam kondisi cuaca yang cerah, kipas dihidupkan selama 10 jam yaitu antara pukul 7.00-17.00 WIB. Udara lingkungan sebagai media pengering dihembuskan ke dalam bangunan ISD dengan menggunakan kipas axial berukuran 15” sebagai pendistribusi udara, kipas ini digerakkan oleh motor dengan daya 2 hp (1.5 kW), perhitungan pressure drop dan tekanan statis kipas disajikan pada Lampiran 12. Ruangan ISD yang berkapasitas 7500 kg dibagi menjadi 5 tingkatan loading sesuai dengan kapasitas proses dari ERK-Hybrid yaitu 1500 kg setiap kali proses. Untuk mendapatkan distribusi udara yang merata, lantai ISD yang terbuat dari plat esser berpori juga dilengkapi dengan pipa-pipa penyalur udara. Susunan pipa-pipa penyalur udara dalam model ISD 3 dimensi dapat dilihat pada Gambar 8.
Gambar 8 Model ISD 3 Dimensi untuk simulasi CFD 3.3.5 Pengukuran Kecepatan Udara dan Suhu Kajian terhadap distribusi udara dalam ruangan ISD diamati melalui pengukuran parameter suhu, kecepatan udara, RH dan perubahan massa produk dan kadar air produk yang dikeringkan dengan selang waktu tertentu. Pengukuran parameter suhu dan kecepatan dilakukan pada posisi-posisi tertentu yang dianggap mewakili distribusi dalam ruangan ISD. Perubahan massa dan kadar air produk diukur dengan pengambilan sampel produk pada setiap lapisan dalam ruangan ISD dalam jangka waktu tertentu.
38
Pengukuran kecepatan aliran udara dilakukan dengan menggunakan anemometer. Pengukuran ini dilakukan sebanyak tiga kali ulangan pada titik-titik yang sama dengan pengukuran suhu. Titik pengukuran kecepatan udara sesuai titik-titik unit termokopel pengukuran suhu dapat dilihat pada Tabel 10. Pengukuran suhu dilakukan dengan menggunakan thermokopel CA yang dihubungkan dengan hybrid recorder, pencatatan dilakukan setiap 2 jam selama proses pengeringan berlangsung. Lokasi titik-titik pengukuran suhu dan kecepatan dilakukan pada bidang zx untuk ketinggian y dari 0 (lantai ISD/dasar bahan) untuk masing-masing tingkatan tumpukan, dan pada bidang zy untuk jarak dari sisi dinding ISD untuk masing-masing tingkatan tumpukan. Tebal setiap tumpukan untuk setiap kali loading adalah 0.50 m. Pembentukan posisi pengukuran dilakukan dengan membuat grid dari kawat halus, gambar titik-titik pengukuran sebanyak 20 titik dapat dilihat pada Lampiran 13, koordinat lokasi titik-titik pengukuran tersebut dapat dilihat pada Tabel 10. Tabel 10 Koordinat lokasi titik-titik pengukuran suhu, kecepatan udara dan RH dalam ISD dengan termokopel X(m) Y (m) 0.8 0.75 - 0.8 0.75 0.2 0.75 - 0.2 0.75 0.8 1.25 - 0.8 1.25 0.2 1.25 - 0.2 1.25 0.8 1.75 - 0.8 1.75 1.05 0.75 1.07 0.75 -1.05 1.75 -1.07 1.75 Jumlah
Z (m) 0.2 - 0.2 - 0.8 0.8 0.2 - 0.2 - 0.8 0.8 0.2 - 0.2 0 0 0 0
Unit 1 1 1 1 1 1 1 1 1 1 1 bk (RH1) 1 bb (RH1) 1 bk (RH2) 1 bb (RH2) 14
X(m) 0.2 - 0.2 0.8 - 0.8 0.2 - 0.2 0.8 - 0.8 0.2 - 0.2 0 0 -
Y (m) 1.75 1.75 2.25 2.25 2.25 2.25 2.75 2.75 2.75 2.75 3.25 3.25 -
Z (m) - 0.8 0.8 0.2 - 0.2 - 0.8 0.8 0.2 - 0.2 - 0.8 0.8 0 0 -
Unit 1 1 1 1 1 1 1 1 1 1 1 bk (RH3) 1 bb (RH3) 12
3.3.6 Pengukuran RH RH udara pengering diukur dengan menggunakan termokopel yang dibuat sebagai bola basah dan bola kering. Pengukuran dilakukan pada 4 titik yaitu: RH lingkungan 1 titik dan 3 titik di dalam bangunan ISD dengan jarak masing-masing kearah sumbu y dari bidang zx (lantai pengering). Pengukuran RH dilakukan setiap 2 jam selama proses pengeringan berlangsung. Koordinat titik pengukuran RH dalam ruang ISD dapat dilihat pada Tabel 10.
39
Untuk data pengukuran yang digunakan pada validasi terhadap hasil simulasi CFD, pengukuran RH menggunakan bola kering dan bola basah dilakukan pada 20 titik sesuai dengan titik pengukuran suhu yang hanya dilakukan satu kali pada kondisi steady. 3.3.7 Pengukuran Kadar Air Jagung Kadar air jagung diukur dengan menggunakan moisture tester dan oven drying. Metode pengukuran kadar air dengan oven drying adalah berdasarkan metode primer (Henderson & Perry 1976), persentase kadar air dihitung dengan Persamaan (5) dan (6). Pengambilan sampel ukur dilakukan pada 6 titik pada garis tengah diameter ISD (Lampiran 14) koordinat lokasi titik-titik pengambilan sampel pada bidang xz dengan ketinggian y dapat dilihat pada Tabel 11. Pengukuran dilakukan setiap 2 jam, selama masa pengeringan dan 1 kali setelah selesai penyimpanan. Tabel 11 Koordinat lokasi titik-titik pengambilan sampel untuk pengukuran kadar air jagung x(m) 0 0 0 0 0 0 Jumlah
y (m) 0.1 0.1 0.1 0.1 0.1 0.1
z (m) 0.9 0.55 0.2 -0.2 -0.55 - 0.9
Sampel 1 1 1 1 1 1 6
x(m) 0 0 0 0 0 0
y (m) 0.4 0.4 0.4 0.4 0.4 0.4
z (m) 0.9 0.55 0.2 -0.2 -0.55 - 0.9
Sampel 1 1 1 1 1 1 6
3.3.8 Validasi Model Simulasi Validasi dilakukan untuk membandingkan hasil pengukuran dan hasil simulasi yang didapat, pada titik-titik dan lokasi tertentu. Kriteria hasil validasi dianalisis dengan metode curve-fitting dan standar deviasi. Sementara besarnya error dalam validasi dihitung dengan persamaan mean absolute error (MAE) berikut ini: ∑
|
| .......................................................... (46)
Validasi data pada CFD meliputi suhu, aliran udara dan RH, sementara validasi data kadar air jagung dilakukan berdasarkan hasil simulasi model pengeringan tumpukan pada Visual Basic 6.0.
40
3.3.9 Uji Mutu Jagung Hasil Pengeringan dan Penyimpanan di dalam InStore Dryer Pegujian mutu jagung hasil pengeringan menurut SNI dilakukan di Laboratorium Balai Penelitian Tanaman Obat dan Aromatik (BALITRO) untuk pengujian proksimat meliputi: kadar protein (Metode SNI 2356-1991), kadar abu (Metode AOAC 942-05), kadar lemak (Metode SNI 01-2363-1991), kadar serat (AOAC 962.09 tahun 1982). Sementara pengujian kontaminasi aflatoxin (menggunakan metode ELISA) dilakukan di Laboratorium Balai Besar Penelitian Veteriner (BBALITVET). Sampel untuk pengujian diambil pada awal dan akhir proses pengeringan dan setelah penyimpanan selama 30 hari, dengan cara sampling pada layer yang mewakili. Disamping itu juga dilakukan uji viabilitas (Metode SNI 01-6944-2003) dengan cara mengkecambahkan biji.
IV HASIL DAN PEMBAHASAN 4.1 Simulasi Aliran Udara pada ISD 4.1.1 Pembentukan Grid Perhitungan Untuk memecahkan persoalan aliran fluida salah satu yang penting diperhatikan dalam simulasi CFD adalah pembentukan grid. Pada permasalahan aliran fluida dalam penelitian ini, ukuran grid yang digunakan adalah 0.001 m untuk seluruh bentuk bangunan geometri ISD. Volume grid sebesar 72108 volume dan jumlah titik (node) sebanyak 17747. Grid dalam batas volume ISD dibuat menggunakan element yang tidak terstruktur dengan model tet/hybrid dan tipe Tgrid, selanjutnya grid tersebut dilengkapi dengan kondisi batas. Hasil pembentukan grid untuk model ISD yang disimulasikan disajikan pada Gambar 9.
Gambar 9 Pembentukan grid pada domain perhitungan
42
4.1.2 Distribusi Suhu dan Kecepatan Udara Pengering Simulasi 1 Distribusi suhu dan kecepatan udara pengering dalam ISD pada Simulasi 1 ditunjukkan pada Gambar 10, 11,12 dan 13. Nilai Hasil Simulasi 1 selengkapnya dapat dilihat pada Lampiran 15. Udara lingkungan yang dihembuskan oleh kipas pendistribusi terdorong ke depan sepanjang lubang inlet dan melewati lantai berpori kemudian disebarkan ke seluruh ruangan ISD. Suhu udara yang paling besar berada di depan inlet dan juga pada bagian bawah ISD yaitu pada kisaran 33.6-34 oC. Setelah menyebar ke seluruh ruangan ISD suhu mulai berkurang, namun sebaran suhu ruangan di bagian atas dan tengah cenderung seragam yang ditunjukkan oleh warna orange dengan nilai kisarannya 33.3-33.6 oC dan nilai rata-rata 33.5 oC. Sebaran suhu dalam ISD ini ternyata juga dipengaruhi oleh efek pori sebagian pada pipa input. Pipa input dengan sebagian berpori ini ternyata menghalangi sebaran udara yang membawa suhu masuk, sehingga sebaran suhu di sekitar pipa input menjadi sangat bervariasi yaitu antara 27-33.7 oC. variasi ini ditunjukkan oleh warna biru, hijau dan kuning pada Gambar 10.
Bagian tengah dan atas ISD, kisaran suhu 33.3-33.6oC Suhu di sekitar pipa input antara 27-33.7 oC
Suhu di depan Inlet antara 33.6-34 oC
Bagian bawah ISD, kisaran suhu 33.6-34oC
Gambar 10 Distribusi suhu udara di dalam ISD Simulasi 1
43
Pengaruh pipa yang berpori sebagian terhadap variasi sebaran suhu di sekitarnya juga terjadi pada pipa-pipa di lokasi lainnya, sehingga secara akumulasi pengaruhnya terhadap ketidakseragaman sebaran suhu di sekitarnya juga semakin besar. Kondisi ini dapat dilihat pada Gambar 11.
Suhu disekitar pipa inlet antara 27-33.7oC
Gambar 11 Pengaruh pipa input berpori sebagian terhadap sebaran suhu ISD pada Simulasi 1. Kecepatan aliran udara tertinggi terdapat di depan kipas yang masuk melalui inlet, yaitu pada kisaran 7.79-8.19 m/dtk yang ditunjukkan oleh warna merah. Aliran udara ini melewati lantai pengering yang berpori dan sebagian lainnya juga masuk melalui pipa-pipa input dan pipa-pipa output menuju outlet ISD. Saat mencapai di tengah ruangan kecepatan udara mulai menurun berkisar pada 0.410.81 m/dtk ditunjukkan oleh warna biru, pada gambar juga terlihat adanya kecepatan aliran dengan kisaran 0-0.41 m/dtk yang disebabkan oleh halangan pori-pori lantai ISD. Secara visual sebaran kecepatan aliran udara dapat dilihat pada Gambar 12. Pipa dengan berpori setengah juga mempengaruhi vektor aliran udara. Posisi pori yang tertutup dan menghadap ke bagian dalam ruangan menjadi penghalang bagi udara untuk menembus ke arah dinding ataupun sebaliknya, sehingga membuat udara harus memutar melalui pori di depannya dan menuju ke arah belakang pipa. Hal ini membuat kecepatan udara menjadi sangat berkurang
44
dan bahkan sampai kisaran 0-0.4 m/dtk di lokasi-lokasi antara dinding dan bagian pipa yang tidak berpori, kondisi ini dapat dilihat pada Gambar 13.
0-0.41 m/dtk 0.41-0.81 m/dtk Kecepatan aliran udara di depan inlet antara 7.79- 8.19 m/dtk
0.81-1.23 m/dtk
Gambar 12 Distribusi kecepatan udara didalam ISD pada Simulasi 1
Lokasi antara pipa setengah berpori dan dinding, kecepatan aliran 0-0.41 m/dtk. Aliran udara tidak bisa menembus langsung ke arah dinding
Aliran udara memutar di depan bidang pipa setengah berpori untuk berbalik menuju arah dinding
Gambar 13 Pengaruh pipa input berpori sebagian terhadap vektor aliran udara dalam ISD pada Simulasi 1
45
4.1.3 Distribusi Suhu dan Kecepatan Udara Pengering Simulasi 2 Pada Simulasi 2, keseluruhan pipa input dibuat berpori, untuk melihat kemungkinan pengurangan hambatan dari sebaran suhu dan kecepatan aliran udara yang terdapat pada ISD di lapangan seperti yang telah terlihat pada Simulasi 1. Distribusi suhu dan kecepatan udara pengering dalam ISD pada Simulasi 2 ditunjukkan pada Gambar 14, 15, 16 dan 17. Nilai hasil Simulasi 2 selengkapnya dapat dilihat pada Lampiran 16. Kipas pendistribusi mendorong udara lingkungan sepanjang lubang inlet dan dihembuskan ke dalam ruang ISD. Kecepatan aliran udara yang membawa suhu udara lingkungan melewati lantai berpori untuk kemudian disebarkan ke seluruh ruangan ISD. Suhu udara yang paling besar berada di depan inlet, bagian bawah, dan sebagian sisi ISD di depan kipas (velocity inlet) yaitu pada kisaran suhu 33.7-34 oC. Setelah menyebar ke seluruh ruangan ISD suhu mulai berkurang, namun sebaran suhu ruangan dibagian atas dan tengah cenderung seragam antara range 33.4-33.7 oC yang ditunjukkan oleh warna orange, dengan nilai rata-rata 33.5 oC. Sebaran suhu pada Simulasi 2 dapat dilihat pada Gambar 14.
kisaran suhu 33.4-33.7oC
kisaran suhu di depan inlet 33.7 -34oC
kisaran suhu 33.7-34oC
Gambar 14 Distribusi suhu udara di dalam ISD Simulasi 2
46
Sebaran suhu dalam ISD pada Simulasi 2 ternyata tidak dipengaruhi oleh efek pipa input. Hal ini dikarenakan pipa input yang berpori keseluruhan mampu memperlancar sebaran udara yang membawa suhu udara secara lebih merata. Sebaran suhu di sekitar pipa input terlihat masih sama dengan sebaran suhu secara keseluruhan yaitu berkisar 33.4-34 oC, kondisi tersebut dapat dilihat pada Gambar 15.
Pipa inlet yang berpori seluruhnya tidak mempengaruhi sebaran suhu
Gambar 15 Pengaruh pipa input berpori sebagian terhadap sebaran suhu ISD pada Simulasi 2 Kecepatan aliran udara tertinggi pada Simulasi 2 terdapat di depan kipas yang masuk melalui inlet, yaitu pada kisaran 7.79-8.19 m/dtk yang ditunjukkan oleh warna merah. Aliran udara ini melewati lantai pengering yang berpori dan sebagian lainnya juga masuk melalui pipa-pipa input dan pipa-pipa output menuju outlet ISD. Pada bagian bawah ruang ISD setelah melalui lantai, kecepatan udara berkisar antara 0.81-1.23 m/dtk, saat mencapai ditengah ruangan kecepatan udara mulai menurun berkisar pada 0.41-0.81 m/dtk ditunjukkan oleh warna biru, pada Simulasi 2 ini juga menunjukkan bahwa masih ada lokasi-lokasi yang mempunyai kecepatan pada kisaran 0-0.41 m/dtk akibat pengaruh tahanan lantai ISD yang berpori. Secara visual, profil sebaran kecepatan aliran udara dapat dilihat pada Gambar 16.
47
0-0.41 m/dtk 0.41-0.81 m/dtk Kecepatan aliran udara di depan inlet antara 7.79-8.19 m/dtk
0.81-1.23 m/dtk
Gambar 16 Distribusi kecepatan udara di dalam ISD pada Simulasi 2 Pipa input yang seluruhnya berpori pada Simulasi 2 ini ternyata mampu mengurangi halangan pergerakan aliran udara di dalam ISD, dibandingkan ketika menggunakan pipa input dengan setengah berpori yang mempengaruhi sebaran aliran udara pada Simulasi 1. Posisi pori di keseluruhan permukaan pipa mampu membantu aliran udara menembus ke segala arah. Khusus untuk bagian yang menghadap ke dalam ruangan ISD yang tadinya terhalangi oleh bidang tidak berpori sehingga menghalangi udara untuk menembus ke arah dalam ruangan ataupun ke arah dinding, pada Simulasi 2 tidak terjadi lagi. Karena udara bisa menembus langsung ke arah dalam ruangan maupun ke arah dinding tanpa harus memutar melalui bidang yang berpori. Hal ini ditunjukkan dengan ruang yang memiliki kecepatan aliran udara antara 0-0.41 m/dtk antara pipa input dan dinding ISD menjadi berkurang, walaupun masih ada namun hal ini bukanlah halangan yang berasal dari pipa saja, tetapi merupakan akumulasi oleh halangan lantai dan pipa yang menyebabkan berkurangnya tekanan aliran udara. Namun demikian, pemakaian pipa yang berpori seluruhnya ternyata dapat membantu pergerakan aliran udara ke segala arah. Pengaruh pipa input terhadap vektor aliran udara pada Simulasi 2 dapat dilihat pada Gambar 17.
48
Lokasi antara pipa inlet dan dinding, berkurangnya ruang dengan kecepatan aliran 0-0.41 m/dtk
Aliran udara dapat menembus langsung kearah dinding melalui pori-pori pipa input
Gambar 17 Pengaruh pipa input berpori sebagian terhadap vektor aliran udara dalam ISD pada Simulasi 2 4.1.4 Distribusi RH Udara Hasil Simulasi RH di dalam bangunan ISD ditentukan berdasarkan perhitungan dengan menggunakan Persamaan (3). Pada kondisi suhu lingkungan 34 oC dan suhu bola basah 30 oC, diperoleh RH lingkungan sebesar 75%. Secara keseluruhan RH hasil simulasi untuk bangunan ISD di lapangan (Simulasi 1) rata-rata RH berkisar antara 67.1% sampai dengan 68%, data selengkapnya disajikan pada Lampiran 15. Sedangkan sebaran RH rata-rata hasil Simulasi 2 berkisar antara 60.3% sampai dengan 60.5%, yang disajikan pada Lampiran 16. 4.1.5 Keragaman Kecepatan Udara, Suhu dan RH Tingkat keragaman kecepatan aliran udara, suhu dan RH pada kedua simulasi didapatkan dari nilai rata-rata hasil simulasi. Nilai keragaman dan masing-masing standar deviasi tersebut digunakan untuk mengevaluasi pengaruh pemakaian jenis pipa input terhadap keseragaman kecepatan aliran udara, suhu dan RH pada kondisi operasi yang sama.
49
Keragaman sebaran suhu ditunjukkan oleh profil suhu pada 5 ketinggian dalam ruang ISD hasil simulasi. Untuk Simulasi 1 dan 2 dapat dilihat pada Gambar 18 dan 19.
Rata-rata 31.8oC Rata-rata 31.9oC Rata-rata 31.9 oC Rata-rata 32oC Rata-rata 31.8oC
Gambar 18 Profil suhu pada 5 ketinggian Simulasi 1
Rata-rata 33.6oC Rata-rata 33.6oC Rata-rata 33.7oC Rata-rata 33.7oC Rata-rata 33.7oC
Gambar 19 Profil suhu pada 5 ketinggian Simulasi 2
50
Pada Simulasi 1 (ISD lapangan), rata-rata suhu terendah terdapat pada ketinggian 0.75 dan 2.75 m yaitu 31.8 oC sementara suhu tertinggi terdapat pada ketinggian 1.25 m sebesar 32 oC, rata-rata suhu untuk kelima ketinggian sebesar 32 oC. Pada Simulasi 1 didapatkan standar deviasi untuk sebaran suhu sebesar 2.6 o
C. Pada Simulasi 2 rata-rata suhu terendah terdapat pada ketinggian 2.25 dan
2.75 m sebesar 33.6 oC, sedangkan suhu tertinggi terdapat pada ketinggian 0.75, 1.25, dan 1.75 m sebesar 33.7 oC. Nilai standar deviasi sebaran suhu pada Simulasi 2 sebesar 0.06 oC. Perbandingan nilai keragaman suhu pada kedua simulasi secara grafis dapat dilihat pada Gambar 20, sementara datanya disajikan pada Lampiran 17.
34.0 33.5
Suhu (oC)
33.0 32.5 Simulasi-1
32.0
Simulasi-2 31.5 31.0 30.5 0.75
1.25
1.75
2.25
2.75
Ketinggian (m)
Gambar 20 Keragaman suhu pada kedua simulasi Profil kecepatan aliran udara pada ke-5 ketinggian hasil simulasi untuk Simulasi 1 dan 2 dapat dilihat pada Gambar 21 dan 22. Kecepatan aliran udara pada kedua simulasi juga mengalami keragaman dan perbedaan. Pada Simulasi 1 kecepatan aliran udara terkecil terdapat pada ketinggian 2.75 m dengan nilai 0.23 m/dtk, sementara yang kecepatan terbesar terdapat pada ketinggian 1.25 m sebesar 0.40 m/dtk, sedangkan rata-rata untuk seluruh ketinggian sebesar 0.32 m/dtk. Standar deviasi untuk sebaran kecepatan aliran udara pada Simulasi 1 sebesar 0.22 m/dtk.
51
Rata-rata 0.23 m/dtk Rata-rata 0.27 m/dtk Rata-rata 0.32 m/dtk Rata-rata 0.40 m/dtk Rata-rata 0.39 m/dtk
Gambar 21 Profil kecepatan aliran udara pada 5 ketinggian Simulasi 1
Rata-rata 0.37 m/dtk Rata-rata 0.36 m/dtk Rata-rata 0.44 m/dtk Rata-rata 0.53 m/dtk Rata-rata 0.59 m/dtk
Gambar 22 Profil kecepatan aliran udara pada 5 ketinggian Simulasi 2
52
Untuk Simulasi 2, kecepatan aliran udara terbesar berada pada ketinggian 0.75 m sebesar 0.59 m/dtk, sedangkan kecepatan terkecil berada pada ketinggian 2.25 m sebesar 0.36 m/dtk. Kecepatan aliran udara rata-rata pada Simulasi 2 adalah sebesar 0.46 m/dtk dengan standar deviasi sebesar 0.20 m/dtk. Perbandingan keragaman kecepatan aliran udara pada kedua simulasi disajikan secara grafis pada Gambar 23, sementara penyajian datanya dapat dilihat pada Lampiran 17.
Kecepatan aliran udara (m/dtk)
0.70 0.60 0.50 0.40 Simulasi-1 0.30
Simulasi-2
0.20 0.10 0.00 0.75
1.25
1.75
2.25
2.75
Ketinggian (m)
Gambar 23 Keragaman kecepatan aliran udara pada kedua simulasi Keragaman RH pada kedua simulasi juga menunjukkan perbedaan yang signifikan. Data perbandingan keragaman RH udara pada kedua simulasi disajikan pada Lampiran 17 dan secara grafis dapat dilihat pada Gambar 24. Pada Simulasi 1, didapatkan RH terendah sebesar 67.1% yang berada pada ketinggian 1.25 m, sementara nilai RH tertinggi terdapat pada ketinggian 2.75 m sebesar 68%. Nilai RH rata-rata untuk Simulasi 1 adalah 67.6% dengan standar deviasi bernilai 10.7%. Untuk Simulasi 2, RH terendah berada pada ketinggian 0.75 m sebesar 60.3% sementara RH tertinggi sebesar 60.5% yang terdapat pada ketinggian 1.75, 2.25 dan 2.75 m. Rata-rata RH pada Simulasi 2 adalah 60.4% dengan standar deviasi sebesar 0.2%.
53
70.0 68.0
RH udara (%)
66.0 64.0 Simulasi-1 62.0
Simulasi-2
60.0 58.0 56.0 0.75
1.25
1.75
2.25
2.75
Ketinggian (m)
Gambar 24 Keragaman RH udara pada kedua simulasi Berdasarkan perbandingan nilai keragaman dan standar deviasi untuk parameter sebaran suhu, kecepatan aliran udara dan RH pada kedua simulasi, terlihat bahwa Simulasi 2 memiliki keragaman dengan nilai standar deviasi yang lebih kecil. Hal tersebut menunjukkan bahwa Simulasi 2 memiliki tingkat keseragaman suhu, kecepatan aliran udara dan RH yang lebih baik dibandingkan Simulasi 1. 4.2 Validasi Suhu, Kecepatan Aliran Udara dan RH Hasil simulasi model aliran udara menggunakan CFD berupa kontur distribusi suhu dan kecepatan aliran udara, menunjukkan besar dan arah aliran udara sebagai media pengering yang digunakan dalam ISD. Validasi model dilakukan dengan membandingkan data ukur kecepatan aliran udara dan suhu pada 20 titik hasil pengukuran dengan hasil simulasi. Nilai hasil simulasi dan hasil pengukuran suhu pada bidang xz dengan ketinggian y dapat dilihat pada Lampiran 18. Gambar 25 menunjukkan penyajian grafis hasil validasi suhu dengan membandingkan suhu hasil pengukuran dan hasil simulasi. Perbedaan antara data suhu hasil pengukuran dan suhu hasil simulasi CFD diberikan dalam nilai standar deviasi sebesar 0.45 oC. Penyimpangan dari validasi ini dinyatakan dalam total error sebesar 12.64 oC, dengan rata-rata error 0.63 oC pada range 0.22-1.19 oC.
54
Hubungan antara suhu hasil simulasi CFD dengan suhu hasil pengukuran memiliki nilai korelasi sebesar 0.66. 40 SD = 0.45 oC Total error = 12.64 oC Rata-rata error= 0.63 oC R = 0.66
39 38
Suhu (C)
37 36 35 34 33 32 31 30 0
1
2
3
4
5
y=0.75 m
6
7
8
y=1.25 m
9 10 11 12 13 14 15 16 17 18 19 20 21 y=1.75 m
y=2.25 m
Point pengukuran pada bidang xz
y=2.75 m T-Ukur T-CFD
Gambar 25 Validasi suhu udara hasil simulasi terhadap suhu pengukuran Validasi kecepatan aliran udara antara hasil pengukuran dan hasil simulasi disajikan pada Gambar 26, sementara nilai hasil simulasi dan hasil pengukurannya dapat dilihat pada Lampiran 18.
Kecepatan aliran udara (m/dtk)
4.0
SD = 0.20 m/dtk Total error = 5.71 m/dtk Rata-rata error = 0.29 m/dtk R= 0.73
3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 0
1
2
3
y=0.75 m
4
5
6
7
y=1.25 m
8
9 10 11 12 13 14 15 16 17 18 19 20 21
y=1.75 m
y=2.25 m
Titik pengukuran pada bidang xz
y=2.75 m V Ukur V Cfd
Gambar 26 Validasi kecepatan aliran udara hasil simulasi terhadap data pengukuran.
55
Dari grafik validasi terlihat bahwa kecepatan aliran udara hasil simulasi telah mengikuti data pengukuran, walaupun pada bebarapa titik terlihat jelas perbedaannya. Perbedaan ini disebabkan oleh penentuan jarak grid yang tidak persis sama antara pengukuran dan simulasi. Pada ketinggian 0.75 m sampai 1.25 m terlihat perbedaan yang besar antara nilai ukur dan simulasi, namun pada ketinggian 1.75 m sampai 2.75 m perbedaan berkurang dan nilai simulasi mendekati nilai pengukuran. Perbedaan nilai pengukuran dan hasil simulasi dinyatakan dalam nilai standar deviasi sebesar 0.20 m/dtk. Sementara penyimpangannya dinyatakan dalam total error sebesar 5.71 m/dtk dengan ratarata error 0.29 m/dtk pada range 0.04-0.97 m/dtk. Hubungan antara kecepatan aliran udara hasil simulasi dengan hasil pengukuran memiliki nilai korelasi sebesar 0.73. 85 SD = 0.91% Total error = 25.85% Rata-rata error= 1.29% R=0.66
80 75
RH (%)
70 65 60 55 50 45 0
1
2
3
y=0.75 m
4
5
6
7
y=1.25 m
8
9 10 11 12 13 14 15 16 17 18 19 20 21 y=1.75 m
y=2.25 m
Titik pengukuran pada bidang xz
y=2.75 m RH ukur RH-Hitung
Gambar 27 Validasi RH hasil perhitungan terhadap RH hasil pengukuran Validasi RH dilakukan dengan membandingkan RH hasil pengukuran dengan RH perhitungan berdasarkan suhu hasil simulasi CFD Fluent 6.1. Hasil validasi ini dapat dilihat pada Gambar 27, nilai hasil perhitungan dan hasil pengukuran RH disajikan pada Lampiran 19. Secara umum terlihat bahwa adanya kesamaan antara RH hasil pengukuran dan RH hasil perhitungan dari suhu hasil simulasi CFD. Perbedaan dalam validasi ini dinyatakan dengan standar deviasi sebesar 0.91%, sementara penyimpangannya dinyatakan dalam total error sebesar
56
25.85% dengan rata-rata error 1.29% pada range 0.03-2.66%. Hubungan antara RH hasil perhitungan dengan RH hasil pengukuran pada ISD lapangan memiliki nilai korelasi sebesar 0.66. 4.3 Perubahan Kadar Air Jagung pada Simulasi 4.3.1 Kadar Air Jagung Simulasi 1 Simulasi 1 merupakan simulasi proses pengeringan pada kondisi musim hujan dengan ketebalan tumpukan 2.50 m, dengan kadar air awal jagung 18% b.k, suhu jagung 29.5 oC, suhu udara masuk 31 oC dengan RH 73% dan laju massa udara 12.7 kg/mnt-m2. Kondisi tersebut merupakan kondisi lapangan hasil pengukuran. Perubahan kadar air setelah simulasi pengeringan selama 150 jam selengkapnya dapat dilihat pada Lampiran 20 dan penyajian grafisnya dapat dilihat pada Gambar 28.
20.0
Suhu udara = 31oC, RH=73 %, Laju massa udara=12.7 kg/mnt-m2, Me=16.6% b.k
19.5
Kadar air (% b.k)
19.0
Layer 1
18.5
Layer 10 Layer 40
18.0
Layer 50
17.5
Layer 100 17.0
Layer 150
16.5
Layer 200 Layer 250
16.0 0
10
20
30
40
50
60
70
80
90 100 110 120 130 140 150
Waktu (jam ke-)
Gambar 28 Perubahan kadar air pada Simulasi 1 Secara keseluruhan pada awal-awal pengeringan berlangsung, pada seluruh layer terjadi peningkatan kadar air dalam waktu yang bervariasi. Pada Layer 1 peningkatan terjadi hanya pada 1 jam pertama dari kadar air sebesar 0.8%, peningkatan ini dikarenakan terjadinya kondensasi uap air dari lapisan sebelumnya. Ketebalan setiap layer sebesar 0.01 m memungkinkan adanya dua
57
lapisan biji yang bertumpuk dan ukuran biji yang bervariasi menyebabkan asumsi satu layer hanya ada satu lapisan biji menyebabkan hasil tidak persis tepat. Sehingga pada Layer 1 pun masih mungkin terjadi kenaikan kadar air akibat uap air yang dibawa udara dari lapisan sebelumnya. Namun hal tersebut tidak berlangsung lama, karena setelah 1 jam pertama kadar air pada Layer 1 terus menurun menuju kadar air keseimbangannya. Hal yang sama juga terjadi pada Layer 10 (meningkat sebesar 1.41% b.k), Layer 40 (meningkat sebesar 1.48% b.k), Layer 50 dan 100 (masing-masing meningkat sebesar 1.50% b.k). Peningkatan kadar air pada ke-4 layer tersebut juga terjadi pada satu jam pertama pengeringan berlangsung, namun setelah satu jam kadar air pada ke-4 layer tersebut mengalami penurunan. Pada Layer 150, peningkatan kadar air masih terjadi sampai jam ke-3 proses pengeringan berlangsung sebesar 1.52% b.k. Setelah jam ke-3 kadar air pada layer ini menurun. Pada Layer 200 terjadi peningkatan sampai jam ke-4 sebesar 1.52% b.k dan pada Layer 250 terjadi peningkatan kadar air sampai jam ke-5 dengan nilai sebesar 1.59% b.k. Secara umum terlihat bahwa trend penurunan kadar air baru terjadi dengan baik setelah melewati 10 jam proses pengeringan, penurunan ini juga terjadi bervariasi sesuai dengan lapisannya masing-masing. Karena pada prinsipnya pengeringan tumpukan tebal ini adalah proses pengeringan lapisan tipis yang dikondisikan bertumpuk, maka uap air dari lapisan sebelumnya masih mungkin terkondensasi pada lapisan setelahnya tergantung pada suhu, RH dan kecepatan udara yang membawa uap air tersebut. Pada Layer 1, 10, 40 dan 50 terlihat penurunan tajam terjadi setelah jam ke-3
sampai jam ke-40, selanjutnya
penurunan kadar air berlangsung secara perlahan-lahan. Sementara pada layer setelahnya yaitu Layer 100, 150, 200 dan 250 laju penurunan kadar air lebih lambat dari keempat layer sebelumnya. Hal ini dikarenakan tumpukan biji yang semakin tebal, sehingga udara dan RH yang melewati layer-layer tersebut telah terlebih dahulu terjenuhkan oleh uap air yang dibawa dari layer-layer sebelumnya. Namun demikian, penurunan kadar air pada layer-layer tersebut terus berlangsung menuju kadar air keseimbangan.
58
Pada Simulasi 1 dengan suhu 31 oC dan RH 73% didapatkan Me sebesar 16.6% b.k. Dengan kondisi ini, maka dibutuhkan waktu selama 150 jam pengeringan untuk mengurangi kadar air dari 18% b.k menuju kadar air keseimbangan. Pada akhir Simulasi 1 (jam ke-150) didapat kadar air rata-rata sebesar 16.6% b.k, nilai tersebut telah mencapai kadar air keseimbangannya. 4.3.2 Kadar Air Jagung Simulasi 2 Simulasi 2 merupakan simulasi dengan kondisi proses pengeringan pada musim kemarau. Kondisi yang diberikan adalah ; ketebalan tumpukan jagung 2.50 m dengan kadar air awal 18% b.k, laju massa udara 12.7 kg/mnt-m2, suhu udara pengering yang digunakan adalah 33 oC dengan RH 59.8%. Hasil simulasi Simulasi 2 selengkapnya dapat dilihat pada Lampiran 21, sementara penyajian grafisnya pada Gambar 29. 21.5
Suhu udara = 33 oC, RH=59.8 %, Laju massa udara=12.7 kg/mnt-m2, Me=13.5% b.k
Kadar air (% b.k)
20.5 19.5
Layer 1
18.5
Layer 10 Layer 40
17.5
Layer 50
16.5
Layer 100
15.5
Layer 150 Layer 200
14.5
Layer 250
13.5 12.5 0
10
20
30
40
50
60
70
80
90
100
110
120
Waktu (jam ke-)
Gambar 29 Perubahan kadar air pada Simulasi 2 Pada Simulasi 2 terlihat bahwa pada semua lapisan tetap mengalami peningkatan kadar air akibat kondensasi uap air dari lapisan sebelumnya. Pada Layer 1 terjadi peningkatan kadar air sebesar 0.11% b.k pada jam ke-1 proses pengeringan. Setelah jam ke-2 kadar air pada Layer 1 terus menurun secara drastis hingga jam ke-30, selanjutnya penurunan berlangsung lambat. Pada Layer 10 terjadi peningkatan sebesar 1.12% b.k pada jam ke-1, sementara pada Layer 40
59
dan 50 juga terjadi peningkatan kadar air pada jam ke-1 (masing-masing sebesar 1.40% b.k dan 0.49% b.k). Sementara pada Layer 100 peningkatan kadar air terjadi hingga jam ke-2 sebesar 0.48% b.k. Memasuki jam ke-3 hingga jam ke-30 penurunan kadar air pada Layer 10, 40 dan 50 berlangsung secara drastis hingga melawati jam ke-30 penurunan kadar air mulai berlangsung secara perlahan. Pada Layer 100 laju penurunan yang tinggi terjadi antara jam ke-3 hingga jam ke-50 dan selanjutnya berlangsung lambat. Pada Layer 150 dan 200 kadar air meningkat hingga jam ke-3 (sebesar masing-masing 1.49% b.k dan 1.5% b.k) dan selanjutnya penurunan kadar air berlangsung dengan laju yang tinggi setelah jam ke-3 hingga jam ke-60. Untuk Layer 250 peningkatan kadar air hingga jam ke-6 pengeringan berlangsung dengan nilai 1.54% b.k. Selanjutnya perubahan kadar air berlangsung cepat pada jam ke-20 hingga jam ke-80 untuk Layer 200 dan jam ke30 hingga jam ke-90 untuk Layer 250. Selanjutnya penurunan terus berlangsung secara perlahan-lahan menuju kadar air keseimbangan. Perubahan RH dari 73.2% pada Simulasi 1 menjadi 59.8% pada Simulasi 2 ternyata memberikan pengaruh yang cukup besar pada proses pengeringan. Terjadinya peningkatan kadar air yang lebih kecil pada layer-layer yang ada merupakan pengaruh RH udara yang kering, karena semakin kering RH udara, maka uap air
yang diserap oleh udara akan semakin banyak. Hal tersebut
menyebabkan kemungkinan kondensasi pada layer-layer bagian atas menjadi lebih kecil karena kapasitas tampung uap air dalam udara yang kering menjadi lebih besar. Sehingga kalaupun terjadi kondensasi maka pengaruhnya sangat kecil bila dibandingkan kenaikan kadar air yang terjadi pada Simulasi 1. Kondensasi uap air dari RH kering pada Simulasi 2 juga mengakibatkan perlambatan laju penurunan kadar air pada awal-awal proses pengeringan, namun demikian penurunan kadar air terus berlangsung menuju kadar air keseimbangan. Pada Simulasi 2 dengan suhu 33 oC dan RH 59.8% didapatkan kadar air keseimbangan sebesar 13.5% b.k. Dengan kondisi ini, dibutuhkan waktu selama 120 jam pengeringan untuk mengurangi kadar air dari 18% b.k menuju kadar air keseimbangan tersebut. Pada akhir proses Simulasi 2 didapatkan kadar air ratarata sebesar 13.5% b.k, nilai tersebut memperlihatkan bahwa proses pengeringan pada simulasi ini telah mencapai kadar air keseimbangannya. Bila dilihat dari segi
60
waktu pengeringan dan perubahan kadar air bijian, maka udara lingkungan dengan RH 59.8% memiliki potensi yang sangat baik untuk proses pengeringan. 4.4 Perubahan Kadar Air Jagung Percobaan 4.4.1 Kadar Air Jagung Percobaan 1 Data pengukuran perubahan kadar air bijian dalam ISD selama 50 jam pengeringan menggunakan udara lingkungan disajikan pada Lampiran 22. secara grafis dapat dilihat pada Gambar 30. Hasil pengukuran data perubahan kadar air pada ISD di lapangan menunjukkan terjadinya fluktuasi kadar air pada biji-bijian. Fluktuasi ini disebabkan oleh suhu dan RH udara lingkungan yang dihembuskan ke dalam ruang ISD sangat bervariasi dan fluktuatif. Sehingga saat udara lembab masuk ke dalam ruangan ISD maka biji-bijian akan menyerap kandungan air yang berada di udara, dan juga sebaliknya ketika udara kering masuk, maka udara dengan RH yang kering tersebut akan kembali menampung dan membawa uap air dari biji-bijian sehingga terjadi penurunan kadar air pada bijian. Selama 50 jam pengeringan dengan RH udara rata-rata sebesar 73% dan suhu rata-rata 31 oC, diperoleh kadar air keseimbangan 16.6% b.k.
19.0
Kadar air (%.b.k)
18.5 18.0 17.5 17.0 16.5 16.0 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50
Waktu (jam ke-)
Layer 10-Ukur Layer 40-Ukur
Gambar 30 Kadar air hasil pengukuran selama 50 jam pengeringan pada Percobaan 1
61
Fluktuasi kadar air jagung yang terjadi di dalam ISD selama 50 jam pengeringan memperlihatkan trend menurun. Fluktuasi ini disebabkan oleh adanya kondensasi uap air yang dibawa udara dari layer-layer sebelumnya. Hal ini dapat dilihat dari penggambaran grafis antara jam ke-14 sampai jam ke-50, pada waktu-waktu tersebut ketika kadar air pada Layer 10 menurun maka pada Layer 40 kadar air nya meningkat. Di sini terlihat adanya perpindahan kandungan air yang dibawa udara ke Layer 40 dari layer-layer sebelumnya. Namun demikian fluktuasi tersebut tetap mengikuti trend penurunan menuju kadar air keseimbangannya. Pengaruh fluktuasi RH terhadap perubahan kadar air dapat dilihat pada
19.0
120
18.5
110 100
18.0 90 17.5 80 17.0 70 16.5
RH udara (%)
Kadar air (%.b.k)
Gambar 31, penyajian data nya disajikan pada Lampiran 23.
60
16.0
50 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50
Waktu (jam ke-)
RH udara Layer 10-Ukur Layer 40-Ukur
Gambar 31 Pengaruh fluktuasi RH terhadap perubahan kadar air bijian pada Percobaan 1 Secara umum ketika RH udara yang masuk ke ruangan ISD tinggi maka kadar air bijian juga menjadi bertambah, dan sebaliknya ketika RH udara yang masuk rendah maka kadar air juga akan berkurang. Kemampuan udara menyerap uap air dari bijian dan memindahkannya sangat penting dalam proses pengeringan, dan hal tersebut tentunya harus didukung oleh tingkat kekeringan pada RH udara.
62
4.4.1 Kadar Air Jagung Percobaan 2 Data pengukuran perubahan kadar air bijian dalam ISD selama 40 jam pengeringan menggunakan udara lingkungan pada Percobaan 2 disajikan pada Lampiran 24, secara grafis dapat dilihat pada Gambar 32. Hasil pengukuran kadar air pada Percobaan 2 memperlihatkan terjadinya fluktuasi kadar air pada bijibijian. Selama 40 jam pengeringan RH udara yang masuk ke ISD rata-rata sebesar 59.8% dan suhu rata-rata 33 oC diperoleh kadar air keseimbangan 13.5% b.k.
20.5 19.5
Kadar air (%.b.k)
18.5 17.5 16.5 15.5 14.5 13.5 12.5 0
2
4
6
8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40
Waktu (jam ke-)
Layer 10-Ukur Layer 40-Ukur
Gambar 32 Kadar air hasil pengukuran selama 40 jam pengeringan pada Percobaan 2 Fluktuasi kadar air jagung yang terjadi di dalam ISD selama 40 jam pengeringan memperlihatkan trend menurun. Fluktuasi tersebut terjadi karena kondensasi uap air yang dibawa udara dari layer-layer sebelumnya. Dari penggambaran grafis terlihat bahwa ketika kadar air pada Layer 10 menurun maka pada Layer 40 kadar air nya meningkat. Kondisi ini memperlihat adanya perpindahan kandungan air yang dibawa udara ke Layer 40 dari layer-layer sebelumnya. Secara umum fluktuasi kadar air pada Percobaan 2 mengikuti trend penurunan menuju kadar air keseimbangannya. Pengaruh fluktuasi RH terhadap perubahan kadar air dapat dilihat pada Gambar 33, penyajian data nya disajikan pada Lampiran 25. Dari penyajian grafis terlihat bahwa terjadi kecenderungan fluktuasi kadar air terhadap fluktuasi RH.
63
Pada Percobaan 2 dengan RH udara yang relatif lebih rendah (antara 55.5-72.7%) menyebabkan fluktuasi kadar air cenderung menurun. 20
85
19
80 75
17 16
70
15 65
14 13
60
RH udara (%)
Kadar air (%.b.k)
18
12 55
11 10
50 0
2
4
6
8
10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40
Waktu (jam ke-)
RH udara Layer 10-Ukur Layer 40-Ukur
Gambar 33 Pengaruh fluktuasi RH terhadap perubahan kadar air bijian pada Percobaan 2 4.5 Validasi Perubahan Kadar Air Jagung Pendugaan penurunan kadar air bijian memakai model pengeringan tumpukan tebal (bed) dimana tumpukan dibagi menjadi 50 layer dengan tinggi layer 0.5 m dan ketebalan setiap layer sebesar 0.01 m. Validasi ini dilakukan untuk data pada Layer 10 dan 40 yang diukur di lapangan terhadap data hasil simulasi. 4.5.1 Validasi Kadar Air Jagung Percobaan 1 Data hasil validasi antara pengukuran dan simulasi pada Percobaan 1 dapat dilihat pada Lampiran 26, sementara penyajian grafisnya dapat dilihat pada Gambar 34 dan 35. Hasil validasi antara kadar air pengukuran dan kadar air simulasi pada Layer 10 didapatkan standar deviasi sebesar 0.31% b.k, dengan total error sebesar 11.57% b.k dan rata-rata error sebesar 0.44% b.k pada range 0.07-0.92% b.k. Hasil validasi ini menunjukkan bahwa antara data hasil pengukuran dan data simulasi tidak jauh berbeda, walaupun pada beberapa waktu perbedaan tersebut terlihat sangat jelas. Namun dengan nilai error rata-rata sebesar 0.44% b.k menunjukkan bahwa antara kedua data tersebut memiliki
64
penyimpangan yang relatif kecil. Penurunan kadar air pada Layer 10 memiliki korelasi sebesar 0.90.
20.0
Suhu udara = 31oC, RH=73 %, Me=16.6 b.k, Laju massa udara=12.7 kg/dtk-m2
19.5
SD = 0.31% b.k Total error = 11.57% b.k Rata-rata error= 0.44% b.k R = 0.90
Kadar air (%.b.k)
19.0 18.5 18.0 17.5 17.0 16.5 16.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50
Waktu (jam ke-)
L10-simulasi L10-ukur
Gambar 34 Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 10 pada Percobaan 1
20.0
Suhu udara = 31oC, RH=73 %, Me=16.6 % b.k, Laju massa udara=12.7 kg/dtk-m2
19.5
SD = 0.31% b.k Total Error = 11.23% b.k Rata-rata error= 0.43% b.k R = 0.35
Kadar air (%.b.k)
19.0 18.5 18.0 17.5 17.0 16.5 16.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50
L40-simulasi
Waktu (jam ke-)
L40-ukur
Gambar 35 Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 40 pada Percobaan 1
65
Validasi kadar air pada Layer 40 menunjukkan penyimpangan yang tidak jauh berbeda dari Layer 10. Besarnya penyimpangan yang terjadi antara data hasil pengukuran dan data hasil simulasi diberikan dalam nilai total error sebesar 11.23% b.k dengan rata-rata error 0.43% b.k pada range 0.01-1.89% b.k dan standar deviasi sebesar 0.31% b.k. Pada Layer 40 penurunan kadar air mempunyai korelasi sebesar 0.35. Nilai korelasi yang kecil pada penurunan kadar air di Layer 40 ini disebabkan oleh adanya variasi kadar air awal yang masuk ke dalam ISD. Variasi tersebut menyebabkan hasil pengukuran pada percobaan tidak sepenuhnya bisa mengikuti trend simulasi. Penyimpangan tersebut dapat dilihat pada jam ke-6 sampai jam ke-12, antara rentang waktu tersebut ketika trend pada simulasi menurun hasil pengukuran justeru menunjukkan peningkatan kadar air. Disamping variasi kadar air awal, peningkatan kadar air tersebut juga dapat disebabkan oleh kondensasi uap air yang di bawa udara dari layer-layer sebelumnya. 4.5.2 Validasi Kadar Air Jagung Percobaan 2 Data hasil validasi antara pengukuran dan simulasi pada Percobaan 2 dapat dilihat pada Lampiran 27, sementara penyajian grafisnya dapat dilihat pada Gambar 36 dan 37. Hasil validasi antara kadar air pengukuran dan kadar air simulasi Layer 10 pada Percobaan 2 didapatkan standar deviasi sebesar 0.32% b.k, dengan total error sebesar 9.49% b.k dan rata-rata error sebesar 0.45% b.k pada range 0.06-3.12% b.k. Hasil validasi ini menunjukkan bahwa antara data hasil pengukuran telah mengikuti trend hasil simulasi dengan baik, walaupun pada beberapa waktu terjadi perbedaan yang sangat jelas. Perbedaan jelas terlihat antara jam ke-8 sampai jam ke-16, pada rentang waktu tersebut data pengukuran mengalami peningkatan, sementara pada simulasi telah memperlihatkan trend penurunan. Namun dengan nilai error rata-rata sebesar 0.45% b.k menunjukkan penyimpangan yang relatif kecil antara data pengukuran dan simulasi. Penurunan kadar air pada Layer 10 memiliki nilai korelasi sebesar 0.88.
66
Suhu udara = 33oC, RH=59.8 %, Me=13.5 b.k, Laju massa udara=12.7 kg/dtk-m2 20.0 SD = 0.32% b.k Total error = 9.49% b.k Rata-rata error= 0.45% b.k R = 0.88
19.0
Kadar air (%.b.k)
18.0 17.0 16.0 15.0 14.0 13.0 12.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50
Waktu (jam ke-)
L10-simulasi L10-ukur
Gambar 36 Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 10 pada Percobaan 2
Suhu udara = 34oC, RH=59.8 %, Me=13.5 % b.k, Laju massa udara=12.7 kg/dtk-m2 20.0
SD = 0.30% b.k Total Error = 9.06% b.k Rata-rata error= 0.43% b.k R = 0.84
19.0
Kadar air (%.b.k)
18.0 17.0 16.0 15.0 14.0 13.0 12.0
0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50
Waktu (jam ke-)
L40-simulasi L40-ukur
Gambar 37 Validasi kadar air bijian hasil simulasi terhadap hasil pengukuran Layer 40 pada Percobaan 2 Validasi perubahan kadar air hasil pengukuran dan simulasi Layer 40 Percobaan 2 menunjukkan penyimpangan yang tidak jauh berbeda dari Layer 10. Pada Layer 40 didapatkan total error sebesar 9.06% b.k dengan rata-rata error
67
0.43%b.k pada range 0.07-2.81% b.k dan standar deviasi sebesar 0.30% b.k. Penurunan kadar air pada Layer 40 memiliki nilai korelasi sebesar 0.84. Perbedaan antara kadar air hasil pengukuran dan hasil simulasi pada Percobaan 1 dan 2 yang telah divalidasi diduga karena tidak tepatnya penerapan kelembaban udara pengering pada tumpukan. Seperti asumsi yang digunakan bahwa pengeringan tumpukan atau lapisan tebal ini merupakan kumpulan dari banyak lapisan-lapisan tipis. Namun bila pada pengeringan lapisan tipis pengeringan terjadi pada kelembaban yang relatif sama dengan udara pengering, asumsi ini akan membuat penyimpangan pada kasus lapisan tebal. Hal ini dikarenakan pada lapisan tebal uap air yang keluar dari bijian tidak semuanya langsung terbawa oleh udara pengering, sehingga lapis batas bijian (boundary layer) lebih tebal. Kondisi ini akan menambah kesalahan (error) yang cukup berarti apa lagi ketika kadar air bijian masih tinggi. Namun ketika kadar air telah berkurang, maka penerapan nilai RH pada simulasi lebih mendekati kondisi yang sebenarnya. Hal ini secara umum terlihat setelah 12 jam waktu pengeringan, pada kedua percobaan model simulasi telah dapat mengikuti data pengukuran. Menurut Brooker et al. (1974), perbedaan antara hasil simulasi dan percobaan dipengaruhi oleh kurang akuratnya model pengeringan lapisan tipis, kurang tepatnya persamaan kadar air isotermis biji pada RH tinggi dan tidak tepatnya nilai parameter input model. Asumsi bahwa tidak terjadinya penyusutan volume (shrinkage) juga turut menyumbangkan kesalahan dalam perhitungan. Perbedaan lainnya yang cukup berarti adalah bahwa penyusunan model lapisan tipis adalah berdasarkan pada kondisi suhu dan RH tetap, sedangkan kenyataannya pada percobaan sangat berfluktuasi. 4.6 Analisis Mutu Jagung Hasil Pengeringan dan Penyimpanan 4.6.1 Mutu Perdagangan Proses pengeringan selama 40 jam dalam In-Store Dryer (ISD) menggunakan udara lingkungan, telah mampu mengurangi kadar air jagung dari rata-rata 16.40% b.b (range 16.02-16.77% b.b) pada awal proses pengeringan menjadi 13.40% b.b (range 11.88-13.14% b.b) pada akhir proses pengeringan. Selama masa penyimpanan kadar air pada jagung terus turun menuju kadar air keseimbangannya sesuai dengan suhu dan RH lingkungan, setelah 30 hari
68
penyimpanan kadar air rata-rata jagung dalam ISD menjadi 12.33% b.b (range 11.0-14.0% b.b). Perbandingan parameter mutu hasil percobaan dengan persyaratan mutu SNI 01-03920-1995 (standar mutu jagung untuk perdagangan) disajikan pada Gambar 38 sementara data selengkapnya dapat dilihat pada Lampiran 28.
Persyaratan mutu (% max)
18 16 14
Syarat Mutu I
12
Syarat Mutu II
10
Syarat Mutu III
8
Syarat Mutu IV
6
Sebelum Pengeringan
4
Setelah Pengeringan
2
Setelah Penyimpanan
0
Parameter mutu
Gambar 38 Perbandingan parameter mutu hasil percobaan dengan SNI jagung untuk perdagangan Perbandingan antara parameter mutu hasil percobaan dengan persyaratan mutu SNI 01-03920-1995 menunjukkan bahwa mutu jagung selama proses pengeringan dan penyimpanan berada pada tingkat mutu IV. Perubahan yang signifikan terjadi pada kadar air, sementara parameter yang sangat mempengaruhi tingkat mutu pada jagung hasil percobaan adalah butir rusak sebesar 7.03% sehingga mutu jagung hasil percobaan berada pada tingkat mutu IV. Kerusakan ini terjadi pada saat pemindahan jagung oleh Bucket Elevator setelah proses pengeringan tahap awal pada ERK-Hybrid, disamping itu kerusakan juga dapat terjadi akibat benturan-benturan saat loading dan unloading serta proses perputaran silinder ketika berlangsungnya pengeringan pada ERK-Hybrid. Sementara selama proses pengeringan dan penyimpanan pada ISD butir rusak mengalami peningkatan sebesar 0.01%, sementara parameter lain yang meningkat adalah butir warna lain sebesar 0.08%, dan kotoran sebesar 0.04%. Peningkatan persentase pada parameter-parameter tersebut terjadi karena adanya fluktuasi suhu
69
dan RH dalam ISD, perbedaan suhu dan kelembaban yang terjadi terus menerus dalam waktu lama dapat menyebabkan keretakan pada bijian dan juga perubahan warna bijian. Bijian yang retak ataupun pecah menjadi bagian-bagian yang halus selama penyimpanan akan menjadi kotoran dari bijian itu sendiri, namun secara umum peningkatan nilai parameter-parameter tersebut tidak terlalu signifikan, sehingga tidak memberikan pengaruh besar pada tingkat mutu dari hasil proses pengeringan dan penyimpanan di dalam ISD. 4.6.2 Mutu Benih Persyaratan SNI untuk benih jagung hibrida di dasarkan pada SNI 01-69442003. Hasil pengujian tingkat kemurnian benih dan viabilitas sebelum dan setelah proses dalam ISD dan perbandingan dengan SNI disajikan pada Gambar 39 sementara penyajian data selengkapnya dapat dilihat pada Lampiran 28.
120
Nilai (%)
100 80
Persyaratan SNI Sebelum proses Setelah proses
60 40 20 0 Kadar air (% maks)
Benih murni Daya Kotoran Benih (% min) berkecambah (% maks) (% min) Parameter persyaratan
Gambar
39 Hasil pengujian tingkat kemurnian benih dan viabilitas sebelum dan setelah proses dalam ISD dan perbandingan dengan SNI
Kadar air rata-rata yang terukur setelah proses pengeringan dan penyimpanan dalam ISD sebesar 12.33% b.b masih melebihi persyaratan SNI sebesar maksimum 12% b.b. Namun kalau melihat dari range kadar air yang terukur yaitu antara 11 sampai dengan 14% b.b, maka sebenarnya kadar air setelah proses dalam ISD tersebut telah masuk ke dalam keriteria SNI, hanya saja variasi kadar
70
air yang terjadi menyebabkan nilai kadar air tersebut belum dapat memenuhi persyaratan SNI. Benih murni yang didapatkan dari hasil percobaan sebesar 85.8% belum memenuhi standar yang disyaratkan SNI sebesar 98%. Hal tersebut dikarenakan jagung yang diamati tercampur oleh kotoran benih sebesar 14.5%, berupa : butir pecah, butir rusak dan butir warna lain seperti yang telah paparkan pada persyaratan mutu menurut SNI 01-03920-1995. Sementara untuk daya berkecambah sebelum proses pengeringan dan penyimpanan dalam ISD didapatkan nilai sebesar 85%, nilai tersebut telah memenuhi standar SNI minimum,
namun
setelah
proses
pengeringan
dan
penyimpanan
daya
berkecambah menjadi 84.7% atau turun sebesar 0.3%. penurunan tersebut memang tidak signifikan namun menyebabkan batas minimum dari SNI menjadi tidak tercapai. Prosentase daya berkecambah yang hanya bisa mencapai batas minimum dari yang disyaratkan SNI disebabkan oleh suhu pengeringan pada ERK-Hybrid yang berkisar antara 50 sampai dengan 60 oC. Suhu pengeringan tersebut telah melampaui suhu udara maksimum untuk pengeringan benih jagung sebesar 43.3 oC (Hall 1970). Tidak terpenuhinya persyaratan mutu dari hasil percobaan bukanlah efek langsung dari proses pengeringan dan penyimpanan selama di dalam ISD, namun hal tersebut terjadi akibat pengeringan tahap awal dengan suhu tinggi pada ERK-Hybrid dan juga proses pemindahan bijian dari ERK-Hybrid ke ISD yang mengakibatkan terjadinya kerusakan pada jagung. Sebagai unit penyimpan, ISD telah berfungsi dengan baik yang ditunjukkan oleh tidak terjadinya peningkatan kerusakan bijian selama pengeringan dan penyimpanan, disamping itu pengurangan kadar air pada pengeringan dengan suhu rendah juga telah menunjukkan kinerja yang baik dengan penurunan kadar air yang signifikan. 4.6.3 Mutu Pakan Ternak Untuk dijadikan pakan ternak, maka mutu jagung harus memenuhi persyaratan SNI 01-4483-1998. Perbandingan kadar proksimat jagung antara hasil percobaan dan persyaratan mutu SNI dapat dilihat pada Gambar 40 sementara data selengkapnya disajikan pada Lampiran 28.
71
18
Nilai (%)
16 14
Persyaratan SNI
12
Sebelum Proses
10
Setelah Proses
8 6 4 2 0 Kadar air (% maks)
Kadar protein Kadar serat kasar kasar (% min) (% maks)
Kadar abu (% maks)
Kadar lemak (% min)
Komposisi Nutrisi (kimia) jagung
Gambar 40 Perbandingan kandungan nutrisi (kimia) jagung antara hasil percobaan dan persyaratan SNI Komposisi nutrisi (kimia) jagung sebelum dan setelah proses pengeringan dan penyimpanan di dalam ISD mengalami perubahan pada masing-masing parameter. Kadar air mengalami penurunan dari rata-rata 16.40% b.b sebelum pengeringan dan penyimpanan menjadi rata-rata 12.33% b.b setelah proses pengeringan dan penyimpanan, penurunan tersebut melebihi standar yang disyaratkan yaitu sebesar maksimum 14% b.b. Sementara untuk nilai protein, sebelum pengeringan dan penyimpanan didapat rata-rata sebesar 7.76% nilai ini lebih baik dari standar SNI sebesar minimum 7.5%, namun setelah proses pengeringan dan penyimpanan nilai protein jagung berkurang menjadi rata-rata 6.62% sehingga tidak memenuhi lagi syarat minimum SNI. Untuk kadar serat didapatkan rata-rata 2.56% sebelum proses dan rata-rata 2.53% setelah proses berlangsung, nilai tersebut lebih baik dari standar SNI yang mensyaratkan kadar serat maksimum 3%. Kadar abu juga menunjukkan bahwa jagung hasil percobaan memiliki kadar abu yang lebih baik dari standar yang ditetapkan SNI, yaitu rata-rata sebesar 1.24% sebelum proses dan 1.23% setelah proses berlangsung, sementara SNI menetapkan kadar abu maksimum sebesar 2%. Untuk kadar lemak, nilai pada jagung hasil percobaan juga masih lebih baik dari standar SNI yaitu rata-rata sebesar 3.80% dan 3.24% untuk sebelum dan setelah proses, sementara SNI menetapkan nilai kadar lemak minimum sebesar 3%.
72
Kontaminasi mikotoksin jenis aflatoxin B1 pada jagung percobaan sebelum pengeringan dan penyimpanan didapatkan nilai rata-rata sebesar 18.48 ppb, sementara setelah proses pengeringan dan penyimpanan menjadi rata-rata 21.10 ppb. Hal tersebut menunjukkan bahwa selama penyimpanan 30 hari dalam ISD terjadi peningkatan kontaminasi aflatoxin sebesar 2.63 ppb. Namun nilai kontaminasi aflatoxin selama proses dalam ISD tersebut masih dibawah ambang batas toleransi untuk konsumsi ternak yang ditetapkan SNI sebesar maksimum 50 ppb. Untuk komponen butir pecah, butir warna lain, benda asing dan kepadatan, masih dalam batas persyaratan SNI baik sebelum maupun setelah proses pengeringan dan penyimpanan dalam ISD. Sehingga berdasarkan kandungan gizi pada analisa proksimat dan kontaminasi aflatoxin serta keseluruhan parameter dalam SNI 01-4483-1998, maka jagung hasil pengeringan dan penyimpanan di dalam ISD masih memenuhi syarat yang baik untuk dijadikan pakan ternak.
V SIMPULAN DAN SARAN 5.1 Simpulan 1. Distribusi suhu dan aliran udara di dalam ISD dapat disimulasikan dengan menggunakan teknik CFD. Hasil validasi kecepatan aliran udara dan suhu antara simulasi terhadap pengukuran didapatkan koefisien korelasi 0.73 dan 0.66. Sebaran suhu dan kecepatan aliran udara lebih seragam dengan jenis pipa input yang berpori seluruhnya dibandingkan jenis pipa input setengah berpori. 2. Validasi sebaran RH antara hasil perhitungan dan pengukuran didapatkan korelasi 0.66. Pipa input yang berpori seluruhnya memberikan sebaran RH yang lebih seragam dibandingkan jenis pipa input setengah berpori. 3. Dari kapasitas total ISD sebesar 7500 kg, dilakukan validasi penurunan kadar air pada kapasitas 1500 kg. Hasil validasi antara simulasi terhadap pengukuran pada layer 10 dan 40, didapatkan nilai korelasi 0.90 dan 0.35 pada musim hujan, 0.88 dan 0.84 pada musim kemarau. 4. Dari analisis mutu jagung menunjukkan bahwa proses pengeringan dan penyimpanan di dalam ISD mampu mempertahankan mutu jagung dengan baik sebagai bahan pakan ternak. 5.2 Saran 1. Untuk pengembangan ISD sebaiknya digunakan pipa input dengan pori seluruhnya, mengingat sebaran aliran udara, suhu dan RH hasil simulasi CFD untuk jenis pipa tersebut lebih seragam. 2. Untuk kinerja pengeringan yang lebih baik, perlu pengontrolan terhadap input aliran udara ke dalam ISD dengan menyalakan kipas pada range RH yang lebih rendah (55-70%), dan juga perlu dilakukan pengkondisian suhu udara sebelum masuk ke dalam ISD dengan cara pemanasan udara. 3. Perlu dilakukan validasi pada kapasitas total ISD serta analisis ekonomi untuk pengembangan ISD.
DAFTAR PUSTAKA Abdullah K. 1993. Optimization of Solar Drying System. Proc. of the 5th International Energy Conference. Seoul, October 18-22, 1993. Abdullah K. 1995. Optimasi Dalam Perencanaan Alat Pengering Hasil Pertanian Dengan Energi Surya. Laporan Akhir Penelitian Hibah Bersaing. Direktorat Pembinaan Penelitian dan Pengabdian Pada Masyarakat. Direktorat Jendral Pendidikan Tinggi. Kontrak No. 039/P4M/DPPM/PHB/95. Abdullah K. 1998. Greenhouse Effect Solar Dyer for Coffe and Cocoa beans. Final Report. University Research for Graduate Education. Contract No. 032/HTPP-II/URGE/1996. Directorate General of Higher Education, Indonesia. [Anonim]. 1994. ASAE Standard. USA. [Anonim]. 1998. SNI 01-4483-1998 Jagung Bahan Baku Pakan. Jakarta: BSN. [Anonim]. 2003. SNI 01-6944-2003 Benih Jagung Hibrida. Jakarta: BSN. Bala BK. 1997. Drying and Storage of Cereal Grains. New Delhi : Oxford & IBH Publishing Co. Pvt. Ltd. Brooker DB, Bakker-Arkema FW, Hall CW. 1974. Drying Cereal Grain. Connecticut : The AVI Publishing Company Inc. Wesport. Brooker DB, Bakker-Arkema FW, Hall CW. 1992. Drying and Storage of Grains and Oilseeds. New York: An Avi Book, Van Nostrand Reinhold. Cengel YA. 2003. Heat Transfer. New York : Mc.Graw-Hill,Inc. Chikubu S. 1974. Characteristic of Japanese Rice and Storage Principle of Brown Rice. National Food Research Institut, Ministry of Agriculture and Forestry. Fiscal. Devahastin S. 2000. Panduan Praktis Mujumdar untuk Pengeringan Industrial. Tambunan AH, Wulandani D, Hartulistiyoso E, Nelwan LO, Penerjemah. Bogor : IPB Press. Terjemahan dari: Handbook of Industrial Drying 2nd Edition. Dharmaputra O, Retnowati S, Sunjaya, Ambarwati S. 1993. Populasi A.flavus dan kandungan aflatoksin pada jagung ditingkat petani dan pedagang di Propinsi Lampung. Yogyakarta: Makalah Kongres Nasional XII dan Seminar Ilmiah Perhimpunan Fitopatologi Indonesia. Francis BJ, Wood JF. 1982. Changes in The Nutritive Value of Feed Concentrate during Storage; Hand Books of Nutritive Value of Processed Food Vol. II Di dalam: Milolov Rechyl Jr, Editor. CRC Series in Nutrition and Food. Florida: CRC Press Inc.
75
Garraway MO, Evans QC. 1984. Fungal Nutritions and Physiology. New York: John Wiley and Sons Inc. Hall C.W. 1970. Handling and Storage Food Grain in Tropical and Subtropical Areas. Roma: FAO. Hall CW, Davis DC. 1979. Processing Equipment for Agricultural Product. Connecticut: The AVI Publishing Company,Inc. Westport. Hall CW. 1980. Drying and Storage of Agricultural Crops. Connecticut: The AVI Publishing Company,Inc. Westport. Harlos LS, Jeon YW, Bockhop CW. 1983. Design and performance of multipurpose dryer using non-conventional energy sources. Manila: Agricultural Engineering Departement. IRRI. Hartono. 1978. Teknologi Lepas Panen Jagung dan Palawija Serupa Lainnya. Jakarta: Lokakarya jagung, Bulog. Henderson SM, Perry RL. 1976. Agricultural Process Engineering. 3rd Edition. Conecticut: The AVI Publishing Company, Inc. Wesport. Holman JP. 1995. Perpindahan Kalor. Jakarta : Erlangga. Jeon WY, Harlos LS, Bockhop CW.1983. Evaluation of a multipurpose dryer using non-conventional energy sources. Manila: Agricultural Engineering Departement. IRRI. Kim KS, Shin MG, Kim BC, Thim JH, Cheigh HS, Muhlbauer W, Kwon TW. 1989. An Ambient-air In-Storage Paddy Drying System for Korean Farm. Agricultural Machanization in Asia, Africa and Latin America., 20:2. Komar N. 1988. Mempelajari Sistem Lumbung Pengering Gabah Bahan Bakar Sekam [tesis]. Bogor: Program Pascasarjana, Institut Pertanian Bogor. Koto HA. 1984. Simulasi Penyimpanan Gabah dalam Silo Besi Kedap Udara [tesis]. Bogor: Program Pascasarjana, Institut Pertanian Bogor. Kristanto A. 2007. Teknologi Pascapanen untuk Peningkatan Mutu Jagung. http://www.tanindo.co.id/abdi11/hal0901.hmt. [30 Januari 2007]. Manalu LP. 1999. Pengering Energi Surya dengan Pengaduk Mekanis untuk Pengeringan Kakao[tesis]. Bogor: Program Pascasarjana, Institut Pertanian Bogor. Nagler MJ, Buangsuwon D, Jewer K, Siriacha P. 1986. Production and quality controll of low aflatoxin maize in rainy seasons. Bangkok: Proceding of Departement of Agriculture Research Conference.
76
Napitupulu VM. 1993. Rancangan dan Uji Kinerja Kipas Untuk Pengeringan [tesis]. Bogor: Program Pascasarjana, Institut Pertanian Bogor. Nelwan LO. 2005. Study on Solar Assisted Dryer with Rotating Rack for Cocoa Beans [disertasi]. Bogor: Sekolah Pascasarjana, Institut Pertanian Bogor. Nelwan LO. 1997. Pengering Kakao dengan Energi Surya Menggunakan Rak Pengering dengan Kolektor Tipe Efek Rumah Kaca [tesis]. Bogor: Program Pascasarjana, Institut Pertanian Bogor. Nugroho E. 1986. Simulasi Pengeringan Gabah [skripsi]. Bogor: Fateta IPB. Soemangat, Syarief AM, Subekti D, Purwadaria HK. 1987. Studies on the implementation of the pit dryer at the village level in Yogyakarta, Indonesia. Bangkok: Asean Seminar on Grain Post Harvest Technology. Soemartono. 1968. Teknik Pengolahan Padi. Jakarta: Departemen Pertanian RI. Soesarsono W. 1977. Teknik Pengolahan dan Penyimpanan Hasil Panen. Bogor: Departemen Teknologi Hasil Petanian Fatemeta IPB. Suseno H. 1974. Fisiologi Tumbuhan. Bogor: Departemen Botani Faperta IPB. Thahir R, Sudaryono, Soemardi, Soeharmadi. 1988. Teknologi Pascapanen Jagung. Bogor: Pusat Penelitian dan Pengembangan Tanaman Pangan. Versteeg HK, Malalasekera W. 1995. An introduction to computational fluid dynamic. The finite volume methode. Malaysia: Longman Sc & Technical. Widodo TW, Harjono, Tokumoto O, Matsumoto. 1994. Hasil Analisis Teknis dan Ekonomis Proses Pengeringan Padi dengan Menggunakan “DS System”. Enjiniring Pertanian 1:1-12. Wijandi S. 1988. Teknologi Penyimpanan Komoditas Pertanian. Bogor: Fatemeta Institut Pertanian Bogor. Wilcke WF, Morey RV, Hansen DJ. 1993. Reducing Energy Use for Ambientair Corn Drying. Applied Engineering in Agriculture 9:2. William PC. 1991. Storage of Grain and Seeds. London: CRC Press Inc. Wulandani D. 2005. Kajian Distribusi Suhu, RH dan Aliran Udara Pengering untuk Optimasi Disain Pengering Efek Rumah Kaca [disertasi]. Bogor: Sekolah Pascasarjana, Institut Pertanian Bogor.
LAMPIRAN
78
Sistem Pengerinng ERK-Hybridd dan In-Store Dryer Terintegrassi Lampiran 1 S
ANGAN KETERA In-Store Dryer : 1 2 3 4 5 6 7 8
aluran inlet Sa Sa aluran outlet Ou utlet udara ISD Kip pas ISD Ka atup penutup Lantai pengering berlu ubang Sa aluran outlet biji-bijian Pin ntu kontrol
ng ERK-Hybrid : Pengerin 1 2 3 4 5 6 7 8 9
Silinder pengering Kip pas Pe enukar panas Kip pas ekshaus Mo otor penggerak silind der Tangki air Ho opper bahan bakar Mo otor dan screw feede er bahan bakar Tungku pembakaran
79
Lampiran 2 Profil alirann udara padaa pipa setenggah berpori (Brooker et al. a 1992)
80
Lampiran 3 Arah alirann udara pada pipa input ddan output (B L Brooker et all. 1992)
81
Lampiran 4 Susunan pippa-pipa penyyalur udara dalam L d ISD
82
Lampiran 5 Standar ASAE untuk ukuran dan kapasitas Silo
83
Lampiran 5 (lanjutan)
84
Lampiran
6 Algoritma numerik volume hingga dengan metode SIMPLE (Versteeg & Malalasekera 1995)
85
Lampiran
7 Asumsi, kondisi awal dan kondisi batas yang digunakan pada simulasi CFD
Asumsi 1. Udara tidak termanpatkan (incompressible), ρ konstan. 2. Panas jenis, konduktivitas dan viskositas udara konstan (bilangan Prandtl udara konstan). 3. Udara bergerak dalam kondisi steady 4. Udara lingkungan dianggap konstan selama simulasi, yaitu pada 34oC. 5. Kecepatan udara pada kipas dianggap konstan. 6. Aliran udara dalam ISD dengan bukaan yang sangat lebar dianggap sebagai aliran udara di atas bidang datar (sepanjang bidang datar) dan laminer apabila nilai bilangan Reynolds kurang dari 500 000 (Cengel, 2003). Perhitungan udara laminer: Sifat thermofisik udara yang digunakan dalam model simulasi CFD adalah: Suhu inlet : 34oC Suhu outlet : 34oC Suhu fluida operasi : (34+34)/2 = 34oC --- > = 307.15 K Diameter spesifik = 2.5 m Sifat Massa jenis Panas jenis Konduktivitas panas Visikositas dinamik Bilangan Prandtl
Simbol ρ Cp k μ Pr
Nilai 1.225 1.0057 0.02624 1.7894 x 10-5 0.71
Satuan kg/m3 kJ/kgoC W/m.K kg/m.s -
Bilangan Reynold (Re) : 1.225 0.2 2.5 1.7894 x10 Sehingga aliran dianggap laminer.
34229.35
86
Lampiran 7 (lanjutan) Kondisi Awal : Untuk kedua simulasi, menggunakan kondisi awal sebagai berikut: 1. Kecepatan aliran udara awal pada semua koordinat (x,y dan z) = 0 m/detik. 2. Suhu dinding = suhu lingkungan 3. Tekanan udara 1 atm = 101.325 kPa Kondisi Batas : 1. Inlet sekaligus kipas dianggap sebagai velocity inlet dengan kecepatan udara masuk 8.2 m/detik (Hasil pengukuran). 2. Outlet dianggap sebagai outflow dengan ratio bukaan 1. 3. Dinding bersifat adiabatis (tidak ada pertukaran panas) 4. Ketebalan dinding 0.002 m 5. Dinding dan seluruh pipa penyalur udara terbuat dari alumanium (plat esser) dengan propertis (Holman 1995): Density (ρ) : 2719 kg/m3 Specific Heat (Cp) :871 j/kg-k Thermal Conductivity (k) : 202.4 w/m-k 6. Porositas lantai dan pipa-pipa penyalur udara adalah 40%, dalam simulasi ditetapkan sebagai porous jump, dengan parameter input: a. Permeabilitas permukaan (α) dihitung dengan persamaan (FLUENT ver.6.1): 150 1
0.0045 150
0.432 1 0.432
3.374 10
b. Koefisien porous jump (C2) dihitung dengan persamaan (FLUENT ver.6.1): 3.5 1
3.5 1 0.432 0.004 0.432
Keterangan : Dp : diameter lubang pada plat = 0.004 m ε
: persen lubang pada plat = 0.4
6 164.60 1/
87
Lampiran 8 Algoritma simulasi pengeringan tumpukan tebal Mulai
Masukkan data kadar air awal, suhu jagung, suhu udara, kelembaban mutlak udara, laju massa udara, tebal tumpukan, total waktu simulasi
Hitung kadar air keseimbangan Loop waktu Loop lapisan Hitung tekanan uap udara Hitung RH udara Hitung panas laten penguapan jagung
Hitung suhu jagung
Hitung kadar air jagung
Hitung kelembaban mutlak udara
tdk
Total Kedalaman tercapai ?
ya Hitung kadar air rata-rata
Waktu tercapai ?
ya Cetak hasil perhitungan: kadar air, suhu udara,kelembaban mutlak, RH
Selesai
tdk
88
Lampiran 9 Parameter yang digunakan dalam simulasi pengeringan tumpukan Parameter Panas jenis biji jagung Panas jenis udara Panas jenis air Panas jenis uap air Bulk density jagung
Simbol
Nilai
Satuan
Sumber
Cpg
1.527
kJ/kg K
Bala 1997
Cpa Cpl Cpw
1.008 4.186 1.890
kj/kg C kj/kg K kj/kg K
ρd
745
kg/m3
Bala 1997 Bala 1997 Bala 1997 Brooker 1992
Massa jenis udara
ρa
1.11
kg/m
Jansen terjemahan (Arismunandar 1995)
Tekanan atmosfer Laju udara Koeffisien pindah panas konveksi
Patm Ga
101.325 12.7 25.4433 G1.3
kPa kg/mnt-m2
Bala 1997 Pengukuran
kJ/mnt/m3/K
Bala 1997
hcv
3
Panas laten penguapan jagung (Lg) dalam satuan kJ/kg, dihitung dengan persamaan (Strohman&Yeoger 1967 dalam Bala 1997): 1
0.8943 exp
0.1232
.........................................(1)
Kelembaban mutlak (H) dalam satuan kg/kg didapatkan dengan persamaan (Brooker 1974): .
....................................................................................(2)
Kelembaban relatif (RH) dalam satuan %, dihitung dengan persamaan (Brooker 1974): ...........................................................................................(3) Kadar air keseimbangan (Me) dalam satuan % b.k, dihitung dengan persamaan (Brooker et al. 1992) : .
.
.
.
.....................................(4)
Persamaan empiris pengeringan untuk jagung page equation (Hal 1970; Van Rest & Isaacs 1968 dalam Bala 1997): ................................................................................(5)
89
.......................................................................(6) ln
ln ...................................................(7)
Persamaan ini adalah dari ekspresi persamaan garis : .....................................................................................(8) Dimana :
ln
ln Jika persamaan (6) diplot pada pada sumbu x dan y, maka akan didapatkan sebuah garis lurus dengan nilai positif slope u , selanjutnya jika data eksperimen diplotkan dari persamaan (8) maka u dan k didapat dari grafik dengan menarik garis lurus. Hal tersebut ditunjukkan pada gambar berikut:
y2
Slope = u
y1
x1
x2
ln t
Gambar 1 penentuan k dan u dari plot Dalam persamaan yang digunakan, u adalah exponen yang nilainya diberikan dengan persamaan (Bala 1997) : 0.1248
1.2011
1.0685
0.00702
..................(9)
90
Lampiran 10 Kode program Visual Basic untuk simulasi pengeringan tumpukan 'Simulasi Model Persamaan Differensial Parsial 'untuk Pengeringan Biji-bijian Lapisan Tebal pada In-Store Dryer 'Oleh : Diswandi Nurba 'Program Studi Ilmu Keteknikan Pertanian - Sekolah Pasca Sarjana IPB '2008 Dim TG(260), MC(260), AT(260), HM(260), RH(260)
Private Sub Command1_Click() 'Input data AC = Val(Text1.Text) GC = Val(Text2.Text) LC = Val(Text3.Text) WC = Val(Text4.Text) BD = Val(Text5.Text) AP = Val(Text6.Text) Z = Val(Text7.Text) LX = Val(Text8.Text) DT = Val(Text9.Text) DZ = Val(Text10.Text) RP = Val(Text11.Text) EL = 0.01 'batas nilai max absolute selisih RH dg RH udara jenuh menju hitung Suhu Bijian ESS = 0.01 'batas nilai max absolut selisih RH dgn RH udara jenuh menuju hitung delta KA EX = 0.001 'batas min absolut selisih waktu iterasi dg waktu proses menuju waktu total TN = Val(Text15.Text) MO = Val(Text16.Text) / 100 NNi = Val(Text20.Text) * 60 TAd = Val(Text21.Text) RXd = Val(Text22.Text) / 100 GAd = Val(Text23.Text) 'Process TI = TN + 273.15 'SUB HITUNG TEKANAN UAP PS = Exp(52.575 - 6790 / TI - 5.0281 * Log(TI)) AV = 0 For L = 1 To LX MC(L) = MO MEQ = MC(L) RX = Exp(-C0 * (TI ^ C1) * Exp(C2 * (TI ^ C3) * MEQ)) RH(L) = RX P = RH(L) * PS HN = (0.622 * P) / (AP - P) HM(L) = HN AV = AV + MC(L) * DZ Next L AV = AV / Z T = 0: TPR = 60: DTPR = 60 FileNumber = FreeFile() Open "D:\Jagung Drying Result.txt" For Output As #FileNumber For NN = 1 To NNi 'suhu udara dan RH udara lingkungan TA = TAd: RX = RXd TI = TA + 273.15 'Perhitungan tekanan uap dan Kelembaban PS = Exp(52.575 - 6790 / TI - 5.0281 * Log(TI)) H = (0.622 * RX * PS) / (AP - RX * PS) HN = H T = T + DT 'Perhitungan Kadar Air Keseimbangan MEQ = ((Log(1 - RX)) / ((-8.6541 * 10 ^ -5) * (TA + 49.81) * (100 ^ 1.8634))) ^ (1 / 1.8634) 'looping untuk iterasi pada setiap lapisan
91
AV = 0 For L = 1 To LX If RX > 0.98 Then RX = 0.98 TI = TA + 273.15 DC = -0.0098 - 0.19 * RX + 0.01242 * TA U = 0.1248 + 1.2011 * RX - 1.0685 * (RX ^ 2) + 0.00702 * TA
100 200
250
300
If MO = MEQ Then GoTo 100 If MC(L) = MEQ Then GoTo 100 If DC = 0 Then GoTo 100 AZ = Log(Abs(MC(L) - MEQ)) BZ = Log(Abs(MO - MEQ)) TE1 = Abs(BZ - AZ / DC) TE2 = 1 / U TE = TE1 ^ TE2 'Perubahan kadar air DM = -(MC(L) - MEQ) * DC * U * ((TE + (DT / 60)) ^ (U - 1)) * (DT / 60) GoTo 200 DM = 0 A2 = 2 * (TA - TG(L)) B3 = GC + LC * MC(L) F1 = WC * TA + 2501.61 - TG(L) * LC GA = GAd 'Laju massa udara masuk 'Perhitungan Heat Transfer Coeffisien pada tumpukan biji HT = 25.4433 * (GA ^ 1.3) 'Perhitungan Panas laten Penguapan bijian LG = 2.32 * ((1094 - 0.026 * (TA + 17.78)) * (1 + 2.4962 * Exp(-21.73 * MC(L)))) YY = AC * TA + LG - LC * TG(L) XB = DM E1 = AC + WC * (HN - (XB * BD * DZ / (GA * DT))) GE = GA * E1 TP = (BD / DT) * XB TP = A2 + (TP * ((2 * YY / HT) + (F1 * DZ / GE))) BB = B3 + LC * XB BT = 1 + (BD / DT) * (2 * B3 / HT + DZ * BB / GE) DG = TP / BT DA = -(BD * DZ / ((GA * DT) * E1)) * ((DG * BB) - (XB * F1)) TI = TA + DA + 273.15 PS = Exp(52.575 - 6790 / TI - 5.0281 * Log(TI)) H = HN - (XB * BD * DZ / (GA * DT)) P = (H * AP) / (0.622 + H) RX = P / PS FXX = RX - RP If FXX <= EL Then GoTo 1000 If FXX < 0 Then GoTo 1000 If FXX = 0 Then GoTo 1000 XD = 0.0006000001 IT = 0 CX = DM FX = RX - RP CXX = CX + XD XB = CXX 'PERHITUNGAN SUHU UDARA DAN SUHU BIJIAN E1 = AC + WC * (HN - (XB * BD * DZ / (GA * DT))) GE = GA * E1 TP = (BD / DT) * XB TP = A2 + (TP * ((2 * YY / HT) + (F1 * DZ / GE))) BB = B3 + LC * XB BT = 1 + (BD / DT) * (2 * B3 / HT + DZ * BB / GE) DG = TP / BT DA = -(BD * DZ / ((GA * DT) * E1)) * ((DG * BB) - (XB * F1)) TI = TA + DA + 273.15 PS = Exp(52.575 - 6790 / TI - 5.0281 * Log(TI)) H = HN - (XB * BD * DZ / (GA * DT)) P = (H * AP) / (0.622 + H) RX = P / PS FXX = RX - RP If (FX * FXX) < 0 Then GoTo 400 If (FX * FXX) = 0 Then GoTo 900
92
400
500 550 600
CX = CX + XD FX = FXX IT = IT + 1 If (IT - 600) <= 0 Then GoTo 300 For K = 1 To 5 XV = (CX + CXX) / 2 XB = XV E1 = AC + WC * (HN - (XB * BD * DZ / (GA * DT))) GE = GA * E1 TP = (BD / DT) * XB TP = A2 + (TP * ((2 * YY / HT) + (F1 * DZ / GE))) BB = B3 + LC * XB BT = 1 + (BD / DT) * (2 * B3 / HT + DZ * BB / GE) DG = TP / BT DA = -(BD * DZ / ((GA * DT) * E1)) * ((DG * BB) - (XB * F1)) TI = TA + DA + 273.15 PS = Exp(52.575 - 6790 / TI - 5.0281 * Log(TI)) H = HN - (XB * BD * DZ / (GA * DT)) P = (H * AP) / (0.622 + H) RX = P / PS FVE = RX - RP If (FX * FVE) < 0 Then GoTo 500 If (FX * FVE) = 0 Then GoTo 800 CX = XV FX = FVE GoTo 550 CXX = XV FXX = FVE Next K XV = (CX * FXX - CXX * FX) / (FXX - FX) XB = XV E1 = AC + WC * (HN - (XB * BD * DZ / (GA * DT))) GE = GA * E1 TP = (BD / DT) * XB TP = A2 + (TP * ((2 * YY / HT) + (F1 * DZ / GE))) BB = B3 + LC * XB BT = 1 + (BD / DT) * (2 * B3 / HT + DZ * BB / GE) DG = TP / BT DA = -(BD * DZ / ((GA * DT) * E1)) * ((DG * BB) - (XB * F1)) TI = TA + DA + 273.15 PS = Exp(52.575 - 6790 / TI - 5.0281 * Log(TI)) H = HN - (XB * BD * DZ / (GA * DT)) P = (H * AP) / (0.622 + H) RX = P / PS
FVE = RX - RP If Abs(FVE) < ESS Then GoTo 900 If (FX * FVE) < 0 Then GoTo 700 If (FX * FVE) = 0 Then GoTo 900 CX = XV FX = FVE GoTo 600 700 CXX = XV FXX = FVE 800 CXX = XV 900 DM = CXX GoTo 250 1000 TG(L) = TG(L) + DG TA = TA + DA AT(L) = TA MC(L) = MC(L) + DM RH(L) = RX HN = H HM(L) = HN AV = AV + MC(L) * DZ Next L AV = AV / Z TIM = T If Abs(TIM - TPR) >= EX Then GoTo 1100 'Progress Bar
93
baris = 2000 ProgressBar1.Max = baris ProgressBar1.Min = 0 For nilai = 1 To baris ProgressBar1.Value = nilai Next nilai 'Show results in a table With Me.ListView1 .ListItems.Clear .View = lvwReport .ColumnHeaders.Clear .ColumnHeaders.Add , , " Layer ", 700 .ColumnHeaders.Add , , " Moisture Content (%b.k) ", 2100 .ColumnHeaders.Add , , " Air Temp (C) ", 1400 .ColumnHeaders.Add , , " Humidity (kg/kg) ", 1500 .ColumnHeaders.Add , , " RH (%) ", 1300 For i = 1 To LX .ListItems.Add , , i .ListItems(.ListItems.Count).ListSubItems.Add , , Round(MC(i) * 100, 5) .ListItems(.ListItems.Count).ListSubItems.Add , , Round(AT(i), 5) .ListItems(.ListItems.Count).ListSubItems.Add , , Round(HM(i), 5) .ListItems(.ListItems.Count).ListSubItems.Add , , Round(RH(i) * 100, 5) Write #FileNumber, NN, i, MC(i), AT(i), HM(i), RH(i) Next i End With TPR = TPR + DTPR Text17.Text = Str(T) Text18.Text = Str(AV * 100) Text19.Text = Str(MEQ * 100) 'Refresh Text17.Refresh Text18.Refresh 'Text19.Refresh ListView1.Refresh 1100 Next NN Close #FileNumber End Sub Private Sub Command2_Click() Text1.Text = "" Text2.Text = "" Text3.Text = "" Text4.Text = "" Text5.Text = "" Text6.Text = "" Text7.Text = "" Text8.Text = "" Text9.Text = "" Text10.Text = "" Text11.Text = "" Text15.Text = "" Text16.Text = "" Text17.Text = "" Text18.Text = "" Text19.Text = "" Text20.Text = "" Text21.Text = "" Text22.Text = "" Text23.Text = "" Me.ListView1.ListItems.Clear End Sub Private Sub Command3_Click() End End Sub Private Sub Command4_Click() Form2.Show
94
End Sub Private Sub Command5_Click() Form4.Show End Sub Private Sub Command6_Click() 'Reset Text15.Text = "" Text16.Text = "" Text17.Text = "" Text18.Text = "" Text19.Text = "" Text20.Text = "" Text21.Text = "" Text22.Text = "" Text23.Text = "" End Sub
95
Lampiran 11 Interface program simulasi pengeringan tumpukan
96
Lampiran 12 Perhitungan pressure drop dan tekanan statis kipas Tekanan statis pada ruang ISD (plenum chamber) diperlukan untuk memaksa udara mengalir melewati lantai pengering dan massa jagung. Tekanan statis merupakan gaya tegak lurus terhadap dinding saluran dan besarnya tidak tergantung pada kecepatan aliran. Penentuan besarnya tekanan statis meliputi: a. Tahanan dari massa jagung terhadap aliran udara b. Tahanan lantai berlubang terhadap aliran udara c. Penurunan tekanan pada saluran d. Penurunan tekanan karena perubahan luas penampang e. Pengaruh tekanan dinamis Dari pengukuran didapatkan kecepatan aliran udara pada kipas axial yang dipakai pada ISD sebesar 8.20 m/dtk, maka debit dapat dihitung : . Luas kipas dengan diameter 15 inchi adalah : 3.14
3.14 0.19
0.114
Sehingga : 8.20 0.114
0.93
/
--- > 56.09 m3/menit --- > 1981.91 cfm
Bila aliran udara dialirkan melalui lapisan bijian (massa jagung), tahanan dari aliran tersebut menurunkan tekanan (pressure drop) disebabkan oleh kehilangan energi karena gesekan dan turbulensi aliran. Tahanan ini dapat diatasi dengan menambah tekanan pada saluran masuk atau membuat kondisi hampa udara pada saluran keluar. Dengan persamaan (Hukill dan Shedd 1955 dalam Brooker 1992) maka dapat dihitung penurunan tekanan sebagai berikut: ∆
....................................................................... (2)
dimana : ΔP
: penurunan tekanan per ft, (inchi air/ ft)
a dan b
: tetapan, tergantung dari macam bijian.
Qa
: aliran massa udara (cfm/ft2)
97
Tabel 1 Tetapan untuk persamaan (2) Biji-bijian Nilai a Nilai b Barley 0.00069 0.070 Oats 0.00080 0.074 Jagung pipil 0.00065 0.154 Kacang kedelai 0.00036 0.092 Gandum 0.00092 0.050 Sumber : Hukill dan Shedd (1955) dalam Brooker (1992) Luas lantai ISD: Lantai berbentuk lingkaran r = 1.25 m 3.14 1.25
4.91
52.81
sehingga : . .
37.53
/
didapat prossure drop : ∆
0.00065 37.53 ln 1 0.154 37.53
0.48
/
Dengan tebal lapisan jagung = 2.50 m = 8.20 ft, maka tahanan oleh lapisan jagung adalah : ∆
∆
∆
∆
............................................................................ (3) 8.20 0.48
3.94
atau 978.09 Pa (0.98 kPa) Udara yang dihembuskan ke dalam tumpukan biji juga melalui lantai dan plat besi berlubang, maka terjadi penurunan tekanan setelah plat besi tersebut terlewati, Brooker (1992) memberikan persamaan: ∆
.......................................................................... 4
Dari pengukuran didapat pori lantai (Of) = 40% Dari Brooker (1992) Porositas tumpukan jagung ( ε ) = 43% Sehingga tahanan lantai dapat dihitung: ∆
10 9
37.53 0.43 0.4
0.01
Kehilangan tekanan oleh pembebasan saluran ditetapkan dengan persamaan (ASHRAE , 1969) :
98
.............................................................................. 5
∆ dimana :
c1: tetapan pembesaran penampang V: kecepatan sebelum memasuki perubahan penampang (fpm) Luas permukaan kipas, diameter 15 inchi adalah 0.114 1981.91 1.23
1616
1.23
/
dari gambar bentuk saluran berikut (ASHRAE 1969 dalam Brooker 1992) maka c1 = 1.00
Gambar 1. Bentuk saluran pada ISD sehingga : ∆
1
1616 4005
0.16
Maka tekanan total adalah penjumlahan dari ketiga tahanan tersebut: ∆
∆ 3.94
∆
0.01
0.16
4.11
Tekanan dinamis adalah tekanan yang hanya tergantung dari gesekan dan kerapatan udara. Tekanan dinamis dapat dihitung dengan persamaan
.................................................................................... 6 dimana: Pdin
: tekanan dinamis (tekanan kecepatan) (in. air)
V
: kecepatan aliran udara (ft/menit) 1616 4005
0.16
Sehingga tekanan statis dihitung dengan persamaan :
............................................................................. 7 4.11
0.16
Pstatis = 0.98 kPa
3.95
99
Tekanan statis yang dibutuhkan untuk pengeringan jagung dengan tumpukan setinggi 2.50 m (7500 kg) dalam ISD adalah 0.98 kPa atau 982.2 Pa. Sehingga tekanan kipas axial sebesar 90 Pa masih kurang untuk melakukan kerjanya. Namun penataan pipa saluran udara sebanyak 13 pipa yang menyediakan ruang bebas untuk pergerakan udara dalam ISD dapat membantu kinerja kipas axial tersebut. Tabel 2 Faktor konversi (Brooker 1992) 1 cfm 1 cfm/bu 1 cfm/cwt 1 cfm/ft2 1 cfm/ft3 1 cfm/ton 1 lb/hr ft2 1 ft3/hr 1 atm 1 bar 1 in. Hg 1 in. water 1 in.water/ft 1 lb/ft2 1 lb in2 1 mm Hg
Airflow = 0.0283 m3/min = 1.114 m3/min.t = 0.804 m3/m3.min = 0.6238 m3/min.t = 0.3048 m3/m2.min = 1 m3/m3.min = 0.0311 m3/min.t = 4.883 kg/hr.m2 = 0.472 x 10-3 m3/min Pressure = 101.3 kPa = 100 kPa = 3.376 kPa = 248.8 Pa = 816.40 Pa/m = 47.88 Pa = 6.895 kPa = 0.133kPa
Area and Length = 0.3048 m = 0.0254 m = 4046.87 m2 = 0.404687 ha 1 ft2 = 0.0929 m2 1 in.2 = 6.452 cm2 1 mi = 1.6093 km Volume 1 ft3 = 0.02832 m3 1 gal (U.S) = 0.00379 m3 1 in.3 = 16.387 cm3 1L = 0.001 m3 1 bu = 0.03524 m3 1 bu/acre = 0.0871 m3/ha 1 bu/hr = 0.03524 m3/hr 1 yd3 = 0.76455 m3
1 ft 1 in 1 acre
100
Lampiran 13 Lokasi titik pengukuran suhu dan kecepatan udara di dalam ruangan ISD
101
Lampiran 14 Lokasi titik pengambilan sampel pengukuran kadar air
102
Lampiran 15 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 41 42 43 44 45 46 47 48 49 50 51 52 53 54
X (m) 0 0.6 0.8 0.4 0.6 0.2 -0.6 -0.8 -0.4 -0.6 -0.2 -0.76 0 0.76 0.76 0.76 0 -0.76 -0.76 -0.2 0.8 0.2 -0.8 -0.2 0.2 0.2 -0.2 0 0.6 0.8 0.4 0.6 0.2 -0.6 -0.8 -0.4 -0.6 -0.2 -0.76 0 0.76 0.76 0.76 0 -0.76 -0.76 -0.2 0.8 0.2 -0.8 -0.2 0.2 0.2 -0.2
Koordinat Y (m) 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25 1.25
Z (m) 0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2 0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2
T(C)
V (m/s)
RH
33.7 33.5 33.5 33.6 33.6 32.2 33.6 33.6 33.6 33.6 33.4 27.5 28.3 27.5 28.2 27.5 28.3 27.5 28.3 32.7 33.3 33.6 33.5 33.7 33.6 33.6 33.6 33.7 33.5 33.5 33.5 33.6 33.6 33.6 33.6 33.6 33.6 33.7 27.7 28.8 27.9 28.5 27.9 29.0 27.7 28.5 33.5 33.6 33.6 33.5 33.6 33.6 33.6 33.7
1.02 0.41 0.49 0.43 0.36 0.36 0.47 0.51 0.36 0.59 0.96 0.15 0.09 0.15 0.09 0.11 0.20 0.10 0.10 0.33 0.48 0.33 0.45 0.35 0.75 0.65 0.34 0.90 0.41 0.52 0.48 0.28 0.35 0.37 0.50 0.36 0.56 1.43 0.13 0.05 0.12 0.10 0.07 0.32 0.09 0.12 0.35 0.36 0.33 0.59 0.42 0.53 0.56 0.62
60.5 60.9 61.1 60.8 60.7 65.5 60.6 60.6 60.6 60.6 61.5 85.9 81.8 85.9 82.4 86.0 81.9 86.0 82.1 63.7 61.5 60.7 60.8 60.5 60.5 60.6 60.6 60.5 60.8 61.1 61.0 60.7 60.8 60.5 60.6 60.6 60.5 60.5 85.0 79.8 83.9 80.9 84.0 78.9 85.1 81.1 60.8 60.8 60.7 60.9 60.5 60.6 60.7 60.5
Rata-rata & St. Deviasi
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
67.9 10.7 31.8 2.6 0.39 0.25
% % C C m/dtk m/dtk
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
67.1 10.2 32.0 2.5 0.40 0.29
% % C C m/dtk m/dtk
103
Lampiran 15 (lanjutan) 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
0 0.6 0.8 0.4 0.6 0.2 -0.6 -0.8 -0.4 -0.6 -0.2 -0.76 0 0.76 0.76 0.76 0 -0.76 -0.76 -0.2 0.8 0.2 -0.8 -0.2 0.2 0.2 -0.2 0 0.6 0.8 0.4 0.6 0.2 -0.6 -0.8 -0.4 -0.6 -0.2 -0.76 0 0.76 0.76 0.76 0 -0.76 -0.76 -0.2 0.8 0.2 -0.8 -0.2 0.2 0.2 -0.2
1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25
0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2 0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2
33.5 33.4 33.5 33.6 33.6 33.4 33.6 33.5 33.6 33.6 33.5 27.6 28.3 27.6 28.3 27.6 28.4 27.6 28.4 33.6 33.6 33.4 33.6 33.6 33.6 33.5 33.6 33.5 33.4 33.5 33.6 33.6 33.5 33.6 33.5 33.6 33.5 33.5 27.6 28.4 27.6 28.5 27.6 28.5 27.6 28.5 33.6 33.6 33.5 33.5 33.6 33.6 33.5 33.5
0.41 0.42 0.41 0.36 0.26 0.20 0.31 0.42 0.62 0.45 0.54 0.11 0.03 0.12 0.07 0.04 0.11 0.05 0.10 0.32 0.22 0.33 0.53 0.54 0.34 0.29 0.98 0.27 0.14 0.36 0.37 0.25 0.23 0.42 0.37 0.64 0.41 0.32 0.09 0.05 0.10 0.05 0.03 0.11 0.02 0.06 0.38 0.10 0.38 0.46 0.49 0.33 0.32 0.45
61.1 61.3 60.9 60.7 60.8 61.4 60.5 60.9 60.5 60.6 61.1 85.3 81.9 85.3 81.9 85.4 81.7 85.4 81.7 60.7 60.7 61.2 60.8 60.7 60.8 61.1 60.6 61.1 61.3 60.8 60.6 60.7 60.8 60.5 61.0 60.7 60.8 61.2 85.4 81.3 85.4 80.9 85.5 81.3 85.5 81.2 60.6 60.6 61.1 61.0 60.8 60.8 60.9 60.9
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
67.6 10.6 31.9 2.6 0.32 0.22
% % C C m/dtk m/dtk
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
67.5 10.5 31.9 2.6 0.27 0.17
% % C C m/dtk m/dtk
104
Lampiran 15 (lanjutan) 109 0 110 0.6 111 0.8 112 0.4 113 0.6 114 0.2 115 -0.6 116 -0.8 117 -0.4 118 -0.6 119 -0.2 120 -0.76 121 0 122 0.76 123 0.76 124 0.76 125 0 126 -0.76 127 -0.76 128 -0.2 129 0.8 130 0.2 131 -0.8 132 -0.2 133 0.2 134 0.2 135 -0.2 Rata-rata Standar Deviasi
2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75 2.75
0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2
33.5 33.5 33.6 33.6 33.6 33.6 33.6 33.5 33.5 33.5 33.5 27.5 27.8 27.5 27.9 27.5 27.8 27.5 27.8 33.7 33.6 33.5 33.6 33.5 33.6 33.5 33.5 31.9 2.60
0.47 0.25 0.18 0.31 0.41 0.26 0.40 0.24 0.45 0.21 0.44 0.05 0.03 0.05 0.02 0.01 0.03 0.01 0.03 0.29 0.20 0.41 0.29 0.21 0.42 0.32 0.37 0.32 0.22
61.0 61.2 60.6 60.6 60.7 60.6 60.6 60.9 60.8 61.1 61.1 86.0 84.4 86.0 84.2 86.1 84.5 86.2 84.4 60.4 60.6 61.0 60.8 60.9 60.6 60.9 60.9 67.6 10.7
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
68.0 11.4 31.8 2.8 0.23 0.16
% % C C m/dtk m/dtk
105
Lampiran 16 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 41 42 43 44 45 46 47 48 49 50 51 52 53 54
Koordinat y 0 0.75 0.6 0.75 0.8 0.75 0.4 0.75 0.6 0.75 0.2 0.75 -0.6 0.75 -0.8 0.75 -0.4 0.75 -0.6 0.75 -0.2 0.75 -0.76 0.75 0 0.75 0.76 0.75 0.76 0.75 0.76 0.75 0 0.75 -0.76 0.75 -0.76 0.75 -0.2 0.75 0.8 0.75 0.2 0.75 -0.8 0.75 -0.2 0.75 0.2 0.75 0.2 0.75 -0.2 0.75 0 1.25 0.6 1.25 0.8 1.25 0.4 1.25 0.6 1.25 0.2 1.25 -0.6 1.25 -0.8 1.25 -0.4 1.25 -0.6 1.25 -0.2 1.25 -0.76 1.25 0 1.25 0.76 1.25 0.76 1.25 0.76 1.25 0 1.25 -0.76 1.25 -0.76 1.25 -0.2 1.25 0.8 1.25 0.2 1.25 -0.8 1.25 -0.2 1.25 0.2 1.25 0.2 1.25 -0.2 1.25 x
z 0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2 0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2
T(C )
V (m/s)
RH
33.7 33.6 33.6 33.7 33.7 33.7 33.7 33.7 33.7 33.7 33.8 33.8 33.7 33.8 33.6 33.7 33.7 33.8 33.7 33.7 33.6 33.7 33.7 33.7 33.7 33.7 33.7 33.7 33.7 33.6 33.6 33.7 33.7 33.7 33.7 33.7 33.7 33.7 33.8 33.7 33.8 33.6 33.7 33.7 33.7 33.7 33.7 33.6 33.7 33.7 33.7 33.6 33.6 33.7
1.08 0.40 0.42 0.45 0.35 0.35 0.41 0.48 0.43 0.58 1.07 1.14 0.42 0.99 0.41 0.77 0.77 0.90 0.49 0.34 0.45 0.30 0.48 0.64 0.65 0.52 0.62 0.74 0.43 0.44 0.47 0.31 0.26 0.25 0.50 0.46 0.59 1.42 0.88 0.30 0.84 0.39 0.56 0.43 0.58 0.54 0.28 0.38 0.32 0.63 0.55 0.43 0.50 0.86
60.2 60.5 60.6 60.5 60.5 60.3 60.2 60.2 60.3 60.3 60.1 60.0 60.2 60.1 60.7 60.3 60.2 60.1 60.3 60.4 60.7 60.4 60.3 60.2 60.2 60.4 60.2 60.3 60.4 60.6 60.7 60.3 60.4 60.2 60.3 60.2 60.2 60.2 60.1 60.3 60.1 60.7 60.4 60.3 60.2 60.4 60.4 60.7 60.5 60.4 60.2 60.6 60.6 60.2
Rata-rata & St. Deviasi
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
60.3 0.2 33.7 0.1 0.59 0.25
% % C C m/dtk m/dtk
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
60.4 0.2 33.7 0.1 0.53 0.25
% % C C m/dtk m/dtk
106
Lampiran 16 (lanjutan) 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108
0 0.6 0.8 0.4 0.6 0.2 -0.6 -0.8 -0.4 -0.6 -0.2 -0.76 0 0.76 0.76 0.76 0 -0.76 -0.76 -0.2 0.8 0.2 -0.8 -0.2 0.2 0.2 -0.2 0 0.6 0.8 0.4 0.6 0.2 -0.6 -0.8 -0.4 -0.6 -0.2 -0.76 0 0.76 0.76 0.76 0 -0.76 -0.76 -0.2 0.8 0.2 -0.8 -0.2 0.2 0.2 -0.2
1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 1.75 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25 2.25
0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2 0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2
33.6 33.6 33.6 33.6 33.7 33.7 33.7 33.7 33.7 33.7 33.6 33.8 33.7 33.7 33.6 33.6 33.6 33.7 33.7 33.7 33.7 33.5 33.7 33.6 33.6 33.6 33.7 33.6 33.5 33.6 33.6 33.7 33.7 33.7 33.7 33.7 33.6 33.5 33.7 33.7 33.7 33.6 33.5 33.5 33.6 33.6 33.7 33.7 33.5 33.7 33.6 33.7 33.6 33.6
0.44 0.42 0.42 0.33 0.30 0.20 0.33 0.43 0.65 0.43 0.67 0.83 0.12 0.84 0.40 0.27 0.28 0.35 0.52 0.30 0.40 0.40 0.53 0.44 0.35 0.30 0.96 0.19 0.30 0.24 0.24 0.29 0.24 0.41 0.38 0.63 0.34 0.54 0.74 0.20 0.76 0.24 0.25 0.15 0.14 0.36 0.29 0.40 0.43 0.43 0.33 0.36 0.20 0.60
60.7 60.7 60.7 60.8 60.3 60.3 60.2 60.3 60.2 60.3 60.6 60.1 60.3 60.2 60.7 60.5 60.8 60.3 60.4 60.4 60.5 60.9 60.4 60.5 60.5 60.7 60.3 60.8 60.9 60.6 60.6 60.2 60.3 60.2 60.5 60.3 60.6 60.9 60.1 60.3 60.2 60.6 61.1 60.9 60.7 60.5 60.2 60.5 60.9 60.5 60.5 60.5 60.7 60.6
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
60.5 0.2 33.7 0.1 0.44 0.20
% % C C m/dtk m/dtk
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
60.5 0.3 33.6 0.1 0.36 0.17
% % C C m/dtk m/dtk
107
Lampiran 16 (lanjutan) 109 0 2.75 110 0.6 2.75 111 0.8 2.75 112 0.4 2.75 113 0.6 2.75 114 0.2 2.75 115 -0.6 2.75 116 -0.8 2.75 117 -0.4 2.75 118 -0.6 2.75 119 -0.2 2.75 120 -0.76 2.75 121 0 2.75 122 0.76 2.75 123 0.76 2.75 124 0.76 2.75 125 0 2.75 126 -0.76 2.75 127 -0.76 2.75 128 -0.2 2.75 129 0.8 2.75 130 0.2 2.75 131 -0.8 2.75 132 -0.2 2.75 133 0.2 2.75 134 0.2 2.75 135 -0.2 2.75 Rata-rata Rata standar deviasi
0.4 0.6 0.2 0 -0.6 -0.8 -0.6 -0.2 0 0.6 0.8 -0.76 -0.76 -0.76 0 0.76 0.76 0.76 0 -0.8 -0.2 0.8 0.2 -0.2 -0.2 0.2 0.2
33.6 33.6 33.6 33.6 33.7 33.7 33.7 33.6 33.6 33.6 33.5 33.7 33.7 33.7 33.7 33.5 33.6 33.5 33.7 33.7 33.7 33.5 33.7 33.7 33.7 33.6 33.7 33.7 0.10
0.68 0.30 0.38 0.44 0.45 0.26 0.39 0.27 0.38 0.20 0.37 0.58 0.24 0.57 0.33 0.18 0.37 0.12 0.19 0.28 0.29 0.45 0.25 0.38 0.36 0.58 0.60 0.46 0.20
60.6 60.5 60.5 60.6 60.2 60.2 60.3 60.5 60.5 60.8 60.8 60.2 60.2 60.2 60.4 61.0 60.8 61.0 60.5 60.2 60.2 60.9 60.4 60.5 60.4 60.6 60.4 60.4 0.20
Rata RH St Dev RH Rata T St.Dev T Rata V St.Dev V
: : : : : :
60.5 0.2 33.6 0.1 0.37 0.14
% % C C m/dtk m/dtk
108
Lampiran 17 Perbandingan keragaman suhu, kecepatan udara dan RH di dalam ISD pada kedua simulasi CFD.
Keragaman suhu Ketinggian (m)
0.75
Rata-rata (oC) Standar deviasi (oC)
31.8 2.6
Rata-rata (oC) Standar deviasi (oC)
33.7 0.05
1.25 1.75 Skenario 1 32 31.9 2.5 2.6 Skenario 2 33.7 33.7 0.05 0.06
2.25
2.75
Rata-rata
31.9 2.6
31.8 2.8
31.9 2.6
33.6 0.08
33.6 0.07
33.7 0.06
Keragaman kecepatan aliran udara Ketinggian (m) Rata-rata (m/dtk) Standar deviasi (m/dtk ) Rata-rata (m/dtk) Standar deviasi (m/dtk )
0.75 1.25 Skenario 1 0.39 0.40 0.2 0.29 Skenario 2 0.59 0.53 0.25 0.25
1.75
2.25
2.75
Rata-rata
0.32 0.22
0.27 0.17
0.23 0.16
0.32 0.22
0.44 0.20
0.36 0.17
0.37 0.14
0.46 0.20
Keragaman RH Ketinggian (m)
0.75
Rata-rata (% ) Standar deviasi ( % )
67.9 10.7
Rata-rata (% ) Standar deviasi ( % )
60.3 0.2
1.25 1.75 Skenario 1 67.1 67.6 10.2 10.6 Skenario 2 60.4 60.5 0.2 0.2
2.25
2.75
Rata-rata
67.5 10.5
68 11.4
67.6 10.7
60.5 0.3
60.5 0.2
60.4 0.2
109
Lampiran 18 Data validasi suhu dan kecepatan udara hasil pengukuran dan simulasi CFD serta nilai error dan standar deviasinya Poin
Simbol
1 A1 2 A2 3 A3 4 A4 5 B1 6 B2 7 B3 8 B4 9 C1 10 C2 11 C3 12 C4 13 D1 14 D2 15 D3 16 D4 17 E1 18 E2 19 E3 20 E4 max min Rata-rata Total Error Korelasi ( R )
x -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2
Koordinat y 0.75 0.75 0.75 0.75 1.25 1.25 1.25 1.25 1.75 1.75 1.75 1.75 2.25 2.25 2.25 2.25 2.75 2.75 2.75 2.75
z -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8
T-ukur 34.5 33.4 33.9 33.6 34.5 34.3 33.7 34.7 34.1 33.7 34.2 33.8 34.0 34.2 34.3 33.7 34.2 34.6 34.4 34.0 34.7 33.4 34.1
Suhu ( C ) T-CFD 33.62 32.22 33.48 33.35 33.61 33.56 33.46 33.65 33.51 33.37 33.53 33.48 33.49 33.54 33.55 33.45 33.52 33.63 33.60 33.48 33.65 32.22 33.46
Std 0.60 0.84 0.29 0.15 0.61 0.53 0.19 0.74 0.39 0.22 0.49 0.24 0.38 0.49 0.53 0.18 0.45 0.69 0.55 0.36 0.84 0.15 0.45
Error 0.85 1.19 0.41 0.22 0.86 0.76 0.27 1.04 0.55 0.32 0.69 0.35 0.53 0.70 0.75 0.26 0.64 0.97 0.78 0.50 1.19 0.22 0.63 12.64 0.66
v-ukur 0.25 0.23 0.18 0.43 0.03 0.25 0.06 0.45 0.10 0.08 0.03 0.36 0.11 0.05 0.03 0.21 0.19 0.05 0.05 0.23 0.45 0.03 0.17
Kecepatan udara (m/s) v-CFD SD 0.51 0.18 0.36 0.09 0.49 0.22 0.96 0.38 0.50 0.33 0.35 0.07 0.52 0.33 1.43 0.69 0.42 0.23 0.20 0.08 0.41 0.27 0.54 0.13 0.37 0.19 0.23 0.13 0.36 0.23 0.32 0.08 0.24 0.03 0.26 0.15 0.18 0.09 0.44 0.15 1.43 0.69 0.18 0.03 0.45 0.20
Error 0.26 0.13 0.31 0.53 0.47 0.10 0.46 0.97 0.32 0.12 0.38 0.18 0.26 0.18 0.33 0.12 0.04 0.21 0.13 0.21 0.97 0.04 0.29 5.7 0.73
110
Lampiran 19 Data validasi RH udara hasil pengukuran dan perhitungan serta nilai error dan standar deviasinya. Poin
Simbol
1 A1 2 A2 3 A3 4 A4 5 B1 6 B2 7 B3 8 B4 9 C1 10 C2 11 C3 12 C4 13 D1 14 D2 15 D3 16 D4 17 E1 18 E2 19 E3 20 E4 max min Rata-rata Total Error Korelasi
x -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2
Koordinat y 0.75 0.75 0.75 0.75 1.25 1.25 1.25 1.25 1.75 1.75 1.75 1.75 2.25 2.25 2.25 2.25 2.75 2.75 2.75 2.75
z -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8 -0.2 -0.8 0.2 0.8
RH-ukur 61.00 64.72 63.03 64.15 61.00 61.53 63.59 60.27 62.41 63.73 61.87 63.24 62.55 61.80 61.60 63.66 62.07 60.60 61.33 62.69 64.72 60.27 62.34
RH ( % ) RH-hitung 60.57 65.52 61.07 61.49 60.61 60.77 61.12 60.47 60.94 61.43 60.89 61.06 61.02 60.85 60.82 61.16 60.91 60.56 60.65 61.06 65.52 60.47 61.15
SD 0.30 0.57 1.39 1.88 0.27 0.54 1.74 0.14 1.04 1.63 0.69 1.54 1.08 0.67 0.55 1.77 0.82 0.02 0.48 1.15 1.88 0.02 0.91
Error 0.42 0.80 1.96 2.66 0.39 0.76 2.47 0.20 1.47 2.30 0.98 2.18 1.52 0.95 0.77 2.50 1.16 0.03 0.68 1.63 2.66 0.03 1.29 25.85 0.66
111
Lampiran 20 Perubahan kadar air setelah simulasi pengeringan selama 150 jam pada Simulasi 1. T in RH Waktu Me Laju massa udara Waktu (jam ke-) 0 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 41 42 43 44 45
Layer 1 18.00 18.85 18.76 18.67 18.59 18.51 18.43 18.36 18.29 18.23 18.17 18.11 18.05 18.00 17.94 17.89 17.85 17.80 17.76 17.72 17.67 17.64 17.60 17.56 17.53 17.50 17.46 17.43 17.41 17.38 17.35 17.32 17.30 17.28 17.25 17.23 17.21 17.19 17.17 17.15 17.13 17.11 17.10 17.08 17.06 17.05
: : : : : Layer 10 18.00 19.41 19.30 19.20 19.10 19.00 18.91 18.82 18.74 18.66 18.58 18.51 18.44 18.37 18.31 18.24 18.18 18.13 18.07 18.02 17.97 17.92 17.87 17.82 17.78 17.74 17.70 17.66 17.62 17.59 17.55 17.52 17.49 17.46 17.43 17.40 17.37 17.35 17.32 17.30 17.27 17.25 17.23 17.21 17.19 17.17
Layer 40 18.00 19.48 19.40 19.32 19.24 19.16 19.09 19.01 18.94 18.87 18.80 18.73 18.66 18.60 18.54 18.47 18.42 18.36 18.30 18.25 18.20 18.15 18.10 18.05 18.00 17.96 17.91 17.87 17.83 17.79 17.75 17.72 17.68 17.64 17.61 17.58 17.55 17.52 17.49 17.46 17.43 17.40 17.38 17.35 17.33 17.31
31 oC 73 % 150 jam 16.6 % b.k 2 12.7 kg/menit-m Kadar air (% b.k) Layer Layer Layer 50 100 150 18.00 18.00 18.00 19.50 19.50 18.00 19.44 19.52 19.49 19.37 19.49 19.52 19.29 19.44 19.50 19.22 19.39 19.46 19.14 19.34 19.43 19.07 19.28 19.39 19.00 19.23 19.35 18.94 19.18 19.31 18.87 19.13 19.28 18.80 19.08 19.24 18.74 19.03 19.20 18.68 18.97 19.16 18.61 18.92 19.12 18.55 18.87 19.08 18.50 18.82 19.04 18.44 18.77 19.00 18.38 18.72 18.96 18.33 18.67 18.92 18.28 18.63 18.88 18.23 18.58 18.84 18.18 18.53 18.80 18.13 18.48 18.76 18.08 18.44 18.72 18.04 18.39 18.68 17.99 18.35 18.64 17.95 18.30 18.60 17.91 18.26 18.56 17.86 18.22 18.52 17.83 18.17 18.48 17.79 18.13 18.44 17.75 18.09 18.40 17.71 18.05 18.36 17.68 18.01 18.32 17.64 17.98 18.29 17.61 17.94 18.25 17.58 17.90 18.21 17.55 17.87 18.17 17.52 17.83 18.14 17.49 17.80 18.10 17.46 17.76 18.07 17.44 17.73 18.03 17.41 17.70 18.00 17.38 17.67 17.96 17.36 17.64 17.93
Layer 200 18.00 18.00 19.53 19.52 19.52 19.50 19.48 19.45 19.42 19.39 19.37 19.34 19.31 19.28 19.25 19.22 19.19 19.15 19.12 19.09 19.06 19.03 18.99 18.96 18.93 18.89 18.86 18.83 18.79 18.76 18.73 18.69 18.66 18.62 18.59 18.55 18.52 18.48 18.45 18.42 18.38 18.35 18.31 18.28 18.25 18.21
Layer 250 18.00 18.00 18.05 19.54 19.56 19.59 19.58 19.56 19.54 19.52 19.49 19.47 19.45 19.42 19.40 19.38 19.35 19.33 19.30 19.27 19.25 19.22 19.20 19.17 19.14 19.11 19.08 19.06 19.03 19.00 18.97 18.94 18.91 18.88 18.85 18.82 18.79 18.76 18.73 18.70 18.66 18.63 18.60 18.57 18.54 18.51
Rata-rata 18.00 18.84 19.19 19.33 19.28 19.23 19.17 19.12 19.07 19.01 18.96 18.91 18.86 18.81 18.76 18.71 18.67 18.62 18.58 18.53 18.49 18.45 18.41 18.37 18.33 18.29 18.25 18.21 18.18 18.14 18.10 18.07 18.03 18.00 17.97 17.94 17.90 17.87 17.84 17.81 17.78 17.75 17.73 17.70 17.67 17.65
112
Lampiran 20 (lanjutan) 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99
17.03 17.02 17.00 16.99 16.98 16.97 16.95 16.94 16.93 16.92 16.91 16.90 16.89 16.88 16.87 16.86 16.85 16.85 16.84 16.83 16.82 16.82 16.81 16.80 16.80 16.79 16.78 16.78 16.77 16.77 16.76 16.75 16.75 16.74 16.74 16.74 16.73 16.73 16.72 16.72 16.71 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.68 16.68 16.68 16.67 16.67
17.15 17.13 17.11 17.10 17.08 17.06 17.05 17.03 17.02 17.01 16.99 16.98 16.97 16.96 16.94 16.93 16.92 16.91 16.90 16.89 16.88 16.87 16.86 16.86 16.85 16.84 16.83 16.82 16.82 16.81 16.80 16.80 16.79 16.78 16.78 16.77 16.77 16.76 16.76 16.75 16.75 16.74 16.74 16.73 16.73 16.72 16.72 16.72 16.71 16.71 16.70 16.70 16.69 16.69
17.28 17.26 17.24 17.22 17.20 17.18 17.16 17.14 17.13 17.11 17.09 17.08 17.06 17.05 17.03 17.02 17.01 16.99 16.98 16.97 16.96 16.95 16.94 16.93 16.91 16.90 16.90 16.89 16.88 16.87 16.86 16.85 16.84 16.84 16.83 16.82 16.81 16.81 16.80 16.79 16.79 16.78 16.78 16.77 16.77 16.76 16.75 16.75 16.74 16.74 16.74 16.73 16.72 16.72
17.33 17.31 17.29 17.27 17.25 17.23 17.21 17.19 17.17 17.15 17.13 17.12 17.10 17.08 17.07 17.05 17.04 17.03 17.01 17.00 16.99 16.98 16.96 16.95 16.94 16.93 16.92 16.91 16.90 16.89 16.88 16.87 16.86 16.86 16.85 16.84 16.83 16.83 16.82 16.81 16.80 16.80 16.79 16.79 16.78 16.77 16.77 16.76 16.76 16.75 16.75 16.74 16.73 16.73
17.61 17.58 17.55 17.52 17.50 17.47 17.44 17.42 17.40 17.37 17.35 17.33 17.31 17.28 17.26 17.24 17.23 17.21 17.19 17.17 17.15 17.14 17.12 17.11 17.09 17.08 17.06 17.05 17.03 17.02 17.01 17.00 16.98 16.97 16.96 16.95 16.94 16.93 16.92 16.91 16.90 16.89 16.88 16.87 16.86 16.86 16.85 16.84 16.83 16.83 16.82 16.81 16.80 16.80
17.90 17.86 17.83 17.80 17.77 17.74 17.71 17.68 17.65 17.62 17.60 17.57 17.54 17.52 17.49 17.47 17.45 17.42 17.40 17.38 17.36 17.33 17.31 17.29 17.27 17.25 17.24 17.22 17.20 17.18 17.17 17.15 17.13 17.12 17.10 17.09 17.07 17.06 17.05 17.03 17.02 17.01 17.00 16.98 16.97 16.96 16.95 16.94 16.93 16.92 16.91 16.90 16.88 16.88
18.18 18.15 18.11 18.08 18.05 18.02 17.99 17.95 17.92 17.89 17.86 17.83 17.80 17.78 17.75 17.72 17.69 17.67 17.64 17.61 17.59 17.56 17.54 17.51 17.49 17.47 17.44 17.42 17.40 17.38 17.36 17.34 17.32 17.30 17.28 17.26 17.24 17.22 17.21 17.19 17.17 17.16 17.14 17.13 17.11 17.10 17.08 17.07 17.06 17.04 17.03 17.02 16.99 16.99
18.47 18.44 18.41 18.38 18.35 18.32 18.28 18.25 18.22 18.19 18.16 18.13 18.10 18.07 18.04 18.01 17.98 17.95 17.92 17.89 17.86 17.83 17.81 17.78 17.75 17.72 17.70 17.67 17.65 17.62 17.60 17.57 17.55 17.52 17.50 17.48 17.46 17.43 17.41 17.39 17.37 17.35 17.33 17.31 17.29 17.27 17.26 17.24 17.22 17.20 17.19 17.17 17.14 17.14
17.62 17.59 17.57 17.54 17.52 17.50 17.47 17.45 17.43 17.41 17.39 17.37 17.35 17.33 17.31 17.29 17.27 17.25 17.23 17.22 17.20 17.18 17.17 17.15 17.14 17.12 17.11 17.09 17.08 17.07 17.05 17.04 17.03 17.02 17.00 16.99 16.98 16.97 16.96 16.95 16.94 16.93 16.92 16.91 16.90 16.89 16.88 16.88 16.87 16.86 16.85 16.84 16.83 16.83
113
Lampiran 20 (lanjutan) 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150
16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.65 16.65 16.65 16.65 16.65 16.64 16.64 16.64 16.64 16.64 16.64 16.63 16.63 16.63 16.63 16.63 16.63 16.63 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.60 16.60 16.60 16.60 16.60
16.69 16.69 16.69 16.68 16.68 16.68 16.67 16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.65 16.65 16.65 16.65 16.65 16.65 16.64 16.64 16.64 16.64 16.64 16.64 16.63 16.63 16.63 16.63 16.63 16.63 16.63 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61 16.61
16.72 16.71 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.68 16.68 16.68 16.68 16.67 16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.65 16.65 16.65 16.65 16.65 16.64 16.64 16.64 16.64 16.64 16.64 16.63 16.63 16.63 16.63 16.63 16.63 16.63 16.63 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.61
16.73 16.73 16.72 16.72 16.71 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.68 16.68 16.68 16.68 16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.66 16.65 16.65 16.65 16.65 16.65 16.64 16.64 16.64 16.64 16.64 16.64 16.63 16.63 16.63 16.63 16.63 16.63 16.63 16.62 16.62 16.62 16.62 16.62 16.62 16.62 16.62
16.79 16.79 16.78 16.78 16.77 16.76 16.76 16.75 16.75 16.74 16.74 16.74 16.73 16.73 16.72 16.72 16.72 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.68 16.68 16.68 16.68 16.67 16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.66 16.65 16.65 16.65 16.65 16.65 16.64 16.64 16.64 16.64 16.64 16.64 16.63 16.63
16.88 16.87 16.86 16.85 16.84 16.84 16.83 16.82 16.82 16.81 16.80 16.80 16.79 16.78 16.78 16.77 16.77 16.76 16.76 16.75 16.75 16.74 16.74 16.73 16.73 16.73 16.72 16.72 16.71 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.68 16.68 16.68 16.68 16.67 16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.65 16.65
16.98 16.97 16.96 16.95 16.94 16.93 16.92 16.91 16.90 16.89 16.88 16.88 16.87 16.86 16.85 16.84 16.84 16.83 16.82 16.82 16.81 16.80 16.80 16.79 16.78 16.78 16.77 16.77 16.76 16.76 16.75 16.75 16.74 16.74 16.73 16.73 16.73 16.72 16.72 16.71 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.69 16.68 16.68
17.13 17.11 17.10 17.08 17.07 17.06 17.04 17.03 17.02 17.01 17.00 16.98 16.97 16.96 16.95 16.94 16.93 16.92 16.91 16.90 16.89 16.89 16.88 16.87 16.86 16.85 16.85 16.84 16.83 16.82 16.82 16.81 16.81 16.80 16.79 16.79 16.78 16.78 16.77 16.77 16.76 16.76 16.75 16.75 16.74 16.74 16.73 16.73 16.72 16.72 16.72
16.82 16.82 16.81 16.80 16.80 16.79 16.79 16.78 16.77 16.77 16.76 16.76 16.75 16.75 16.74 16.74 16.74 16.73 16.73 16.72 16.72 16.71 16.71 16.71 16.70 16.70 16.70 16.69 16.69 16.69 16.68 16.68 16.68 16.68 16.67 16.67 16.67 16.67 16.66 16.66 16.66 16.66 16.66 16.65 16.65 16.65 16.65 16.65 16.64 16.64 16.64
114
Lampiran 21 Perubahan kadar air setelah simulasi pengeringan selama 120 jam pada Simulasi 2. T in RH Waktu Me Laju massa udara Waktu (jam Layer ke-) 1 0 18.00 1 18.11 2 17.75 3 17.43 4 17.13 5 16.86 6 16.62 7 16.39 8 16.18 9 15.99 10 15.82 11 15.66 12 15.51 13 15.37 14 15.24 15 15.13 16 15.02 17 14.92 18 14.83 19 14.74 20 14.66 21 14.58 22 14.51 23 14.45 24 14.39 25 14.33 26 14.28 27 14.23 28 14.19 29 14.14 30 14.10 31 14.07 32 14.03 33 14.00 34 13.97 35 13.94 36 13.91 37 13.89 38 13.87 39 13.84 40 13.82 41 13.80 42 13.78 43 13.77 44 13.75 45 13.74
: : : : : Layer 10 18.00 19.12 18.74 18.38 18.04 17.73 17.44 17.18 16.93 16.70 16.48 16.28 16.10 15.92 15.76 15.61 15.48 15.35 15.23 15.12 15.01 14.92 14.83 14.74 14.66 14.59 14.52 14.46 14.40 14.34 14.29 14.24 14.20 14.15 14.11 14.08 14.04 14.01 13.98 13.95 13.92 13.90 13.87 13.85 13.83 13.81
Layer 40 18.00 19.40 19.21 18.98 18.76 18.54 18.32 18.10 17.89 17.68 17.48 17.28 17.08 16.90 16.72 16.55 16.38 16.22 16.07 15.93 15.79 15.66 15.54 15.42 15.31 15.20 15.10 15.01 14.92 14.84 14.76 14.69 14.62 14.55 14.49 14.43 14.38 14.32 14.28 14.23 14.19 14.15 14.11 14.07 14.04 14.01
33 oC 59.8 % 120 jam 13.5 % b.k 2 12.7 kg/mnt-m Kadar air (% b.k) Layer Layer 50 100 18.00 18.00 19.39 19.42 19.26 19.48 19.08 19.43 18.89 19.34 18.71 19.26 18.52 19.17 18.32 19.07 18.13 18.98 17.94 18.88 17.75 18.77 17.57 18.66 17.38 18.55 17.20 18.43 17.03 18.31 16.85 18.19 16.69 18.07 16.53 17.94 16.37 17.81 16.22 17.68 16.08 17.55 15.94 17.41 15.81 17.28 15.68 17.14 15.56 17.01 15.45 16.87 15.34 16.74 15.23 16.61 15.13 16.48 15.04 16.35 14.95 16.23 14.87 16.10 14.79 15.98 14.72 15.87 14.65 15.76 14.58 15.65 14.52 15.54 14.46 15.44 14.40 15.34 14.35 15.25 14.30 15.16 14.25 15.07 14.21 14.99 14.17 14.91 14.13 14.83 14.09 14.76
Layer 150 18.00 18.00 19.47 19.49 19.47 19.42 19.37 19.33 19.27 19.22 19.16 19.11 19.05 18.98 18.92 18.85 18.78 18.70 18.63 18.55 18.47 18.38 18.29 18.20 18.11 18.01 17.91 17.81 17.71 17.60 17.50 17.39 17.27 17.16 17.05 16.93 16.82 16.70 16.59 16.47 16.36 16.25 16.14 16.03 15.92 15.82
Layer 200 18.00 18.00 19.46 19.50 19.49 19.48 19.45 19.42 19.39 19.36 19.32 19.29 19.25 19.22 19.18 19.14 19.10 19.05 19.01 18.96 18.91 18.86 18.81 18.76 18.70 18.64 18.58 18.52 18.45 18.39 18.32 18.25 18.17 18.09 18.01 17.93 17.85 17.76 17.67 17.58 17.49 17.39 17.30 17.20 17.10 17.00
Layer 250 18.00 18.00 18.03 19.52 19.54 19.53 19.54 19.53 19.50 19.48 19.46 19.44 19.41 19.39 19.36 19.34 19.31 19.28 19.25 19.22 19.19 19.16 19.13 19.09 19.06 19.02 18.98 18.94 18.90 18.86 18.81 18.77 18.72 18.67 18.62 18.57 18.52 18.46 18.40 18.34 18.28 18.21 18.15 18.08 18.01 17.94
Ratarata 18.00 18.68 18.92 18.98 18.83 18.69 18.55 18.42 18.28 18.16 18.03 17.91 17.79 17.68 17.57 17.46 17.35 17.25 17.15 17.05 16.96 16.86 16.77 16.69 16.60 16.51 16.43 16.35 16.27 16.20 16.12 16.05 15.97 15.90 15.83 15.76 15.70 15.63 15.57 15.50 15.44 15.38 15.32 15.26 15.20 15.14
115
Lampiran 21 (lanjutan) 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99
13.72 13.71 13.69 13.68 13.67 13.66 13.65 13.64 13.63 13.62 13.62 13.61 13.60 13.59 13.59 13.58 13.58 13.57 13.57 13.56 13.56 13.55 13.55 13.54 13.54 13.54 13.53 13.53 13.53 13.52 13.52 13.52 13.52 13.51 13.51 13.51 13.51 13.51 13.50 13.50 13.50 13.50 13.50 13.50 13.50 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49
13.79 13.77 13.76 13.74 13.73 13.71 13.70 13.69 13.68 13.66 13.65 13.64 13.64 13.63 13.62 13.61 13.60 13.60 13.59 13.58 13.58 13.57 13.57 13.56 13.56 13.55 13.55 13.55 13.54 13.54 13.53 13.53 13.53 13.53 13.52 13.52 13.52 13.52 13.51 13.51 13.51 13.51 13.51 13.50 13.50 13.50 13.50 13.50 13.50 13.49 13.49 13.49 13.49 13.49
13.98 13.95 13.92 13.90 13.87 13.85 13.83 13.81 13.79 13.78 13.76 13.74 13.73 13.72 13.70 13.69 13.68 13.67 13.66 13.65 13.64 13.63 13.62 13.61 13.61 13.60 13.59 13.59 13.58 13.57 13.57 13.56 13.56 13.55 13.55 13.55 13.54 13.54 13.54 13.53 13.53 13.53 13.52 13.52 13.52 13.52 13.51 13.51 13.51 13.51 13.51 13.50 13.50 13.50
14.06 14.03 14.00 13.97 13.94 13.91 13.89 13.87 13.84 13.82 13.81 13.79 13.77 13.75 13.74 13.72 13.71 13.70 13.69 13.67 13.66 13.65 13.64 13.64 13.63 13.62 13.61 13.60 13.60 13.59 13.58 13.58 13.57 13.57 13.56 13.56 13.55 13.55 13.55 13.54 13.54 13.53 13.53 13.53 13.53 13.52 13.52 13.52 13.52 13.51 13.51 13.51 13.51 13.51
14.69 14.63 14.56 14.51 14.45 14.40 14.35 14.30 14.25 14.21 14.17 14.13 14.10 14.06 14.03 14.00 13.97 13.95 13.92 13.90 13.87 13.85 13.83 13.81 13.79 13.78 13.76 13.75 13.73 13.72 13.70 13.69 13.68 13.67 13.66 13.65 13.64 13.63 13.62 13.62 13.61 13.60 13.59 13.59 13.58 13.58 13.57 13.57 13.56 13.56 13.55 13.55 13.54 13.54
15.71 15.61 15.52 15.42 15.33 15.24 15.16 15.07 14.99 14.92 14.85 14.78 14.71 14.65 14.58 14.53 14.47 14.42 14.37 14.32 14.28 14.23 14.19 14.15 14.12 14.08 14.05 14.02 13.99 13.96 13.94 13.91 13.89 13.87 13.85 13.83 13.81 13.79 13.77 13.76 13.74 13.73 13.71 13.70 13.69 13.68 13.67 13.66 13.65 13.64 13.63 13.62 13.61 13.61
16.89 16.79 16.69 16.58 16.48 16.37 16.27 16.16 16.06 15.96 15.86 15.76 15.67 15.57 15.48 15.39 15.30 15.22 15.14 15.06 14.98 14.91 14.84 14.77 14.71 14.64 14.59 14.53 14.47 14.42 14.37 14.33 14.28 14.24 14.20 14.16 14.12 14.09 14.06 14.03 14.00 13.97 13.94 13.92 13.89 13.87 13.85 13.83 13.81 13.79 13.78 13.76 13.73 13.73
17.86 17.78 17.70 17.62 17.54 17.45 17.36 17.27 17.18 17.09 16.99 16.90 16.80 16.70 16.60 16.50 16.40 16.30 16.21 16.11 16.01 15.91 15.82 15.72 15.63 15.54 15.45 15.36 15.28 15.20 15.12 15.04 14.97 14.90 14.83 14.76 14.70 14.64 14.58 14.52 14.47 14.42 14.37 14.33 14.28 14.24 14.20 14.16 14.13 14.09 14.06 14.03 13.97 13.97
15.09 15.03 14.98 14.93 14.88 14.82 14.78 14.73 14.68 14.63 14.59 14.54 14.50 14.46 14.42 14.38 14.34 14.30 14.27 14.23 14.20 14.16 14.13 14.10 14.07 14.04 14.02 13.99 13.97 13.94 13.92 13.90 13.87 13.85 13.84 13.82 13.80 13.78 13.77 13.75 13.74 13.72 13.71 13.70 13.69 13.67 13.66 13.65 13.64 13.64 13.63 13.62 13.60 13.60
116
Lampiran 21 (lanjutan) 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120
13.49 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48
13.49 13.49 13.49 13.49 13.49 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48 13.48
13.50 13.50 13.50 13.50 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.48 13.48 13.48 13.48 13.48 13.48 13.48
13.50 13.50 13.50 13.50 13.50 13.50 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.49 13.48 13.48 13.48 13.48 13.48
13.54 13.53 13.53 13.53 13.52 13.52 13.52 13.52 13.51 13.51 13.51 13.51 13.51 13.50 13.50 13.50 13.50 13.50 13.50 13.50 13.49
13.60 13.59 13.59 13.58 13.58 13.57 13.56 13.56 13.56 13.55 13.55 13.54 13.54 13.54 13.53 13.53 13.53 13.52 13.52 13.52 13.52
13.72 13.71 13.69 13.68 13.67 13.66 13.65 13.64 13.63 13.62 13.62 13.61 13.60 13.60 13.59 13.58 13.58 13.57 13.57 13.56 13.56
13.94 13.92 13.90 13.87 13.85 13.83 13.81 13.80 13.78 13.76 13.75 13.73 13.72 13.71 13.69 13.68 13.67 13.66 13.65 13.64 13.63
13.60 13.59 13.58 13.58 13.57 13.57 13.56 13.56 13.55 13.55 13.55 13.54 13.54 13.53 13.53 13.53 13.52 13.52 13.52 13.52 13.51
117
Lampiran 22 Hasil pengukuran nilai kadar air bijian selama 50 jam pada Percobaan 1 Jam ke-
0
2
4
6
8
10
12
14
16
18
Point A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6
% b.k 18.24
16.02 19.51
Layer 10 Rata-rata % b.b 15.43 17.93
17.67
17.99 19.91
18.52
18.77
18.47
18.56
18.32
18.74
17.80
16.43 14.42
15.64
15.47 15.34
15.49
15.72 16.25
15.78
15.01 16.35
15.10
16.09 17.60
16.58
18.18 17.13
15.59
13.95
19.17
16.87 16.76
16.18 15.56
15.37
16.21
17.65 19.55
15.80
15.64
18.16
18.65 19.41
15.85 16.17
16.08
18.54
18.30 18.12
15.62
15.02
19.16
19.66 16.85
15.25 16.60
14.43 14.36
14.96
14.23 17.30
15.38 14.62
% b.k 18.30
17.80 18.50
Layer 40 Rata-rata % b.b 15.47 18.20
18.70
15.38
17.68
19.30 18.42
15.19
15.02
18.18
18.83 19.29
13.81 16.33
Rata-rata
14.74
18.90 19.00 17.82 17.40 17.49 16.76 17.57 17.09 17.82 17.50 17.97 17.95 17.34 17.03 17.95 17.75 18.10 18.11 17.34 17.07 17.88 17.91 18.23 18.06 17.94 17.57 18.74 18.12 19.10 18.55 18.64 18.12 18.24 18.10 18.65 18.21 18.14 18.45 18.48 19.47 18.33 18.26 18.12 18.38 18.06 18.30 18.64 18.57 18.08
15.11 15.61
Rata-rata
15.40
15.75 18.87
17.35
17.60
17.72
17.93
18.55
18.27
18.52
18.34
15.90 15.97 15.12 14.82 14.88 14.35 14.94 14.60 15.12 14.90 15.23 15.22 14.78 14.55 15.22 15.08 15.33 15.33 14.78 14.58 15.16 15.19 15.42 15.30 15.21 14.95 15.78 15.34 16.04 15.65 15.71 15.34 15.43 15.32 15.72 15.41 15.36 15.57 15.60 16.30 15.49 15.44 15.34 15.53 15.30 15.47 15.71 15.66 15.31
15.87
14.79
14.97
15.05
15.20
15.64
15.45
15.62
15.50
118
Lampiran 22 (lanjutan)
20
22
24
26
28
30
32
34
36
38
A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6
17.48
18.11 17.02
14.88 17.54
16.79
16.24 19.25
17.42
17.07
17.53
17.17
17.74
17.01
16.89
14.62 15.20
14.65
14.34 15.75
15.06
14.58 14.04
14.53
14.53 14.48
14.45
14.53 17.61
16.48
16.19 16.84
14.92
14.34
17.01
16.60 19.23
14.71 14.98
14.97
16.73
17.00 16.93
14.58
15.10
17.61
17.07 16.33
14.45 14.23
14.13
17.78
16.74 18.70
14.83
15.06
16.45
17.12 17.92
13.97 16.14 15.06
17.74
17.24 17.62
14.92
14.37
17.74
16.89 16.60
15.34 14.54
14.24 16.13
14.97
14.15 16.50
13.93 14.41
14.16
18.02 17.95 18.59 18.17 18.37 17.94 18.21 18.22 18.30 18.40 18.12 17.80 17.47 18.19 18.24 18.23 17.79 17.64 17.42 17.96 15.07 18.37 18.14 18.36 17.85 18.48 18.55 18.30 18.08 17.69 18.45 18.27 18.39 18.58 18.13 17.91 17.91 18.51 18.08 18.50 18.67 17.63 18.58 17.79 17.68 17.20 17.35 17.14 17.71 17.95 18.04 17.73 18.11 17.09 17.37 17.97 16.80 18.22 18.27 17.83
18.17
18.18
17.93
17.55
18.16
18.29
18.22
17.62
17.77
17.75
15.27 15.22 15.68 15.38 15.52 15.21 15.41 15.41 15.47 15.54 15.34 15.11 14.87 15.39 15.43 15.42 15.10 14.99 14.84 15.22 13.10 15.52 15.35 15.51 15.14 15.60 15.65 15.47 15.31 15.03 15.58 15.45 15.53 15.67 15.35 15.19 15.19 15.62 15.31 15.61 15.73 14.99 15.67 15.10 15.02 14.68 14.78 14.63 15.04 15.22 15.28 15.06 15.34 14.59 14.80 15.23 14.39 15.41 15.45 15.14
15.38
15.38
15.20
14.93
15.37
15.46
15.41
14.98
15.09
15.07
119
Lampiran 22 (lanjutan)
40
42
44
46
48
50
A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6
16.71
16.27 16.07
14.32 16.35
17.43
16.06 16.14
16.54
16.80
16.70
14.22 14.45
14.39
14.22 14.20
14.31
14.68 16.98
16.83
16.78 16.79
14.19
14.49
17.20
17.01 16.74
13.84 13.90 14.50
16.95
16.58 16.55
14.05
14.84
16.95
16.57 16.88
13.99 13.85
14.54 14.34
14.52
14.40 16.80
14.37 14.38
14.38
17.19 17.55 17.73 17.87 17.49 16.40 17.70 17.77 18.55 17.66 17.43 16.89 16.28 18.05 17.83 18.48 17.77 16.78 18.02 18.09 18.60 18.65 17.55 18.03 16.27 17.87 17.71 17.82 17.56 16.82 17.32 18.43 18.35 18.01 17.50 16.95
17.37
17.67
17.53
18.16
17.34
17.76
14.67 14.93 15.06 15.16 14.89 14.09 15.04 15.09 15.65 15.01 14.84 14.45 14.00 15.29 15.14 15.60 15.09 14.37 15.27 15.32 15.68 15.72 14.93 15.27 14.00 15.16 15.05 15.12 14.94 14.40 14.76 15.56 15.51 15.26 14.90 14.50
14.80
15.01
14.91
15.36
14.78
15.08
120
Lampiran 23 Pengaruh fluktuasi RH terhadap perubahan kadar air pengukuran pada Percobaan 1 Jam ke 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46 48 50 Rata - rata
Layer 10-Ukur % b.k % b.b 17.93 15.19 18.52 15.62 18.77 15.80 18.47 15.59 18.56 15.64 18.32 15.49 18.74 15.78 17.80 15.10 17.60 14.96 17.30 14.74 17.54 14.92 17.42 14.83 17.07 14.58 17.53 14.92 17.17 14.65 17.74 15.06 17.01 14.53 16.89 14.45 17.61 14.97 16.50 14.16 16.35 14.05 16.54 14.19 16.80 14.39 16.70 14.31 16.98 14.52 16.80 14.38
Layer 40-Ukur % b.k % b.b 18.20 15.40 18.87 15.87 17.35 14.79 17.60 14.97 17.72 15.05 17.93 15.20 18.55 15.64 18.27 15.45 18.52 15.62 18.34 15.50 18.17 15.38 18.18 15.38 17.93 15.20 17.55 14.93 18.16 15.37 18.29 15.46 18.22 15.41 17.62 14.98 17.77 15.09 17.75 15.07 17.37 14.80 17.67 15.01 17.53 14.91 18.16 15.36 17.34 14.78 17.76 15.08
RH lingk. (%)
T in (oC)
95.0 78.8 65.9 68.2 67.6 88.0 79.4 60.5 64.0 75.0 95.0 67.3 60.0 59.0 60.0 91.0 75.4 68.3 60.0 75.8 93.0 79.5 69.1 59.0 70.0 73
23.5 31.0 33.0 33.5 33.5 25.0 31.5 32.0 32.0 33.4 24.5 31.2 32.2 32.3 32.3 24.3 32.0 32.3 32.4 32.5 24.5 31.7 32.0 32.5 33.0 31
121
Lampiran 24 Hasil pengukuran nilai kadar air bijian selama 40 jam pada Percobaan 2 Jam ke-
0
2
4
6
8
10
12
14
16
18
Point A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6
% b.k 20.48
18.60 18.18 17.86 19.54 18.39 18.60 19.10 19.32 18.82 18.82 17.65 19.05 18.82 18.60 17.24 16.28 19.05 18.18 18.39 18.60 18.60 17.24 18.18 19.05 18.39 18.39 17.58 19.05 19.05 18.18 18.39 18.60 19.05 19.05 19.05 18.82 19.05 17.24 17.65 20.00 19.05 20.24 17.65 17.65 18.82 18.60 19.05 17.24 19.05 18.82 17.24
16.09 17.78
Layer 10 Rata-rata % b.b 17.00 19.09
18.80
18.63
17.96
18.31
18.45
18.71
18.70
18.60
17.04
15.69 15.38 15.15 16.35 15.53 15.69 16.04 16.19 15.84 15.84 15 16.00 15.84 15.69 14.71 14 16 15.38 15.53 15.69 15.69 14.71 15.38 16.00 15.53 15.53 14.95 16.00 16.00 15.38 15.53 15.69 16.00 16.00 16.00 15.84 16.00 14.71 15.00 16.67 16.00 16.83 15.00 15.00 15.84 15.69 16.00 14.71 16.00 15.84 14.71
13.86 15.09
Rata-rata
16.02
15.82
15.70
15.22
15.47
15.59
15.76
15.75
15.68
14.55
% b.k 21.18
20.24 19.05 18.39 17.24 20.24 18.39 18.60 18.18 19.32 18.60 16.09 20.59 18.18 17.44 18.60 18.60 17.98 20.00 19.05 17.24 18.39 19.10 18.75 19.54 18.82 18.39 18.82 17.44 20.24 16.09 18.82 17.65 17.65 18.60 19.05 18.60 17.65 17.24 16.28 18.82 18.82 15.91 17.44 18.82 17.44 19.05 17.65 18.60 18.82 17.65 15.05
16.28 16.67
Layer 40 Rata-rata % b.b 17.48 20.15
18.51
18.37
18.58
18.83
18.18
18.13
17.68
18.20
16.00
16.83 16.00 15.53 14.71 16.83 15.53 15.69 15.38 16.19 15.69 13.86 17.07 15.38 14.85 15.69 15.69 15.24 16.67 16.00 14.71 15.53 16.04 15.79 16.35 15.84 15.53 15.84 14.85 16.83 13.86 15.84 15.00 15.00 15.69 16.00 15.69 15.00 14.71 14.00 15.84 15.84 13.73 14.85 15.84 14.85 16.00 15.00 15.69 15.84 15.00 13.08
14.00 14.29
Rata-rata
16.77
15.61
15.51
15.66
15.85
15.37
15.35
15.02
15.40
13.79
122
Lampiran 24 (lanjutan)
20
22
24
26
28
30
32
34
36
38
40
A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6 A1 A2 A3 A4 A5 A6
15.38
14.89
13.33 15.71
16.85 14.77
17.24
15.65
15.93
13.87
14.52
14.65
13.77
14.34
14.46
13.48
12.67
12.75
11.76
12.26
11.93
13.88
11.43
12.78
12.87 11.88
16.16
15.56
16.48
16.09
12.10
17.24
17.23
16.30
16.23
15.73
16.72
16.48
16.99
18.18 17.24
13.46
14.67
14.15
13.97
13.86
14.32
14.71
14.52
15.38 14.71 16.93
14.02
14.48
14.71 12.62 15.13
13.59
13.14
13.21 14.29 16.43
16.13 15.56 11.88
13.91
15.24 13.46
15.22 16.67 12.18
13.59
13.86 13.86
17.24 14.44 12.63
14.24
13.89 13.89
18.18 17.24 12.54
14.15
14.29 16.67
17.98 15.56
12.38 10.89 13.49
15.73
14.10
13.46 13.86
16.09 16.09
13.21 12.75
14.13 12.22
14.77
11.65
16.61
16.13 16.13
12.75 12.75
15.22 14.61
12.90
12.08
12.50 12.62
14.61 14.61
13.54
11.54
16.48
13.73 14.29 15.09
16.67 20.00
13.21 12.04
14.29 14.44
13.98
13.73
13.86 12.38
15.22 13.68
13.33
14.56
16.41
15.56 16.09
12.50
16.09 14.13
14.61
13.53
12.62 12.62
14.13 14.29
13.19
14.71
15.91
14.29
16.67 17.78
13.00 14.02
14.44 14.44
13.04
13.57
14.42 12.87
14.94 16.30
17.05
12.96
16.67
14.15
14.11
13.89 13.46 16.99
15.38 14.71
14.52
123
Lampiran 25 Pengaruh fluktuasi RH terhadap perubahan kadar air pengukuran pada Percobaan 2
0
Layer 10-Ukur % b.k % b.b 19.09 16.02
Layer 40-Ukur % b.k % b.b 20.15 16.77
2
18.80
15.82
18.51
4
18.63
15.70
18.37
6
17.96
15.22
8
18.31
15.47
10
18.45
12
18.71
14 16
Jam ke
RH lingk. (%)
T in (oC)
-
-
15.61
68.4
29.7
15.51
58.4
32.8
18.58
15.66
58.4
34.2
18.83
15.85
58.4
34.5
15.59
18.18
15.37
58.1
33.5
15.76
18.13
15.35
72.7
28.2
18.70
15.75
17.68
15.02
63.9
33.4
18.60
15.68
18.20
15.40
55.9
34.3
18
17.04
14.55
16.00
13.79
61.0
33.5
20
15.71
13.57
16.41
14.10
55.5
33.7
22
15.65
13.53
16.61
14.24
58.4
33.6
24
15.93
13.73
16.16
13.91
60.3
33.4
26
13.87
12.08
17.23
14.67
56.8
33.6
28
14.52
12.67
16.23
13.97
58.4
34.0
30
14.65
12.78
16.72
14.32
55.9
33.7
32
13.77
12.10
16.99
14.52
62.5
29.5
34
14.34
12.54
16.93
14.48
58.1
33.5
36
14.46
12.63
15.13
13.14
58.4
34.1
38
13.88
12.18
16.43
14.11
58.4
34.0
40
13.49
11.88
16.99
14.52
57.2
33.7
59.8
33
Rata-rata
124
Lampiran 26 Validasi perubahan kadar air pengukuran dan simulasi pada Percobaan 1. Kadar air (% b.k)
Jam ke-
Layer 10 Ukur
Simulasi
Layer 40 SD
Error
Ukur
Simulasi
SD
Error
0
17.93
18.00
0.05
0.07
18.20
18.00
0.14
0.20
2
18.52
19.30
0.55
0.78
18.87
19.40
0.38
0.54
4
18.77
19.10
0.23
0.33
17.35
19.24
1.33
1.89
6
18.47
18.91
0.31
0.44
17.60
19.09
1.05
1.48
8
18.56
18.74
0.13
0.18
17.72
18.94
0.86
1.22
10
18.32
18.58
0.18
0.26
17.93
18.80
0.61
0.86
12
18.74
18.44
0.21
0.30
18.55
18.66
0.08
0.12
14
17.80
18.31
0.35
0.50
18.27
18.54
0.19
0.27
16
17.60
18.18
0.41
0.58
18.52
18.42
0.07
0.10
18
17.30
18.07
0.55
0.77
18.34
18.30
0.03
0.04
20
17.54
17.97
0.30
0.43
18.17
18.20
0.02
0.02
22
17.42
17.87
0.31
0.45
18.18
18.10
0.06
0.08
24
17.07
17.78
0.50
0.71
17.93
18.00
0.05
0.07
26
17.53
17.70
0.12
0.17
17.55
17.91
0.25
0.36
28
17.17
17.62
0.32
0.46
18.16
17.83
0.23
0.33
30
17.74
17.55
0.13
0.19
18.29
17.75
0.38
0.54
32
17.01
17.49
0.34
0.48
18.22
17.68
0.38
0.54
34
16.89
17.43
0.38
0.54
17.62
17.61
0.01
0.01
36
17.61
17.37
0.17
0.24
17.77
17.55
0.16
0.22
38
16.50
17.32
0.58
0.82
17.75
17.49
0.18
0.26
40
16.35
17.27
0.65
0.92
17.37
17.43
0.04
0.06
42
16.54
17.23
0.48
0.68
17.67
17.38
0.20
0.29
44
16.80
17.19
0.27
0.38
17.53
17.33
0.14
0.20
46
16.70
17.15
0.32
0.45
18.16
17.28
0.62
0.87
48
16.98
17.11
0.09
0.13
17.34
17.24
0.07
0.10
50
16.80
17.08
0.20
0.28
17.76
17.20
0.40
0.56
max
18.77
19.30
0.65
0.92
18.87
19.40
1.33
1.89
min
16.35
17.08
0.05
0.07
17.34
17.20
0.01
0.01
Rata-rata Total Error Korelasi ( R )
0.31
0.44
0.31
0.43
11.57
11.23
0.90
0.35
125
Lampiran 27 Validasi perubahan kadar air pengukuran dan simulasi pada Percobaan 2 Kadar air (% b.k)
Jam ke-
Layer 10 Ukur
Simulasi
Layer 40 SD
Error
Ukur
Simulasi
SD
Error
0
19.09
18.00
0.77
1.09
20.15
18.00
1.52
2.15
2
18.80
18.74
0.05
0.06
18.51
19.21
0.49
0.70
4
18.63
18.04
0.41
0.58
18.37
18.76
0.28
0.39
6
17.96
17.44
0.36
0.51
18.58
18.32
0.18
0.26
8
18.31
16.93
0.98
1.38
18.83
17.89
0.67
0.94
10
18.45
16.48
1.39
1.97
18.18
17.48
0.50
0.70
12
18.71
16.10
1.85
2.61
18.13
17.08
0.74
1.05
14
18.70
15.76
2.08
2.94
17.68
16.72
0.68
0.96
16
18.60
15.48
2.21
3.12
18.20
16.38
1.29
1.82
18
17.04
15.23
1.28
1.81
16.00
16.07
0.05
0.07
20
15.71
15.01
0.49
0.70
16.41
15.79
0.44
0.62
22
15.65
14.83
0.58
0.83
16.61
15.54
0.76
1.07
24
15.93
14.66
0.90
1.27
16.16
15.31
0.61
0.86
26
13.87
14.52
0.46
0.65
17.23
15.10
1.50
2.13
28
14.52
14.40
0.09
0.12
16.23
14.92
0.93
1.31
30
14.65
14.29
0.25
0.36
16.72
14.76
1.39
1.96
32
13.77
14.20
0.30
0.43
16.99
14.62
1.68
2.38
34
14.34
14.11
0.16
0.23
16.93
14.49
1.73
2.44
36
14.46
14.04
0.29
0.41
15.13
14.38
0.53
0.76
38
13.88
13.98
0.07
0.10
16.43
14.28
1.52
2.15
40
13.49
13.92
0.30
0.43
16.99
14.19
1.98
2.81
max
19.09
18.74
2.21
3.12
20.15
19.21
1.98
2.81
min
13.49
13.92
0.05
0.06
15.13
14.19
0.05
0.07
0.32
0.45
0.30
0.43
Rata-rata Total Error
9.49
9.06
Korelasi ( R )
0.88
0.84
126
Lampiran 28 Mutu jagung pada percobaan a. Perbandingan parameter mutu hasil percobaan dengan persyaratan mutu SNI 01-03920-1995 Persyaratan Mutu SNI (% Maks) No
Komponen Utama
1 2 3 4 5
Kadar air Butir rusak Butir warna lain Butir pecah Kotoran
I
II
III
IV
14 2 1 1 1
14 4 3 4 1
15 6 7 3 2
17 8 10 5 2
Mutu Jagung Percobaan Sebelum Setelah Setelah Pengeringan Pengeringan Penyimpanan 16.40 13.20 12.33 7.03 7.03 7.04 0.43 0.43 0.51 4.10 4.10 4.10 1.07 1.07 1.11
b. Perbandingan parameter mutu hasil percobaan dengan persyaratan mutu SNI 01-6944-2003 No
Komponen
Persyaratan SNI
1 2 3 5
Kadar air (% maks) Benih murni (% min) Daya berkecambah (% min) Kotoran Benih (% maks)
12 98 85 2
Sebelum proses 16.40 85.8 85.0 14.50
Jagung Percobaan Setelah Pengeringan 13.20 85.8 84.7 14.50
Setelah proses 12.33 85.5 84.7 14.50
c. Perbandingan parameter mutu hasil percobaan dengan persyaratan mutu SNI 01-6944-2003 No 1 2 3 4 5 6 7 8 9 10
Komponen Kadar air (% maks) Kadar protein kasar (% min) Kadar serat kasar (% maks) Kadar abu (% maks) Kadar lemak (% min) Mikotoksin: Aflatoksin (ppb maks) Okratoksin (ppb maks) Butir pecah (% maks) Warna lain (% maks) Benda asing (% maks) Kepadatan (min) kg/m3
Persyaratan SNI 14 7.5 3 2 3
Sebelum Proses 16.40 7.76 2.56 1.24 3.80
Setelah Proses 12.33 6.62 2.53 1.23 3.24
50 5 5 5 2 700.0
18.48 4.10 0.43 0 700.0
21.10 4.10 0.51 0 700.0
127
Lampiran 29 Kontaminasi aflatoxin pada jagung hasil uji laboratorium sebelum proses pengeringan dan penyimpanan
128
Lampiran 30 Kandungan abu, lemak dan serat pada jagung hasil uji laboratorium sebelum proses pengeringan dan penyimpanan.
129
Lampiran 31 Kandungan protein pada jagung hasil uji laboratorium sebelum proses pengeringan dan penyimpanan
130
Lampiran 32 Kontaminasi aflatoxin pada jagung hasil uji laboratorium setelah proses pengeringan dan penyimpanan
131
Lampiran 33 Kandungan abu, lemak dan serat pada jagung hasil uji laboratorium setelah proses pengeringan dan penyimpanan
132
Lampiran 34 Kandungan protein pada jagung hasil uji laboratorium setelah proses pengeringan dan penyimpanan
133
Lampiran 35 L 3 Sistem Peengering Efeek Rumah Kaca K Hybrid (ERK-Hybrrid) dan InStore Dryyer Terintegrrasi yang di ujicoba dalaam penelitiann
a Pengerinng ERK-Hybbrid dan In-SStore Dryer tterintegrasi a)
a) In-Store D Dryer dan Bucket Elevatoor
134
Lampiran 36 Alat-alat ukur L u yang diigunakan dallam penelitiaan
a) Term mokopel untuuk mengukurr suhu dalam m tump pukan biji jaggung
i
ii
b) Termokkopel untuk mengukur bola basahh; i) untuk menguk kur bola bassah dalam tuumpukan bijji, ii) botol tempat air a sebagai bbola basah
c. Moissture tester digital d
135
Lampiran 36 6 (lanjutan)
a) H Hybrid Recorrder
b) Chino recoorder
c)) Anemometer
136
Lampiran 36 6 (lanjutan)
a) Tiimbangan diigital
b) Jangka sorong
ng c) Oven dryin