SIMULATION OF PARTICLE SIZE EFFECT ON LOW RANK COAL GASIFICATION Name Department Advisors
: Nuur Faridatun H (2305 100 004) Ida Meipurwati (2305 100 016) : Chemical Engineering FTI-ITS : Prof. Dr. Ir Sugeng Winardi, M.Eng Dr. Tantular Nurtono, ST., M Eng
Abstract This research is conducted by simulation method based on CFD (Computational Fluid Dinamycs) because the experimental method is very difficult to be conducted and considered to be inneffctive due to the complexity of coal gasification process. The difficulties are multiphase reaction involved in the process, very high operating temperature, instantaneous reaction, and the difficulty to obtain estimation of gasification mechanism in the gasifier which depends on the equipment design and its fluid flow pattern. Gasification is the process of converting various carbon-based feedstocks to clean synthetic gas (syngas), which is primarily a mixture of hydrogen (H2) and carbonmonoxide (CO) as fuels, through an incomplete combustion. This gasification process is very interesting because it produces energy and syngas that can be used for industry necessity. Numerical simulations of the oxygen-blown coal gasification process inside an entrained flow gasifier involve chemical reaction, mass and heat transfer, and simultaneously fluid flow pattern can display detail information about what is occurred in the gasifier. The detail display of temperature distribution profile, reaction efficiency, and gas composition as the result of the gasification process could provide alternative solution to optimize the operating conditions of the low rank coal gasification. To increase the performance of gasification, the
3
study is required to analyze the conditions of low rank coal especially about granule size or coal particle size. Cylinder geometry is used for gasifier with 3 m of length and the diameter is 20 cm. The material used in this research is pulverized coal and pure oxygen as an oxidizer. Pulverized coal and oxygen are injected into gasifier as the primary inlet at 1200 K. Oxygen mass flow is 0,0066 kg/s. Variabel which applied is various of particle size. The methodology that will be conducted includes some stages, the first step is to document physical and chemical properties datas of low rank coal into the CFD shoftware, display the gasification equipment by the preprocessing shoftware, simulate the gasification process by variating operating condition and various of particle size. The bigger the particle size, the higher the temperature coming out of the gasifie . Mass fraction of CO decrease with the bigger the particle size. But mass fraction of H2 increase with the bigger the particle size. Key words : Gasification, coal, syngas, particle size of coal
4
SIMULASI PENGARUH UKURAN PARTIKEL PADA GASIFIKASI BATUBARA KUALITAS RENDAH TUGAS AKHIR Diajukan Untuk Memenuhi Salah Satu Syarat Memperoleh Gelar Sarjana Teknik pada Program Studi S-1 Jurusan Teknik Kimia Fakultas Teknologi Industri Institut Teknologi Sepuluh Nopember Oleh: NUUR FARIDATUN H NRP 2305 100 004
IDA MEIPURWATI NRP 2305 100 016
Disetujui oleh Tim Penguji Tugas Akhir: 1. Prof. Dr. Ir. Sugeng Winardi, M.Eng
…....................
(Pembimbing) 2. Dr. Tantular Nurtono, ST. M.Eng
........................
(Pembimbing) 3. Dr. Ir. Heru Setyawan, M.Eng
........................
(Penguji) 4. Ir. Minta Yuwana, MS.
........................
(Penguji) 5. Dr. Ir. Arief Widjaja, M.Eng (Penguji) SURABAYA AGUSTUS, 2009
5
........................
KATA PENGANTAR Puji syukur ke hadirat Allah SWT atas rahmat dan ridhoNya sehingga kami dapat menyelesaikan skripsi yang berjudul :
“SIMULASI PENGARUH UKURAN PARTIKEL PADA GASIFIKASI BATUBARA KUALITAS RENDAH” Pada kesempatan ini penulis mengucapkan terima kasih atas segala sesuatu yang diberikan sehingga skripsi ini dapat kami selesaikan, yaitu kepada : 1. Bapak Prof. Dr. Ir. Sugeng Winardi, M.Eng selaku Dosen Pembimbing dan Kepala Laboratorium Mekanika Fluida dan Pencampuran atas bimbingan, doa, dukungan dan fasilitas yang telah diberikan sehingga kami dapat menyelesaikan skripsi ini. 2. Bapak Dr. Tantular Nurtono, ST. M.Eng selaku Dosen Pembimbing di Laboratorium Mekanika Fluida dan Pencampuran atas bimbingan, masukan, dan sarannya yang amat bermanfaat sehingga skripsi ini bisa terselesaikan. 3. Ibu Dr.Widiyastuti,ST. MT selaku dosen di Laboratorium Mekanika Fluida dan Pencampuran atas doa dan dukungannya kepada kami. 4. Dr. Ir. Tri Widjaja, M.Eng, selaku ketua Jurusan Teknik Kimia FTI – ITS Surabaya atas fasilitas yang diberikan. 5. Dr. Ir. Kusno Budikarjono, MT selaku koordinator Skripsi dan Tugas Akhir Jurusan Teknik Kimia FTI – ITS. 6. Bapak Prof. Dr. Ir. Ali Altway, MS yang telah memberikan saran dan masukan yang amat berharga terhadap skripsi kami. 7. Bapak Dr. Ir. Heru Setyawan, M.Eng, Bapak Ir. Minta Yuwana, MS. dan Bapak Dr. Ir. Arief Widjaja, M.Eng selaku dosen penguji atas saran dan bimbingannya yang sangat bermanfaat demi kelancaran skripsi ini.
6
8. Ibu Ir. Pantjawarni P dan Bapak Ir. Samsudin Affandi, MS selaku dosen wali kami atas bimbingan yang diberikan selama kami menjadi mahasiswa. 9. Dosen – dosen pengajar serta seluruh karyawan Jurusan Teknik Kimia FTI – ITS. 10. Orang tua dan keluarga kami tercinta yang telah memberi dukungan moral, doa, semangat, nasehat, dan material yang tiada hentinya kepada kami. 11. Saudara-saudara kami Mixing Crew 2008 – 2009, Bang Farid dan teman-teman K-45 di Jurusan Teknik Kimia FTI-ITS atas kerja sama dan bantuannya selama proses penyelesaian skripsi ini. 12. Pihak-pihak yang tidak dapat kami sebutkan satu per satu yang telah membantu kelancaran proses penyelesaian skripsi kami. Kami menyadari bahwa masih terdapat kesalahan dalam pengerjaan skripsi ini. Oleh karena itu, saran dan kritik yang bersifat membangun sangat kami harapkan demi peningkatan kualitas skripsi ini. Semoga skripsi ini dapat bermanfaat dikemudian hari. Aamiin. Surabaya,Agustus 2009
Penyusun
7
DAFTAR ISI
HALAMAN JUDUL LEMBAR PENGESAHAN ABSTRAK i v KATAPENGANTAR DAFTAR ISI ............................................................................. vii DAFTAR GAMBAR.................................................................. ix DAFTAR TABEL ....................................................................... x DAFTAR NOTASI .................................................................... xi BAB I PENDAHULUAN I.1 Latar Belakang ....................................................................... 1 I.2 Perumusan Masalah ............................................................... 2 I.3 Batasan Masalah..................................................................... 2 I.4 Tujuan Penelitian ................................................................... 2 I.5 Manfaat Penelitian ................................................................. 3 BAB II TINJAUAN PUSTAKA II.1 Batubara................................................................................ 5 II.2 Proses Gasifikasi................................................................... 7 II.2.1 Tipe Gasifier.............................................................. 10 II.2.2 Keuntungan Teknologi Gasifikasi Pada Batubara..... 11 II.3 Peneliti Terdahulu............................................................... 12 II.4 Model Matematika Proses Gasifikasi ................................. 14 II.4.1 Transport Senyawa.................................................... 14 II.4.2 Turbulensi.................................................................. 16 II.4.3 Persamaan Energi...................................................... 17 II.4.4 Discrete Phase Model ............................................... 18 II.4.5 Particle Size Distribution .......................................... 24 II.5 Computational Fluid Dynamics (CFD) .............................. 25 BAB III METODOLOGI PENELITIAN III.1 Sistem yang dipelajari ....................................................... 29 III.2 Bahan yang Digunakan...................................................... 30 III.3 Simulasi ............................................................................. 32 III.3.1 Model Sistem yang Dipelajari................................ 32
8
III.3.2 Parameter Devolatilisasi......................................... 33 III.3.3 Parameter Reaksi Heterogen (Char-Gas) .............. 33 III.3.4 Parameter Reaksi Homogen ................................... 34 III.3.5 Parameter Ukuran Partikel Batubara...................... 35 III.4 Kondisi Batas................................................................... 36 III.5 Prosedur Simulasi ............................................................ 36 III.6 Variabel Penelitian........................................................... 36 BAB IV HASIL PENELITIAN DAN PEMBAHASAN IV.1 Pengaruh Ukuran Partikel Batubara Terhadap Distribusi Temperatur...................................................... 38 IV.2 Pengaruh Ukuran Partikel Batubara Terhadap Fraksi Massa CO ............................................................. 41 IV.3 Pengaruh Ukuran Partikel Batubara Terhadap Fraksi Massa H2 .............................................................. 45 IV.4 Pengaruh Ukuran Partikel Batubara Terhadap Fraksi Massa CH4............................................................ 49 IV.5 Pengaruh Ukuran Partikel Batubara Terhadap Fraksi Massa H2O ........................................................... 52 IV.6 Pengaruh Ukuran Partikel Batubara Terhadap Fraksi Massa CO2 ............................................................ 55 IV.7 Pengaruh Ukuran Partikel Batubara Terhadap Karakteristik dan Distribusi Medan Kecepatan Aliran.... 59 BAB V PENUTUP V.1 Kesimpulan...................................................................... 63 V.2 Saran ................................................................................ 63 DAFTAR PUSTAKA
9
DAFTAR GAMBAR Gambar II.1. Gambar II.2.
Kondisi batas untuk pemodelan reflect ........... 21 A Reacting Particle in the Multiple Surface Reactions Model.............................................. 23 Gambar III.1. Dimensi Peralatan Gasifier ............................. 30 Gambar III.2. Grid Gasifier.................................................... 32 Gambar IV.1. Distribusi Temperatur (K) Untuk Berbagai Variasi Ukuran Batubara ............................... 38 Gambar IV.2. Distribusi Fraksi Massa CO............................. 41 Gambar IV.3. Hubungan Antara Tinggi Gasifier [m] dengan Fraksi Massa CO ............................................. 42 Gambar IV.4. Distribusi Fraksi Massa H2 ............................. 45 Gambar IV.5. Hubungan Antara Tinggi Gasifier [m] dengan Fraksi Massa H2 .............................................. 46 Gambar IV.6. Distribusi Fraksi Massa CH4 ........................... 49 Gambar IV.7. Hubungan Antara Tinggi Gasifier [m] dengan Fraksi Massa CH4 .......................................... 50 Gambar IV.8. Distribusi Fraksi Massa H2O ........................... 52 Gambar IV.9. Hubungan Antara Tinggi Gasifier [m] dengan Fraksi Massa H2O........................................... 53 Gambar IV.10. Distribusi Fraksi Massa CO2 ........................... 55 Gambar IV.11. Hubungan Antara Tinggi Gasifier [m] dengan Fraksi Massa CO2 ........................................... 56 Gambar IV.12. Vektor medan kecepatan aliran (m/s) untuk variasi ukuran partikel batubara ...................... 59
10
DAFTAR TABEL Tabel II.1 Analisis Komposisi Batubara di Indonesia............. 6 Tabel II.2 Komposisi Volatile Matter...................................... 7 Tabel II.3 Standard Baku Mutu Udara Ambient ................... 12 Tabel III.1 Komposisi Batubara menurut proximate dan elemental analysis ......................................... 31 Tabel III.2 Data Kinetika untuk devolatilisasi kobayashi model................................................... 33 Tabel III.3 Data Kinetika untuk devolatilisasi single rate model ................................................ 33 Tabel III.4 Data Kinetika untuk reaksi heterogen .................. 34 Tabel III.5 Data Kinetika untuk reaksi homogeny ................. 35 Tabel III.6 Parameter Rosin-Rammler Size Distribution........ 35 Tabel IV.1 Suhu Gas Keluar Gasifier..................................... 40 Tabel IV.2 Komposisi CO (fraksi massa) Keluar Gasifier untuk Berbagai Ukuran Partikel.......................... 42 Tabel IV.3 Komposisi H2 (fraksi massa) Keluar Gasifier untuk Berbagai Ukuran Partikel........................... 46 Tabel IV.4 Komposisi CH4 (fraksi massa) Keluar Gasifier untuk Berbagai Ukuran Partikel ............ 50 Tabel IV.5 Komposisi H2O (fraksi massa) Keluar Gasifier untuk Berbagai Ukuran Partikel ............ 53 Tabel IV.6 Komposisi CO2 (fraksi massa) Keluar Gasifier untuk Berbagai Ukuran Partikel ............ 56
11
DAFTAR NOTASI
NOTASI cmf f g k P T D L mf ma x, y, z’ u, v, w’ Wi e ρ me tw vx v σ Ji Yi Di Sct
KETERANGAN Koefisien konstanta Mixture fraction Gravitasi Energi kinetik turbulen Tekanan Temperature Diameter gasifier Tinggi Gasifier Mass flow batubara Mass flow oksigen Komponen koordinat Komponen kecepatan Berat komponen i Energi Dissipasi Densitas Fluida Viskositas efektif Shear stress di dinding Kecepatan fluida arah x Viscositas kinematik Konstanta pemodelan Fluks difusi spesies i Fraksi massa local masing-masing spesies Diffusivitas Turbulent schmidt number
12
SATUAN [-] [-] [m/s2] [m2/s2] [bar] [K] [m] [m] [kg/s] [kg/s] [-] [-] [gr] [W/kg] [kg/m3] [Pa.s] [Pascal] [m/s] [m2/s] [-] [kgmol/s.m2] [-] [m2/s] [dimensionless]
µt Ri,r Mw,i Er Ar n µj,r Nr YR YP YM σk σε C2 C1ε Sh mj’ hoj Rj’ CD σ H T∞ Ap
Viskositas turbulen Arehenius molar rate untuk pembentukan/ hilangnya spesies i dalam reaksi r Berat molekul spesies i Activation energy untuk reaksi r Pre-exponensial factor Temperatur exponent Eksponent untuk spesies j dalam rekasi r Jumlah spesies dalam reaksi r Fraksi massa reaktan R Fraksi massan produk P Kontribusi dari dilatasi yang berfluktuasi dalam compressible Bilangan Prandtl turbulen untuk k Bilangan Prandtl turbulen untuk ε Konstanta Konstanta Panas yang berasal dari reaksi kimia dan sumber panas lainnya Fraksi massa dari spesies j' Entalpi pembentukan spesies j' Rate volumetrik dari pembentukan spesies j' Koefisien seret Konstanta Stefan-Boltzmann Koefisien transfer panas konveksi Temperatur lokal dari fase kontinyu Luas permukaan partikel
13
[Pa.s] [-] [kg/kgmol] [J/kg] [dimensionless] [-] [-] [-] [-] [-] [-] [dimensionless] [dimensionless] [-] [-] [....] [-] [J/kg] [m3/s] [dimensionless] [W/m2-K4] [W/m2.K] [K] [m2]
Cp mp G Dp,0 Nr D0,r Pr Rj,r µ
Kapasitas panas dari partikel Massa partikel Radiasi insiden Diameter partikel pada awal devolatiliassi Orde reaksi dari r Laju reaksi dari partikel surface spesies per unit area Effectiveness factor Laju reaksi dari particle surface per unit area Laju deplesi dari partikel surface spesies
14
[J/kg.K] [kg] [W/kg] [µm] [-] [kg/m2.s] [dimensionless] (kg/m2-s) [kg/s]