13.01.2013 Views

tc süleyman demirel üniversitesi fen bilimleri enstitüsü metal hidrür ...

tc süleyman demirel üniversitesi fen bilimleri enstitüsü metal hidrür ...

tc süleyman demirel üniversitesi fen bilimleri enstitüsü metal hidrür ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

T.C.<br />

SÜLEYMAN DEMĐREL ÜNĐVERSĐTESĐ<br />

FEN BĐLĐMLERĐ ENSTĐTÜSÜ<br />

METAL HĐDRÜR REAKTÖR TASARIM PARAMETRELERĐNĐN<br />

HĐDROJEN DEPOLAMA ÖZELLĐKLERĐNE ETKĐSĐNĐN<br />

DENEYSEL OLARAK ĐNCELENMESĐ<br />

Muhammet KAYFECĐ<br />

Danışman: Doç. Dr. Fevzi BEDĐR<br />

II. Danışman: Doç. Dr. Hüseyin KURT<br />

DOKTORA TEZĐ<br />

MAKĐNA MÜHENDĐSLĐĞĐ ANABĐLĐM DALI<br />

ISPARTA – 2011


ĐÇĐNDEKĐLER<br />

i<br />

Sayfa<br />

ĐÇĐNDEKĐLER.............................................................................................................. i<br />

ÖZET .. ........................................................................................................................ iv<br />

ABSTRACT ................................................................................................................. v<br />

TEŞEKKÜR ................................................................................................................ vi<br />

ŞEKĐLLER DĐZĐNĐ ....................................................................................................vii<br />

ÇĐZELGELER DĐZĐNĐ................................................................................................. x<br />

SĐMGELER VE KISALTMALAR DĐZĐNĐ................................................................xi<br />

1.GĐRĐŞ......................................................................................................................... 1<br />

1.1. Hidrojen Üretimi.................................................................................................... 3<br />

1.1.1. Fosil yakıt yöntemleri......................................................................................... 5<br />

1.1.2. Elektroliz ............................................................................................................ 5<br />

1.1.3. Isıl kimyasal yöntem........................................................................................... 6<br />

1.2. Hidrojenin Depolanması........................................................................................ 6<br />

1.2.1. Sıkıştırılmış gaz olarak depolama....................................................................... 7<br />

1.2.2. Sıvı olarak depolama .......................................................................................... 9<br />

1.2.3. Metal <strong>hidrür</strong> esaslı depolama............................................................................ 10<br />

1.2.4. Karbon nanotüplerde depolama........................................................................ 11<br />

1.2.5. Cam mikro kürelerde depolama ....................................................................... 12<br />

1.2.6. Bor esaslı depolama.......................................................................................... 13<br />

1.2.7. Mağaralarda depolama ..................................................................................... 14<br />

1.3. Metal Hidrürler .................................................................................................... 14<br />

1.3.1. Metal <strong>hidrür</strong> termodinamiği.............................................................................. 22<br />

1.3.1.1. Aktivasyon..................................................................................................... 22<br />

1.3.1.2. Histeresiz ....................................................................................................... 23<br />

1.3.1.3. Plato eğimi..................................................................................................... 24<br />

1.3.1.4. Tersinir kapasite ............................................................................................ 25<br />

1.3.1.5. Hidrürleşme kinetikleri.................................................................................. 25<br />

1.3.1.6. Metal <strong>hidrür</strong> ömrü.......................................................................................... 26


1.3.1.7. Kristal çatlaması ............................................................................................ 27<br />

1.3.2. Metal <strong>hidrür</strong>ler arası bağlar .............................................................................. 27<br />

1.3.2.1. Đyonik <strong>hidrür</strong>ler.............................................................................................. 28<br />

1.3.2.2. Kovalent <strong>hidrür</strong>ler.......................................................................................... 29<br />

1.3.2.3. Metalik <strong>hidrür</strong>ler............................................................................................ 30<br />

1.3.3. Metal <strong>hidrür</strong> Çeşitleri........................................................................................ 31<br />

1.3.3.1. AB5 inter<strong>metal</strong>ik bileşikleri........................................................................... 32<br />

1.3.3.2. AB2 inter<strong>metal</strong>ik bileşikleri ........................................................................... 34<br />

1.3.3.3. AB inter<strong>metal</strong>ik bileşikleri ............................................................................ 35<br />

1.3.3.4. Diğer inter<strong>metal</strong>ik bileşikler.......................................................................... 36<br />

1.3.3.5. Katı çözelti alaşımları.................................................................................... 36<br />

2.KAYNAK ÖZETLERĐ............................................................................................ 38<br />

3.MATERYAL VE YÖNTEM................................................................................... 53<br />

3.1. Materyal............................................................................................................... 53<br />

3.1.1. Metal <strong>hidrür</strong> reaktör.......................................................................................... 53<br />

3.1.2. Mekanik öğütücü .............................................................................................. 55<br />

3.1.3. Split fırın........................................................................................................... 57<br />

3.1.4. Vakum pompası................................................................................................ 58<br />

3.1.5. Sıcaklık tarayıcı cihaz....................................................................................... 59<br />

3.1.6. Deney tesisatı.................................................................................................... 59<br />

3.1.7. Hassas terazi ..................................................................................................... 60<br />

3.1.8. Hidrojen gazı .................................................................................................... 61<br />

3.1.9. Azot gazı........................................................................................................... 61<br />

3.2. Yöntem ................................................................................................................ 62<br />

3.2.1. Metal <strong>hidrür</strong> reaktör kanat etkenliği ................................................................. 62<br />

3.2.1.1. Kanat verimi .................................................................................................. 63<br />

3.2.1.2. Toplam yüzey verimi..................................................................................... 64<br />

3.2.2. Deneyler............................................................................................................ 66<br />

3.2.3. Yapay sinir ağları.............................................................................................. 67<br />

ii


3.2.3.1. Yapay sinir ağlarının özellikleri .................................................................... 68<br />

3.2.3.2. Yapay sinir alanlarının kullanım alanları ...................................................... 69<br />

3.2.3.3. Yapay sinir ağı yapıları.................................................................................. 70<br />

3.2.3.4. Yapay sinir ağı ile modelleme....................................................................... 71<br />

4.ARAŞTIRMA BULGULARI.................................................................................. 72<br />

4.1. Deneysel Çalışmalardan Elde Edilen Bulgular ................................................... 72<br />

4.2. YSA Kullanılarak Yapılan Modellemeden Elde Edilen Bulgular....................... 91<br />

5. TARTIŞMA VE SONUÇ..................................................................................... 100<br />

5.1. Sonuçlar............................................................................................................. 100<br />

5.2. Öneriler.............................................................................................................. 102<br />

6.KAYNAKLAR...................................................................................................... 104<br />

EKLER ..................................................................................................................... 111<br />

EK-1 Kullanılan Alaşım ve Gazların Kimyasal ve Fiziksel Özellikleri................... 111<br />

EK-2 Deneysel Sistem Resimleri ............................................................................. 113<br />

ÖZGEÇMĐŞ.............................................................................................................. 117<br />

iii


ÖZET<br />

Doktora Tezi<br />

METAL HĐDRÜR REAKTÖR TASARIM PARAMETRELERĐNĐN<br />

HĐDROJEN DEPOLAMA ÖZELLĐKLERĐNE ETKĐSĐNĐN DENEYSEL<br />

OLARAK ĐNCELENMESĐ<br />

Muhammet KAYFECĐ<br />

Süleyman Demirel Üniversitesi<br />

Fen Bilimleri Enstitüsü<br />

Makina Mühendisliği Anabilim Dalı<br />

Danışman: Doç. Dr. Fevzi BEDĐR<br />

Bu çalışmada, <strong>metal</strong> <strong>hidrür</strong> reaktörlerde hidrojen şarj/deşarj işlemi deneysel olarak<br />

incelenmiş ve yapay sinir ağları (YSA) kullanılarak modellenmiştir. Deneysel<br />

çalışmalarda kanatçıksız ve kanatçıklı olmak üzere iki farklı geometride <strong>metal</strong> <strong>hidrür</strong><br />

reaktör tasarlanmış ve imal edilmiştir. Hidrojen şarj/deşarj işlemi üzerine etkili olan<br />

çalışma şartları, reaktör geometrisi ve alaşım gibi parametreler araştırılmıştır. Bu<br />

amaçla yapılan deneysel çalışmalarda reaktöre beslenen hidrojen şarj basıncı 2, 4, 6<br />

ve 8 bar değerlerine ayarlanarak reaktörlerin farklı noktalarındaki sıcaklıklarının<br />

zamanla değişimi incelenmiştir.<br />

Yapılan deneylerde maksimum sıcaklıklar z=60 mm noktasında, 2, 4, 6 ve 8 bar şarj<br />

basınçlarında Reaktör 1’de sırasıyla 88.77, 102.92, 116.24 ve 128.72 °C olurken,<br />

Reaktör 2’de 81.33, 98.1, 113.27 ve 121.52 °C olarak ölçülmüştür. Depolanan<br />

hidrojen kütlesini belirlemek için yapılan çalışmalarda ise 4, 6 ve 8 bar basınçlarda<br />

sırasıyla Reaktrör 1’de 0.88, 1.17, 1.32 g, Reaktör 2’de ise 1.13, 1.25, 1.55 g olarak<br />

ölçülmüştür. Oluşturulan YSA modelinde ise, Reaktör 1 ve Reaktör 2 için depolanan<br />

hidrojen kütlesinin yaklaşık % 99 doğrulukta tahmin edilebileceği hesaplanmıştır.<br />

YSA modelinde R 2 , RMSE ve MAPE değerleri Reaktör 1 için sırasıyla 0.9991,<br />

0.000062 ve 0.02146 olarak bulunurken, Reaktör 2 için sırasıyla 0.9984, 0.01258 ve<br />

0.00424 olarak bulunmuştur. Ayrıca LaNi5 alaşımlarına Al ilavesi reaksiyonu<br />

hızlandırmış ve <strong>hidrür</strong> reaktörlerde şarj süresini kısaltmıştır.<br />

Anahtar Kelimeler: Hidrojen depolama, <strong>metal</strong> <strong>hidrür</strong> reaktör, LaNi4.75Al0.25<br />

alaşımları, reaktör tasarımı, YSA<br />

2011, 118 sayfa<br />

iv


ABSTRACT<br />

Ph. D. Thesis<br />

EXPERIMENTAL STUDY FOR EFFECT OF DESIGN PARAMETERS OF<br />

METAL HYDRIDE REACTOR ON HYDROGEN STORAGE PROPERTIES<br />

Muhammet KAYFECĐ<br />

Süleyman Demirel University<br />

Graduate School of Applied and Natural Sciences<br />

Department of Mechanical Engineering<br />

Supervisor: Assoc. Prof. Dr. Fevzi BEDĐR<br />

In this study, hydrogen charge/discharge process in <strong>metal</strong> hydride reactors has been<br />

investigated experimentally and modeled by using Artificial Neural Network (ANN).<br />

Metal hydride reactors have two different geometry with finned and not finned were<br />

designed and manufactured at the experimentally study. Parameters affecting the<br />

hydrogen charge/discharge process, such as operation conditions, reactor geometry<br />

and alloys were investigated. For this purpose, hydrogen feeding pressure to the<br />

reactor was adjusted to 2, 4, 6 and 8 bar, and the variation of temperature on different<br />

points of the reactor and absorbed hydrogen mass by the time was investigated.<br />

The maximum temperature values on z=60 mm point of the Reactor 1 for 2, 4, 6 and<br />

8 bar charging pressures were 88.77, 102.92, 116.24 and 128.72 °C respectively<br />

while the same values for the same point of Rector 2 were measured as 81.33, 98.1,<br />

113.27 and 121.52 °C respectively. On the studies, to determine absorbed hydrogen<br />

mass was measured as 0.88, 1.17, 1.32 g in Reactor 1 and 1.13, 1.25, 1.55 g in<br />

Reactor 2 for the charging pressures of 4, 6 and 8 bar respectively. It was calculated<br />

that the amount of absorbed hydrogen mass in Reactor 1 and 2 can be estimated by<br />

the developed ANN model with an accuracy of about 99%. In ANN model R 2 ,<br />

RMSE and MAPE values for Reactor 1 were found as 0.999, 0.00006 and 0.0214<br />

respectively while for the Reactor 2 these values were found as 0.998, 0.01258 and<br />

0.0042 respectively. Also, Al addition to the LaNi5 alloys have accelerated the<br />

reaction and decreased the charging time on hydride reactors.<br />

Key Words: Hydrogen storage, <strong>metal</strong> hydride reactor, LaNi4.75Al0.25 alloys, reactor<br />

design, ANN.<br />

2011, 118 pages<br />

v


TEŞEKKÜR<br />

Bu araştırma için beni yönlendiren, karşılaştığım zorlukları bilgi ve tecrübesi ile<br />

aşmamda yardımcı olan değerli danışmanlarım Doç. Dr. Fevzi BEDĐR ve Doç. Dr.<br />

Hüseyin KURT hocalarıma teşekkür ederim. Deneysel çalışmaların yapılmasında<br />

imkânlarını ve desteklerini esirgemeyen Niğde Üniversitesi Mühendislik-Mimarlık<br />

Fakültesi Makine Mühendisliği Bölümü Öğretim Üyeleri Prof. Dr. Mahmut D. MAT,<br />

Doç. Dr. Yüksel KAPLAN hocalarıma ve HYTEM laboratuarı çalışanlarına<br />

teşekkürlerimi sunarım. Ayrıca araştırmanın yürütülmesi aşamasında maddi ve<br />

manevi desteklerini esirgemeyen değerli arkadaşlarım Arş. Gör. Dr. Ali KEÇABAŞ,<br />

Arş. Gör. Engin GEDĐK, Öğr. Gör. Ali UYSAL, Öğr. Gör. Şafak ATAŞ ve adını<br />

burada sayamadığım bütün arkadaşlarıma ve hocalarıma teşekkür ederim.<br />

1710-D–08 No`lu Proje ile tezimi maddi olarak destekleyen Süleyman Demirel<br />

Üniversitesi Bilimsel Araştırma Projeleri Koordinasyon Birimine teşekkür ederim.<br />

Ayrıca bugünlere gelmemde çok büyük emeği bulunan annem, babam ve<br />

kardeşlerime tezimin her aşamasında beni yalnız bırakmayan eşim Nurcan ve kızım<br />

Zeynep Nisa’ya sonsuz sevgi ve saygılarımı sunarım.<br />

vi<br />

Muhammet KAYFECĐ<br />

ISPARTA, 2011


ŞEKĐLLER DĐZĐNĐ<br />

Şekil 1.1. Hidrojen üretim zinciri................................................................................. 4<br />

Şekil 1.2. Sıkıştırılmış hidrojen gazı depolama tankı (kompozit silindir) ................... 9<br />

Şekil 1.3. Sıvı hidrojen depolama tankı ..................................................................... 10<br />

Şekil 1.4. Tek ve çok duvarlı karbon nanotüpler ....................................................... 11<br />

Şekil 1.5. Mikro cam kürelerde hidrojen depolama................................................... 12<br />

Şekil 1.6. Sodyum bor <strong>hidrür</strong> esaslı hidrojen depolama ............................................ 13<br />

Şekil 1.7. Bor esaslı hidrojen depolama..................................................................... 14<br />

Şekil 1.8.a) Metal üzerinde hidrojen kemisorpsiyonu, b) Moleküler potansiyel<br />

kuyusu ve atomik hidrojen ...................................................................... 16<br />

Şekil 1.9. Metal <strong>hidrür</strong>lerde faz dönüşüm şeması...................................................... 17<br />

Şekil 1.10. LaNi5 hidrojen absorpsiyon hızı a) Polikristalize, b) Nano-kristalize,<br />

c) katalizörlü nanokristalize .................................................................... 19<br />

Şekil 1.11. Metal <strong>hidrür</strong> teorik P-C eğrisi.................................................................. 20<br />

Şekil 1.12. P-C izotermleri gerçek diyagramı............................................................ 20<br />

Şekil 1.13. Şematik van’t Hoff grafiği....................................................................... 21<br />

Şekil 1.14. Đkili <strong>hidrür</strong>ler tablosu ............................................................................... 28<br />

Şekil 1.15. Bazı doğal elementler için van’t Hoff diyagramı .................................... 32<br />

Şekil 1.16. Çeşitli AB5 alaşımları için van’t Hoff diyagramı .................................... 33<br />

Şekil 1.17. Bazı AB2 alaşımları için van’t Hoff diyagramı ....................................... 34<br />

Şekil 1.18. Bazı AB alaşımları için van’t Hoff diyagramı......................................... 35<br />

Şekil 3.1. Reaktörlerin kesit görünüşleri, a) Reaktör 1 ve b) Reaktör 2 .................... 54<br />

Şekil 3.2. Reaktör 1 ve Reaktör 2 kesit görünüşleri ve imal ölçüleri......................... 54<br />

Şekil 3.3. Termokuplların <strong>metal</strong> <strong>hidrür</strong> reaktörler üzerine yerleşimi......................... 55<br />

Şekil 3.4. Metal <strong>hidrür</strong> tozlarını öğütmek için kullanılan mekanik öğütücü ............. 56<br />

Şekil 3.5. Split fırın.................................................................................................... 57<br />

Şekil 3.6. Vakum pompası ......................................................................................... 58<br />

Şekil 3.7. Sıcaklık tarayıcı cihaz................................................................................ 59<br />

Şekil 3.8. Deney tesisatı............................................................................................. 60<br />

vii


Şekil 3.9. Hassas terazi............................................................................................... 60<br />

Şekil 3.10. Dikdörtgen profilli dairesel kanat verimleri ............................................ 63<br />

Şekil 3.11. Dikdörtgen ve dairesel kanat dizileri ....................................................... 65<br />

Şekil 3.12. Đleri beslemeli YSA yapısı....................................................................... 71<br />

Şekil 4.1. 2 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2 ............................................... 74<br />

Şekil 4.2. 4 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2 ............................................... 75<br />

Şekil 4.3. 6 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2 ............................................... 77<br />

Şekil 4.4. 8 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2 ............................................... 78<br />

Şekil 4.5. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z = 5 mm noktasında<br />

zamana bağlı sıcaklık değişimleri. .......................................................... 79<br />

Şekil 4.6. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z = 20 mm noktasında<br />

zamana bağlı sıcaklık değişimleri ........................................................... 80<br />

Şekil 4.7. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z = 40 mm noktasında<br />

zamana bağlı sıcaklık değişimleri ........................................................... 81<br />

Şekil 4.8. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z = 60 mm noktasında<br />

zamana bağlı sıcaklık değişimleri ........................................................... 82<br />

Şekil 4.9. Reaktör 2’ de 2, 4, 6 ve 8 bar basınçlarda ve z = 5 mm noktasında<br />

zamana bağlı sıcaklık değişimleri. .......................................................... 83<br />

Şekil 4.10. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z = 20 mm noktasında<br />

zamana bağlı sıcaklık değişimleri ........................................................... 84<br />

Şekil 4.11. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z = 40 mm noktasında<br />

zamana bağlı sıcaklık değişimleri ........................................................... 85<br />

Şekil 4.12. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z = 60 mm noktasında<br />

zamana bağlı sıcaklık değişimleri ........................................................... 86<br />

Şekil 4.13. Reaktör 1 ve Reaktör 2 reaktörlerde 4 bar basınçta depolanan hidrojen<br />

kütlesinin zamana bağlı değişimi ............................................................ 87<br />

Şekil 4.14. Reaktör 1 ve Reaktör 2 reaktörlerde 6 bar basınçta depolanan hidrojen<br />

kütlesinin zamana bağlı değişimi ............................................................ 88<br />

Şekil 4.15. Reaktör 1 ve Reaktör 2 reaktörlerde 8 bar basınçta depolanan hidrojen<br />

kütlesinin zamana göre değişimi............................................................. 88<br />

Şekil 4.16. Reaktör 1’de 4, 6 ve 8 bar basınçlarda depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi............................................................................. 89<br />

viii


Şekil 4.17. Reaktör 2’de 4, 6 ve 8 bar basınçlarda depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi............................................................................. 90<br />

Şekil 4.18. Reaktör 1 için oluşturulan YSA yapısı .................................................... 92<br />

Şekil 4.19. Reaktör 2 için oluşturulan YSA yapısı .................................................... 93<br />

Şekil 4.20. Reaktör 1 için YSA eğitim verileri ile deneysel verilerin<br />

karşılaştırılması ....................................................................................... 97<br />

Şekil 4.21. Reaktör 2 için YSA eğitim verileri ile deneysel verilerin<br />

karşılaştırılması ....................................................................................... 97<br />

Şekil 4.22. Reaktör 1’de depolanan hidrojen kütlesinin zamana ve basınca bağlı<br />

değişimi ................................................................................................... 98<br />

Şekil 4.23. Reaktör 1’de depolanan hidrojen kütlesinin zamana ve basınca bağlı<br />

değişimi ................................................................................................... 98<br />

ix


ÇĐZELGELER DĐZĐNĐ<br />

Çizelge 1.1. Farklı depolama yöntemlerinde elde edilebilecek hidrojen miktarı<br />

ve enerji yoğunluk değerleri...................................................................... 7<br />

Çizelge 1.2 Bazı önemli inter<strong>metal</strong>ik bileşiklerin protatipi ve yapısı ........................ 31<br />

Çizelge 4.1. Reaktör 1 YSA modeli nöron 1–6 için LM algoritması denklem<br />

sabitleri.................................................................................................... 94<br />

Çizelge 4.2. Reaktör 2 YSA modeli nöron 1–6 için LM algoritması denklem<br />

sabitleri.................................................................................................... 94<br />

Çizelge 4.3. Reaktör 1 YSA modeli nöron 7–12 için LM algoritması denklem<br />

sabitleri.................................................................................................... 95<br />

Çizelge 4.4. Reaktör 2 YSA modeli nöron 7–13 için LM algoritması denklem<br />

sabitleri.................................................................................................... 95<br />

Çizelge 4.5. Normalizasyon değerleri ........................................................................ 96<br />

Çizelge 4.6. Đstatiksel veriler ...................................................................................... 96<br />

x


Ab<br />

Ac<br />

At<br />

cH<br />

DH<br />

Asal yüzey alanı (m 2 )<br />

Kesit alanı (m 2 )<br />

Toplam yüzey alanı (m 2 )<br />

SĐMGELER VE KISALTMALAR DĐZĐNĐ<br />

Hidrojen konsantrasyonu<br />

Difüzyon katsayısını (cm 2 /s)<br />

F Faraday sabiti<br />

h Isı taşınım katsayısı (W/m 2 K)<br />

H Reaksiyon entalpisi (kJ/kg)<br />

H Silindir yüksekliği (m)<br />

H/M Hidrojenin <strong>metal</strong>e oranı<br />

L Kanat uzunluğu (m)<br />

m Kütle (kg)<br />

N Kanat sayısı<br />

n Mol sayısı (mol)<br />

NA<br />

Avagadro sayısı (6.22x10 -23 )<br />

P Basınç (bar)<br />

PA<br />

PD<br />

qt<br />

qwo<br />

Ru<br />

Absorpsiyon basıncı (bar)<br />

Desorpsiyon basıncı (bar)<br />

Kanatlı durumda geçen ısı (W/K)<br />

Kanatsız durumda geçen ısı (W/K)<br />

Evrensel gaz sabiti (8.314 J/molK)<br />

S Entropi (kJ/K)<br />

t Kanat et kalınlığı (m)<br />

T Sıcaklık (°C)<br />

t Zaman (s)<br />

V Hacim (m 3 )<br />

α Küre yarıçapı (m)<br />

εf<br />

ηf<br />

θb<br />

Kanat etkenliği<br />

Kanat verimi<br />

Sıcaklık farkı (°C)<br />

xi


1. GĐRĐŞ<br />

Ekonomik kalkınmanın temel öğelerinden biri olan enerji, insanlığın vazgeçilmez<br />

gereksinimlerinden biridir. Dünya nüfusu ve endüstriyel gelişmelere paralel olarak<br />

enerji gereksinimi giderek artmakta buna karşın fosil enerji kaynaklarının rezervleri<br />

hızla tükenmektedir. En son istatiksel değerlendirmelere göre; dünya enerji<br />

ihtiyacının %38.5’ini karşılayan petrolün 41 yıl, %23.7’sini karşılayan doğal gazın<br />

62 yıl, %24.7’sini karşılayan kömürün ise 230 yıl rezerv kullanım süresi<br />

bulunmaktadır. 1900 yılında nüfusu 1.6 milyar, birincil enerji tüketimi yaklaşık 1000<br />

Mtep olan dünyamızda 1997 yılında nüfus 6.5 milyara ulaşmış, birincil ticari enerji<br />

tüketimi 11700 Mtep düzeyine çıkmıştır (www.iea.org, 2011). Böylece bir yüzyıl<br />

içinde dünyanın birincil enerji tüketimi 8 katın üzerinde artış göstermiştir.<br />

Dünyada enerji talebinin karşılanmasında ana kaynakların fosil yakıtlar olması, fosil<br />

yakıtların yanma reaksiyonu ile değerlendirilmesi ve bu reaksiyonda karbondioksit<br />

(CO2), karbon monoksit (CO), kükürt dioksit (SO2) ve azot oksit (NOx) gibi zararlı<br />

emisyonların ortaya çıkması, çok büyük çevre sorunları yol açmaktadır. Bugün<br />

dünyanın en önemli çevre sorunu olan küresel ısınmanın ana nedeni, artan CO2<br />

emisyonu ile atmosferin sera etkisinin güçlenmesidir. Birincil enerji kaynaklarının<br />

rezervlerinin kısıtlı olması, yakıt fiyat artışları, nüfus artışı, endüstrileşme, ulusal<br />

kaynaklarının değerlendirilmesi zorunluluğu, 21. Yüzyılın sosyo-ekonomik<br />

yapılanması, mevcut yakıtların çevre üzerindeki olumsuz etkileri (sera etkisi, küresel<br />

ısınma, iklim değişiklikleri, yağış anormallikleri, asit yağmurları, sağlık problemleri<br />

gibi), yeni enerji teknolojileri kapsamında, depolanabilir, yüksek kalorifik değeri<br />

taşınabilir bir enerji kaynağı olan hidrojen enerjisini ön plana çıkarmıştır.<br />

Hidrojen enerjisi önümüzdeki yüzyılın en önemli enerji kaynaklarından biri olarak<br />

kabul edilmektedir. Bu enerji, sudan elde edilebilmekte ve yüksek verimlilikle, çevre<br />

üzerinde hiçbir olumsuz etki yaratmadan yararlı bir enerjiye dönüştürülebilmektedir<br />

(Altan ve Yörükoğulları, 1997). Dünyanın enerji sorununu çözmek için kullanılacak<br />

hidrojen enerjisi ile milyarlarca yıl yetecek enerji devamlı olarak üretilebilecektir.<br />

Hidrojen gazı bir enerji olarak avantajlarından en önemlisi toksik etki içermeyen ve<br />

1


korozyona neden olmayan bir gaz oluşudur. Bu yüzden gerekli önlemler alındığında<br />

hidrojen enerjisinin kullanılması tehlikesiz ve çok basittir. Kirliliğe neden olmayan<br />

hidrojen yakıt olarak yakılması sonucunda atık olarak su buharı ortaya çıkmaktadır<br />

(Pottier and Bailleux, 1986). Atık madde olarak su buharı oluşması doğal çevreye<br />

zarar vermemesinin yanında, diğer enerji atıklarının verdiği zararların giderilmesine<br />

yardımcı olabilmekte ve teknolojinin birçok alanında hidrojen atıklarından etkin bir<br />

şekilde doğal geri dönüşümlü olarak yararlanılabilmektedir (Altan ve Yörükoğulları,<br />

1997).<br />

Bir enerji kaynağı olarak hidrojenin sahip olduğu en önemli özellik depolanabilir<br />

olmasıdır. Ancak bilinen en hafif gaz olması nedeniyle (yoğunluğu 0.0838 kg/m 3 )<br />

depolanmasında bazı problemler ortaya çıkmaktadır. Geleceğin enerjisi olan<br />

hidrojenin günümüzde yaygın olarak kullanımı için depolanmasında verimli ve etkin<br />

yöntemlerin geliştirilmesi gerekir. Hidrojeni depolamada kullanılan birçok yöntem<br />

vardır. Hidrojen sıkıştırılarak gaz halinde ya da çok düşük sıcaklıklara soğutulup<br />

sıvılaştırılarak sıvı halde depolanabilmektedir. Fakat yapılan araştırmalar, en etkili<br />

depolama yönteminin katı olarak depolama olduğunu göstermiştir. Hidrojen, fiziksel<br />

olarak karbon nanotüplerde, kimyasal olarak da <strong>hidrür</strong>lerde katı halde<br />

depolanabilmektedir. Depolama yöntemleri içinde yüksek depolama kapasiteleri,<br />

uygun çalışma ortamları ve düşük basınçta çalıştıkları için daha güvenli olmaları<br />

nedeniyle en etkili yöntemin <strong>hidrür</strong> şeklinde depolama olduğunu göstermiştir<br />

(Barkhordarian et al., 2006). Metal <strong>hidrür</strong> olarak depolamada hidrojen <strong>metal</strong>ler arası<br />

boşluklarda depolanmaktadır. Hidrür bileşikleri belli bir sıcaklıkta hidrojenle<br />

reaksiyona girerek hidrojeni absorbe eden, ısıtıldıklarında ise depoladıkları hidrojeni<br />

serbest bırakan bileşiklerdir. Hidrojen <strong>hidrür</strong> formunda yükseltilmiş sıcaklıklarda<br />

birçok geçiş <strong>metal</strong>leri ve alaşımlarıyla reaksiyona girer. Bu alaşımlardan<br />

elektropozitif olanlar, örneğin Sc, Yt, lantanitler, aktinitler ve Ti, Va grubu<br />

elementler oldukça reaktif <strong>metal</strong> <strong>hidrür</strong>lerdir (Züttel, 2003).<br />

Yapılan bu çalışmanın amacı, daha yüksek kapasitede hidrojen depolayabilen,<br />

depoladığı bu hidrojenin daha uygun sıcaklıklarda serbest bırakabilen, hafif ve<br />

ekonomik olan yeni <strong>hidrür</strong> alaşımlarının araştırılmasıdır. Depolama <strong>metal</strong>leri ile ilgili<br />

2


araştırılan karakteristikler ise şarj/deşarj edilen hidrojen miktarı, <strong>metal</strong> <strong>hidrür</strong>ün<br />

termal stabilitesi, şarj/deşarj kinetiği, termodinamik ve termofiziksel özellikler,<br />

kristal yapı ve segregasyon veya karbonizasyon gibi yüzey işlemleri olarak<br />

sıralanabilir. Özellikle reaksiyon hızını önemli ölçüde etkileyen parametrelerin<br />

başında alaşım özellikleri ve reaktör tasarımı gelmektedir. Bu çalışmada LaNi-Al<br />

alaşımlarının hidrojen depolama özellikleri deneysel ve teoriksel olarak<br />

incelenmiştir.<br />

Bu amaçla Reaktör 1 ve Reaktör 2 olarak adlandırılan kanatçıksız ve kanatçıklı<br />

olmak üzere iki adet reaktör tasarlanmıştır ve imal edilmiştir. Her iki reaktörde 2, 4,<br />

6 ve 8 bar olarak hidrojen şarj basıncı farklı değerlere ayarlanmış, reaktörlerin z-<br />

ekseni boyunca 5, 20, 40 ve 60 mm noktalarına yerleştirilen termokupllar<br />

kullanılarak reaktör içerisindeki sıcaklık dağılımları belirlenmiştir. Ayrıca şarj<br />

basıncı 4, 6 ve 8 bar arasında değiştirilerek reaktörlerde depolanan hidrojen<br />

kütlesine; şarj basıncının, sıcaklığın ve reaktör tasarımının etkisi belirlenmiştir. Diğer<br />

yandan yapay sinir ağları kullanılarak oluşturulan model ile depolanan hidrojen<br />

kütlesi tahmin edilmiş ve deneysel verileri karşılaştırılarak oluşturulan modelin<br />

geçerliliği ortaya konmuştur.<br />

1.1. Hidrojen Üretimi<br />

Hidrojen doğal bir yakıt olmayıp, birincil enerji kaynaklarından yararlanılarak su,<br />

fosil yakıtlar ve biokütle gibi değişik hammaddelerden üretilebilen sentetik bir<br />

yakıttır. Üretim kaynakları son derece bol ve çeşitlidir. Bunlar arasında su, hava,<br />

kömür ve doğal gaz sayılabilir. Ancak, sayılan bu kaynaklardan kömür ve doğal gaz<br />

fosil yakıt olup, sınırlı rezerve sahiptir. Ayrıca fosil yakıtların hem birincil enerji<br />

kaynağı olarak kullanılması hem de hidrojen üretiminde kullanılması çok büyük<br />

çevre sorunlarına yol açmaktadır. Bu nedenle, hidrojenin temiz enerji kaynakları ile<br />

sudan üretimi en doğru seçim olacaktır. Her türlü birincil enerji yardımıyla üretilen<br />

hidrojen, günümüzde suni gübreden, nebati yağlara, roket yakıtlarına kadar çeşitli<br />

alanlarda kullanılmaktadır. Günümüzde çeşitli uygulamalar için yaklaşık 600 milyar<br />

metreküp hidrojen üretilmektedir. Hidrojen üretimi için en önemli yöntem suyun<br />

3


elektrolizle oksijen ve hidrojene ayrıştırılmasıdır. Diğer yandan hidrojen üretiminde<br />

fosil yakıtların yanında güneş, rüzgâr, dalga enerjileri, jeotermal enerji ve biokütle<br />

gibi birincil enerji kaynaklarından da yararlanılmaktadır (Türe, 2001).<br />

Hidrojenin üretilmesi aşamasında Şekil 1.1’de gösterilen buhar iyileştirme, atık<br />

gazların saflaştırılması, elektroliz, foto süreçler, termokimyasal süreçler, radyoliz<br />

gibi alternatif birçok hidrojen üretim teknolojileri mevcuttur. Üretilen hidrojen boru<br />

hatları veya tankerler ile büyük mesafelere taşınabilir. Hidrojen diğer yakıtlara göre<br />

pahalı olmasına rağmen uzun dönemde teknolojik ilerlemelerle enerji kullanımında<br />

önemli rol oynayacaktır. Pazarın bölgesine ve boyutuna bağlı olarak hidrojenin kg<br />

başına maliyeti 3.52–10.5 TL arasındadır. Ancak bu maliyet göreceli olup, hidrojen<br />

çağına adım atılmakla hızlı düşüşü beklenmektedir. Çevresel zararlar ve yüksek<br />

kullanma verimi dikkate alındığında solar hidrojen enerji sistemleri en düşük etkin<br />

maliyete sahiptir.<br />

UYGULAMA DAĞITIM ÜRETĐM<br />

Şekil 1.1. Hidrojen üretim zinciri (Gosselink, 2002)<br />

4


1.1.1. Fosil yakıt yöntemleri<br />

Günümüzde sanayide kullanılan hidrojen büyük miktarlarda, doğal gaz, petrol<br />

ürünleri veya kömür gibi fosil yakıtlardan elde edilmektedir. En çok kullanılan<br />

yöntemler, doğal gazın katalitik buhar ıslahı, petrolün kısmi oksidasyonu, buhar<br />

demir işlemi ve kömür gazlaştırılması şeklindedir. Bunlardan başka, temel amacı<br />

hidrojen üretimi olmamakla birlikte diğer sanayi maddelerinin üretimi sırasında, yan<br />

ürün olarak hidrojen elde edilmektedir. Bu yöntemler arasında, klor-alkaliden karşıt<br />

klor üretimi, ham petrolün rafineri işleminde hafif gazların üretimi, kok fırınlarında<br />

kömürden kok üretimi ve margarin sanayisinde kimyasal hidrojenerasyon işlemleri<br />

sayılabilir (Türe, 2001).<br />

1.1.2. Elektroliz<br />

Suyun doğru akım kullanılarak hidrojen ve oksijenlerine ayrılması elektroliz olarak<br />

tanımlanmaktadır. Bu yöntem hidrojen üretimi için en basit yöntem olarak<br />

bilinmektedir. Đlke olarak, bir elektroliz hücresi içinde, genelde düzlem bir <strong>metal</strong><br />

veya karbon plakalar olan, iki elektrot ve bunların içine daldırıldığı, elektrolit olarak<br />

adlandırılan iletken bir sıvı bulunmaktadır. Doğru akım kaynağı bu elektrotlara<br />

bağlandığında akım iletken sıvı içinde, pozitif elektrottan negatif elektroda doğru<br />

akacaktır. Bunun sonucu olarak da, elektrolit içindeki su, katottan çıkan hidrojen ve<br />

anottan çıkan oksijene ayrışacaktır. Burada yalnız suyun ayrışmasına karşılık, su iyi<br />

bir iletken olmadığı için elektrolit in içine iletkenliği artırıcı olarak genelde potasyum<br />

hidroksit gibi bir madde eklenmektedir. Suyun elektrolizi için, normal basınç ve<br />

sıcaklıkta, ideal olarak 1.23 Volt yeterlidir. Ancak bu değerlerde tepkimenin oldukça<br />

yavaş olması ve oluşan hidrojenin oksijenden ayrıştırılması gibi nedenlerden dolayı,<br />

elektroliz işleminde daha yüksek gerilimler kullanılmaktadır. Hidrojen üretim hızı,<br />

gerçek akım şiddeti ile orantılı olduğundan, ekonomik nedenlerle yüksek akım<br />

yoğunlukları tercih edilmektedir. Bundan dolayı pratikte suyun ayrıştırılması için<br />

hücre başına uygulanan gerilim genelde 2.0 Volt’dur. Teorik olarak, her metreküp<br />

oksijen için 2.8 kWh elektrik enerjisi yeterli olmakla birlikte, yukarıda belirtilen<br />

nedenlerden dolayı pratikte kullanılan elektrik enerjisi miktarı bir metreküp hidrojen<br />

5


üretimi için 3.9–4.6 kWh arasında değişmektedir. Buna göre elektroliz işleminin<br />

verimi yaklaşık %70 civarında olmaktadır. Ancak, son yıllarda bu alanda yapılan<br />

çalışmalar ve gelişen teknoloji sayesinde yaklaşık %90 verim elde edilmiştir. Pratikte<br />

kullanılan elektroliz hücrelerinde, nikel kaplı çelik elektrotlar kullanılmaktadır (Türe,<br />

2001).<br />

1.1.3. Isıl kimyasal yöntem<br />

Fosil veya nükleer yakıtlardan elde edilen birincil enerjilerin elektroliz yolu ile<br />

hidrojene dönüştürülmesinin verimi, ilk başlarda bu yakıtlardan elde edilecek<br />

elektrik enerjisinin verimine bağlıdır. Elektrik üretim verimi, modern fosil yakıt<br />

santralleri için yaklaşık %38 ve nükleer tesisler için %32 seviyesindedir. Elektroliz<br />

hücresinin ticari olarak %80 verim de çalıştığı düşünüldüğünde, fosil yakıtlardan<br />

elektroliz yoluyla hidrojen elde etmede toplam verim %25–30 olmaktadır. Elektrik<br />

üretimi sırasında oluşan ısı enerjisi, suyun ayrıştırılması için kullanıldığında, daha<br />

yüksek verim elde etmek mümkündür. Ancak, suyun ısı enerjisi ile ayrıştırılması için<br />

en az 2500 °C’lik bir sıcaklık gerekmektedir. Burada, tek basamakta termo-kimyasal<br />

işlem yerine, birkaç basamaklı işlemler önerilmektedir. Bu alanda yapılan çalışmalar<br />

sonucu, çok basamaklı ısıl kimyasal işlemlerde gerekli sıcaklık 950 °C’ye kadar<br />

indirilmiş, toplam verim ise %50 olarak bulunmuştur. Isıl-kimyasal yöntem<br />

üzerindeki çalışmalar yoğun bir şekilde sürmektedir (Türe, 2001).<br />

1.2. Hidrojenin Depolanması<br />

Hidrojenin yaygın şekilde kullanımında başlıca engel depolanmasıdır. Hidrojen gaz<br />

formunda oda sıcaklığı ve basıncında aynı eşdeğer enerji miktarına sahip bir gazdan<br />

3000 kat daha fazla yer kaplar. Bu nedenle de hidrojenin araçlarda ve mobil<br />

uygulamalarda kullanımı için sıkıştırma, sıvılaştırma veya diğer yöntemlerinde<br />

kullanılması gereklidir. Depolamada dört ana teknik kullanılmaktadır. Bunlar;<br />

sıkıştırılmış gaz, kriyojenik sıvı, <strong>metal</strong> <strong>hidrür</strong> ve karbon adsorpsiyonudur. Kısa<br />

dönemde en uygulanabilir olanları kriyojenik sıvı ve <strong>metal</strong> <strong>hidrür</strong>lerde depolamadır.<br />

Metal <strong>hidrür</strong> yöntemi gelişmiş bir yöntem olsa da rekabet edebilir olması için daha<br />

6


fazla araştırma yapılması gereklidir. Karbon adsorpsiyonu ise henüz olgunlaşmış bir<br />

teknik değildir, ancak araştırma-geliştirme çalışmalarının sonunda hedefler<br />

gerçekleştirilirse uygulanabilir yöntem olarak görülmektedir. Hidrojenin son<br />

kullanımda depolama teknikleri her bir uygulama için farklıdır. Çizelge 1.1’de bazı<br />

depolama teknikleriyle depolanabilecek hidrojen miktarları ve enerji yoğunlukları<br />

verilmiştir.<br />

Çizelge 1.1. Farklı depolama yöntemlerinde elde edilebilecek hidrojen miktarı ve<br />

enerji yoğunluk değerleri (Selvam et al. 1986)<br />

Depolama Hidrojen miktarı Hacimce yoğunluk* Enerji Yoğunluğu<br />

Ortamı (ağırlıkça %) (H/l ) (x10 25 ) MJ/kg MJ/l<br />

Gaz halde 100 0.50 141.9 1.20<br />

Sıvı H2 (-253 °C) 100 4.20 141.9 9.92<br />

MgH2 7.65 6.70 9.92 14.32<br />

VH2 2.10 11.4 - -<br />

Mg2NiH4 3.60 5.90 4.48 11.49<br />

TiFeH1,95 1.95 5.50 2.47 13.56<br />

LaNi5H6,7 1.50 7.60 1.94 12.77<br />

NaAlH4 7.40 - - 8.25<br />

NaBH4 (Katı) 10.6 6.80 - -<br />

NaBH4 -20 sol. 4.40 - 44.0 -<br />

NaBH4 -35 sol. 7.70 - 77.0 -<br />

Nanotüpler 1–10 - - -<br />

Benzin - - 47.27 6.6–9.9<br />

Metanol - - 22.69 5.9–8.9<br />

* Bu değerlere tank ağırlığı dâhil edilmemiştir.<br />

1.2.1. Sıkıştırılmış gaz olarak depolama<br />

Bu depolama oda sıcaklığında yüksek basınca dayanıklı tankta yapılmaktadır.<br />

Sıkıştırılmış gaz depolama tankı Şekil 1.2’de gösterilmiştir. Sıkıştırılmış gaz<br />

depolamada tankın ağırlığına dolayısıyla tankın tipine bağlı olarak ağırlıkça %1–7<br />

hidrojen depolanmaktadır. Daha hafif, dayanıklı ve ağırlıkça daha fazla hidrojen<br />

depolayabilen tanklar daha pahalıdır. Doldurma istasyonunda hidrojen gazının<br />

sıkıştırılması için yakıtın enerji içeriğinin %20’si kadarı harcanır.<br />

Sıkıştırılmış gaz olarak basınçlı tanklarda depolama yöntemi iyi bilinen bir depolama<br />

yöntemidir. Hidrojen yüksek basınç tanklarının içine sıkıştırılır. Bu süreci<br />

7


gerçekleştirmek için enerji gereklidir ve sıkıştırılan gazın doldurduğu bu hacim<br />

genellikle oldukça geniştir. Bu hidrojenin geleneksel benzin tanklarına göre düşük<br />

enerji yoğunluğuna sahip olması ile sonuçlanır. Bir hidrojen gaz tankı bir benzin<br />

tankının depoladığı enerjiye eşit miktarda enerji ihtiva ettiğinde benzin tankından<br />

3000 kez daha büyük olacaktır (www.fuelcellstore.com, 2010).<br />

Hidrojen günümüzde genellikle 50 litrelik silindirik depolarda 200–250 bar basınç<br />

altında depolamaktadır. Ancak gaz olarak depolama basıncı 600–700 bar’a kadar<br />

çıkabilmektedir. 50 litrelik tanklarda bile depolansalar hidrojen çok hafif olduğu için<br />

hacimsel enerji yoğunluğu çok düşüktür. Diğer taraftan yüksek basınçtan dolayı<br />

depolama tankları ağır olmaktadır (www.e-sources.com, 2010).<br />

Hidrojen çok yüksek basınç ve düşük sıcaklıklarda ideal gaz özelliği göstermektedir.<br />

Dolayısıyla belli bir basınç ve sıcaklıktaki molekül sayısı ve kütlesi ideal gaz<br />

denkleminden hesaplanabilir.<br />

P V<br />

n = ⇒ m = n M<br />

(1.1)<br />

R T<br />

u<br />

Bu denklemde, n hidrojenin mol sayısı, M; molekül ağırlığı, R;evrensel gaz sabitidir.<br />

Hidrojenden elde edilecek enerji;<br />

E = m ΔH<br />

(1.2)<br />

Bağıntısı ile hesaplanabilir. Bu bağıntıdaki Δ H hidrojen gazının yanma ısısını<br />

vermektedir ve bu değer hidrojen için 140 MJ/kg’dır. Hidrojenin hacimsel enerji<br />

yoğunluğu ise aşağıdaki denklemden hesaplanabilir.<br />

W<br />

hacimsel<br />

E ΔH<br />

P P ΔH<br />

= = =<br />

(1.3)<br />

V n R T R T<br />

Bağıntısı ile hesaplanır. Bu bağıntıda görüldüğü gibi enerji yoğunluğu basınçla<br />

artmaktadır. Fakat basınçtaki artış depo malzemesinin dayanımı ile sınırlıdır.<br />

8<br />

u


Uygulamada basınçlı depo malzemesi olarak ostenitik çelik ve alüminyum alaşımları<br />

kullanılmaktadır. Fakat bu depoların dezavantajı çok ağır olmalarıdır. Depolanan<br />

hidrojenin tüm deponun ağırlığına göre oranı %2–3 civarındadır. Depoların bu<br />

dezavantajları kompozit malzeme kullanılarak giderilebilir. Bu depolarda hidrojen ile<br />

temas eden ince bir <strong>metal</strong> tabakası ve bunu takiben karbon fiberi ile güçlendirilmiş<br />

polimer zarf gelmektedir. Böyle bir kapta hidrojen ağırlık oranı %5’e çıkmaktadır.<br />

Daha gelişmiş depolarda <strong>metal</strong> zarf yerine çok tabakalı polimer kullanılmaktadır.<br />

Şekil 1.2. Sıkıştırılmış hidrojen gazı depolama tankı (kompozit silindir)<br />

(Sirosh, 2002)<br />

1.2.2. Sıvı olarak depolama<br />

Bu teknikte hidrojen atmosfer basıncında -253 °C’de Şekil 1.3’te gösterildiği gibi<br />

oldukça iyi yalıtılmış tanklarda depolanmaktadır. Hidrojen sıvı şekilde olduğu için,<br />

eşdeğer ağırlıktaki gazolinden 3 kat fazla enerji içerir ve eşdeğer enerji içerdiği<br />

durumda da 2.7 kat fazla hacim gerektirir. Bu teknik tank ve yalıtım dâhil ağırlıkça<br />

%16 hidrojen depolar. Ayrıca, sıvılaştırma yakıtın enerji içeriğinin %28’i kadarını<br />

gerektirir. Diğer bir dezavantaj yalıtıma rağmen tanka ısı transferi olmasıdır. Bu ısı<br />

transferi sonucunda hidrojen buharlaşmaktadır. Ancak basınçlı tank kullanılarak bu<br />

problem çözülebilir. Ama bu çözüm de ağırlığı ve boyutu arttırmaktadır.<br />

9


Şekil 1.3. Sıvı hidrojen depolama tankı (Linde gas, 2010)<br />

Hidrojeni sıvılaştırma için harcanan enerji fazla olsa bile, uzay araçları ve<br />

roketlerdeki sıvılaştırma masrafları göz ardı edilmektedir. Sıvı hidrojen büyük<br />

tanklarda depolanmışsa günlük %0.06’sı, eğer küçük tanklarda depolanmışsa günlük<br />

%3’ü buharlaşarak kaybolmaktadır. Bu oranın azaltılması yalıtıma bağlıdır.<br />

1.2.3. Metal <strong>hidrür</strong> esaslı depolama<br />

Metal <strong>hidrür</strong>lerde hidrojen depolama tekniğinde hidrojen granülleri <strong>metal</strong>lerin<br />

atomları arasındaki boşluğa depolanır. Bu amaçla çeşitli <strong>metal</strong>ler kullanılmaktadır.<br />

Depolanan hidrojenin kullanımı için ısı gereklidir. Metal <strong>hidrür</strong> sistemleri güvenilir<br />

ve az yer kaplar, ancak ağırdır ve pahalıdır. Araştırma aşamasında olan<br />

uygulamalarda ağırlıkça yaklaşık %7 oranında hidrojen depolanabilmektedir.<br />

Sıkıştırılmış gaz veya kriyojenik sıvı depolamanın aksine <strong>metal</strong> <strong>hidrür</strong> yeniden<br />

doldurulmada çok az enerji gerektirir. Ancak yakıtın dışarıya salınımı için enerji<br />

harcanır. Düşük sıcaklıkta <strong>metal</strong> <strong>hidrür</strong> depolanmasında bu enerji yakıt hücresinin<br />

veya motorun atık ısısından sağlanabilir. Yüksek sıcaklık <strong>metal</strong> <strong>hidrür</strong> depolaması<br />

10


daha ucuz olmasına rağmen, aracın enerji tüketiminin yarısı <strong>metal</strong>den hidrojeni açığa<br />

çıkarmak için harcanır.<br />

1.2.4. Karbon nanotüplerde depolama<br />

Karbon özellikle yüksek oranda gözenekli çok küçük parçacıklar haline<br />

getirilebilmesi ve karbon atomları ile gaz molekülleri arasında oluşan çekim kuvveti<br />

nedeniyle gaz depolamaya en elverişli maddelerden biridir. Nanotüpler birçok özel<br />

üstün özelliklere sahiptirler. Örneğin elastiklik modülü çelikten beş kat daha fazladır.<br />

Ayrıca tüpün yapısına bağlı olarak bazıları yarı iletken bazıları da iletken olarak<br />

davranırlar. Bu özellikleri dolayısıyla nanotüp kullanarak elektronik cihazları mikro<br />

ve nano boyutlara indirebilmek mümkündür. Hidrojen nanotüplerin içerisine<br />

kimyasal veya fiziksel yollarla depolanmaktadır. Karbon nanotüpleri 1991’de Lijima<br />

tarafından keşfedilmiştir.<br />

Karbon nanotüpler kısaca grafit tabakaların tüp sekline dönüşmüş halidir. Çapları<br />

birkaç nanometre veya 10–20 nanometre mertebesinde, boyları ise mikron<br />

civarındadır. Nanotüpler tek duvarlı olarak üretilebileceği gibi çok-duvarlı tüplerde<br />

üretilebilmektedir. Tek ve çok duvarlı nantüpler ve nanotüplerden oluşan iplik yapı<br />

Şekil 1.4’te verilmiştir. Çeşitli ilavelerle oluşturulan, örneğin alkali ilaveli (Li-K),<br />

nanotüpler de mevcuttur (Hirscher et al., 2002).<br />

Şekil 1.4. Tek ve çok duvarlı karbon nanotüpler (Lee et al., 2002)<br />

11


Karbon nano-tüplerde hidrojen absorpsiyonu oldukça gözenekli süper aktif grafit<br />

yüzeylerde gerçekleşmektedir. Absorbe işlemi karbon atomlarının hidrojen<br />

moleküllerine uyguladığı Van Der Waal’s kuvveti ile gerçekleşmektedir. Bu nedenle<br />

absorbsiyon işlemi kimyasal değil fiziksel bir olaydır. Verilen bir sıcaklıkta absorbe<br />

edilen hidrojen sadece basıncın bir fonksiyonudur. Basınç düşürüldüğünde istenilen<br />

miktarda hidrojen açığa çıkar. Nano yapı basınçlı tanklarda kullanılmakta dolayısıyla<br />

absorbe edilen hidrojen gaz halinde depolanan hidrojene katkı sağlamaktadır. Bu<br />

sistemler ağırlıkça yaklaşık %4 hidrojen depolamaktadır.<br />

1.2.5. Cam mikro kürelerde depolama<br />

Bu yöntemle hidrojen depolamada küçük, içi boş, çapları 25 ile 500 µm arasında<br />

değişen ve duvar kalınlıkları ~1 µm olan cam küreler kullanılmaktadır. Şekil 1.5’te<br />

gösterilen cam mikro küreler yüksek basınç ve 200–400 °C gibi yüksek sıcaklıklarda<br />

hidrojen gazı ile doldurulur. Yüksek sıcaklıkta cam duvarlar geçirgenleşir ve gaz<br />

mikro kürelerin içine dolar. Cam oda sıcaklığına soğutulduğunda, hidrojen kürelerin<br />

içine hapsolur. Hidrojen kullanılacağı zaman küreler tekrar ısıtılarak hidrojen gazının<br />

serbest bırakılması sağlanır. Cam kürelerin depolama kapasitesi 200–490 bar basınç<br />

altında %5–6 civarındadır (Kaya, 2005).<br />

Şekil 1.5. Mikro cam kürelerde hidrojen depolama (Teitel, 1981)<br />

12


1.2.6. Bor esaslı depolama<br />

Sodyum bor <strong>hidrür</strong> esaslı depolamada, alanatlar ve katı haldeki depolama<br />

yöntemlerinden farklı olarak sıvı halde kullanımı esas almaktadır. Çözelti halinde,<br />

sodyum bor <strong>hidrür</strong>, aşağıdaki reaksiyona göre hidrojenini vermekte ve sodyum meta<br />

borata dönüşmektedir. Bu dönüşüm Şekil 1.6’da gösterilmiştir.<br />

NaBH 4<br />

2<br />

2<br />

2<br />

( s)<br />

+ 2H<br />

O → 4H<br />

+ NaBO ( katalizör)<br />

(1.4)<br />

Denklem (1.4)’te görüleceği üzere reaksiyon sonucu açığa çıkan hidrojen miktarı<br />

<strong>hidrür</strong> şeklinde bağlı olan hidrojenin iki katıdır. NaBH4’de mevcuda eşit miktarda<br />

hidrojen suyun parçalanması ile açığa çıkmaktadır. Denklemde verilen reaksiyon<br />

ekzotermiktir. Bunun bir sonucu sistemden elde edilen hidrojenin nemli olmasıdır.<br />

Şekil 1.6. Sodyum bor <strong>hidrür</strong> esaslı hidrojen depolama (Güvendiren vd. 2002)<br />

Sodyum bor <strong>hidrür</strong>de hidrojen depolamanın en önemli üstünlüğü depolanan<br />

hidrojenin oda sıcaklığında geri alınabilmesi ve geri alımın katalizör yardımı ile<br />

kolaylıkla kontrol edilebilmesidir. Nitekim sıvı halde çözelti alevle temas halinde<br />

bile güvenli olmakta, ancak katalizörün çözeltiyle teması durumunda hidrojen çıkışı<br />

sağlanmaktadır. Sistemin en büyük dezavantajı açığa çıkan NaBO2’nin başka bir<br />

ortamda NaBH4’e geri dönüştürülmesi gerekliliğidir. Şekil 1.7’de bor esaslı hidrojen<br />

depolama sistemi verilmiştir.<br />

13


1.2.7. Mağaralarda depolama<br />

Şekil 1.7. Bor esaslı hidrojen depolama<br />

Hidrojen gazını depolamanın en ucuz yöntemi, doğalgaza benzer şekilde, yeraltında,<br />

tükenmiş petrol veya doğal gaz rezervuarlarında depolamaktır. Maliyeti biraz yüksek<br />

olan diğer bir depolama şekli ise, hidrojeni maden ocaklarındaki mağaralarda<br />

saklamaktır. Örneğin Almanya’nın Kiel şehrinde 1971’den beri yerin 1330 m<br />

altındaki bir mağarada hidrojen depolanmaktadır. Ancak mağaralarda saklanan<br />

hidrojenin yılda %1–3 miktarı sızıntı nedeniyle kaybolmaktadır.<br />

1.3. Metal Hidrürler<br />

Metal <strong>hidrür</strong>ler, kafes kusurları bulunan <strong>metal</strong> atomları ve kafes boşlukları ara<br />

bölgeleri içine yerleşmiş hidrojen atomlarından oluşmuştur. Ara bölge boşluk<br />

şeklinde ya da çizgisel hata şeklinde olabilmektedir. Çizgi şeklinde hata olması<br />

durumunda, hidrojen atomlarının dizilimi hata boyunca toplanması şeklinde olur.<br />

Öyle ki, özellikle iki komşu atom moleküler hidrojen formunda yeniden birleşirse<br />

kafes gerilimini arttırmaktadır. Hidrojen adsorpsiyonu esnasında kafes boyutunu<br />

arttırmaktadır (Carter and Cornish, 2001; Fisher et al., 1988). Hidrürleşmede iki<br />

14


farklı yol vardır, bunlardan biricisi kimyasal tutunma ikincisi ise suyun<br />

elektrokimyasal ayrılmasıdır. Bu reaksiyonlar sırasıyla şöyledir (Sandrock, 1999);<br />

ve<br />

x<br />

2<br />

x<br />

+ H 2 ↔ MH + Q<br />

(1.5)<br />

2<br />

M x<br />

x<br />

2<br />

−<br />

M + H 2 O + e−<br />

↔ MH x + OH<br />

(1.6)<br />

Denklemlerde M <strong>metal</strong>i ve Q ekzotermik reaksiyonla oluşan ısıyı göstermektedir.<br />

Elektrokimyasal ayrılmada paladyum gibi bir katalizör olması gerekir.<br />

Hidrojenin kimyasal tutunması işlemi şematik olarak Şekil 1.8a’da verilmiştir.<br />

Şekilde gösterildiği gibi, moleküler hidrojen yüzeye yakın minimum sığ bir<br />

potansiyele ulaşır ve atomik hidrojen yüzeyden daha derinde ama hemen hemen<br />

yüzeye yakın bir bölgededir. Metal kafes içinde hidrojen <strong>metal</strong> kafes çatlak bölgeleri<br />

içinde periyodik potansiyel minimuma sahiptir. Bu durum aşağıda açıklanmış ve<br />

Şekil 1.8b’de verilmiştir. Hidrojen molekülleri <strong>metal</strong> yüzeylere yaklaşırken, zayıf<br />

van Der Waal’s kuvvetlerinin etkisiyle birbirine yaklaşırlar. Molekül zp mesafesinde<br />

potansiyel kuyusuna Ep noktasına ulaşır ve çok geniş kuvvetler moleküler formda<br />

yüzeye yaklaşır. Böylece, hidrojen moleküllerinin ayrılma enerjisi kemisorpsiyon<br />

enerjisi tarafından aşılır. Daha sonra hidrojen molekülü ayrılır ve bireysel hidrojen<br />

atomları kemisorptif kuvvetlerin etkisiyle yüzeye doğru çekilir ve ECH potansiyeline<br />

ulaşır. Bu noktadan sonra çevre havası ısıl enerjisi bile hidrojen atomlarının<br />

titreşimsel genişliğini arttırmaya yeterlidir, böylece hidrojen <strong>metal</strong> yüzeylere<br />

ulaşabilir (Carter and Cornish, 2001).<br />

Metal ve hidrojen genellikle α-fazı ve β-fazı olmak üzere iki farklı <strong>hidrür</strong><br />

formundadır. α-fazında sadece biraz hidrojen absorbe edilmiştir ve β-fazında <strong>hidrür</strong><br />

tamamen oluşmuştur. Bunlara örnek olarak Mg2Ni formu <strong>hidrür</strong>lerde Mg2NiH0.3 ve<br />

Mg2NiH4 gösterilebilir.<br />

15<br />

x<br />

2


Şekil 1.8. a) Metal üzerinde hidrojen kemisorpsiyonu b) Moleküler potansiyel<br />

kuyusu ve atomik hidrojen (Christmann, 1981; Schlapbach, 1988).<br />

Başlangıç şarjında α-fazı oluşmaya başlar ve bundan sonra şarj ve deşarj edilmiş<br />

<strong>hidrür</strong> genellikle faz dönüşümü altındadır. Şöyle ki (Cui et al., 2001);<br />

Mg NiH + 3.<br />

7H<br />

↔ Mg NiH<br />

(1.7)<br />

2<br />

0.<br />

3<br />

Faz dönüşümü şeması Şekil 1.9’da verilmiştir. Şarj esnasında, hidrojen partikül<br />

yüzeyinden β-fazında faz-geçişi ara yüzeyine doğru dağıtılır ve bu form klasik β-fazı<br />

<strong>hidrür</strong>dür. Deşarj esnasında hidrojen faz değişim ara yüzeyinden α-fazına doğru<br />

partikül yüzeyinde dağılır burada hidrojen moleküler formunda yeniden birleşir.<br />

Modelde çözümü basitleştirmek amacıyla alaşım parçacıklarının küresel geometriye<br />

sahip oldukları, α ve β fazlarının her birinde difüzyonun eşit olduğu varsayılmıştır.<br />

Bu varsayımla partükül içerisine difüzyon olan hidrojen küresel koordinatlarda<br />

aşağıdaki denklemde gösterilen Fick ikinci yasası kullanılarak modellenebilir;<br />

∂<br />

( r,<br />

t)<br />

= D<br />

∂t<br />

⎡ 2<br />

∂ cH<br />

( r,<br />

t)<br />

2 ∂c<br />

⎢ +<br />

2<br />

⎣ ∂r<br />

r<br />

cH H<br />

H<br />

16<br />

2<br />

4<br />

( r,<br />

t)<br />

⎤<br />

⎥<br />

∂r<br />

⎦<br />

(1.8)


Burada cH partikül içerisindeki hidrojen konsantrasyonu, t zaman ve r küresel<br />

koordinatı göstermektedir. cH (r0)=cs sınır şartları için denklem<br />

(1.8)’in çözümü;<br />

c<br />

c<br />

H<br />

s<br />

− c<br />

− c<br />

0<br />

0<br />

2α<br />

= 1+<br />

πr<br />

∑<br />

n 1<br />

∞<br />

=<br />

( −1)<br />

n<br />

n<br />

nπr<br />

sin<br />

α<br />

17<br />

2<br />

D Hn<br />

π<br />

− 2<br />

α<br />

2<br />

t<br />

(1.9)<br />

Burada α küre çapı ve DH hidrojenin alaşım partikülleri arasındaki difüzyon<br />

katsayısını (cm 2 /s) göstermektedir.<br />

Deşarj<br />

Şarj<br />

Şekil 1.9. Metal <strong>hidrür</strong>lerde faz dönüşüm şeması (Cui et al., 2001)<br />

Merkezdeki hidrojen konsantrasyonu şöyle bulunabilir;<br />

c<br />

c<br />

H<br />

s<br />

− c<br />

− c<br />

0<br />

0<br />

= 1+<br />

2<br />

∑<br />

n 1<br />

∞<br />

=<br />

D<br />

n<br />

n H<br />

( −1)<br />

− 2<br />

α<br />

2<br />

π<br />

2<br />

t<br />

(1.10)<br />

Mt ve M∞, t sürede küresel partiküllere difüzyonla giren ve çıkan hidrojen miktarları<br />

denklem (1.11) kullanılarak bulunabilir.


M<br />

M<br />

t<br />

∞<br />

6<br />

1−<br />

π<br />

= ∑ ∞<br />

2<br />

n=<br />

1<br />

1<br />

2<br />

n<br />

Böylece difüzyon akışı şöyle ifade edilir;<br />

2<br />

D Hn<br />

π<br />

− 2<br />

α<br />

n=<br />

1<br />

18<br />

2<br />

t<br />

(1.11)<br />

∑ ∞<br />

2 2<br />

dM t 6FD<br />

HM<br />

∞ D Hn<br />

π t<br />

Ι = −F<br />

=<br />

−<br />

(1.12)<br />

2<br />

2<br />

dt α<br />

α<br />

Burada F Faraday sabitidir. Küresel partiküllere difüzyonla giren ve çıkan toplam<br />

hidrojen miktarı M∞ varsayılırsa;<br />

4 3<br />

∞ = ± πα ( c 0 − c )<br />

(1.13)<br />

3<br />

M s<br />

Denklemde ± simgesinde (-) şarj ve (+) deşarjı göstermektedir. Buradan denklem<br />

(1.12) şöyle yazılabilir;<br />

2 2<br />

D Hn<br />

π t<br />

Ι = ± 8 παFD<br />

H ( c 0 − cs<br />

) −<br />

(1.14)<br />

2<br />

α<br />

Uzun süre bölge için, difüzyon yaklaşık olarak aşağıdaki şekilde verilebilir;<br />

∑ ∞<br />

n=<br />

1<br />

2<br />

π D H<br />

Ι = log{<br />

± 8πα<br />

FD<br />

H ( c0<br />

− cs<br />

) } − t<br />

(1.15)<br />

2.<br />

303 α<br />

log 2<br />

Burada, hidrojen difüzyon katsayısı DH; partikül yarıçapı bilinmek koşuluyla<br />

konsantrasyon bilgileri olmadan log Ι eğiminin t zamanına göre değişiminden<br />

bulunabilir.<br />

Bazı araştırmalarda hidrojen konsantrasyonu için α-fazı yayılımının önemli olmadığı<br />

gösterilmekle beraber(Beltowska et al., 2001). Diğer yandan bazı çalışmalar ise β-<br />

fazı difüzyon katsayısının α-fazından daha düşük olduğunu göstermiştir. Buna göre<br />

β-fazında hidrojen difüzyonu daha yavaştır (Cui et al., 2001). Bu kolaylıkla büyük<br />

partiküllerle eksik şarja neden olabilir. Bu durum özellikle geniş partiküllerde


kusurlu şarj işlemine neden olmaktadır. Nano-ölçekli partiküllerle yapılan bir<br />

çalışmada <strong>metal</strong> taneciklerinin boyutu 5-50 µm olduğunda, absorpsiyon ve<br />

desorpsiyon kinetiklerinin büyüklükle geliştirildiği için, ısı iletimini de<br />

geliştirildiğini göstermiştir (Zaluska et al., 2001).<br />

Hidrürleşme kinetikler katalizörlerle geliştirilebilmektedir (Zaluska et al., 2001). Bu<br />

katalizörler sıvı ya da katı fazında olabilir, ancak katalizörün bütün reaksiyonları<br />

etkilemeyeceğinden, depolama miktarının yeterli olması için katalizör miktarının<br />

düşük tutulması gerekir. Şekil 1.10’da LaNi5 hidrojen adsorpsiyonunda nano yapı ve<br />

katalizörlerin etkisi verilmiştir. Burada X ağırlıkça % hidrojen miktarını<br />

göstermektedir.<br />

Şekil 1.10. LaNi5 hidrojen absorpsiyon hızı a) Polikristalize, b) Nano-kristalize ve<br />

c) katalizörlü nanokristalize (Zaluska et al., 2001)<br />

Metal <strong>hidrür</strong>lerin en genel karakteristik yöntem PCT (basınç-konsantrasyon-sıcaklık)<br />

P-C izotermleri grafiğindeki eğridir. Teorik P-C izotermleri α- ve β- fazlarıyla Şekil<br />

1.11’de verilmiştir. Konsantrasyon, genellikle birim <strong>metal</strong> molekülde hidrojen atomu<br />

19


miktarı H/M olarak belirlenir. Metal <strong>hidrür</strong>leri tanımlamak için en uygun parametre<br />

(H/M)mak maksimum hidrojen kapasitesidir.<br />

Şekil 1.11. Metal <strong>hidrür</strong> teorik P-C eğrisi (Schlapbach, 1988)<br />

Şekil 1.12. Gerçek P-C izotermleri diyagramı (Sandrock, 1999)<br />

20


Tersinir kapasite Şekil 1.12’de gösterildiği gibi ∆(H/M) plato genişliği kullanılarak<br />

belirlenir ve <strong>metal</strong> <strong>hidrür</strong>lerin teorik kapasitesini belirlemede kullanılır.<br />

Konsantrasyon özellikle enerji yoğunluğu karşılaştırmalarında toplam kütle içerinde<br />

yüzde kütle olarak ifade edilebilir. Gerçek <strong>metal</strong> <strong>hidrür</strong>ler P-C eğrisinde,<br />

absorpsiyon/desorpsiyon arasında biraz histeresiz gösterirler. Ayrıca plato biraz<br />

eğime sahiptir. Bu özellikler Şekil 1.13’te gösterilmiştir. Termodinamik reaksiyon<br />

denge denkleminde, reaksiyon sabiti K ile şöyle ifade edilebilir (Hangström, 1999);<br />

RT ln K = ΔH<br />

− ΔS<br />

(1.16)<br />

Burada ∆H reaksiyon entalpisi ve ∆S reaksiyon entropisidir. Katı-gaz reaksiyonu için<br />

denge sabitinden gaz basıncına indirgenir. Böylece van’t Hoff denklemi elde<br />

edilmiştir (Hangström, 1999);<br />

ΔH ΔS<br />

ln P = − (1.17)<br />

R T R<br />

lnP’ye karşı 1/T, denge diyagramı (P, T) değerleri ölçeği vant Hoff diyagramı olarak<br />

adlandırılır. Reaksiyon entalpisi denklem (1.17) ile beraber diyagramda açısal<br />

katsayıdan belirlenebilir ve diyagramda pratik uygulamalar için <strong>hidrür</strong><br />

davranışlarının P-T uyumu ile adlandırılır. Teorik van’t Hoff diyagramı genellikle<br />

<strong>metal</strong> <strong>hidrür</strong>lerin gerçek özelliklerini çok iyi şekilde tanımlar. Şematik van’t Hoff<br />

diyagramı Şekil 1.13’te verilmiştir (Sandrock, 1999).<br />

Şekil 1.13. Şematik van’t Hoff grafiği<br />

21


Hidrür oluşumu reaksiyon entalpisi önemli bir niceliktir. Bu miktar genellikle<br />

reaksiyon ekzotermik ise negatiftir ve böylece <strong>hidrür</strong> oluşumu enerjiyi serbest<br />

bırakır. Bu nedenle hidrojenin tekrar <strong>hidrür</strong> tarafından bırakılması için enerji<br />

verilmesi gerekir. Uygulamaların çoğunda çevre havası kullanılmıştır veya en<br />

azından 0–100 °C aralığında bir sıcaklık olmalıdır. Hidrürün çevre havasından ısı<br />

alarak hidrojenin serbest kalması için reaksiyon entalpisinin oldukça küçük olması<br />

gerekir.<br />

1.3.1. Metal <strong>hidrür</strong> termodinamiği<br />

1.3.1.1. Aktivasyon<br />

Metal <strong>hidrür</strong>lerin pratik uygulamaları için çok önemli bir sorun <strong>metal</strong> yüzeylerinin<br />

genellikle her bir <strong>metal</strong> oluşum sürecine bağlı olarak çeşitli kalınlıklarda oksit<br />

tabakası ile kaplı olmasıdır. Endüstriyel anlamda <strong>metal</strong> <strong>hidrür</strong> üretiminde bu olay her<br />

zaman gerçekleşmektedir. Oksit tabakası bir hidrojen bariyeri gibi hareket ederek,<br />

gaz hidrojenin saf <strong>metal</strong>e ulaşmasını engellemektedir. Bu nedenle ilk hidrojen şarjı<br />

genellikle (genel alaşımlar için) oksit tabakası boyunca hidrojeni zorlamak için<br />

yüksek sıcaklık ve basınçta yapılır. Hidrojen geri alınırken (dehidrojenasyon) kafes<br />

özgün boyutuna geri dönerken, hidrojen şarjı (hidrojenasyon) esnasında kafes hacmi<br />

önemli ölçüde artar. Bu genişleme-daralma <strong>metal</strong> partiküllerini kırarak taze <strong>metal</strong><br />

yüzeyini ortaya çıkarır ve partikül boyutunu küçültür (Yamaguchi and Akiba, 1994).<br />

Aktivasyon işlemiyle önemli anizotropik kafes gerilmeleri ve dislokasyonlar<br />

(çizgisel hatalar) ortaya çıkmaktadır (Guegan et al., 1980; Nakamura and Akiba,<br />

2000). LaNi5 alaşımlarında aktivasyonda izotropik ve anizotropik kafes gerilmeleri<br />

<strong>hidrür</strong> fazında katı çözelti fazı gerilmelerinin küçük olduğu durumlarda sadece <strong>hidrür</strong><br />

fazında gerçekleşmektedir. Dahası, dislokasyonların geneli (1/3 tipi Burger<br />

vektörleriyle) α-tipi kenar dislokasyonlarından oluşmaktadır (Yamamoto et al.,<br />

2001).<br />

Mg2Ni için mekanik öğütme işleminden sonra özgül yüzey alanının ve yeni kafes<br />

çatlakları oluşmasından dolayı aktivasyon işleminin kolaylaştığı belirlenmiştir.<br />

22


Ancak La2Mg18Ni alaşımları için aktivasyonda mekanik öğütmenin hiçbir pozitif<br />

katkı sağlamadığı görülmüştür (Chen et al., 1993). Đri taneli malzemelerde<br />

aktivasyon işlemini kolaylaştırmada birçok yöntem test edilmiştir. Örneğin, hidrojen<br />

depolama alaşımları yüzeylerine uygulanan florlama işlemiyle reaksiyon flor içerikli<br />

çözelti ile modifiye edilmiştir (Liu et al., 1992–1996). Yüzey oksitlerini<br />

kaldırıldığında özgül yüzey alanı artar ve yüzeyde daha fazla katalitik alanlar ortaya<br />

çıkararak hidrojen aktivasyon ve absorpsiyon kinetiklerinin geliştirilmesini<br />

sağlamaktadır.<br />

Pratik uygulamalar için tank sistemini dizayn ederken aktivasyonu dikkate almak<br />

önemli bir faktördür. Tankta aktivasyon işlemi geliştirilebilirse, çalışma şartlarından<br />

daha büyük sıcaklık ve basınç aralığında tank dizaynları yapılabilir. Diğer taraftan<br />

tank aktive edilmiş <strong>metal</strong> <strong>hidrür</strong>lerle doldurulması esnasında, hava ile temas sonucu<br />

kirlenme oluşması yönünden aktive olmuş <strong>metal</strong> <strong>hidrür</strong>ler, aktive olmamış <strong>metal</strong><br />

<strong>hidrür</strong>ler göre daha hassastır.<br />

1.3.1.2. Histeresiz<br />

Histeresiz, karmaşık bir olay olup, <strong>hidrür</strong>leşme ve de<strong>hidrür</strong>leşme esnasında<br />

absorpsiyon platosu desorpsiyon platosundan daha yüksek basınçtadır. Bunun sonucu<br />

şekilde (Bkz. Şekil 1.12.) gösterilen histeresiz döngüsü oluşur. Her bir plato iki<br />

termodinamik parametreye karşılık gelmektedir (Flanagan and Oates, 1988).<br />

Histeresiz olayının bir açıklaması <strong>hidrür</strong>leşme ve de<strong>hidrür</strong>leşme işlemlerinde elastik<br />

ve plastik enerjileri bir birine eşit değildir (Balasubramaniam, 1997; Flanagan, et al.,<br />

1987). Diğer bir modelde ise histeresiz koherent gerilimin bir terimi olarak<br />

açıklanmıştır (Schwarz and Khachaturyan, 2006). Bu modelde <strong>metal</strong>-hidrojen sistemi<br />

geniş çatlak kaynakları ve tutarlı zorlanma ile iki faz arasında üretilmiş makroskopik<br />

bariyerle dengede, kısmen açık iki fazlı bir yapı gibi kabul edilmiştir. Bu<br />

makroskopik değer ısıl dalgalanmayla aşılamaz, bu bir metastabil fazdır ve kimyasal<br />

potansiyel engeli aşmak için yeterli değere ulaşıncaya kadar bir bariyer gibi kilitlidir.<br />

Bu makroskopik bariyer, alaşımın niteliğine bağlıdır ve bu yol dönüşüm<br />

başlangıcından bağımsızdır (Schwarz and Khachaturyan, 2006).<br />

23


Tipik <strong>metal</strong> <strong>hidrür</strong>lerdeki kafes genişlemesi %20 ve üzerinde iken burada<br />

de<strong>hidrür</strong>leşme ve <strong>hidrür</strong>leşme fazları arasında kafes genişlemesi bilinenin aksine çok<br />

büyük değildir. Son yıllarda Rabkin ve Skripyunk de<strong>hidrür</strong> matrislerinin ve <strong>hidrür</strong><br />

partiküllerinin bir arada olabildiği sözde-iki fazlı denge kavramını önermişlerdir. Bu<br />

model toz partikül boyutunu küçülterek histeresizi düşürmeyi öngörmektedir.<br />

Nicelik olarak histeresiz serbest enerji kullanılarak ifade edilebilir;<br />

Δ<br />

⎛ P<br />

( histeresiz)<br />

= RT ln ⎜<br />

⎝ P<br />

24<br />

⎞<br />

⎟<br />

⎠<br />

G H2<br />

A<br />

D<br />

(1.18)<br />

Burada PA ve PD sırasıyla absorpsiyon ve desorpsiyon basınçlarıdır. Denklemde PD<br />

ile gösterilen gerçek denge basıncıdır, ancak Flanagan vd. Denge basıncının PA ve PD<br />

arasında bir yerde olduğunu ispatlamışlardır (Flanagan and Oates, 1988). Sandrock<br />

vd. bu parametrenin tek bir malzeme özelliği olmadığını ve örneklerin geçmişine ve<br />

kullanılan test yöntemlerine bağlı olduğunu belirtmişlerdir (Sandrock et al., 1992).<br />

Pratik uygulamalarda depolama tankı servis basıncı üzerine önemli bir etkiye sahip<br />

olduğundan dolayı, histeresiz önemli bir özelliktir.<br />

Aslında hidrojen absorpsiyonu ve desorpsiyonu esnasında tersinir deformasyondan<br />

dolayı malzeme verimi kaybına neden olur. Uygulamaların birçoğunda histeresizin<br />

olabildiğince küçük olması istenmektedir. Histeresizi küçük değerlerde tutma işlemi<br />

elementlerin yer değiştirilmesiyle ve ısıl işlemle gerçekleştirilebilir.<br />

1.3.1.3. Plato eğimi<br />

Gerçek PCI eğrisinin önemli bir özelliği olan plato eğimi genellikle aşağıdaki<br />

eşitlikle gösterilir (Sandrock et al., 1992);<br />

d(ln<br />

P)<br />

Eğim = (1.19)<br />

d(<br />

H / M)


Plato eğimi, plato uzunluğu boyunca sabit olması gerekmez. Plato eğiminin<br />

oluşmasının en önemli nedeni olarak homojen olmayan kompozisyon gösterilebilir<br />

(Park et al., 2004). Erimiş alaşımlar için ısıl işlem plato eğimini azaltabilir.<br />

Genellikle Plato ortasındaki basınç, alaşımı karakterize etmek için kullanılmaktadır.<br />

Histeresizde olduğu gibi plato eğimleri bazı uygulamalarda önemli teknolojik etkiye<br />

sahiptir ve plato eğiminin azaltılması tavsiye edilmektedir.<br />

1.3.1.4. Tersinir kapasite<br />

Pratik uygulamalar için plato genişliği gibi belirlenen tersinir kapasite ∆(H/M)r<br />

önemli bir faktördür. Tersinir kapasite maksimum kapasiteden önemli ölçüde az<br />

olamamalıdır (H/M)mak. Bunlara örnek alarak vanadyum verilebilir. 353 K<br />

sıcaklıkta 10 -6 -10 -7 MPa hidrojen basıncında vanadyum mono<strong>hidrür</strong> VHx (x≈1)<br />

olarak adlandırılır. Daha yüksek basınçlarda (yaklaşık 1 MPa) di<strong>hidrür</strong> VHx (x≈2)<br />

formundadır. Pratik sistemler için mono<strong>hidrür</strong> dehidrojenasyon imkansızdır ve<br />

vanadyumun tersinir kapasitesi VH≈1 ve VH≈2 arasındadır şeklinde söylenebilir.<br />

Hidrojen kapasitesi her bir terimin atomik oranı (H/M) veya yüzde ağırlık olarak<br />

(%ağırlık) belirtilebilir. Atomik oran <strong>metal</strong>ik atom sayısına göre depolanan hidrojen<br />

atomlarının sayısıdır. Yüzde ağırlık, <strong>hidrür</strong> fazının ağırlığındaki depolanan hidrojen<br />

ağırlığıdır. Hidrojen kapasitesini ifade etmek için çok fazla kullanılmamış fakat<br />

mobil ve portatif uygulamalar için önemli teknolojik etkiye sahip diğer bir yol birim<br />

hacimdeki hidrojen atomu sayısıyla, volumetrik kapasitenin belirtilmesidir.<br />

Genellikle tersinir kapasite (∆NH) ve kristal hacim bu sayıyı hesaplamak için<br />

kullanılır. Bu tanımlama kristal içindeki boş hacimleri içermez (her zaman aynı<br />

ölçüde bulunan) ve bu sonuçta en üst hacmi gösterir.<br />

1.3.1.5. Hidrürleşme kinetikleri<br />

Hidrojenasyon işleminde gaz fazındaki bir molekül atomlarına ayrılır ve bu atomlar<br />

daha sonra <strong>metal</strong> <strong>hidrür</strong> kristal kafesine bağlanır. Bütün kinetikler hız-limiti adımı<br />

olarak adlandırılan, yavaş bir adımla sınırlandırılmıştır. Laboratuar işlerinde ve yeni<br />

alaşımların geliştirilmesinde ısı ve kütle transferi tavsiye edilir ancak hız limiti adımı<br />

25


istenmez (Gerard and Ono, 1992). Bunun anlamı maksimum ısı iletimi ve<br />

minimumum hidrojen akış direncini sağlamak için örnek dizaynına çok büyük önem<br />

verilmesi gerekir. Diğer yandan, pratik uygulamalarda, ısı ve kütle transferi son<br />

derece önemlidir. Örnek olarak magnezyum hidrojenasyonunu ele alınırsa; şu andaki<br />

en hızlı reaksiyon kinetikleri 150 ve 250 °C’de %1 mol Nb2O5 ilave edilerek çelik<br />

toplarla öğütülmüş 30 s’de %5’ten daha fazla absorbe edebilmiş MgH2<br />

alaşımlarından elde edilmiştir (Hanada et al., 2006). Bunun anlamı sıcaklığı sabit bir<br />

değerde tutmak birim mol başına 2 kW’lık bir gücün alaşımdan uzaklaştırılması<br />

gerekir. Bu büyük miktardaki ısının depolama tankından dışarıya transfer edilmesi<br />

gerekir. 500 kW ısı gücüne sahip 5 kg’lık bir hidrojen depolama tankı için her<br />

defasında absorpsiyon süresi 300 s olup, bu süre endüstri tarafından da kabul edilmiş<br />

depolama süresidir. Bu basit hesap hidrojenasyon kinetikleri veya en azından büyük<br />

oluşum ısılı <strong>hidrür</strong>lerin ısı transferiyle kontrol edilebileceğini göstermektedir.<br />

Metal <strong>hidrür</strong> tanımlaması için; verilen bir reaksiyon için kinetik eğri zamana göre<br />

değişen bir denkleme dönüştürülmüştür. Reaksiyon hız limit adımlarını anlamak için<br />

bu eğrinin basınç ve sıcaklığa bağımlılığının araştırılması gerekir. Hidrojen <strong>metal</strong><br />

reaksiyonunda reaksiyon ısısından dolayı, oldukça hızlı reaksiyon, <strong>hidrür</strong> fazında<br />

düşük ısıl iletkenlik ve ürün gevrekliği gibi nedenlerden dolayı durum özellikle<br />

karmaşıktır (Gerard and Ono, 1992). Numune boyutu, partikül boyutu, yüzey<br />

özellikleri, katı ve gaz fazının saflığı verilen numunenin kinetikleri üzerine önemli<br />

etkiye sahip faktörlerdir. Dahası iç hız adımlarını belirleme, sadece tozlar üzerinde<br />

değil, ayrıca iri taneli numuneler üzerinde de gerçekleştirilmelidir (Mintz and Bloch,<br />

1985). Tüm bu işlemler farklı deneysel çalışmalar karşılaştırılarak büyük bir dikkatle<br />

yapılmalıdır.<br />

1.3.1.6. Metal <strong>hidrür</strong> ömrü<br />

Pratik uygulamalar için hidrojen depolama sisteminin dayanım özellikleri (kapasite,<br />

kinetikler, tersinirlik, plato basıncı vb.) nispeten cihazın tüm ömrü boyunca sabit<br />

kalması esastır. Farklı uygulamalarda farklı miktarlarda çevrime ihtiyaç vardır ancak<br />

genellikle uygulamaların çoğunda çevrim sayısının birkaç yüzden birkaç bin arasında<br />

26


olması gerekir. Çevrim ömrü genel olarak AB5 ve magnezyum esaslı bileşikleri için<br />

incelendiğinde çevrim ömrünün birçok parametreye bağlı olduğu görülmektedir.<br />

Bunlardan en önemlilerinden birisi inter<strong>metal</strong>ik bileşiklerin <strong>hidrür</strong>leri kimyasal<br />

tepkimelerin etkisiyle yarı kararlı durumdadır ve bu da çevrim ömrünün kısalmasına<br />

neden olmaktadır. Çevrim ömrünü kısaltan diğer sebepler ise: hidrojen gazındaki<br />

kirliliğin zehirlenmeye neden olması, tanecik yığılması, yapısal gevşemeler ve kristal<br />

büyümesi, faz değişimi ve <strong>hidrür</strong> olmayan faz oluşumudur. Yinede AB5 ve<br />

magnezyum esaslı bileşiklerde hidrojen absorpsiyon özelliklerinde köklü<br />

değişiklikler olmadan çevrim bir kaç binin üzerinde sürdürebilmektedir. Çevrim<br />

kararlığını devam ettirebilmek için kullanılan en önemli yol element değişimi<br />

yöntemidir(Hirscher, 2010).<br />

1.3.1.7. Kristal çatlaması<br />

Kristal çatlaması, hidrojenasyon ve <strong>hidrür</strong> alaşımlarının kırılgan yapısından dolayı<br />

hacim değişimi nedeniyle küçük partiküllü tozlar içerisindeki alaşım parçacıklarının<br />

kendi kendine toz haline gelmesidir. Bu da depolama tankı içerisindeki ısı transferi<br />

ve gaz akışını değiştirecek bir malzeme yoğunluğuna neden olabilmektedir<br />

(Hirscher, 2010). Bu etki <strong>hidrür</strong>leşme ve de<strong>hidrür</strong>leşme durumları için genellikle<br />

farklı iç ısı iletkenlikleriyle geliştirilmiştir. Kristal çatlaması olayı <strong>hidrür</strong> reaktör için<br />

ciddi bir problem oluşturabilmektedir, çünkü küçük partiküller tankın alt bölümlerine<br />

düşerek, hidrojenasyonla birlikte kafes genişlemesinden dolayı tankta çatlamalara yol<br />

açacak boyutta genişleyebilmektedir.<br />

1.3.2. Metal <strong>hidrür</strong>ler arası bağlar<br />

Metal <strong>hidrür</strong>ler bir ana <strong>metal</strong> ve hidrojen arasında yoğunlaşmış tek fazlı bileşik<br />

olarak tanımlanabilir (Hirscher, 2010). Đlk <strong>metal</strong> <strong>hidrür</strong>, Graham tarafından bulunmuş<br />

olan ve büyük miktarlarda hidrojen absorbe edebilen paladyumdur. Basit ikili <strong>metal</strong><br />

<strong>hidrür</strong>ler <strong>metal</strong>-hidrojen bağının niteliğine göre iyonik <strong>hidrür</strong>ler, kovalent <strong>hidrür</strong>ler<br />

ve <strong>metal</strong>ik <strong>hidrür</strong>ler olmak üzere üç temel gruba ayrılabilir (Libowitz, 1991). Şekil<br />

1.14’deki periyodik tabloda <strong>hidrür</strong>ler arası bağlara göre elementler verilmiştir.<br />

27


1.3.2.1. Đyonik <strong>hidrür</strong>ler<br />

Şekil 1.14. Đkili <strong>hidrür</strong>ler tablosu<br />

Elektron verme eğiliminde olan bir atomla elektron alma eğilimindeki bir atom<br />

reaksiyon alanına girdiklerinde, aralarında bir elektron alışverişi olur. Bunun<br />

sonucunda elektron alan atom negatif iyon (anyon), elektron veren atom ise pozitif<br />

iyon (katyon) haline gelecektir. Bu şekilde aralarında elektrostatik çekme kuvveti<br />

yaratılan atomlar iyonik bağla bir bileşik oluştururlar. Đyon bağına elektro valans bağ<br />

da denilmektedir. Đyonik bağ oluşumunda, <strong>metal</strong>, düşük elektronegatifliği nedeniyle<br />

bir elektron vererek pozitif bir iyon (katyon) oluşturur. Normal sofra tuzunda,<br />

sodyum ile klor iyonları birbirlerine iyonik bağ ile bağlıdır. Đyonik bağ genellikle<br />

<strong>metal</strong>lerle a<strong>metal</strong>ler arasında gerçekleşir. A<strong>metal</strong> atomlarının elektronegatifliği<br />

yüksektir ve kolayca elektron alıp negatif iyon (anyon) oluşturabilirler. Dolayısıyla,<br />

iki veya daha fazla iyon, elektrostatik kuvvetlerin etkisiyle birbirlerini çekerler. Bu<br />

tür bağlar, hidrojen bağından daha kuvvetli fakat kovalent bağ ile hemen hemen aynı<br />

kuvvettedir.<br />

28


Bu grup tüm alkali <strong>metal</strong>lerin ikili <strong>hidrür</strong>lerini ve kalsiyumdan Baryum’uda içine<br />

alan toprak alkali <strong>metal</strong>leri içermektedir (Libowitz, 1965). Bu bileşiklerde, hidrojen<br />

negatif yüklü (H - ) iyonlardan oluşmaktadır ve halojen serisinin bir üyesi olarak kabul<br />

edilebilir. Bu nedenle, sertlik, kırılganlık, optik özellikleri ve kristal yapıları gibi<br />

birçok fiziksel özellikleri ilgili halojenürlere benzerler. Toprak alkali <strong>hidrür</strong>ler<br />

baryum klorür yapıya sahipken, alkali <strong>metal</strong> <strong>hidrür</strong>ler sodyum klorür yapıya sahiptir<br />

(Ivey and Northwood, 1983). Sodyum <strong>hidrür</strong> (NaH) ve kalsiyum <strong>hidrür</strong> (CaH2) tipik<br />

ikili iyonik <strong>hidrür</strong>lerdir. Đyonik <strong>hidrür</strong>ler ergime sıcaklıklarında veya hemen altında<br />

yüksek iletkenliğe sahiptir. Lityum alüminyum <strong>hidrür</strong> LiAlH4 ve sodyum bor <strong>hidrür</strong><br />

NaBH4 gibi kompleks iyonik <strong>hidrür</strong>ler ticari olarak indirgeyici madde olarak<br />

kullanılmaktadır. Genel olarak ikili iyonik <strong>hidrür</strong>ler magnezyum dışında hidrojen<br />

depolama uygulamalarında oldukça kararlıdır (Ivey and Northwood, 1983). Buna<br />

karşın magnezyum <strong>hidrür</strong>ler gerçek iyonik <strong>hidrür</strong> değildir. Magnezyum <strong>hidrür</strong>ler<br />

hidrojen ve magnezyumun etkileşimiyle oluşan reaksiyon kısmen iyonik kısmen de<br />

kovalenttir. Bu nedenle magnezyum <strong>hidrür</strong>ler iyonik ve kovalent <strong>hidrür</strong>ler arasında,<br />

bir geçiş <strong>hidrür</strong>ü olarak kabul edilmektedir.<br />

1.3.2.2. Kovalent <strong>hidrür</strong>ler<br />

Kovalent bağ, iki atom arasında, bir veya daha fazla elektronun paylaşılmasıyla<br />

karakterize edilen kimyasal bağdır. Genellikle bağ, ortaya çıkan molekülü bir arada<br />

tutan ortak çekim gücü olarak tanımlanabilir. Paylaşılan elektron ya da elektronlar,<br />

her iki çekirdek etrafında dolanacaklar, iki çekirdek arasındaki bölgede daha uzun<br />

süre bulundukları için bu bölgede (-) yüklü bir alan yaratacaklardır. Bu alan, her iki<br />

çekirdeğe bir çekme kuvveti uygulayarak bir bağ yaratır.<br />

Kovalent bağ, genellikle benzer elektronegatifliğe sahip atomlar arasında gerçekleşir.<br />

Bu nedenle a<strong>metal</strong>ler, daha kolaylıkla kovalent bağı tercih eder ve <strong>metal</strong>ler de<br />

kolayca yerlerinden oynatılabilen elektronların daha serbestçe dolaşabildiği <strong>metal</strong>ik<br />

bağ yaparlar. Burada benzer elektronegatif atomlar elektronların çiftlerini paylaşırlar.<br />

Genellikle kovalent <strong>hidrür</strong>ler düşük erime ve kaynama noktalarına sahiptir.<br />

Moleküller arasında zayıf van der Waals kuvvetlerinden dolayı çoğu kovalent<br />

29


<strong>hidrür</strong>ler oda sıcaklıklarında sıvı veya gaz halindedir ve katı halde ısıl dengede<br />

değildir (Maeland, 1977). Kovalent <strong>hidrür</strong>lere su (H2O), hidrojen sülfür (H2S), Silan<br />

(SiH4), alüminyum bor <strong>hidrür</strong> (AlBH4), metan (CH4) ve diğer hidrokarbonlar örnek<br />

olarak verilebilir. Kovalent <strong>hidrür</strong>ler hidrojen gazı ve elementinin diğer reaksiyonu<br />

ile oluşturulamazlar. Sentezlenmeleri için kompleks kimyasal reaksiyonlar<br />

uygulanmalıdır. Sentezlemedeki bu zorluktan dolayı, kovalent <strong>hidrür</strong>ler hidrojen<br />

depolama uygulamaları için uygun değildirler.<br />

1.3.2.3. Metalik <strong>hidrür</strong>ler<br />

Metalik bağ, esas olarak <strong>metal</strong>ler arasındaki, bir ya da daha fazla atomu bir arada<br />

tutan bir kimyasal bağ türüdür. Metal atomlarının lattisindeki serbest elektronların<br />

yer değiştirmiş olarak paylaşılması esasına dayanır. Metal atomları valans bandında<br />

dış yörüngelerinde periyotlarına veya enerji seviyelerine oranla yüksek sayıda<br />

elektron içerirler. Bunlar atomdan ayrılarak bir pozitif iyonun, lattis etrafını saran bir<br />

elektron zarfı, gazı oluştururlar. Negatif elektron gazı ile pozitif çekirdekler arasında<br />

oluşan elektrostatik çekme kuvvetleri <strong>metal</strong> atomlarının bir arada kalmalarını sağlar.<br />

Burada, kovalent ve iyonik bağlardaki gibi merkezi bir bağ söz konusu değildir.<br />

Metalik bağ, düzgün pozitif iyon yığını ile bu yığını kuşatan elektron gazı arasında<br />

ortaya çıkar. Metalik bağ, <strong>metal</strong> atomları veya iyonlar ile yer değiştirmiş elektronlar<br />

arasındaki elektrostatik çekimdir. Metal atomlarının en az bir valans elektronu vardır<br />

ve bu elektronu ne bir komşu atomla paylaşır ne de iyon oluşturmak üzere verirler.<br />

Bunun yerine, <strong>metal</strong> atomlarının dış enerji seviyeleri üst üste biner ve bu anlamda<br />

kovalent bağa benzerlik gösterirler.<br />

Doğada <strong>metal</strong>ik olarak bulunan <strong>hidrür</strong>lerin çoğu hidrojen depolama için<br />

kullanılabilir. Hidrürlerin çoğu nadir <strong>metal</strong>ler ve aktinit serisi de dahil olmak üzere<br />

geçiş <strong>metal</strong>lerinden oluşmaktadır. Bu <strong>hidrür</strong>ler <strong>metal</strong>ler gibi hidrojen depolar ve<br />

<strong>metal</strong>ik bağ formundadır. Bunlar yüksek ısıl ve elektriksel iletkenliğe sahiptir.<br />

Ancak <strong>metal</strong>lerin aksine oldukça kırılgandırlar (Maeland, 1977).<br />

30


1.3.3. Metal <strong>hidrür</strong> çeşitleri<br />

Hidrojen son derece reaktif bir elementtir. Binlerce <strong>metal</strong> ve alaşımlarla katı çözelti<br />

şeklinde <strong>hidrür</strong> biçiminde olduğu bilinir. Birçok doğal element uygun şartlar altında<br />

hidrojen absorbe etmektedir. Fakat Şekil 1.15’de gösterildiği gibi bunların<br />

birçoğunun pratik uygulamalarda kullanılabilmeleri için PCT özellikleri uygun<br />

değildir. Şekildeki grafikte kutucuk içerisinde gösterilen sıcaklık aralığı 0–100 °C ve<br />

basınç aralığını 1–10 atm olarak göstermiştir. Đstenilen PCT özelliklerini elde<br />

edebilmek amacıyla bazı güçlü <strong>hidrür</strong> elementleri (A olarak adlandırılmıştır) ile daha<br />

güçsüz <strong>hidrür</strong>leme özelliklerine sahip olan (B olarak adlandırılmıştır) elementleri<br />

birleştirmek gerekmektedir. Đnter<strong>metal</strong>ik bileşiklerinin çeşitli tiplerinde bileşikler<br />

farklı miktarlarda A ve B elementinden içerirler, bileşikler daha çok katı solüsyon<br />

şeklindedir ve <strong>hidrür</strong> bileşikleri geçiş <strong>metal</strong>leri formundadır.<br />

Çizelge 1.2. Bazı önemli inter<strong>metal</strong>ik bileşiklerin protatipi ve yapısı (Züttel, 2001)<br />

Đnter<strong>metal</strong>ik<br />

Bileşikler<br />

Prototip Hidrür Yapı<br />

AB5 LaNi5 LaNiH6 Hook faz, altıgen<br />

AB2 ZrV2, ZrMn2, TiMn2 ZrV2H5.5 Lave fazı, altıgen veya kübik<br />

AB3 CeNi3, YFe3 CeNi3H4 Altıgen, PuNi3-tipi<br />

A2B7 Y2Ni7, Th2Fe7 Y2Ni7H3 Altıgen, Ce2Ni7-tipi<br />

A6B23 Y6Fe23 Ho6Fe23H12 Kübik, Th6Mn23-tipi<br />

AB TiFe, ZrNi TiFeH2 Kübik, CsCl- veya CrB-tipi<br />

A2B Mg2Ni, Ti2Ni Mg2NiH4 Kübik, MoSi2- veya Ti2Ni-tipi<br />

A elementi genellikle nadir elementler veya kararlı <strong>hidrür</strong> oluşturabilen toprak alkali<br />

<strong>metal</strong>lerdir. B elementi çoğunlukla kararsız <strong>hidrür</strong> oluşturabilen bir geçiş <strong>metal</strong>idir.<br />

Đnter<strong>metal</strong>ik bileşiklerde B ile A arasındaki bazı belirlenmiş oranlar x=0.5, 1, 2, 5<br />

şeklindedir. Farklı tiplerdeki <strong>hidrür</strong> aileleri Çizelge 1.2’de verilmiştir.<br />

31


Şekil 1.15. Bazı doğal elementler için van’t Hoff diyagramı (Sandrock, 1999)<br />

1.3.3.1. AB5 inter<strong>metal</strong>ik bileşikleri<br />

AB5 bileşikleri 100 °C sıcaklıkta ve birkaç atmosfer denge basınçlarında genellikle<br />

<strong>hidrür</strong> formundadır. Diğer önemli özellikleri düşük histeresiz, gaz kirliliğindeki<br />

toleranslar ve başlangıç çevriminde aktivasyon kolaylığıdır. Böylelikle, bu<br />

alaşımların çoğu ana <strong>metal</strong> olarak Ni içerir. AB5 ailesi CaCu5 tipi kafesle altıgen<br />

yapıya sahiptir. AB5 sınıfı alaşımların en önemli örneklerinden birisi LaNi5<br />

alaşımıdır. Bu sınıf alaşımlar olağanüstü çok yönlülüğe sahiptir. Çünkü birçok farklı<br />

element türü A ve B kafes boşlukları yerine geçebilir. A <strong>metal</strong>i için genel örnekler<br />

Mm, Ca, Y, Zr ve B için Al, Mn, Si, Zn, Cr, Fe, Cu ve Co verilebilir. A ve B’nin<br />

yerine geçen bileşikler alaşım makro yapısında önemli değişiklere neden olur.<br />

Örneğin, LaNi5 alaşımında Ni Co ve Fe ile yer değiştirirse <strong>hidrür</strong>leşme üzerinde<br />

32


hacim genişlemesi oluşur, korozyon hızı düşer ve <strong>hidrür</strong> çevrim ömründe gelişmeler<br />

sağlar (Hottinen, 2001).<br />

AB5 alaşımlarının PCT özellikleri Şekil 1.16’da verilmiştir. AB5 alaşımlarında<br />

MmNi5 hariç, histeresiz oldukça düşüktür (Sandrock, 1999). AB5 alaşımlarının en<br />

önemli avantajı koruyucu oksit formunda olmayışları ve H2 içerisindeki küçük O2 ve<br />

H2O kirliliğine karşı dayanıklı olmalarıdır. AB5 bileşiklerinin hidrojen kapasiteleri<br />

düşüktür. Yaklaşık %1 nominal miktarla LaCo5 ve Ca0.7Mm0.3Ni5 ile maksimum<br />

kapasite ağırlıkça %1.90 olarak gerçekleşmiştir. Tersinir kapasite ise LaNi5 ile<br />

ağırlıkça %1.28 olarak belirlenmiştir, nominal miktar ise ağırlıkça yaklaşık %0.7–1<br />

olmuştur (Sandrock, 1999).<br />

Şekil 1.16. Çeşitli AB5 alaşımları için Van’t Hoff diyagramı (Sandrock, 1999)<br />

Diğer yandan AB5 alaşımlarının maliyeti de oldukça yüksektir. Reaktör dahil<br />

edilmeden sadece hammadde göz önüne alındığında alaşım maliyetleri (1999) 8–11<br />

33


€/kg arasındadır. Bu hidrojen depolamanın tersinir kapasitesi de eklenerek hidrojen<br />

depolama maliyetine çevrildiğinde maliyet yaklaşık 0.7–2.0 €/gH olmaktadır<br />

(Sandrock, 1999).<br />

1.3.3.2. AB2 inter<strong>metal</strong>ik bileşikleri<br />

AB5 alaşımları gibi AB2 alaşımları da çevre sıcaklıklarında geniş ve çok yönlü<br />

olumlu PCT özelliklerine sahiptir. A elementi genellikle Ti, Zr, Hf, TH veya bir<br />

lantanittir (atom numarası 57–71). B elementi V, Cr, Mn ve Fe gibi geçişli ya da<br />

geçişsiz <strong>metal</strong>lerden biri olabilir. Bazı AB2 alaşımlarının PCT özellikleri Şekil<br />

1.17’de verilmiştir. AB2 alaşımlarının hidrojen kapasiteleri tersinir esaslı olarak AB5<br />

ile karşılaştırılabilir fakat genellikle yüksek kapasiteye sahiptir. Bu alaşımlarda<br />

özellikle geniş basınç ve sıcaklıklar olanaklıdır.<br />

Şekil 1.17. Bazı AB2 alaşımları için van’t Hoff diyagramı (Sandrock, 1999)<br />

34


AB2 alaşımlarının maksimum kapasiteleri genellikle ağırlıkça %1.5–2 arasında ve<br />

TiCr1.2V0.8 çifti ile ağırlıkça %3.4 kapasiteye ulaşılabilmektedir. Tersinir kapasiteler<br />

ise ağırlıkça %0.9–1.3 arasındadır (Sandrock, 1999).<br />

Hammaddelerin maliyetleri 2.5–12 €/kg arasındadır. AB2 alaşımlarının hidrojen<br />

esaslı maliyetleri yaklaşık 0.4–1.3 €/gH’dir. Böylece AB2 alaşımlarının maliyetleri<br />

AB5 alaşımlarının maliyetlerine göre bir miktar daha düşüktür. Örneğin<br />

Ti0.92Zr0.02V0.43Fe0.09Cr0.05Mn1.5 için maksimum ve tersinir kapasitesi sırasıyla<br />

ağırlıkça %1.9-%1.3 ve maliyeti 0.41 €/gH’dır (Hottinen, 2001).<br />

1.3.3.3. AB inter<strong>metal</strong>ik bileşikleri<br />

AB alaşımları genel olarak TiFe esaslıdır. Bunlar iki platoya sahip olma<br />

eğilimindedir, bunlardan üstte olan çok kararlı değildir ve histeresiz oldukça geniştir.<br />

Şekil 1.18. Bazı AB alaşımları için van’t Hoff diyagramı (Sandrock, 1999)<br />

35


PCT özellikleri Şekil 1.18’de gösterildiği gibi kısmi olarak Fe yerine Mn veya Ni<br />

kullanılarak düzeltilebilir. Şekilde U üst platoyu ve L alt platoyu göstermektedir.<br />

Oksit filmleri genellikle TiFe esaslı alaşımları biçimindedir. Bu hidrojen içindeki<br />

kirlilik duyarlılığını azaltır. Diğer yandan oksitlenme yanma davranışları için bir<br />

eğilim oluşturmaz.<br />

AB alaşımları için maksimum hidrojen kapasitesi ağırlıkça %1-2’den daha düşük<br />

miktarda değişmektedir. Burada ilgi çeken iki alaşım, TiFe ve TiFe0.85Mn0.15,<br />

maksimum kapasitesi ağırlıkça yaklaşık %1.9 ve tersinir kapasitesi %1.5’dir. Bu<br />

alaşımların maliyetleri sırasıyla 0.34 ve 0.35 €/gH olarak oldukça düşüktür. AB<br />

alaşımlarının ticari olarak kullanılamamasının ana sebebi üst plato ve kirlilik<br />

duyarlılığının kararsız olmasıdır. Bazı çalışmalarda (ZrCo)1-x(TiNi)x gibi sözde ikili<br />

<strong>hidrür</strong>ler kullanılarak bu problemlerin aşılması amaçlanmıştır, fakat çok açık verim<br />

artışları elde edilememiştir. Örneğin (ZrCo)0.7(TiNi)0.3 alaşımında H/M=1.6 ve<br />

hidrojen kapasitesi ağırlıkça %1.16’dır (Shimizu et al., 1999; Sandrock, 1999).<br />

1.3.3.4. Diğer inter<strong>metal</strong>ik bileşikler<br />

Yukarıda gösterilen alaşımlara ilaveten, burada birçok inter<strong>metal</strong>ik aileleri hidrojen<br />

absorpsiyon kapasitesine sahiptir, fakat bunlar ticari olarak ilgi çekmemektedir.<br />

Bunlar A2B, AB3, A2B7, A3B alaşımlarını içine alır. Bu alaşımlardan bazıları iyi<br />

hidrojen kapasitelerine sahiptir fakat PCT özellikleri uygun değildir. PCT<br />

özelliklerinin yanında pratikte uygulama alanlarında fakir hidrojen kapasitesine<br />

sahiptir veya çok dar plato eğilimine sahiptir.<br />

1.3.3.5. Katı çözelti alaşımları<br />

Bir katı çözelti alaşımında yöntem birinci elementin bir minör element içerisinde<br />

çözünmesidir. Katı çözelti alaşımları Pd esaslı alaşımlarda iyi PCT özelliklerine<br />

sahiptir. Ancak bu alaşımlar genellikle ağırlıkça yaklaşık %1 olarak çok düşük<br />

hidrojen depolama kapasitesine sahiptir. Bunlar içerisinde vanadyum esaslı alaşımlar<br />

daha yüksek depolama kapasitesine sahiptirler, fakat bunlar oldukça pahalıdır.<br />

36


Özellikle V-Ti-Fe kombinasyonu hidrojen depolama için önemli bir alternatiftir.<br />

Örneğin %3.7 hidrojen depolama kapasitesiyle (V0.9Ti0.1)0.95Fe0.05 maksimum<br />

hidrojen kapasitesine sahiptir ve tersinir hidrojen kapasitesi ağırlıkça %1.8 ve hala en<br />

iyi PCT özelliklerine sahip alaşımdır. Vanadyum esaslı (V0.9Ti0.1)0.95Fe0.05<br />

alaşımlarının maliyeti yaklaşık 11.76 €/kg veya 0.65 €/gH’dir (Sandrock, 1999).<br />

37


2. KAYNAK ÖZETLERĐ<br />

Metal <strong>hidrür</strong>lere olan ilgi 1970 yılında, LaNi5 ve FeTi alaşımlarının hidrojen<br />

depoladığının keşfi ile başlamıştır (Reilly and Wiswall 1974). Bu dönemlerde ülkeler<br />

dünyanın artan enerji ihtiyacı ve petrol rezervlerindeki azalma yüzünden yeni<br />

kaynakların arayışı içerisindeydi. Bu arayış sırasında hidrojen enerjisinin alternatif<br />

bir yakıt olarak kullanılacağının anlaşılması ile birlikte bu konu üzerine yapılan<br />

çalışmalar hız kazandı. Fakat yapılan çalışmalar hidrojen enerjisinin etkin ve verimli<br />

bir şekilde kullanılabilmesinin önündeki en büyük engelin hidrojenin depolanma<br />

problemi olduğunu gösterdi. Hidrojenin bilinen en hafif gaz olması, depolanmasını<br />

büyük bir probleme dönüştürmekteydi. Nitekim ilk petrol krizinin ardından 1975<br />

yılında Amerika, Japonya ve Almanya’da pek çok araştırma programı başlatılmıştır<br />

(Mendelsohn et al., 1979).<br />

Ayrıca sera gazlarının olumsuz etkileri ve fosil yakıtların tükenmesi ile ilgili<br />

endişeler alternatif bir enerji kaynağı kullanımının ne kadar önemli olduğunu ortaya<br />

koymaktadır. Yüksek kalorifik değeri ve temiz bir yakıt olması sebebiyle hidrojen,<br />

üzerinde en çok çalışılan enerji taşıyıcı olarak dikkat çekmektedir. Fakat hidrojenin<br />

depolama problemleri gerek uygulama alanlarını gerekse de ticarileşmesini<br />

kısıtlamaktadır. Son zamanlarda küçük hacimlerde yüksek miktarlarda hidrojen<br />

depolayabilen ve yüksek basınçlarda çalışma olanağı sağlayan <strong>metal</strong> <strong>hidrür</strong> reaktörler<br />

büyük ilgi görmektedir. Bu yüzden bilim adamları çalışmalarını nadir <strong>metal</strong> alaşımlar<br />

üzerinde yoğunlaştırmıştır. Literatürde <strong>metal</strong> <strong>hidrür</strong> tanklarda hidrojen depolama ile<br />

ilgili birçok deneysel ve nümerik çalışma mevcuttur. Hidrojenin <strong>metal</strong> <strong>hidrür</strong><br />

reaktörlerde kimyasal olarak depolanması göz önüne alındığında Mg, La, Na ve Li<br />

gibi <strong>metal</strong>ler dikkat çekmektedir.<br />

Lucas ve Richard (1984), ısı taşınım denklemini ve hidrojen depolama sisteminin<br />

performansını, tek boyutlu olarak ele almışlardır. Geliştirdikleri matematiksel modeli<br />

nümerik olarak çözmüşlerdir. Fakat geliştirilen modelde reaksiyonla açığa çıkan ısı<br />

hesaba katılmamıştır.<br />

38


Mayer vd. (1987), hidrojenin depolanması esnasında oluşan ısı ve kütle transferini<br />

deneysel ve teorik olarak incelenmişlerdir. Çalışmalarında <strong>hidrür</strong><br />

konsantrasyonunun, basınç ve sıcaklığın, zamana ve yatak çapına bağlı olarak<br />

değişimlerini elde etmişlerdir.<br />

Bernauer (1988), çalışmasında, <strong>metal</strong> <strong>hidrür</strong> teknolojisi alanında son on yıldır<br />

yapılan çalışmalar neticesinde meydana gelen gelişmeleri özetlemiş, hidrojenin<br />

depolanmasıyla <strong>hidrür</strong>'ün hareketli ve hareketsiz uygulamalarda, emniyetli bir<br />

şekilde kolaylıkla kullanılabileceğini göstermiştir.<br />

Sun ve Deng (1990), <strong>metal</strong> <strong>hidrür</strong> yataklarında iki boyutlu zamana bağlı ısı ve kütle<br />

transferi için nümerik bir model sunmuşlardır. Silindirik ve Kartezyen koordinatlarda<br />

ısı ve kütle transferini bir boyutlu ve iki boyutlu olarak incelemişlerdir. Bu çalışmada<br />

diferansiyel denklemler sonlu farklar kullanılarak impilisit olarak çözülmüştür. Bu<br />

çözümden elde ettikleri nümerik sonuçlarla deneysel sonuçların uyum içinde<br />

olduğunu göstermişlerdir.<br />

Gopal ve Murthy (1992), silindirik biçimli bir <strong>metal</strong> <strong>hidrür</strong> yatak içerisinde ısı ve<br />

kütle transferi için bir boyutlu bir matematiksel model geliştirmişler. Yapılan çalışma<br />

sonucunda alaşımlar ve çalışma şartları için geniş bir aralık içeren boyutsuz formda<br />

veriler elde edilmiştir. Bu veri yatak özellikleri ve çalışma parametrelerinin bir terimi<br />

olarak boşluk, ortalama sıcaklık ve konsantrasyon için korelasyon elde etmede<br />

kullanılmıştır. Bu korelasyon LaNi4.7Al0.3 yatakta örnek olarak gösterilmiştir.<br />

Jemni ve Nasrallah (1995), <strong>metal</strong>-hidrojen reaktörlerinde absorbsiyon işlemi için iki<br />

boyutlu bir matematiksel model geliştirmişlerdir. Yatak geometrisi, ısı transfer<br />

katsayısı ve şarj basıncı gibi proses parametrelerinin <strong>hidrür</strong> oluşumu üzerindeki<br />

etkilerini araştırmışlardır. Diğer bir çalışmalarında ise hidrojen akışını hesaba katarak<br />

katı ve gazın termal dengesini incelemişler ve reaktörde oluşan ısı ve kütle transferini<br />

karakterize eden bir matematiksel model geliştirmişlerdir.<br />

39


Nasrallah ve Jemni (1997), <strong>metal</strong>-hidrojen reaktöründe matematiksel olarak sıcaklık<br />

ve kütle analizi yapmışlar ve çalışmalarında şu hipotezleri kullanmışlarıdır; Katı ve<br />

gaz, aynı sıcaklıktadır. Reaktörde basınç değişmesinin etkisi ihmal edilebilirdir.<br />

Denge basıncı değişimine hidrojen konsantrasyonunun etkisi ihmal edilebilirdir.<br />

Burada yazarlar bu hipotezler ve hipotezler olmadan elde edilen nümerik sonuçlara<br />

göre hipotezlerin geçerliliğini test etmişlerdir. LaNi5 hidrojen reaktörü için yapılan<br />

analizler sonucunda ortaya atılan hipotezlerin geçerli olduğunu ortaya koymuştur.<br />

Jemni vd. (1999), yaptıkları çalışmada <strong>metal</strong>-hidrojen (LaNi5-H2) reaktörü deneysel<br />

ve teorik olarak incelenmiştir. Çalışmada öncelikle deneysel olarak ısıl iletkenlik,<br />

<strong>hidrür</strong> yatak ve reaktör çevresindeki akışkan arasındaki iletim, denge basıncı ve<br />

reaksiyon kinetikleri, başlangıç koşullarını hesaba katarak, sıcaklık ve zamana bağlı<br />

hidrojen basınçlarını belirlemeyi hedeflemişlerdir. Đkinci hedefi ise teorik ve<br />

deneysel sonuçları karşılaştırarak teorik modelin geçerliliğini belirlemektir. Deneysel<br />

çalışmada denge basıncı, reaksiyon kinetikleri ve LaNi5-H2 efektif ısıl iletkenliği<br />

belirlenmiştir. Reaktör içerisinde ölçülen sıcaklık değerleriyle, efektif ısıl iletkenlik<br />

(λc) ve <strong>hidrür</strong> yatağı reaktör çevresindeki akışkan arasındaki ısı iletimi<br />

hesaplanmıştır. Isıl iletkenlik 1.087 W/mK olarak belirlenmiştir.<br />

Johari vd. (2001), yaptıkları çalışmada MmNi5 ve MmNi4.85Fe0.15 için hidrojen gazı<br />

depolama karakteristikleri gösterilmiştir. Metallerin sıcaklık ve basınç etkisi altındaki<br />

absorpsiyon kapasitesi gösterilmiştir. Çalışma şartlarının kısıtlamalarından basınç 2–<br />

30 bar arasında, sıcaklık 298–323 K aralığında seçilmiştir. Çalışmada MmNi5<br />

alaşımının hidrojen depolama kapasitesinin çalışma sıcaklığı ile ters orantılı olarak<br />

değiştiği belirlenmiştir. Bu çalışmada MmNi5 için maksimum absorpsiyon 298 K<br />

sıcaklığında başarılmıştır ve plato basıncı 10 bar olarak ölçülmüştür. Hidrojen içeriği<br />

kütlesel ve hacimsel olarak %1.20 olarak hesaplanmıştır. MmNi5 için 313 ve 323 K<br />

sıcaklığında hacimsel ve kütlesel sonuçlar %0.9 ve %0.8 bulunmuş, plato basıncı<br />

sırasıyla 20 ve 24 bar olarak hesaplanmıştır.<br />

Mat ve Kaplan (2001), yaptıkları çalışmada LmNi5 <strong>hidrür</strong> reaktörlerde hidrojen<br />

absorpsiyonunu matematiksel olarak incelemişlerdir. LmNi5 depolama tanklarında<br />

40


<strong>metal</strong> <strong>hidrür</strong> formunda sürekli matematiksel model ile nümerik olarak çözülmüştür.<br />

Model kompleks ısı/kütle transferi ve <strong>hidrür</strong> yatakdaki kimyasal reaksiyonları<br />

içermektedir. Düşük denge basıncının <strong>hidrür</strong> formunu geliştirdiği belirlenmiştir.<br />

Hidrür işleminin başlangıç zamanlarında absorbe edilen hidrojen kütlesi hızlı şekilde<br />

artarken daha sonra reaksiyon ısısı nedeniyle yatak sıcaklığının artmasından sonra<br />

depolama hızı yavaşlamıştır. Literatür değerleriyle karşılaştırıldığında bulunan<br />

sonuçların uyum içerisinde olduğu görülmüştür. Sonuçlar seçilen bölgelerdeki<br />

izoterm, izokonsantrasyon ve konsantrasyon histeresizlerini göstermiştir. Hidrojen<br />

formasyonunda sıcaklığın güçlü bir fonksiyonu olduğu için ana parametre denge<br />

basıncı olduğu belirlenmiştir. Denge basıncı düşük olduğunda <strong>hidrür</strong>leşme hızlı<br />

gerçekleşmektedir. Başlangıçta <strong>hidrür</strong>leşme çok hızlı gerçekleşmekle beraber daha<br />

sonra ekzotermik reaksiyon sonucu oluşan ısının yükselmesiyle reaksiyon hızı da<br />

düşmektedir. Bu esnada hidrojen girişinde ve soğuyan duvarlarda reaksiyon devam<br />

etmektedir. Nümerik sonuçlar deneysel verilerle karşılaştırıldığında modelin<br />

başarılığı olduğu, deneysel verilerle karşılaştırıldığında sıcaklık profilleri ve <strong>hidrür</strong><br />

formasyonunun uyum içinde olduğu göstermişlerdir.<br />

Güvendiren vd. (2002), yaptıkları çalışmada magnezyuma değişik katkı maddeleri<br />

eklenmiş ve bu maddelerin atritörle yapılan öğütme sonrası oluşan mekanik<br />

alaşımlandırmadaki etkinliği incelenmiştir. Etkinlik yapısal incelme olarak ele<br />

alınmış ve öğütme sonrası incelme, gerek <strong>metal</strong>ografik gerekse X-ışınları kırınım<br />

yöntemleri ile tespit edilmiştir. Bu amaçla sisteme sabit hacim oranında (%5) farklı<br />

sertlik değerlerinde SiC, Al2O3, V, grafit ilaveleri yapılmıştır. Yapılan çalışma<br />

yapısal incelmede katkı maddesinin sertliğinin tek başına belirleyici olmadığını,<br />

katkı maddesinin öğütme sırasında büyüklüğünü muhafazasının daha etkin olduğunu<br />

göstermiştir. Bu açıdan yapılan ilaveler içerisinde en etkininin grafit ve Al2O3 olduğu<br />

gözlenmiştir.<br />

Aldaş vd. (2002), yaptıkları çalışmada <strong>metal</strong>-<strong>hidrür</strong> yatakta hidrojen depolanması<br />

sayısal olarak incelenmiştir. Kullanılan matematiksel model <strong>hidrür</strong> oluşumu sırasında<br />

gerçekleşen kompleks ısı-kütle transferi ve kimyasal reaksiyon içermektedir. Hidrür<br />

oluşumu, <strong>metal</strong> <strong>hidrür</strong> yatağın içindeki sıcaklık ve H/M oranı dağılımları ve sıcaklık-<br />

41


zaman ve H/M oranı-zaman grafikleri yardımıyla incelenmiştir. Nümerik sonuçlar,<br />

<strong>hidrür</strong> oluşumunda en önemli parametrenin denge basıncı olduğunu göstermektedir.<br />

Denge basıncı sıcaklıkla artmaktadır. Kimyasal reaksiyon sonucu oluşan ısı yatakta<br />

sıcaklığın artmasına neden olmaktadır. Hidrür oluşumunun soğutulan yüzeylere<br />

yakın bölgelerde (denge basıncı düşük olduğu için) yüksek, merkeze yakın<br />

bölgelerde de yavaş gerçekleşmektedir. Metal <strong>hidrür</strong> yatağının ısı iletim katsayısı<br />

düşük olduğu için özellikle kimyasal reaksiyonun hızlı olduğu ilk zamanlarda yatakta<br />

merkezle cidar arasında önemli ölçüde sıcaklık gradyanı oluştuğu görülmüştür.<br />

Sayısal olarak hesaplanan H/M oranı, literatürdeki deneysel sonuçlarla<br />

karşılaştırılmış ve bulunan sonuçların uyum içinde olduğu görülmüştür.<br />

Aldaş vd. (2002), yaptıkları çalışmada <strong>hidrür</strong> yataklarda ısı ve kütle transferi, akışkan<br />

akışı ve kimyasal reaksiyonların <strong>hidrür</strong> işlemi üzerine etkisini 3 boyutlu nümerik<br />

olarak incelenmiştir. Hidrürleşmenin, soğutulan duvar sınırlarına yaklaştıkça daha<br />

fazla hızlandığını belirlemişlerdir. Bunların sonucunda sistemde sıcaklık dağılımı,<br />

üzerindeki akışkan akışı nedeniyle hidrojen absorpsiyonunda önemli olmayan<br />

miktarlarda gelişmeler sağlanmıştır. Hidrürleşme önemli derecede sıcaklığın bir<br />

fonksiyonudur. Hidrürleşme başlangıçta yatağın bütününde uniformdur. Reaksiyon<br />

ısısından dolayı sistemde sıcaklık yükselir, <strong>hidrür</strong> işlemin merkeze doğru gidildikçe<br />

yavaşlamaktadır ancak duvar sınırlarında ve giriş bölümünde gelişmiştir. Nümerik<br />

sonuçlar sıvı akışının sistem sıcaklık profili üzerinde önemli etkisinin olduğunu<br />

göstermiştir. Ancak toplam <strong>hidrür</strong> formasyonunun akışkan akışı tarafından<br />

etkilenmediği bulunmuştur. Bazı bölgelerde akışkan akışı nedeniyle sıcaklık düşüşü<br />

olmuştur. Bu düşüş yatağın diğer tarafındaki sıcaklık artışı ile dengelediği<br />

göstermişlerdir.<br />

Bilgili (2003), yapılan çalışmada <strong>metal</strong> <strong>hidrür</strong>lerde hidrojen depolaması ile ilgili bir<br />

model oluşturulmuş, bu modelin depolayabileceği hidrojen oranları, bu oranların<br />

seçilen geometriye bağımlılığı ve reaktör içindeki sıcaklık dağılımlarını<br />

belirlenmiştir. Bu noktada <strong>metal</strong> <strong>hidrür</strong> ortamların hidrojen depolama<br />

karakteristiklerinin bağlı olduğu değişkenlerden biri de reaktör yatak sıcaklığıdır. Bu<br />

nedenle çalışmada, farklı bölge kesitlerde sıcaklık değişimleri incelenmiştir. Çalışma<br />

42


sonucunda H/M oranları soğutulan, düşük sıcaklıkta bulunan bölgelerde, daha<br />

yüksek değerlere sahip olduğu belirlenmiştir. Diğer bir ifadeyle soğutulan bölgelerde<br />

hidrojen soğurulması daha hızlı gerçekleşmiştir. Çalışmada üç farklı geometri ele<br />

alınarak matematiksel modellemeler yapılmıştır. Geometri-1 iç çapı 5 mm dış çapı<br />

25 mm ve yüksekliği 60 mm; Geometri-2 iç çapı 10 mm, dış çapı 20 mm ve<br />

yüksekliği 30 mm; Geometri-3 ise 10 mm iç çapında 17 mm dış çapında ve 50 mm<br />

yüksekliğe sahiptir. Üç farklı geometri için yapılan çözümlerde elde edilen soğurma<br />

olayının gerçekleştiği ilk 650 s zaman diliminde en hızlı soğurma işlemi geometri<br />

3’te gerçekleşmiştir. Bununla birlikte, hidrojenin <strong>metal</strong> yatak içine tamamen şarj<br />

olması için geçen süreler karşılaştırıldığında geometri 3’ün geometri 1’e oranla<br />

yaklaşık 2.8 kat daha hızlı sürede şarj olduğu belirlenmiştir. Geometri-2 ve<br />

Geometri-3 için bu fark ise çok düşüktür. Fakat <strong>metal</strong> <strong>hidrür</strong> yatak hacmi açısından<br />

bakıldığında Geometri-1’in Geometri-2’ye göre 4 kat daha büyük bir hacme sahip<br />

olması bu süreyi etkileyen önemli bir parametredir.<br />

Demiralp vd. (2004), yaptıkları çalışmada <strong>metal</strong> <strong>hidrür</strong> yataklarda ısı ve kütle<br />

transferini deneysel ve teorik olarak incelemişlerdir. Deney esnasında, LaNi5 alaşımı<br />

ile dolu reaktöre sabit basınç altında hidrojen şarj edilmiştir. Reaktördeki sıcaklık<br />

değişimleri farklı noktalarda ölçülerek, bilgisayara kaydedilmiştir. Deneyler sıcaklık<br />

değerleri oda sıcaklıklarına düşünceye kadar sürdürülmüştür. Teorik kısımda ise bu<br />

işlem esnasında kompleks ısı ve kütle transferini ve sıvı akışını göz önünde<br />

bulunduran iki boyutlu bir matematiksel model geliştirilerek nümerik olarak<br />

çözülmüştür. Elde edilen deneysel verilerle nümerik sonuçlar karşılaştırıldığında<br />

verilerin uyum içerisinde olduğu ve geliştirilen matematiksel modelin geçerliliği test<br />

edilmiştir. Deneylerde %99.999 saflıkta hidrojen ve oksidasyonu önleyebilmek için<br />

argon, istenmeyen gazları sistemden atmak için vakum pompası, öğütücü, reaktör,<br />

vakum ve basınç manometresi, termokupl ve veri kaydedici kullanılmıştır. Sistemde<br />

kullanılan reaktör 40 mm iç çapında (1 mm et kalınlığında), 120 mm yüksekliğinde<br />

olup, St 60 malzemeden yapılmıştır. 5000 s süre ile reaktördeki sıcaklık dağılımının<br />

değişik şarj basınçlarında zamana göre değişimi ölçülmüştür. Reaktörde ölçülen<br />

sıcaklıklar ilk 500 s içinde çok farklı değerlerde olmakla beraber zaman ilerledikçe<br />

birbirlerine çok yaklaşmakta ve ortam sıcaklığına ulaşmaktadır. Yüksek şarj<br />

43


asınçlarında ölçülen sıcaklık değerleri, düşük şarj basınçlarında ölçülen sıcaklık<br />

değerlerinden daha büyüktür. Buna sebebi, şarj basıncındaki artışın reaksiyon hızını<br />

arttırıcı yönde etki etmesidir. Dolayısıyla yüksek şarj basınçlarında ekzotermik<br />

reaksiyon daha hızlı olmakta, şarj basıncı azaldıkça reaksiyon hızı da azalmaktadır.<br />

Reaksiyon sonucunda oluşan ısının merkezde daha yüksek olduğu, merkezden cidara<br />

doğru gelindikçe ise azaldığı görülmektedir. Ayrıca z düzleminde artıkça sıcaklığın<br />

artığı, maksimum sıcaklık eğrileri ise z=60 mm’de olduğunu belirlemişlerdir.<br />

Oi vd. (2004), plaka-kanatçık tipi ısı değiştirici tabanlı <strong>metal</strong> <strong>hidrür</strong> reaktörün ısı<br />

transferi karakteristiğini deneysel ve nümerik olarak araştırmıştır. Hidrojen şarj<br />

işlemi sırasında deneysel ve nümerik sonuçlar yakın değerler ortaya koyarken, deşarj<br />

işlemi sırasında önemli farklılıklar tespit edilmiştir. Bu durum nümerik çalışmada<br />

kullanılan <strong>metal</strong> <strong>hidrür</strong> malzemenin ısıl iletkenlik katsayısı ve porozite değerlerinin<br />

gerçek değerlerinden uzak olarak kullanılmasına bağlanmıştır.<br />

Demircan vd. (2005), çalışmada iki adet LaNi5–H2 reaktöre hidrojen absorpsiyonu<br />

deneysel ve teorik olarak araştırılmıştır. Deneysel çalışmada iki tank LaNi5<br />

alaşımlarıyla dolduruluş ve sabit basınçta hidrojen şarj edilmiştir. Tanklardaki<br />

sıcaklık değişimi farklı noktalardan ölçülmüş ve bilgisayara kaydedilmiştir. Hidrür<br />

işlemi sıcaklık ölçülerek belirlenmiştir. Teorik çalışmada ise, iki boyutlu<br />

matematiksel model kompleks ısı ve kütle transferi ve akışkan akışı varsayımıyla<br />

nümerik olarak çözülmüştür. Matematiksel modelin geçerliliğini belirlemek için<br />

nümerik sonuçlar deneysel verilerle karşılaştırılmıştır. Deneysel veriler ve nümerik<br />

sonuçların karşılaştırıldığında kabul edilebilir sonuçlar elde edilmiştir.<br />

MacDonald ve Rowe (2006), yaptıkları çalışmalarında <strong>metal</strong> <strong>hidrür</strong> reaktör ile çevre<br />

arasında ısı transferine yoğunlaşmış ve <strong>metal</strong> <strong>hidrür</strong> reaktöre eklenen kanatçıkların<br />

etkisini incelemişlerdir. Ayrıca istenilen hidrojen miktarını elde edebilmek amacıyla<br />

kanatçıkların önemi, geliştirilen tek boyutlu rezistans analizi ve iki boyutlu transient<br />

model ile ele alınmıştır. Eklenen kanatçıkların reaktör içerisindeki hidrojen basıncı<br />

üzerinde önemli bir etkisinin olduğu saptanmıştır.<br />

44


Da Li vd. (2007), LaNi5−xAlx alaşımlarına Al ilavesinin hidrojen depolama hızına ve<br />

reaksiyon kinetiklerine etkisini belirlemek amacıyla deneysel bir çalışma<br />

yapmışlardır. Dört farklı LaNi5−xAlx (x=0.25, 0.50, 0.75, 1.00) alaşımının hidrojen<br />

depolama özelliklerini belirlemek amacıyla sistematik çalışmalar gerçekleştirilmiştir.<br />

Hidrojenerasyon öncesi ve sonrası mikro yapı özellikleri, aktivasyon, kinetikler ve<br />

termodinamik özellikleri, yanmazlık özellikleri gibi özellikler dört farklı alaşım<br />

yapısı için araştırılmıştır. Deney sonuçları alaşımlar mükemmel aktivasyon<br />

özellikleri ve kinetik özelliklere sahip olduğunu ortaya koymuştur. Alaşımların<br />

kristal yapısının x-ray kırınım analizleri alaşımlara 0≤x≤1 arasında Al ilavesinin<br />

alaşımların kristal yapısını değiştirmediğini göstermiştir, ancak kafes sabiti biraz<br />

artmıştır. Alaşım-H2 sisteminin termodinamik özellikleriyle entalpi mutlak değeri<br />

artarken, x değeriyle denge basıncı artarken, hidrojen depolama yeteneği ve<br />

histeresiz düşmüştür. Ayrıca x değerinin artmasıyla alaşımların hidrojen absorbe<br />

etme hızı düştüğünü belirtmişlerdir.<br />

Li vd. (2007), çalışmalarında dört farklı özellikte LaNi5−xAlx (x=0.25, 0.50, 0.75,<br />

1.00) alaşımlarının hidrojen depolama özellikleri incelenmiştir. Dört farklı Al<br />

ilavesinde hidrojenasyon öncesi ve sonrasında mikro yapı özellikleri, aktivasyon<br />

kinetikleri ve termodinamik özellikleri, anti yanma özellikleri de araştırılmıştır.<br />

Sonuçlar bu alaşımların çok iyi aktivasyon ve kinetik özelliklerine sahip olduklarını<br />

ortaya çıkarmıştır. Alaşımların kristal yapısının x-ray kırınım analizi alaşım içerisine<br />

0 ≤ x ≤ 1 aralığında Al ilavesinin alaşımların kristal yapısında değişikliğe neden<br />

olmadığını göstermiştir, fakat kafes sabiti biraz artmıştır. Alaşım-H2 sisteminin<br />

termodinamik özelliklerinin ölçülmesiyle, x değerinin artmasıyla entalpi mutlak<br />

değerinin artmasına rağmen; denge basıncı, hidrojen absorbe etme yeteneği ve<br />

histeresiz düştüğü belirlenmiştir. Ayrıca x değerinin artmasıyla alaşımların hidrojen<br />

absorbe etme hızının düştüğü bulunmuştur.<br />

Melloulia vd. (2007), deneysel çalışmalarında spiral ısı değitirici şeklinde bir reaktör<br />

geliştirmişlerdir. Yerleştirilen bu ısı değiştirici toplam depolama problemlerini<br />

azaltmada önemli katkılar sağlamıştır. Performans çalışmalarında basınç, tank hacmi,<br />

absorpsiyon sıcaklığı ve toplam ısı transfer katsayısı değerlerinde önemli gelişmeler<br />

45


sağlanmıştır. Yapılan deneysel çalışmalar sonucunda ısı değiştirici kullanımıyla<br />

hidrojen şarj/deşarj süresinde önemli ölçüde düştüğü görülmüştür. Bunlara ilave<br />

olarak soğutucu akışkan sıcaklığı ve debisi, uygulanan basınç ve hidrojen tankı<br />

hacmi gibi parametrelerin etkisi incelenmiş ve bu parametrelerin seçiminin önemi<br />

üzerinde durulmuştur.<br />

MacDonald vd. (2007), çalışmalarında <strong>metal</strong> <strong>hidrür</strong>lerin hidrojen depolama<br />

özellikleri deneysel ve sayısal modelleme kullanılarak belirlenmiştir. Özel olarak<br />

imal edilen <strong>metal</strong> <strong>hidrür</strong> tankında hidrojen giriş ve çıkış debi değerleri belirlenmiştir.<br />

Metal <strong>hidrür</strong> yatağında sıcaklık ölçümleri için radyal ve eksenel olarak termokupllar<br />

yerleştirilmiştir. Bu veriler önceden geliştirilmiş 2 boyutlu matematiksel modelde<br />

kullanılmıştır. Model parametrik çalışmada <strong>metal</strong> <strong>hidrür</strong> sistemlerinin bazı anahtar<br />

değişkenlerinin belirlemede kullanılmıştır. Bu değişkenler denge basıncı, yapılan<br />

çalışmaya bağlı olarak alaşım konsantrasyonu ve efektif ısı iletkenlik değeridir. Bu<br />

değerler konsantrasyona ve denge basıncına göre analiz edilmiştir. Denge basıncı<br />

eğimi, özellikle iç çevrimin doğruluğunu belirlemede önemli bir özelliktir olduğunu<br />

belirlemişlerdir.<br />

Mohan vd. (2007), yaptıkları çalışmada bir pratik <strong>metal</strong> <strong>hidrür</strong> esaslı hidrojen<br />

depolama cihazı tasarlamışlar. Bu cihaz hidrojen gazını dağıtmak için birçok filtre ve<br />

<strong>hidrür</strong> yatakta hidrojen şarj veya deşarjına bağlı olarak ısıtıcı veya soğutucu ısı<br />

değiştirici tüplerden oluşmaktadır. Bu çalışmada LaNi5 <strong>hidrür</strong> alaşımları kullanılan<br />

bir modelleme gösterilmiştir. Çalışmada geometrik ve çalışma parametrelerinin<br />

depolama cihazlarının hidrojen depolama verimlerine etkileri belirlenmiştir. Sonuçta<br />

tasarım parametrelerinin <strong>hidrür</strong> performansını önemli oranda etkilediği<br />

belirlemişlerdir.<br />

Kınacı (2007), yaptıkları çalışmada AB (A=Al, Be, Cu, Fe, Ni, Sb, V ve B=Ti) tipi<br />

ara <strong>metal</strong>lerin, toplam enerjileri, hidrojen depolama kapasiteleri ve kararlılıkları olası<br />

bir depolama malzemesi bulmak amacıyla incelenmiştir. Ayrıca, araştırılan<br />

malzemelerin termodinamik özellikleri ve atom-elektron yapıları arasındaki ilişkiler<br />

gösterilmiştir. Bu iş için bir temel prensipler psödo-potansiyel yöntemi olan<br />

46


yoğunluk fonksiyonlu teori (DFT) ve bu teori kapsamındaki genelleştirilmiş değişim<br />

ölçüsü yaklaşımı (GGA) kullanılmıştır. Hesaplamalar, CuTiH dışında, deneyle<br />

gözlemlenen yapıları doğru olarak tahmin etmiş atom ve hücre parametrelerini DFT<br />

için kabul edilen hata sınırları içinde bulmuştur. CuTi ara <strong>metal</strong> sisteminde CuTiH<br />

yapısında deneysel olarak gözlemlenenden belirgin derecede düşük oluşum<br />

enerjisine sahip bir yapı belirlenmiştir. Bu çelişki, deneysel yapının yarı-kararlı ve<br />

hidrojenin yapı içindeki hareketi için gereken aktivasyon enerjisinin yüksek<br />

oluşundan kaynaklanıyor olabilir. Çalışma sonucunda AlTi ve SbTi ara <strong>metal</strong>lerinin,<br />

bu malzemelerin <strong>hidrür</strong>lerinin kararsız olması nedeniyle hidrojen depolamaya uygun<br />

olmadığı görülmüştür. Araştırılan diğer sistemlerde <strong>hidrür</strong> kararlığı VTi, CuTi, NiTi,<br />

BeTi, FeTi sırasıyla azalmaktadır. VTi, FeTi ve NiTi sistemlerinde dört-<strong>hidrür</strong><br />

(MTiH4) oluşumunda, hidrojen¸ çevresindeki <strong>metal</strong> atomu düzeninde altı-kenarlıdan<br />

dört-kenarlıya geçiş öngörülmüştür. Bu kompozisyonda FeTi ve NiTi için belirlenen<br />

<strong>hidrür</strong>ler kararlı değildir. Ancak oluşum enerjilerinin sıfıra yakınlığı nedeniyle<br />

kimyasal kompozisyonda veya depolama parametrelerinde yapılacak uygun bir<br />

değişiklik <strong>hidrür</strong>lerin oluşmasını sağlayabilir. VTi depolama kapasitesine<br />

bakıldığında VTiH6’nın bile negatif oluşum enerjisine sahip olması nedeniyle umut<br />

vaat etmektedir. Diğer yandan oluşan <strong>hidrür</strong>lerin kararlılığı genellikle yüksektir ve<br />

bu nedenle hidrojeni serbest bırakması daha zor olacaktır.<br />

Muthukumar vd. (2007), çalışmada, içerisine MmNi4.6Al0.4 alaşımları doldurulmuş<br />

silindirik biçimli bir hidrojen depolama cihazında birleştirilmiş ısı ve kütle transfer<br />

işlemleri sayısal olarak ticari yazılım Fluent 6.1.22 kullanılarak gösterilmiştir. Cihaz<br />

için hidrojen depolama verimi, hidrojen şarj basıncı ve absorpsiyon sıcaklıkları gibi<br />

çalışma parametreleri değiştirilerek belirlenmiştir. Ayrıca, <strong>hidrür</strong> yatak kalınlıkları ve<br />

toplam ısı transfer katsayısı gibi farklı yatak parametrelerinin değiştirilmesiyle cihaz<br />

depolama verimi üzerine etkileri belirlenmiştir. Hidrür yatak ortalama depolama<br />

sıcaklığı ve farklı kaynak basınçlarında hidrojen depolama kapasiteleri literatürde<br />

bulunan değerlerlerle karşılaştırıldığında iyi sonuçlar alınmıştır. Başlangıçta<br />

<strong>hidrür</strong>leme işlemi esnasında hidrojen absorpsiyonu hızlı bir artış göstermiş daha<br />

sonra reaksiyon ısısından dolayı <strong>hidrür</strong> yatak sıcaklığının artmasıyla yavaş yavaş<br />

düşmüştür. Verilen herhangi bir absorpsiyon sıcaklığında hidrojen absorpsiyon oranı<br />

47


ve hidrojen depolama kapasitesi, şarj basıncının artmasıyla bulunmuştur. Hidrojen<br />

absorpsiyon kapasitesindeki değişim, kaynak basıncının 5–35 bar arasında 5 bar<br />

arttırılarak farklı değerlerde reaksiyon oranı ve sıcaklığı gösterilmiştir. Ayrıca toplam<br />

ısı transfer katsayısı 750 W/m 2 den 1250 W/m 2 ’ye ve soğutucu akışkan sıcaklığı<br />

288–298 K arasında değiştirilerek hidrojen depolama kapasitesi üzerindeki etkisi<br />

araştırılmıştır. Sonuçta ise ısı transfer oranı <strong>hidrür</strong>leşme oranını hızlandırdığı<br />

belirlenmiştir. 35 bar ve 298 K kaynak şartlarında MmNi4.6Al0.4 birim kg alaşımda<br />

yaklaşık 13.1 g hidrojen depolamıştır.<br />

Melnichuk vd. (2008), çalışmalarında, hidrojen absorpsiyonu performans testleri için<br />

bir küçük <strong>metal</strong> <strong>hidrür</strong> hidrojen depolama prototipi geliştirmişlerdir. Çalışmalarında<br />

<strong>metal</strong> <strong>hidrür</strong>lerin aktivasyon davranışı, desorpsiyon kapasitesi, desorpsiyon akış<br />

oranı ve ısı transferi üzerine yoğunlaşmıştır. AB5 inter<strong>metal</strong>ik tipi alaşımlar<br />

kullanılarak bir <strong>hidrür</strong> yatak yerleştirilmiştir. MmNi4.7Al0.3 kimyasal<br />

kompozisyonunda oda sıcaklıklarında denge basıncı yaklaşık 300 kPa olarak<br />

belirlenmiştir ve bu değer çalışma şartları için yeterlidir. Diğer yandan basınç<br />

değişimi, farklı debiler ve başlangıç şarj basınçları için sabit sıcaklıklarda deşarj<br />

testleri esnasında görüntülenmiştir. Deşarj esnasında dinamik basınç ölçümlerindeki<br />

azalma zamanın bir fonksiyonu olarak ölçülmüştür. Aynı zamanda, <strong>hidrür</strong> ayrışma<br />

kinetikleri ve ısı transfer şartlarıyla ilişkili olarak deşarj debisindeki artışla dinamik<br />

deşarj basıncındaki azalma bulunmuştur.<br />

Kim vd. (2008), çalışmalarında geliştirilmiş yüksek ısıl iletkenliğe sahip, gözenekli<br />

kompakt <strong>metal</strong> <strong>hidrür</strong> teknolojisi kullanarak, hidrojen kompresörlerinin<br />

performanslarını belirlemeyi amaçlamışlardır. Hidrojen depolama için soğutma<br />

ortamına girilen sıcaklıklar (20 ve 30 °C), sıkıştırma için ısıtma ortamı (60, 70, 80 ve<br />

90 °C) ve <strong>metal</strong> <strong>hidrür</strong> tipi (LaNi5 ve LaNi4.75Al0.25) anahtar parametreler olarak<br />

dikkate alınmıştır. Her bir koşulda sıcaklık dağılımı ve sıkıştırma basıncı<br />

kaydedilmiştir. Araştırmada LaNi4.75Al0.25 alaşımları LaNi5 alaşımından daha hızlı<br />

absorpsiyon hızı göstermiştir. Buna ilaveten sıkıştırma oranı esasen soğutma ve<br />

ısıtma arasındaki sıcaklık farkına bağlı olduğu; LaNi5 reaktörde LaNi4.75Al0.25<br />

reaktörden daha az sıkıştırma işi ve verimi gösterdiği bulunmuştur. Şartlar dikkate<br />

48


alındığında, 20 °C soğutma sıcaklığı ve 90 °C ısıtma sıcaklığında LaNi5<br />

çalıştırıldığında maksimum iş (2.44 kJ) ve verim (%5.91) gerçekleşmiştir.<br />

Cheng vd. (2008), yaptıkları bu çalışmalarında hidrojen şarj/deşarj çevriminin<br />

LaNi4.25Al0.75 hidrojen depolama özelliklerine etkisi incelenmiştir. Şarj/deşarj ve şarj<br />

kinetikleri için başlangıç aktivasyonu ve 1000 şarj/deşarj çevriminden sonra her ikisi<br />

için basınç izotermleri 363, 383 ve 403 K sıcaklıklarında ölçülmüştür. Uzun süre<br />

çevrimlerden sonra, şarj kinetiklerinde çok önemli bozulmalar gözlenirken plato<br />

basıncı ve depolama kapasitesi kinetiklerinde çok küçük değişimler gözlenmiştir.<br />

1000 çevrim sonunda başlangıç aktivasyonu sonrası (0.78–301 µm) ile<br />

karşılaştırıldığında tane büyüklüğü dağılımı daralmıştır (0.36–103 µm) bunun anlamı<br />

da hacim 34.81’den 11.34’e düşmüştür. 1000 çevrimden sonra LaNi4.25Al0.75 hidrojen<br />

depolama özelliklerindeki değişim LaNi4.25Al0.75 özellikleri, püskürtme olgusu ve<br />

çevrimlerde hidrojen şarjından dolayı kirlenme disproporsiyonu ile açıklanmıştır.<br />

Muthukumar ve Manvendra (2009), MmNi4·6Al0·4 içeren silindirik <strong>metal</strong> <strong>hidrür</strong><br />

reaktörde birleştirilmiş ısı ve kütle transferi bir iki boyutlu nümerik analizle<br />

göstermişlerdir. Hidrojen absorpsiyonunu anlama, enerji, momentum ve kütlenin<br />

korunumu ve reaksiyon kinetikleri denklem için ana denklemler, sonlu hacimler<br />

yöntemi eş zamanlı olarak kullanılarak çözülmüştür. MmNi4·6Al0·4 esaslı hidrojen<br />

depolama cihazları üzerindeki verim araştırmaları, hidrojen kaynak basıncı,<br />

absorpsiyon (soğutucu akışkan) sıcaklığı, toplam ısı transfer katsayısı ve <strong>hidrür</strong> yatak<br />

kalınlığına bağlı değiştiği ortaya çıkmıştır. Hidrojen depolama verimi üzerinde enerji<br />

denklemindeki konveksiyon etkisi teriminin göz ardı edilebilir olduğu bulunmuştur.<br />

Sonuçlar bilgisayar simülasyonundan elde edilen değerlerle deneysel verilerin<br />

uyumlu olduğunu göstermiştir. 30 bar ve 298 K kaynak şartlarında, MmNi4·6Al0·4<br />

depolama yaklaşık ağırlıkça %1.28 değerindedir, bu değer deneysel %1.3 verilere<br />

çok yakındır. Toplam ısı transfer katsayısının absorpsiyon zamanını önemli oranda<br />

düşürdüğünü belirlemişlerdir.<br />

Park vd. (2009), çalışmalarında ısıl enerji uygulamalarında birleştirilmiş AB5 tipi<br />

yüksek güçlü <strong>metal</strong> <strong>hidrür</strong> reaktörlerini araştırmışlardır. Her biri sırasıyla düşük ve<br />

49


yüksek basınçlı <strong>metal</strong> <strong>hidrür</strong> içeren LaNixAl5-x veya CayMm1-yNi5 alaşımlarından<br />

birini içeren iki reaktörün birleştirilmesi ile test reaktörleri oluşturulmuştur. Yüksek<br />

güçlü reaktörleri geliştirmek amacıyla, <strong>metal</strong> <strong>hidrür</strong>ler çok ince (yaklaşık 1–2 mm<br />

kalınlığında) bakırla kaplanmıştır ve gözenekli <strong>metal</strong> <strong>hidrür</strong>ler yoğunlaşıncaya kadar<br />

sıkıştırılmıştır. Deneyler esnasında, sistem reaktif kinetiklerinin dinamik davranışları<br />

gözlenmiştir. Test edilen sistemler arasında, birleştirilmiş Ca0.6Mm0.4Ni5 ve<br />

LaNi4.75Al0.25 ısıl enerji uygulamaları için çok verimli olduğu görülmüştür. Hidrür ve<br />

de<strong>hidrür</strong>leme işleminin son anında denge durumuna (yaklaşık olarak 250 s) ve<br />

üretim/absorpsiyon ısısının en yüksek miktarına ulaşmıştır. Alaşımda y’nin en düşük<br />

değerinde Ca1-yMmyNi5 alaşımları, daha hızlı reaksiyon kinetiklerinde daha fazla<br />

kalsiyum içeriğine neden olur. LaNixAl5-x reaktöründeki durumda ise, alüminyum<br />

ilavesi reaksiyon kinetiklerinin iyileşmesine neden olur. Ayrıca, düşük basınç <strong>hidrür</strong><br />

reaktör LaNixAl5-x, yüksek basınçlı <strong>metal</strong> <strong>hidrür</strong> Ca1-yMmyNi5 kullanıldığı zaman<br />

denge durumuna daha uzun zamanda ulaşıldığı sonucuna ulaşmışlardır.<br />

Satheesh vd. (2009), çalışmalarında MmNi4.6Al0.4/ MmNi4.6Fe0.4 <strong>hidrür</strong> çifti ile<br />

çalışan bir <strong>metal</strong> <strong>hidrür</strong> soğutma sisteminin sayısal analizi gösterilmiştir.<br />

Birleştirilmiş reaktör çiftinde zamana bağlı ısı ve kütle transferleri için bir silindirik<br />

model oluşturularak, kararsız, iki boyutlu matematiksel modelde nümerik olarak<br />

çözülmüştür. Denklem sistemleri ful implisit sonlu hacimler yöntemi (SHM)<br />

tarafından çözülmüştür. Reaktörlerde, reaksiyon yatak sıcaklık dağılımı, hidrojen<br />

konsantrasyonu ve denge basıncının, sabit ve değişken duvar sıcaklıkları sınır<br />

şartlarına etkisi araştırılmıştır. Basınç-konsantrasyon-sıcaklık çizimi bir dinamik<br />

korelasyon gösterilmiştir. Verilen 363/298/278 K sıcaklıklarında (TH/TM/TC), çift ve<br />

tek etkili <strong>metal</strong> <strong>hidrür</strong> sistemlerinde sabit ve değişken duvar sıcaklık sınır şartlarında,<br />

çevrim zamanı sırasıyla 1470 ve 1765.6 s olarak bulunmuştur. Sayısal sonuçlar<br />

LaNi4.61Mn0.26Al0.13/La0.6Y0.4Ni4.8Mn0.2 <strong>hidrür</strong> çifti için literatürde gösterilen<br />

sonuçlarlarla karşılaştırılmış ve uyum içerisinde olduğunu belirtmişlerdir.<br />

Kaplan (2009), yaptığı çalışmada farklı dizaynlardaki <strong>metal</strong> <strong>hidrür</strong> reaktörlerde şarj<br />

işlemi üzerindeki ısı transfer mekanizmasının etkisi farklı hidrojen şarj basınçları<br />

altında araştırılmıştır. Üç farklı silindirik reaktör aynı temel boyutlarda dizayn<br />

50


edilmiş ve üretilmiştir. Reaktör-1 20 mm iç çapında, 122 mm dış uzunluk ve 1 mm<br />

kalınlıktadır. Reaktör-2; Reaktör-1 ile aynı boyutlardadır fakat reaktörden daha iyi ısı<br />

uzaklaştırabilmek için 42 mm çapında ve 4.5 mm aralıklarla yerleştirilmiş 1 mm<br />

kalınlığında 22 dairesel kanatçık bulunmaktadır. Reaktör-3’te ise zorlanmış<br />

konveksiyonla ısı taşınımı sağlamak için reaktör çevresinde su sirkülâsyonu<br />

yapılmıştır. Şarj esnasında farklı bölgelerden 1–10 bar basınç aralığında reaktör<br />

sıcaklıkları ölçülmüştür. Üçüncü reaktör bütün çalışma basınçlarında hızlı şarj<br />

zamanıyla en düşük sıcaklık artışı göstermiştir. Reaktörlerin efektif ısı transfer<br />

katsayıları deneysel sonuçlara göre hesaplanmıştır ve bunlar sırasıyla 5.5 ±1 W/m 2 K,<br />

35 ± 1 W/m 2 K ve 113 ± 1 W/m 2 K olarak bulunmuştur. Deneysel sonuçlar <strong>hidrür</strong><br />

reaktör şarjında genel olarak ısı transferine bağlı olduğu ve daha iyi soğutulmuş<br />

reaktörlerin hızlı şarj özellikleri gösterdiğinin ortaya koymuştur.<br />

Askri vd. (2009), <strong>metal</strong> <strong>hidrür</strong> depolama tanklarında hidrojen depolama zamanı<br />

üzerinde en büyük faktörün reaksiyon yatağından ısının hızlı şekilde taşınması<br />

olduğunu gösterdiği çalışmasında; Đki boyutlu bir model geliştirerek bu modelin<br />

geçerliliğini deneysel çalışmalarla göstermiştir. Bu model ilk olarak absorpsiyon<br />

işlemi esnasında tank yüzeylerinin ısıl yoğunluğunu hesaplamada kullanılmıştır.<br />

Yatak duvarında çelik ve pirinç malzeme kullanılması durumları test edilmiş ve bu<br />

deneyler depolama zamanı üzerinde önemli bir etkisinin olmadığını göstermiştir.<br />

Daha sonra kurulmuş model farklı <strong>metal</strong> <strong>hidrür</strong> reaktör dizaynı için dinamik<br />

davranışlar araştırılmıştır. i) Bir silindirik tank, ii) dıştan kanatçıklı silindirik tank,<br />

iii) konsantrik boru ile içten soğutulmuş bir silindirik tank ve iv) kanatçıklı<br />

konsantrik borularla donatılmış silindirik tank olmak üzere dört farklı tasarım<br />

yapılmıştır. Yapılan optimizasyon sonucunda hidrojen depolama zamanında yaklaşık<br />

%80 gelişme elde edilmiştir.<br />

Baldissin ve Lombardo (2009), LaNi4.8Al0.2 alaşımlarında hidrojen absorpsiyonu ve<br />

desorpsiyonu için iki boyutlu bir matematiksel model geliştirmişler ve modelin<br />

geçerliliğini deneysel olarak ispatlamışlardır. Model bir enerji, kütle ve bir<br />

momentum dengesinden oluşmuştur. Bu diferansiyel denklemler sonlu elemanlar<br />

yöntemi kullanılarak Comsol yazılımında nümerik olarak çözülmüştür. Teorik ve<br />

51


deneysel sonuçlar karşılaştırıldığında <strong>metal</strong> <strong>hidrür</strong>lerde hidrojen depolamanın en<br />

önemli tarafı kullanılan reaktörlerin ısıl karakteristikleri, ısı değiştirici özellikleri:<br />

toplam ısı transfer katsayısı h ve efektif ısıl etkenliği k en önemli parametrelerdir.<br />

Deneysel çalışmada kullanılan cihaz paslanmaz çelikten imal edilmiş ve sıcaklığın<br />

sabit bir değerde tutulması için su banyosu içerisindedir. Reaktör içerisine<br />

LaNi4.8Al0.2 malzeme doldurulmuştur. Đki farklı reaktör kullanılarak yapılan<br />

deneylerde Reaktör 1’de ekzotermik reaksiyon sonucu oluşan maksimum sıcaklık<br />

yaklaşık 26 °C’ye ulaşmıştır. Kanatçıklı olarak imal edilen Reaktör 2 ise su<br />

banyosuna daldırılmış durumdadır ve burada da maksimum sıcaklık yaklaşık 24 °C<br />

olarak ölçülmüştür. Deneylerde ölçülen değerler matematiksel modelle<br />

karşılaştırıldığında uyum içerisinde oldukları belirlenmiştir.<br />

Mellouli vd. (2010), yaptıkları çalışmada, yakıt pili kullanan araçlar için, kanatlı<br />

spiral borularla donatılmış ısı değiştiricili <strong>metal</strong> <strong>hidrür</strong> hidrojen depolama tanklarında<br />

ısı ve kütle transferinin optimizasyonu iki boyutlu bir matematiksel model<br />

geliştirilmiştir. Model nümerik olarak çözülmüş kompleks ısı ve kütle transferini de<br />

hesaba katarak deneysel verilerle karşılaştırılmış ve iyi sonuçlar elde edilmiştir.<br />

Uzunluk, kalınlık ve tasarlanan kanatçıkların <strong>metal</strong> <strong>hidrür</strong> reaktörlerin şarj işlemi<br />

üzerine etkisi ortaya çıkarılmıştır. Buna ilaveten, oluşturulan model farklı şekillerde<br />

dizayn edilmiş <strong>metal</strong> <strong>hidrür</strong> reaktörlerin dinamik davranışlarını belirlemek için<br />

kullanılmıştır. Ayrıca spiral kanatçıklı borulara ısı değiştirici ilave edilerek soğutma<br />

dizayn seçenekleri araştırılmıştır.<br />

Dhaou vd. (2010), çalışmalarında <strong>metal</strong> <strong>hidrür</strong> tanklarda hidrojen reaksiyon zamanı<br />

(kısaca MHV) önemli derecede <strong>hidrür</strong> yataktan ısı transferi yapılmasına bağlı<br />

olduğunu göstermişlerdir. Bu çalışmada, <strong>metal</strong> <strong>hidrür</strong> tanklara hidrojen şarj işlemi<br />

performansı üzerine kanatçıklı spiral ısı değiştiricilerin geometrisi ve çalışma<br />

parametrelerinin etkisi deneysel olarak araştırılmıştır. Deneysel sonuçlar kanatçıklı<br />

spiral ısı değiştirici kullanıldığı zaman reaktör şarj zamanının oldukça düştüğünü<br />

göstermiştir. Ayrıca farklı parametrelerin etkisi (soğutucu akışkanın sıcaklığı ve<br />

debisi, uygulanan basınç ve hidrojen tank hacmi) tartışılmıştır ve elde edilen<br />

sonuçlarla seçilen bu parametrelerin çok iyi sonuçlar verdiğini belirlemişlerdir.<br />

52


3. MATERYAL VE YÖNTEM<br />

Yapılan bu çalışmada hidrojen şarj/deşarj işlemlerine etki eden proses<br />

parametrelerinden; hidrojen basıncı, debisi ve <strong>metal</strong> <strong>hidrür</strong> reaktör geometrisinin,<br />

reaktör ile çevre arasındaki ısı transferine ve hidrojen depolama özellikleri<br />

üzerindeki etkisinin belirlenmesi amaçlanmıştır. Diğer yandan kullanılan malzeme ve<br />

çalışma şartları gibi parametrelerinde depolanan hidrojen kütlesine etkisi deneysel<br />

olarak araştırılacaktır. Yapılan literatür araştırmasında LaNi5 alaşımlarına Al ilavesi<br />

alaşımın ısı iletim katsayısını ve gözenekliliğini arttırdığı, bunun sonucu olarak<br />

ekzotermik reaksiyonları hızlandırdığı için yapılan çalışmada <strong>hidrür</strong> alaşımı olarak<br />

LaNi4.75Al0.25 alaşımı tercih edilmiştir. Ayrıca deneylerde şarj basıncı 2, 4, 6 ve 8 bar<br />

basınçlara ayarlanarak hidrojen şarj basıncın reaksiyonlara ve depolanan hidrojen<br />

kütlesine etkisi ortaya konacaktır. Daha sonra ileri beslemeli yapay sinir ağları<br />

kullanılarak oluşturulan YSA modeli depolanan hidrojen kütlesinin basınç sıcaklık<br />

ve zamana bağlı olarak değişimi belirlenecektir. Son olarak deneysel sonuçlarla YSA<br />

kullanılarak elde edilen veriler karşılaştırılarak oluşturulan modelin geçerliliği ortaya<br />

konacaktır.<br />

3.1. Materyal<br />

3.1.1. Metal <strong>hidrür</strong> reaktör<br />

Reaktör geometrisinin hidrojen şarj/deşarj proseslerine etkisini belirlemek amacıyla<br />

Şekil 3.1’de gösterildiği gibi Reaktör 1 ve Reaktör 2 olarak adlandırılan iki farklı<br />

reaktör dizayn ve imal edilmiştir. Her iki reaktörde St 42 çelikten imal edilmiştir.<br />

Reaktör 1 kanatçıksız olarak dizayn edilmiş, 20 mm iç çapında, 125 mm dış<br />

uzunluğunda ve 2 mm et kalınlığındadır. Reaktör 2, Reaktör 1 ile aynı boyutlara<br />

sahiptir. Fakat reaktör içerisinde oluşan reaksiyon ısısını daha hızlı uzaklaştırmak<br />

için 4 mm aralıklarla yerleştirilmiş 40 mm dış çapında 1 mm et kalınlığında 23 adet<br />

dairesel kanatçık yerleştirilmiştir (Bkz. Şekil 3.2). Her iki reaktörde hidrojen<br />

şarjından sonra <strong>hidrür</strong> alaşımlarının hacimleri artacağı için yaklaşık %70 oranında<br />

135 g ağırlığında LaNi4.75Al0.25 malzeme ile doldurulmuştur.<br />

53


Şekil 3.1. Reaktörlerin kesit görünüşleri, a) Reaktör 1 ve b) Reaktör 2<br />

Şekil 3.2. Reaktör 1 ve Reaktör 2 kesit görünüşleri ve imal ölçüleri<br />

54


Şekil 3.3. Termokuplların <strong>metal</strong> <strong>hidrür</strong> reaktörler üzerine yerleşimi<br />

Metal <strong>hidrür</strong> reaktörlerde sıcaklık dağılımını belirlemek amacıyla reaktörlerin farklı<br />

noktalarına termokupllar yerleştirilmiştir. Termokupllar reaktör merkezinden z<br />

ekseni boyunca 5, 20, 40 ve 60 mm noktalarına yerleştirilmiştir ve reaktör üzerinde<br />

yerleştirildiği noktalar Şekil 3.3’te gösterilmiştir. Ayrıca termokuplların reaktör<br />

üzerine bağlanışı ve ölçümler Ek-2’de verilmiştir.<br />

3.1.2. Mekanik öğütücü<br />

Mekanik öğütücü 60 mm çapında ve 80 mm yüksekliğinde olup Şekil 3.4 ’te<br />

gösterilmiştir. Alaşımlar reaktör içerisine doldurulmadan önce tane boyutlarını<br />

küçültmek için öğütme işlemine tabi tutulmuştur. Öğütme işleminde her bir<br />

öğütmede karıştırıcı içerisine 1/10 oranında çelik bilye ve yaklaşık 100 g<br />

LaNi4.75Al0.25 ilave edilmiştir. Deneylerde kullanılan LaNi4.75Al0.25 alaşımları<br />

Ergenics Inc. Firmasından temin edilmiştir. LaNi alaşımlarının fiziksel ve kimyasal<br />

özellikleri Ek-1.1’de verilmiştir. Literatürde yapılmış çalışmalarda incelenerek<br />

öğütme işleminde öğütücü hızı 500 d/d olarak ayarlanmıştır. Öğütme işlemi 4 saat<br />

boyunca yapılmış ve her iki reaktör için aynı tane boyutlarında (yaklaşık 40 µm) aynı<br />

55


sürelerde öğütülmüş alaşım kullanılmıştır. Ayrıca öğütücü çevresinde soğutucu<br />

akışkan olarak su dolaştırılarak öğütücünün soğutulması ve <strong>hidrür</strong> yapısının<br />

bozulmaması sağlanmıştır.<br />

Şekil 3.4. Metal <strong>hidrür</strong> tozlarını öğütmek için kullanılan mekanik öğütücü<br />

Öğütme hızı hem <strong>hidrür</strong> miktarını hem de <strong>hidrür</strong>leşme hızını etkilemektedir. Hidrür<br />

miktarına olan etkisi, hazne içindeki sıcaklıkla ilgilidir. Hidrürleşmenin oluşabilmesi<br />

için yüksek bir sıcaklık gerekmektedir. LaNi4.75Al0.25 havayla temas eder etmez<br />

oksitlendiği için, <strong>hidrür</strong>leşmenin gerçekleşeceği aktif yüzey yok olmaktadır. Bu<br />

yüzden öncelikle bu oksit tabakasının kırılması gerekmektedir. Bunun için de yüksek<br />

sıcaklığa ihtiyaç vardır. Kullanılan yöntemin en önemli özelliklerinden biri<br />

<strong>hidrür</strong>leşme için gereken bu yüksek sıcaklığı elde etmek için herhangi bir ısıl işlem<br />

uygulanmamasıdır. Hidrür oluşumu için gerekli olan sıcaklık, öğütme sırasında<br />

hazne içerisindeki çarpışmalardan elde edilmektedir. Bu yüzden içerideki sıcaklığı<br />

artıran her etken <strong>hidrür</strong> oluşumunu hızlandırır. Hidrürleşmenin ilk başladığı yerler<br />

kusurlar ve boşluklar olduğundan, deformasyonla oluşan kusurların ve boşlukların<br />

sayısı arttığında <strong>hidrür</strong>leşme hızı da artmış olur. Bunun yanı sıra hız arttığında<br />

56


malzeme daha küçük parçalara ayrılır yani parçacık boyutu küçülür. Boyutu<br />

küçüldükçe yüzey/hacim oranı artar.<br />

Öğütme işleminde diğer bir önemli bir parametre de öğütme süresidir. Öğütme süresi<br />

arttıkça parçacık büyüklüğü azalır. Ancak kullanılan ve sonuçta elde edilmek istenen<br />

malzemenin yapısına uygun olan bir öğütme süresi vardır. Bu süre malzemenin<br />

hazırlanma koşularına da bağlıdır. Uygun süreden daha kısa ya da daha uzun yapılan<br />

öğütmeler sonucunda istenen yapı elde edilemez. Ayrıca uzun süreli öğütme, enerji<br />

tüketimine ve öğütme aletlerinin aşınması sonucu meydana gelebilecek kirlenme<br />

riskinin artmasına neden olmaktadır.<br />

3.1.3. Split fırın<br />

Reaktörlerde hidrojen şarjından önce aktivasyon işleminin gerçekleşmesi için reaktör<br />

içerisinin tamamen nemden ve havadan arındırılması gerekmektedir. Hidrür<br />

reaktörlerini aktive etmek için öncelikle Şekil 3.5’te gösterilen fırın içerisinde 220 °C<br />

sıcaklığa ısıtılarak 2 saat boyunca vakum yapılmıştır.<br />

Şekil 3.5. Split fırın<br />

57


Aktivasyon işleminde her bir hidrojen şarjından sonra reaktör, fırın içerisinde<br />

ısıtılarak vakum yapılmıştır. Deneyler esnasında da hidrojen şarjından sonra<br />

endotermik reaksiyon için gerekli ısı (deneylerde 200–220 °C arasında<br />

belirlenmiştir) fırın kullanılarak sağlanmıştır.<br />

3.1.4. Vakum pompası<br />

Reaktörlerde <strong>hidrür</strong> tozlarının aktivasyonu için reaktörün yaklaşık 220 °C ye kadar<br />

ısıtılıp vakum yapılması gerekir. Fırın içerisinde ısıtılan reaktörden çıkan nem ve<br />

hava vakum yapılarak atmosfere atılmıştır. Bu işlem aktivasyonun tam olarak<br />

gerçekleşebilmesi için deneylerden önce 8 bar basınç altında sekiz defa tekrar<br />

edilmiştir. Hidrojen şarjından sonra reaktör tekrar fırın içerisine yerleştirilerek Şekil<br />

3.6’da gösterilen vakum pompasıyla vakumlama işlemi yapılmıştır. Aktivasyon<br />

işleminin sonunda reaktör hidrojen depolayabilme yeteneğine sahip olmaktadır.<br />

Aktivasyon işleminden sonra her bir deneyin sonunda depolanan hidrojen<br />

boşaltılmasında da yine vakum pompası kullanılmıştır. Vakum pompasının teknik<br />

özellikleri Ek-1.4’te verilmiştir.<br />

Şekil 3.6. Vakum pompası<br />

58


3.1.5. Sıcaklık tarayıcı<br />

Sıcaklık tarayıcı farklı noktalardan aynı anda sıcaklık kaydı yapmak için<br />

kullanılmıştır. Şekil 3.7’de verilen sıcaklık tarayıcı aynı anda 8 farklı noktadan<br />

sıcaklık verisi alabilmekte ve aldığı veriyi bir USB kablosu ile veri kaydedici<br />

bilgisayara aktarabilmektedir. Yapılan deneylerde ekzotermik reaksiyon sonucu<br />

oluşan sıcaklık değişimi reaktörlerin z=5, 20, 40 ve 60 mm noktalarına yerleştirilmiş<br />

4 adet K tipi termokupl kullanılarak 4000 s boyunca kaydedilmiştir.<br />

3.1.6. Deney tesisatı<br />

Şekil 3.7. Sıcaklık tarayıcı<br />

Metal <strong>hidrür</strong> reaktörlerde tasarımın, reaksiyon hızının, depolanan hidrojen miktarının<br />

ve LaNi4.75Al0.25 alaşımlarının hidrojen absorpsiyon özelliklerinin belirlenmesi için<br />

oluşturulan deney düzeneği Şekil 3.8’de gösterilmiştir. Düzenek bir reaktör, ince<br />

partikül üretmek için çelik bilyeli mekanik öğütücü, deneylerden önce reaktör<br />

içerisindeki nem ve havayı boşaltmak için bir vakum pompası, reaktör içerisinde<br />

kalan nemi ve havayı buharlaştırmak için fırın, manometreler, termokupllar, veri<br />

toplama sistemi, hidrojen basıncını ayarlamak için regülatör, %99.999 saflıkta<br />

hidrojen içeren bir tank ve atmosfer şartlarının altında çalışma sağlamak için azot<br />

tankı bulunmaktadır.<br />

59


3.1.7. Hassas terazi<br />

Şekil 3.8. Deney tesisatı<br />

Deneylerde zaman göre depolanan hidrojen miktarının belirlenmesi amacıyla şekil<br />

3.9’da gösterlen hassas terazi kullanılmıştır. Terazi maksimum 3100 g ağırlık<br />

tartabilmekte, hassasiyeti 0,01 g, doğrusallık ± 0.02 g, stabilizasyon süresi 2 s ve 10–<br />

30 °C sıcaklık aralığında hassas ölçüm yapabilmektedir. Terazi bir RS–232<br />

bağlantısı ile bilgisayara bağlanabilmekte ve hassas olarak veriler alınabilmektedir.<br />

Şekil 3.9. Hassas terazi<br />

60


3.1.8. Hidrojen gazı<br />

Hidrojen 1500’lü yıllarda keşfedilmiş, 1700’lü yıllarda yanabilme özelliğinin farkına<br />

varılmış, evrenin en basit ve en çok bulunan elementi olup, renksiz, kokusuz,<br />

havadan 14.4 kez daha hafif ve tamamen zehirsiz bir gazdır. Atomik sembolü H olan<br />

hidrojenin atom ağırlığı 1.00797, atom sayısı 1 olup en basit ve en hafif elementtir.<br />

Hidrojen doğada en çok bulunan element olmasına rağmen, hafifliği sebebi ile<br />

atmosfere yükselip orada serbest kaldığından, yeryüzünde serbest halde çok az<br />

bulunur. Atmosferin çok küçük bir bileşenidir ve yeryüzünde yaklaşık 0.5 ppm<br />

bolluktadır ve 2500 km’nin üzerindeki atmosferin hemen hemen tamamı hidrojendir.<br />

Güneş ve diğer yıldızların termonükleer tepkimeyle vermiş olduğu ısının yakıtı<br />

hidrojen olup, evrenin temel enerji kaynağıdır. Evrenin tamamındaki atomların<br />

%90’ı ve evren kütlesinin %75’i hidrojendir. Hidrojen gazının termofiziksel<br />

özellikler ayrıntılı olarak Ek-1.3’te verilmiştir.<br />

3.1.9. Azot gazı<br />

Azot, daha çok havanın sıvılaştırılması ve ayrıştırılması yöntemi ile elde edilir. Sıvı<br />

ve/veya gaz fazlarında ticari olarak bulunur. Sıvı fazda; çift cidarlı, vakum ve perlit<br />

malzemesi ile yalıtılmış özel kriyojenik tanklar içinde depolanır ve nakledilir. Gaz<br />

fazda; basınç altında boru hattı ile veya başınca dayanıklı, dikişsiz çelik tüpler içinde<br />

sıkıştırılmış olarak tedarik edilir. TSE standartlarına göre sıvı/gaz azotun birinci<br />

sınıfının safiyeti minimum %99.999, ikinci sınıfının safiyeti minimum %99.99<br />

olmalıdır. Deneylerde <strong>metal</strong> <strong>hidrür</strong> tozları öğütülürken ve reaktörlere doldurma<br />

işlemi esnasında öğütücü içerisine sürekli olarak azot gazı gönderilerek alaşımların<br />

havayla temas etmesi ve oksidasyonu önlenmiştir. Bütün öğütme süresince ve<br />

alaşımın reaktör içerisine doldurulması esnasında inert gaz olarak azot gazı<br />

kullanılmıştır. Azot gazının termofiziksel özellikleri Ek-1.2’de ayrıntılı olarak<br />

verilmiştir.<br />

61


3.2. YÖNTEM<br />

3.2.1. Metal <strong>hidrür</strong> reaktör kanat etkenliği<br />

Kanat kullanımı, bir yüzeyden ısı geçişini artırmak içim etkin yüzey alanını artırmayı<br />

amaçlar. Bununla birlikte kanadın kendisi orijinal yüzeyden ısı geçişine bir iletim<br />

direnci gösterir. Bu nedenle, kanat kullanımının ısı geçişini mutlaka artıracağı<br />

önceden söylenemez. Bu husus kanat etkenliği tanımı tanımlanarak<br />

değerlendirilebilir. Kanat etkenliği, εf, kanatlı halde geçen ısının kanatsız halde<br />

geçebilecek ısıya oranı olarak tanımlanır. Böylece,<br />

olup Ac,b tabandaki kanat kesit alanıdır.<br />

q<br />

f<br />

ε f =<br />

(3.1)<br />

hA cθ<br />

b<br />

Gerçekçi her tasarımda εf’nin değeri mümkün olduğunca büyük olmalıdır ve genel<br />

olarak, εf ≥ 2 olmadıkça kanat kullanımı uygun değildir. Kanatların yerleştirilme<br />

düzeni ısı taşınım katsayısını değiştirebilir ancak bu etki genellikle göz ardı edilir. Bu<br />

nedenle, kanatlı yüzeyin taşınım katsayısının kanatsız yüzeyinkine eşit sayılması<br />

durumunda sonsuz kanat yaklaşımı için;<br />

k P<br />

ε f =<br />

(3.2)<br />

h A<br />

Aynı denklem, ısı taşınım katsayısının küçük olduğu durumlarda kanat kullanmanın<br />

yararını da göstermektedir. Akışkan gaz olduğunda ve özellikle yüzeyden ısı geçişi<br />

doğal taşınımla oluğunda kanatlara daha çok gerek duyulmaktadır. Kanatlar bir gaz<br />

ile bir sıvıyı ayıran yüzeylerde kullanılacaksa, genellikle düşük ısı taşınım<br />

katsayısının bulunduğu taraf olan gaz tarafına yerleştirilir. Bilinen bir örnek otomobil<br />

radyatörünün borularıdır. Kanatlar, içinde su akışı olan (büyük h) boruların iç<br />

yüzeyine değil, üzerinden ortam havasının aktığı (küçük h) dış yüzeyin uygulanır.<br />

εf>2, kanat uygulamasını haklı kılacak bir kıstas olarak kullanılacaksa, (kP/hAc) > 4<br />

c<br />

62


olması gerekecektir. Kanadın L=2.3/m’den fazla uzun olması anlamsızdır (Incropera<br />

and Dewitt, 2006).<br />

3.2.1.1. Kanat verimi<br />

Kanat ısıl performansının bir diğer ölçüsü kanat verimi ηf’dir. Taşınım için sıcaklık<br />

farkı, dip (x=0) ve akışkan arasındaki sıcaklık farkıdır, θb=Tb - T∞. Bu nedenle bir<br />

kanadın yayabileceği enerjinin en yüksek değeri bütün kanat yüzeyi taban<br />

sıcaklığında olduğu zaman gerçekleşecektir. Ancak bu ideal bir durumdur ve kanat<br />

içinde bir sıcaklık eğişimi her zaman vardır. Buradan yola çıkarak kanat verimi;<br />

q<br />

f<br />

f<br />

η f = =<br />

(3.3)<br />

q max h A f θb<br />

biçiminde tanımlanabilir. Burada Af kanadın yüzey alanıdır.<br />

Şekil 3.10. Dikdörtgen profilli dairesel kanat verimleri<br />

63<br />

q


Adyabatik uçlu ve sabit kesit alanlı düz bir kanat için kanat verimi aşağıda<br />

verilmiştir.<br />

hPkA c θb<br />

tanh mL tanh m L<br />

η f =<br />

=<br />

(3.4)<br />

h P Lθ<br />

m L<br />

b<br />

Kanat veriminin en yüksek ve en düşük değerleri olan 1 ve 0 değerlerine sırasıyla, L<br />

uzunluğu 0 ve ∞’a yaklaşıldığında erişilmektedir.<br />

Dikdörtgen profilli dairsel kanat verimleri Şekil 3.10’da verilmiştir. Ucunda ısı kaybı<br />

olan düz dikdörtgen bir kanattan ısı geçişini veren denklemi kullanmak oldukça<br />

zordur. Bu denklemi kullanmak yerine, dikdörtgen kanatlar için Lc=L + (t/2) ve iğne<br />

kanat için Lc=L + (D/4) biçiminde düzeltilmiş kanat yüksekliği tanımları yapılarak,<br />

adyabatik uç için geçerli olan denklem kullanılabilir. Düzeltme, taşınım uçlu gerçek<br />

kanattan ısı geçişi ile adyabatik uçlu daha uzun bir kanadın eşdeğer davranış<br />

gösterdiği varsayımına dayanmaktadır. Bu nedenle taşınım uçlu kanattan ısı geçişi;<br />

şeklindedir. Buradan kanatçık verimi;<br />

bağıntılarıyla hesaplanabilir.<br />

3.2.1.2. Toplam yüzey verimi<br />

q = h P k A θ tanh L<br />

(3.5)<br />

f<br />

c<br />

b<br />

64<br />

c<br />

tanh mL<br />

f<br />

mL<br />

= η (3.6)<br />

Tek bir kanadın ısıl davranışını gösteren kanat verimi ηf yerine, bir kanat dizisi ve<br />

üzerine yerleştirildiği yüzeyin ısıl davranışını gösteren toplam yüzey verimi ηo<br />

kullanılabilir. Şekil 3.11’de verilen örnek kanat dizilerinde S kanat hatvesidir.


Şekil 3.11. Dikdörtgen ve dairesel kanat dizileri (Incropera and Dewitt, 2006).<br />

Toplam yüzey verimi,<br />

q<br />

t<br />

t<br />

η o = =<br />

(3.7)<br />

q mak hA tθ<br />

b<br />

olarak tanımlanır. Burada At, kanatların ve üzerlerine yerleşikleri yüzeyin (asal<br />

yüzey) toplam alanı, Qt ise kanatlar ve asal yüzeyden olan toplam ısı geçişini<br />

göstermektedir.<br />

N adet kanat varsa ve asal yüzey alanı Ab ile gösterilirse, toplam yüzey alanı;<br />

t<br />

f<br />

65<br />

q<br />

A = NA + A<br />

(3.8)<br />

A t f 1<br />

b<br />

= NA + 2πr<br />

( H − Nt)<br />

(3.9)<br />

şeklindedir. Olabilecek en faza ısı geçişi, kanat yüzeylerinin ve asal yüzeyin Tb<br />

sıcaklığında olmaları durumunda gerçekleşir. Toplam yüzeyden taşınımla ısı geçişi;<br />

Q = Nη<br />

hA θ + hA θ<br />

(3.10)<br />

t<br />

f<br />

f<br />

b<br />

olarak gösterilebilir. Burada, ısı taşınım katsayısı h, kanat yüzeylerinde ve asal<br />

yüzeyde eşit kabul edilmiş olup, ηf bir kanadın verimidir. Buradan;<br />

b<br />

b


⎡ NA ⎤<br />

f<br />

q t = h[<br />

Nηf<br />

A f + ( A t − NA f ) ] θb<br />

= hA t ⎢1<br />

− ( 1−<br />

ηf<br />

) ⎥θ<br />

b (3.11)<br />

⎣ A t ⎦<br />

NA f<br />

η o = 1− ( 1−<br />

ηf<br />

)<br />

(3.12)<br />

A<br />

A<br />

p<br />

t<br />

2 2<br />

= 2π(<br />

2 − r )<br />

(3.13)<br />

2c<br />

olup kanatsız durumda taşınılma geçen ısı;<br />

wo<br />

1<br />

66<br />

1<br />

q = h(<br />

2πr<br />

H)<br />

θ<br />

(3.14)<br />

denklemi kullanılarak bulunabilir (Incropera and Dewitt, 2006).<br />

3.2.2. Deneyler<br />

Deneysel çalışmalarda öncelikle Reaktör 1 ve Reaktör 2 olmak üzere iki adet reaktör<br />

tasarlanmış ve imal edilmiştir. Ardından deney tesisatı oluşturulmuştur. Bu aşamadan<br />

sonra deneysel çalışmalara geçilmiştir. Deneysel çalışmalarda aşağıdaki prosedürler<br />

göz önüne alınmıştır;<br />

Deneylerde önce içerisine 1/10 oranında çelik bilye yerleştirilmiş 60 mm çapında ve<br />

80 mm yüksekliğindeki mekanik öğütücüde 4 saat süreyle öğütülen LaNi4.75Al0.25<br />

malzeme azot gazı altında <strong>hidrür</strong> reaktöre doldurulmuştur. Kullanılan azot gazı,<br />

alaşımda olası oksidasyonu önlemek için kullanılmıştır. Öğütücü hızı 500 d/d olarak<br />

ayarlanmıştır(Kaplan, 2009).<br />

4 saat boyunca öğütülen LaNi4.75Al0.25 malzeme azot gazı altında reaktöre<br />

doldurulmuştur. Her bir reaktöre 135 g ağırlığında, reaktör hacminin %70 oranında<br />

alaşım doldurulmuştur. Bu miktar her iki reaktörde de muhafaza edilmiştir.<br />

Reaktör düşük basınç altında (~10 -4 mmHg) vakum pompasına bağlanarak fırın<br />

içerisine yerleştirilmiştir. Yaklaşık 220 °C ye ısıtılan Reaktör 2 saat boyunca<br />

b


vakumlanmıştır. Bu işlemden sonra aktivasyon işlemini gerçekleştirebilmek için<br />

hidrojen tesisatına bağlanan reaktör literatür (Kaplan 2009) incelenerek 8 bar<br />

hidrojen basıncı altında şarj edilmiştir. Bu işlemler her iki reaktör için tam olarak<br />

aktivasyon işleminin gerçekleşmesi için 8 defa tekrar edilmiştir.<br />

Aktivasyon işlemi tam olarak gerçekleştikten sonra deneysel işlemlere geçilmiş,<br />

öncelikle Reaktör 2 ile deneylere başlanmıştır. Reaktör üzerine 4 adet termokupl z-<br />

ekseni boyunca 5, 20, 40 ve 60 mm konumlarına yerleştirilmiştir. Hızlı bağlantı<br />

parçasıyla tesisata bağlanan reaktör hassas terazi üzerine yerleştirilmiştir. Sıcaklık<br />

tarayıcı cihaz ve hassas terazi veri kaydedici kablolarla bilgisayara bağlanmıştır.<br />

Hidrojen şarj basıncı bir basınç regülâtörü yardımıyla 2, 4, 6 ve 8 bar gibi farklı<br />

basınçlara ayarlanmıştır. Deneylerde literatüre göre yaklaşık olarak 4000 s sürede<br />

reaktör sıcaklıkları oda şartlarına düştüğü için yapılan deneylerde de deney süresi<br />

4000 s olarak belirlenmiştir.<br />

Deneysel çalışmalarda şarj basıncı 8 bar olarak ayarlanmış ve 4000 s boyunca<br />

sıcaklık değişimi ve kütle değişimi kaydedilmiştir. Deney sonunda Reaktör 220 °C<br />

ye ısıtılarak depolanan hidrojen tamamen boşaltılmıştır. Daha sonra basınç 2, 4 ve 6<br />

bar olarak değiştirilerek deneylere devam edilmiştir.<br />

Reaktör 2’de bütün basınç değerlerinde deneyler gerçekleştirildikten sonra Reaktör 1<br />

deneylerine geçilmiştir. Reaktör 1’de öncelikle aktivasyon işlemi gerçekleştirişmiş<br />

ve ardından 2, 4, 6 ve 8 bar basınçlarda Reaktör 2’de olduğu gibi aynı şartlar altında<br />

aynı sürelerde deneyler tekrarlanmıştır.<br />

3.2.3. Yapay sinir ağları<br />

Yapay sinir ağları (YSA) kavramı, beynin çalışma ilkelerinin bilgisayarlar üzerinde<br />

benzetimlerinin yapılması fikri ile ortaya çıkmıştır. Yapılan ilk çalışmalar, beyni<br />

oluşturan nöronların matematik modellerinin kurulması şeklindedir. Günümüzde<br />

YSA olarak isimlendirilen alan, birçok nöronun belirli biçimlerde bir araya getirilip<br />

67


ir işlevin gerçeklenmesi üzerindeki yapısal olduğu kadar matematiksel sorulara<br />

yanıt arayan bir bilim dalı olmuştur. YSA, karar hızı açısından insan beyni ile<br />

yarışabilecek aşamayı henüz kat etmemiş olmalarına rağmen, karmaşık<br />

eşleştirmelerin hassas bir biçimde gerçekleyebilmeleri sayesinde uygulama alanlarını<br />

zamanla genişlemektedir. Bu bağlamda, bir YSA yapısının çözebileceği problem<br />

uzayı, insan beyninin çözebildiği problem uzayının oldukça kısıtlı bir kümesi<br />

olacaktır. YSA yaklaşımını uygulanabilir kılan bazı temel özelliklerin<br />

değerlendirilmesinde, bunun dikkate alınması gerekir.<br />

3.2.3.1. Yapay sinir ağlarının özellikleri<br />

Temel özelliklerden birincisi; sistemin paralelliği ve toplumsal işlevin yapısal olarak<br />

dağılmışlığıdır (Haykin, 1994). Bu ağ yapısı içerisinde birçok nöronun eşzamanlı<br />

olarak çalıştığı ve karmaşık bir işlevin çok sayıda nöron aktivitesinin bir araya<br />

gelmesinden oluştuğu anlamına gelir.<br />

Đkinci özellik ise genelleme yeteneğidir. Ağ yapısının, eğitim sırasında kullanılan<br />

sayısal bilgilerden eşleştirmeyi tanımlayan kaba özellikleri belirlemesi ve dolayısıyla<br />

eğitim sırasında kullanılmayan girdiler için de anlamlı yanıtlar üretilmesi, genelleme<br />

yeteneği olarak açıklanabilir.<br />

Bir başka özellik ise ağ fonksiyonunun doğrusal olmayan oluşudur. Yapı üzerinde<br />

dağılmış belli tipteki doğrusal olmayan alt birimler özellikle, istenen eşleştirmenin<br />

denetim ya da tanılama işlemlerinde olduğu gibi doğrusal olmayan olması<br />

durumunda işlevin doğru biçimde yerine getirilebilmesini matematik yönden olası<br />

kılar. Farklı nöron aktivasyon fonksiyonlarının irdelenmesi, işlevin doğru biçimde<br />

gerçeklenebilmesi için yapısal esnekliğe gerekliliği olan ağ parametrelerinin başarımı<br />

arttıracak biçimde güncellenebilmesi içindir.<br />

Son özellik ise; sayısal ortamda tasarlanan sinir ağı yaklaşımlarının tüm devre<br />

gerçeklenebilirliklerinin olmasıdır. Bu sistemler, robotik uygulamaları ile birlikte<br />

68


düşünüldüğünde, günlük hayatta yaşam kalitesinin arttırılmasında önemli rol<br />

oynayacaklardır.<br />

Günümüzde YSA’nın pek çok alanda uygulamasına rastlamak mümkündür.<br />

Özellikle; örüntü tanıma, işaret işleme, sistem tanılama ve doğrusal olmayan denetim<br />

alanlarında YSA değişik modelleri başarı ile kullanılmıştır (Narenda and Pathasaraty,<br />

1990; Efe, Kaynak, 1999; Efe et al. 1999).<br />

Ele alınan problemin YSA ile çözümünde karşılaşılan ilk seçenek öğrenme<br />

mekanizması üzerindedir. Literatürde öğreticili öğrenme, öğreticisiz öğrenme ve<br />

takviyeli öğrenme olmak üzere üç tip öğrenme stratejisinden bahsedilmektedir<br />

(Chen, 1996). Yaklaşımlar arasındaki temel farklılık; istenen çıkış değerlerinin var<br />

olup olmamasından kaynaklanmaktadır. Eğer eğitici, sistem çıkışlarının istenilen<br />

değerlerini sağlayıp YSA çıktıları ile doğrudan karşılaştırabiliyorsa, bu tip öğrenme,<br />

öğreticili öğrenme olarak adlandırılır. Bu genelde; veri seti üzerinden elde edilen<br />

istenilen çıkış değerleri ve YSA çıktıları arasında hesaplanan ortalama karesel hata<br />

gibi bir hata fonksiyonunun en aza indirgenmesi şeklinde formüle edilebilir.<br />

Takviyeli öğrenme; tam olarak istenilen çıkış değerlerinin bilinmiyor olduğu<br />

öğreticili öğrenmenin özel bir halidir. Yalnızca, YSA çıktısının doğru olup<br />

olmamasına ilişkin bilgiye göre belirlenir. Ağı yapılandırma koşulları, istenilen çıktı<br />

değerleri sağlayamıyorsa ve öğrenme, yalnızca girdi verilerine ilişkin korelasyonlarla<br />

gerçekleniyorsa, bu tip öğrenme öğreticisiz öğrenme olarak adlandırılır. Öğreticisiz<br />

öğrenme algoritmaları daha çok, sistemin geçmişte karşı karşıya kaldığı veri<br />

kümesinin içerdiği istatistiksel bilgilerin çıkarılmasını amaçlar. Böylelikle çok<br />

elemanlı veri kümeleri içerisinde, deneyim yoluyla bilgi genelleştirmesi yapılabilir.<br />

3.2.3.2. Yapay sinir ağlarının kullanım alanları<br />

YSA bilgi sınıflama ve bilgi yorumlamanın da içinde bulunduğu çok değişik<br />

problemlerin çözümünde ve de özellikle çok sayıda veriyi değerlendirme gerektiren<br />

işlemlerde oldukça yaygın olarak kullanılmaktadır. Bir problemin YSA ile<br />

69


çözülebilmesi için aşağıdaki şartlardan birinin sağlanması gerekmektedir (Öztemel,<br />

2003).<br />

• Problemlere, sadece yapay sinir ağları yönteminin pratik çözümler<br />

üretebilmesinin söz konusu olması,<br />

• Başka çözüm yolları olmasına rağmen, YSA daha kolay ve daha etkin<br />

çözümler üretebilecek olması.<br />

Uygulamalar incelendiğinde YSA, doğrusal olamayan, karmaşık, kesin olmayan, çok<br />

boyutlu, eksik ve problemin çözümü için matematiksel bir modelin olmadığı<br />

durumlarda yaygın olarak kullanıldığı görülmektedir. YSA mühendislikten tıbba,<br />

endüstriden iktisada kadar yaygın kullanımı olmakla birlikte bunlar sınıflandırma,<br />

tahmin ve modelleme olmak üzere üç grupta toplanabilir. En çok kullanıldığı<br />

alanların bazıları şu şekildedir:<br />

• Mekanik parçaların ömürlerini ve kırılmalarının tahmin edilmesi,<br />

• Veri madenciliği<br />

• Zeki araçlar ve robotlar için optimum rota belirleme,<br />

• Robot hareket mekanizmalarının kontrol edilmesi<br />

• Üretim planlama ve çizelgeleme, kalite kontrolü<br />

• Elektrik işareti tanıma<br />

• Radar ve sonar sinyalleri sınıflandırma<br />

• Optik karakter tanıma vs. (Öztemel, 2003; Elmas, 2003)<br />

3.2.3.3. Yapay sinir ağı yapıları<br />

YSA sınıflandırmak için birçok yaklaşım vardır. Örneğin, YSA öğrenme<br />

yöntemlerine veya eğitilme yöntemlerine göre; eğitimde kullanılan tekrarlanan veya<br />

tekrarlanmayan temel fonksiyonun aktivasyon fonksiyonu ile değişimine göre<br />

sınıflandırılabilmektedir (Ham et al., 2001). YSA hücrelerin birbirleri ile çeşitli<br />

şekillerde bağlanmalarından oluşmaktadır. Hücre çıkışları, ağırlıklar üzerinden diğer<br />

hücrelere ya da kendisine giriş olarak bağlanabilmekte ve bağlantılarda gecikme<br />

70


irimi de kullanılabilmektedir. Hücrelerin bağlantı şekillerine, öğrenme kurallarına<br />

ve aktivasyon fonksiyonlarına göre çeşitli YSA yapıları geliştirilmiştir.<br />

Đleri beslemeli YSA’da, hücreler katmanlar seklinde düzenlenir ve bir katmandaki<br />

hücrelerin çıkışları bir sonraki katmana ağırlıklar üzerinden giriş olarak<br />

verilmektedir. Giriş katmanı, dış ortamlardan aldığı bilgileri hiçbir değişikliğe<br />

uğratmadan orta (gizli) katmandaki hücrelere iletmektedir. Bilgi, orta ve çıkış<br />

katmanında işlenerek ağ çıkışı belirlenmektedir. Bu yapısı ile ileri beslemeli ağlar<br />

doğrusal olmayan statik bir işlevi gerçekleştirir. Đleri beslemeli 3 katmanlı YSA’nın,<br />

orta katmanında yeterli sayıda hücre olmak kaydıyla, herhangi bir sürekli fonksiyonu<br />

istenilen doğrulukta yaklaştırabileceği gösterilmiştir. En çok bilinen geriye yayılım<br />

öğrenme algoritması, bu tip YSA’nın eğitiminde etkin olarak kullanılmakta ve bazen<br />

bu ağlara geriye yayılım ağları da denmektedir. Şekil 3.12’de giriş, orta ve çıkış<br />

katmanı olmak üzere üç katmanlı ileri beslemeli YSA yapısı verilmektedir.<br />

3.2.3.4. Yapay sinir ağı ile modelleme<br />

Şekil 3.12. Đleri beslemeli YSA yapısı<br />

Depolama deneylerinden elde edilen sonuçları değerlendirmek ve depolanabilecek<br />

hidrojen miktarını önceden tahmin edebilmek amacıyla Pythia - Neural Network<br />

aracından yararlanılmıştır. Deneyler sonucunda kütle değerlerinin %80’i ağı eğitmek<br />

ve %20’si test amaçlı olarak kullanılmıştır.<br />

71


4. ARAŞTIRMA BULGULARI<br />

4.1. Deneysel Çalışmalardan Elde Edilen Bulgular<br />

Bu çalışmada Reaktör 1 ve Reaktör 2 olmak üzere iki farklı reaktör imal edilmiş ve<br />

her iki reaktöre aynı sürelerde öğütülmüş, aynı tane boyutlarında 135 g LaNi4.75Al0.25<br />

<strong>metal</strong> <strong>hidrür</strong> tozu doldurulmuştur. Yapılan deneylerde LaNi5 alaşımlarına Al<br />

ilavesinin <strong>hidrür</strong> reaksiyonlarına, depolama hızına, sıcaklık profiline, depolama<br />

kinetiklerine etkisi ve ayrıca <strong>hidrür</strong> reaktör tasarımının bu parametreler üzerindeki<br />

etkisi araştırılmıştır. Deneylerde hidrojen şarj basıncı 2, 4 6 ve 8 bar olarak<br />

değiştirilerek her basınç değeri için 4000 s boyunca işlem sürdürülmüş ve bu basınç<br />

değerlerinde reaktörlerin z-ekseni boyunca 5, 20, 40 ve 60 mm noktalarına bağlanan<br />

termokupllar vasıtasıyla zamana bağlı sıcaklık ve kütle değişimleri ölçülmüştür.<br />

Reaktörlerde yüksek şarj basınçlarında ölçülen sıcaklık değerleri, düşük şarj<br />

basınçlarında ölçülen sıcaklık değerlerinden daha yüksek değerlere çıkmıştır. Bunun<br />

sebebi şarj basıncındaki artışın reaksiyon hızını arttırıcı yönde etki etmesidir.<br />

Dolayısıyla yüksek şarj basınçlarında ekzotermik reaksiyon daha hızlı olmakta, şarj<br />

basıncı azaldıkça reaksiyon hızı da azalmaktadır. Reaksiyon sonucunda oluşan ısının<br />

merkezde daha yüksek olduğu, merkezden cidara doğru gelindikçe azaldığı<br />

görülmektedir. Ayrıca z-ekseni boyunca mesafe arttıkça sıcaklığın arttığı, maksimum<br />

sıcaklık eğrilerinin ise z=60 mm noktasında olduğu görülmektedir. Buna ilaveten<br />

reaktör uzunluğu boyunca farklı konumlar için maksimum sıcaklığa ulaşma süreleri<br />

de farklıdır. Reaktörün girişe yakın bölümü hidrojenle daha önce karşılaştığı için,<br />

reaktörün alt bölümüne göre burada maksimum sıcaklığa daha önce ulaşılmıştır.<br />

Giriş bölgesine yaklaştıkça hidrojen absorpsiyonu tamamlanırken, z-ekseni boyunca<br />

gidildikçe hidrojenle daha fazla kontak kurma şansı artmaktadır. Böylece reaktör<br />

tabanına doğru sıcaklık zamanla artarak maksimum değerlere ulaşmıştır. Ayrıca<br />

hidrojen absorpsiyon reaksiyonunun doğasında olan, ekzotermik reaksiyonla ısı<br />

üretiminden dolayı oluşan sıcaklık artışı sonucu reaktör girişine doğru reaksiyon<br />

sayısında düşmeler meydana gelmektedir.<br />

72


Bunların yanında reaktörlerde depolanan hidrojen miktarının zamana göre değişimini<br />

belirlemek için yapılan deneylerde ise hidrojen şarj basıncının ve reaktör tasarımının<br />

depolanan hidrojen miktarını etkilediği belirlenmiştir. Bu amaçla reaktörlere 4, 6 ve<br />

8 bar şarj basınçlarda hidrojen şarjı yapılmış ve depolanan hidrojen miktarı zamana<br />

bağlı olarak ölçülmüştür. Reaktörlerde hidrojen şarj basıncı arttıkça depolanan<br />

hidrojen miktarı da artmıştır. Ayrıca Reaktör 2’de ısı transferinin hızlı gerçekleşmesi<br />

sonucu depolanan hidrojen miktarı Reaktör 1’e göre daha fazla olmuştur. Hidrojen<br />

depolama hızı, ilk 500 s sürede hızlı ekzotermik reaksiyonlar sonucu çok hızlı<br />

artarken daha sonra reaksiyonların yavaşlamaya başlamasıyla düşmüştür.<br />

Diğer yandan LaNi5 alaşımlarına Al ilavesinin absorpsiyon hızı, aktivasyon<br />

özellikleri, kinetik özellikleri ve termodinamik özelliklerini geliştirdiği gözlenmiştir.<br />

LaNi5-xAlx alaşımlarında x değeri 0.25, 0.50, 0.75 ve 1 olabilmektedir. Ancak bu<br />

değerin çok fazla artması reaksiyon hızını yavaşlattığı bu nedenle çalışmada<br />

kullanılan x=0.25 değerinin optimum değer olduğu literatür çalışmalarından da<br />

anlaşılmaktadır.<br />

Reaktör 1 ve Reaktör 2’ye 2 bar basınçta 4000 s boyunca hidrojen şarj edilmiş ve<br />

deneylerde reaktörün z=5, 20, 40 ve 60 mm konumlarında zamana göre sıcaklık<br />

değişimleri Şekil 4.1a ve b’de verilmiştir. Buna göre reaktör merkezinden tabana<br />

doğru gidildikçe maksimum sıcaklığa ulaşma süresi artmaktadır. Bunun nedeni de<br />

ekzotermik reaksiyonun reaktör girişinde başlamasıdır. Reaktörün bütün<br />

konumlarında maksimum sıcaklığa yaklaşık 300 s süre içerisinde ulaşılmıştır. Bu<br />

sürede reaksiyon hızlı şekilde gerçekleşmekte ve sıcaklık maksimum değere<br />

ulaşmaktadır. Sıcaklık maksimum değere ulaştıktan sonra reaksiyonların azalmaya<br />

başlamasıyla, sıcaklık azalmaya ve 4000 s sürenin sonunda yaklaşık oda<br />

sıcaklıklarına düşmüştür. Reaktör 1’de maksimum sıcaklık z=60 mm konumunda<br />

209 s sürede 88.77 °C olarak ölçülmüştür. z=5 mm noktasında maksimum sıcaklık<br />

199 s sürede 83.59 °C, z=20 mm noktasında 206 s sürede 87.79 °C ve z=40 mm<br />

noktasında 224 s sürede 87.78 °C olarak ölçülmüştür. Reaktör 2’de ise maksimum<br />

sıcaklık z=60 mm konumunda 225 s sürede 81.33 °C olarak ölçülmüştür. Bunun<br />

yanında z=5 mm noktasında maksimum sıcaklık 219 s sürede 75.71 °C, z=20 mm<br />

73


noktasında 216 s sürede 79.7 °C ve z=40 mm noktasında 224 s sürede 80.7 °C olarak<br />

ölçülmüştür.<br />

Şekil 4.1. 2 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2<br />

74


Şekil 4.2. 4 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2<br />

Reaktör 1 ve Reaktör 2’ye 4 bar basınçta 4000 s boyunca hidrojen şarj edilmiş ve<br />

deneylerde reaktörlerin z=5, 20, 40 ve 60 mm konumlarında zamana göre sıcaklık<br />

75


değişimleri Şekil 4.2a ve b’de verilmiştir. Reaktör merkezinden tabana doğru<br />

gidildikçe maksimum sıcaklığa ulaşma süresi artmakta, reaktörün bütün<br />

konumlarında 300 s’de maksimum ulaşılmıştır. Bu sürede reaksiyon hızlı şekilde<br />

gerçekleşmekte ve sıcaklık maksimum değere ulaşmaktadır. Sıcaklık maksimum<br />

değere ulaştıktan sonra reaksiyonların azalmaya başlamasıyla sıcaklıkta azalmaya ve<br />

4000 s sonunda yaklaşık oda sıcaklıklarına düşmüştür. Reaktörlerde maksimum<br />

sıcaklıklar 2 bar şarj basınçtakine göre daha yüksek değerlere çıkmıştır. Reaktör 1’de<br />

maksimum sıcaklık z=60 mm noktasında 238 s sürede 102.92 °C olarak ölçülmüştür.<br />

z=5 mm noktasında maksimum sıcaklık 229 s sürede 94.94 °C, z=20 mm noktasında<br />

228 s sürede 100.66 °C ve z=40 mm noktasında 226 s sürede 103.25 °C olarak<br />

ölçülmüştür. Reaktör 2’de ise maksimum sıcaklık z=60 mm noktasında 258 s sürede<br />

98.1 °C olarak ölçülmüştür. Bunun yanında z=5 mm noktasında maksimum sıcaklık<br />

250 s sürede 84.49 °C, z=20 mm noktasında 251 s sürede 91.66 °C ve z=40 mm<br />

noktasında 256 s sürede 96.41 °C olarak ölçülmüştür.<br />

Reaktör 1 ve Reaktör 2’ye 6 bar basınçta 4000 s boyunca hidrojen şarj edilmiş ve<br />

deneylerde reaktörün z=5, 20, 40 ve 60 mm konumlarında zamana göre sıcaklık<br />

değişimleri Şekil 4.3a ve b’de verilmiştir. Buna göre reaktör merkezinden tabana<br />

doğru gidildikçe maksimum sıcaklığa ulaşma süresi artmaktadır. Reaktörün bütün<br />

konumlarında maksimum sıcaklığa yaklaşık 300 s süre içerisinde ulaşılmıştır. Bu<br />

sürede reaksiyon hızlı şekilde gerçekleşmekte ve sıcaklık maksimum değere<br />

ulaşmaktadır. Sıcaklık maksimum değere ulaştıktan sonra reaksiyonların azalmaya<br />

başlamasıyla sıcaklık azalmaya ve 4000 s’nin sonunda yaklaşık oda sıcaklıklarına<br />

düşmüştür. Reaktörlerde maksimum sıcaklıklar 4 bar basınca göre daha yüksek<br />

değerlere çıkarken; Reaktör 1’de z=60 mm konumunda 242 s sürede 116.24 °C<br />

olarak ölçülmüştür. z=5 mm noktasında maksimum sıcaklık 237 s sürede 95.8 °C,<br />

z=20 mm noktasında 245 sürede 107.12 °C ve z=40 mm noktasında 251 s sürede<br />

114.11 °C olarak ölçülmüştür. Reaktör 2’de ise maksimum sıcaklık z=60 mm<br />

noktasında 212 s sürede 113.27 °C olarak ölçülmüştür. Bunun yanında z=5 mm<br />

noktasında maksimum sıcaklık 210 s sürede 96.92 °C, z=20 mm noktasında 213 s<br />

sürede 107.45 °C ve z=40 mm noktasında 211 s sürede 111.08 °C olarak<br />

ölçülmüştür.<br />

76


Şekil 4.3. 6 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2<br />

77


Şekil 4.4. 8 bar basınçta reaktör uzunlukları boyunca zamana bağlı sıcaklık<br />

değişimleri, a) Reaktör 1 ve b) Reaktör 2<br />

Reaktör 1 ve Reaktör 2’ye 8 bar basınçta 4000 s boyunca hidrojen şarj edilmiş ve<br />

deneylerde reaktörün z=5, 20, 40 ve 60 mm konumlarında zamana göre sıcaklık<br />

değişimleri Şekil 4.4a ve b’de verilmiştir. Buna göre reaktör merkezinden tabana<br />

78


doğru gidildikçe maksimum sıcaklığa ulaşma süresi artmaktadır. Reaktörün bütün<br />

konumlarında maksimum sıcaklığa yaklaşık 300 s süre içerisinde ulaşılmıştır. Bu<br />

sürede reaksiyon hızlı şekilde gerçekleşmekte ve sıcaklık maksimum değere<br />

ulaşmaktadır. Sıcaklık maksimum değere ulaştıktan sonra reaksiyonların azalmaya<br />

başlamasıyla sıcaklık azalmaya ve 4000 s’nin sonunda yaklaşık oda sıcaklıklarına<br />

düşmüştür. Reaktörlerde maksimum sıcaklıklar 6 bar basınca göre daha yüksek<br />

değerlere çıkmıştır. Reaktör 1’de maksimum sıcaklık z=60 mm konumunda 256 s<br />

sürede 128.72 °C olarak ölçülmüştür. z=5 mm noktasında maksimum sıcaklık 246 s<br />

sürede 118.58 °C, z=20 mm noktasında 255 s sürede 126.13 °C ve 40 mm noktasında<br />

258 s sürede 128.12 °C olarak ölçülmüştür. Reaktör 2’de ise maksimum sıcaklık<br />

z=60 mm konumunda 244 s sürede 121.52 °C olarak ölçülmüştür. Diğer yandan z=5<br />

mm noktasında maksimum sıcaklık 248 s sürede 107.58 °C, z=20 mm noktasında<br />

258 s sürede 118.5 °C ve z=40 mm noktasında 252 s sürede 120.37 °C olarak<br />

ölçülmüştür.<br />

Sıcaklık (°C)<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

79<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.5. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z=5 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

2 bar


Reaktör 1’de z=5 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana göre sıcaklık değişimi Şekil 4.5’te verilmiştir. Bu konumda<br />

her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa ulaşılmış daha<br />

sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı olarak düşmeye<br />

başlamış 4000 sonunda oda sıcaklıklarına gelmiştir. Maksimum sıcaklık 8 bar<br />

basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta düşmüştür. Bunların<br />

yanında aynı şartlarda Reaktör 2’ye göre buradaki sıcaklıklar daha yüksek değerlere<br />

çıkmıştır. Ayrıca basınç yükseldikçe ekzotermik reaksiyon sayısının artmasına bağlı<br />

olarak hidrojen absorpsiyon hızı ve depolanan hidrojen miktarı artmıştır. Reaktör<br />

1’de z=5 mm konumunda; 2, 4, 6 ve 8 bar basınçlarda maksimum sıcaklıklara erişim<br />

süreleri 199, 229, 237, 246 s ve bu sürelerde gerçekleşme maksimum sıcaklıklar<br />

sırasıyla 83.59, 94.94, 95.8 ve 118.58 °C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

130<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

80<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.6. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z=20 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

Reaktör 1’de z=20 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.6’da verilmiştir. Bu<br />

konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

2 bar


ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

olarak düşmeye başlamış 4000 s sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

düşmüştür. Bunların yanında aynı şartlarda Reaktör 2’ye göre buradaki sıcaklıklar<br />

daha yüksek değerlere çıkmıştır. Ayrıca basınç yükseldikçe ekzotermik reaksiyon<br />

sayısının artmasıyla depolanan hidrojen miktarı da artmıştır. Reaktör 1’de z=20 mm<br />

konumunda; 2, 4, 6 ve 8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 206,<br />

228, 245, 255 s ve bu sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 87.79,<br />

100.66, 107.12 ve 126.13 °C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

130<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

81<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.7. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z=40 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

Reaktör 1’de z=40 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.7’de verilmiştir. Bu<br />

konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

olarak düşmeye başlamış 4000 s sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

2 bar


düşmüştür. Bunların yanında aynı şartlarda Reaktör 2’ye göre buradaki sıcaklıklar<br />

daha yüksek değerlere çıkmıştır. Ayrıca basınç yükseldikçe ekzotermik reaksiyon<br />

sayısının artmasıyla depolanan hidrojen miktarı da artmıştır. Reaktör 1’de z=40 mm<br />

konumunda; 2, 4, 6 ve 8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 224,<br />

226, 251, 258 s ve bu sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 87.78,<br />

103.25, 114.11 ve 128.12 °C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

130<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

82<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.8. Reaktör 1’de 2, 4, 6 ve 8 bar basınçlarda ve z=60 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

Reaktör 1’de reaktörde z=60 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6<br />

ve 8 bar 4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.8’de verilmiştir.<br />

Bu konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

olarak düşmeye başlamış 4000 s sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

düşmüştür. Bunların yanında aynı şartlarda Reaktör 2’ye göre buradaki sıcaklıklar<br />

daha yüksek değerlere çıkmıştır. Ayrıca basınç yükseldikçe ekzotermik reaksiyon<br />

sayısının artmasıyla depolanan hidrojen miktarı da artmıştır. Reaktör 1’de z=60 mm<br />

2 bar


konumunda; 2, 4, 6 ve 8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 209,<br />

238, 242, 256 s ve bu sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 88.77,<br />

102.92, 116.24 ve 128.72 °C olarak ölçülmüştür.<br />

Reaktör 2’de z=5 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.9’da verilmiştir. Bu<br />

konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

olarak düşmeye başlamış 4000 sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

düşmüştür. Ayrıca basınç yükseldikçe ekzotermik reaksiyon sayısının artmasıyla<br />

depolanan hidrojen miktarı da artmıştır. Reaktör 2’de z=5 mm konumunda; 2, 4, 6 ve<br />

8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 199, 250, 190, 188 s ve bu<br />

sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 75.71, 84.49, 96.92 ve 107.58<br />

°C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

83<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.9. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z=5 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

2 bar


Reaktör 2’de z=20 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.10’da verilmiştir. Bu<br />

konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

olarak düşmeye başlamış 4000 s sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

düşmüştür. Ayrıca basınç yükseldikçe ekzotermik reaksiyon sayısının artmasıyla<br />

depolanan hidrojen miktarı da artmıştır. Reaktör 2’de z=20 mm konumunda; 2, 4, 6<br />

ve 8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 216, 251, 213, 258 s ve<br />

bu sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 79.7, 91.66, 107.45 ve<br />

118.5 °C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

84<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.10. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z=20 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

Reaktör 2’de z=40 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.11’de verilmiştir. Bu<br />

konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

2 bar


olarak düşmeye başlamış 4000 s sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

düşmüştür. Ayrıca basınç yükseldikçe ekzotermik reaksiyon sayısının artmasıyla<br />

depolanan hidrojen miktarı da artmıştır. Reaktör 2’de z=40 mm konumunda; 2, 4, 6<br />

ve 8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 224, 256, 211, 244 s ve<br />

bu sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 80.7, 96.41, 111.08 ve<br />

120.37 °C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

85<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.11. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z=40 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

Reaktör 2’de z=60 mm konumunda farklı hidrojen şarj basınçlarda 2, 4, 6 ve 8 bar<br />

4000 s boyunca zamana bağlı sıcaklığın değişimi Şekil 4.12’de verilmiştir. Bu<br />

konumda her bir basınç değerinde yaklaşık 300 s sürede maksimum sıcaklığa<br />

ulaşılmış daha sonra reaksiyon sayısı azalmaya ve sıcaklıkta reaksiyonlara bağlı<br />

olarak düşmeye başlamış 4000 s sonunda oda sıcaklıklarına gelmiştir. Maksimum<br />

sıcaklık 8 bar basınçta gerçekleşmiş ve basınç düştükçe maksimum sıcaklıkta<br />

düşmüştür. Ayrıca basınç yükseldikçe ekzotermik reaksiyon sayısının artmasıyla<br />

depolanan hidrojen miktarı da artmıştır. Reaktör 2’de z=60 mm konumunda; 2, 4, 6<br />

2 bar


ve 8 bar basınçlarda maksimum sıcaklıklara erişim süreleri 225, 258, 212, 244 s ve<br />

bu sürelerde gerçekleşen maksimum sıcaklıklar sırasıyla 81.33, 98.1, 113.27 ve<br />

121.52 °C olarak ölçülmüştür.<br />

Sıcaklık (°C)<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

0 800 1600 2400 3200 4000<br />

Zaman (s)<br />

86<br />

8 bar<br />

6 bar<br />

4 bar<br />

Şekil 4.12. Reaktör 2’de 2, 4, 6 ve 8 bar basınçlarda ve z=60 mm noktasında zamana<br />

bağlı sıcaklık değişimleri<br />

Hidrojen şarj basıncının absorpsiyon hızına ve depolanan hidrojen kütlesine etkisinin<br />

belirlenmesi için yapılan deneylerde hidrojen şarj basıncı 4, 6 ve 8 bar basınçlara<br />

ayarlanarak kütle ölçümleri yapılmıştır. Deneylerde şarj basıncı arttıkça reaksiyon<br />

sayısının arttığı ve buna bağlı olarak depolanan hidrojen miktarının arttığı, ayrıca<br />

Reaktör 2’de oluşan reaksiyon ısısının kanatçıklar yardımıyla daha hızlı<br />

uzaklaştırılmasına bağlı olarak depolanan hidrojen miktarının Reaktör 1’e göre daha<br />

fazla olduğu görülmüştür.<br />

Reaktör 1 ve Reaktör 2 için 4 bar hidrojen şarj basıncında yapılan deneylerde zamana<br />

bağlı kütle değişimleri Şekil 4.13’te verilmiştir. Deneylerde depolanan hidrojen<br />

kütlesinin 2400 s süreden sonra değişmediği gözlendiği için deney süreleri 2400 s<br />

olarak belirlenmiştir. Đlk 300 s içerisinde reaktörlerde sıcaklık hızlı şekilde<br />

2 bar


yükselirken depolanan hidrojen kütlesi de hızlı şekilde artmıştır. 1000 s süreden<br />

sonra reaksiyonların azalmasıyla sıcaklık düşmeye başlamış ve hidrojen absorpsiyon<br />

hızı yavaşlamıştır. 4 bar basınçta 135 g LaNi4.75Al0.25 alaşım içerisinde Reaktör 1’de<br />

depolanan hidrojen miktarı 1.01 g olurken, Reaktör 2’de 1.13 g hidrojen<br />

depolanmıştır.<br />

Depolanan hidrojen kütlesi (g)<br />

1,2<br />

1,0<br />

0,8<br />

0,6<br />

0,4<br />

0,2<br />

0,0<br />

0 500 1000 1500 2000 2500<br />

Zaman (s)<br />

87<br />

Reaktör 2<br />

Reaktör 1<br />

Şekil 4.13. Reaktör 1 ve Reaktör 2’de 4 bar basınçta depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi<br />

Reaktör 1 ve Reaktör 2 için 6 bar şarj basıncında yapılan deneylerde zamana bağlı<br />

kütle değişimleri Şekil 4.14’te verilmiştir. Deneylerde depolanan hidrojen kütlesinin<br />

2400 s süreden sonra değişmediği gözlendiği için deney süreleri 2400 s olarak<br />

belirlenmiştir. Đlk 300 s içerisinde reaktörlerde sıcaklık hızlı şekilde yükselirken<br />

depolanan hidrojen kütlesi de hızlı şekilde artmıştır. 1000 s süreden sonra<br />

reaksiyonların azalmasıyla sıcaklık düşmeye başlamış ve hidrojen absorpsiyon hızı<br />

yavaşlamıştır. 6 bar basınçta 135 g LaNi4.75Al0.25 alaşım içerisinde Reaktör 1’de<br />

depolanan hidrojen miktarı 1.17 g olurken, Reaktör 2’de 1.25 g hidrojen<br />

depolanmıştır.


Depolanan hidrojen kütlesi (g)<br />

1,4<br />

1,2<br />

1,0<br />

0,8<br />

0,6<br />

0,4<br />

0,2<br />

0,0<br />

0 500 1000 1500 2000 2500<br />

Zaman (s)<br />

88<br />

Reaktör 2<br />

Reaktör 1<br />

Şekil 4.14. Reaktör 1 ve Reaktör 2’de 6 bar basınçta depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi<br />

Depolanan hidrojen kütlesi (g)<br />

1,6<br />

1,4<br />

1,2<br />

1,0<br />

0,8<br />

0,6<br />

0,4<br />

0,2<br />

Reaktör 2<br />

Reaktör 1<br />

0,0<br />

0 500 1000 1500 2000 2500<br />

Zaman (s)<br />

Şekil 4.15. Reaktör 1 ve Reaktör 2’de 8 bar basınçta depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi


Reaktör 1 ve Reaktör 2 için 8 bar şarj basıncında yapılan deneylerde zamana bağlı<br />

kütle değişimleri Şekil 4.15’te verilmiştir. Deneylerde depolanan hidrojen kütlesinin<br />

2400 s süreden sonra değişmediği gözlendiği için deney süreleri 2400 s olarak<br />

belirlenmiştir. Đlk 300 s içerisinde reaktörlerde sıcaklık hızlı şekilde yükselirken<br />

depolanan hidrojen kütlesi de hızlı şekilde artmıştır. 1000 s süreden sonra<br />

reaksiyonların azalmasıyla sıcaklık düşmeye başlamış ve hidrojen absorpsiyon hızı<br />

yavaşlamıştır. 8 bar basınçta 135 g LaNi4.75Al0.25 alaşım içerisinde Reaktör 1’de<br />

depolanan hidrojen miktarı 1.39 g olurken, Reaktör 2’de 1.55 g hidrojen<br />

depolanmıştır.<br />

Depolanan hidrojen kütlesi (g)<br />

1,4<br />

1,2<br />

1,0<br />

0,8<br />

0,6<br />

0,4<br />

0,2<br />

0,0<br />

0 500 1000 1500 2000 2500<br />

Zaman (s)<br />

89<br />

8 bar<br />

6 bar<br />

Şekil 4.16. Reaktör 1’de 4, 6 ve 8 bar basınçlarda depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi<br />

Reaktör 1’de farklı hidrojen şarj basınçlarında 4, 6 ve 8 bar depolanan hidrojen<br />

kütlelerinin zamana göre değişimleri karşılaştırılmalı olarak Şekil 4.16’da verilmiştir.<br />

Reaktör 1’de hidrojen şarj basıncı arttıkça ekzotermik reaksiyon artmakta dolayısıyla<br />

depolanan hidrojen miktarı ve şarj hızı artmaktadır. Đlk 500 s sürede reaksiyon hızlı<br />

bir şekilde gerçekleşmekte buna paralel olarak hidrojen kütlesi de artmaktadır. Bu<br />

süreden sonra kütle artış hızı yavaşlamakta ve maksimum kapasiteye geldikten sonra<br />

4 bar


sabit bir değerde kalmaktadır. Reaktörde sıcaklık artışıyla kütle artışı bir birine<br />

paralel olarak gerçekleşmektedir. Reaktör 1’de şarj basıncı 4, 6 ve 8 bar olarak<br />

değiştirilmiştir. Bu basınçlarda depolanan hidrojen miktarı sırasıyla 1.01, 1.17 ve<br />

1.39 g olarak gerçekleşmiştir. Hidrojen şarj basıncıyla orantılı olarak ekzotermik<br />

reaksiyon hızlanmakta ve dolayısıyla depolanan hidrojen miktarı da artmıştır.<br />

Reaktör 2’de göre ısı transferinin yavaş gerçekleşmesinden dolayı reaksiyon hızı<br />

daha yavaştır. Buna rağmen yaklaşık 500 s sürede reaktör tam olarak şarj edilmiştir.<br />

Ancak depolanan hidrojen miktarları Reaktör 2’ye göre daha az gerçekleşmiştir.<br />

Ayrıca LaNi5 alaşımlarına Al ilavesi reaksiyon hızını arttırmıştır.<br />

Depolanan hidrojen kütlesi (g)<br />

1,6<br />

1,4<br />

1,2<br />

1,0<br />

0,8<br />

0,6<br />

0,4<br />

0,2<br />

0,0<br />

0 500 1000 1500 2000 2500<br />

Zaman (s)<br />

90<br />

6 bar<br />

Şekil 4.17. Reaktör 2’de 4, 6 ve 8 bar basınçlarda depolanan hidrojen kütlesinin<br />

zamana bağlı değişimi<br />

Reaktör 2’de farklı hidrojen şarj basınçlarında 4, 6 ve 8 bar depolanan hidrojen<br />

kütlelerinin zamana bağlı değişimleri karşılaştırılmalı olarak Şekil 4.17’de<br />

verilmiştir. Reaktörde hidrojen şarj basıncı arttıkça ekzotermik reaksiyon artmakta,<br />

dolayısıyla depolanan hidrojen miktarı ve şarj hızı artmaktadır. Đlk 500 s sürede<br />

reaksiyon hızlı bir şekilde gerçekleşmekte buna paralel olarak hidrojen kütlesi de<br />

artmaktadır. Bu süreden sonra kütle artış hızı yavaşlamakta ve maksimum kapasiteye<br />

geldikten sonra sabit bir değerde kalmaktadır. Reaktörde sıcaklık artışıyla kütle artışı<br />

8 bar<br />

4 bar


ir birine paralel olarak gerçekleşmektedir. Reaktörde şarj basıncı 4, 6 ve 8 bar<br />

olarak değiştirilmiştir. Bu basınçlarda depolanan hidrojen miktarı sırasıyla 1.13, 1.25<br />

ve 1.55 g olarak gerçekleşmiştir. Hidrojen şarj basıncıyla orantılı olarak ekzotermik<br />

reaksiyon hızlanmış ve dolayısıyla depolanan hidrojen miktarı da artmıştır. Reaktör<br />

1’e göre ısının daha hızlı uzaklaştırılmasıyla yaklaşık 500 s sürede reaktör tam olarak<br />

şarj edilmiştir. Ayrıca LaNi5 alaşımlarına Al ilave edilmesi reaksiyon hızını<br />

arttırmıştır.<br />

4.2. YSA Kullanılarak Yapılan Modellemeden Elde Edilen Bulgular<br />

Metal <strong>hidrür</strong> reaktörlerde depolanan hidrojen miktarını tahmin edebilmek amacıyla<br />

iki farklı ağ oluşturulmuştur. YSA bazı deney sonuçlarıyla eğitilerek lineer olmayan<br />

problemlerin çözümünde kullanılabilmektedir. YSA’nda veriler iki kısma ayrılır.<br />

Birinci kısım ağın eğitiminde kullanılırken ikinci kısım modelin geçerliliğini<br />

göstermek için kullanılır. Ağın eğitilmesiyle bağlantı ağırlıklarını belirlenir, buradan<br />

çıkış değerlerini tahmin ederek hedef değerlere yaklaşmayı amaçlamaktadır. Ağın<br />

verimi, istenen ve tahmin edilen değerler arasındaki % hata hesaplanarak<br />

belirlenebilir (Satish and Setty, 2005).<br />

Geliştirilen YSA modelinin verimini belirlemek için farklı hata analiz yöntemleri<br />

kullanılmaktadır. Genellikle, determinasyon katsayısı (R 2 ), karelerin karekökü<br />

ortalama hatası (RMSE) ve ortalama mutlak hata yüzdesi (MAPE) yöntemleridir.<br />

R 2 ’nin en yüksek değeri, RMSE ve MAPE’nin en düşük değerleri, geliştirilen YSA<br />

modelinin en iyi verimini göstermektedir. Bu parametreler aşağıdaki denklemlerden<br />

hesaplanmaktadır;<br />

RMSE<br />

⎡<br />

2 ⎤<br />

2 ∑(<br />

mden<br />

, i − m YSA,<br />

i )<br />

R = 1−<br />

⎢<br />

⎥<br />

(4.1)<br />

2<br />

⎢⎣<br />

∑(<br />

m YSA ) ⎥⎦<br />

⎡1<br />

⎢<br />

⎣n<br />

n<br />

= ∑ i=<br />

1<br />

( m<br />

YSA,<br />

i<br />

91<br />

− m<br />

den,<br />

i<br />

)<br />

2<br />

⎤<br />

⎥<br />

⎦<br />

1/<br />

2<br />

(4.2)


m YSA − mden<br />

MAPE = × 100<br />

(4.3)<br />

m<br />

YSA<br />

YSA yönteminin kullanılabilmesi için bütün giriş ve çıkış değerleri 0 ile 1 arasında<br />

normalize edilmiştir;<br />

V<br />

N<br />

VR<br />

− Vmin<br />

= (4.4)<br />

V − V<br />

mak<br />

Burada VR, Vmin ve Vmak sırasıyla verilerin gerçek, minimum ve maksimum<br />

değerlerini göstermektedir.<br />

Kullanılan modelde her iki ağ için giriş verileri; deney süresi, şarj basıncı ve sıcaklık<br />

değişimi, çıkış değeri ise; depolanan hidrojen kütlesidir. Çıkış değeri olan kütle<br />

değişimini belirlemede, genelleme ve enterpolasyon yapabilmek için çok katmanlı<br />

ileri beslemeli ağ eğitilmiş ve test edilmiştir. Kullanılan ağ yapıları Şekil 4.18 ve<br />

4.19’da verilmiştir. Gizli tabaka Reaktör 1’de 12, Reaktör 2’de 13 nörona sahiptir.<br />

Giriş tabakası zaman, sıcaklık ve basınç olmak üzere 3 giriş nöronuna sahiptir.<br />

Çıkışta ise kütle değişimini veren bir nöron vardır. Yapılan modellemede 2400 giriş<br />

verisi kullanılmıştır. Alınan verilerin tamamı YSA’nı eğitmek için kullanılmamış,<br />

verilerin bazıları test amacıyla kullanılmıştır.<br />

92<br />

min<br />

Şekil 4.18. Reaktör 1 için oluşturulan YSA yapısı


Şekil 4.19. Reaktör 2 için oluşturulan YSA yapısı<br />

Yapılan modellemede geri yayılımlı öğrenme algoritmalı tek gizli katmanla çalışan<br />

12 (Reaktör 1) ve 13 (Reaktör 2) nöronlu ileri beslemeli ağ kullanılmıştır. Algoritma<br />

değişimi için Levenberg-Markuardt yöntemi kullanılmıştır. LM yöntemiyle gerçekte<br />

Newton’un yaklaşma yöntemi kullanmaktadır. Kullanılan FERMI fonksiyonu;<br />

F<br />

z<br />

1<br />

= (4.5)<br />

−4<br />

( z−0.<br />

5)<br />

1+<br />

e<br />

Burada z giriş ağırlıkları toplamıdır. Depolanan hidrojen kütlesi m, LM algoritması<br />

kullanılarak aşağıdaki denklemden bulunabilir;<br />

Reaktör 1 için;<br />

Reaktör 2 için;<br />

1<br />

F ( m)<br />

= −4<br />

( 0.<br />

804230 F7<br />

−0.<br />

484840 F8<br />

−1.<br />

223922 F9<br />

−2.<br />

513080 F10<br />

+ 1.<br />

883421F11<br />

) −0.<br />

5<br />

(4.6)<br />

1+<br />

e<br />

1<br />

F ( m)<br />

= −4<br />

( 0.<br />

374875 F7<br />

−1.<br />

40497 F8<br />

+ 0.<br />

910134 F9<br />

+ 0.<br />

313569 F10<br />

−1.<br />

471658 F11<br />

+ 2.<br />

171917 F12<br />

) −0.<br />

5 (4.7)<br />

1+<br />

e<br />

Her bir giriş değeri bağlantı ağırlıklarıyla çarpılır. Ağırlıklandırılmış giriş değerleri<br />

doğrusal değerlere eklenerek çıkış değerlerine dönüştürülür. Bu çıkış değerleri ikinci<br />

93


katmandaki nöronlar için giriş değerleri olarak kullanılır. Buradan gizli katman<br />

seviyesindeki her bir Ei hesaplanmıştır.<br />

1<br />

Fi = (4.8)<br />

−4<br />

( E I −0.<br />

5)<br />

1+<br />

e<br />

Nöron 1-8 için Ei denklem (4.9) kullanılarak hesaplanmıştır. Burada i nöron sayısını<br />

göstermektedir. Elde edilen sabitler Çizelge 4.1 ve 4.2’de verilmiştir.<br />

= C t + C P + C T<br />

(4.9)<br />

Ei 1i<br />

2i<br />

3i<br />

Çizelge 4.1. Reaktör 1 YSA modeli nöron 1–6 için LM algoritması denklem sabitleri<br />

i<br />

C1i<br />

Sabitler<br />

C2i C3i<br />

1 2.206519 -2.093989 1.045560<br />

2 1.020514 -0.901983 0.111571<br />

3 0.670883 0.230205 -1.553855<br />

4 1.825874 0.059120 0.282857<br />

5 -1.288623 -0.341080 0.045053<br />

6 -1.483582 -0.360242 -1.554945<br />

Çizelge 4.2. Reaktör 2 YSA modeli nöron 1–6 için LM algoritması denklem sabitleri<br />

i<br />

C1i<br />

Sabitler<br />

C2i C3i<br />

1 1.098551 -0.335540 0.582720<br />

2 0.199532 -0.584620 0.025831<br />

3 -0.152180 0.738283 -0.267530<br />

4 2.026555 -1.797250 1.008123<br />

5 -1.773360 0.278109 0.113874<br />

6 -1.431860 0.053772 -1.217070<br />

94


Çizelge 4.3. Reaktör 1 YSA modeli nöron 7–12 için LM algoritması denklem<br />

sabitleri<br />

i<br />

C1i C2i<br />

Sabitler<br />

C3i C4i C5i C6i<br />

7 0.364408 0.787332 0.088055 0.085690 0.524498 0.394243<br />

8 0.654251 -0.567383 0.504829 -0.703947 -0.266245 -0.846468<br />

9 -0.768490 0.041987 1.057702 -1.173076 -0.648012 -0.478241<br />

10 -2.171420 -0.234487 1.413316 -0.995119 1.213655 2.626653<br />

11 -0.189810 -1.778331 1.266436 0.637397 -1.268316 -0.839176<br />

12 0.804230 -0.484840 -1.223922 -2.513080 1.883421<br />

Çizelge 4.4. Reaktör 2 YSA modeli nöron 7–13 için LM algoritması denklem<br />

sabitleri<br />

i<br />

C1i C2i<br />

Sabitler<br />

C3i C4i C5i C6i<br />

7 -0.980440 -0.167450 -0.08593 -0.977760 0.841411 -0.864720<br />

8 0.349570 0.705810 1.171809 -2.350020 1.228656 1.152940<br />

9 1.091870 -1.548970 -0.212430 0.219092 -1.029880 0.467993<br />

10 -0.862360 -0.871070 -0.339940 -0.915530 0.271038 0.130278<br />

11 -0.606370 0.039687 0.669292 -0.743100 0.070864 1.167493<br />

12 0.045887 -1.283120 1.860205 -0.715690 0.563241 -1.475930<br />

13 0.374875 -1.404970 0.910134 0.313569 -1.471660 2.171917<br />

Reaktör 1 ve Reaktör 2 için sırasıyla nöron 7–11, nöron 7–12 denklem (4.10)’da<br />

yerine yazılarak Ei hesaplanmıştır. Modelden elde edilen sabitler Çizelge 4.3 ve<br />

4.4’de verilmiştir.<br />

E = C F + C F + C F + C F + C F + C F<br />

(4.10)<br />

i<br />

1i<br />

1<br />

2i<br />

2<br />

3i<br />

3<br />

Reaktör 1’de nöron 12 ve Reaktör 2’de nöron 13 denklem (4.11) ve (4.12)<br />

kullanılarak Ei hesaplanmıştır.<br />

Reaktör 1 için;<br />

i<br />

7i<br />

7<br />

8i<br />

8<br />

9i<br />

9<br />

95<br />

4i<br />

4<br />

E = C F + C F + C F + C F + C F<br />

(4.11)<br />

10i<br />

10<br />

5i<br />

5<br />

11i<br />

6i<br />

11<br />

6


Reaktör 2 için;<br />

E = C F + C F + C F + C F + C F + C F<br />

(4.12)<br />

i<br />

7i<br />

7<br />

8i<br />

8<br />

9i<br />

9<br />

10i<br />

Reaktör 1 ve Reaktör 2’de depolanan hidrojen miktarını YSA kullanarak tahmin<br />

edebilmek için geliştirilen model eğitim verilerinin deneysel verilerle<br />

karşılaştırılması sırasıyla Şekil 4.20 ve 4.21’de verilmiştir.<br />

96<br />

10<br />

11i<br />

Çizelge 4.5. Normalizasyon değerleri<br />

11<br />

12i<br />

Reaktör 1 Reaktör 2<br />

Vmin Vmak Vmin Vmak<br />

t 0 2400 0 2400<br />

P 4 8 4 8<br />

T 36.01 128.72 37.30 120.38<br />

Çizelge 4.6. Đstatiksel veriler<br />

R 2 RMSE MAPE<br />

Reaktör 1 0.999 0.00006 0.0214<br />

Reaktör 2 0.998 0.01258 0.0042<br />

YSA eğitim verileri ile deneyler sonucu elde edilen verilerin her ikisi arasındaki<br />

sapmalar çok küçük ve kabul edilebilir seviyede olduğu görülmektedir. Buna göre<br />

depolanan hidrojen kütlesinin oluşturulan ağ yardımıyla Reaktör 1’de %99.9 ve<br />

Reaktör 2’de %99.8 doğrulukta tahmin edilebileceği görülmüştür.<br />

12


YSA Verileri<br />

1.6<br />

1.4<br />

1.2<br />

1.0<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0.0<br />

97<br />

R 2 = 0.9991<br />

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6<br />

Deneysel veriler<br />

Şekil 4.20. Reaktör 1 için YSA eğitim verileri ile deneysel verilerin karşılaştırılması<br />

YSA Verileri<br />

1.6<br />

1.4<br />

1.2<br />

1.0<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0.0<br />

R 2 = 0.9984<br />

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6<br />

Deneysel veriler<br />

Şekil 4.21. Reaktör 2 için YSA eğitim verileri ile deneysel verilerin karşılaştırılması<br />

Reaktör 1 ve Reaktör 2 için oluşturulan YSA modeli ile deneyler sonucunda elde<br />

edilen kütle değişim verileri karşılaştırmalı olarak sırasıyla Şekil 4.22 ve 4.23’te<br />

verilmiştir.


Şekil 4.22. Reaktör 1’de hidrojen kütlesinin zamana ve basınca bağlı değişimi<br />

Şekil 4.23. Reaktör 2’de hidrojen kütlesinin zamana ve basınca bağlı değişimi<br />

98


Hidrojen şarjında ilk 500 s süre içerisinde reaksiyonlara bağlı olarak depolanan<br />

hidrojen kütlesi hızlı şekilde artmıştır. Buna bağlı olarak YSA verilerinin hata oranı<br />

artarken, reaksiyonların yavaşlamaya başlamasıyla kütle değişimi de yavaşlamış ve<br />

YSA verileri deneysel sonuçlara yaklaşmıştır. Şekillerde de görüldüğü gibi sıcaklık,<br />

basınç ve zamana bağlı olarak geliştirilen YSA modeli ile elde edilen kütle<br />

verileriyle deneysel olarak elde edilen kütle verilerinin uyum içinde olduğu<br />

görülmüştür.<br />

99


5. TARTIŞMA VE SONUÇ<br />

5.1. Sonuçlar<br />

Bu çalışmada, aynı iç hacmine sahip iki farklı reaktör tasarlanmış ve imal edilmiştir.<br />

Reaktör 1 temel esaslarda, Reaktör 2 kanatçıklı olarak imal edilerek kanatçıkların<br />

reaktör içerisindeki ısı transferine ve reaksiyon kinetiklerine etkisi incelenmiştir.<br />

Ayrıca LaNi5 alaşımlarına Al ilavesinin reaksiyon hızına, aktivasyon özelliklerine ve<br />

reaksiyon kinetiklerine etkisi belirlenmiştir. Reaktörlerde sıcaklık değişim<br />

profillerini görmek amacıyla 2–8 bar arasında farklı hidrojen şarj basınçları altında,<br />

reaktörlerin farklı konumlarından sıcaklık ölçümleri yapılmıştır. Bunun yanında<br />

hidrojen şarj basıncının ve reaktör tasarımının depolanan hidrojen miktarına etkisini<br />

belirlemek için her iki reaktörle 4, 6 ve 8 bar basınçlarda kütle ölçümleri yapılmıştır.<br />

Bunların yanında, depolanan hidrojen miktarını teorik olarak tahmin edebilmek<br />

amacıyla Pythia-Neural Networks programı kullanılarak bir YSA modeli<br />

geliştirilmiştir. Yapılan bu çalışmalardan şu sonuçlar elde edilmiştir;<br />

LaNi4.75Al0.25 alaşımlarında aktivasyon hızı, depolama hızı ve miktarını arttırmak<br />

için 4 saat boyunca öğütme işlemi yapılmıştır. Böylece alaşım tane boyutları<br />

küçültülerek yüzey alanının, dolayısıyla depolanan hidrojen miktarının arttırılması<br />

sağlanmıştır.<br />

LaNi4.75Al0.25 alaşımları kullanılması LaNi5 alaşımlarına göre reaksiyon hızında<br />

önemli artışlar sağlanmıştır. Ayrıca reaksiyon kinetiklerinde ve termodinamik<br />

özeliklerinde önemli gelişmeler sağlamıştır. Bu nedenle reaktörlerin her bir basınç<br />

değerinde 400 s sürede şarj edilmesi sağlanmıştır.<br />

Reaktörlere hidrojen şarjı esnasında ilk 300 s sürede, bütün basınç değerlerinde<br />

ekzotermik reaksiyon sonucu sıcaklık hızlı bir şekilde yükselmiştir. Reaktörlerde<br />

oluşan maksimum sıcaklık değerleri hidrojen şarj basıncına ve kanatçıksız-kanatçıklı<br />

olması durumlarına bağlı olarak değişmekle beraber basınç arttıkça maksimum<br />

sıcaklık değerlerinin de arttığı görülmüştür. Ayrıca Reaktör 2’nin kanatçıklı<br />

100


olmasından dolayı ısı hızlı şekilde reaktörden uzaklaştırılmış ve bu nedenle sıcaklık<br />

değerleri Reaktör 1’e göre daha düşük değerlere ulaşmıştır.<br />

Reaktörlere hidrojen şarjında, <strong>metal</strong> <strong>hidrür</strong> alaşımlarıyla hidrojen öncelikle reaktör<br />

girişlerine yakın bölgelerle temas ettiğinden dolayı, ilk önce bu bölgede sıcaklık<br />

artışı maksimum değerlere ulaşmıştır. Reaktör z-ekseni boyunca gidildikçe<br />

maksimum sıcaklıklara ulaşma süreleri de artmaktadır. Reaktörün alt kısımlarında<br />

reaksiyon devam ederken reaktör girişinde reaksiyonlar azalmaya başlamıştır.<br />

Sıcaklık değerleri yaklaşık olarak 300 s sürede bütün noktalarda maksimum değere<br />

ulaştıktan sonra 4000 s sonunda oda sıcaklıklarına düşmektedir.<br />

Reaktör 1’de 2 bar şarj basıncında maksimum sıcaklıklar z=5 mm noktasında 199 s<br />

sürede 75.71 °C, z=20 mm noktasında 206 s sürede 79.7 °C, z=40 mm noktasında<br />

224 s sürede 80.7 °C ve z=60 mm noktasında 209 s sürede 88.77 °C olarak<br />

ölçülmüştür. 4 bar şarj basıncında maksimum sıcaklıklar z=5 mm noktasında 329 s<br />

sürede 94.94 °C, z=20 mm noktasında 368 s sürede 100.66 °C, z=40 mm noktasında<br />

396 s sürede 103.25 °C ve z=60 mm noktasında 398 s sürede 102.92 °C olarak<br />

ölçülmüştür. 6 bar şarj basıncında maksimum sıcaklıklar z=5 mm noktasında 237 s<br />

sürede 95.8 °C, z=20 mm noktasında 245 s sürede 107.12 °C, z=40 mm noktasında<br />

251 s sürede 114.11 °C ve z=60 mm konumunda 292 s sürede 116.24 °C olarak<br />

ölçülmüştür. 8 bar şarj basıncında maksimum sıcaklıklar z=5 mm noktasında 276 s<br />

sürede 118.58 °C, z=20 mm noktasında 285 s sürede 126.13 °C, z=40 mm<br />

noktasında 306 s sürede 128.12 °C ve z=60 mm konumunda 356 s sürede 128.72 °C<br />

olarak ölçülmüştür.<br />

Reaktör 2’de 2 bar şarj basıncında maksimum sıcaklıklar 5 mm noktasında 199 s<br />

sürede 75.71 °C, 20 mm noktasında 206 s sürede 79.7 °C, 40 mm noktasında 224 s<br />

sürede 80.7 °C ve 60 mm noktasında 225 s sürede 81.33 °C olarak ölçülmüştür. 4<br />

bar şarj basıncında maksimum sıcaklıklar z=5 mm noktasında 250 s sürede 84.49<br />

°C, z=20 mm noktasında 251 s sürede 91.66 °C, 40 mm noktasında 256 s sürede<br />

96.41 °C ve z=60 mm noktasında 268 s sürede 98.1 °C olarak ölçülmüştür. 6 bar şarj<br />

basıncında maksimum sıcaklıklar z=5 mm noktasında 190 s sürede 96.92 °C, z=20<br />

101


mm noktasında 193 s sürede 107.45 °C, z=40 mm noktasında 201 s sürede 111.08<br />

°C ve z=60 mm noktasında 202 s sürede 113.27 °C olarak ölçülmüştür. 8 bar şarj<br />

basıncında maksimum sıcaklıklar z=5 mm noktasında 188 s sürede 107.58 °C, z=20<br />

mm noktasında 188 s sürede 118.5 °C, z=40 mm noktasında 272 s sürede 120.37 °C<br />

ve z=60 mm noktasında 284 s sürede 121.52 °C olarak ölçülmüştür.<br />

Reaktör 1’de şarj basıncı 4, 6 ve 8 bar olarak değiştirilmiştir. Bu basınçlarda<br />

depolanan hidrojen miktarı sırasıyla 1.01, 1.17 ve 1.39 g olarak gerçekleşmiştir.<br />

Hidrojen şarj basıncıyla orantılı olarak ekzotermik reaksiyon hızlandığından dolayı<br />

depolanan hidrojen miktarı da artmıştır.<br />

Depolanan hidrojen miktarının belirlenmesi amacıyla yapılan deneylerde Reaktör<br />

2’de şarj basıncı 4, 6 ve 8 bar olarak değiştirilmiştir. Bu basınçlarda depolanan<br />

hidrojen miktarı sırasıyla 1.13, 1.25 ve 1.55 g olarak gerçekleşmiştir. Hidrojen şarj<br />

basıncıyla orantılı olarak ekzotermik reaksiyon hızlandığından dolayı depolanan<br />

hidrojen miktarı da artmıştır. Hidrür oluşumu soğutulan reaktörün soğutulan<br />

kısımlarında daha hızlı gerçekleşmektedir. Bu nedenle kanatçıklı olarak imal edilen<br />

Reaktör 2’de depolanan hidrojen miktarı Reaktör 1’e göre daha yüksek değerlere<br />

çıkmıştır.<br />

Reaktörlerde depolanan hidrojen kütlesinin YSA kullanılarak tahmin edilebilmesi<br />

için oluşturulan modelde Reaktör 1’de %99.9 ve Reaktör 2’de %99.7 doğrulukta<br />

tahmin edilebileceği görülmüştür. R 2 , RMSE ve MAPE değerleri Reaktör 1 için<br />

sırasıyla 0.999, 0.00006 ve 0.0214 olarak bulunurken, Reaktör 2 için sırasıyla 0.998,<br />

0.01258 ve 0.0042 olarak bulunmuştur. Her iki reaktör için oluşturulan modelden<br />

elde edilen bulgular deneysel verilerle karşılaştırıldığında uyum içerisinde oldukları<br />

görülmüştür.<br />

5.2. Öneriler<br />

LaNi4.75Al0.25 alaşımları doldurulmuş <strong>metal</strong> <strong>hidrür</strong> reaktörlerde hidrojen depolama<br />

parametrelerine <strong>hidrür</strong> reaktör tasarımının etkisini incelemek amacıyla yapılan bu<br />

102


çalışmada, aktivasyon hızı, reaksiyon hızı ve termodinamik özelliklerde çok iyi<br />

sonuçlar alınmakla beraber, daha sonraki çalışmalarda aşağıda önerilen çalışmalar<br />

yapılabilir.<br />

• Deneylerin yapıldığı ortam sıcaklığı sabit tutulmalı ve çevre sıcaklığının<br />

deney sonuçlarını etkilemesi önlenmelidir.<br />

• Reaktör, hidrojen depolanan malzemeye (LaNi5, MgNi2, CaNi5 vb) yakın<br />

veya aynı malzemeden imal edilmeli.<br />

Daha kısa sürede doldurulup boşaltılabilecek reaktör tasarlanması ve imalatı için;<br />

Efektif ısı transfer katsayısını arttırmak için farklı ısı değiştiricisi tasarımları<br />

yapılabilir. Bu tasarımlar aşağıdaki gibi olabilir.<br />

• Đçten veya dıştan kanatçık,<br />

• Đçten veya dıştan soğutucu akışkan sirkülasyonu,<br />

• Yukarıdaki tasarımların birkaçının veya tamamının birlikte uygulanması ile<br />

reaktör içerisinde oluşan ısı hızlı şekilde transfer edilebilecektir.<br />

Reaktörden istenilen debide hidrojenin geri alınması için ise reaktörün ısıtılması<br />

gerekir. Bunun için;<br />

• PCI kontrollü elektrikli ısıtıcı kullanarak optimum ısıl kontrol sağlanabilir.<br />

• Mobil uygulamalarda, taşıtların motor veya egzoz ısısından yararlanarak<br />

reaktörden hidrojeni geri almak için gerekli ısı (endotermik reaksiyon)<br />

sağlanabilir.<br />

103


6. KAYNAKLAR<br />

Aldaş K., Mat, M.D., Kaplan, Y., 2002. A three-dimensional mathematical model for<br />

absorption in a <strong>metal</strong> hydride bed. International Journal of Hydrogen Energy,<br />

27, 1049 – 1056.<br />

Aldaş, K., Mat, M.D., 2002. Metal-<strong>hidrür</strong> yataklarda hidrojen depolanmasının sayısal<br />

analizi. Turkish Journal of Engineering and Environmental Sciences, 26, 201-<br />

207.<br />

Altan, M., Yörükoğulları, E., 1997. Hidrojen Zeolit Sisteminin Enerji<br />

Teknolojisindeki Önemi. Dünya Enerji Konseyi Türk Milli Komitesi, Türkiye<br />

7. Enerji Kongresi.<br />

Askri, F., Ben Salah, M., Jemni, A., Ben Nasrallah, S., 2009. Optimization of<br />

hydrogen storage in <strong>metal</strong>-hydride tanks. International Journal of Hydrogen<br />

Energy, 34(2), 897–905.<br />

Balasubramaniam, R. 1997. Hysteresis in <strong>metal</strong>–hydrogen systems. Journal of Alloys<br />

and Compounds, 253–254, 203–206.<br />

Baldissin, D., Lombardo, D., 2009. Thermofluidynamic modelling of hydrogen<br />

absorption and desorption in a LaNi4.8Al0.2 hydride bed. Excerpt from the<br />

proceedings of the COMSOL Conference, Milan.<br />

Barkhordarian, G., Klassen, T., Bormann, R., 2006. Catalytic mechanisms of<br />

transition <strong>metal</strong> oxides on hydrogen sorption reaction of magnesium. Journal<br />

of Physical Chemistry, B, 110 (22), 110-120.<br />

Beltowska, B. M., Kleperis, J., Wojcik, G., Czerwinski, A., Skowronski, J.,<br />

Kopczyk, M., 2001. Electrochemical behavior of <strong>metal</strong> Hydrides. Journal of<br />

Solid State Electrochemical, 5, 229–249.<br />

Bernauer, O., 1988. Metal hydride technology. International Journal of Hydrogen<br />

Energy, 13(3), 181–190.<br />

Bilgili, M., 2003. Metal <strong>hidrür</strong> yataklarda hidrojen soğurulmasının sayısal analizi,<br />

Gazi Üniversitesi Fen Bilimleri Enstitüsü, Yüksek Lisans Tezi, 144 s,<br />

Ankara.<br />

Carter, T.J., Cornish L.A., 2001. Hydrogen in Metals, Engineering Failure Analysis,<br />

8, 113-121.<br />

Chen, C. H., 1996. Fuzzy Logic and Neural Network Handbook. IEEE Pres.<br />

104


Cheng, H.H., Yang, H.G., Li, S.L., Denga, X.X., Chena, D.M., Yanga, K., 2008.<br />

Effect of hydrogen absorption/desorption cycling on hydrogen storage<br />

performance of LaNi4.25Al0.75. Journal of Alloys and Compounds, 453, 448–<br />

452.<br />

Christmann K., 1981. Hydrogen Adsorption on Metal Surfaces. In: Atomistics of<br />

Fracture Conference Proceedings, Eds. Latanision and Pickensr, Plenum, NY,<br />

USA.<br />

Cui, N., Luo, J. L., Chuang, 2001. Study of hydrogen diffusion in α- and β-phase<br />

hydrides of Mg2Ni alloy by microelectrode technique, Journal of<br />

Electroanalytical Chemistry, 503(1–2), 92–98.<br />

Da Li, C., Hong-hui, C., Lei, M., De-min, C., Man-qi, L., Ke, Y., 2007. Effects of Al<br />

partial substitution for Ni on properties of LaNi5−xAlx. Transactions of<br />

Nonferrous Metals Society of China, 17, 967–971.<br />

Demiralp, M., Mat, M.D., Kaplan, Y. Bayrak, M., Veziroğlu, T.N., 2004. Metal<br />

<strong>hidrür</strong> yataklarda ısı ve kütle transferinin deneysel ve teorik olarak<br />

incelenmesi. Niğde Üniversitesi Müh. Bilimleri Dergisi, 8(1- 2), 13–22.<br />

Demircan, A., Demiralp, M., Kaplan,Y., Mat, M.D., Veziroglu, T.N., 2005.<br />

Experimental and theoretical analysis of hydrogen absorption in LaNi5–H2<br />

reactors. International Journal of Hydrogen Energy, 30, 1437 – 1446.<br />

Dhaou, H., Souahlia, A., Mellouli, S., Askri, F., Jemni, A., Nasrallah, S.B., 2010.<br />

Experimental study of a <strong>metal</strong> hydride vessel based on a finned spiral heat<br />

exchanger. International journal of hydrogen energy, 35, 1674–1680.<br />

E4 Tech, Đnternet Sitesi. www.e-sources.com, Erişim Tarihi:05.01.2010.<br />

Efe, M. Ö., Kaynak, O., 1999. A comparative study of neural network structures in<br />

ıdentification of nonlinear systems. Mechatronics, 9(3), 287–300.<br />

Efe, M. Ö., Abadoğlu, E., Kaynak, O., 1999. Analysis and design of a neural network<br />

assisted nonlinear controller for a bioreactor. International Journal of Robust<br />

and Nonlinear Control, 9(11), 799–815.<br />

Elmas, Ç., 2003. Yapay Sinir Ağları (Kuram, Mimari, Eğitim, Uygulama). Seçkin<br />

Yayıncılık, 192s. Ankara.<br />

Fischer P., Yvon K., 1988. Crystal and magnetic structures of ternary <strong>metal</strong> hydrides:<br />

a comprehensive review. Hydrogen in Inter<strong>metal</strong>lic Compounds I, Springer,<br />

63, 87-138.<br />

Flanagan, T.B., Oates, W.A. 1988. Thermodynamics of Inter<strong>metal</strong>lic Compound-<br />

Hydrogen Systems. Hydrogen in Inter<strong>metal</strong>lic Compounds I (editor L.<br />

Schlapbach), Springer Verlag, p. 49, Berlin.<br />

105


Flanagan, T.B., Park, C.N., Everett, D.H. 1987. Hysteresis in <strong>metal</strong> hydrides: an<br />

illustration of entropy production. Journal of Chemical Education, 64, (11),<br />

944-946.<br />

Fuel Cell Store, Đnternet Sitesi. www.fuelcellstore.com, Erişim Tarihi:01.01.2010.<br />

Gerard, N., Ono, S. 1992. Hydride formation and decomposition kinetics, in<br />

Hydrogen in Inter<strong>metal</strong>lic Compounds II, Vol. 67, Springer- Verlag, Berlin,<br />

Chapter 4.<br />

Gopal, M.R., Murthy, S.S, 1992. Prediction of heat and mass transfer in annular<br />

cylindrical <strong>metal</strong> hydride beds. International journal of Hydrogen Energy,<br />

17(10), 795–805.<br />

Gosselink, J.W., 2002. Pathways to a more sustainable production of energy:<br />

sustainable hydrogen-a research objective for shell. International Journal of<br />

Hydrogen Energy, 27(11-12), 1125-1129.<br />

Mendelsohn, M., Gruen, D., Dwight, A., 1979. Group 3A and 4A substituted AB5<br />

hydrides. Inorganic Chemistry, 18(12), 3343–3345.<br />

Güvendiren, M., Ünalan, H.E., Öztürk, T., 2002. Hidrojen depolama amacıyla<br />

magnezyum tozlarının öğütülmesinde katkı maddelerinin etkisi. Toz<br />

Metalürjisi Konferansı.<br />

Hanada, N., Ichikawa, T., Hino, S., Fujii, H., 2006. Remarkable improvement of<br />

hydrogen sorption kinetics in magnesium catalyzed with Nb2O5. Journal of<br />

Alloys and Compounds, 420(1–2), 46–49.<br />

Hangström, M., 1999. Improved <strong>metal</strong> hydrides for energy applications based on<br />

gas-solid reactions. Doctor’s dissertation, Helsinki University of Technology<br />

Publication in Engineering Physics, TKK-F-A793, Finland.<br />

Ham, M.F., Kostanic, I., 2001. Principles of Neurocomputing for Science &<br />

Engineering. 643p. McGraw-Hill, Inc. New York, U.S.A.<br />

Haykin, S., 1994. Neural Networks. Macmillan College Printing Company, New<br />

Jersey.<br />

Hirscher, M., Becher, M., Haluska, M., Quintel, A., Skakalova, V., Choi, Y.M.,<br />

Dettlaff-Weglikowska, U., Roth, S., Stepanek, I., Bernier, P., Leonhardt, A.,<br />

Fink, J., 2002. Hydrogen storage in carbon nanostructures. Journal of Alloys<br />

and Compounds, 330, 654–658.<br />

Hirscher, M., 2010. Handbook of Hydrogen Storage: New Materials for Future<br />

Energy Storage. Wiley-Vch, Verlag, 353pp., Weinheim, Germany.<br />

106


Hoffman, K.C., 1973. Economics of Hydrogen Energy Systems. Brookhaven<br />

National Laboratory Upton, New York.<br />

Hottinen, T., 2001. Technical Review and Economic Aspects of Hydrogen Storage<br />

Technologies, Helsinki University of Technology, Department of<br />

Engineering Physics and Mathematics, Master's thesis, 81 p., Helsinki.<br />

International Energy Agency, Đnternet Sitesi. www.iea.org, Erişim Tarihi:16.02.2011<br />

Ivey, D.G., Northwood, D.O., 1983. Storing energy in <strong>metal</strong> hydrides: a review of<br />

the physical <strong>metal</strong>lurgy. Journal of Materials Science, 18 (2), 321-347<br />

Jemni A., Ben Nasrallah S., 1995. Study of two dimensional heat and mass transfer<br />

during absorption in a <strong>metal</strong>- hydrogen reactor. International Journal of<br />

Hydrogen Energy, 20(1), 43–52.<br />

Jemni, A., Nasrallah, S.B., Lamloumi, J., 1999. Experimental and theoretical study<br />

of a <strong>metal</strong>-hydrogen reactor. International Journal of Hydrogen Energy, 24,<br />

631–644.<br />

Johari, A., Yunus, R.M., Noordin, M.S., Ghazali, M.N.F.Md., Abdullah, T.A.T.,<br />

2001. MmNi5 and MmNi4.85Fe0.15 as hydrogen storage medium. Journal<br />

Technology, 35(F) 1–10.<br />

Kaplan, Y., 2009. Effect of design parameters on enhancement of hydrogen charging<br />

in <strong>metal</strong> hydride reactors. Int. Journal of Hydrogen Energy, 34, 2288 – 2294.<br />

Kaya, S., 2005. Magnezyum Bor Hidrür Sentezi ve Hidrojen Çeviriminde Kullanımı,<br />

Gazi Üniversitesi Fen Bilimleri Enstitüsü Yüksek Lisans Tezi, Ankara.<br />

Kınacı, A, First principles investigation of hydrogen storage in inter<strong>metal</strong>lic systems.<br />

Graduate School of the Natural and Applied Sciences, Middle East Technical<br />

University, Master Thesis, 89p., Ankara.<br />

Kim, J.K., Park, I.S., Kim, K.J., Gawlik, K., 2008. A hydrogen-compression system<br />

using porous <strong>metal</strong> hydride pellets of LaNi5-xAlx. International journal of<br />

hydrogen energy, 33, 870–877.<br />

Lee, S., Lee, Y., 2002. Hydrogen Storage in single-walled carbon nanotubes. Applied<br />

Physics Letters, 76 (20), 2877–2899.<br />

Libowitz, G.G., 1991. Presented at the hydrogen storage materials, batteries and<br />

electrochemistry. Electrochemical Society Proceedings, 92–5, 3.<br />

Libowitz, G.G., 1965. The Solid-State Chemistry of Binary Metal Hydrides. W.A.<br />

Benjamin, New-York.<br />

Linde Gas, www.linde-gas.com, Erişim Tarihi:16.10.2010.<br />

107


Liu, F.J., Sandrock, G., Suda, S., 1992. Activation characteristics of chemically<br />

treated LaNi4.7Al0.3. Journal of Alloys and Compounds, 190(1), 57–60.<br />

Liu, F.J., Suda, S., Sandrock, G., 1996. Effects of Ni-substitution and F-treatment on<br />

the hydriding behaviors and microstructures of AB2-compound<br />

(Ti,Zr)(Mn,Cr)2. Journal of Alloys and Compounds, 232(1–2), 232–237.<br />

Lucas, G.G., Richards, W.L. 1984. Mathematical modeling of hydrogen storage<br />

systems. International Journal of Hydrogen Energy, 9(3), 225–231.<br />

MacDonald, B.D., Rowe, A.M., 2006. A thermally coupled <strong>metal</strong> hydride hydrogen<br />

storage and fuel cell system. Journal of Power Sources, 161(1), 346–355.<br />

MacDonald, B.D., Rowe, A.M., 2007. Experimental and numerical analysis of<br />

dynamic <strong>metal</strong> hydride hydrogen storage systems. Journal of Power Sources<br />

174, 282–293.<br />

Maeland, A.J. 1977. Presented at the Hydrides for Energy Storage. Geilo, Norway,<br />

Pergamon.<br />

Mat, M.D., Kaplan, Y., 2001. Numerical study of hydrogen absorption in an Lm−Ni5<br />

hydride reactor. International Journal of Hydrogen Energy, 26, 957–963.<br />

Mayer, U., Groll, M., Supper, W., 1987. Heat and mass transfer in <strong>metal</strong> hydride<br />

reaction beds: Experimental and theoretical results. Journal of Less-Common<br />

Metals, 131, 235–244.<br />

Mellouli, S., Askri, F., Dhaou, H., Jemni, A., Ben Nasrallah, S., 2010. Numerical<br />

simulation of heat and mass transfer in <strong>metal</strong> hydride hydrogen storage tanks<br />

for fuel cell vehicles. Int. Journal of Hydrogen Energy, 35(4), 1693–1705.<br />

Melloulia, S., Askri, F., Dhaou, H., Jemni, A., Nasrallah, S.B. 2007. A novel design<br />

of a heat exchanger for a <strong>metal</strong>-hydrogen reactor. International Journal of<br />

Hydrogen Energy, 32, 3501–3507.<br />

Melnichuk, M., Andreasen, G., Corso, H.L., Visintin, A., Peretti H.A., 2008. Study<br />

and characterization of a <strong>metal</strong> hydride container. International Journal of<br />

Hydrogen Energy, 33, 3571–3575.<br />

Mintz, M.H., Bloch, J., 1985. Evaluation of the kinetics and mechanisms of<br />

hybriding reactions. Progress in Solid State Chemistry, 16(3), 163–194.<br />

Mohan, G., Maiya, M.P., Murthy, S.S., 2007. Performance simulation of <strong>metal</strong><br />

hydride hydrogen storage device with embedded filters and heat exchanger<br />

tubes. International Journal of Hydrogen Energy, 32, 4978–4987.<br />

108


Muthukumar, P., Manvendra M. U, 2009. Study of coupled heat and mass transfer<br />

during absorption in a MmNi4.6Al0.4 based hydrogen storage device. Sadhana,<br />

34, 1–16.<br />

Muthukumar, P., Madhavakrishna, U., Anupam D., 2007. Parametric studies on a<br />

<strong>metal</strong> hydride based hydrogen storage device. International Journal of<br />

Hydrogen Energy, 32, 4988–4997.<br />

Nakamura, Y., Akiba, E. 2000. In-situ X-ray diffraction study on LaNi5 and<br />

LaNi4.75Al0.25 in the initial activation process. Journal of Alloys and<br />

Compounds, 308(1–2), 309–318.<br />

Narenda, K.S., Pathasaraty, K., 1990. Identification and control of dynamical<br />

systems using neural networks. IEEE Transactions on Neural Networks, 1(1),<br />

4–27.<br />

Nasrallah S.B., Jemni, A. 1997. Heat and mass transfer models in <strong>metal</strong> hydrogen<br />

reactor. International Journal of Hydrogen Energy, 22(1), 67–76.<br />

Oi, T, Maki, K, Sakaki, Y, 2004. Heat transfer characteristics of the <strong>metal</strong> hydride<br />

vessel based on the plate-fin type heat exchanger. Journal of Power Sources,<br />

125, 52–61.<br />

Öztemel, E., 2003. Yapay Sinir Ağları. Papatya Yayıncılık, 232s, Đstanbul.<br />

Park, S., Kim, J.K., Kim, K.J., Zhanga,J., Park, C., Gawlik K., 2009. Investigation of<br />

coupled AB5 type high-power <strong>metal</strong> hydride reactors. International journal of<br />

hydrogen energy, 34, 5770–5777.<br />

Park, C.N., Luo, S., Flanagan, T.B., 2004. Analysis of sloping plateaux in alloys and<br />

inter<strong>metal</strong>lic hydrides: I. Diagnostic features. Journal of Alloys and<br />

Compounds, 384(1–2), 203–207.<br />

Percheron-Guegan, A., Lartigue,C., Achard, J.C., Germi, P., Tasset, F., 1980.<br />

Neutron and X-ray diffraction profile analyses and structure of LaNi5,<br />

LaNi5−xAlx and LaNi5−xMnx inter<strong>metal</strong>lics and their hydrides (deuterides).<br />

Journal of the Less Common Metals, 74(1), 1–12.<br />

Pottier, J., Bailleux, C., 1986. Hydrogen: A Gas Of The Past, Present and Future A<br />

Drama In 3 Acts. Hydrogen Energy Progress, 1, 197-216.<br />

Reilly J. J., Wiswall, R. H. 1974. Formation and properties of iron titanium hydride.<br />

Journal of Inorganic Chemical, 13, 218–222.<br />

Sandrock G., 1999. A panoramic overview of hydrogen storage alloys from a gas<br />

reaction point of view. Journal of Alloys and Comp., 293–295, 877- 888.<br />

109


Sandrock, G., Suda, S., Schlapbach, L. 1992. Applications, Topics in Applied<br />

Physics. Hydrogen in Inter<strong>metal</strong>lic Compounds II, Springer-Verlag, Berlin, p.<br />

197–258.<br />

Satheesh, A., Muthukumar, P., Dewan, A., 2009. Computational study of <strong>metal</strong><br />

hydride cooling system. International Journal of Hydrogen Energy, 34(7),<br />

3164–3172.<br />

Schlapbach L., 1988. Introduction. Hydrogen in Inter<strong>metal</strong>lic Compounds I, Vol.63,<br />

pp.1-10, Springer, Germany.<br />

Schwarz, R.B., Khachaturyan, A.G. 2006. Thermodynamics of open two-phase<br />

systems with coherent interfaces: Application to <strong>metal</strong>–hydrogen<br />

systems. Acta Materialia, 54(2), 313–323.<br />

Selvam P., Viswanathan B., Swamy, C.S., Srinivasan, V., 1986. Magnesium and<br />

Magnesium alloy hydrides. International Journal of Hydrogen Energy, 11(3),<br />

169-192.<br />

Shimizu, E., Aoki, K., Masumoto, T., 1999. Hydrogen absorption properties of<br />

amorphous and crystalline alloys in the pseudobinary ZrCo–TiNi system.<br />

Journal of Alloys and Compounds, 293–295, 526–530.<br />

Sirosh, N., 2002. Hydrogen Composite Tank Program. Proceedings of the U.S. DOE<br />

Hydrogen Program Review, NREL/CP–610–32405.<br />

Sun, D.V., Deng, S.J., 1990. Numerical solution of the two-dimensional non-steady<br />

heat and mass transfer problem in <strong>metal</strong> hyride beds. International Journal of<br />

Hyrogen Energy, 15(11), 807–816.<br />

Teitel, R. 1981. Hydrogen storage in glass microspheres. Report BNL 51439,<br />

Brookhaven National Laboratories.<br />

Türe, Đ. E., 2001. Hidrojen Enerjisi, ISBN-975-8547-00-3. Temiz Enerji Vakfı<br />

Yayınları, TUBĐTAK Matbaası.<br />

Yamaguchi, M., Akiba, E., 1994. Ternary hydrides, in electronic and magnetic<br />

properties of Metals and Ceramics. Part II, vol. 3 (ed. K.H. Jürgen Buschow),<br />

VCH, Weinheim, p. 333.<br />

Zaluska, A., Zaluski, L., Ström-Olsen J.O., 2001. Structure, catalysis and atomic<br />

reactions on the nano-scale: a systematic approach to <strong>metal</strong> hydrides for<br />

hydrogen storage. Applied Physics A Materials Science & Processing, 72,<br />

157–165.<br />

Zuttel, A., 2003. Materials for hydrogen storage. Materials Today, pp.24-33.<br />

110


EKLER<br />

EK–1 Kullanılan Alaşım ve Gazların Kimyasal-Fiziksel Özellikleri<br />

Çizelge 1. Lantanyum nikelin (LaNi) fiziksel ve kimyasal özellikleri<br />

Özellik Değer<br />

Yoğunluk, ρ 4230 kg/m 3<br />

Özgül ısı, cp<br />

111<br />

419 J/kgK<br />

Isı iletim katsayısı, k, 1.6 W/mK<br />

Aktivasyon enerjisi, Ea 25817 J/mol H2<br />

Reaksiyon enerjisi, H<br />

Soğutucu akışkan ile ısı taşınım katsayısı, h<br />

10 J/mol H2<br />

500 W/m 2 K<br />

Hidrojen ile reaktör yüzeyi arasındaki ısı 50 W/m 2 K<br />

taşınım katsayısı, h<br />

Gözeneklilik, ε, (%) 50.5<br />

Gaz sabiti, (RH2) (J/kgK) -<br />

Referans kinematik viskozite, υ (m 2 /s) -<br />

Dinamik viskozite, µ, (Ns/m 2 ) -<br />

Isı yayılım sabiti, α, (m 2 /s) 87.6x10 -6<br />

Çizelge 2. Azot gazının fiziksel özellikleri<br />

Özellik Değer<br />

Erime Noktası (1 atm) -209.9 °C<br />

Kaynama Noktası (1 atm) -195.8 °C<br />

Yoğunluk (gaz, 0°, 1 atm) 1.2506 g/m 3<br />

Moleküler Ağırlık 28.01 g/mol-g<br />

Yoğunluk (likit, 1 atm) 0.8082 kg/l<br />

Spesifik Isı 0.249<br />

Buharlaşma Isısı 199 kJ/kg<br />

Yoğunluk, gaz (15 °C, 1 atm) 1.185 kg/m 3<br />

Özgül ağırlık, gaz (Hava=1) 0.967<br />

Kritik Sıcaklık -147 °C<br />

Kritik Basınç 33.5 atm


Çizelge 3. Hidrojen gazının fiziksel özellikleri<br />

Özellik Değer<br />

Yoğunluk 0.084 kg/m 3<br />

Erime Noktası -259.14 °C<br />

Kaynama Noktası -252.87 °C<br />

Hava içindeki difüzyonu, 0.61 cm 2 /s<br />

Sabit basınçta özgül ısısı, 14.89 J/gK<br />

Havada ateşlenme sınırı, (%hacim) 4.0–75.0<br />

Havada ateşlenme enerjisi, 0.02 mJ<br />

Ateşlenme sıcaklığı, 585 °C<br />

Havada alev sıcaklığı, 2045 °C<br />

Patlama enerjisi, ( 0.17 g TNT k/j<br />

Alev hızı 2.91 m/s<br />

Çizelge 4. Vokum pompası teknik özellikleri<br />

Özellik Değer<br />

Silindir hacmi (Süpürme hacmi) 50 Hz 1.8 m 3 /h<br />

Hız 50 Hz 1.6 m 3 /h<br />

Maksimum vakum (toplam basınç) 3.10 -3 mbar<br />

Bağlantı 11 mm dış çapında (3/8)<br />

Hızlı sökülebilir bağlantı<br />

Maksimum çıkış basıncı 0.5 bar<br />

Maksimum giriş basıncı (su buharı için) 15 mbar<br />

Maksimum su buharı pompalama hızı 16 g/h<br />

Motor gücü 50/60 Hz 0.16 kW<br />

Çalışma sıcaklıkları 12–40 °C<br />

Ağırlık 10 kg<br />

Ses 54 dB<br />

Yağ kapasitesi-maksimum<br />

Yağ kapasitesi-minimum<br />

112<br />

0.28 lt<br />

0.20 lt


EK-2 Deneysel Sistem Resimleri<br />

Şekil 1. LaNi4.75Al0.25 alaşımı (toz)<br />

Şekil 2. Reaktör 1 deney tesisatı<br />

113


Şekil 3. Reaktör 2 deney tesisatı<br />

Şekil 4. Mekanik öğütücü<br />

114


Şekil 5. H2 ve N2 manometre ve kontrol vanaları<br />

115


Şekil 6. Termokuplların reaktör üzerine yerleşimi<br />

116


ÖZGEÇMĐŞ<br />

Adı Soyadı : Muhammet KAYFECĐ<br />

Doğum Yeri ve Yılı : Anamur, 1981<br />

Medeni Hali : Evli<br />

Yabancı Dili : Đngilizce<br />

Eğitim Durumu (Kurum ve Yıl)<br />

Lise :Bozyazı Lisesi, Mersin-Bozyazı, 1996-1999<br />

Lisans : Zonguldak Karaelmas Üniversitesi Karabük Teknik Eğitim Fakültesi<br />

Makina Eğitimi Bölümü, 1999-2003<br />

Yüksek Lisans :Zonguldak Karaelmas Üniversitesi Fen Bilimleri Enstitüsü Makina<br />

Eğitimi Anabilim Dalı, 2003-2005<br />

Çalıştığı Kurum/Kurumlar ve Yıl: Karabük Üniversitesi, (2008-Halen)<br />

Yayınları (SCI ve diğer makaleler):<br />

1. Kurt, H. and Kayfeci, M., 2009. Prediction of thermal conductivity of<br />

ethylene glycol–water solutions by using artificial neural networks, Applied<br />

Energy, 86(10), 2244–2248.<br />

2. Keçebaş, A. and Kayfeci, M., 2010. Effect on optimum insulation thickness,<br />

cost and saving of storage design temperature in cold storage in Turkey,<br />

Energy Education Science and Technology, Part A, Energy Science and<br />

Research, 25(1-2), 117–127.<br />

3. Keçebaş A., Kayfeci, M. and Gedik, E., Performance investigation of the<br />

Afyon geothermal district heating system for building applications: Exergy<br />

analysis, Applied Thermal Engineering, 31, 1229-1237, 2011.<br />

4. Keçebaş A. ve Kayfeci, M., 2009. Mekanik Beslemeli Kazanların Otomatik<br />

Kontrolü, 5. Uluslararası Đleri Teknolojiler Sempozyumu (IATS’09), s.1705-<br />

1709, 13-15 Mayıs, Karabük, Türkiye.<br />

5. Kayfeci, M., Gedik, E., Sağıroğlu, S. ve Kurt, H., 2009. Taşıt<br />

Đklimlendirmesinde Alternatif Klima Sistemlerinin Kullanımının<br />

Araştırılması, 5. Uluslararası Đleri Teknolojiler Sempozyumu (IATS’09),<br />

s.1882–1886, 13–15 Mayıs, Karabük, Türkiye.<br />

117


6. Gedik, E., Kayfeci, M., Keçebaş A. and Kurt, H., 2009. Magnetic<br />

Refrigeration Technology Applications On Near-room Temparature, 5th<br />

International Advanced Technologies Symposium (IATS’09), pp.1887-1893.<br />

May 13-15, Karabuk, Turkey.<br />

7. Keçebaş, A., Gedik, E. and Kayfeci, M., 2010. Investigation of the energy<br />

and exergy performance of a solar ejector cooling system using various<br />

working fluids. 5th International Ege Energy Symposium and Exhibition<br />

(IEESE-5), 27-30 June, Pamukkale University, Denizli, Turkey.<br />

8. Kayfeci, M. ve Kurt, H., 2007. Sıvılar Đçin Isı Đletim Katsayısı Ölçüm<br />

Cihazının Tasarımı, Đmali ve Test Edilmesi, Gazi Üniversitesi. Müh. Mim.<br />

Fak. Der. 22(4), 907–915.<br />

9. Gürbüz, H. ve Kayfeci, M., 2008. Otomobillerde Metal Hidrid Esaslı Klima<br />

Sistemlerinin Kullanılabilirliğinin Araştırılması, Mühendis ve Makina,<br />

49(584), 20-24.<br />

10. Yılmaz, S., Kayfeci, M. ve Keçebaş, A., 2010. Çift Fazlı (Isı Borulu) Güneş<br />

Kollektörlerini Konut Isıtmasında Kullanımının Deneysel Olarak<br />

Đncelenmesi, Tesisat Mühendisliği Dergisi, Sayı 116 Mart/Nisan, Sayfa:5-10.<br />

11. Keçebaş, A., Gedik, E. ve Kayfeci, M., 2010. Fosil Yakıtların Kullanımından<br />

Kaynaklanan Hava Kirliliği Üzerine Jeotermal Enerji ve Doğalgaz<br />

Kullanımının Etkisi: Afyon Örneği, Makine Teknolojileri Elektronik Dergisi,<br />

7(3), 23-30.<br />

12. Keçebaş, A., Kayfeci, M. and Gedik, E., 2010. Afyonkarahisar’da ölçülen<br />

sülfür dioksit ve partikül madde değişimleri üzerine bir değerlendirme,<br />

Mühendis ve Makine Dergisi, 51(610), 14-19.<br />

13. Kayfeci, M., Kurt, H. ve Bedir F., 2009. Metal Hidrür Esaslı Isı Pompalarının<br />

Đklimlendirmede Kullanımının Araştırılması, IV. Ulusal Hidrojen Enerjisi<br />

Kongresi ve Sergisi (Uluslararası Katılımlı) Bildiriler Kitabı, Sayfa 57-64,<br />

15-16 Ekim, Kocaeli, Türkiye.<br />

14. Keçebaş, A. ve Kayfeci, M., Alternatif Evsel Klima Sistemlerinin<br />

(Termoelektrik, Absorpsiyonlu ve Metal Hidrid) Buhar Sıkıştırmalı Sistemle<br />

Karşılaştırılması, 2008. VII. Ulusal Temiz Enerji Sempozyumu Bildiriler<br />

Kitabı, Sayfa: 743-751, 17–19 Aralık, Đstanbul, Türkiye.<br />

118

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!