PENENTUAN PARAMETER KINETIK DEKOMPOSISI KATALITIK METANA DENGAN KATALIS Ni-Cu-Al UNTUK PRODUKSI NANO KARBON SKRIPSI

UNIVERSITAS INDONESIA

PENENTUAN PARAMETER KINETIK DEKOMPOSISI KATALITIK METANA DENGAN KATALIS Ni-Cu-Al UNTUK PRODUKSI NANO KARBON

SKRIPSI

MARCHO RIZAL 0806333272

FAKULTAS TEKNIK TEKNIK KIMIA DEPOK JUNI 2012

Penentuan parameter..., Marcho Rizal, FT UI, 2012

i

UNIVERSITAS INDONESIA

PENENTUAN PARAMETER KINETIK DEKOMPOSISI KATALITIK METANA DENGAN KATALIS Ni-Cu-Al UNTUK PRODUKSI NANO KARBON

SKRIPSI Diajukan sebagai salah satu syarat untuk memperoleh gelar Sarjana Teknik

MARCHO RIZAL 0806333272

FAKULTAS TEKNIK TEKNIK KIMIA DEPOK JUNI 2012


HALAMAN PERNYATAAN ORISINALITAS

Skripsi ini adalah hasil karya saya sendiri, dan semua sumber baik yang

dikutip maupun dirujuk telah saya nyatakan dengan benar.

Nama

: Marcho Rizal

NPM

: 0806333272

Tanda Tangan : Tanggal

: 26 Juni 2012

ii Penentuan parameter..., Marcho Rizal, FT UI, 2012

HALAMAN PENGESAHAN

Skripsi ini diajukan oleh Nama NPM Program Studi Judul Skripsi

: : Marcho Rizal : 0806333272 : Teknik Kimia : Penentuan Parameter Kinetik Dekomposisi Katalitik Metana dengan Katalis Ni-Cu-Al untuk Produksi Nano Karbon

Telah berhasil dipertahankan di hadapan Dewan Penguji dan diterima sebagai bagian persyaratan yang diperlukan untuk memperoleh gelar Sarjana Teknik pada Program Studi Teknik Kimia, Fakultas Teknik, Universitas Indonesia

DEWAN PENGUJI

Pembimbing 1

: Dr. Ir. Praswasti P.D.K Wulan, M.T.

Pembimbing 2

: Dr . rer. nat. Ir. Yuswan Muharam, MT

Penguji

: Dr. Ir. Tania Surya Utami, M.T.

Penguji

: Ir. Mahmud Sudibandriyo MSc, PhD

Penguji

: Ir. Rita Arbianti MSi

Ditetapkan di Tanggal

: Departemen Teknik Kimia FTUI, Depok : 26 Juni 2012

iii Penentuan parameter..., Marcho Rizal, FT UI, 2012

KATA PENGANTAR

Puji syukur penulis panjatkan kepada Tuhan Yang Maha Esa atas berkat dan rahmat-Nya sehingga penulis dapat menyelesaikan skripsi ini dengan sebaik baiknya. Skripsi dengan judul “Metode Alternatif Penentuan Angka Oktana

Bahan Bakar Komersial Dengan Menggunakan Model Kinetika Pembakaran Hidrokarbon Multikomponen”

ini disusun sebagai salah satu syarat untuk

mencapai gelar Sarjana Teknik Jurusan Teknik Kimia pada Fakultas Teknik Universitas Indonesia. Dalam penulisan skripsi ini, penulis berterima kasih atas bantuan dan dukungan yang telah diberikan oleh berbagai pihak. Secara khusus, penulis mengucapkan terima kasih yang sedalam-dalamnya kepada : 1. Prof. Dr. Ir. Widodo Wahyu Purwanto, DEA selaku ketua Departemen Teknik Kimia Fakultas Teknik Universitas Indonesia. 2. Dr. Ir. Praswast PDK Wulan selaku pembimbing skripsi 3. Dr. rer. nat. Ir. Yuswan Muharam, MT selaku pembimbing skripsi 4. Orang tua tercinta atas semua dukungan doa dan kasih sayang serta segala hal yang telah diberikan. 5. Rekan-rekan mahasiswa Teknik Kimia angkatan 2008 atas semua kerja sama dan bantuannya selama ini serta Seluruh teman-teman saya terutama Chandra Hadiwijaya, Kenny Viriya, Shelly Apsari, dan Fiona Wijaya yang sangat membantu dalam penyelesaian skripsi ini. 6. Semua staf dan karyawan di Departemen Teknik Kimia Universitas Indonesia atas segala bantuannya. 7. Serta berbagai pihak yang tidak dapat disebutkan satu persatu untuk segala kontribusinya. Harapan penulis skripsi ini dapat bermanfaat secara keilmuan bagi setiap orang yang membacanya. Kritik dan saran yang membangun sangat diharapkan oleh penulis demi tercapainya hasil yang lebih baik. Depok, 2012 Penulis

iv Penentuan parameter..., Marcho Rizal, FT UI, 2012

HALAMAN PERNYATAAN PERSETUJUAN PUBLIKASI TUGAS AKHIR UNTUK KEPENTINGAN AKADEMIS Sebagai sivitas akademik Universitas Indonesia, saya yang bertanda tangan di bawah ini:

Nama

: Marcho Rizal

NPM

: 0806333272

Program Studi

: Teknik Kimia

Departemen

: Teknik Kimia

Fakultas

: Teknik

Jenis karya

: Skripsi

demi pengembangan ilmu pengetahuan, menyetujui untuk memberikan kepada Universitas Indonesia Hak Bebas Royalti Noneksklusif (Non-exclusive RoyaltyFree Right) atas karya ilmiah saya yang berjudul : “Penentuan Parameter Kinetik Dekomposisi Katalitik Metana dengan Katalis Ni-

Cu-Al untuk Produksi Nano Karbon” beserta perangkat yang ada (jika diperlukan). Dengan Hak Bebas Royalti Noneksklusif

ini

Universitas

Indonesia

berhak

menyimpan,

mengalihmedia/formatkan, mengelola dalam bentuk pangkalan data (database), merawat, dan memublikasikan tugas akhir saya selama tetap mencantumkan nama saya sebagai penulis/pencipta dan sebagai pemilik Hak Cipta. Demikian pernyataan ini saya buat dengan sebenarnya.

Dibuat di

: Depok

Pada tanggal : 26 Juni 2012 Yang menyatakan

(Marcho Rizal)

v Penentuan parameter..., Marcho Rizal, FT UI, 2012

ABSTRAK Nama

: Marcho Rizal

Program Studi : Teknik Kimia Judul

: Penentuan Parameter Kinetik Dekomposisi Katalitik Metana

dengan Katalis Ni-Cu-Al untuk Produksi Nano Karbon

Nano Teknologi semakin berkembang sekarang ini. Reaksi dekomposisi katalitik metana merupakan salah satu cara untuk memproduksi nano karbon. Kendala dari proses produksi nano karbon dengan menggunakan reaksi dekomposisi katalitik metana berada pada parameter kinetika yang belum diketahui terutama jika memperhitungkan deaktivasi katalis. Penelitian ini dilakukan untuk mengembangkan model serta menentukan parameter kinetik dari reaksi-reaksi elementer dekomposisi katalitik metana dengan mengikutsertakan reaksi deaktivasi katalis yang akan diselesaikan dengan menggunakan regresi non linear, metode Marquadt. Model menghasilkan nilai pra eksponensial yang berkisar antara 0.00042 hingga 34.3, dengan kisaran nilai energy aktivasi antara 50,274 kJ/mol hingga 104,673kJ/mol

Kata kunci : nano karbon, dekomposisi katalitik metana, kinetik, regresi non linear, Marquadt

vi Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

ABSTRACT Name

: Marcho Rizal

Major of Study : Chemical Engineering

Title

: Determination of Kinetic Parameter for Catalytic Decomposition

of Methane with Ni-Cu-Al Catalyst for Carbon Nano Production

Nano technology is advancing right now. Catalytic Decomposition of Methane is one of many ways to produce Nano Carbon. The problem with this method is the kinetic parameter especially if we include the deactivation of the catalyst used. This research is done to develop a model to determine the parameter kinetic of catalytic decomposition of methane that includes the deactivation of catalyst and will be solved by using a non linear regression, the Marquadt method. This research gives the pre exponensial number ranging from 0.00042 to 34.3, with activation energy ranging from 50,274 kJ/mol to 104,673kJ/mol

Keywords : Nano carbon, catalytic decomposition of methane, kinetic, non linear regression, marquadt

vii Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

DAFTAR ISI HALAMAN PERNYATAAN ORISINALITAS .................................................... ii

HALAMAN PENGESAHAN ................................................................................ iii KATA PENGANTAR ........................................................................................... iv

HALAMAN PERNYATAAN PERSETUJUAN PUBLIKASI .............................. v ABSTRAK ............................................................................................................. vi

DAFTAR ISI ........................................................................................................ viii DAFTAR TABEL .................................................................................................. ix DAFTAR GAMBAR .............................................................................................. x DAFTAR SIMBOL................................................................................................ xi DAFTAR LAMPIRAN ........................................................................................ xiii BAB 1 PENDAHULUAN ...................................................................................... 1 1.1 1.2 1.3 1.4

Latar Belakang .......................................................................................... 1 Rumusan Masalah ..................................................................................... 3 Tujuan Penelitian ...................................................................................... 3 Ruang Lingkup Masalah ........................................................................... 3

BAB 2 TINJAUAN PUSTAKA ............................................................................. 4 2.1 2.2 2.3 2.4 2.5 2.6

Nano Karbon ............................................................................................. 4 Cara Memproduksi CNT........................................................................... 5 Reaksi Dekomposisi Metana..................................................................... 6 Kondisi Operasi ...................................................................................... 10 Matrix Laboratory (MatLab)................................................................... 11 Regresi Non Linear ................................................................................. 12

BAB 3 METODE PENELITIAN.......................................................................... 15 3.1 3.2

Diagram Alir Penelitian .......................................................................... 15 Rincian Kegiatan Penelitian .................................................................... 16

BAB 4 HASIL DAN PEMBAHASAN................................................................. 21 4.1 Mekanisme Reaksi .................................................................................. 21 4.2 Model Kinetika dan Fitting Model terhadap Data Eksperimen .............. 28 4.3 Penentuan Parameter Kinetik (Nilai Pra Eksponensial Konstanta Laju Reaksi dan Energi Aktivasi).................................................................................. 32 BAB 5 KESIMPULAN ......................................................................................... 36 DAFTAR PUSTAKA ........................................................................................... 37

viii Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

DAFTAR TABEL Tabel 3. 1 Tabel Model Persamaan Laju Reaksi Masing-Masing spesies zat ...... 17 Tabel 4. 1 Tabel nilai Konstanta Laju Reaksi untuk nilai k1 hingga k11

……32

Tabel 4. 2 Tabel nilai K dan Deviasi Standar pada suhu 750oC .......................... 33

Tabel 4. 3 Tabel nilai K dan Deviasi Standar pada suhu 700oC .......................... 33

Tabel 4. 4 Tabel nilai K dan Deviasi Standar pada suhu 650oC .......................... 33 Tabel 4. 5 Tabel nilai K dan Deviasi Standar pada suhu 600oC .......................... 33

Tabel 4. 6 Tabel nilai K dan Deviasi Standar pada suhu 550oC .......................... 33 Tabel 4. 7 Tabel nilai K dan Deviasi Standar pada suhu 500oC .......................... 33 Tabel 4. 8 Tabel Nilai Pra Eksponensial dan Energi Aktivasi dari Masing-Masing k ............................................................................................................................. 34

ix Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

DAFTAR GAMBAR Gambar 2. 1 Logo MatLab .................................................................................. 12 Gambar 2. 2 Tampilan Antarmuka MatLab ........................................................ 12

Gambar 3. 1 Diagram Alir Penelitian…………………………………………...15

Gambar 3. 2 Algoritma Program Regresi Non Linear......................................... 19 Gambar 4. 1 Ilustrasi Keseluruhan Reaksi ……………………………………21

Gambar 4. 2 Skema Mekanisme Detail Dari Reaksi yang Digunakan Dalam Model .................................................................................................................... 22 Gambar 4. 3 Reaksi Pemutusan Hidrogen dari Metana Tahap 1......................... 23 Gambar 4. 4 Reaksi Pemutusan Hidrogen dari Metana Tahap 2......................... 23 Gambar 4. 5 Reaksi Pemutusan Hidrogen dari Metana Tahap 3......................... 24 Gambar 4. 6 Reaksi Pemutusan Hidrogen dari Metana Tahap 4......................... 24 Gambar 4. 7 Pembentukan Gas Hidrogen dari Atom H yang Terikat dengan Inti Aktif ...................................................................................................................... 25 Gambar 4. 8 Segregasi Karbon yang Terikat pada Katalis ke Bagian Depan Katalis ................................................................................................................... 25 Gambar 4. 9 Difusi Karbon Melalui Nikel .......................................................... 26 Gambar 4. 10 Pembentukan Nano Karbon .......................................................... 27 Gambar 4. 11 Deaktivasi Katalis karena Enkapsulasi Katalis .............................. 28 Gambar 4. 12 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=750 oC ......................................................... 29 Gambar 4. 13 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=700 oC ......................................................... 29

Gambar 4. 14 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=650 oC ......................................................... 30 Gambar 4. 15 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=600 oC ......................................................... 30 Gambar 4. 16 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=550 oC ......................................................... 31 Gambar 4. 17 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=500 oC ......................................................... 31

x Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

DAFTAR SIMBOL o ANi o

luas permukaan spesifik partikel nikel [m2/gcat]

CCNi,F konsentrasi karbon yang terlarut dalam nikel pada bagian depan

3

partikel [mol C/m ] o CCNi,R

konsentrasi karbon yang terlarut dalam nikel pada bagian belakang

partikel [mol C/m3] o DNi

Diameter partikel nikel [mm]

o DCNi

Difusivitas karbon di nikel [m2/h]

o

Perubahan energy pada segregasi karbon [J/Mol]

o (g)

fasa gas

o k1

konstanta laju reaksi pada tahap (1) [molekul/(site.sec.atm)]

o k2

konstanta

laju

reaksi

balik

pada

reaksi

tahap

(1)

[molekul/(site.sec)] o k3

konstanta laju reaksi pada tahap (2) [molekul/(site.sec.)]

o k4

konstanta

laju

reaksi

balik

pada

reaksi

tahap

(2)



o k6

konstanta

laju

reaksi

balik

pada

reaksi

tahap

(3)



o k8

konstanta

laju

reaksi

balik

pada

reaksi

tahap

(4)



o k10

konstanta

laju

reaksi

balik

pada

reaksi

tahap

(5)

[molekul/(site.sec.atm)] o k11

laju reaksi enkapsulasi pada tahap (9) [molekul/(site.sec)]

o L

panjang jalur difusi [m]

o PCH4

tekanan parsial metana [atm]

o PH2

tekanan parsial hydrogen [atm]

o R

konstanta gas universal (J/(mol.kelvin)]

o S

sisi aktif logam

xi Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

o T

suhu

o xf

konsentrasi karbon pada bagian depan pada partikel nikel [g C/g

Ni] o xR

konsentrasi karbon pada bagian belakang dari partikel nikel [g C/g

Ni] o θi.s

konsentrasi permukaan dari sisi yang diisi oleh spesies i [I = CH3,

CH2, CH, C, H]

o θt

total sisi aktif

o θv

sisi aktif yang tidak diisi spesies apapun

o ρNi

massa jenis nikel

o η

jumlah atom di permukaan logam

o Cv

jumlah sisiaktif kosong setiap unit massa katalis dibagi dengan

bilangan Avogadro [gmol/gcat] o Ci.s

konsentrasi permukaan yang diisi oleh spesies i [gmol/gcat]

xii Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

DAFTAR LAMPIRAN Lampiran A Penurunan persamaan model……………………………….. .39 Lampiran B Data Eksperimen yang Digunakan…………………………...47

Lampiran C Plot ln k Terhadap 1/T ………………………………………..55

Lampiran B Program yang digunakan untuk penelitian………………… 61

xiii Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

BAB 1 PENDAHULUAN 1.1

Latar Belakang

Nanoteknologi dan bahasa pemrograman komputer adalah dua hal yang semakin berkembang sekarang ini. Salah satu riset yang menarik tentang nanoteknologi adalah carbon nanotube (CNT) dari divisi nanokarbon. CNT

menjadi pusat sorotan media setelah ditemukan oleh Iijima pada tahun 1991. CNT sangat berpotensi untuk digunakan pada berbagai aplikasi seperti nanoteknologi, elektronik, optik, dan lainnya. CNT memiliki kekuatan yang luar biasa dan sifat elektrik yang unik serta merupakan konduktor panas yang efisien. Sifat-sifat tersebut membuat CNT dapat digunakan sebagai penyimpan hidrogen, nanoscale transistor, flat-panel display, superkapasitor, nanoprobes dan sensor (Daenan et al., 2003). Reaksi dekomposisi katalitik metana merupakan salah satu reaksi yang menjanjikan untuk mengkonversi unsur karbon menjadi material adi karbon nanotube melalui reaksi CH4  C + 2H2

(Chen et al, 2001).

Proses ini juga merupakan salah satu bagian dari chemical vapor deposition (CVD). Selain cara ini juga ada cara lain untuk memproduksi CNT yaitu arc discharge dan laser ablation. Kelebihan reaksi dekomposisi metana dibandingkan dengan kedua cara diatas adalah biaya yang lebih ekonomis. Dalam penelitian ini kinetika dari reaksi dekomposisi katalitik metana akan diteliti. Kinetika berhubungan dengan laju dari reaksi kimia dan bagaimana laju reaksi dipengaruhi oleh berbagai faktor seperti konsentrasi dan temperatur. Kinetika memberikan bukti mengenai mekanisme dari proses kimia yang terjadi. Hal ini akan membantu untuk meprediksi hasil dan jumlah dari reaksi kimia yang terjadi. Oleh sebab inilah analisa kinetika mikro digunakan untuk penelitian ini. Kinetika dari reaksi dekomposisi katalitik metana akan berbeda tergantung dari berbagai macam faktor seperti penggunaan jenis katalis, mekanisme yang ditinjau, dan dasar dari persamaan yang digunakan untuk menurunkan persamaan laju reaksi. Penelitian tentang produksi CNT dengan menggunakan reaksi dekomposisi metana sudah banyak dilakukan sebelumnya. Salah satu contohnya adalah penggunaan katalis berbeda akan menyebabkan perubahan pada kinetika,

INDONESIA Penentuan parameter..., Marcho1Rizal, FT UI,UNIVERSITAS 2012

2

seperti yang sudah dibandingkan antara alumina dengan karbon aktif (Zongqing, 2005). Peninjauan mekanisme yang digunakan juga akan mempengaruhi kinetika

dan penurunan persamaan laju reaksi. Dalam tim riset CNT Universitas Indonesia, sejak tahun 2005 telah

dilakukan penelitian mengenai reaksi dekomposisi katalitik metana dengan menggunakan katalis terstruktur gauze. Dalam perkembangannya akan dilakukan menggunakan katalis terstruktur pelat sejajar menggunakan katalis Ni-Cu-Al

untuk memperoleh luas permukaan yang tinggi dan menghasilkan CNT yang lebih seragam (alligned CNT) . Terkait dengan hal tersebut, maka dilakukan studi kinetika sebagai salah satu dasar perancangan, optimasi dan scale up reaktor. Kinetika reaksi dekomposisi metana yang akan dilakukan berbasis mekanisme reaksi (kinetika intrinsik) yang melibatkan tahap reaksi elementer (Chen et.al., 2001). Internal reaktor menggunakan bentuk granular berukuran 0,125 mm untuk mempermudah limitasi tahanan internal dan eksternal (Wulan, 2011). Parameterparameter yang tidak diketahui diestimasi melewati prosedur estimasi parameter dengan menggunakan data eksperimen dimana set parameter terbaik digunakan untuk meminimalisir hasil dari error jumlah kuadrat antara model dan data dipertimbangkan. Untuk mengestimasi parameter yang digunakan, digunakan metode regresi non linear. Dalam menyelesaikan masalah pemodelan dengan regresi non linear, terdapat berbagai metode yang dapat digunakan. Salah satunya adalah metode Marquadt. Metode Marquadt ini menggunakan teknik interpolasi antara metode Gauss-Newton dan metode steepest descent. Interpolasi didapat dengan menambahkan matriks diagonal, dimana nilai dari λ dipilih pada setiap iterasi

sehingga vektor parameter terkoreksi akan menghasilkan jumlah kuadrat yang lebih rendah pada iterasi. Metode Marquadt sangat efektif dalam penggunaan untuk pengujian dan validasi data, karena ketepatan, dimana jumlah dari kuadrat bilangan residual dapat dicari hingga nilai yang paling kecil. Penggunaan metode Marquadt pada penelitian ini karena metode Marquadt merupakan salah satu metode regresi non linear yang memiliki tingkat ketepatan yang sangat tinggi.


3

Pada penelitian ini, fenomena yang terjadi pada kinetika reaksi akan diekspresikan dalam persamaan matematis yang akan diselesaikan untuk mendapatkan informasi mengenai parameter kinetika dari seluruh reaksi sehingga bisa digunakan sebagai basis untuk desain reactor.

1.2

Rumusan Masalah

Bagaimana mendapatkan parameter kinetik dari reaksi dekomposisi

katalitik metana untuk produksi CNT dengan menggunakan katalis Ni-CuAl.

1.3

Tujuan Penelitian a.

Mengembangkan

model

kinetika

transien

dengan

mempertimbangkan deaktivasi katalis dari reaksi dekomposisi katalitik metana dengan menggunakan dasar mekanisme milik Reyyan et.al (2009) b.

Mendapatkan parameter kinetika dari reaksi dekomposisi katalitik metana dengan katalis Ni-Cu-Al

1.4

Ruang Lingkup Masalah a.

Kinetika mikro berdasarkan Mekanisme reaksi didasarkan pada mekanisme yang melibatkan tahap reaksi elementer (Reyyan et.al., 2009).

b.

Mempertimbangkan deaktivasi katalis

c.

Temperatur yang digunakan pada rentang 500 – 750 oC.

d.

Dilakukan pada tekanan 1 atm.

e.

Katalis yang digunakan adalah Ni/Cu/Al dalam bentuk granular.

f.

Metode regresi non-linear menggunakan metode Marquadt

g.

Perangkat lunak yang digunakan untuk penelitian ini adalah MatLab

h.

Data eksperimen yang digunakan adalah jumlah karbon di setiap

waktu i.

Data eksperimen yang digunakan adalah data eksperimen milik Wulan (2011)


BAB 2 PUSTAKA TINJAUAN

Pada Bab II ini akan dilakukan tinjauan pusataka dari apa yang akan dilakukan dan apa yang ada dari penelitian ini. Yang dibahas dalam tinjauan

pustaka ini adalah Carbon Nanotube (CNT), produksi CNT, reaksi dekomposisi metana, perangkat lunak MatLab, dan regresi non linear.

2.1

Nano Karbon Nano karbon didefinisikan sebagai material karbon yang bukan hanya

ukuran partikelnya saja berukuran nanometer akan tetapi struktur dan teksturnya juga berukuran nanometer (Michio, 2004). Carbon Nanotube ditemukan oleh Iijima dengan transmission electron microscopy (TEM) pada tahun 1991 (Iijima, 2002). Sejak saat iu banyak sekali penelitian tentang nanokarbon dilakukan, seperti pada struktur, sifat, dan metode preparasinya. Pada intan, masing-masing atomnya tersusun secara rapat dan secara simetris terkoordinasi memenuhi ruang tiga dimensi. Sedangkan grafit tersusun dari karbon heksagonal yang membentuk lembar atom dua dimensi dimana diantara keduanya tersebut terdapat jarak yang cukup panjang. selain kedua bentuk tersebut, ada beberapa bentuk susunan unsur karbon lainnya yaitu (Petterson,1999): Fullerenes dan Nanotubes adalah lembaran grafit yang terbungkus baik satu lapis maupun lebih yang membentuk bola yang stabil, disebut fullerenes, atau membentuk tabung molekul, disebut nanotube. Nanopores adalah fiber yang sangat kuat yang tersusun dari ikatan nanotubes

1) Nanofibers adalah susunan dari berlapis-lapis grafit yang membentuk fiber berdiamaeter kurang dari 1 mikrometer. 2) Karbon Aktif biasanya berarti granula atau partikel kecil grafit yang dimurnikan. 3) Carbon Fibers yaitu polimer karbon berantai dengan kandungan karbon tinggi Carbon Nanotubes sendiri bisa


5

2.1.1

Mekanisme Pertumbuhan Nano karbon

Beberapa model pertumbuhan dari nano karbon sudah di ajukan yang berasal dari pirolisis hidrokarbon pada permukaan logam. Model pertama menyarankan

mekanisme

4

tahap.

Pada

tahap

pertama,

hidrokarbon

terdekomposisi pada permukaan logam, melepas hydrogen dan karbon yang kemudian larut dalam partikel. Tahap kedua meliputi difusi karbon melalui partikel logan dan pengendapan ada bagian belakang atau yang berlawanan dari

muka atau permukaan untuk membentuk badan filament. Supply karbon ke bagian muka lebih cepat dibandingkan difusi melalui sebagian besar atau logam, sehingga menyebabkan akumulasi karbon pada bagian muka, dimana hal ini harius dihindari untuk menanggulangi penghalangan dari sisi aktif. Hal ini terjadi pada difusi pada permukaan dan karbon membentuk kulit atau lapisan pada badan utama filamen pada tahap ketiga. Pada tahap keempat, penutupan dan deaktivasi katalis dan terminasi pertumbuhan dari nano karbon terjadi. Model lain menyatakan bahwa terdapat mekanisme dimana difusi yang terjadi pada sebagian besar karbon tidak terlalu signifikan dan karbon ditransport mengitari pertikel melalui difusi pada permukaan. Model lain mengemukakan sebuah mekanisme dimana karbon membentuk sebuah hemispherical graphene cap pada partikel katalis dan nano karbon tumbuh dari hal tersebut.

2.2

Cara Memproduksi CNT 2.2.1. Arc Discharge Metode Arc Discharge ini awalnya digunakan untuk memproduksi C60

fullerenes, adalah cara yang paling umum dan mudah untuk memproduksi CNT.

Tetapi teknik ini menghasulkan campuran komponen yang kompleks dan membutuhkan purifikasi lebih lanjut, untuk memisahkan CNT dari jelaga dan residu dari katalis logam. Metode ini menghasilkan CNT dari busur penguapan dari dua tongkat karbon yang diletakkan di masing-masing ujung, sejauh 1m dalam sebuah ruang tertutup yang berisikan gas inert pada tekanan rendah. Aliran listrik langsung sebesar 50 sampai 100 A,akan menghasilkan temperatur pelepasan yang tinggi antara kedua elektroda. Pelepasan ini menguapkan


6

permukaan dari salah satu elektroda karbon dan membentuk tumpukan kecil berbentuk tongkat pada elektroda yang lainnya

2.2.2.

Laser Ablation

Metode ini dilakukan dengan menguapkan tongkat grafit dengan campuran katalis kobalt dan nikel pasa temeratur 1200

O

C dala argon yang mengalir,

dilanjutkan dengan pemanasan dalam vakum pada 1000 OC untuk menyingkirkan

C60 dan fullerenes lain. Penguapan awal dengan laser dan diikuti dengan tembakan ke dua, menguapkan target secara seragam. Penggunaan 2 tembakan meminimalisir tumpukan karbon sebagai jelaga. Tembakan laser kedua akan menghancurkan partikel besar yang terbentuk karena tembakan pertama dan memasukkan mereka kedalam struktur nanotube. Material yang dihasilkan dari metode ini berbentuk tambang, berdiameter 10-20 nm, dan memiliki panjang sekitar 100 μm atau lebih. Tiap tambang terdiri dari nanotube berdinding tunggal. Dengan memvariasikan komposisi katalis dan parameter lainnya, diameter dan besar distribusi dapat di variasikan

2.3

Reaksi Dekomposisi Metana Dekomposisi didefinisikan sebagai salah satu dari reaksi kimia yang

menguraikan atau memutuskan ikatan-ikatan suatu senyawa menjadi unsur-unsur atau senyawa yang lebih sederhana. Salah satu contoh adalah dekompsisi metana yang merupakan hidrokarbon yang paling stabil. Sehingga yang terjadi pada reaksi ini adalah pemutusan ikatan H-C dari metana yang menjadi komponen yang lebih sederhana yaitu hydrogen dan karbon. Pemilihan metana sebagai reaktan

untuk produksi nanokarbon dan hydrogen karena metana merupakan hidrokarbon dengan perbandingan karbon dengan metana yang paling tinggi. Selain itu metana bisa di dapat langsung dari alam tanpa harus diolah terlebih dahulu, sehingga mengurangi biaya produksi. Pada reaksi dekomposisi metana atau Methane Decomposition Reaction (MDR) terjadi reaksi (Chen, 2001) : CH4C+2H2

ΔH298 = +75KJ/mol


(2.1)

7

Reaksi perngkahan metana dapat dibedakan menjadi 2 jenis yaitu Cracking) dan perengkahan secara tidak perengkahan langsung (Direct Methane

langsung (Indirect Methane Cracking). Reaksi dekomposisi metana tergolong reaksi perengkahan langsung dan saat ini sudah diketahui terdapat 2 golongan untuk reaksi perengkahan ini yaitu: 

Dekomposisi metana secara thermal (Thermal Cracking) yang akan menghasilkan karbon hitan dan hydrogen,



Dekomosisi katalitik metana yang menggunakan katalis berdasar logam transisi (Fe, Co, dan Ni) yang akan menghasilkan hydrogen dan nanokarbon.

Analisis termodinamika dari reaksi dekomposisi metana menyatakan bahwa nilai untuk energi bebas Gibbs

dan energi penguraian metana

pada

suhu 198 K, masing-masing sebesar 50,8 kJ/mol dan 75 kJ/mol (Song L, 2005). Nilai

yang positif ini menunjukkan bahwa reaksi tidak dapat berjalan dengan

spontan. Meskipun reaksi dapat berjalan, konversi yang dihasilkan reaksi tersebut tidak akan maksimal. Sedangkan nilai

yang positif menandakan reaksi

bersifat endotermis. Hal ini akan mempengaruhi konversi yang akan meningkat seiring dengan meningkatnya suhu reaksi sehingga reaksi ini harus dilakukan pada temperatur sangat tinggi. Temperatur reaksi memegang peranan penting dalam tinjauan ekonomi suatu proses industri. Proses yang dilakukan pada temperatur yang lebih rendah akan lebih menguntungkan jika kita bandingkan dengan proses bertemperatur tinggi. Hal ini disebabkan energi yang dibutuhkan lebih sedikit. Oleh karena itu, untuk menurunkan temperatur reaksi perlu ditambahkan katalis yang dapat

menurunkan energi aktivasi.

2.3.1 Katalis Pada Reaksi Dekomposisi Metana Masalah yang sering timbul pada reaksi dekomposisi katalitik metana adalah deaktovasi katalis yang terjadi akibat pembentukan karbon pada permukaan katalis (Muradov, 2001). Karbon yang dihasilkan dari perengkahan metana akan membentuk deposit karbon yang dapat menutupi inti aktif katalis sehingga terjadi deaktivasi katalis dan menurunnya lifetime katalis. Katalis yang


8

digunakan saat ini sangat mudah terdeaktivasi pada suhu tinggi. Untuk mengatasi maka kondisi operasi dijaga pada suhu hal ini dan menjaga stabilitas katalis,

rendah (<873 K) dengan memperkaya aliran gas metana, dan senantiasa mengeluarkan hydrogen dari hasil produk sehingga konversi metana menjadi lebih

tinggi.

Penggunaan katalis berbasis nikel (Ni) pada penelitian dekompsisi metana karena nikel merupakan logam yang paling aktif untuk reaksi dekomposisi

katalitik metana disbanding logam-logam lainnya. Tetapi Ni mudah terdeaktivasi dan sintering, sehingga menyebabkan diameter partikel Ni membesar. Hal ini akan mempengaruhi kualitas karbon nanotube. Penambahan promoter baik promoter structural maupun tekstural dapat mencegah terjadinya sintering. Promoter structural yang biasa digunakan adalah Cu (Chen, 2001; Chen, 2004) karena dapat menurunkan titik leleh katalis dan meningkatkan aktivitas katalis. Titik leleh yang rendah diperlukan karena kondisi katalis dalam keadaan liquid adalah syarat terbentuknya karbon nanotube (Chen,

2001). Cu mempunyai

afinitas tinggi terhadap struktur grafit dan menghambat pertumbuhan deposit karbon pada permukaan Ni. Sedangkan penambahan promoter tekstural akan memperbaiki tekstur dari katalis sehingga baik untuk digunakan. Pada penelitian ini digunakan katlis Ni-Cu-Al, dimana Cu digunakan sebagai promotor untuk mencegah sintering dan alumina digunakan sebagai yang menguatkan ikatan antar logam sehingga dapat menjaga diameter Ni tetap kecil.

2.3.2 Kinetika Pada Reaksi Dekomposisi Metana Di dalam studi kinetika dikenal adanya teknik penentuan kinetika keadaan

tunak (steady state kinetik) dan kinetika dalam keadaan dinamik (transient kinetik). Untuk mendapatkan persamaan laju reaksi katalitik dapat dilakukan dengan dua metode, yaitu: 

Metode analisa kinetika makro



Metode analisa kinetika mikro

Pada metode analisa kinetika makro, penyusunan persamaan laju reaksi menggunakan prinsip ‘hukum pangkat’ yang kemudian dievaluasi dengan data percobaan. Metode ini etbagi menjadi dia yaitu hukum pangkat sederhana (smiple


9

power law) dan hukum pangkat kompleks (complex power law). Pengolahan data kinetika dengan metode hukum pangkat atau pendekatan empiris dimaksudkan

untuk mendapatkan persamaan laju reaksi atau menentukan harga konstanta laju reaksi dan orde reaksi terhadap konsentrasi reaktan atau tekanan parsial reaktan.

Pada reaksi katalitik, orde reaksi biasanya berupa bilangan pecahan dan pesarnya dipengaruhi oleh konsentrasi.

Persamaan yang diperoleh dari metode ini dalam perancangan reactor hanya

terbatas pada konsidi pecobnaan penentuan persamaan tersebut karena metode ini berdasarkan persamaan empiris. Oleh karena itu, pada perkembangan selanjutnya sering persamaan laju reaksi katalitik ditentikan berdasarkan mekanisme reaksinya. Persamaan kinetika mekanistik ini diharapkan dapat berlaku pada kondisi yang lebih lebar dibandingkan dengan persamaan kinetika dari metode hukum pangkat (Slamet, 1999). Sedangkan pada metode analisa kinetika mikro, penyusunan persamaan lajunya bertumpu pada mekanisme reaksi yang terjadi. Analisa kinetika mikro adalah pengujian reaksi-reaksi katalitik berdasarkan pada reaksi-reaksi kimia elementer yang terjadi pada permukaan katalitik dan hubungannya satu sama lain dengan permukaan katalitik selama siklus katalitik terjadi. Strategi yang dipakai pada metode analisa kinetika mikro adalah dengan cara memadukan antara data eksperimen yang ada dengan konsep-konsep teoritis dari system katalitik yang terkait. Salah satu contoh analisa kinetika mikro yang sering digunakan adaah metode

isotherm

adsopsi,

dibutuhkan

hubungan

kesetimbangan

antara

konsentrasi/tekanan parsial reaktan dan produk di fasa fluida dengan banyaknya reaktan dan produk yang teradsorpsi di permukaan katalis. Hubungan kesetimbangan pada temperatur konstan ini disebut isotherm adsorpsi.

Di dalam adsorpsi kimiawi, ada beberapa jenis isotherm adsorpsi, diantaranya adalah isotherm Langmuir. Isoterm inibiasa di sebut dengan isotherm Langmuir-Hinshelwood-Hougen-Watson. Isotherm ini diturunkan berdasarkan asumsi bahwa keaktifan inti katalis seragam dan adsorpsi yang terjadi adalah reaksi reversible elementer antara adsorbat dengan inti aktif katalis. Tahap penentu laju reaksinya adalah reaksi permukaan antara reaktan-reaktan

yang

teradsorpsi yang berlangsung sangat lambat, sehingga mengendalikan laju reaksi.


10

Mekanisme Langmuir ini sering dipakai karena bentuknya sederhana dan hasil yang memadai (Slamet, 1996).

Jenis isotherm adsorpsi yang lain adalah Eley-Rideal.Mekanisme ini terjadi pada kondisi yang sangat vakum, dengan asumsi bahwa salah satu reaktar

teradsorpsi secara kimiawi pada inti aktif katalis, sementara reaktan kedua berada pada fasa gas atau fasa teradsorpsi fisik leman (Hinink, 1995; Slamet 1996). 2.4

Kondisi Operasi

Kondisi operasi yang baik akan mengahasilkan produk akhir dari dekomposisi metana yang sesuai harapan. Untuk dekomposisi metana, terdapat dua kondisi operasi yang sangat mempengaruhi hasil akhir, yaitu temperatur dan tekanan.

2.4.1. Temperatur Operasi Dekomposisi katalitik matanamembutuhkan temperatur yang cukup tinggi. Semakin tinggi temperatur operasi, maka gas metana akan semakin cepat tedekomposisi. Hal ini akan mempengaruhi reaksi ini baik dari segi kualitas maupun kuantitas. Namun di sisi lain, suhu yang tinggi akan mengakibatkan katalis mudah mengalami sintering dan karbon yang dihasilkan akan semakin banyak sehingga deaktivasi katlis juga semakin cepat terjadi. Meskipun memiliki dampak negatif, penggunaan temperatur tinggi menyebabkan jumlah karbon yang terbentuk akan gmnejadi lebih banyak. Hasil uji produk membuktikan bahwa CNT hanya akan terbentuk pada temperatur lebih tinggi dari pada temperatur terbentuknya karbon nanofiber pada umumnya.

Karbon yang terbentuk pada temperatur tinggi lebih teratur dengan ketebalan seragam (Resketenko, 2003), tetapi penambahan temperatur juga akan menyebabkan diameter dan ketebalan dinding dari nanotube yang telah terbentk semakin mengecil (Sinnots, 1999), sedangkan ng diharapkan adalah nanotube karbon dengan dinding yang tebal. Semuahal ini membuat tidak ada temperatur tertentu yang mutlak untuk reaksi ini. Temperatur optimal untuk reaksi dekompsisi metana agar terbentuk CNT yang diharapkan berada pada temperatur 650-750oC.


11

2.4.2. Tekanan Operasi

Pada reaksi dekomposisi metana, tekanan operasinya adalah tekanan rendah atau tekanan atmosferik. Hal ini dikarenakan reaksi dekomposisi metana adalah

reaksi endotermis dimana penambahan tekanan malah akan menggeser kesetimbangan reaksi ke arah kiri (menambah reaktan). 2.5

Matrix Laboratory (MatLab) Matrix Laboratory (MatLab) adalah sebuah bahasa programing generasi

keempat yang memberikan pengguna sebuah lingkungan untuk perhitungan numeris. MatLab mengintegrasikan komputasi perhitungan, visualisasi dan programing dalam lingkungan yang mudah untuk digunakan dimana masalah dan penyelesaian diekspresikan dalam notasi matematis. MatLab juga merupakan system interaktif yang data dasarnya tidak membutuhkan dimensi. Hal ini memberikan kesempatan pada pengguna untuk menyelesaikan berbagai perhitungan teknis terutama yang memiliki matriks dan vector. Sistem dari MatLab terdiri dari bagian-bagian utama yaitu: 1. Desktop Tools and Development Environment 2. Mathematical Function Library 3. The Language 4. Graphics 5. External Interfaces Walaupun matlab bertujuan utama untuk perhitungan numeris, dengan adanya tambahan dari simulink, maka matlab mampu menambahkan simulasi

daerah secara grafis, dan desain berbasis dasar model. MatLab digunakan oleh berbagai jenis pengguna berlatar belakang berbeda-beda seperti insinyur, peneliti, dan ekonom. Penelitian ini menggunakan MatLab sebagai perangkat untuk Validasi data kinetik dan proses curve fitting.


12

Gambar 2. 1 Logo MatLab (www.Matlab.com)

Gambar 2. 2 Tampilan Antarmuka MatLab (www.Matlab.com)

2.6

Regresi Non Linear Analisa Regresi atau regresi adalah aplikasi dari metode matematis dan

statistik untuk menganalisa data eksperimen dan mem-fit model matematis dari


13

data dengan mengestimasi parameter-parameter yang tidak diketahui dari model. menemani analisa regresi, memberikan Serangkaian tes statistik yang biasanya

identifikasi model, verifikasi model dan efisiensi desain dari program eksperimen. sebuah penyelesaian, akan ada lebih dari Pada regresi non-linear, dalam

satu parameter yang tidak diketahui dan haris dihitung dengan mem-fit model dengan data eksperimen. Dalam regresi non linear, mungkin saja dilakukan fitting untuk beberapa variable terikat secara bersamaan dengan menggunakan beberapa

persamaan terhadap data. (Constantinides, 2000) Pada regresi non linear, terdapat beberapa teknik untuk meminimalisir jumlah dari kuadrat residu atau nilai sisa. Beberapa metode yang ada akan membantu kita untuk memfit model yang mengandung beberapa variable, terhadap data eksperimen untuk mendapatkan nilai dari parameter model yang memiliki nilai jumlah dari kuadrat residual keseluruhan yang kecil. Ditambah lagi dengan adanya analisa statistik dari hasil regresi, akan mambantu kita untuk : 

Menentukan apakah model memuaskan atau tidak dalam eksperimen



Membedakan model yang bersaing



Menentukan ketepatan dari estimasi parameter dengan membangun daerah terpercaya pada ruang parameter



Mengukur korelasi antara parameter dengan meninjau matrix korelasi koefisien



Melakukan tes untuk memverifikasi data eksperimen yang data dari populasi eksperimen yang sama



Melakukan tes untuk memverifikasi apakah residu antara data dan model terdistribusi secara acak.

2.6.1.

Metode Marquadt

Marquadt telah mengembangkan sebuah teknik interpolasi antara metode Gauss-Newton dan metode steepest descent. Interpolasi ini didapat dengan menambahkan matriks diagonal (λI) ke matriks (J’J): (2.2) Nilai dari λ dipilih dari setiap iterasi sehingga parameter vektor terkoreksi akan menghasilak nilai jumlah dari kuadrat yang lebih kecil dari iterasi yang


14

bersangkutan. Dapat mudah dilihat bahwa nilai λ yang kecil dibanding elemen matriks (J’J), metode Marquadt akan mendekati metode Gauss-Newton. Sedangkan ketika nilai λ sangat besar, maka metode ini akan mendekati metode steepest descent, dengan pengecualian faktor skala yang tidak mempengaruhi arah dari parameter koreksi vektor tetapi member nilai step yang kecil. Berdasarkan Marquadt, nilai Φ yang kecil diinginkan untuk di daerah sekitar fungsi yang dilinearkan untuk memberikan representasi dari fungsi non linear

yang memenuhi syarat. Oleh sebab itu metode untuk menentukan λ harus memberikan nilai λ kecil ketika metode Gauss-Newton akan konvergen secara efisien dan nilai λ ketika metode steepest descent dibutuhkan. Metode Marquadt dapat diaplikasikan pada metode Newton. Dalam kasus ini matriks diagonal λI ditambahkan pada matriks Hessian: (2.3) Metode Marquadt memiliki langkah-langkah sebagai berikut: 1. Asumsi tebakan awal parameter vektor b. 2. Menentukan nilai yang besar pada λ, misalkan 1000. Hal ini menunjukkan bahwa iterasi pertama akan didominasi oleh metode steepest descent dan memastukan metode berjalan kenilai jumlah kuadrat residu yang rendah. 3. Mengevaluasi

matriks

jacobian

J

dari

persamaan

model.

Juga

mengevaluasi matriks Hessian H jika menggunakan metode Newton. 4. Mencari vektor koreksi Δb. 5. Evaluasi estimasi parameter vektor baru dengan persamaan: q

(2.4)

6. Menghitung nilai baru dari Φ. Jika Φ(m+1) < Φ(m), kurangi nilai λ. Jika Φ(m+1) > Φ(m) biarkan parameter yang lama [

] dan tambahkan

nilai λ. 7. Ulangi langkah 3 sampai 6 hingga salah satu atau kedua kondisi ini terpenuhi: a. Φtidak berubah secara signifikan. b.

Δb menjadi sangat kecil (sesuai dengan nilai lriteria konvergensi).


BAB 3 METODE PENELITIAN

3.1

Diagram Alir Penelitian

Menentukan Mekanisme Reaksi

Membuat Persamaan laju reaksi dari masing-masing spesies yang terlibat

Regresi non linear dengan metode Marquadt pada MatLab

Data eksperimen

Konstanta laju reaksi (k1-k11) pada masing-masing suhu yang terdapat pada data percobaan

Plot ln k terhadap 1/T untuk setiap nilai k (k1-k11)

Nilai pra eksponensial energi aktivasi untuk masing-masing nilai k (k1-k11)

Gambar 3. 1 Diagram Alir Penelitian


16

3.2

Rincian Kegiatan Penelitian

a. Menentukan mekanisme reaksi Tahap ini diawali dengan menentukan reaksi yang terjadi pada metana. Pada penelitian ini, terdapat tiga percobaan dekomposisi katalitik

tahap utama dalam reaksinya yaitu Reaksi di permukaan katalis (Surface Reaction), pembentukan nano karbon (Nano Carbon Forming), dan deaktivasi katalis (Catalyst Deactivation). Pada tahap ini, juga di turunkan

persamaan laju yang terjadi pada setiap tahap dalam percobaan. Persamaan laju diturunkan dari reaksi-reaksi yang terjadi pada setiap tahapannya.

b. Pemodelan kinetik Permodelan dari tahap ini dimaksudkan untuk merumuskan beberapa model persamaanlaju reaksi dekomposisi katalitik metana menjadi hydrogen dan nano karbon. Penyusunan persamana kinetika ini dilakukan berdasarkan metode analisa kinetika yang bertumpu pada metode analisa kinetika mikro yang dijbarkan menurut mekanisme reaksi katalitik yang mungkin terjadi: 

Reaksi dekomposisi metana secara keseluruhan (Chen et al, 2001) (3.1)



Tahapan

reaksi

dekomposisi

katalitik

metana

berdasarkan

mekanisme (Reyyan, 2009):

I. Reaksi pada permukaan (Surface Reaction) 1. CH4 + 2S

CH3.S + H.S

2. CH3.S + S

CH2.S + H.S

3. CH2.S + S

CH.S + H.S

4. CH.S + S 5. 2H.S

II.

C.S + H.S

H2 + 2S

Pembentukan nano karbon 6. C.S 7. CNi,F

CNi,F + S CNi,R


17

8. CNi,R  Cf

III.

Deaktivasi katalis

9. nC.S  nCP

Penjabaran dari penurunan model persamaan kinetika reaksi dekomposisi katalitik metana dapat dilihat pada lampiran A. Tabel 3.1 adalah ringkasan model persamaan laju reaksi dekomposisi katalitik metana menjadi hidrogen dan nano

karbon berdasarkan mekanisme reaksi: Tabel 3. 1 Tabel Model Persamaan Laju Reaksi Masing-Masing spesies zat

Nama spesies zat

Model Persamaan laju reaksi

Persamaan

CH4

RCH4 =

(3.2)

CH3.S

RCH3.S=

(3.3)

CH2.S

CH.S C.S

H.S

(3.4)

RCH2.S =

(3.5)

RCH.S =

(3.6)

RC.S =

RH.S =

(3.7)

H2

(3.8) (3.9)

(3.10)

CP

(3.11)


18

Seluruh persamaan yang terdapat pada tabel 3.1, akan dimasukkan kedalam program yang sudah di buat terlebih dahulu di Matlab.

Persamaan-persamaan inilah yang akan di cari solusinya setelah parameter-parameter (k1 – k11) di tebak nilai awalnya terlebih dahulu.

c. Regresi non linear dengan menggunakan metode Marquadt Pada tahap ini, dilakukan regresi untuk menentukan parameter

parameter yang masih belum diketahui. Tahap ini diawali dengan melakukan regresi non linear terhadap data yang sudah di dapat dengan menggunakan metode Marquadt. Kemudian dilakukan trial and error untuk mendapatkan nilai parameter-parameter kinetik yang paling kecil nilai error nya atau dalam grafik akan menghasilkan bentuk yang paling sesuai dengan persebaran data yang ada. Tahap ini dilakukan pada perangkat lunak Matrix Laboratory. Fitting model dilakukan dengan cara menyelesaikan persamaan-persamaan model secara bersamaan. Kemudian hasil utama dari model yang berupa jumlah karbon yang terbentuk selama rentang waktu tertentu. Hasil simulasi kemudian akan dicocokan dengan data eksperimen yang ada. Berikut

ini

adalah

algoritma

dari

regresi

yang

dilakukan:


19

Memasukkan persamaan-

Mengevaluasi Matriks

jacobian model

persamaan model yang

sudah di buat

Mencari Vektor Koreksi Δb Memasukkan data eksperimen (t, gC) Tidak

Φ(m+1+)>Φ(m) Asumsi tebakan awal vektor b (nilai k)

Tidak

Ya

Kurangi nilai λ

Tambahkan nilai λ

Menentukan nilai awal λ (Damping Factor) (>1000) ΔΦ < 10% Memastikan iterasi pertama berjalan

Ya Selesai Gambar 3. 2 Algoritma Program Regresi Non Linear


20

d. Penentuan nilai parameter kinetik Pada tahap ini, setelah didapatkan nilai konstanta laju reaksi dari masing-masing suhu, nilai parameter kinetik akan didapatkan dengan memplot nilai ln k (konstanta laju reaksi) yang didapat terhadap 1/T,

sehingga nanti akan didapatkan nilai pra eksponensial dan nilai dari energi

aktivasi dari masing-masing konstanta laju reaksi.


BAB 4

HASIL DAN PEMBAHASAN

Hasil dari penelitian ini akan ditampilkan, dianalisis dengan jelas. Hasil

penelitian ini berupa nilai parameter kinetik dari reaksi dekomposisi katalitik

metana yang dilakukan dalam eksperimen. Mekanisme reaksi kimia akan menghasilkan parameter kinetik setelah hasil dari model dibandingkan dengan

data eksperimen. Dengan mensimulasikan dengan variasi suhu, akan didapatkan juga nilai pra eksponensial dan energi aktivasi dari reaksi yang dimodelkan pada penelitian ini.

4.1

Mekanisme Reaksi Mekanisme reaksi yang digunakan pada penelitian ini didasarkan pada

mekanisme reaksi milik Reyyan (2009), dimana pada model reaksi dekomposisi katalitk metana ini, reaksi dibagi menjadi tiga bagian utama yaitu reaksi di permukaan katalis yaitu reaksi pemutusan – pemutusan atom hidrogen dari metana, pembentukan nano karbon yang terdiri dari larutnya atom karbon dalam katalis, difusi karbon dari bagian depan katalis ke bagian belakang katalis, serta pembentukan nano karbon itu sendiri, dan reaksi yang terakhir adalah enkapsulasi katalis, dimana sebagian atom karbon akan membentuk sebuah lapisan karbon yang akan menutupi katalis. Berikut ini adalah ilustrasi dari keseluruhan reaksi yang terjadi: H2 CH4

Pembentukan Nano Karbon(Tahap

Tahap 5

6 sampai 8) Katalis Nano CH4

C.S

Karbon

Reaksi Permukaan (tahap 1 Tahap 9

sampai tahap 4) Enkapsulasi

Gambar 4. 1 Ilustrasi Keseluruhan Reaksi

21 Penentuan parameter..., Marcho Rizal, FT UI, 2012

22

Pada reaksi ini, hidrogen yang sudah terlepas dari karbon dapat berikatan kembali dengan karbon pada permukaan katalis sehingga dapat terjadi proses desorpsi dari karbon terhadap katalis. Proses hidrogenasi ini akan terjadi pada setiap intermediate selama di reaksi permukaan katalis.. Gambar 4.2 adalah skema

mendetail mekanisme yang digunakan:

CH4 H.S

H2

(1) CH3 . S

H.S

(2)

(5)

CH2 . S

H.S

(3) CH . S H.S (4) Lapisan karbon

(9)

C.S (6)

enkapsulasi

C NiF (xF)

(7) C NiR (xR)

(8)

(8) Nano Karbon

Gambar 4. 2 Skema Mekanisme Detail Dari Reaksi yang Digunakan Dalam Model

Universitas Indonesia Penentuan parameter..., Marcho Rizal, FT UI, 2012

23

4.1.1 Reaksi di permukaan

Pada model, reaksi di permukaan terdiri dari reaksi elementer yang terjadi di permukaan nikel atau katalis seperti yang tersaji pada Gambar 4.3 hingga 4.7 yaitu:

CH4 + 2S

CH3.S + H.S

(4.1)

CH4

CH3 H

Ni

Ni

Gambar 4. 3 Reaksi Pemutusan Hidrogen dari Metana Tahap 1 Pada tahap ini, metana yang berada di permukaan katalis akan mengalami dehidrogenasi, dimana salah satu atom hidrogen yang berikatan pada karbon akan diputus atom hidrogen tersebut akan menempel pada sisi aktif katalis dan intermediate CH3 (CH3.S) terbentuk.

CH3.S + S

CH2.S + H.S

(4.2)

CH3

CH2

H Ni

H Ni

H

Gambar 4. 4 Reaksi Pemutusan Hidrogen dari Metana Tahap 2


24

Pada tahap ini, intermediate CH3 (CH3.S) yang berada di permukaan katalis akan mengalami dehidrogenasi, dan menghasilkan intermediate CH2 ( CH2.S).

CH2.S + S

CH.S + H.S

(4.3)

CH2

CH

H

H

Ni

Ni

H

H

H

Gambar 4. 5 Reaksi Pemutusan Hidrogen dari Metana Tahap 3 Pada tahap ini, intermediate CH2 ( CH2.S) yang berada di permukaan katalis akan mengalami dehidrogenasi, dan menghasilkan intermediate CH (CH.S).

CH.S + S

C.S + H.S

(4.4)

CH

C H

H H

Ni

Ni

H

H

H

H

Gambar 4. 6 Reaksi Pemutusan Hidrogen dari Metana Tahap 4 Pada tahap ini, intermediate CH (CH.S) yang berada di permukaan katalis akan mengalami dehidrogenasi, dan menghasilkan intermediate C (C.S).

2H.S

H2 + 2S

(4.5)


25

H

H

H2

Ni

Ni

H

H

H

2

Gambar 4. 7 Pembentukan Gas Hidrogen dari Atom H yang Terikat Intiyang Aktif Pada tahap ini, seluruh atomdengan hidrogen sudah terlepas dari karbon akan berikatan satu sama lain, dan membentuk gas hidrogen 4.1.2 Pembentukan Nano Karbon Reaksi pada permukaan katalis akan menghasilkan atom karbon yang teradsorb (C.S) dimana spesies ini akan terlarut kedalam partikel nikel pada fasa gas dan akan berdifusi dari bagian depan ke bagian belakang dari partikel nikel, kemudian membentuk nano karbon. Pada model ini, tahapan yang meliputi bagian pembentukan nano karbon adalah sebagai berikut:

Segregasi C.S

CNi,F + S

(4.6)

C

C

Gambar 4. 8 Segregasi Karbon yang Terikat pada Katalis ke Bagian Depan Katalis


26

Pada bagian depan nikel, sebuah laspisan dengan konsentrasi karbon yang tinggi akan terbentuk karena sifat segregasi karbon dalam nikel. Permukaan yang kaya akan karbon, dan konsentrasi karbon di permukaan akan berkurang. Konsentrasi dari karbon yang terlarut dapat digambarkan dengan isotherm

segregasi. Untuk atom karbon kompleks dengan sisi aktif, “C.S”, dapat

digambarkan sebagai atom karbon yang teradsorb. Terdapat persamaan Langmuir untuk segregasi digunakan untuk menentukan nilai konsentrasi karbon pada

lapisan segregasi atau konsentrasi karbon pada bagian depan katalis. (4.7) (4.8) Dimana Dimana θC.S adalah fraksi dari katalis yang tertutup karbon dan = -45,217-14.24 T (j/mol) (Snoeck et al, 1997).

4.1.4Difusi karbon melewati nikel CNi,F

CNi,R

(4.9)

C

Front

2/3 DNi

C

Rear

Gambar 4. 9 Difusi Karbon Melalui Nikel Perbedaan dari potensi kimia dari material karbon akan menghasilkan perbedaan dalam tingkat kelarutan pada bagian depan dan belakang katalis (nikel). Hal ini akan menghasilkan gaya dorong utuk difusi karbon melewati nikel. Laju dari difusi dapat ditulis dengan hukum Fick: (4.10)


27

Dengan jarak rata-rata yang dilewati oleh karbon diasumsikan sebesar 2/3 diameter nikel,

= difusivitas karbon dalam nikel,

karbon terlarut di nikel pada bagian depan Ni,

= konsentrasi dari

= konsentrasi dari karbon

terlarut di nikel pada bagian belakang Ni dan ANi = Luas Spesifik dari partikel nikel

4.1.5 Pembentukan nano karbon CNI,R  Cf

(4.11)

Ni

Nano Karbon

Gambar 4. 10 Pembentukan Nano Karbon Reaksi pembentukan nano karbon terjadi terlalu cepat sehingga diasumsikan bahwa seluruh atom karbon yang mencapai bagian belakang dari katalis membentuk nano karbon, Sehingga tahap ini tidak dimasukkan dalam model.

4.1.6 Deaktivasi Katalis Pada bagian ini, atom karbon yang sudah teradsorb, tidak seluruhnya terlarut dalam katalis. Sebagian dari karbon yang teradsorb akan membentuk sebuah lapisan sehingga mengenkapsulasi permukaan aktif dari katalis. Pada akhirnya, sisi aktif yang tertutup akan meghalangi interaksi antarmuka dari reaktan berfasa gas dalam hal ini metana dan intermediet-intermediet yang ada terhadap sisi aktif


28

dari nikel atau katalis. Sehingga diasumsikan bahwa perubahan pada penutupan permukaan karena enkapsulasi adalah penyebab utama dari deaktivasi. Pembentukan karbon yang mengenkapsulasi katalis dapat ditulis sebagai berikut: nC.S  nCp

(4.12)

Cenkapsulasi

Gambar 4. 11 Deaktivasi Katalis karena Enkapsulasi Katalis Berdasarkan Chen et al. nilai n yang terbaik adalah tiga (3) dimana nilai ini akan memberikan kecocokan paling tinggi terhadap data yang ada.

4.2

Model Kinetika dan Fitting Model terhadap Data Eksperimen Model kinetika dari dekomposisi katalitik metana untuk menghasilkan

nano karbon terdiri dari tiga bagian besar yaitu reaksi di permukaan katalis, pembentukan nano karbon, dan deaktifasi dari katalis (enkapsulasi). Model kinetika reaksi dekomposisi katalitik metana terdiri dari 7 persamaan diferensial serta 2 persamaan aljabar linear. Dalam model kinetic, persamaan untuk nilai segregasi didapatkan dari data yang diberikan oleh Chen et al (2001). Yaitu:

nilai dari simulasi akan merepresentasikan nilai dari eksperimen (terutama nilai dari jumlah karbon yang terbentuk). Model kinetika ini diselesaikan tanpa menentukan tahap penentu laju reaksi. Model ini juga menggunakan asumsi bahwa seluruh atom karbon yang berdifusi melalui katalis akan membentuk nano karbon hingga katalis terdeaktivasi. (Reyyan, 2009)


29

Model akan disimulasikan dengan menggunakan parameter kinetik yang belum diketahui dan kemudian ditebak dengan menggunakan dasar reaksi maju dan reaksi balik dari dekomposisi metana berdasarkan Snoeck et al (1997) yaitu 8:1 yang merupakan rasio antara reaksi dekomposisi metana dengan reaksi

hidrogenasi karbon sehingga menghasilkan metana dengan menggunakan katalis

berbasis nikel. Kemudian hasil simulasi akan dibandingkan dengan data. Berikut

Jumlah carbon (gram)

adalah perbandingan antara model dengan data yang ada: 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 -0.05 0

datal model

1000

2000

3000

4000

5000

Waktu (detik)

Gambar 4. 12 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=750 oC

Jumlah Carbon (gr)

0.3 0.25 0.2

0.15

data

0.1

model

0.05 0 0.00

1000.00

2000.00

3000.00

4000.00

Waktu (detik)



30

0.35

Jumlah carbon (gr)

0.3

0.25 0.2

0.15

data model

0.1

0.05 0 -0.05

0

1000

2000

3000

4000

5000

Waktu (detik)


0.3

Jumlah Carbon (s)

0.25 0.2 0.15 data 0.1

model

0.05 0 0 -0.05

1000

2000

3000

4000

5000

Waktu (detik)



31

0.35

Jumlah carbon (gr)

0.3 0.25

0.2

0.15

data model

0.1

0.05 0 -0.05

0

1000

2000

3000

4000

Waktu (detik)


0.18 0.16

Jumlah carbon

0.14 0.12 0.1 0.08

data

0.06

model

0.04 0.02 0 0

1000

2000

3000

4000

Waktu (detik)

Gambar 4. 17 Fitting Hasil Simulasi Terhadap Data Untuk Karbon yang Terbentuk Terhadap Waktu pada T=500 oC Simulasi dilakukan dengan menggunakan kriteria konvergensi (Δb) sebesar 10-10 dengan nilai kepercayaan sebesar 95%. Dapat dilihat pada gambar 4.12 hingga 4.17 bahwa secara umum model dan data sudah cukup sesuai dan dapat merepresentasikan apa yang terjadi dalam reaksi tersebut. Asumsi yang


32

digunakan untuk membuat model kinetik akan mempengaruhi efisiensi dari model, seperti asumsi bahwa konsentrasi karbon di bagian belakang katalis diasumsikan nol (0) karena atom karbon yang membentuk nano karbon terlalu cepat dalam penggabungannya. Selain itu penyebaran ukuran partikel nikel dan

pembentukan nano karbon dianggap seragam (Reyyan, 2009).

4.3

Penentuan Parameter Kinetik (Nilai Pra Eksponensial Konstanta

Laju Reaksi dan Energi Aktivasi) Parameter kinetik dari masing-masing persamaan dalam model didapatkan dengan mencocokan data dengan hasil simulasi dan data eksperimen pada suhu 500oC hingga 750oC. Pada simulasi, akan didapatkan nilai konstanta laju reaksi (k1 hingga k11) untuk suhu yang ada dengan interval 50oC. Berikut ini adalah nilai dari masing-masing konstanta laju reaksi pada masing-masing suhu: Tabel 4. 1 Tabel nilai Konstanta Laju Reaksi untuk nilai k1 hingga k11 Suhu (oC) 750 700 650 600 550 500 Suhu (oC)

750 700 650 600 550 500

k1 8,47E-06 6,12E-06 4,34E-06 3,21E-06 2,00E-06 1,00E-06

k2 1,05E-06 8,19E-07 5,13E-07 4,17E-07 3,04E-07 1,37E-07

k3 8,75E-05 6,33E-05 4,12E-05 3,28E-05 2,46E-05 1,09E-05

k4 1,12E-05 7,17E-06 5,32E-06 4,12E-06 2,65E-06 1,26E-06

k5 8,43E-04 6,24E-04 4,13E-04 3,21E-04 2,31E-04 1,08E-04

k6 k7 k8 k9 k10 k11 1,02E-04 8,12E-03 1,00E-03 8,08E-02 1,00E-02 2,00E-02 7,55E-05 6,23E-03 7,40E-04 6,19E-02 7,50E-03 9,00E-03 5,00E-05 4,42E-03 5,12E-04 4,09E-02 5,00E-03 7,00E-03 4,54E-05 3,20E-03 4,08E-04 3,31E-02 4,00E-03 3,00E-03 2,66E-05 2,00E-03 2,82E-04 2,12E-02 2,50E-03 8,00E-04 1,28E-05 1,20E-03 1,15E-04 1,08E-02 1,25E-03 4,00E-04

Deviasi rata-rata dari keseluruhan nilai K yang didapatkan dari penelitian kurang dari 5%. Berikut ini adalah nilai deviasi dari masing-masing konstanta laju reaksi untuk setiap suhu yang ada:


Tabel 4. 2 Tabel nilai K dan Deviasi Standar pada suhu 750oC

k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11

Nilai k 8.47E-06 1.05E-06 8.75E-05 1.12E-05 8.43E-04 1.02E-04 8.12E-03 1.00E-03 8.08E-02 1.00E-02 2.00E-02

Deviasi 2.54E-07 5.25E-08 3.50E-06 3.58E-07 4.38E-05 2.24E-06 2.92E-04 4.30E-05 4.04E-03 2.70E-04 0.0006

Tabel 4. 5 Tabel nilai K dan Deviasi

Standar pada suhu 600oC

k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11


nilai k 3.21E-06 4.17E-07 3.28E-05 4.12E-06 3.21E-04 4.54E-05 3.20E-03 4.08E-04 3.31E-02 4.00E-03 3.00E-03

deviasi 9.63E-08 2.09E-08 1.31E-06 1.32E-07 1.67E-05 9.99E-07 1.15E-04 1.75E-05 1.66E-03 1.08E-04 0.00009


o

Standar pada suhu 700 C




k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11


k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11



o


Standar pada suhu 650 C

k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11


deviasi 6E-08 1.52E-08 9.84E-07 8.48E-08 1.20E-05 5.85E-07 7.20E-05 1.21E-05 1.06E-03 6.75E-05 0.000024


k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11



deviasi 3E-08 6.85E-09 4.36E-07 4.03E-08 5.62E-06 2.82E-07 4.32E-05 4.95E-06 5.40E-04 3.38E-05 0.000012

Dengan persamaan Arrhenius , dapat diperoleh faktor tumbukan (nilai pra eksponensial), A (molekul/sisi aktif. Detik) dan nilai energi aktivasi, Ea (kJ/mol). Untuk nilai A dan Ea diperoleh dari hasil plot nilai ln (k) dengan (1/T). Hasil pengujian antara ln (k) dengan (1/T) dilakukan pada setiap

nilai k agar menghasilkan energi aktivasi untuk masing-masing reaksi pemutusan

serta reaksi deaktivasi yang terjadi. Tabel 4.1 menunjukkan hasil pengujian untuk nilai pra eksponensial masing-masing konstanta laju reaksi dan energi aktivasi yang diperoleh:

Tabel 4. 8 Tabel Nilai Pra Eksponensial dan Energi Aktivasi dari Masing-Masing k

k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11

A (molekul/sisi aktif.detik) 0,005253 0,000426 0,035295 0,006054 0,353101 0,046468 3,092562 0,540101 34,39805 5,139727 0,000219535

Energi Aktivasi (J/mol) 54.398,5 50.781,91 51.047,96 53.798,44 51.230,87 51.746,34 50.274,76 53.068,26 51.180,98 52.785.59 104.673,3

. Untuk hasil plot antara ln k terhadap 1/T untuk setiap nilai k dapat dilihat pada bagian lampiran D. Dapat dilihat pada tabel 4.1, bahwa nilai pra-eksponensial paling rendah terdapat pada nilai k1 dan k2. Kedua nilai k ini digunakan pada reaksi pemutusan

hidrogen dari gas metana pada tahap satu (1). Reaksi ini merupakan reaksi yang paling lambat diantara seluruh reaksi yang terjadi di permukaan. Sedangkan reaksi paling cepat terjadi pada pemutusan hidrogen yang paling akhir (CH.S C.S). Nilai Energi aktivasi yang didapat dari hasil simulasi berkisar antara 50.274 j/mol hingga 54.398 J/mol untuk reaksi di permukaan. Nilai ini berbeda dengan nilai energi aktivasi milik Wulan (2011) yang sebesar 67,76 kJ/mol. Perbedaan nilai ini kemungkinan berasal dari model dan asumsi yang digunakan, dimana pada model ini asumsi Pseudo-Steady-State-Hypothesis (PSSH) tidak digunakan. 34 Penentuan parameter..., Marcho Rizal, FT UI, 2012

35

Selain itu pada model ini tahap pembatas laju tidak ditentukan dan seluruh model di selesaikan secara simultan untuk masing-masing spesies yang telibat dalam reaksi ini.

Nilai energi aktivasi pada penilitian ini menghasilkan nilai yang lebih

rendah dibandingkan dengan Zein (2004) yang menyatakan bahwa energi aktivasi

reaksi ini adalah 60 kJ/mol. Nilai yang dihasilkan penelitian ini cukup masuk akal jika dibandingkan dengan hasil Zein (2004) karena perbedaan nilai yang didapat

tidak besar dan juga karena perbedaan dari katalis yang digunakan. Data eksperimen yang digunakan untuk penelitian ini menggunakan katalis Ni-Cu-Al sedangkan zein menggunakan katalis Ni/TiO2.


BAB 5 KESIMPULAN

Kesimpulan dari Penelitian ini adalah sebagai berikut:

1. Model kinetika reaksi dekomposisi metana dengan katalis Ni-Cu-Al telah

berhasil dikembangkan dari mekanisme reaksi elementer milik Reyyan (2009)

2. Model Valid karena memiliki simpangan yang kecil (± 5%) 3. Model kinetika yang dikembangkan adalah model transien yang mempertimbangkan deaktivasi katalis 4. Berikut adalah nilai parameter kinetik yang didapat dari hasil penelitian:

k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11

A (molekul/sisi aktif.detik) 0,005253 0,000426 0,035295 0,006054 0,353101 0,046468 3,092562 0,540101 34,39805 5,139727 0,000219535

Energi Aktivasi (J/mol) 54.398,5 50.781,91 51.047,96 53.798,44 51.230,87 51.746,34 50.274,76 53.068,26 51.180,98 52.785.59 104.673,3


DAFTAR PUSTAKA Chen J., Li Y., Qin Y., Chang L. Formation of Bamboo-Shaped Carbon Filaments and Dependence of Their Morphology on Catalys Composition and

Reaction Conditions. Carbon 2001;39.1467-1475.

Chen, Jiuling, et al. Production of Cox-free Hydrogen and Nanocarbon by Direct Decomposition of Undiluted Methane on Ni-Cu-Alumina Catalyst. Applied

Catalysis: General 2004; 269:178-186. Constantinides A., Mostoufi N. 2000. Numerical Methods for Chemical Engineers with MATLAB Applications. New Jersey; Prentice Hall Dai L., 2006. Carbon Nanotechnology. Elsevier B.V. Hinink, J.P. A Quantitative of Transient Kinetic Experiment: Oxidation of CO by O2/NO on Pt. PhD. Thesis, The Hague, The Netherlands, 1995. Iijima, S, Carbon Nanotubes : Past, Present, Future. Physica B 2002; 323:1-5. Muradov, N. (2001). Catalysis of Methane Decomposition over Elemental Carbon. Catalysis Communication, 2, 89 – 94. Petterson, Joakim, Ove Hjortsberg. Hydrogen storage alternatives: A Technological and Economics Assessment. KFB Reports 1990; Vol 27. Resketenko, T.V., dkk. 2003. Carbon capacious Ni-Cu-Al2O3 catalyst for high temperature methane decompotition. Applied Catalyst A: General, 247, 51-63. Reyyan, et al. 2009. Partial Regeneration on Ni-based Catalyst for Hydrogen Production via Methane Cracking Part II: Modelling and Optimization. Turk J Chem 2009: 33(@ TUBITAK); p. 825-841. S.B. Sinnots, dkk, Chem. Phys. Lett. 315 (1999) 25

Slamet, Studi Kinetika Reaksi Reformasi Metana dengan Karbon DIoksida Menggunakan Katalis Ni/Al2O3. Tesis Magister UI 1996. Snoeck, et al. 1997. Kinetic Study of the Carbon Filament Formation by Methane Cracking on a Nickel Catalyst. England Song L. Pengaruh keasaman katalis berbasis Ni-Cu terhadap kinerja reaksi dekomposisi katalitik metana menjadi hydrogen dan nanokarbon. Skripsi. Departemen Teknik Gas dan Petrokimia FTUI; 2005.


38

Wulan, Praswasti PDK, Reaksi Dekomposisi Metana Dengan Katalis Ni-Cu-Al Untuk Produksi Carbon Nanotube : Kinetika Reaksi dan Pemodelan Reaktor. Disertasi. Departemen Teknik Kimia FTUI; 2011

Zein, S. H. S., Mohamed, A. R., and Sai, P. S. T., (2004), Kinetic Studies on

Catalytic Decomposition of Methane to Hydrogen and Carbon over Ni/Tio2

Catalyst,

Ind.

Eng.

Chem.

Res.,

43,

pp.

4864-4870


LAMPIRAN A PENURUNAN PERSAMAAN MODEL 1. CH4 + 2S  CH3.S + H.S

dikalikan

Ctotk’1 = k1 dan Ctotk’2 = k2

r1 = 2. CH3.S + S  CH2.S + H.S

dikalikan

Ctotk’3 = k3 dan Ctotk’4 = k4 39 Penentuan parameter..., Marcho Rizal, FT UI, 2012

40

r2 =

3. CH2.S + S CH.S + H.S

dikalikan


r3 =

4. CH.S + S C.S + H.S

dikalikan



41

r3 =

5. 2H.S H2 + 2S

dikalikan


r3 = 6. C.S  CNi,F + S

K11 dan k12 tidak dikonversi menjadi k’11 dan k’12 karena tidak peril mengalikan dengan Ctot/Ctot. satuannya berbeda


42

r6 =

Sehingga menghasilkan:


43

Persamaan Langmuir sederhana bisa digunakan :

Dimana θC.S adalah penutupan permukaan karbon dengan

= -45,217-

14.24 T (j/mol) Maka, konsesntrasi dari karbon pada bagian depan nikel x F (g carbon/ g Ni) bisa ditentukan dari persamaan:

(chen et al; snoeck et al) 7. CNi,F  CNi,R Laju difusi dari karbon melewati partikel nikel digunakan untuk menentukan berat spesifik karbon dalam model.


44

r7 = laju difusi karbon melalui nikel (mol C/ gcat h) l = rata panjang jalur, yang diasumsikan sebesar 2/3 dari diameter nickel

= difusivitas karbon dalam nikel = = konsentrasi dari karbon terlarut di nikel pada bagian depan Ni = konsentrasi dari karbon terlarut di nikel pada bagian belakang Ni ANi = Luas Spesifik dari partikel nikel (0.15m2/gcat)

8. Reaksi ini terjadi terlalu cepat, itulah mengapa diasumsikan bahwa setiap mol karbon yang mencapai bagian belakang partikel Ni berkontribusi pada struktur filament/nanotube. Sehinga CNi,R bisa diabaikan dalam model.

9.

-rC.S = rCp Konstanta laju dari reaksi enkapsulasi adalah parameter yang tidak diketahui dan diestimasi sepanjang model.


45

Laju net masing-masing spesies o CH4

Terdapat hanya pada reaksi tahap 1 saja (bolak balik)

-r1 =

(3.1)

o CH3.S

Terdapat pada reaksi tahap 1 dan 2 maka RCH3.S

=

r1

–

r2

= (3.2)

o CH2.S Terdapat pada reaksi tahap 2 dan 3 maka RCH2.S

= r2

– r3 = (3.3)

o CH.S Terdapat pada reaksi tahap 3 dan 4 maka RCH.S = r3-r4 = (3.4) o C.S RC.S = r4 =

(3.5)

o H.S Terdapat pada semua tahap reaksi. Maka: RH.S = r1 + r2 + r3 + r4 - 2r5 =

(3.6) o H2 Terdapat pada reaksi tahap 5 saja maka : RH2 = r5 =

(3.7)

o xF (3.8)


46

o xR

o CP

(3.9)

r8 =

dimana n = 3

(3.10)


LAMPIRAN B DATA EKSPERIMEN YANG DIGUNAKAN Seluruh Data pada lampiran ini merupakan data eksperimen milik Wulan (2011) Tabel 1 Data eksperimen jumlah karbon yang terbentuk pada suhu 750oC Waktu (detik) 0 334 335 335 336 336 338 338 338 346 1066 1067 1067 1067 1224 1225 1225 1225 1226 1226 1227 2021 2234 2234 2234

Jumlah Karbon yang terbentuk (gram) 0 0.075916 0.075916 0.075815 0.075714 0.075613 0.075613 0.075613 0.075411 0.086308 0.168037 0.168037 0.168037 0.168239 0.154819 0.154719 0.154517 0.154517 0.154517 0.193868 0.193868 0.260571 0.269248 0.269046 0.269046

2235 2235 2413 2413 2414 2890 2891 2891 3048 3049 3049 3049 3050 3085 3085 3085 3086 3086 3087 3087 3087 3088 3904 3904

0.269046 0.269147 0.340282 0.340282 0.34008 0.271165 0.271468 0.27167 0.366516 0.366516 0.366516 0.366415 0.366314 0.329284 0.329183 0.329183 0.329183 0.329385 0.329485 0.329485 0.329485 0.329385 0.401427 0.401427


Tabel 2 Data eksperimen jumlah karbon yang terbentuk pada suhu 700oC Waktu (detik) 0.00 70.00 100.00 150.00 279.00 279.00 279.00 306.00 306.00 307.00 307.00 308.00 308.00 404.00 611.00 612.00 612.00 612.00 613.00 1143.00 1143.00 1548.00 1548.00 2285.00 2285.00 2508.00 3543.00 3543.00 3544.00

Jumlah Karbon yang terbentuk (gram)

0 0.04 0.09 0.12 0.144536 0.144536 0.144738 0.153516 0.153617 0.153819 0.153919 0.15402 0.154121 0.161173 0.170849 0.171051 0.171152 0.171353 0.171454 0.18363 0.183832 0.19873 0.19888 0.208708 0.208506 0.227825 0.248253 0.248253 0.248253


Tabel 3 Data eksperimen jumlah karbon yang terbentuk pada suhu 650oC waktu (detik) 0 113 451 452 461 462 464 465 465 466 466 899 899 900 900 901 1054 1054 1055 1093 1121 1121 1122 1122 1357 1710 1710 1792 1792 1793


1794 1794 1794 1795 1795 1796 1796 2070 2230 2230 2231 2509 2509 3003 3004 3004 3051 3057 3057 3058 3058 3059 3407 3427 3427 3839 3839 3852 3853

0 0.065018 0.151809 0.152412 0.139909 0.140242 0.147486 0.147679 0.14754 0.147456 0.147456 0.150933 0.150782 0.150816 0.150665 0.150699 0.172781 0.172627 0.17204 0.193274 0.195079 0.195079 0.19511 0.19511 0.200193 0.215439 0.215764 0.219285 0.219606 0.219628

0.220292 0.220613 0.220934 0.221277 0.221438 0.22178 0.22178 0.238707 0.255881 0.245881 0.255896 0.259433 0.259433 0.261195 0.262422 0.26259 0.278669 0.274013 0.274013 0.274027 0.273692 0.273539 0.320585 0.297449 0.297787 0.318793 0.319134 0.326607 0.312278


Tabel 4 Data eksperimen jumlah karbon yang terbentuk pada suhu 600oC waktu (detik) 0 42 42 43 43 44 44 44 45 45 46 46 100 100 101 101 174 174 174 175 176 176 177 178 179 220 250 288 300 300 400 600 700 800 900 1008 1009 1009

Jumlah Karbon yang terbentuk (gram) 0 0.025676 0.025575 0.025575 0.025474 0.025474 0.025575 0.025575 0.025676 0.025776 0.025776 0.025978 0.042425 0.042526 0.042627 0.042829 0.054533 0.054533 0.054533 0.054331 0.054331 0.054331 0.054331 0.054634 0.054836 0.065 0.066 0.07 0.087 0.09 0.11 0.13 0.14 0.135 0.137 0.142921 0.142921 0.143022

1157 1157 1220 1220 1221 1221 1221 1222 1222 1223 1223 1223 1224 1224 1225 1225 1226 1368 1369 1369 1370 1370 1370 1371 1371 1376 1384 1429 1430 1437 1450 1450 1503 1504 1504 1504 1505 1505 1506 1506 1506


0.152002 0.152204 0.1529 0.1529 0.1529 0.1529 0.1529 0.1529 0.153 0.153101 0.153202 0.153303 0.153505 0.153505 0.153606 0.153505 0.153202 0.157844 0.157945 0.158045 0.158045 0.158045 0.158045 0.158045 0.158247 0.161577 0.166622 0.168438 0.153337 0.15142 0.149604 0.149705 0.163326 0.163427 0.163427 0.163326 0.163225 0.163024 0.163024 0.163024 0.163125

51

1507 1525 1526 1538 1538 1564 1564 1564 1710 1805 1805 1805 1854 1855 1855 1882 1883 1883 1901 1901 1901 1902 1902 1903 1903 1903 2011 2011 2012 2105 2106 2266 2267 2394 2395 2395 2396 2396 2396 2397 2398 2398

0.163225 0.164582 0.164784 0.160546 0.160546 0.163999 0.163797 0.163696 0.167015 0.165087 0.164885 0.164784 0.159213 0.159112 0.159112 0.167083 0.166982 0.166982 0.16002 0.160121 0.160222 0.160423 0.160423 0.160524 0.160524 0.160524 0.169483 0.169483 0.169483 0.15678 0.156679 0.181602 0.181702 0.18029 0.180189 0.180189 0.180189 0.180189 0.180189 0.179987 0.179785 0.179684

2398 2399 2399 2400 2400 2540 2555 2556 2556 2557 2557 2563 2563 2564 2564 2564 2568 2568 2569 2569 2569 2651 2697 2697 2697 2698 2803 2804 2804 2804 2805 2806 2807 2807 2808 2808 2808 2808 2809 2809 2810 2939

0.179483 0.179382 0.179281 0.179281 0.17918 0.198452 0.195122 0.195425 0.195627 0.195728 0.195728 0.193508 0.193407 0.193306 0.193104 0.193003 0.193003 0.192902 0.192902 0.192902 0.192902 0.196333 0.206927 0.207129 0.20723 0.20723 0.210762 0.210762 0.210863 0.211064 0.211165 0.211266 0.211064 0.210963 0.210762 0.210661 0.204102 0.204102 0.204102 0.204102 0.204001 0.222971


52

2939 2940 2954 2954 2975 3010 3010 3011 3011 3012 3012 3012 3100 3100 3101 3159 3160 3160 3160 3279 3280 3490 3490 3496 3523 3524 3524 3525 3525

0.222971 0.22287 0.220347 0.220246 0.214596 0.219944 0.219944 0.219944 0.219843 0.219742 0.21954 0.21954 0.228823 0.228924 0.229125 0.233767 0.233868 0.233969 0.23417 0.233666 0.233666 0.242444 0.242545 0.241738 0.242949 0.24305 0.243151 0.243151 0.243251

3525 3526 3526 3527 3527 3527 3528 3534 3595 3596 3596 3597 3597 3597 3656 3657 3657 3658 3677 3753 3835 3836 3836 3837 3837 3837 3918

0.243352 0.243352 0.243352 0.243352 0.243352 0.243352 0.243554 0.245774 0.246783 0.246783 0.246884 0.246985 0.246985 0.246985 0.248599 0.248902 0.249104 0.249104 0.238812 0.191893 0.245774 0.245875 0.246077 0.246178 0.246278 0.246379 0.246581


Tabel 5 Data eksperimen jumlah karbon yang terbentuk pada suhu 550oC Waktu (Detik) 0 0 143 143 144 299 468 468 890 1235 1445 1636 1636 2209 2698 3002 3002 3003 3003 3450 3451 3451 3451 3532 3668 3667 3667


0 0 0.038893 0.054937 0.055037 0.055239 0.10598 0.168709 0.15998 0.1775 0.199295 0.21567 0.215771 0.232327 0.263753 0.275021 0.255122 0.255324 0.255526 0.27389 0.273789 0.273688 0.273486 0.273082 0.3123 0.300231 0.300231


Tabel 6 Data eksperimen jumlah karbon yang terbentuk pada suhu 550oC Waktu (Detik) 0 63 63 63 66 66 67 67 67 68 68 127 127 127 128 431 432 433 500 500 765 777 998 1188 1188 1188 1189 1189 1190 1190 1250 1252 1252 1255 1371 1374 1390 1390 1390

Jumlah Karbon yang terbentuk (gram) 0 0.057762 0.057964 0.058064 0.058165 0.058266 0.058367 0.058569 0.05867 0.058771 0.058771 0.060284 0.060385 0.060385 0.060284 0.060183 0.060082 0.060082 0.074915 0.075116 0.080342 0.08524 0.09023 0.094915 0.095116 0.095117 0.095217 0.095318 0.095419 0.095419 0.09552 0.09552 0.115196 0.10783 0.108032 0.108334 0.109444 0.110453 0.121149

1391 1391 1392 1392 1392 1393 1400 1821 1821 1835 1835 1848 1848 1849 1849 1850 1850 1851 2004 2005 2005 2005 2006 2006 2008 2009 2039 2077 2077 2312 2438 2438 2520 2520 2520 2521 3488 3488 3489 3634 3635


0.12125 0.121451 0.121552 0.121754 0.121754 0.121754 0.121754 0.104075 0.104277 0.108111 0.108212 0.107808 0.107808 0.10801 0.10801 0.108212 0.108212 0.108212 0.109624 0.109725 0.109725 0.109725 0.109624 0.109423 0.120925 0.120824 0.121934 0.125567 0.125667 0.130208 0.142619 0.142518 0.145948 0.145948 0.145948 0.145847 0.15281 0.15281 0.15281 0.148067 0.148067

LAMPIRAN D PLOT ANTARA LN K TERHADAP 1/T -11.5 0.0009 -12

0.001

0.0011

0.0012

0.0013

0.0014

-12.5 ln k

-13

-13.5

y = -6543.5x - 5.2493 R² = 0.9909

-14

1/T

Dari plot nilai ln k1 terhadap 1/T, didapatlah nilai untuk pra eksponensial untuk k1 adalah sebesar0,005253 dengan nilai energi aktivasi sebesar 54.398,5 J/mol. -13.5 0.0009 -14

0.001

0.0011

0.0012

0.0013

0.0014

ln k

-14.5 -15 -15.5 -16

y = -6108.9x - 7.7621 R² = 0.9741 1/T

Dari plot nilai ln k2 terhadap 1/T, didapatlah nilai untuk pra eksponensial untuk k2 adalah sebesar 0,000426 dengan nilai energi aktivasi sebesar 50.781,91 J/mol.


56

ln k

-8 0.0009 -8.5

0.001

0.0011

-9

-9.5

-10

-10.5

0.0012

0.0013

0.0014

-11

y = -6140.6x - 3.3445 R² = 0.9743

-11.5 -12

1/T


-11 0.0009

0.001

0.0011

0.0012

0.0013

0.0014

-11.5

ln k

-12 -12.5 -13 -13.5 -14

y = -6460.7x - 5.1077 R² = 0.9826

1/T



57

-5 0.0009 -5.5

0.001

-6

0.0012

0.0013

0.0014

-6.5 ln k

0.0011

-7

-7.5

-8 -8.5

y = -6162.3x - 1.0417 R² = 0.9818

-9 -9.5

1/T


-8 0.0009 -8.5

0.001

0.0011

0.0012

0.0013

0.0014

ln k

-9 -9.5 -10 -10.5 -11 -11.5

y = -6224.3x - 3.0695 R² = 0.9712

1/T



58

-4 0.0009

0.001

0.0011

0.0012

0.0013

0.0014

-4.5

ln k

-5

-5.5

-6

-6.5

y = -6047.4x + 1.1298 R² = 0.9981

-7

1/T

Dari plot nilai ln k7 terhadap 1/T, didapatlah nilai untuk pra eksponensial untuk k7 adalah sebesar 0,000426 dengan nilai energi aktivasi sebesar 50.781,91 J/mol. -6 0.0009 -6.5

0.001

0.0011

0.0012

0.0013

0.0014

-7

ln k

-7.5 -8 -8.5 -9 -9.5

y = -6383x - 0.6167 R² = 0.9661

1/T



59

ln k

-1 0.0009 -1.5

0.001

0.0011

-2

-2.5

-3

-3.5

0.0012

0.0013

0.0014

-4 y = -6156.8x + 3.5381 R² = 0.9867

-4.5 -5

1/T


-3 0.0009 -3.5

0.001

0.0011

0.0012

0.0013

0.0014

-4

ln k

-4.5 -5 -5.5 -6 -6.5 -7

y = -6349.8x + 1.6371 R² = 0.9872

1/T

Dari plot nilai ln k10 terhadap 1/T, didapatlah nilai untuk pra eksponensial untuk k10 adalah sebesar 5,139727 dengan nilai energi aktivasi sebesar 52.785.59 J/mol.


60

-3 0.0009

0.001

0.0011

0.0012

0.0013

0.0014

-4

ln k

-5

-6

-7

y = -12590x + 8.4243 R² = 0.9817

-8 -9

1/T



LAMPIRAN C PROGRAM YANG DIGUNAKAN UNTUK PENELITIAN function [x,y] = CBT (k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,x0,y0) k1 = input('nilai k1 = '); k2 = input('nilai k2 = '); k3 = input('nilai k3 = '); k4 = input('nilai k4 = '); k5 = input('nilai k5 = '); k6 = input('nilai k6 = '); k7 = input('nilai k7 = '); k8 = input('nilai k8 = '); k9 = input('nilai k9 = '); k10 = input('nilai k10 = '); [x,y] = ode23s(@persamaan,[O 3600],[.2668 0 0 0 0 0 .7332 1]); function dy = persamaan (x,y,k) dy = zeros (8,1); dy(1) = k2*y(2)*y(7)-k1*y(1)*y(8)^2 ; dy(2) = k1*y(1)*y(8)^2 - k2*y(2)*y(6) - k3*y(2)*y(8) + k4*y(3)*y(6); dy(3) = k3*y(2)*y(8)^2-k4*y(3)*y(6) - k5*y(3)*y(8) + k6*y(4)*y(6); dy(4) = k5*y(3)*y(8)^2-k6*y(4)*y(6) - k7*y(4)*y(8)+k8*y(5)*y(6); dy(5) = k7*y(4)*y(8) - k8*y(5)*y(6); dy(6) = k1*y(1)*y(8)^2-k2*y(2)*y(6)+k3*y(2)*y(8)^2k4*y(3)*y(6)+k5*y(3)*y(8)^2-k6*y(4)*y(6)+k7*y(4)*y(8) k8*y(5)*y(6)+2*(k9*y(7)^2-k10*y(7)*y(8)^2); dy(7) = k9*y(7)^2-k10*y(7)*y(8)^2; dy(8) = 1-y(2)-y(3)-y(4)-y(5)-y(6); end disp(' Nonlinear Regression Analysis'), disp(' ') % Input experimental data disp(' ') datain = 3; disp(' ') switch datain case 1 fdata = input(' Name of file for storing the data = '); if exist ([fdata,'.mat']) disp(' Data file already exists. Use another name.') fdata = input(' Name of file for storing the data = '); end v = input(' Number of dependent variables = '); for m = 1:v countpoints=0; fprintf('\n Variable %2d\n',m) datasets=input(' How many data sets for this variable? = '); for nset = 1:datasets fprintf('\n Data set %2d\n',nset) npoints = input(' How many points in this set? = '); disp(' ') for np = 1:npoints countpoints=countpoints+1; fprintf(' Point %2d\n',np) xm = input(' Enter independent value = '); ym = input(' Enter dependent value = ');


62

x(countpoints,m)=xm; y(countpoints,m)=ym; end end n(m) = countpoints;

end eval(['save ',fdata ' n x y']) case 2 fdata = input(' Name of file for storing the data = '); if exist ([fdata,'.mat']) disp(' Data file already exists. Use another name.') fdata = input(' Name of file for storing the data = '); end v = input(' Number of dependent variables = '); for m = 1:v fprintf('\n Variable %2d\n',m) xm = input(' Vector of independent variable = '); x(1:length(xm),m) = (xm(:).')'; ym = input(' Vector of dependent variable = '); y(1:length(ym),m) = (ym(:).')'; n(m) = length(ym); end eval(['save ',fdata ' n x y']) case 3 fname = 'data.mat'; eval(['load ', fname]) [r,v] = size(y); end disp(' ') namex = 'time,seconds'; for m = 1:v nam = ['ch4' 'ch3s' 'ch2s' 'chs' 'cs' 'hs' 'h2' 'tetav']; namey(m,1:length(nam))= nam; end % Input type of equation(s) eqin = 2; disp(' ') switch eqin case 1 % Algebraic equation fnctn = input(' Name of M-file containing algebraic equation(s) = '); ODEfile = []; x0 = []; y0 = []; case 2 % ODE ODEfile = 'persamaan'; disp(' ') x0 = 0; y0 = [y(1,1) y(2,1) y(3,1) y(4,1) y(5,1) y(6,1) y(7,1) y(8,1)]; fnctn = []; otherwise% Algebraic and ODE fnctn = input(' Name of M-file containing algebraic equation(s) = '); ODEfile = input(' Name of M-file containing differential equation(s) = '); disp(' ')


63

x0 = input(' Value of independent variable at boundary condition = '); y0 = input(' Value(s) of dependent variable(s) at boundary condition = '); end % Input method of solution

method = 1;

disp(' ') b0 = [k1,k2,k3,k4,k5,k6,k7,k8,k9,k10]; disp(' ') trace = 0; if trace == 0 disp(' ') disp(' Please wait for final results'), disp(' ') end % Regression [b, yc, w, JTJ] = NLR(b0, n, x, y, fnctn, ODEfile, x0, y0, ... method, [], trace); if trace == 1 disp(' ') disp(' Please wait for final results') end % Statistical properties [sd, cl] = statistics(b, n, x, y, yc, w, JTJ, 95, fout); % Displaying final results disp(' ') disp(' **********************************************************') disp(' Final Results') disp(' **********************************************************') disp(' ') disp(' No. Parameter Standard 95% Confidence interval ') disp(' deviation for the parameters ') disp(' lower value upper value') for m = 1:length(b) fprintf(' %2d %10.4e %10.4e %10.4e %10.4e\n',... m,b(m),sd(m),cl(m,:)) end disp(' ') disp(' **********************************************************') % Plotting the results for m = 1:v figure(m) [xx,loc] = sort(x(1:n(m),m)); plot(x(1:n(m),m),y(1:n(m),m),'o',xx,yc(loc,m))


64

xlabel(namex) ylabel(namey(m,:)) end fclose('all'); end


PENENTUAN PARAMETER KINETIK DEKOMPOSISI KATALITIK METANA DENGAN KATALIS Ni-Cu-Al UNTUK PRODUKSI NANO KARBON SKRIPSI

Recommend Documents