)( ) TR - Gıda Teknolojileri Elektronik Dergisi

1. Giriş
www.teknolojikarastirmalar.com
ISSN:1304-4141
Makine Teknolojileri Elektronik Dergisi
2007 (4) 31-37
Kısa Makale
Van Der Waals Gazlarına Ait Genel Đş Formülleri
TEKNOLOJĐK
ARAŞTIRMALAR
Zekai Kazım TELLĐ 1 , Đbrahim ÜÇGÜL 2 , Murat ÖZTÜRK 2
1 Süleyman Demirel Üniversitesi Mühendislik Mimarlık Fakültesi ISPARTA
2 Süleyman Demirel Üniversitesi Yenilenebilir Enerji Kaynakları Araştırma ve Uygulama Merkezi ISPARTA
Özet
Gerçek gazların yüksek basınçlara sıkıştırılması, kimya sanayisinde bazı yüksek basınçlı reaksiyonlar için
gereklidir. Sıkıştırma olaylarında iş ve teknik iş formülleri izoterm, izentropik ve izobar durum değişimleri için
türetilmesi gerekmektedir.
Sunulan çalışmada, Van der Waals denklemi ve Termodinamiğin I. Yasası kullanılarak, gerçek gazların
izoterm, izentropik ve izobar sıkıştırma işlemleri için geçerli iş ve teknik iş bağıntıları türetilmiştir. Bilindiği
üzere mol kütlesi büyük gazlar, yüksek basınç değerlerinde, ideal tutumdan önemli ölçüde farklı olmakta ve
belirgin gerçek gaz özelliği göstermektedir.
Anahtar Kelimeler: Van der Waals gazları, Đş formülleri, Durum değişimleri
Basınç, sıcaklık ve hacim arasındaki bağıntının deneysel incelenmesinde; ideal gaz denkleminin birçok gaz için
gerçek şartlara uymadığı görülmüştür. Basınç ne kadar büyük ve sıcaklık ne kadar küçük olursa, sapma o derece
büyük olmaktadır. Ancak basıncın düşük ve sıcaklığın büyük olması durumunda ideal gaz denklemlerin yaklaşma
görülmektedir. Đdeal şartları, sapma individüel olduğundan belli basınç ve sıcaklık değerleriyle belirlemek mümkün
değildir. Newton’un çekim yasası göz önüne alındığında, moleküller arası çekim kuvvetlerin uzaklığın karesine ters
orantılı ve moleküllerin mol kütlelerine doğru orantılı olduğu görülür. Gaz ne kadar yoğun ve molekül kütlesi ne
kadar büyük olursa, moleküller arası çekim kuvvetleri o kadar etkindir. Dolayısıyla gaz gerçek tutumludur [1].
Endüstriyel proseslerin çoğunda gazların yüksek basınç değerlerine sıkıştırılması gerekmektedir. Sıkıştırma çok
kademeli olarak gerçekleştirilmekte, ileri kademelerde gazlar ideal gaz özelliğini kaybetmekte ve gerek gaz tutumu
göstermektedirler. Gazların ideal tutumdan sapmasında iki özellik etkilidir.
a-) Moleküller arası çekim kuvvetleri
b-) Đndividüel mol hacmi özellikleri
Basıncın yüksek olması; belirli hacimdeki molekül sayısının çok olduğunu göstermektedir. Dolayısıyla moleküller
arasındaki uzaklık küçük, çekim kuvvetleri etkin olmaktadır. Aynı zamanda düşük sıcaklıklarda moleküllerin
kinetik enerjisi düşüktür ve hareket çekim kuvvetlerinin etkisi altında gerçekleşmektedir. Gerçek gazlarda
gözlemlenen bu özellikleri kullanarak Van der Waals, ideal gaz denklemini aşağıdaki gibi tahsis etmiştir.
( P + Pi
)( V − b)
= RT
Burada, Pi çekim kuvvetlerinden ileri gelen moleküler iç basınçtır ve mol hacminin karesi ile ters orantılıdır.
(1)

Teknolojik Araştırmalar: MTED 2007 (4) 31-37 Van Der Waals Gazlarına Ait Genel Đş Formulleri
a
Pi =
V
2
Denklem (1)’in sol tarafında ki ikinci terim mol hacminin b büyüklüğü ile oluşumunu belirlemektedir. Denklem
(1)’de ki b, Denklem (2)’deki a katsayısı Van der Waals sabitlerini göstermektedir ve her gaz için farklı değerleri
alırlar [2]. Sonuç olarak Van der Waals denklemi;
( V − b)
RT
⎛ a ⎞
⎜ P ⎟
⎜
+
⎟
=
2
⎝ V ⎠
şeklini alır. Çok yüksek basınç değerlerinde a ve b basıncın fonksiyonudur (a,b = f(P)). Gerçek gazların izoterm
sıkıştırılması fugasiteler üzerinden belirlenir. Teknik iş;
Wt12 2
∫
= − VdP
1
şeklinde olup, sürekli akış ve sürekli açık sistemler için SASA işidir. Serbest entalpi farkı alınarak hesaplanmalıdır.
Sarf edilen teknik iş aşağıdaki ifade ile belirlenir.
Wt12 = -∆GT (5)
Ancak gerçek gazların sıkıştırılması zamanın küçük kesrinde olması nedeniyle, gerçek gazların durum değişimi
belirleyen denklemlerin bilinmesine ihtiyaç vardır [2].
Sunulan çalışmada, Van der Waals denklemi ve Termodinamiğin I. yasası kullanılarak gerçek gazların izoterm,
izentropik ve izobar durum değişimleri ayrıntılı olarak incelenmiştir. Bu durum değişimleri için geçerli iş ve teknik
iş bağıntıları kullanımları son derece basit olacak şekilde nümeriksel olarak türetilmiştir.
2. Đş ve Teknik Đş Formüllerinin Türetilmesi
2.1. Đzoterm Durum Değişimi
Herhangi bir termodinamik işlem sırasında sıcaklık değişmiyor ise bu işleme sabit sıcaklıkta (izoterm) durum
değişimi denir. Termodinamikte doğrudan izoterm durum değişimi mevcut değildir. Đzoterm sıkıştırma işlemi,
adyabatik sıkıştırma ve ardından soğutma şeklinde olmaktadır. Bir kompresörde izoterm sıkıştırmanın sağlanması,
ancak durum değişimi çok yavaş gerçekleştiği zaman sağlanabilir.
2.1.1. Đzoterm Durum Değişimde Đş Formülü
Literatürde işin tanımını veren pek çok ifade mevcuttur. [3] nolu literatürde işin tanımı,
W
= ∫ PdV
şeklinde verilmiştir. (3) nolu Van der Waals denkleminden basınç;
RT
a
P = −
V − b V
2
şeklinde ifade edilebilir. Bu ifadenin V ’ye göre türevi alınırsa;
32
(2)
(3)
(4)
(6)
(7)

Telli, Z.K., Üçgül, Đ., Öztürk, M.
PdV
=
RT
−1
−2
( V − b)
dV
− aV
dV
33
Teknolojik Araştırmalar: MTED 2007 (4) 31-37
elde edilir. Türetilen bu ifade Denklem (6)’da yerine yazılır ve gerekli işlemler yapılarak düzenlenirse;
⎛ V − b ⎞ a ⎛V
o ⎞
W = RT
ln⎜
⎟ + ⎜ ⎟
⎜ ⎟ ⎜
−1
⎟
⎝V
o − b ⎠ V o ⎝ V ⎠
V o =
ifadesi elde edilir. Burada molar hacimler
ifade edilir [2].
RT
P
o
ve
Z RT
V =
P
2.1.2. Đzoterm Durum Değişiminde Teknik Đş Formülü
Teknik iş [3];
W t
= −∫V
dP
(8)
(9)
şeklinde tanımlanmıştır. Z, real faktör olarak
şeklinde verilmiştir. Van der Waals denkleminden hareketle, Denklem (7) ifadesinin diferansiyeli alınrsa,
dP = −RT
−2
−3
( V − b)
dV
+ 2aV
dV
elde edilir. Bu ifadeyi Denklem (10)’da verilen teknik iş formuna dönüştürmek için denklemin her iki tarafı ( −V
)
ile çarpılarak;
VdV
2adV
−V
dP = RT
− 2 2
V
( V − b)
ifadesi türetilir. Bu ifadenin integralini almak için literatürde [5] verilen Bronstein integral kalıbından hareketle
aşağıdaki ifade türetilir.
−
∫
VdP
=
⎛
RT
⎜
⎝
− b
( V − b)
+ ln
( V − b)
⎞
⎟
⎠
V
Vo
1
+ 2a
V
Bu denklem düzenlenirse real gazlarda teknik iş formülü;
W
t
⎛ V − b ⎞ bRT
⎛V
o − b ⎞ 2a
⎛V
o ⎞
= RT
ln⎜
⎟ − ⎜ ⎟ + ⎜ ⎟
⎜ ⎟ ⎜
−1
⎟ ⎜
−1
−
⎟
⎝V
o − b ⎠ V o b ⎝ V − b ⎠ V o ⎝ V ⎠
şeklinde bulunur.
2.2. Đzentropik Durum Değişimi
V
Vo
Isı alış verişi olmayan, s=sabit bağıntısına uygun olarak değişen işleme izentropik (veya tersinir adyabatik) işlem
denilmektedir.
(10)
(11)
(12)
(13)
(14)

Teknolojik Araştırmalar: MTED 2007 (4) 31-37 Van Der Waals Gazlarına Ait Genel Đş Formulleri
2.2.1. Đzentropik Durum Değişiminde Đş Formülü
(7) noluVan der Waals denklemine, literatürden [4] alınan,
χ′
⎛V
o − b ⎞
T = T ⎜ ⎟
o ⎜ ⎟
⎝ V − b ⎠
denklemi yerine yazılır ve denklemde,
A o
( ) χ′
V −
= RTo
b
(16)
χ ′ =
tanımlaması yapılarak,
A a
P = − χ′
+ 1
V
( V − b)
2
R
C v
ile birlikte,
basınç ifadesi elde edilir. Bu denklemin her iki tarafı d V ile çarpılıp düzenlenirse;
PdV
= A
( V − b)
−
( χ′
+ )
1 −2
dV
− aV
dV
W =
oluşturulur. Đş ifadesi ∫ PdV
’yi türetebilmek için yukarıdaki ifadenin integrali Brostein integral kalıbı
kullanılarak alınırsa aşağıdaki sonuç elde edilir.
∫
PdV
= A
[ − ( χ ′ + 1)
+ 1]
−(
′ + 1)
( V − b)
1 χ
V
+ 1
Bu eşitlik düzenlendiği takdirde aşağıdaki ifade bulunur.
a ⎛V
o ⎞
W = ⎜ ⎟
⎜
−1
⎟
−
V o ⎝ V ⎠ χ ′
A
⎡⎛
V
⎜
Burada A ifadesi yerine yerleştirilirse;
V o
( ) ⎥ ⎥
o
⎢ −
⎢
⎜ ⎟
1
χ′
V o − b ⎝ V − b ⎠
a ⎛V
⎡
o ⎞ RT
⎛ ⎞
o V o − b
W = ⎜ ⎟ − ⎢⎜
⎟
⎜
−1
V
⎟ ′ ⎢
⎜ ⎟
o ⎝ V ⎠ χ
⎣⎝
V − b ⎠
⎣
χ′
− b ⎞
⎟
⎤
−1⎥
⎥
⎦
şeklinde ki izentropik durum değişimindeki iş ifadesi türetilmiş olur.
2.2.2. Đzentropik Durum Değişiminde Teknik Đş Formülü
χ′
Van der Waals denkleminde, Denklem (15) ve (16)’dan türetilen Denklem (17)’nin diferansiyeli;
dP
A
−
( χ ′ + )
= χ′
+ 1
( V − b)
1 dV
−3
+ 2aV
dV
⎤
⎦
34
(15)
(17)
(18)
(19)
(20)
(21)
(22)

Telli, Z.K., Üçgül, Đ., Öztürk, M.
şeklinde bulunur. Teknik iş için eşitliğin her iki tarafı V ile çarpılarak;
VdV
dV
V ( χ ′ + 1)
− 2a
2
V
− dP = A χ′
+ 2
( V − b)
ifadesi elde edilir. bu ifadenin integral çözümü için literatürden [5] alınan,
∫ ( ) ( )( ) ( )( ) ⎟ xdx 1 ⎛ −1
b ⎞
= ⎜
+
n 2
n−2
n−1
ax + b a ⎝ n − 2 ax + b n −1
ax + b ⎠
ifadesinin açılımını kullanarak;
W
t
−1
b
1
= −∫V
dP = A(
χ′
+ 1) −
− 2a
−
χ′
+ 1
χ ′
V
elde edilir. Gerekli düzenlemeler yapılarak;
W
t
2a
⎛V
= ⎜
V
⎜
o ⎝ V
o
χ′
( V − b)
( χ ′ + 1)(
V − b)
⎞ χ ′ + 1 ⎡ ⎛V
o − b ⎞
−1⎟
+ ⎢ − ⎜ ⎟
⎟
RTo
1
′ ⎢
⎜ ⎟
⎠ χ
⎣ ⎝ V − b ⎠
izentropik durum değişimindeki teknik iş eşitliği elde edilir.
2.3. Đzobar Durum Değişimi
35
V
V o
V
V o
⎤ RT
⎡ ⎛ ⎞
ob
V o − b
⎥ + ⎢1
− ⎜ ⎟
⎥ − ⎢
⎜ ⎟
⎦
V o b
⎣ ⎝ V − b ⎠
χ′
χ′
+ 1
Teknolojik Araştırmalar: MTED 2007 (4) 31-37
Teknik pratik de genellikle sürekli çalışan sistemlerde basınç sabittir. Örnek olarak gaz ısıtıcısı veya gaz
soğutucusu verilebilir. Đzobar durum değişiminde sadece hacimsel debi ve sıcaklık değişmektedir.
2.3.1. Đzobar Durum Değişiminde Đş Formülü
P = sabit, dP = 0 olmak üzere izobar durum değişiminde iş;
( V V o )
W = P −
ifadesinden hesaplanır. Burada basınç;
A a
P = − χ′
+ 1
V
( V − b)
2
şeklindedir. Literatürden [4] alınan Denklem (15) ve ayrıca hesaplamalarda kısaltma sağlayan Denklem (16)
3
kullanılabilir. Ayrıca
V o = 22,
414m
kmol
ve V ≈ b [2] arasında incelenen gerçek gazlar için değerler alınarak
tablo oluşturulur ve gerekli değerler bu tablodan alınarak W işi hesaplanır. Tablo 1’de örnek olarak hidrojen
gazının molar hacminin basınç ve sıcaklıkla değişimi verilmiştir.
⎤
⎥
⎥
⎦
(23)
(24)
(25)
(26)
(27)
(28)

Teknolojik Araştırmalar: MTED 2007 (4) 31-37 Van Der Waals Gazlarına Ait Genel Đş Formulleri
V (m 3 /kmol)
Tablo 1. Hidrojen gazı için molar hacmin basınç ve sıcaklıkla değişimi
P(bar) T(K) V (m 3 /kmol)
36
P(bar) T(K)
22,414 1 273 8 433,4 415,6
20 119 286 6 650,6 467,4
18 138,1 298,5 4 1154,9 551,7
16 163 313,2 2 3091,2 733,6
14 196,8 330,7 1 8356 978,1
12 244,6 352,2 0,5 23043 1311,7
10 316,3 379,4 0,25 66307,6 1780,8
Buraya kadar türetilen iş formülleri tüm basınç değerleri için geçerlidir.
2.3.2. Đzobar Durum Değişiminde Teknik Đş Formülü
Đzobar durum değişiminde dP=0 olduğu için
olduğu görülür.
3. Örnek Uygulama
W t
2
∫
= − VdP
1 ifadesinden teknik işin izobar durum değişiminde sıfır
Örnek uygulama olarak izoterm sıkıştırma seçilmiştir. Havanın 1 bardan 300 bara, 300 K sıcaklıkta izoterm olarak
sıkıştırıldığı kabul edilmiştir. Literatürde mevcut olan teknik fugasiteler ve bu çalışmada izoterm sıkıştırma işi için
türetilen kesin formüller yardımıyla hesap uygulaması yapılmış ve sonuçlar karşılaştırılmıştır.
T = 300 K için havanın kritik sıcaklığı (Tkr,Hava) ve kritik basıncı (Pkr,Hava) aşağıdaki şekilde hesaplayabiliriz.
T = T xn + T
xn
kr,
Hava Kr,
N N Kr,
O O
2 2
2
Tkr,H = T kr,
N . 0,
79 + T . 0,
21
2
kr,
O2
= 126.0,79 + 155.0,21
= 132 K
P = P xn + P
xn
kr,
Hava Kr,
N N Kr,
O O
2 2
2
Pkr,H = P kr,
N . 0,
79 + P . 0,
21
2
kr,
O2
= 33,93.0,79 + 50,8.0,21
= 37,5 bar
T* = T/Tkr,H=300/132=2,27
*
P 2 = P Pkr,
H =
P1 = 1 barda, 300 37,
5 8
2
2
=
.
Teknik iş fugasiteler üzerinden W t = − n RT
ln f 2 f1
şeklinde hesaplanır[6]. Burada f i = z fi.
Pi
olup burada i f
.
z i bileşenine ait fugasite ve fugasite kofisyenidir. = 1
ve fi
z 1,
0,
95
1 = f
z f alınarak 1 1 1 P z f = f ve 2 2 2 P z f f
2 =
W t = −8,
315.
300.
ln 285 1 = −14100
kJ
olarak teknik iş fugasiteler üzerinden hesaplanır.
n kmol, R = 8,
315 kj kmolK , literatürden [6],
= için f 1 bar, 285 300 . 95 , 0 f = = ve,
Kesin formüllerle hesaplama yapmak için, bu çalışmada türetilen Denklem (14) kullanılmıştır.
W t
V 2 − b bRT
⎛ V 1 − b ⎞ 2a
⎛ V 1 ⎞
= RT
ln − ⎜ ⎟ + ⎜ ⎟
⎜
−1
⎟ ⎜
−1
V − −
⎟
1 b V 1 b ⎝V
2 − b ⎠ V 1 ⎝V
2 ⎠
1 =
2

Telli, Z.K., Üçgül, Đ., Öztürk, M.
Real faktörler literatürden [2], Z1=1 alınarak V 1 Z RT
V
1
= Z RT
1
8315 J kmolK.
300K
3
P1
=
≅ 25m
5 2
10 N m
2 1,
1 = Z alınarak T R Z V P 2 2
V 2 = Z 2 RT
1,
1.
8315.
300
P2
=
≅ 0,
09
5
300.
10
bulunur. Van der Waals sabitleri;
1 = ifadesinden 2
m
3
P1 = 1 ifadesinden 1
V ;
kmol
kmol
bhava = 0,0394.0,79 + 0,0312.0,21 = 0,0377 m 3 /kmol
ahava = 1,398.0,79 + 1,3375.0,21 = 1,39 bar m 6 /kmol 2 ’dir.
Tüm veriler yerine yazılarak;
W
t
W t
37
V ;
Teknolojik Araştırmalar: MTED 2007 (4) 31-37
5
0,
09 − 0,
0377 0,
0377.
8,
315.
300 ⎛ 25 − 0,
09 ⎞ 2.
1,
39.
10 ⎛ 25 ⎞
= 8,
315.
300ln
−
⎜
−1⎟
+ ⎜ −1⎟
25 − 0,
0377 25 − 0,
09 ⎝ 0,
09 − 0,
041 ⎠ 25 ⎝ 0,
09 ⎠
= −14224
kJ
kmol
Hava 1 kmol alındığın için W t = −14224kJ
elde edilir.
Görüleceği üzere fugasitelerden hesaplanan teknik iş ve kesin formüllerle hesaplanan teknik iş aynı sonucu
vermektedir. Ayrıca şu da ortaya çıkmaktadır; 1 kmol hava içerisinde 300 bar basınçta 14100 kJ’lük enerji
depoludur.
4. Sonuçlar
Gerçek gazların izoterm, izentropik ve izobar durum değişimleri için Van der Waals denklemi ve Termodinamiğin
I. Yasası kullanılarak iş ve teknik bağıntıları elde edilmiştir. Đzobar durum değişimi için nümeriksel olarak türetilen
ve uygulaması çok basit olan iş formülünün örneklendirilmesinde özellikle hidrojen seçilmiştir. Bunun temel
nedeni son yıllarda temiz enerji taşıyıcısı olarak hidrojenin ön plana çıkmasıdır. Ayrıca çeşitli enerji
kaynaklarından (doğal gaz reformundan, kömürün gazlaştırılmasından, güneşten, rüzgardan vb.) üretilen hidrojenin
sabit basınç altında depolanmasında ne kadarlık bir işe ihtiyaç duyulacağı belirtilmiştir. Bununla beraber izoterm
durum değişimi için yapılan örnek bir uygulamada gerçek gaz olarak hava alınmıştır. Türetilen kesin formüllerle
yapılan hesaplama sonucunun; fugasiteler üzerinden yapılan hesaplamayla aynı olduğu görülmüş ve türetilen kesin
formüllerin kullanılabilirliği ortaya konmuştur. Bu çalışmada elde edilen izoterm, izentropik ve izobar durum
değişimleri için iş ve teknik iş formülleri daha sonra yapılacak olan deneysel çalışmalarda rahatlıkla kullanılacaktır.
5. Kaynaklar
1. Bejan, A., “Advanced Engineering Thermodynamics”, John Wiley, 1988.
2. Telli, Z.K., “Termodinamik Problemleri Đle Birlikte”, Palme Yayıncılık, Ankara, 1998.
3. Çengel, Y.A., “Mühendislik Yaklaşımıyla Termodinamik”, Literatür Yayıncılık, Ankara, 1996.
4. Telli, Z.K., “Termodinamik Semineri”, Palme Yayıncılık, Ankara, 1998.
5. Bronstein,I.,N., Semendjajew, K.,A., “Taschenbuhder Mathematik, Verlagtlarri Deutsch, 1968, Zürrch und
Frankfurt, 8 Auflage, p.283-330.
6. Telli, Z.K, “Yakıtlar ve Yanma”, Süleyman Demirel Üniversitesi, Mühendislik Mimarlık Fakültesi Makine
Mühendisliği Bölümü, Isparta, 1996.