LAPORAN HASIL PENELITIAN SESUAI PRIORITAS NASIONAL

Perubahan iklim dan pelestrarian lingkungan (biodiversity)

LAPORAN HASIL PENELITIAN SESUAI PRIORITAS NASIONAL

PENINGKATAN REAKTIFITAS ABSORBENT Ca(OH)2 UNTUK PENYISIHAN SO2 HASIL PEMBAKARAN SAMPAH MENGGUNAKAN BAG FILTER REKTOR: PERCOBAAN DAN MODELING Oleh: Dr. Ir. Mariana, M.Si Dr. Ir. Syahrul, M.Sc. Dr. Ir. Farid Mulana, M.Eng

Dibiayai oleh Direktorat Jenderal Pendidikan Tinggi, Departemen Pendidikan Nasional, sesuai dengan Surat Perjanjian Pelaksanaan Hibah Kompetitif Penelitian Sesuai Prioritas Nasional Nomor : 212/SP2H/PP/DP2M/V/2009, Tanggal 30 Mei 2009

UNIVERSITAS SYIAH KUALA NOPEMBER, 2009

RINGKASAN SO2 merupakan jenis gas yang sangat berbahaya bagi kesehatan makhluk hidup. Kehadiran SO2 di dalam udara dapat berekasi dengan radikal bebas (seperti O • dan OH • ) membentuk asam dan dapat menyebabkan hujan asam. Hujan asam bila sampai dibumi akan menyebabkan berbagai macam penyakit seperti batuk, iritasi tenggorokan, gatal dan lain-lain yang merupakan penyakit tropis. Disamping hal tersebut diatas, reaksi antara SO2 dengan radikal bebas menyebabkan penipisan lapisan ozone, sehingga akan menyebabkan timbulnya panas global (global warming) bumi. Teknologi penghilangan SO2 dari hasil pembakaran sampah dapat dilakukan dengan menggunakan proses kering maupun proses basah.

Penghilangan SO2 dengan menggunakan proses kering lebih ekonomis bila

dibandingkan dengan proses basah karena sederhana, mudah pemeliharaan dan tidak menghasilkan limbah cair. Namun demikian, kelemahan proses kering adalah konversi absorsi rendah dan efisiensi penyisihan

SO2 relatif kecil. Salah satu cara mengatasi

masalah tersebut di atas adalah dengan menggunakan absorbent yang mempunyai reaktifitas yang tinggi serta pemahaman terhadap mekanisme penyerapan SO2 di dalam absorbent dengan pembentukan model matematik. Penelitian ini diharapkan dapat menghasilkan: (1) sorbent dengan reaktifitas yang tinggi dan (2) model matematik yang dapat menggambarkan fenomena absorbsi SO2 menggunakan absorbent. Penelitian ini ditargetkan dapat digunakan untuk menyerap gas hasil pembakaran secara efektif, mudah dan murah.

ii

SUMMARY SO2 is a type of gas which is very dangerous for health of living beings. The presence of SO2 in the air can react with free radicals (such as O • and OH • and) to form acidic and can cause acid rain. The acid rain in the earth will cause various diseases such as cough, throat irritation, itching and other tropical diseases and etc. In addition, the reaction of SO2 with free radicals cause depletion of ozone layer, which will give rise to a global heat (global warming) of the earth. Technology of the SO2 removal can be done by using a dry process and wet process. The removal of SO2 using the dry process is more economical when compared with wet processes because of simple, easy maintenance and absence liquid waste. However, the weakness of the dry process is the low efficiency of gas removal and sorbent utilization. One way to overcome the above problems is to use absorbent having high reactivity and an understanding of the mechanism of SO2 absorption in the solid sorbent with the formation of mathematical models. This research is expected to result in: (1) sorbent with a high reactivity and (2) mathematical models that can describe the phenomenon of SO2 absorption using solid sorbent. The tarfet of the research is to be used for absorbing gas pollutants from incinerator in effective, easy and inexpensively.

iii

PRAKATA Segala Puji dan rasa syukur kita panjatkan kehadirat Allah SWT, dengan berkat rahmat, dan hidayahNya kami telah dapat melaksanakan penelitian yang berjudul “Peningkatan reaktifitas absorbent Ca(OH)2 untuk penyisihan SO2 hasil pembakaran sampah menggunakan bag filter reaktor: Percobaan dan Modeling”.

Penelitian ini dapat terlaksana dengan dukungan dana dari program hibah kompetitif penelitian dengan Dibiayai oleh Direktorat Jenderal Pendidikan Tinggi, Departemen Pendidikan Nasional Sesuai dengan Surat Perjanjian Pelaksanaan Hibah Kompetiotif Penelitian Sesuai Prioritas Nasional Nomor: 212/SP2H/PP/DP2M/V/2009, Tanggal 30 Mei 2009.

Sebagai suatu penelitian yang berhubungan dengan perubahan iklim dan pelestarian lingkungan akibat adanya proses pembakaran yang terjadi di alam, pengamatan terhadap variabel proses yang berhubungan dengan peningkatan reaktifitas sorbent sehingga dapat menyerap gas-gas hasil pembakaran serta model matematik yang dapat menggambarkan fenomena pemisahan gas menggunakan sorben padat telah disusun.

Pada kesempatan ini, kami juga menyampaikan penghargaan dan terima kasih yang tak terhingga kepada semua pihak yang telah mendukung kegiatan penelitian ini, terutama jurusan Teknik Kimia Universitas Syiah Kuala, para laboran, dan teknisi di laboratorium Proses dan Operasi Teknik Kimia, Balai Riset dan Standarisasi (Baristan) Provinsi NAD serta para mahasiswa.

Banda Aceh, Nopember 2009

Dr. Ir. Mariana, M.Si

iv

DAFTAR ISI HALAMAN PENGESAHAN...............................................................................................i A. LAPORAN HASIL PENELITIAN RINGKASAN .......................................................................................................................ii SUMMARY ..........................................................................................................................iii PRAKATA............................................................................................................................ivi DAFTAR ISI ........................................................................................................................v DAFTAR TABEL.................................................................................................................vi DAFTAR GAMBAR ............................................................................................................vii DAFTAR LAMPIRAN .........................................................................................................ix BAB I PENDAHULUAN ....................................................................................................1 BAB II TINJAUAN PUSTAKA...........................................................................................3 2.1 Persiapan Sorbent Ca(OH)2 yang Mempunyai Reaktifitas yang Tinggi dengan Penambahan Additive DE ...............................................................................................3 2.2 Model Matematik ............................................................................................................4 2.2.1 Penurunan model matematik terhadap absorpsi di dalam reaktor bag filter..........5 BAB III TUJUAN DAN MANFAAT PENELITIAN ..........................................................12 3.1 Tujuan Penelitian ......................................................................................................12 3.2 Manfaat Penelitian ....................................................................................................12 BAB IV METODE PENELITIAN .......................................................................................13 4.1 Persiapan Sorbent Ca(OH)2 yang mempunyai Reaktifitas yang Tinggi dengan Penambahan Additive DE...............................................................................................13 4.1.1 Penghalusan sampel tanah diatomeae....................................................................15 4.1.2 Pengayakan sampel tanah diatomeae.....................................................................15 4.2 Tes Aktifitas (Penyisihan SO2 menggunakan sorbent Ca(OH)2/DE ...............................16 4.3 Tahap Pengujian Model Matematik ................................................................................18 BAB V HASIL DAN PEMBAHASAN................................................................................20 5.1 Tanah Diatomeae ............................................................................................................20 5.2 Peningkatan Reaktifitas Sorbent .....................................................................................21 5.2.1 Pengaruh reaksi terhadap densitas .........................................................................21 5.2.2 Pengaruh reaksi terhadap pH .................................................................................22 5.2.3 Pengaruh reaksi terhadap berat sampel..................................................................22 5.2.4 Pengaruh temperatur terhadap konversi ................................................................23 5.2.5 Pengaruh laju putaran pengaduk terhadap konversi ..............................................24 5.2.6 Pengaruh waktu reaksi terhadap konversi .............................................................24 5.2.7 Analisa TG-DTA ...................................................................................................26 5.2.8 Analisa XRD..........................................................................................................29 5.3 Pengujian Reaktifitas Sorbent .........................................................................................33 5.4 Pengembangan Model Matematik...................................................................................33 5.5 Pengujian Keabsahan Model ..........................................................................................33 BAB VI KESIMPULAN DAN SARAN ..............................................................................35 DAFTAR PUSTAKA ...........................................................................................................39 LAMPIRAN..........................................................................................................................42 B. DRAF ARTIKEL ILMIAH

v

DAFTAR TABEL

Tabel 5.1 Analisa hasil tanah diatomea Ujung Batee, Aceh Besar ...........................19

vi

BAB I PENDAHULUAN

Jutaan sampah dibuang ke lingkungan setiap harinya. Sebagian besar sampah ditumpuk dalam bentuk landfill di permukaan tanah, dibuang ke sungai atau laut tanpa penenganan terlebih dahulu. Sampah merupakan masalah yang paling serius terutama untuk daerah perkotaan yang berpenduduk pesat. Salah satu penanganan sampah dapat dilakukan melalui proses pembakaran di dalam sebuah incinerator (alat pembakaran sampah).

Gas-gas hasil pembakaran sampah dalam incinerator unumnya mengandung: (1) debu; (2) gas-gas seperti sulfur dioksida (SO2), nitrogen oksida (NOx), carbon monoksida (CO), asam hydroklorik; (3) dan logam-logam. Diantara gas-gas tersebut di atas, yang paling berbahaya adalah SO2 karena dapat berekasi dengan radikal bebas di udara membentuk asam yang dapat menyebabkan terjadinya hujan asam dan penipisan lapisan ozone sehingga menyebabkan terjadinya panas global bumi serta dapat menyebabkan berbagai penyakit seperti iritasi pernapasan, batuk dan sebagainya yang dapat meningkatkan angka kematian (Pearlman dkk., 1971; Goldstein dkk., 1973; Bates, 1972; Fairchild dkk, 1972).

Berbagai teknologi proses penyisihan SO2 di dalam gas telah dilakukan. Proses yang umum digunakan adalah proses kering dan proses basah. Penggunaan proses kering mempunyai beberapa keunggulan dari proses basah antara lain: (1) murah; (2) mudah penanganan; dan (3) tidak menghasilkan limbah cair. Salah satu proses pemisahan SO2 secara kering yang efektif, sederhana, mudah penanganan dan murah adalah menggunakan reaktor bag filter dengan menggunakan absorbent Ca(OH)2. Namun demikian proses ini juga mempunyai beberapa kelemahan antara lain adalah rendahnya konversi sorbent dan effisiensi pemisahan gas (Mariana dkk., 2003; Seno dkk., 1996; Uchida dkk., 1979).

Untuk meningkatkan efisiensi pemisahan gas dan meningkatkan reaktifitas sorbent (konversi kalsium), berbagai variasi bahan tambahan (additive) dapat ditambahkan ke dalam sorbent. Sorbent yang diperoleh dengan mencampur kalsium (Ca(OH)2) di dalam air dengan bahan mengandung silika dapat meningkatkan konversi kalsium dibandingkan dengan koversi kalsium yang diperoleh tanpa campuran additive (Garea dkk., 1996). Tanah

1

diatome (DE) atau fly ash (FA) merupakan material pozzolonic yang mengandung SiO2, Al2O3, dan CaO (Jozewicz dkk., 1988). penelitian FA sebagai sumber silika telah banyak dilakukan tetapi penelitian penggunaan DE sebagai sumber silika hampir tidak ada.

Disamping hal tersebut di atas, peningkatan efesiensi gas dapat juga dilakukan dengan pemahaman terhadap mekanisme penyerapan SO2 di dalam absorbent dengan pembentukan model matematik. Berbagai model matematik absorpsi gas di dalam absorbent telah dikembangkan. Model matematik untuk peneyarapan SO2 dalam batu kapur pertama sekali dikembangkan oleh Klingspor dkk, (1983, 1984). Pengembangan model matematik tersebut didasarkan pada shrinking core model dengan order reaksi satu untuk reaktor fixed bed. Mura dkk. (1991) mempelajari kenetika reaksi pemisahan sulfur dengan CaO pada sebuah fixed bed reactor menggunakan model grain. Georgakis dkk. (1979) mengembangkan model perubahan ukuran grain untuk reaksi gas – padat. Mariana dkk. (2003) telah mengembangkan model matenatik untuk penyerapan SO2 di dalam Ca(OH)2. Model ini belum mempertimbangkan pengaruh reaktifitas Ca(OH)2 terhadap penyisihan SO2. Penelitian ini difokuskan pada dua kajian utama yaitu: (1) persiapan absorbent Ca(OH)2 yang mempunyai reaktifitas yang tinggi dengan penambahan additive tanah diatome (DE); dan (2) pengembangan model matematika yang telah dikembangkan oleh Mariana dkk. (2003) untuk memprediksi penyisihan SO2 di dalam gas hasil pembakaran dengan menggunakan absorben Ca(OH)2/DE dengan mempertimbangkan pengaruh reaktifitas sorbent terhadap penyisihan gas. Untuk membuktikan keabsahan model yang akan dikembangkan, maka akan dibandingkan dengan data percobaan menggunakan gas SO2 yang disetimbangkan dengan gas N2.

2

BAB II TINJAUAN PUSTAKA

2.1 Persiapan Sorbent Ca(OH)2 yang Mempunyai Reaktifitas yang Tinggi dengan Penambahan Additive DE Berbagai literature mengemukakan bahwa reaksi pozzolanic antara silika dengan Ca(OH)2 di dalam larutan (slurry) membentuk kalsium silikat hidrat (hydrated calcium silicates (CSHs)) dan kalsium alumina (CAHs) yang tinggi yang dapat meningkatkan reaktivitas dari sorbent Ca(OH)2 (Liu dkk., 2002; Ho dan Shih, 1992; Ishizuka dkk., 2000). Gambar 2.1 menunjukkan model reaksi pozzolanic Ca(OH)2 dengan silika yang mengandung DE.

Hydrated Calcium Silicate (CSHs) Ca(OH) 2

D.E. Slurring

Gambar 2.1 Model reaksi pozzolanic Ca(OH)2 dengan silika yang mengandung DE (Mariana, 2003)

Beberapa literature juga menyebutkan bahwa reaktifitas sorbent sangat tergantung pada luas permukaan kontak sorbent (specific surface area) (Klinspor dkk., 1983; 1984). Secara jelas pengaruh luas permukaan kontak terhadap laju reaksi dengan menggunakan Ca(OH)2 sebagai sorbent padat telah direport oleh beberapa literatur (Irabien dkk., 1992; Ortiz dkk., 1993). Tsuchiai dkk. (1995) menyebutkan bahwa penghilangan SO2 maksimum diperoleh ketika pada waktu sebuah sorbent mencapai diameter pori rata-rata (mean pore diameter) maksimum, tetapi pada saat itu luas permukaan kontak spesifik (specific surface area) belum mencapai maksimum. Renado dkk. (1999) menunjukkan bahwa volume mesopore dan macropore merupakan komponen yang paling utama yang berhubungan dengan reaktifitas sorbent dan konversi penyisihan gas.

3

Lin dan Shih (2003) menyebutkan bahwa peningkatan luas permukaan kontak spesifik (specific surface area) sangat dipengaruhi oleh peningkatan volume total pori. Jozewicz dkk. (1987) menyebutkan bahwa peningkatan spesies hidrat (hydrated species) yang terdapat di dalam permukaan sorbent

dapat meningkatkan reaktifitas padatan dan

diharapkan untuk menghilangkan hidrasi air dengan pengontakan sorbent dan gas panas untuk menghasilkan pori baru yang mempunyai luas permukaan kontak baru. Peristiwa ini umumnya dikenal dengan proses kalsinasi.

Jika reaktifitas sorbent sangat tergantung pada luas permukaan kontak spesifik (specific surface area) dan volume total pori, maka sangatlah penting untuk meningkatkan luas permukaan kontak spesifik dan volume pori dari sebuah sorbent padat melalui proses kalsinasi.

2.2 Model Matematika

Gas hasil pembakaran sampah megandung polutan-polutan yang sangat berbahaya seperti sulfur dioksida (SO2), partikel padatan, dan lain-lain. Perkembangan proses yang ekonomis untuk menghilangkan polutan-polutan tersebut sangat dibutuhkan. Untuk tujuan tersebut di atas, sebuah reaktor bag filter dengan menggunakan sorbent Ca(OH)2 sangat menjanjikan karena merupakan proses yang sederhana, mudah penanganan, ekonomis dan tidak perlu penanganan limbah cair.

Namun demikian, proses ini mempunyai kelemahan yaitu

kecilnya konversi dari sorbent dan rendahnya efisiensi pemisahan gas (Uchida dkk., 1979; Seno dkk., 1996).

Untuk menyelesaikan masalah kecilnya efisiensi pemisahan gas SO2, perlu mempelajari prinsip penyisihan SO2 dan kinetika penyerapan/absorpsi. Berbagai model matematik untuk absorpsi gas di dalam sorbent padat telah dikembangkan. Model matematik yang pertama sekali dikembangkan oleh Klingspor dkk., (1983; 1984) untuk rekasi SO2 dan batu kapur (CaCO3) di dalam sebuah reaktor fixed bed menggunakan shrinking core model. Mura dkk. (1991) mempelajari kinetika penghilangan SO2 proses kering dengan menggunakan sorbent CaO dalam sebuah reaktor fixed bed berdasarkan model grain. Georgakis dkk. (1979) mengembangkan model matematik berdasarkan grain untuk reaksi gas-padat. Ramachandran dan Smith (1977) mempelajari pengaruh perubahan porositas

4

terhadap reaksi gas-padat berdasarkan konsep grain dan pellet (kumpulan grain), sedangkan Sampath dkk. (1975) juga mengembangkan model matematik untuk reaksi gaspadat non-catalytic melalui analisa transient yang berdasarkan model grain-pellet.

Mariana dkk. (2003) telah mengembangkan model matematik gas-padat di dalam reaktor bag filter dengan mengasumsi bahwa reaktor bag filter sama dengan reaktor fixed bed. Model grain digunakan untuk satu dimensi reaktor fixed bed. Kinetika reaksi SO2 dengan Ca(OH)2 dilakukan pada berbagai variabel kelembaban (humidity) dan temperature. Perpindahan gas melalui lapisan produk (product layer) dimodelkan dengan pertimbangan konveksi gas, sedangkan perpindahan gas melalui lapisan (bulk) gas menuju permukaan reaksi grain dimodelkan berdasarkan pertimbangan batasan difusi dan tahanan difusi padatan.

2.2.1 Penurunan Model Matematik terhadap Absorpsi di dalam Reaktor Bag Filter Mari kita pertimbangkan reaksi berikut: SO2 (gas) + Ca(OH)2 (padat)

CaSO3 (padat) + H2O

Untuk menyederhanakan model, maka dilakukan berbagai asumsi antara lain: (1) partikel tersusun dari grain yang sama; (2) baik partikel maupun grain mempunyai ukuran yang sama; (3) diameter grain meningkat dengan pembentukan reaksi dan total porositas partikel menurun; (4) konsentrasi gas dan profil konversi padatan di sekitar partikel berubah hanya ke arah reaktor vertikal; (5) partikel dimodelkan pada koordinat berbentuk bola (spherical) satu dimensi/arah; (6) tidak ada dispersi axial atau difusi sepanjang arah vertikal; (7) model isothermal; dan (8) tidak ada gradien tekanan sepanjang arah vertikal.

5

Grain:

Product

Time Unreacted core

Concentration of gas reactant

Reactant

rG 0

C Ag C AS

CA C AS

rG r rGC rGC r rG Radial Position Gambar 2.2a Schematik yang menunjukkan model grain

Berdasarkan asumsi-asumsi di atas dan dengan asumsi bahwa perkiraan kondisi tunak (steady) untuk phenomena difusi seperti ditunjukkan pada Gambar 2.2a sampai dengan Gambar 2.2c dan dapat dibuat persamaan berikut:

Particle:

CAg =CZ

R , CAg

R0

Gambar 2.2b Schematik yang menunjukkan model partikel/pellet

6

Particle layer: Gas reactant

Z= 0

C=CZ0

Z0

Z = Z0

C=CZ

Gambar 2.2c Schematik yang menunjukkan model partikel layer

2.2.1.1 Grain Jika sebuah grain dikontrol oleh reaksi kimia, reaksi order satu irreversible, lapisan difusi debu dan difusi antar molekul diabaikan, laju reaksi grain dapat ditulis sebagai berikut: 4πrGc k s C A g 2

N GA =

r k 1 + Gc s Ds

 r  1 − Gc rG 

  

.

(2-1)

Menurut Levenspiel (1972), penurunan volume atau diameter un-reacted core dapat menurunkan massa grain, N GA sebagai berikut: N GA = −4πrGc ρ absorbent 2

drGc . dt

(2-2)

Substitusi Persamaan (2-2) ke dalam Persamaan (2-1) memberikan laju reaksi untuk diameter un-reacted core sebagai:

k s C Ag / ρ absorbent drGc =− , dt rGc k s  rGc  1 −  1+ Ds  rG 

(2-3)

dimana nilai awal pada t = 0 , rGc = rG = rGo .

7

Persamaan (2-3) dapat dibuat variabel tak berdimensi dan koefisien sebagai berikut (Mariana dkk. (2001; 2003)):

g1 =

rG r k , g 2 = Gc , τ ′ = s t , rGo rGo rGo

φG =

C Ao

ρ absorbent

dg 2 = −φ G dτ ′

C Ag Ds ,C= , k s rGo C Ao

, β=

βC  g  β + g 2  1 − 2  g1  

(2-3’)

pada τ ′ = 0 , g 1 = g 2 = 1

Karena volume produk padat lebih besar dari volume reaktan padat, maka hubungan antara padat-padat dan gas-padat sangat dibutuhkan (Ulrichson dkk., 1980). Hal ini dapat didefinisikan sebagai parameter α dan dapat dihitung sebagai:

α=

ρ absorbent . ρ solid product

Koreksi antara g1 dan g2 dapat didefinisikan sebagai: g 1 = α + (1 − α )g 2 . 3

3

(2-4)

Konversi padatan sorbent di dalam grain dapat dihitung dengan persamaan:

xG = 1 − g 2 (τ ′, R ′) . 3

(2-5)

2.2.1.2 Partikel/pellet

Jika proses diasumsi isothermal dan pseudo-steady state, distribusi konsentrasi gas di dalam partikel/pellet dapat diuraikan dengan persamaan diferensial berikut dengan kondisi batas (Mura dkk. 1992): 1 ∂  2 ∂C A g  R De ∂R R 2 ∂R 

pada R = 0 ,

∂C Ag ∂R

  − rA = 0 ,  

(2-6)

=0,

pada R = Ro , C Ag = C z . Laju reaksi gas A di dalam partikel ( rA ) dapat dinyatakan dengan menghitung laju reaksi grain dan jumlah grain bola (sphere) per unit volume partikel adalah nol. Jumlah

8

(1 − ε Po ) / (43 πrGo3 ) diasumsi kontant selama reaksi (Ramachandran dan Doraiswamy, 1982). Oleh karena itu laju reaksi gas A di dalam partikel ( rA ) sebagai berikut: rA = N GA

atau

rA =

1 − ε Po 4 3 πrGo 3

3(1 − ε Po ) rGo

3

2

rGc k s C Ag r k 1 + Gc s Ds

 r  1 − Gc rG 

  

.

(2-7)

Dengan menggunakan persamaan tak berdimensi berikut, kita dapat mengurangi Persamaan (2-6) dengan bentuk tak berdimensi berikut: 2 C Ag 3(1 − ε P o )Ro k s De R ,δ = ,C= , φp = , R′ = Ro Deo C Ao Deo rGo

1 ∂  2 ∂C   R′ δ  = φp ∂R ′  R ′ 2 ∂R ′ 

atau

δ (R ′ )

βg 2 2 C

(2-6’)

 g  β + g 2  1 − 2  g1  

∂ 2C 2 ∂C ∂C ∂[δ (R ′)] + δ (R ′) + = φp 2 R ′ ∂R ′ ∂R ′ ∂R ′ ∂R ′

[

]

β g 2 2 (R ′) C 2

 g (R ′)   β + g 2 (R ′) 1 − 2 g 1 (R ′ )  

(2-8)

Kondisi batas Persamaan (2-6’) adalah pada R ′ = 0 ,

∂C =0, ∂R ′

pada R ′ = 1 , C = 1 . ∂[δ (R ′)] Dengan menggunakan korelasi ∂ C = x , = y , dan ∂R ′ ∂R ′ f1 = φ p

β [g 2 (R ′ )]2 C

 g (R ′ )  β + g 2 (R ′)1 − 2  g 1 (R ′ )  

,

Persamaan diferensial berikut diperoleh: ∂x 1  2  = f 1 − δ (R ′) x − xy  .  ∂R ′ δ (R ′)  R′ 

(2-9)

Akhirnya, konversi partikel dapat dihitung dari laju konversi grain berikut 1

x p = 3 ∫ xG R ′ 2 d R ′ .

(2-10)

0

9

2.2.1.3 Partikel layer (lapisan partikel)

Dengan mengasumsi bahwa situasi reaktor bag filter sama dengan reaktor fixed bed, maka neraca massa gas di dalam partikel layer dapat dinyatakan dengan persamaan berikut: u

∂C z ∂C z +εz = −rv , ∂z ∂t

(2-11)

Dengan kondisi awal dan kondisi batas: padat z = 0 , t ≥ 0 , C z = C zo , pada z > 0 , t = 0 , C z = 0 . Analog dengan laju reaksi di dalam partikel, laju reaksi gas A di dalam pertikel layer ( rv ) dapat dirumuskan sebagai rv =

1−εz N PA . 4 3 πR o 3

(2-12)

Laju reaksi di dalam partikel ( N PA ) merupakan fungsi laju reaksi sebuah grain ( N GA ) dan dapat dinyatakan sebagai berikut: N PA =

1 4 πRo 3 3

R

2 ∫ N GA .4πR dR = 0

R

3 Ro

3

∫N

GA

.R 2 dR .

(2-13)

0

Model persamaan di dalam partikel layer dapat dinyatakan sebagai:

∂C ∂C z 9(1 − ε z ) 0 =− u z +εz N GA .R 2 dR . 6 ∫ ∂z ∂t 4πRo 0 R

(2-14)

Untuk penyederhanaan simulasi, Persamaan (2-14) dapat dapat dibuat persamaan tak berdimensi berikut: C Ao

′

9 (1 − ε z )rGo k s z o k zε C z , = z , z′ = , γ = s o z , φ z = − 3 C zo zo rGo u uRo 2

′ 1  ∂C ′  βg 2 2 C ∂C Ao ′ R ′ 2 dR ′ . + γ  Ao  = φ z C Ao ∫   ′ ′ ∂z ∂τ  g  0   β + g 2  1 − 2  g1  

(2-14’)

Dengan kondisi batas ′ pada z ′ = 0 , τ ′ ≥ 0 , C Ao = 1 , ′ pada z ′ > 0 , τ ′ = 0 , C Ao = 0 .

10

1

′ ′ Dengan menggunakan f  C Ao  = φ z C Ao ∫   0

βg 2 2 C  g  β + g 2  1 − 2  g1  

, dan metoda forward difference,

Persamaan (2-14’) menjadi ′ C Ao (i , j + 1) =

∆z ′ ∆z ′ 1+γ ∆τ ′ C ′ (i + 1, j ) , + ∆τ ′ C ′ (i − 1, j ) + φ C ′ (i , j ) f (i . j ) Ao Ao z Ao 2 2

1−γ

(2-15)

dimana i dan j merupakan konstanta untuk metode numeric menggunakan metode perkiraan finite difference (Riggs, 1988).

Konversi sorbent padat di dalam partikel layer merupakan konversi keseluruhan proses dan dapat dihitung sebagai grain dan konversi partikel layer dapat dihitung dengan: 1

x z = ∫ x p dz ′ .

(2-16)

0

Persamaan (2-3’) dan (2-6’) dapat dihitung menggunakan metode Runge-Kutta-Gill method, sedangkan Persamaan (2-14`) dapat diselesaikan dengan perkiraan finite difference (Riggs, 1988).

11

BAB III TUJUAN DAN MANFAAT PENELITIAN

3.1 Tujuan Penelitian

Tujuan dari penelitian ini adalah: 1) Mempersiapkan absorbent Ca(OH)2 yang mempunyai reaktifitas yang tinggi. Persiapan absorben tersebut dilakukan dengan penambahan additive tanah diatome (DE) ke dalam sorben Ca(OH)2 2) Mengembangkan model matematika yang telah dikembangkan oleh Mariana dkk. (2003) untuk memprediksi penyisihan SO2 di dalam gas hasil pembakaran dengan menggunakan absorben Ca(OH)2/DE. Pengembangan model matematik tersebut dengan mempertimbangkan pengaruh reaktifitas sorbent terhadap penyisihan gas. 3) Untuk membuktikan keabsahan model yang dikembangkan, maka hasil perhitungan menggunakan model matematik dibandingkan dengan data percobaan menggunakan gas SO2 yang disetimbangkan dengan gas N2.

3.2 Manfaat Penelitian

Adapun manfaat dari penelitian adalah: 1) Menghasilkan sorbent yang reaktif sehingga dapat digunakan untuk mengabsorp gas-gas yang berbahaya dari proses pembakaran sampah 2) Menghasilkan model matematik yang dapat fenomena pemisahan gas-gas yang berbahaya dengan menggunakan sorbent padat 3) Hasil penelitian ini sangat memberi konstribusi dalam pengembangan ilmu Neraca Massa dan Energi yang menuju pada pengembangan Ipteks, sehingga dapat menambah wawasan dan ketrampilan peneliti seta sebagai wahana pertukaran informasi perkembangan penelitian. 4) Melalui aktivitas penelitian dan pendidikan yang berkesinambungan diharapkan dapat menangani problem-problem yang dihadapi oleh industriindustri kimia terkait yang ada di Indonesia dan Nanggroe Aceh Darussalam khususnya.

12

BAB IV METODE PENELITIAN

4.1. Persiapan Sorbent Ca(OH)2 yang mempunyai Reaktifitas yang Tinggi dengan Penambahan Additive DE

Persiapan bahan yang diperlukan meliputi pengadaan bahan-bahan berikut: 1. Ca(OH)2 2. Tanah diatome 3. Aquades

Tanah diatome dihaluskan dan diayak sesuai dengan variable proses yang dilakukan. Adapun ukuran partikel tanah diatomea yang dipersiapkan adalah: (1) Ukuran 150 (2) Ukuran 250 mesh (3) Ukuran 350 mesh

Sedangkan peralatan yang dipersiapkan meliputi: 1. Erlenmeyer 2. Pipet volum 3. Aluminium foil 4. Corong pemisah 5. Gelas ukur 6. Magnetic stirrer/water batch 7. Termometer

Kapur hidrat (hydrated lime) yang digunakan dalam percobaan ini adalah Ca(OH)2 murni (99.9% - Ca(OH)2), sedangkan tanah diatome yang digunakan berasal dari NAD.

13

1. Ca(OH)2 + additive

1

2. Water 3. Mixer 5

4. Water bath 5. Heater 6. Magnetic stirrer

4 2 3 6

Gambar 4.1a Persiapan sorbent (Skema)

Gambar 4.1b Persiapan sorbent yang sedang dilakukan

Persiapan sorbent ditunjukkan pada Gambar 4.1a dan Gambar 4.1b. Kapur hidrat dan additive DE, bersama-sama dengan 100 g air dimasukkan ke dalam beaker polypropylene conical pada berbagai variasi rasio Ca(OH)2/additives dengan rasio berat air/padatan

(Ca(OH)2+additives) adalah 10. Beaker kemudian ditutup dan dimasukkan ke dalam was water bath pada temperatur 65oC. Kemudian diaduk dengan magnetic stirrer selama 2 jam. Setelah diaduk, slurry dikeringkan selama 24 jam pada temperature 120oC dan solid kemudian dikalsinasi pada berbagai variasi temperature selama 2 jam. Cake hasil pengeringan dihaluskan dan kemudian setelah dilakukan analisa spesifik surface area

14

(SBET), ditabur pada bag filter untuk test aktifitas. Komposisi sorbent dianalisa dengan menggunakan XRD system Rint 2000, Rigaku Co. Japan. SEM, Hitachi S-3000N digunakan untuk melihat marfologi sorbent. Padatan sorbent dianalisa menggunakan analisa thermogravimetric.

4.1.1 Penghalusan sampel tanah diatomea Penghalusan tanah diatomea dilakukan dengan menggunakan Ball Mill di Laboratorium Proses dan Operasi Teknik Kimia Unsyiah seperti ditunjukkan pada Gambar 4.2. Penghalusan ini dilakukan untuk mempersiapkan tanah diatomea dalam berbagai varisi ukuran partikel.

Gambar 4.2 Persiapan penghalusan sampel tanah diatome menggunakan crusher

4.1.2 Pengayakan sampel tanah diatomea

Proses pengayakan sampel diatome dilakukan dengan menggunakan Ayakan Shaker di Laboratorium Proses dan Operasi Teknik Kimia Unsyiah (Gambar 4.3).

Pengayakan ini bertujuan untuk mempersiapkan ukuran partikel tanah diatome seperti yang diharapkan. Ukuran partikel yang dipersiapkan adalah: -

150 mesh

-

250 mesh

-

350 mesh

15

Gambar 4.3. Proses pengayakan menggunakan shaker sedang dijalankan 4.2. Test Aktifitas (Penyisihan SO2 menggunakan Sorbent Ca(OH)2/DE)

Percobaan untu reaksi sorbent dengan SO2 dilakukan dalam sebuah reaktor bag filter (Gambar 4.2a dan Gambar 4.2b). Aliran gas yang mengandung 1000 ppm SO2 dan N2, dengan laju alir 2 liter/min dilewatkan sebuah pemanas listrik untuk menaikkan temperature menjadi 150 oC dan kemudian direaksikan dengan sorbent padat di dalam sebuah filter. Jumlah sorbent dibuat constant sebesar 0.8 g dengan waktu reaksi dilakukan selama 1200 sec.

16

(1) (2) (3) (4)

Gas Heater Reaktor Temperatur controller

NaOH

NaOH 1. 2. 3. 4. 5. 6. 7. 8. 9.

Particle layer

Reaktor Filter Ekeltrik Furnace Micro Feeder Gas Mixer Flow Meter Three Way Cock H2O Absorber SO2 Analyzer

Thimble filter

Gambar 4.2a Perangkat percobaan (Skema)

(4) (1)

1. 2. 3. 4.

Gas Heater Bag Filter Reaktor Temperatur Controler

(3) (2)

Gambar 4.2b Perangkat percobaan tes aktifitas sorbent

17

4.2.1 Analisa Konsentrasi SO2

Analisa konsentrasi gas SO2 dilakukan dengan menggunakan dengan menangkap terlebih dahulu SO2 hasil percobaan dengan menggunakan larutan induk (penangkap) sebagai berikut: Untuk Basis 1000 ml Larutan: 1.

HgCl2 = 10.86 gram

2.

KCl

3.

EDTA = 0.066 gram

= 5.96 gram

+ 1000 ml aquabidest

Hasil tangkapan selanjutnya dianalisa dengan menggunakan metode Parrosanilin di Laboratorium Balai Riset dan Standarisasi Industri (Baristan) Provinsi NAD.

4.3. Tahap Pengujian Model Matematika

Setelah melakukan percobaan dan mengembangkan model matematik, maka proses selanjutnya yang akan dilakukan adalah membandingkan hasil percobaan dengan hasil prediksi menggunakan model matematika yang dikembangkan untuk melihat keabsahan model tersebut. Bila selisih antara hasil prediksi model dan percobaan (relative error) lebih kecil dari 15%, maka dapat disimpulkan bahwa model yang dikembangkan dapat menggambarkan fenomena pemisahan SO2 dengan sorbent Ca(OH)2/DE. Sebaliknya bila relative error lebih besar dari 15%, maka model perlu dikoreksi sampai menghasilkan relative error lebih kecil dari 15%. Relative error dapat dihitung dengan menggunakan

korelasi berikut:

 Penyis ihan SO2 prediksi  % Relative error =  − 1 ×100%  Penyisihan SO2 percobaan 

(4-1)

Adapun diagram alir keseluruhan penelitian yang akan dilakukan baik pada proses batch maupun proses kontinyu ditunjukan pada Gambar 4.3.

18

Ca(OH)2

DE

Slurry

Air

Diaduk

Dikeringkan

Persiapan Sorbent

Dihaluskan Ca(OH)2/DE

SO2

Reaktor

Model matematika

Analisa: Konsentrasi SO2

Prediksi % kehilangan Cr(III)

% kehilangan Cr(III)

Analisa XRD Analisa TG-DTA Analisa SiO2 sisa

Pengetesan Aktifitas

>15%

% Relative error <15% Model dapat menggambarkan fenomena penghilangan Cr(III) menggunakan lautan alkali

Gambar 4.3 Diagram alir keseluruhan penelitian yang dilakukan

19

BAB V HASIL DAN PEMBAHASAN

5.1 Tanah Diatomeae

Adapun tanah diatome yang digunakan berasal dari Desa Ujung Batee, Aceh Besar dan diperlihatkan pada Gambar 5.1. Komposisi hasil analisa tanah diatome tersebut ditunjukkan pada Tabel 5.1. Sedangkan hasil analisa dari tanah diatome dari Dinas Perindustrian dan Perdagangan (Baristan) NAD dapat dilihat pada Lampiran A.

Gambar 5.1 Tanah Diatomea dari Desa Ujung Batee, Aceh Besar

Tabel 5.1 Analisa hasil tanah diatomea Ujung Batee, Aceh Besar No

Parameter Uji

Metode Uji

Satuan

Hasil Uji

1 2 3

SiO2 CaO Al2O3

Gravimetri Titrimetri Gravimetri

% % %

33.94 T.T 4.89

Sumber: Baristan, NAD, 22 Mei 2009

Tabel 5.1 menunjukkan bahwa kandungan SiO2 di dalam tanah diatome adalah sebesar 33,94%. Sehingga dapat disimpulkan bahwa tanah diatome dari Kreung Raya Aceh Besar

20

ini layak digunakan sebagai sumber silika untuk meningkatkan reaktifitas sorbent Ca(OH)2.

5.2 Peningkatan Reaktifitas Sorben 5.2.1 Pengaruh reaksi terhadap densitas

Komposisi suatu larutan sangat mempengaruhi terhadap densitas larutan. Adapun hasil analisa densitas awal dan densitas akhir larutan reaksi antara Ca(OH)2 dan tanah diatomeae (DE) pada sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 350 mesh hasil percobaan diperlihatkan pada Gambar 5.2.

350 mesh

Gambar 5.2 Hasil analisa densitas awal dan densitas akhir hasil percobaan menggunakan sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 350 mesh

Tabel 5.2 menunjukkan bahwa hasil analisa densitas akhir setelah reaksi sedikit lebih tinggi dari densitas awal. Hal ini menjunjukkan bahwa didalam larutan sudah terkandung kalsium silicate hidrate (CSH) hasil reaksi berikut:

Ca(OH)2 + SiO2 + H2O

CaO.SiO2.2H2O

(5.1)

(Kalsium silikat hidrat)

21

5.2.2 Pengaruh reaksi terhadap pH

Komposisi suatu larutan sangat mempengaruhi terhadap pH larutan. Adapun hasil analisa pH awal dan pH akhir larutan reaksi antara Ca(OH)2 dan tanah diatomeae (DE) pada sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 350 mesh hasil percobaan diperlihatkan pada Gambar 5.3.

Gambar 5.3 Hasil analisa pH awal dan pH akhir hasil percobaan menggunakan sorben

Ca(OH)2/DE = 10:1 dan ukuran partikel 350 mesh Hasil analisa pH sebelum dan sesudah reaksi diperlihatkan pada Gambar 5.3. Gambar 5.3 menunjukkan bahwa pH akhir lebih tinggi dari pH awal. Hal ini membuktikan bahwa terbentuknya hasil reaksi antara Ca(OH)2 dengan SiO2 membentuk CaO.SiO2.2H2O.

5.2.3 Pengaruh reaksi terhadap berat sampel

Komposisi suatu larutan sangat mempengaruhi terhadap berat sampel yang dihasilkan. Adapun hasil analisa berat awal dan berat akhir larutan reaksi antara Ca(OH)2 dan tanah diatomeae (DE) pada sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 350 mesh hasil percobaan diperlihatkan pada Gambar 5.4.

Gambar 5.4 menunjukkan bahwa berat sampel akhir sedikit lebih tinggi dari berat sampel awal. Hal ini membuktikan bahwa terbentuknya hasil reaksi antara Ca(OH)2 dengan SiO2 membentuk CaO.SiO2.2H2O yang lebih ringan karena berbentuk hidrate.

22

Gambar 5.4 Hasil analisa berat sampel awal dan densitas akhir hasil percobaan menggunakan sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 350 mesh 5.2.4 Pengaruh temperatur terhadap konversi

Pengaruh temperatur terhadap reaksi menggunakan sorben Ca(OH)2/DE = 3:1, 350 rpm dan ukuran partikel 350 mesh ditunjukkan pada Gambar 5.5.

Gambar 5.5 Pengaruh temperature terhadap konsentrasi Ca(OH)2 pada sorben Ca(OH)2/DE = 10:1, ukuran partikel 200 mesh

Gambar 5.5 menunjukkan bahwa peningkatan temperatur reaksi menyebabkan peningkatan konversi sampai mencapai 65oC dan menurun kembali pada temperature reaksi 80oC. Gambar 3.5 juga menunjukkan bahwa temperatur reaksi 65oC menghasilkan konversi tertinggi. Hal ini sesuai dengan hasil penelitian yang telah dilakukan oleh Mariana dkk.

23

(2003) yang menyatakan bahwa temperature 65oC merupakan temperature optimum untuk reaksi Ca(OH)2 dengan SiO2 dalam DE. 5.2.5 Pengaruh laju putaran pengaduk terhadap konversi

Pengaruh laju putaran pengaduk terhadap konversi menggunakan sorben Ca(OH)2/DE = 3:1, ukuran partikel 350 mesh dan temperature reaksi 65oC ditunjukkan pada Gambar 5.6. Gambar 5.6 menunjukkan bahwa semakin besar laju putaran pengaduk, maka konversi Ca(OH)2 yang bereaksi semakin meningkat.

Gambar 5.6 Pengaruh laju putaran pengaduk terhadap konversi pada sorbent Ca(OH)2/DE = 3:1, ukuran partikel 350 mesh, temp. reaksi 56oC 5.2.6 Pengaruh waktu reaksi terhadap konversi

Pengaruh waktu reaksi terhadap konversi menggunakan sorben Ca(OH)2/DE = 3:1, ukuran partikel 350 mesh dan 350 rpm ditunjukkan pada Gambar 5.7 sampai 5.10.

Gambar 5.7 Pengaruh waktu terhadap konversi pada sorbent Ca(OH)2/DE = 3:1, ukuran partikel 350 mesh, 350 mesh dan T = 50oC

24



Gambar 5.7 sampai Gambar 5.9 menunjukkan bahwa waktu reaksi mempunyai pengaruh yang significant terhadap konversi. Konversi semakin meningkat terhadap waktu reaksi sampai 20 menit, selanjutnya setelah 20 menit konversi mempunyai harga yang konstan. Hal ini menunjukkan bahwa waktu reaksi optimum antara Ca(OH)2 dengan SiO2 adalah 20 menit dalam range variable percobaan yang dilakukan.

25

Gambar 5.10 Pengaruh waktu terhadap konversi pada sorbent Ca(OH)2/DE = 3:1, ukuran partikel 350 mesh, 350 mesh pada berbagai variasi temperature

Gambar 5.10 menunjukkan bahwa temperature reaksi 60oC mempunyai konversi reaksi yang lebih tinggi pada berbagai waktu reaksi. Hal ini relevan dengan Gambar 3.5 yang menunjukkan bahwa konversi tertinggi diperoleh pada temperature reaksi 65oC.

5.2.7 Analisa TG-DTA

Analisa TG-DTA dilakukan Lembaga Ilmu Pengetahuan Indonesia, Pusat Penelitian Fisika Komplek PUSPIPTEK Serpong. Analisa TG-DTA menggunakan sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh pada temperature 65oC, 70oC dan 75oC ditunjukkan pada Gambar 5.11 sampai Gambar 5.13. Gambar 5.11 sampai Gambar 5.13 menunjukkan sorben terdiri dari Ca(OH)2, SiO2 dan CaO.SiO2.2H2O. Hal ini menunjukkan bahwa reaksi antara SiO2 dengan Ca(OH)2 menghasilkan CaO.SiO2.2H2O, sehingga sorbent Ca(OH)2/DE lebih vorus, sehingga lebih reaktif dari sorben Ca(OH)2. Hal ini sesuai dengan hasil penelitian yang telah dilakukan oleh Mariana, dkk (2003) yang untuk Ca(OH)2/FA.

26

CaO.SiO2.2H2O SiO2 Ca(OH)2

Gambar 5.11 TG-DTA sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh pada temperature reaksi 65oC

27

CaO.SiO2.2H2O SiO2 Ca(OH)2

Gambar 5.12 TG-DTA sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh pada temperature reaksi 70oC

28

SiO2

CaO.SiO2.2H2O Ca(OH)2

Gambar 5.13 TG- sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh pada temperature reaksi 75oC

5.2.8 Analisa XRD (X-ray Diffraction Patterns)

Untuk melihat adanya produk CaO.SiO2.2H2O hasil reaksi dapat dilihat dari analisa XRD. Analisa XRD dilakukan Lembaga Ilmu Pengetahuan Indonesia, Pusat Penelitian Fisika Komplek PUSPIPTEK Serpong. Gambar 5.11 sampai Gambar 5.13 menunjukkan analisa XRD sorben Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh untuk variasi temperatur

29

reaksi. Dari Gambar 5.10 sampai Gambar 5.12 dapat dilihat bahwa hasil reaksi antara Ca(OH)2 dengan SiO2 menghasilkan CSH (CaO.SiO2.2H2O).

CSHS (CaO.SiO2.2H2O Ca(OH)2 SiO2

Gambar 5.11 Profil XRD untuk sorbent Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh pada temperature reaksi 65oC

30


Gambar 5.12 Profil XRD untuk sorbent Ca(OH)2/DE = 10:1 untuk produk pada temperature reaksi 70oC

31


Gambar 5.13 Profil XRD untuk sorbent Ca(OH)2/DE = 10:1 untuk produk pada temperature reaksi 75oC

32

5.3 Pengujian Reaktifitas Sorben

Reaktifitas sorben dapat dilihat dengan menggunakan bag filter reaktor. Adapun hasil reaksi antara SO2 dengan sorben Ca(OH)2/DE ditunjukkan pada Gambar 5.14.

Gambar 5.14 Profil XRD untuk sorbent Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh 5.4 Pengembangan Model Matematik

Pengembangan model matematik didasarkan pada model matematik yang telah dikembangkan oleh mariana dkk (2003). Model matematik terdiri dari grain, partikel dan partikel layer. Peningkatan reaktifitas sorbent dilakukan dengan penambahan SiO2 ke dalam sorben Ca(OH)2, maka untuk menyelesaikan model matematik menggunakan sifatsifat fisik sorbent campuran berikut:

ln η m = ∑i x i ln µ i + ∑∑ x i x j Gij

(5-1)

Gij = ∑ ∆ i − ∑ ∆ j + W

(5-2)

i≠ j

W=

(0.3161)(N i − N j )2 Ni + N j

− (0.1188)(N i − N j )

(5-3)

5.5 Pengujian Keabsahan Model

Pengujian model matematik dilakukan dengan data percobaan. Adapun hasil perbandingan antara data percobaan dan perhitungan dengan menggunakan model matematik yang dikembangkan ditunjukkan pada Gambar 5.15.

33

Gambar 5.15 Perbandingan antara hasil percobaan dengan hasil perhitungan untuk sorbent Ca(OH)2/DE = 10:1 dan ukuran partikel 200 mesh

Gambar 5.15 menunjukkan bahwa hasil percobaan mempunyai kecendrungan yang sama dengan hasil perhitungan. Hal ini menunjukkan bahwa model matematik yang dikembangkan dapat menggambarkan fenomena pemisahan gas dengan menggunakan sorbent padat.

34

BAB VI KESIMPULAN DAN SARAN

6.1

Kesimpulan

Dari hasil pembahasan penelitian ini dapat disimpulkan bahwa temperature dan laju putaran pengaduk mempunyai peranan yang penting terhadap konversi reaksi antara SiO2 dengan Ca(OH)2 sehingga dapat meningkatkan reaktifitas sorben Ca(OH)2/DE. Temperatur reaksi 60oC menghasilkan konversi tertinggi, sehingga dapat disimpulkan bahwa temperature reaksi 65oC merupakan temperature optimum untuk reaksi Ca(OH)2 dengan SiO2. Hasil penelitian juga menunjukkan bahwa konversi mencapai konstan pada waktu reaksi 20 menit.

Model matematik yang dikembangkan dapat menggambarkan pemisahan gas menggunakan sorben padat. Hasil perhitungan menggunakan model yang dikembangkan mempunyai kesesuaian dengan data percobaan, menunjukkan bahwa model yang dikembangkan dapat menggambarkan fenomena pemisahan gas menggunakan sorben padat.

6.2

Saran

Berdasarkan hasil dan pembahasan pada penulisan laporan ini, disarankan untuk di lakukan penelitian selanjutnya untuk lebih mendetail pada pemisahan gas menggunakan sorben untuk secara simultan seperti untuk gas HCL, CO dan HCl/CO. Demikian juga perlu pengembangan model yang bisa menggambarkan fenomena pemisahan gas secara simultan menggunakan sorben padat.

Ucapan terima kasih

Dengan selesainya peneltian ini kami mengucapkan terima kasih atas bantuannya sehingga penelitian ini berjalan dengan baik. Terima kasih kepada member Lab Operasi Teknik Kimia, Unsyiah. Terima kasih khusus kepada Direktorat Jenderal Pendidikan Tinggi, Departemen Pendidikan Nasional yang telah membiayai penelitian ini sesuai dengan Surat Perjanjian Pelaksanaan Hibah Kompetitif Penelitian Sesuai Prioritas Nasional Nomor : 212/SP2H/PP/DP2M/V/2009, Tanggal 30 Mei 2009.

35

Daftar Lambang dan Notasi

[cm2]

A

= luas permukaan filter (surface area of filter)

C

= konsentrasi gas A tanpa dimensi, ( C Ag / C Ao )

CA

= konsentrasi gas A (SO2)

[moles/cm3]

CAg

= konsentrasi gas A (SO2) dalam partikel

[moles/cm3]

C Ao

= konsentrasi gas A awal

[moles/cm3]

′

[–]

= koncentrasi gas A tak berdimensi, ( C z / C zo )

[–]

CAS

= koncentrasi gas A (SO2) di pembukaan grain

[moles/cm3]

Cz

= koncentrasi gas A di lapisan partikel (particle layer)

[moles/cm3]

C zo

= koncentrasi gas A di perm. lapisan pertikel (particle layer surface)

[moles/cm3]

dP

= diameter of partikel

D AB

= koefisien difusi gas A melalui gas B

[cm2/s]

De

= koefisian difusi efektif gas A

[cm2/s]

Deo

= koefisian difusi efektif gas A awal ( D NA (ε Po ) )

[cm2/s]

D KA

= koefisien difusi Knudsen gas A

[cm2/s]

D NA

= koefisien difusi gas A

[cm2/s]

Do

= faktor pre-eksponensial Arrhenius untuk koefisien difusi

[cm2/s]

Ds

= koefisien difusi gas dalam padatan grain

[cm2/s]

C Ao

[cm]

2

E r , E d = energi aktivasi reaksi dan difusi

[kJ/mole]

g1 , g 2

= jari-jari grain tak berdimensi, ( rG / rGo , rGc / rGo )

[–]

Gij

= parameter yang didefinisikan dengan Persamaan (5-2)

[-]

th

i

= i step waktu untuk reaksi

[t]

j

= jth step waktu jarak dari permukaan ke lapisan partikel

[z]

ko

= faktor pre-eksponensial Arrhenius untuk reaksi

ks

= laju konstanta reaksi

[cm/s]

N

= jumlah atom karbon

[-]

N GA

= laju reaksi gas A dalam grain

[moles/s]

N PA

= laju reaksi gas A dalam partikel

[moles/s]

rA

= laju reaksi gas A dalam partikel/pellet per unit volume

[m/s]

[moles/cm⋅s]

36

r , rG

= jarak radial dalam grain dan jari-jari grain

[cm]

rGo

= jari-jari awal grain ( 3 /[ρ absorbent SBET ] )

[cm]

rGc

= jari-jari un-reacted core dari grain

[cm]

rv

= laju reaksi gas A dalam lapisan partikel per unit volume

R

= jarak radial dalam partikel

[cm]

Rabsorbent = densitas bulk dari sorbent

[moles/cm3]

[moles/cm3⋅s]

Ro

= jari-jari partikel awal

R′

= jari-jari partikel tak berdimensi, ( R / Ro )

SBET

= luas permukaan partikel

t

= laju reaksi

u

= kecepatan superfisial gas A

W


[-]

Wabsorbent = jumlah sorbent dalam lapisan partikel (particle layer)

[mol]

x, y


[cm] [–] [cm2/g] [s] [cm/s]

[-]

xG , x p , x z = konversi dalam grain, partikel, dan lapisan partikel (particle layer)

[–]

z′

= jarak tak berdimensi lapisan partikel (particle layer), ( z / zo ),

[–]

z

= jarak dari permukaan lapisan partikel (surface of particle layer)

[cm]

zo

= tebal lapisan partikel (particle layer) ( Wabsorbent / [AR absorbent ])

[cm]

Greek letters

α

= koefisian volume ekspansi, ( ρ absorbent / ρ solidproduct )

β

=

ρ

= densitas

∆

= konstanta Grunberg dan Nissan

[-]

εP 

= porositas partikel setiap waktu t

[–]

ε Po

= porositas partikel/pellet awal

[–]

εz 

= porositas lapisan partikel (particle layer) = 1 − [Rabsorbent / ρ absorbent ]

[–]

φG

= konsentrasi tak berdimensi, (C Ao / ρ absorbent )

[–]

Ds k s rGo

[–] [–] [moles/cm3]

37

φP

3(1 − ε Po )Ro k s Deo rGo 2

=

9(1 − ε z )rGo k s z o

[–]

2

φz

=−

γ

=

k s zoε z rGo u

[−]

µ

= viskositas

[cP]

δ

= koefisien difusi efektif tak berdimensi gas A, ( De / Deo )

[–]

τ′

= waktu tak berdimensi, [k s t ] / rGo

[–]

uRo

3

[–]

Subcripts i

= komponen i

j

= komponen j

m

= campuran

38

DAFTAR PUSTAKA

Bates D.V., “Air pollutants and Human Lung: The James Waring Memorial Lecture”, Am. Rev Respir Dis., 105, pp. 111-13 (1972) Fairchild, G.A., J. Roan, and J. McCarroll, “Atmospheric Pollutants and the Pathogenesis of Viral Respiratory Infection, Arch Environ Health, 25, September, pp. 174-182 (1972) Garea, A., J.R. Viguri, and A. Irabien, “ Fly-Ash/Calcium Hydroxide Mixtures for SO2 Removal: Structural Properties and Maximum Yield,” Chem. Eng. J., 66, pp. 171179 (1996). Georgakis, C., C.W. Chang and J. Szekely, “A Changing Grain Size Model for Gas-Solid Reactions,” Chem. Eng. Sci., 34, pp.11072 (1979). Goldstein, E., M.E. Carroll, and P.D. Hoeprich, “Effect of Nitrogen Dioxide on Pulmonary Bacterial Defense Mechanism”, Erch Environ Health, 26, pp. 202-204 (1973) Ho, C.S and S.M. Shih, “Ca(OH)2/Fly Ash Sorbents for SO2 Removal,” Ind. Eng. Chem. Res., 31, pp. 1130-1135 (1992). Irabien, A., F. Cortabitarte, and I. Ortiz, “Kinetics of Flue Gas Desulfurization at Low Temperatutes: Nonideal Surface Adsorption Model,” Chem. Eng. Sci. 47, pp. 15331543 (1992). Ishizuka, T., H. Tsuchiai, T. Murayama, T. Tanaka, and H. Hattori, ”Preparation of Active Absorbent for Dry-Type Flue Gas Desulfurization from Calcium Oxide, Coal Fly Ash, and Gypsum,” Ind. Eng. Chem. Res., 39, pp. 1390-1396 (2000). Jozewicz, W.; Chang, J.; Brna, T.; Sedman, C. Reactivation of Solids from Furnace Injection of Limestone for SO2 Control. Environ. Sci. Technol. 21, 7, pp. 664-670 (1987). Jozewicz, W., J.C.S. Chang, C.B. Sedman, and T.G. Brna, “Silica-Enhanced Sorbents for Dry Injection Removal of SO2 from Flue Gas,” J. Air Pol. Cont. Ass., 38, pp. 10271034 (1988) Klingspor, J., H.T. Karlsson and I. Bjerle, “A Kinetic Study of the Dry SO2-Limestone Reaction at Low Temperature,” Chemical Engineering Communication, 22, pp. 81103 (1983). Klingspor, J., A. Stromberg, H.T. Karlsson and I. Bjerle, “Similarities Between Lime and Limestone in Wet-dry Scrubbing,” Chemical Engineering Processing, 18, pp. 239-

39

247 (1984). Levenspiel, O., “Chemical Reaction Engineering,” pp.363, Second ed., John Wiley & Sons, New York (1972). Lin, R.B. and S.M. Shih, “Characterization of Ca(OH)2/Fly Ash Sorbents for Flue Gas Desulfurization,” Powder Tech., 131, pp. 212-222 (2003). Liu, C.F., S.M. Shih, and R.B. Lin, “Kinetic of the Reaction of Ca(OH)2/Fly Ash Sorbent with SO2 at Low Temperatures,” Chem. Eng. Sci., 57, pp. 93-104 (2002). Mariana, C. Chen, Y. Tsujimura, A. Maezawa and S. Uchida, “Modeling of SO2

Absorption in a Bag Filter Based on Grain Model,” Proc. of Regional Symposium on Chemical Engineering, Bandung, Indonesia, pp. MO10.1-MO10.6 (2001) Mariana, Oobuchi, Y., Maezawa, A. and Uchida, S., “Ca(OH)2/Fly Ash Sorbent for HCl

Removal in a Bag Filter Reactor, “Proc. of Indonesian 12th Student Scientific Meeting in Japan, Osaka, Japan, pp. 809-813 (2003) Mariana, C. Chen, Y. Tsujimura, A. Maezawa and S. Uchida, “Experiments on SO2

Absorption in a Bag Filter Based on the Grain Model,” J. Chin. Inst. Chem. Engrs., 34, 2, pp. 2011-209 (2003).

Mura, G., A. Lallai, and P. Olla, “On the Kinetics of Dry Desulfurization with Calcium Oxide,” Chem. Eng. Journal, 46, pp.119-128 (1991). Mura, G. and A. Lallai, “On the Kinetics of Dry Reaction Between Calcium Oxide and Gas,” Chem. Eng. Sci., 47, 9-11, pp. 2407-2411 (1992) Ortiz, I., F. Cortabitarte, A. Garea, and A. Irabien, “Flue Gas Desulfurization at Low Temperatures: Characterization of the Structural Changes in the Solid Sorbent”, Powder Technol. 75, pp. 167-172 (1993). Pearlman, M.E., J.F. Finklea, J.P. Creason, C.M. Shy, M.M. Young, and J.M. Horton, ”Nitrogen Dioxide and Lower Respiratory Illness”, Pediatrics, 47, 2, pp. 391-398 (1971) Ramachandran, P.A. and J.M. Smith, “Effect of Sintering and Porosity Changes on Rates of Gas-Solid Reactions,” Chem. Eng. Journal, 14, pp.137 (1977). Ramachandran, P.A. and L.K. Doraiswamy, “Modeling of Noncatalytic Gas-Solid Reactions,” AIChE J., 28, 6, pp. 881 (1982) Renedo, M.J., J. Fernandez, A. Garea, A., A. Ayerbe, and J.A. Irabien, J.A. Microstructural Changes in the Desulfurization Reaction at Low Temperature. Ind. Eng. Chem. Res. 38, pp. 1384-1390 (1999)

40

Riggs, J.B., “An Introduction To Numerical Methods For Chemical Engineers,” Texas Tech University Press, USA (1988) Sampath, B.S., P.A. Ramachandran and R. Hughes, “Modeling of Non-Catalytic Gas-Solid Reactions-I. Transient Analysis of the Particle-Pellet Model,“ Chem. Eng. Sci., 30, pp.125 (1975). Seno, T., A. Maezawa, T. Fujiwara and S. Uchida, “Simultaneous Absorption of Acid Gases By Slaked Lime in Bag Filter System,” Proc. of Fourth Japan-Korea Symposium on Separation Technology, Tokyo, Japan, pp. 263-266 (1996). Tsuchiai, H, T. Ishizuka, T. Ueno, H. Hattori, and H. Kita, “Highly Active Sorbent for SO2 Removal Prepared from Coal Fly Ash,” Ind. Eng. Chem. Res. 34, pp 1404-1411 (1995) Uchida, S., S. Kageyama, M. Nogi and H. Karakida, “Reaction Kinetics of HCl and Limestone,” J. Chin. Inst. Chem. Engrs., 10, pp. 45-49 (1979). Ulrichson, L.D. and D.J. Mahoney, “Pore Closure Models for Gas-Solid Reactions: The Effect of Bulk Flow and Reversibility,” Chem. Eng. Sci. 35, pp. 567 (1980)

41

Lampiran 1 Foto Peralatan Penelitian

5 Mei 2009 Gambar L.1 Penghalusan tanah diatomeae (DE) sedang dilakukan (Gambar 1)

5 Mei 2009

Gambar L.2 Penghalusan tanah diatomeae (DE) sedang dilakukan (Gambar 2)

42

7 Mei 2009 Gambar L.3 Penyiapan tanah diatomeae (DE) menurut ukuran partikel yang diinginkan (Gambar 1)

5 Mei 2009 Gambar L.4 Penyiapan tanah diatomeae (DE) menurut ukuran partikel yang diinginkan (Gambar 2)

43

23 Mei 2009 Gambar L.5 Reaksi antara Ca(OH)2 dengan SiO2 menggunakan reaktor (Gambar 1)

25 Mei 2009 Gambar L.6 Reaksi antara Ca(OH)2 dengan SiO2 menggunakan reaktor (Gambar 2)

44

Gambar L.7 Sorben hasil reaksi antara Ca(OH)2 dengan SiO2

Gambar L.8 Reaktor tempat reaksi antara Ca(OH)2 dengan SiO2

45

Gambar L.9 Heater untuk menaikkan temperatur gas sebelum masuk reaktor

Gambar L.10 Alat pengontrol temperatur dalam reaktor

46

Gambar L.10 Alat pengontrol temperatur dalam reaktor

Gambar L.11 Valve pengatur gas N2 dan SO2 masuk reaktor

47

Gambar L.12 Tabung gas yang masuk reaktor

Gambar L.13 Contoh hasil analisa SiO2

48

Gambar L.14 Contoh hasil analisa TG-DTA (1)


49



50

Gambar L.17 Contoh hasil analisa XRD

51

Lampiran 2 Instrumen Penelitian Penelitian ini dilaksanakan di Laboratorium Operasi Teknik Kimia Universitas Syiah Kuala (Unsyiah). Sarana pendukung penelitian ini yang tersedia di Laboratorium Operasi Teknik Kimia Unsyiah dan Balai Riset dan Standarisasi Industri (Baristan) Provinsi NAD ditunjukkan pada Tabel L.13.

Tabel L.13 Sarana pendukung penelitian ini yang tersedia di Laboratorium Operasi Teknik Kimia Unsyiah dan Balai Riset dan Standarisasi Industri (Baristan) Provinsi NAD No

Peralatan yang ada

1

Oven driyer

2

Water batch

3

Heater

4

Erlenmeyer

5

pH meter

6

Analisa SO2

Kegunaan

1. Untuk mengeringkan slurry sorbent setelah dicampur dengan DE 2. Untuk mengkalsinasi sorbent Ca(OH)2/DE setelah dikeringkan 1. Untuk memanaskan slurry (air + Ca(OH)2 + DE) sehingga sesuai dengan yang diinginkan pada proses persiapan sorbent 2. Untuk menjaga temperatur slurry tetap konstant pada proses persiapan sorbent Untuk menaikkan temperatur water batch sehingga dapat meningkatkan temperatur slurry pada proses persiapan sorbent Wadah persiapan sorbent dan analisa hasil percobaan Untuk mengukur pH slurry sorebnt pada proses persiapan sorbent Untuk menganalisa konsentrasi SO2 setiap saat pada proses pengetesan aktifitas

Kemampuan

Keterangan

Dapat digunakan sampai temperature > 1000oC

Lab OTK Unsyiah

Sangat bagus

Lab OTK Unsyiah

Kondisinya bagus

Lab OTK Unsyiah

Jumlahnya cukup

Lab OTK Unsyiah

Ketelitiannya cukup bagus

Lab OTK Unsyiah

Ketelitiannya cukup bagus

Laboratorium Balai Riset dan Standarisasi Industri (Baristan) Provinsi NAD.

52

Disamping sarana pendukung di atas, Fakultas Teknik Unsyiah juga memiliki sebuah bengkel yang dapat dimanfaatkan untuk merangkai peralatan pengetesan keaktifan penyerapan SO2 oleh sorbent Ca(OH)2/DE. Teknik Kimia Unsyiah juga mempunyai Program Magister (S2). Hal ini sangat membantu pelaksanaan penelitian yang akan dilakukan baik dalam membantu mahasiswa S1 dalam percobaan maupun pengolahan dan analisa data dan pengembangan model matematik. Pengetesan XRD dan TG-DTA dilakukan di Lembaga Ilmu Pengetahuan Indonesia, Pusat Penelitian Fisika, Kompleks PUSPIPTEK, Serpong.

53

Lampiran-3 Curiculum Vitae Tim Peneliti

CURRICULUM VITAE KETUA PENELITI 1.

Nama Lengkap dan Gelar

: Dr. Ir. Mariana, M.Si.

2.

Jenis Kelamin

:P

3.

Tempat/Tanggal Lahir

: Aceh Besar/15 Juli 1967

4.

Alamat

: Desa Nusa, Kecamatan Lhoknga, Aceh Besar

5.

NIP

: 19670715 199303 2 003

6.

Pangkat/Golongan

: Penata TK. I/IIId

7.

Jabatan Fungsional

: Lektor

8.

Jabatan Struktural

: Staf Pengajar Fakultas Teknik Jurusan T. Kimia Unsyiah

9.

Jurusan/prodi

: Teknik Kimia

10. Fakultas

: Teknik

11. Riwayat Pendidikan

:

No

PENDIDIKAN

IJAZAH/TAHUN

SPESIALISASI

1

S1 (Teknik Kimia Unsyiah)

1992

Teknik Kimia

2

S2 (Teknik Kimia ITB)

1997

Teknologi Proses, Pemisahan dan Modeling

3

S3 (Material Scinece & Chemical Engineering, Shizuoka University, Hamatasu, Jepang)

2004

Teknologi Proses, Pemisahan dan Modeling

12. Pengalaman Penelitian (5 Tahun Terakhir) No

JUDUL PENELITIAN

1

Treatment of Flue Gas from Incinerators using Dry and Wet Scrubber Systems

2

Waste Water Treatment (Metal Separation from Aqueous Solution by Alkaline Solution)

3

Cr(III) Removal from Aqueous Solution by Alkaline Solution (Batch process): Modeling and Experiment, Fundamental Research, DIKTI, 2008

TAHUN

SUMBER BIAYA

2004

Monbusho Scholarship, Jepang

2005 - 2006

Post Doctoral Shizuoka University, Jepang

2008

Fundamental Research, DIKTI, 2008

54

13. Publikasi (5 Tahun Terakhir) No

JUDUL PENELITIAN

TAHUN

1

Experimental and Modeling Study on CO2 Absorption in a Cyclone Scrubber by Phenomenological Model and Neural Networks, Korean J. Chem. Eng., 21,3, pp. 589-594, 2004

2004

2

Ca(OH)2/Diatomaceous Earth Sorbents for HCl Removal in a Bag Filter Reactor J. Chin. Inst. Chem. Engrs., 35, 4, pp.285-288, 2004

2004

3

Preparation of Ca(OH)2/Diatomaceous Earth Sorbents by Calcination for HCl Removal in a Bag Filter Reactor, The 3rd Hokkaido Indonesian Student Association Scientific Meeting, Hokkaido, Japan, March, 2004

2004

4

Modeling Study of Dust Collection Efficiencies in a Wet Cyclone Scrubber, Annual Meeting of The Society of Separation Process Engineers, Japan (SSPEJ), Tokyo, Japan, June, 2004

2004

5

odeling of Cr(III) Removal from Aqueous Solution by Alkaline Solution, Indonesian 14th Indonesian Student Scientific Meeting in Japan, Nagoya, Japan, September, 2005

2005

6

Cr(III) Removal from Aqueous Solution by Alkaline Solution, J. Chem. Eng. Japan., 39, 7, pp. 724-730, 2006

2006

7

Cr(III) Removal from Aqueous Solution by Alkaline Solution: Modeling and Theoretical Study, Proc. of Rgional Symposium on Computational Mechanics and Numerical Analysis (CMNA-2006), Banda Aceh, Indonesia, May, 2006.

2006

8

Effect of pH on Cr(III) Removal from Aqueous Solution by Alkaline Solution: Modeling and Theoretical Study, Proc. of Annual symposium of Chemical Engineering Syiah Kuala University, Banda Aceh, Indonesia, September, 2006.

2006

9

Cr(III) Removal from Aqueous Solution by Alkaline Solution, J. Chem. Eng. Japan., 39, 7, pp. 724-730, 2006

2006

10

Gas Absoprption by Alkaline Solution in a Cylcone Scrubber : Experimental and Modeling, Jurnal Rekayasa Kimia dan Lingkungan Vol. 6 No. 1 , Hal. 1- 6, 2007

2007

11

HCl Removal By Ca(OH)2/Fly Ash Sorbent in A Bag Filter Reactor, Jurnal Rekayasa Kimia dan Lingkungan Vol. 6 No. 2 , Hal. 60- 63, 2007

2007

12

Pengaruh pH terhadap Penghilangan Cr(III) dari Limbah Cair menggunakan Larutan Alkali: Modeling dan eksperiment Studi, Jurnal Rekayasa Kimia dan Lingkungan Vol. 7 No. 2 , Hal. 67 - 72, 2008

2008

Banda Aceh, 25 Nopember 2009 Ketua Peneliti,

Dr. Ir. Mariana, M.Si NIP 19670715 199303 2 003

55

CURRICULUM VITAE ANGGOTA PENELITI 1. Nama Lengkap dan Gelar

: Dr. Ir. Syahrul, M.Sc.

2. Jenis Kelamin

: L

3. Tempat/Tanggal Lahir

:

4. Alamat

:

5. NIP

: 19631211 199002 1 001

6. Pangkat/Golongan

:

7. Jabatan Fungsional

: Lektor

8. Jabatan Struktural

: Staf Pengajar Fakultas Pertanian Jurusan Teknik Pertanian Unsyiah

9. Jurusan/prodi

: Teknik Pertanian

10. Fakultas

: Pertanian


:

No

PENDIDIKAN

1

S1 (Teknik Pertanian IPB)

2

S2 (

3

S3 (

IJAZAH/TAHUN

SPESIALISASI

) )


JUDUL PENELITIAN

TAHUN

SUMBER BIAYA

1 2 3 13. Publikasi (5 Tahun Terakhir) No

JUDUL PENELITIAN

TAHUN

1 2 3 4 5 6 7

56

8 9 10 11 12 Banda Aceh, 25 Nopember 2009 Ketua Peneliti,

Dr. Ir. Syahrul, M.Sc NIP 19631211 199002 1 001

57

CURRICULUM VITAE ANGGOTA PENELITI 1. Nama Lengkap dan Gelar

: Dr. Farid Mulana, ST., M.Eng.

2. Jenis Kelamin

:L

3. Tempat/Tanggal Lahir

: Aceh Besar, 8 Februari 1972

4. Alamat

: Jl. Prada 1 Lr. Seulanga No. 20, Kp. Pineung, B. Aceh

5. NIP

: 19720208 199702 1 001

6. Pangkat/Golongan

: Penata Muda Tk. 1/IIIb

7. Jabatan Fungsional

:

8. Jabatan Struktural Unsyiah

: Staf Pengajar Fakultas Teknik Jurusan T. Kimia

9. Jurusan/prodi

: Teknik Kimia

10. Fakultas

: Teknik


:

No

PENDIDIKAN

1

S1 (Teknik Kimia Unsyiah)

2

S2 (

3

S3 (

IJAZAH/TAHUN

SPESIALISASI

) )


JUDUL PENELITIAN

TAHUN

SUMBER BIAYA

1 2 3 13. Publikasi (5 Tahun Terakhir) No

JUDUL PENELITIAN

TAHUN

1 2 3 4 5 6 7

58

8 9 10 11 12 Banda Aceh, 25 Nopember 2009 Ketua Peneliti,

Dr. Farid Mulana, ST., M.Eng. NIP 19720208 199702 1 001

59

LAPORAN HASIL PENELITIAN SESUAI PRIORITAS NASIONAL

Recommend Documents