FARKLI LAMEL TİPLERİNİN KANATLI BORULU ISITMA … · 2018. 8. 15. · Grafik 1’deki gibi hava...
Transcript of FARKLI LAMEL TİPLERİNİN KANATLI BORULU ISITMA … · 2018. 8. 15. · Grafik 1’deki gibi hava...
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________1
_______________________________________________________________________________________________Sempozyumu Bildirisi
FARKLI LAMEL TİPLERİNİN KANATLI BORULU ISITMA
BATARYASI PERFORMANSI ÜZERİNE ETKİLERİNİN
SAYISAL YÖNTEMLERLE KARŞILAŞTIRILMASI
Yrd. Doç. Dr. Ayhan ONAT Mert Ekin ÖZKAN Görkem ZENGİN ÖZET
Kanatlı borulu ısı değiştiricilerin uygun tasarımı ve seçimi, enerji verimliliği açısından çok önemlidir. Enerji verimli ürünlerin tasarımı ve seçimi, ekonomik ve ekolojik olarak büyük avantajlar sağlamaktadır.
Bu çalışmanın amacı, değişik lamel tiplerine sahip kanatlı borulu ısıtma bataryalarının kapasite ve basınç kaybı gibi performans değerlerinin sayısal yöntemlerle hesaplanarak karşılaştırma yapılmasına bağlı olarak enerji verimli ürün tasarımına destek vermektir.
Yapılan bu çalışmada, aynı şartlarda çalışan aynı geometrik özelliklerdeki ısıtma bataryalarında farklı lamel tiplerinin kullanılmasıyla oluşan performans farklılıkları sayısal yöntemlerle incelenmiştir. Isı değiştirici kapasite hesapları yapılırken ε-ntu yöntemi kullanılmıştır.
Anahtar Kelimeler: Enerji Verimliliği, Kuru Soğutucu, Isıtma Bataryası, Kanatlı Borulu Isı Değiştirici, ε -Ntu Yöntemi
ABSTRACT
Proper design and selection of fin and tube heat exchangers are very important with respect to energy efficiency. Design and selection of energy efficient products provide advantages economically and ecologically.
The purpose of this study is to support design of energy efficient products based on comparing fin and tube heater coils’ performances which has different fin types as calculating performance values like capacity and pressure drop with numerical methods.
In this study, performance differences based on different fin types on heater coils, which has the same geometrical properties and works under the same ambient conditions, has investigated with numerical methods. ε-NTU method were used calculating heat exchanger capacity.
Key Words: Energy Efficiency, Dry Cooler, Heater Coil, Fin And Tube Heat Exchanger, ε -NTU Method
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________2
_______________________________________________________________________________________________Sempozyumu Bildirisi
1. GİRİŞ
Kanatlı borulu ısıtma bataryaları kuru soğutma prensibine göre çalışan endüstriyel ekipmanlardır. Endüstrinin hemen her alanında, iklimlendirme, enerji ve diğer tüm alanlarda gerek atık ısının uzaklaştırılmasında, gerekse mahal iklimlendirme ihtiyacında kanatlı borulu ısı değiştiriciler kullanılabilmektedir.
Değişen çevre koşulları, su kaynaklarının azalması, enerji krizleri ve diğer nedenler işletmelerin soğutma ekipmanlarından ve üretici firmalardan temel beklentilerini düşük enerji tüketimi, düşük kurulum ve işletme maliyetlerine yoğunlaştırmaktadır.
Isıtma bataryalarında ısıyı ortamdan uzaklaştıran temel akışkan havadır. Havayı ısı transfer yüzeyiyle maksimum seviyede temas etmeye zorlamak, basınç kaybı yüksek bir sonuç ortaya çıkarır. Basınç kaybını azaltmak için havanın ortamda rahat bir akış halinde dolaşımı ise ısı transfer kapasitesini olumsuz etkileyecek sonuçlar doğurur. Bu nedenle hava tarafı dizaynının da olabildiğince verimli olması genel olarak sistemi verimli bir hale getirecektir.
Sektörde çeşit çeşit lamel yapıları bulunmaktadır. Başlıca; dalgalı, düz ve panjurlu lamel tipleri kullanılmaktadır. Bu lamel tiplerinin birbirlerine göre avantajları ve dezavantajları bulunmaktadır.
Bu çalışmada çeşitli notasyon ve lamel yapılarındaki kanatlı borulu ısıtma bataryalarının, farklı dizayn koşullarındaki performanslarının enerji verimliliği yönünden karşılaştırılması literatürde test edilerek doğrulanan sayısal çalışmalar yardımıyla yapılmıştır. Farklı koşullardaki ısıl performansları incelenirken ε-NTU metodu kullanılmıştır. 2. ÇALIŞMA PRENSİBİ
Kanatlı borulu ısı değiştiricilerinin çalışma prensibi oldukça basit olmakla beraber, ısı geçişini etkileyen parametrelerin çokluğu sebebiyle performans analizini yapmak bir o kadar zordur. Isı geçişi birbirine çapraz olarak akan, yani aralarındaki açı 90° olan akışkanlar arasında olur. İç akışkan, ısı değiştiricisi içinde mümkün olduğu kadar gezinir.
Şekil 1. Kanatlı Borulu Isı Değiştiricisinin Sembolik Gösterimi
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________3
_______________________________________________________________________________________________Sempozyumu Bildirisi
Şekil 1’de standart bir kanatlı borulu ısı değiştiricinin çalışma mantığı gösterilmektedir. Soğuk akışkanı hava, sıcak akışkanı da su olarak ele aldığımızda Grafik 1’deki gibi bir sıcaklık – konum grafiği elde ederiz.
Grafik 1. Kuru Soğutucu Sıcaklık Analiz Grafiği
Genelde ısıtma bataryalarında, başlangıçta su ortamdaki sıcak akışkan iken, hava ise soğuk akışkandır. Grafik 1’deki gibi hava suyun giriş sıcaklığına yaklaşırken, su da havanın giriş sıcaklığına yaklaşmıştır. Bu grafik iki akışkanın da tek fazda olduğu, herhangi bir faz değişiminin gerçekleşmediği şartlar için geçerlidir.
Grafik 2. Isıtma Bataryasında Havanın Psikrometrik İncelemesi
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________4
_______________________________________________________________________________________________Sempozyumu Bildirisi
Grafik 2 incelendiğinde standart bir ısıtma bataryasında giren havanın psikrometrik diyagram üzerinde izlediği yol gösterilmiştir. Hava içeriğindeki su miktarı değişmediği için entalpi artışı kuru termometre sıcaklığına bağlıdır. Bu nedenle yüksek ısı transfer kapasitelerinin ısıtma bataryaları ile sağlanması için yüksek hava debilerine ihtiyaç duyulmaktadır. 3. HESAPLAMA ALGORİTMASI
Enerji verimliliği tüm sistemlerde temel olarak elde edilen gücün harcanan güce oranı olarak tanımlanır. Kanatlı borulu ısı değiştiricilerinde elde edilen güç olarak soğutma kapasitesi, harcanan güç olarak ise pompa ve fanların harcadığı güç alınabilir. Bu yüzden enerji verimliliğini doğrudan etkileyen ilk parametre ısı değiştiricinin soğutma kapasitesidir.
Endüstride kanatlı borulu ısıtma bataryası seçimi yapılırken genelde su debisi veya su çıkış sıcaklıkları bilinmemektedir. Bu yüzden klasik yöntemlerle hesaplama yapmak mümkün olmamaktadır.
Bu çalışmada ε -NTU hesaplama yöntemi kullanılarak geliştirilen iteratif algoritmayla ısıtma bataryası performansları sayısal olarak hesaplanmıştır.
Bu algoritma, temel olarak 6 aşamadan meydana gelmektedir. • akışkanların termofiziksel özelliklerinin hesaplanması • geometrik hesaplamalar • su ısıl hesaplamalarının yapılması • hava ısıl hesaplamalarının yapılması • ısı değiştirici ısıl hesaplamalarının yapılması • ısı değiştirici etkenliğinin hesaplanması İlk aşamada akışkanların sıcaklık değerlerine bağlı olarak, özgül ısıları, viskozite değerleri, yoğunlukları ve termal iletkenlikleri hesaplanır.
Ardından akışkan hızlarının ve ısı transfer alanlarının belirlenmesi için temel geometrik hesaplamalar yapılır.
Boru içindeki akışkan hızının hesaplanması için debi ve alan parametrelerine ihtiyaç duyulmaktadır. Boru içindeki suyun debisi bilinmediği için burada iteratif bir yöntem izlenerek afaki bir ε değeri atanır.
ε değeri, ısı değiştiricinin genel etkenliği olarak tanımlanıp
ε = Q / Qmax (1) olarak formülize edilebilir. Buradaki Qmax teorik kapasite olarak tanımlanabilir ve boru içindeki akışkanın boru dışındaki akışkanın giriş sıcaklığına eşitlendiği durumdaki ısı değiştiricinin kapasitesidir.
Qmax = cmin x ΔTmax (W) (2) Afaki ε değeri ve dolasıyla afaki kapasite parametreleri bilindiğinden, boru içi akışkanın giriş ve çıkış sıcaklıkları yardımıyla akışkan debisi belirlenir.
𝑚𝑠𝑢 =𝑄
𝑐𝑚𝑎𝑥 (kg/s) (3)
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________5
_______________________________________________________________________________________________Sempozyumu Bildirisi
Akışkan debisine bağlı olarak hesaplanan hızdan Reynolds sayısı hesaplanarak akış karakteristiği belirlendikten sonra, Nusselt sayısı ve buna bağlı olarak boru içi akışkanının ısı transfer katsayısı hesaplanır.
ν𝑠𝑢 = 𝑚𝑠𝑢
𝐴𝑠𝑢𝑥 𝜌𝑠𝑢 (m/s) (4)
𝑅𝑒𝑠𝑢 =𝜌𝑠𝑢𝑥 ν𝑠𝑢 𝑥 𝐷𝑖ç
μ𝑠𝑢 (5)
Bu çalışmada incelenen bataryalarda türbülanslı akış olduğu için (𝑅𝑒𝑠𝑢 > 2300) türbülanslı akış için Gnielinski’nin korelasyonu seçilmiştir.
𝑁𝑢𝑠𝑢 = ((
𝑓8⁄ )𝑥(𝑅𝑒𝑠𝑢−1000)𝑥 𝑃𝑟𝑠𝑢
1+12.7𝑥√(𝑓
8⁄ )𝑥(𝑃𝑟2/3−1)
) (6)
Burada;
𝑓 = (0,79 𝑥 ln 𝑅𝑒𝑠𝑢 − 1,64)−2 (7)
olarak tanımlanabilir. Su için ısı transfer katsayısı Nu sayısından hareketle hesaplanabilir:
ℎ𝑆𝑢 = 𝑁𝑢𝑠𝑢 𝑥 𝑘𝑠𝑢
𝑑𝑖ç (W/m2K) (8)
Hava debisi doğrudan belirli olduğu için direkt olarak Reynolds sayısı hesaplanabilir.
𝑅𝑒ℎ𝑎𝑣𝑎 =𝜌ℎ𝑎𝑣𝑎𝑥 νℎ𝑎𝑣𝑎 𝑥 𝐷ℎ
μℎ𝑎𝑣𝑎 (9)
Ancak ısı transfer katsayının hesaplanması için gerekli Colburn sayısının hesaplanması için bir çok deneysel çalışma yapılarak sayısal korelasyonlar ortaya konulmuştur. Tüm bu korelasyonlar lamel geomerisine bağlı olarak değişmekledir.
Düz lamel geometrisi Şekil 2.’de gösterilmiş olup bu geometri için korelasyonlar şu şekildedir:
Şekil 2. Düz Lamel Geometrisi
N = 1 için;
j = 0,173 x ReDc0,346x (
Pt
Pl)
P1
x (Dc
ξk)1,161 x (
Dh
ξk)1,035 x (
ξk
Pt)P2 (10)
N > 1 için;
http://en.wikipedia.org/wiki/Mu_(letter)http://en.wikipedia.org/wiki/Mu_(letter)
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________6
_______________________________________________________________________________________________Sempozyumu Bildirisi
j = 0,078 x ReDcP3 x NP4 x (
ξk
Dh)P5 x (
Pt
Fp)1,026 (11)
Bu denklemlerde kullanılan P1, P2, P3, P4 ve P5katsayıları şu şekilde hesaplanır:
P1 = − 0,22 x ln(ReDc) + 1,88 (12)
P2 = 0,106 x ln(ReDc) (13)
P3 = 0,16 x ln[N x (Fp
Dh)0,42] − 0,349 (14)
P4 =− 0,094 X (
plDh
)1,38
ln(ReDc)− 1,405 (15)
P5 = 1,263 x ln (ReDc
N) − 5,97 (16)
Bu denklemdeki Dh şu şekilde hesaplanır:
Dh =4 x Amin x L
Ao (m) (17)
Dalgalı lamel geometrisi ise Şekil 3.’te gösterilmiş olup bu geometri için korelasyonlar şu şekildedir:
Şekil 3. Dalgalı Lamel Geometrisi
Ddış ≥ 8,4 mm ise;
j = 1,79097 x ReDc
0,1707−1,374 x (Plξf
)−0,493
x (ξkDc
)−0,886
x N−0,43x (Pdxf
)−0,0296
x (Pl
ξf)−0,456 x N−0,27 x (
ξk
Dc)−1,34x (
Pd
Xf) (18)
Ddış < 8,4 mm ise;
j = 0,324 x ReDcJ1
x (ξk
Pl)J2 x (tan θ)J3 x (
Pl
Pt)J4 x N0,428 (19)
Bu denklemlerde kullanılan J1, J2, J3, J4 katsayıları şu şekilde hesaplanır:
𝐽1 = −0,229 + 0,115 ( 𝐹𝑝
𝐷𝑐)
0,6
𝑥 ( 𝑃𝑙
𝐷ℎ)
0,54
𝑥 𝑁−0,284𝑥 ln (0,5 tan 𝜃) (20)
𝐽2 = −0,251 + 0,232 𝑥 𝑁1,37
ln(𝑅𝑒𝑑𝑐)−2,303 (21)
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________7
_______________________________________________________________________________________________Sempozyumu Bildirisi
𝐽3 = −0,439 + (𝐹𝑝
𝐷ℎ)
0,09
𝑥 (𝑃𝑙
𝑃𝑡)
−1,75
𝑥 𝑁−0,93 (22)
𝐽4 = 0,502 𝑥 ln(𝑅𝑒𝐷𝑐) − 2,54 (23)
Panjurlu lamel geometrisi de Şekil 4.’te gösterilmiş olup bu geometri için korelasyonlar şu şekildedir:
Şekil 4. Panjurlu Lamel Geometrisi
Redc < 1000 ise;
j = 0,0501 x ReDc0,094 x (
Pl
Pt)1,44 x (
ξk
Dc)0,322x N−0,809 x (
Lp
ξk )−0,177 x (
Lh
Lp)−0,254 (24)
Redc ≥ 1000 ise;
j = 0,962 x ReDcL1 x (
Pl
Pt)−1,51 x (
ξk
Dc)−0,107 x N−1,06 x (
Lh
Lp)−0,433 (25)
Elde edilen Colburn sayısından hareketle havanın ısı taşınım katsayısı belirlenir.
ℎℎ𝑎𝑣𝑎 = (𝑗ℎ𝑎𝑣𝑎∗𝑚ℎ𝑎𝑣𝑎∗𝑐𝑝ℎ𝑎𝑣𝑎
𝑝𝑟ℎ𝑎𝑣𝑎0.66 𝑥 𝐴𝑜
) (W/m2K) (26)
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________8
_______________________________________________________________________________________________Sempozyumu Bildirisi
Kanatlı borulu ısı değiştiricilerinde kanatların amacı ısı transfer yüzeyini arttırmaktır. Fakat bu ısı transfer alanının tamamının üstünden akışkan geçmez. Çünkü kanat alanının bir kısmı borular tarafından kaplıdır ve bir kısmı akışkanla temas halinde değildir. Dolayısıyla ısı değiştiricisinin kapasitesi hesaplanırken bu transfer alanının ne kadarının verimli kullanıldığını hesaplamak gerekir. Kanat verimliliğinin hesaplanması
için μ ve l katsayılarını hesaplanması gerekmektedir. Bu katsayıları aşağıdaki bağıntılar yardımıyla hesaplanır.
𝜇 = √2 𝑥 ℎ𝑑𝚤ş
𝑘𝑑𝚤ş 𝑥 𝛿𝑑𝚤ş (27)
𝑙 =(𝑃𝑡−𝐷𝑐)
2 (28)
Ψ𝑘𝑎𝑛𝑎𝑡 = tanh(𝑚 𝑥 𝑙)
𝑚 𝑥 𝑙 (29)
Bir kanatın verimliliği hesaplandıktan sonra tüm kanatların toplam verimliliğine ulaşmak için şu bağıntı kullanılır.
Ψ𝑡𝑜𝑝 = 𝐴𝑘𝑎𝑛𝑎𝑡 𝑥 Ψ𝑘𝑎𝑛𝑎𝑡 𝑥+ 𝐴𝑡ü𝑝
𝐴𝑡𝑜𝑝𝑙𝑎𝑚 (30)
Su ve hava tarafının ısı taşınım katsayılarından hareketle bataryanın toplam ısı transfer katsayısı şu şekilde hesaplanır.
1
𝑈 𝑥 𝐴=
1
Ψ𝑡𝑜𝑝𝑥 ℎ𝑑𝚤ş 𝑥 𝐴𝑑𝚤ş+
𝜉𝑏𝑜𝑟𝑢
𝑘𝑏𝑜𝑟𝑢 𝑥 𝐴𝑑𝚤ş+
1
ℎ𝑖ç 𝑥 𝐴𝑖ç (K/W) (31)
Hesaplanan ısı transfer katsayılarına bağlı olarak toplam ısı transfer katsayısı hesaplanır. Buna bağlı olarak:
𝑁𝑇𝑈 = 𝑈 𝑥 𝐴
𝑐𝑚𝑖𝑛 (32)
formülüyle ısı transfer ünitesi sayısının (NTU) değeri hesaplanır. Isı değiştiricilerde etkenlik ifadesi akışkanların ısı sığalarının oranları (cmin/cmax) ve ısı transfer ünitesi sayısının bir fonksiyonudur. Etkenlik ifadesi akış biçimine ve kanat formuna göre değişiklik göstermektedir. Geliştirilen hesap yöntemi boru-boru olduğu için her bir boru için ayrı inceleme yapılmıştır ve çapraz akış kabulü yapılmıştır.
ε = f(NTU, cr) (33) Burada hesaplanan ε değeri ilk başta afaki olarak atanan ε değerine göre hesaplandığı için bu iki değerin birbirine yakınsama durumuna bakılır. Yakınsama olmaması durumunda hesaplanan ε değeriyle aynı işlemler yakınsama oluncaya kadar tekrarlanır. Yakınsama sağlandığında hesaplanan ε değerine göre gerçek kapasite hesaplanmış olur.
Enerji verimliliğini doğrudan etkileyen başka bir parametre de basınç kayıplarıdır. Isı değiştiricilerde oluşan basınç kayıpları doğrudan fan ve pompa seçimlerini etkileyecek ve buna bağlı olarak motor seçimlerini de etkileyecektir. Hava tarafı basınç kaybı için ise literatürde deneysel çalışmalara bağlı olarak çeşitli ampirik korelasyonlar mevcuttur. Buradaki korelasyonların lamel geometrisine göre değişmekle birlikte bu çalışmada aşağıdaki korelasyonlar kullanılmıştır: Düz lamel için;
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________9
_______________________________________________________________________________________________Sempozyumu Bildirisi
f = 0,0146 x ReDcP6 x (
Pt
Pl)
1,959
x (ξk
Dc)P7 x N0,021 (34)
Bu denklemlerde kullanılan P6 ve P7 katsayıları şu şekilde hesaplanır:
P6 = − 0,0535 + 0,01166
ln (Ptpl
)+ 0,123 (
Fp
Dc) (35)
P7 = 2,319 − 19,59
ln (ReDc) (36)
Dalgalı lamel için;
Ddış ≥ 8,4 mm ise;
f = 0,05273 x ReDcf1 x (
Pd
Xf)f2 x (
ξk
Pt)f3 x (ln
(Adış
Aiç )−2,276 x (Dh
Dc)0,1325x N0,02305 (37)
Ddış < 8,4 mm ise;
f = 0,01915 x ReDcF1 x (tanθ)F2 x (
ξk
Pl)F3 x (ln (
Adış
Aiç))−5,35x (
Dh
Dc)1,3796x N−0,0916 (36)
Bu denklemlerde kullanılan F1, F2 ve F3 katsayıları şu şekilde hesaplanır:
F1 = 0,4604 – 0,01336 x (ξ𝑘
𝑃𝑡)0,58 𝑥 ln (
𝐴𝑑𝚤ş
𝐴𝑖ç) 𝑥 (𝑡𝑎𝑛𝜃)−1,5 (37)
F2 = 3,247 x (ξ𝑘
𝑃𝑡)1,4 𝑥 ln(
𝐴𝑑𝚤ş
𝐴𝑖ç) (38)
F3 = −20,113
ln (𝑅𝑒𝐷𝑐) (39)
Panjurlu lamel için;
f = 0,72 x ReDcZ1 x (
ξk
Dc)Z2 x N0,322 x (
Pl
Pt)−2,494 x (
Lh
Lp)1,34 x (
Lp
ξk)1,69 (40)
Bu denklemlerde kullanılan L1, Z1 ve Z2 katsayıları şu şekilde hesaplanır:
L1 = −0,501 x (Pl
Pt)−0,424 x N−0,253 x (
Lh
Lp)0,127 (41)
Z1 = −0,392 x (pL
p′)−0,726 x N0,0782 x (
Lh
Lp)0,352 x (
Lp
ξk)0,63 (42)
Z2 = −3,3 + ( 3,72 x (Pl
Pt)
−1
) − (1,29 x (Pl
Pt)
−2
) (43)
Elde edilen f sürtünme faktörüne göre hava tarafı basınç kaybı Darcy-Weishbach eşitliğine göre
hesaplanır.
𝛥𝑃ℎ𝑎𝑣𝑎 = 𝑓 𝑥 𝐿
𝐷 𝑥
𝜌𝑉2
2 (Pa) (44)
Bu çalışmada kullanılan korelasyonlara göre hesaplanmış olan Colburn sayısı ve sürtünme faktörlerindeki sonuçlarla, test sonuçları arasındaki sapma aşağıdaki gibidir:
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________10
_______________________________________________________________________________________________Sempozyumu Bildirisi
Tablo 1. Düz Lamel İçin Colburn Sayısı ve Sürtünme Faktörü Sapma Değerleri
Sapma Colburn
Sayısı Sürtünme Faktörü
±10% 76.3% 67.4%
±15% 87.4% 84.2%
±20% 92.8% 92.5%
Ortalama Sapma 7.64% 8.97% Tablo 2. Dalgalı Lamel İçin Colburn Sayısı ve Sürtünme Faktörü Sapma Değerleri
Sapma Colburn
Sayısı Sürtünme Faktörü
±10% 79.6% 84,9%
±15% 95,1% 97,3%
±20% 97,4% 99,6%
Ortalama Sapma 6,44% 5,01% Tablo 3. Panjurlu Lamel İçin Colburn Sayısı ve Sürtünme Faktörü Sapma Değerleri
Sapma Colburn
Sayısı Sürtünme Faktörü
±10% 79.7% 59%
±15% 93,2% 75,3%
±20% 96,4% 86,9%
Ortalama Sapma 6,83% 10,1% 4. LAMEL TİPLERİNİN ISIL PERFORMANS VE ENERJİ TÜKETİMİ AÇISINDAN KARŞILAŞTIRMASI
Bu bölümde aynı dizayn koşullarında ve soğutma kapasitelerinde farklı lamel tiplerinin performans karşılaştırmaları yapılmıştır. Bu performans karşılaştırmalarında sadece hava tarafı dikkate alınmıştır. Bu nedenle belirlenen akışkanın ısı değiştirici içerisindeki dolanımı için gerekli olan pompa ve pompanın enerji tüketimi dikkate alınmamıştır.
Karşılaştırmaların performans analizleri yapılırken sayısal metotlar kullanılmıştır. Analizler için kanatlı borulu ısı değiştiriciler eş konstrüktif notasyonlarda kullanılmıştır. Bunun nedeni şartları mümkün olduğunca sadeleştirip, farkları azaltıp analizi daha derinlemesine, olabildiğince lamel yapısındaki farklılıkların performansa etkisinin incelenebilmesinden kaynaklıdır.
Tablo 4.’te karşılaştırma yapılacak ısı değiştiricilerin konstrüktif özellikleri verilmiştir.
Tablo 4. Karşılaştırmak İçin Kullanılan Isı Değiştiricinin Konstrüktif Özellikleri
Boru Sayısı 56 Adet
Sıra Sayısı 4 Adet
Devre Sayısı 112 Adet
Geçiş Sayısı 2 Adet
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________11
_______________________________________________________________________________________________Sempozyumu Bildirisi
Borular Arası Mesafe 40 mm
Sıralar Arası Mesafe 35 mm
Boru Dış Çapı 15,8 mm
Boru Et Kalınlığı 0,4 mm
Lamel Dizili Boru Uzunluğu 10.000 mm
Lameller Arası Mesafe (Hatve) 2,1 mm
Lamel Kalınlığı 0,1 mm
Boru Malzemesi Bakır -
Lamel Malzemesi Alüminyum -
Koşul 1
Su giriş sıcaklığı 45 °C,
Su çıkış sıcaklığı 40 °C,
Hava sıcaklığı 35 °C,
Giriş Bağıl Nem %50,
Su debisi 75 m3/h, Tablo 5. Koşul 1 İçin Fan Motoru Enerji Tüketimleri
Lamel
Tipi
Gerekli Hava
Debisi
Hava Direnci Hava Tarafı Isı
Transfer Katsayısı
Hava Çıkış
Sıcaklığı
Fan Motoru
Enerji Tüketimi
m3/s Pa W/(m2K) °C kW
Düz 57,6 67 61 41,4 7,2
Dalgalı 53,7 75 67 41,8 7,5
Panjurlu 50,4 78 76 42,3 7,3
Koşul 2
Su giriş sıcaklığı 45 °C,
Su çıkış sıcaklığı 40 °C,
Hava sıcaklığı 32,5 °C,
Bağıl Nem %50,
Su debisi 75 m3/h, Tablo 6. Koşul 2 İçin Fan Motoru Enerji Tüketimleri
Lamel
Tipi
Gerekli Hava
Debisi
Hava Direnci Hava Tarafı Isı
Transfer Katsayısı
Hava Çıkış
Sıcaklığı
Fan Motoru
Enerji Tüketimi
m3/s Pa W/(m2K) °C kW
Düz 39,0 35 48 41,9 2,5
Dalgalı 37,1 41 54 42,3 2,8
Panjurlu 35,3 43 61 42,8 2,8
Koşul 3
Su giriş sıcaklığı 45 °C,
Su çıkış sıcaklığı 40 °C,
Hava sıcaklığı 30 °C,
Bağıl Nem %50,
Su debisi 75 m3/h,
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________12
_______________________________________________________________________________________________Sempozyumu Bildirisi
Tablo 7. Koşul 3 İçin Fan Motoru Enerji Tüketimleri
Lamel
Tipi
Gerekli Hava
Debisi
Hava Direnci Hava Tarafı Isı
Transfer Katsayısı
Hava Çıkış
Sıcaklığı
Fan Motoru
Enerji Tüketimi
m3/s Pa W/(m2K) °C kW
Düz 29,7 22 40 42,2 1,2
Dalgalı 28,4 26 45 42,7 1,3
Panjurlu 27,4 28 52 43,2 1,4
3,5
Grafik 3. Koşul 1, Koşul 2 ve Koşul 3 İçin Lamel Tiplerinin Enerji Tüketimi Karşılaştırması
Koşul 4
Su giriş sıcaklığı 45 °C,
Su çıkış sıcaklığı 40 °C,
Hava sıcaklığı 35 °C,
Bağıl Nem %50,
Su debisi 100 m3/h, Tablo 8. Koşul 4 İçin Fan Motoru Enerji Tüketimleri
Lamel
Tipi
Gerekli Hava
Debisi
Hava
Direnci
Hava Tarafı Isı
Transfer Katsayısı
Hava Çıkış
Sıcaklığı
Fan Motoru
Enerji Tüketimi
m3/s Pa W/(m2K) °C kW
Düz 87,3 136 79 40,6 22,1
Dalgalı 80,6 149 87 41,1 22,4
Panjurlu 74,4 151 97 41,6 20,9
Koşul 5
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________13
_______________________________________________________________________________________________Sempozyumu Bildirisi
Su giriş sıcaklığı 45 °C,
Su çıkış sıcaklığı 40 °C,
Hava sıcaklığı 32,5 °C,
Bağıl Nem %50,
Su debisi 100 m3/h, Tablo 9. Koşul 5 İçin Fan Motoru Enerji Tüketimleri
Lamel
Tipi
Gerekli Hava
Debisi
Hava
Direnci
Hava Tarafı Isı
Transfer Katsayısı
Hava Çıkış
Sıcaklığı
Fan Motoru
Enerji Tüketimi
m3/s Pa W/(m2K) °C kW
Düz 56,9 67 61 41,1 7,1
Dalgalı 53,5 75 67 41,6 7,4
Panjurlu 50,4 79 76 42,2 7,4
Koşul 6
Su giriş sıcaklığı 45 °C,
Su çıkış sıcaklığı 40 °C,
Hava sıcaklığı 30 °C,
Bağıl Nem %50,
Su debisi 100 m3/h, Tablo 10. Koşul 6 İçin Enerji Tüketimleri
Lamel
Tipi
Gerekli Hava
Debisi
Hava
Direnci
Hava Tarafı Isı
Transfer Katsayısı
Hava Çıkış
Sıcaklığı
Fan Motoru
Enerji Tüketimi
m3/s Pa W/(m2K) °C kW
Düz 42,5 41 51 41,4 3,2
Dalgalı 40,3 47 56 42 3,5
Panjurlu 38,5 50 64 42,6 3,6
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________14
_______________________________________________________________________________________________Sempozyumu Bildirisi
Grafik 4. Koşul 4, Koşul 5 ve Koşul 6 İçin Lamel Tiplerinin Enerji Tüketimi Karşılaştırması
5. SONUÇ
Kanatlı borulu ısı değiştiricilerin ısıl performansı ve enerji tüketimleri direkt olarak rejim ve lamel yapısı ile
değişken olduğu 6 koşulda da görülmektedir. Genel olarak su çıkış sıcaklığının hava giriş sıcaklığına
yakın olduğu dizayn şartlarında daha fazla hava debisi gerektiği için enerji tüketimleri de buna bağlı
olarak yüksek değerlerdedir.
Bu bildiride araştırılan ve bilgi verilmek istenen ana konu, akademisyenler ve araştırmacılar tarafından
termal analizleri yapılmış başlıca lamel yapılarının aynı notasyon ve dizayn koşullarındaki enerji
tüketimlerini karşılaştırmaktır.
Analizlerde de görüldüğü üzere düz lamel, termal performansı sağlamak için en fazla hava debisine
ihtiyaç duyan lamel yapısıdır. Bunun avantajı iki akışkan arasındaki sıcaklık farkının fazla olduğu
durumlarda düşük basınç kaybı değerlerinde çalışması avantajı ile birleşerek düşük enerji tüketimi
değerlerine ihtiyaç duymasını sağlamaktadır. Lakin diğer incelenen lamel yapılarına kıyasla da en düşük
ısı transfer katsayısına sahip olmalarının sonucu olarak, yüksek termal performans gerektiren zor
koşullarda, yani iki akışkan arasındaki sıcaklık farkının düşük olduğu durumlarda, verimsiz bir lamel tipi
olduğu görülmektedir.
Dalgalı lamel yapısı gerek hava direnci olarak gerekse ısı transfer katsayısı bakımından ortalama
değerlere sahip ideal yapıdadır. Aynı termal dizaynda panjurlu ve düz lamele kıyasla daha ortalama hava
debisine ihtiyaç duyması sonuç olarak daha ortalama tüketim değerlerinde çalıştığını göstermiştir.
Panjurlu lamel diğer lamel yapılarına kıyasla daha yüksek ısı transfer katsayılarına ulaşmaktadır. Bunun
sonucunda aynı termal dizayn ve konstrüktif notasyonda daha az hava debisine ihtiyaç duymaktadır.
Lakin havayı diğer lamel yapılarına göre daha fazla zorladıkları için ünite içerisindeki hava direnci (basınç
kaybı) diğer lamel yapılarına kıyasla yüksektir. Bunların en net avantajı zor dizayn şartlarında düşük hava
debisi ihtiyacı ile görülmektedir.
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________15
_______________________________________________________________________________________________Sempozyumu Bildirisi
Bu çalışmada akademisyenler ve araştırmacılar tarafından performans analizi yapılan belirli
notasyonlardaki lamel yapıları enerji tüketimi ve dolayısıyla enerji verimliliği açısından analiz edilmiştir.
Bilindiği üzere her kanatlı borulu ısı eşanjörünün gerek üretim tekniği, gerek kullanılan malzeme gerekse
lamel yapısı ile değişken kendine özgü performans karakteristiği vardır. Buradaki temel husus şudur ki
elde edilen sonuçlar referanslardaki notasyonlar için geçerlidir. Örnek ile açıklamak gerekirse, dalgalı
lamel yapısını farklı notasyonlarda farklı imalat teknikleri ve malzeme kalitesiyle daha verimli bir lamel
haline getirmek mümkündür.
Bu çalışma sonucunda ideal lamel yapısının seçiminde termal dizayn kriterlerinin önemi irdelenmiştir.
Farklı dizayn koşulları için lamel tiplerinin, gerekli hava debisi, hava direnci, ısı transfer katsayısı ve temel
kriter olarak fan motorunun toplam elektrik enerjisi tüketimleri analiz edilmiştir.
6. TERMİNOLOJİ
SEMBOLLER ANLAMLARI A Alan c Isıl Sığa D Çap f Sürtünme Faktörü h Isı Taşınım Katsayısı j Colburn Sayısı k Isı İletim Katsayısı L Uzunluk m Kütlesel Debi N Bir Sıradaki Boru Sayısı
NTU Isı Transfer Ünite Sayısı P Basınç Pr Prandtl Sayısı Pl Borular Arası Boyuna Mesafe Pt Borular Arası Enine Mesafe ε Isı Değiştirici Etkenliği Q Isıtma Bataryası Isı Transfer Kapasitesi T Sıcaklık U Isı Transfer Katsayısı V Hız
Re Reynold Sayısı Nu Nusselt Sayısı 𝜌 Yoğunluk 𝜇 Dinamik Viskozite ξk Hatve Fp θ Dalga Açısı Ψ Verimlilik
KAYNAKLAR
[1] Holman J. P., Thermodynamics, McGraw-Hill, New York. Streeter, Victor, L., Fluid Mechanics, 5th
Edition, McGraw-Hill, New York, ISBN
[2] Xiaokui Maa, Guoliang Dinga, Yuanming Zhanga, Kaijian Wang, (8 November 2006), Airside heat
transfer and friction characteristics for enhanced fin-and-tube heat exchanger with hydrophilic coating
under wet conditions, Science Direct.
-
13. ULUSAL TESİSAT MÜHENDİSLİĞİ KONGRESİ – 19-22 NİSAN 2017/İZMİR_______________________________________16
_______________________________________________________________________________________________Sempozyumu Bildirisi
[3] C.C. Wang, Y.J. Du, Y.J. Chang, W.H. Tao, Airside performance of herringbone fin-and-tube heat
exchangers in wet conditions, Canadian Journal of Chemical Engineering.
[4] C.C. Wang, Y.T. Lin, C.J. Lee, Heat and momentum transfer for compact louvered fin-and-tube heat
exchangers in wet conditions, International Journal of Heat and Mass Transfer
[5] McQuiston, F.C., (1978), Correlation of heat, mass and momentum transport coefficients for plate-fin-
tube heat transfer surfaces with staggered tubes. ASHRAE Transactions, Vol. 84, No. 1, pp. 294-309.
[6] Gesellschaft Verfahrenstechnik und Chemieingenieurwesen, (2010), VDI Heat Atlas, 2. baskı,
Springer-Verlag Berlin Heidelberg, Düsseldorf.
[7] Shah R.K. ve Sekulic D.P., (2003), Fundementals of Heat Exchanger Design, John Wiley & Sons, Inc.,
New Jersey.
[8] Kakaç S., Bergles A.E., Mayinger F., Yüncü H., (1999). Heat Transfer Enhancement of Heat
Exchangers, Series E: Applied Sciences – Vol. 355, Netherlands.
[9] Cengel Y. A., Ghajar A., (2010) , Heat and Mass Transfer: Fundamentals and Applications, 4. Baskı.,
McGraw-Hill Science/Engineering/Math., New York.
[10] İsa K., Onat A., (2012) , İklimlendirme ve Soğutma Sistemlerinde Enerji Verimliliği, 1. Basım., Doğa
Yayıncılık., İSTANBUL.
ÖZGEÇMİŞ
Ayhan ONAT
Ayhan ONAT, 1989 yılında Marmara Üniversitesi’nin Enerji Eğitimi bölümünden mezun olmuştur. 1994
yılında Yüksek Lisansını ve 2002 yılında Doktora öğrenimini aynı bölümde tamamlamıştır.
Akademisyenliğinin yanında birçok sektörel projede yer almış olan Ayhan ONAT halen Marmara
Üniversitesi Makine Mühendisliği bölümünde Yrd. Doç. Dr. olarak görev yapmaktadır
Mert Ekin ÖZKAN
Mert Ekin ÖZKAN, 2016 yılında Marmara Üniversitesi’nin Makine Mühendisliği bölümünden mezun
olmuştur. Mert Ekin ÖZKAN, AIRONN Havalandırma ve Klima Sistemleri Sanayi ve Dış Ticaret Ltd.
Şti.’nde Proje ve Ar-Ge Mühendisi olarak çalışmaktadır
Görkem ZENGİN
Görkem ZENGİN, 2015 yılında Marmara Üniversitesi’nin Makine Mühendisliği bölümünden mezun
olmuştur. Görkem ZENGİN, Cenk Endüstri Tesisleri İmalat Ve Taahhüt A.Ş. bünyesindeki Kapalı ve Kuru
Sistemler Departmanı’nda Makine Mühendisi olarak çalışmaktadır