Akademik araştırmalar, anlaşılır dil

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

27 Eylül 2026, Pazar
VERİANLABağımsız bilim yayıncılığı
Menüyü aç veya kapat
...
Home / Uygulamalı Bilimler / Enerji Araştırmaları / Yüksek Sıcaklıklı Dökme Demir Isıl Depolama Kanallarında Kritik Altı CO₂'nin Zorlanmış Taşınımla Isı Transferinin Deneysel İncelenmesi
Enerji Araştırmaları

Yüksek Sıcaklıklı Dökme Demir Isıl Depolama Kanallarında Kritik Altı CO₂'nin Zorlanmış Taşınımla Isı Transferinin Deneysel İncelenmesi

Araştırma, yüksek sıcaklıklı katı duyulur ısı depolama sistemlerinde depolanmış ısının kritik altı karbondioksit (CO₂) kullanılarak ne kadar etkin çekilebildiğini deneysel olarak incelemektedir.

14/08/2026  Veri Anla 49 görüntüleme
Yüksek Sıcaklıklı Dökme Demir Isıl Depolama Kanallarında Kritik Altı CO₂'nin Zorlanmış Taşınımla Isı Transferinin Deneysel İncelenmesi

Araştırma, yüksek sıcaklıklı katı duyulur ısı depolama sistemlerinde depolanmış ısının kritik altı karbondioksit (CO₂) kullanılarak ne kadar etkin çekilebildiğini deneysel olarak incelemektedir. Araştırmacılar bu amaçla elektrikle ısıtılan dökme demir bir ısıl depolama ortamı, CO₂ dolaşım hattı, değişken frekanslı pompa, doymuş evaporatör, debi ölçer ve sıcaklık-basınç sensörlerinden oluşan kapalı çevrim bir test sistemi kurmuştur. Deneyler 0,2 MPa basınçta ve yüksek Reynolds sayılarında gerçekleştirilmiştir.

Deney sistemindeki CO₂, dökme demir depolama bloğunun içindeki dairesel kanallardan geçerken katı ortamdan ısı almakta ve sıcaklığı akış doğrultusunda yükselmektedir. Isınmış CO₂ daha sonra doymuş evaporatöre ulaşarak taşıdığı enerjiyi besleme suyuna aktarmaktadır. CO₂ çevrim boyunca tek fazlı gaz durumunda kalmakta; faz değişimi yalnız evaporatörün su tarafında gerçekleşmektedir. Bu nedenle araştırma, CO₂'nin dökme demir yüzeyle gerçekleştirdiği zorlanmış taşınımla ısı transferini izole biçimde değerlendirmektedir.

Deneylerde Reynolds sayısı yaklaşık 7,5×104'ten 1,35×105'e yükseldiğinde deneysel Nusselt sayısı yaklaşık 74'ten 97'ye çıkmıştır. Bu sonuç, incelenen çalışma aralığında akış şiddeti ve türbülans arttıkça CO₂ ile sıcak dökme demir yüzey arasındaki taşınımla ısı transferinin güçlendiğini göstermektedir.

Araştırmacılar ölçülen verileri Dittus–Boelter, Gnielinski, Colburn ve Zhou korelasyonlarıyla karşılaştırmış; bu klasik bağıntıların araştırmadaki yüksek sıcaklıklı dökme demir–kritik altı CO₂ sistemi için deneysel Nusselt değerlerini yeterli doğrulukta temsil etmediğini göstermiştir. Bunun üzerine deney verilerinden dört yeni ampirik korelasyon türetilmiş ve mühendislik kullanımı için \(Nu=0.202Re^{0.465}Pr^{-2.19}\) bağıntısı önerilmiştir.

Önerilen korelasyon için belirleme katsayısı R²=0,988, ortalama mutlak yüzde hata (MAPE) %0,715 ve maksimum bağıl sapma %2,404 olarak bildirilmiştir. Araştırmacılar bağıntının yalnız deneyin kapsadığı 70.000≤Re≤140.000, 9≤u≤15 m/s, 450–550 °C ve 0,2 MPa koşullarında kullanılmasını önermektedir. Dolayısıyla sonuç, farklı basınçlar veya süperkritik CO₂ koşullarına doğrudan genellenmemelidir.

Türkiye açısından çalışma; güneş ısısı, elektrikten ısıya dönüşüm, endüstriyel atık ısı veya yenilenebilir kaynaklardan gelen fazla elektriğin yüksek sıcaklıklı duyulur ısı depolama sistemlerinde değerlendirilmesi üzerinde çalışan mühendislik uygulamaları açısından yöntemsel önem taşımaktadır. Ancak araştırma Türkiye'de kurulmuş bir enerji depolama tesisi, yerel enerji maliyeti veya gerçek bir Türk santralinin işletme verilerini incelememektedir. Dökme demir–CO₂ yaklaşımının Türkiye'deki herhangi bir tesis için uygulanabilirliği; sistem sıcaklığı, depolama kapasitesi, basınçlı CO₂ altyapısı, malzeme maliyeti, pompalama ihtiyacı, güvenlik ve bütün çevrim verimi dikkate alınarak ayrıca doğrulanmalıdır.

Araştırmanın temel sorusu nedir?

Çalışmanın temel sorusu, yüksek sıcaklığa kadar ısıtılmış dökme demir katı depolama ortamından kritik altı CO₂ akışıyla ısı çekilirken taşınımla ısı transferinin nasıl değiştiği ve bu davranışın mühendislik tasarımında kullanılabilecek güvenilir bir Nusselt korelasyonuyla ifade edilip edilemeyeceğidir.

Katı duyulur ısı depolama sistemlerinde depolanan enerjinin kullanılabilir hale gelmesi yalnız depolama malzemesinin toplam ısı kapasitesine bağlı değildir. Isının akışkana hangi hızla aktarılabildiği, deşarj süresini ve buna bağlı olarak gerekli depolama malzemesi miktarını da etkiler. Araştırmacılar bu nedenle CO₂ ile dökme demir arasındaki zorlanmış taşınım katsayısının doğru tahmin edilmesini sistem boyutlandırması açısından kritik görmektedir.

Neden dökme demir ısıl depolama ortamı kullanılmıştır?

Çalışmada katı duyulur ısı depolama malzemesi olarak dökme demir seçilmiştir. Araştırmacılar bu seçimi malzemenin yüksek ısıl iletkenliği, mekanik dayanımı ve yüksek sıcaklıklardaki yapısal kararlılığıyla gerekçelendirmektedir.

Taş veya beton gibi geleneksel katı depolama ortamlarına kıyasla daha yüksek ısıl iletkenliğe sahip metalik bir ortamın, depolanan enerjinin daha hızlı alınmasına olanak verebileceği düşünülmektedir. Bununla birlikte çalışma dökme demirin başka depolama malzemelerinden ekonomik olarak daha üstün olduğunu doğrudan test etmemiştir; karşılaştırma deneyinin ana konusu malzeme maliyeti değil, dökme demir içindeki CO₂ kanallarının ısı transferidir.

Deney sistemi nasıl çalışmaktadır?

Deney düzeneği kapalı çevrim olarak tasarlanmıştır. Elektrikli ısıtıcı önce dökme demir depolama ortamını yüksek sıcaklığa getirmektedir. Daha sonra değişken frekanslı pompa CO₂'yi depolama ortamındaki iç kanallardan geçirir. CO₂ sıcak dökme demirden ısı aldıktan sonra doymuş evaporatöre yönlendirilir ve taşıdığı ısıyı besleme suyuna aktarır.

Sistemde ölçüm ve veri toplama için CO₂ hattında debi ölçer, basınç transmitterleri ve sıcaklık transmitterleri bulunmaktadır. Deney düzeneğinin şemasında CO₂, su ve oluşan buhar devreleri ayrı hatlar olarak gösterilmiş; ölçüm elemanları veri toplama sistemi ve bilgisayara bağlanmıştır.

Önemli bir fiziksel ayrım, CO₂'nin bütün dolaşım çevrimi boyunca tek fazlı gaz halinde kalmasıdır. Evaporatördeki faz değişimi su tarafındadır. Bu nedenle araştırmanın ölçtüğü olay CO₂'nin kaynama veya yoğuşması değil, sıcak katı yüzeyden gaz fazındaki CO₂'ye zorlanmış taşınımla ısı aktarımıdır.

Verianla Live: Deneyde ısının izlediği yol

Bu süreç gösterimi yalnız araştırmada tanımlanan kapalı çevrim deney sistemini özetlemektedir. CO₂'nin faz değiştirdiği varsayılmamıştır; kaynak çalışmada CO₂'nin çevrim boyunca tek fazlı gaz olarak kaldığı açıkça belirtilmiştir.

AşamaİşlemKaynak
Isıl şarjDökme demir ısıl depolama ortamı elektrikli ısıtıcıyla yüksek sıcaklığa çıkarılır.Deneysel yöntem, Bölüm 2.1
CO₂ dolaşımıDeğişken frekanslı pompa CO₂'yi kapalı çevrimde hareket ettirir.Şekil 1
Isı çekimiCO₂, dökme demir içindeki dairesel kanallardan geçerken zorlanmış taşınımla ısı alır.Bölüm 2.1–2.3
ÖlçümGiriş/çıkış sıcaklıkları, basınç ve debi izlenir.Şekil 1
Evaporatöre ısı aktarımıIsınan CO₂ taşıdığı enerjiyi evaporatörde besleme suyuna aktarır.Bölüm 2.1
CO₂'nin geri dönüşüSoğuyan CO₂ pompa üzerinden yeniden depolama ortamına döner; CO₂ çevrimi tek fazlı kalır.Bölüm 2.1
 

Dökme demir depolama ünitesinin geometrisi nasıldır?

Isıl depolama ortamı, içinde dairesel akış kanalları bulunan masif bir dökme demir bloktan oluşmaktadır. Kanallar şaşırtmalı düzende ve blok boyunca düzenli biçimde dağıtılmıştır. Ayrı bir iç kanal elektrikli ısıtıcının yerleştirilmesi için kullanılmaktadır.

Kaynak çalışmanın gerçek deney düzeneğine ait fotoğrafları, çok sayıda büyük çaplı metal borunun dökme demir blok içerisinden geçtiğini göstermektedir. Depolama ünitesinin dış tarafında ısıl yalıtım ve ısıtıcı bağlantıları da görülmektedir.

Geometrik parametreSembolDeğer
Isı transfer kanalı çapıd1200 mm
Isıtma kanalı çapıd2253 mm
Isıtma kanalı uzunluğuL1820 mm
Isıl depolama ünitesi yüksekliğih1300 mm
Kısa simetrik delik aralığıδ1600 mm
Uzun simetrik delik aralığıδ2990 mm
Delik–dış duvar mesafesiδ3125 mm
Delik merkezleri arası mesafeS1250 mm
Isıl depolama ünitesi sayısıN32
Toplam etkin ısı transfer uzunluğuL26.240 mm

Birden fazla depolama ünitesi seri bağlanarak gerekli toplam ısı transfer uzunluğu oluşturulmuştur. Çalışmada 32 özdeş ünitenin toplam etkin kanal uzunluğu 26,24 m'dir. Kanallar geometrik olarak özdeş kabul edildiği için ayrıntılı ısı transferi hesabında tek bir kanal temsilî test kesiti olarak kullanılmıştır.

CO₂ özellikleri nasıl belirlenmiştir?

Yoğunluk, dinamik ve kinematik viskozite, sabit basınç özgül ısı kapasitesi ve ısıl iletkenlik gibi CO₂ termofiziksel özellikleri NIST Chemistry WebBook SRD 69 verilerinden alınmıştır. İlgili özellikler hesaplamada giriş, çıkış veya aritmetik ortalama sıcaklığa göre değerlendirilmiştir.

Deney kesitinin iyi yalıtılmış olduğu belirtilmiş ve statik depolama ortamı ile bağlantı borularındaki ısı kayıplarının CO₂ tarafından alınan enerjiye göre ihmal edilebilir olduğu kabul edilmiştir. Bu nedenle veri indirgemesinde test kesiti yaklaşık adyabatik olarak ele alınmıştır.

CO₂'nin aldığı ısı nasıl hesaplanmıştır?

Kararlı durumda CO₂ tarafından alınan ısı hızı, kütlesel debi ile giriş ve çıkış entalpi yaklaşımının özgül ısı–sıcaklık biçiminde ifade edilmesiyle hesaplanmıştır:

\[ Q_c=q_m\left(c_{out}T_{out}-c_{in}T_{in}\right) \]

Burada \(Q_c\) CO₂'nin ısı alma hızını, \(q_m\) kütlesel debiyi, \(c_{in}\) ve \(c_{out}\) ilgili sıcaklıklardaki özgül ısı kapasitelerini, \(T_{in}\) ve \(T_{out}\) ise CO₂ giriş ve çıkış sıcaklıklarını temsil etmektedir.

Konvektif ısı transfer katsayısı genel enerji dengesi:

\[ Q_c=h_fA\Delta t \]

üzerinden belirlenmiştir. Ardından deneysel Nusselt sayısı:

\[ Nu_{exp}=\frac{h_fd_1}{\lambda} \]

şeklinde hesaplanmıştır. Burada \(h_f\) taşınımla ısı transfer katsayısını, \(d_1\) kanal çapını ve \(\lambda\) CO₂'nin ısıl iletkenliğini ifade etmektedir.

Reynolds ve Prandtl sayıları nasıl tanımlanmıştır?

Deneyde Reynolds sayısı:

\[ Re=\frac{ud_1}{\nu} \]

bağıntısıyla hesaplanmıştır. Burada \(u\) kanal içindeki akış hızı ve \(\nu\) CO₂'nin kinematik viskozitesidir.

Prandtl sayısı ise:

\[ Pr=\frac{\mu c_p}{\lambda} \]

olarak tanımlanmıştır. Prandtl hesabındaki termofiziksel özellikler giriş ve çıkış sıcaklıklarının aritmetik ortalamasında değerlendirilmiştir.

Deney hangi çalışma aralığında gerçekleştirilmiştir?

ParametreAraştırmada verilen aralık/değer
CO₂ basıncı0,2 MPa
Kütlesel debi33.994–53.764 kg/h
Reynolds sayısı7,5×104–1,35×105
Akış rejimiYüksek Reynolds sayılı zorlanmış taşınım
Isıl depolama ünitesi32 adet seri bağlı dökme demir ünite

Deneysel belirsizlikler ne düzeydedir?

Araştırmacılar doğrudan ölçümlerde Tip A ve Tip B belirsizliklerini, dolaylı hesaplanan büyüklüklerde ise hata yayılımını dikkate almıştır. Kaynakta verilen ana belirsizlikler şöyledir:

Ölçülen/hesaplanan büyüklükBelirsizlik
Basınç%0,10
CO₂ giriş/çıkış sıcaklığı%0,42
Kütlesel debi%0,16
Duvar sıcaklığı%0,45
Isı akısı%2,61
Isı transfer katsayısı%3,63

Bu belirsizlik analizi deneyin güçlü yönlerinden biridir; çünkü önerilen korelasyon yalnız bir nümerik modelden değil, ölçüm cihazlarının hata sınırları açıkça değerlendirilen fiziksel testlerden türetilmiştir.

Sıcaklık dağılımları ne göstermiştir?

CO₂ sıcaklığı depolama ortamının girişinden çıkışına doğru sürekli artmıştır. Bu eğilim, akışkanın kanal boyunca sıcak dökme demirden enerji aldığını doğrudan göstermektedir. Buna karşılık ortalama duvar sıcaklığının eksen boyunca görece düzgün kaldığı raporlanmıştır.

Kaynak çalışmanın sıcaklık grafiğinde düşük, orta ve yüksek Reynolds sayısını temsil eden üç çalışma vakası boyunca CO₂'nin eksenel sıcaklığının kademeli biçimde yükseldiği görülmektedir. Ayrı duvar sıcaklığı grafiğinde Case A'nın ortalama duvar sıcaklığı yaklaşık 645 °C, Case B'nin yaklaşık 627 °C ve Case C'nin yaklaşık 597 °C civarındadır.

Araştırmacılar duvar ve akışkan sıcaklıklarında ani sıçramaların görülmemesini kararlı çalışma koşullarının ve dökme demirin yüksek ısıl iletkenliğinin göstergesi olarak yorumlamaktadır.

Case B'deki Reynolds sayısı neden ayrıca not edilmelidir?

Kaynağın Tablo 3'ünde Case A için Re=7,6×104, Case B için Re=8,7×105 ve Case C için Re=1,3×105 yazılmıştır. Ancak aynı çalışma deney aralığını 7,5×104–1,35×105 olarak tanımlamaktadır. Case B'nin 8,7×105 değeri bu aralığın çok üzerindedir ve Case C'den de büyüktür.

Bu nedenle Case B değeri kaynak içi sayısal tutarsızlık olarak değerlendirilmelidir. Burada olası bir yazım hatası varsayılarak 8,7×104 biçiminde sessizce düzeltilmemiştir.

Reynolds sayısı arttığında Nusselt sayısı nasıl değişmiştir?

Deneyin en açık bulgularından biri Nusselt sayısının Reynolds sayısıyla monoton biçimde artmasıdır. Reynolds sayısı yaklaşık 75.000'den 135.000'e yükselirken Nusselt sayısı yaklaşık 74'ten 97'ye çıkmıştır.

Düşük Reynolds bölgesinde veri noktalarında belirli saçılmalar bulunurken yüksek Reynolds bölgesinde genel artış eğilimi daha belirgin hale gelmektedir. Araştırmacılar bunu artan akış hızıyla türbülanslı taşınımın güçlenmesine bağlamaktadır.

Verianla Live: Reynolds sayısı arttıkça Nusselt sayısındaki değişim

Aşağıdaki iki değer, çalışmanın Şekil 4 ve sonuç metninde verilen yaklaşık uç değerlerdir. Ara değerler üretilmemiştir.

Çalışma noktasıReynolds sayısı ReNusselt sayısı NuKaynak
İncelenen aralığın düşük ucu≈75.000≈74Şekil 4 ve Bölüm 3.2
İncelenen aralığın yüksek ucu≈135.000≈97Şekil 4 ve Bölüm 3.2
 

Klasik ısı transfer korelasyonları neden yeterli olmamıştır?

Araştırma Dittus–Boelter, Gnielinski, Colburn ve Zhou korelasyonlarını referans olarak kullanmıştır. Bu bağıntıların önemli bölümü düzgün dairesel borularda, belirli termal sınır koşulları ve geleneksel akışkan özellikleri altında geliştirilmiştir.

Çalışmanın karşılaştırma grafiklerinde klasik bağıntılardan hesaplanan Nusselt değerleri deneysel Nusselt değerlerinin önemli ölçüde üzerinde kalmaktadır. Özellikle Dittus–Boelter ve Colburn bağıntıları, deneysel yaklaşık 74–97 aralığına karşı çok daha yüksek tahminler üretmektedir. Zhou bağıntısı klasik bağıntıların bazılarına göre deney verilerine yaklaşsa da yine belirgin sistematik sapma göstermektedir.

Araştırmacılar bu farkı yüksek sıcaklıklı katı depolama ortamının ideal olmayan termal sınır koşulları, yüksek Reynolds sayısı ve CO₂'nin sıcaklığa bağlı termofiziksel özellikleri gibi etkenlerin klasik bağıntıların geliştirildiği koşullardan farklı olmasıyla ilişkilendirmektedir.

Hangi yeni Nusselt korelasyonları geliştirilmiştir?

Deney verilerinden önce klasik kuvvet yasası yapısında üç bağıntı, ardından logaritmik uzayda ikinci dereceden genişletilmiş bir bağıntı türetilmiştir:

Kaynakta verilen denklemAmpirik bağıntı
Denklem (14)\(Nu=0.445Re^{0.458}\)
Denklem (15)\(Nu=0.514Re^{0.456}Pr^{0.4}\)
Denklem (16)\(Nu=0.202Re^{0.465}Pr^{-2.19}\)
Denklem (17)\(\ln Nu=-34.645+6.580\ln Re-0.267(\ln Re)^2+3.766\ln Pr\)

İlk üç ilişki, mühendislikte yaygın kullanılan Reynolds–Prandtl kuvvet yasası biçimini korumaktadır. Dördüncü ilişki ise yüksek Reynolds bölgesindeki eğriliği daha iyi yakalamak amacıyla \((\ln Re)^2\) terimi eklenmiş daha esnek bir ampirik modeldir.

Neden matematiksel olarak en yüksek R²'ye sahip model seçilmemiştir?

Araştırmacılar yalnız belirleme katsayısını maksimize etmeyi yeterli görmemiştir. Genişletilmiş polinom model en yüksek R² değerine sahip olmasına rağmen ortalama ve maksimum tahmin hatalarının daha büyük olduğu belirtilmiştir.

Çalışmanın Tablo 6'sında dördüncü satır için R²=0,994, MAPE=%4,137 ve maksimum sapma=%8,083 verilmiştir. Bu nedenle daha karmaşık modelin eğitim verisine daha iyi uyum göstermesine karşın mühendislik tahmininde hata davranışının daha zayıf olduğu sonucuna varılmıştır.

Önerilen nihai korelasyon nedir?

Araştırmanın önerdiği nihai bağıntı:

\[ \boxed{Nu=0.202Re^{0.465}Pr^{-2.19}} \]

şeklindedir.

Bu korelasyon için çalışma:

  • R²: 0,988
  • MAPE: %0,715
  • Maksimum bağıl sapma: %2,404

değerlerini bildirmektedir. Araştırmacılar bu bağıntının uyum doğruluğu ile en kötü durum hatası arasında daha dengeli sonuç verdiğini belirterek mühendislik hesabı için bunu önermiştir.

Denklem numaralarında hangi tutarsızlık bulunmaktadır?

Kaynağın Bölüm 3.4.2'sinde dört yeni bağıntı açık biçimde Denklem (14), (15), (16) ve (17) olarak numaralandırılmıştır. Ancak hemen sonraki Tablo 6'nın satırları Eq.12, Eq.13, Eq.14 ve Eq.15 biçimindedir. Sayfa 45–48'de ayrı verilen tahmin grafikleri de aynı kaymış Eq.12–Eq.15 etiketlerini kullanmaktadır.

Buna karşılık Bölüm 3.4.4 ve sonuç bölümü nihai bağıntının Denklem (16) olduğunu açıkça belirtmektedir. Ayrıca özet bölümündeki R²=0,988, MAPE=%0,715 ve maksimum sapma=%2,404 değerleri Tablo 6'nın üçüncü satırındaki değerlerle aynıdır.

Dolayısıyla bu Verianla metninde önerilen matematiksel bağıntı kaynakta Denklem (16) olarak verilen haliyle korunmuş, ancak Tablo 6 ve Şekil 9'daki iki basamaklık numaralandırma uyuşmazlığı sessizce düzeltilmemiştir.

Önerilen korelasyon hangi koşullarda kullanılabilir?

Çalışmanın sonuç bölümünde açıkça verilen geçerlilik alanı şöyledir:

DeğişkenGeçerlilik aralığı
Reynolds sayısı70.000 ≤ Re ≤ 140.000
CO₂ akış hızı9 ≤ u ≤ 15 m/s
Sıcaklık450 ≤ T ≤ 550 °C
Basınç0,2 MPa

Bu sınırlar korelasyonun yorumlanmasında kritik önemdedir. Bağıntının 0,2 MPa dışındaki basınçlara, daha farklı sıcaklıklara, farklı kanal geometrilerine veya süperkritik CO₂'ye doğrudan uygulanabileceği çalışma tarafından gösterilmemiştir.

Çalışmanın desteklediği sonuçlar

  • Kritik altı CO₂, incelenen kapalı çevrim deney sisteminde yüksek sıcaklıklı dökme demir depolama ortamından zorlanmış taşınımla ısı çekmiştir.
  • İncelenen Reynolds aralığında Reynolds sayısı arttıkça Nusselt sayısı sistematik biçimde yükselmiştir.
  • Yaklaşık Re=75.000–135.000 aralığında deneysel Nu yaklaşık 74'ten 97'ye yükselmiştir.
  • Klasik Dittus–Boelter, Gnielinski, Colburn ve Zhou korelasyonları mevcut deney koşullarında belirgin tahmin sapmaları göstermiştir.
  • Deney verilerinden dört yeni ampirik Nusselt korelasyonu geliştirilmiştir.
  • \(Nu=0.202Re^{0.465}Pr^{-2.19}\) bağıntısı doğruluk ve maksimum hata arasında en uygun dengeyi sağladığı gerekçesiyle önerilmiştir.
  • Önerilen korelasyon için R²=0,988, MAPE=%0,715 ve maksimum bağıl sapma=%2,404 raporlanmıştır.

Çalışmanın kanıtlamadığı sonuçlar

  • Çalışma dökme demir depolamanın bütün alternatif ısıl depolama malzemelerinden ekonomik olarak üstün olduğunu kanıtlamamaktadır.
  • Deney, süperkritik CO₂ ısı transferini incelememektedir.
  • Önerilen korelasyon farklı basınç ve kanal geometrileri için doğrulanmamıştır.
  • Çalışma tam bir enerji depolama santralinin elektrikten-elektriğe çevrim verimini ölçmemiştir.
  • Pompalama elektriği, toplam sistem ekonomisi ve seviyelendirilmiş enerji depolama maliyeti hesaplanmamıştır.
  • Dökme demirin çok uzun süreli termal çevrimler sonucundaki mekanik veya mikroyapısal yaşlanması araştırılmamıştır.
  • Sonuçlar Türkiye'deki belirli bir tesis veya santral için doğrudan performans garantisi değildir.

Çalışmanın Yöntemi ve Bulguları

Deneysel tasarımın özeti

ÖğeÇalışmada kullanılan yapı
Isıl depolama ortamıDairesel iç kanallara sahip masif dökme demir bloklar
Isı transfer akışkanıKritik altı CO₂, tek fazlı gaz
ÇevrimKapalı döngü
Basınç0,2 MPa
Isıl şarjElektrikli ısıtıcı
CO₂ dolaşımıDeğişken frekanslı pompa
Isı alıcısıDoymuş evaporatörün su tarafı
Isı transfer kanal çapı200 mm
Depolama ünitesi sayısı32
Toplam etkin kanal uzunluğu26,24 m
Kütlesel debi aralığı33.994–53.764 kg/h
Reynolds aralığı7,5×104–1,35×105

Üç temsilî deney vakası

Çalışma sıcaklık profillerini tartışmak için düşük, orta ve yüksek Reynolds bölgelerinden üç temsilî vaka seçmiştir. Kaynak Tablo 3'te değerler aşağıdaki biçimde basılmıştır:

VakaReKütlesel debiBasınçGiriş sıcaklığı
Case A7,6×1043,3×104 kg/h0,2 MPa443,1 °C
Case B8,7×105 — kaynakta bu şekilde yazılmıştır3,6×104 kg/h0,2 MPa430,8 °C
Case C1,3×1055,3×104 kg/h0,2 MPa425,4 °C

Case B'deki Reynolds değeri araştırmanın genel çalışma aralığıyla uyuşmadığından düzeltilmiş bir değer olarak kullanılmamalıdır. Bu tablo yalnız kaynağın basılı halini belgelemek amacıyla verilmiştir.

Nusselt–Reynolds eğiliminin fiziksel yorumu

Reynolds sayısı arttığında akışın atalet etkileri viskoz etkilere göre güçlenmekte ve araştırmanın incelenen türbülanslı çalışma aralığında CO₂ ile kanal duvarı arasındaki taşınım yoğunlaşmaktadır. Deneyde bunun karşılığı Nusselt sayısındaki monoton artıştır.

Araştırmacılar düşük Reynolds tarafındaki saçılmaları duvar sıcaklığının ideal biçimde tamamen sabit olmamasına ve yerel sıcaklık farkındaki küçük değişimlerin hesaplanan Nusselt değerlerini etkileyebilmesine bağlamaktadır. Yüksek Reynolds tarafında ise artış eğilimi daha güçlü olmasına rağmen küçük yerel termal değişiklikler yine veri saçılımına katkıda bulunmaktadır.

Klasik korelasyon karşılaştırmalarının temel mesajı

Dittus–Boelter karşılaştırma grafiğinde klasik modelin deneysel yaklaşık 74–98 aralığındaki Nu değerlerine karşı yaklaşık 160–250 düzeyine kadar tahminler ürettiği görülmektedir. Benzer sistematik yüksek tahmin Gnielinski ve Colburn grafiklerinde de bulunmaktadır.

Zhou bağıntısı sıcaklık oranını içeren daha gelişmiş bir ifade olmasına rağmen deneysel değerlere tam olarak oturmamaktadır. Bu sonuçlar, araştırmacıların mevcut yüksek sıcaklıklı dökme demir–CO₂ sistemine özgü bir regresyon bağıntısı geliştirme gerekçesini oluşturmaktadır.

Regresyon yöntemi

İlk korelasyonlar genel:

\[ Nu=aRe^bPr^c \]

yapısından hareketle geliştirilmiştir. Deney verileri logaritmik uzaya dönüştürülmüş ve katsayılar sıradan en küçük kareler temelli çoklu doğrusal regresyonla belirlenmiştir.

Yüksek Reynolds bölgesindeki eğriliğin klasik kuvvet yasasıyla tam olarak temsil edilemediği görülünce \(\ln(Re)\) değişkeninin ikinci dereceden terimini içeren genişletilmiş bir polinom modeli de oluşturulmuştur. Araştırmacılar bunun deney aralığı içinde daha esnek uyum sağlayabildiğini, ancak aralık dışına ekstrapolasyonda daha düşük sağlamlığa sahip olabileceğini belirtmektedir.

Önerilen korelasyonun istatistiksel değerlendirmesi

Önerilen ilişki için raporlanan göstergeDeğer
Belirleme katsayısı R²0,988
Ortalama mutlak yüzde hata, MAPE%0,715
Maksimum bağıl sapma%2,404

Çalışmanın parity grafiklerinde yeni korelasyonlardan elde edilen tahminlerin deneysel değerler çevresinde toplandığı ve bütün veri noktalarının ±%10 bandı içinde kaldığı belirtilmektedir. Bununla birlikte korelasyonların eğitiminde ve değerlendirilmesinde aynı deney veri havuzunun nasıl ayrıldığına ilişkin bağımsız bir dış test veri seti veya çapraz doğrulama prosedürü açıklanmamıştır.

Deneyin güçlü yönleri

  • Sonuçların yalnız CFD veya teorik hesap yerine fiziksel kapalı çevrim deney sisteminden elde edilmesi.
  • Gerçek dökme demir ısıl depolama ünitesinin ve çoklu kanal mimarisinin kullanılması.
  • CO₂ sıcaklığı, basıncı ve debisinin gerçek zamanlı ölçülmesi.
  • Ölçüm ve türetilmiş büyüklükler için belirsizlik analizinin yapılması.
  • Klasik korelasyonlarla doğrudan karşılaştırma sunulması.
  • Ampirik ilişkinin yalnız uyum katsayısıyla değil MAPE ve maksimum sapmayla da değerlendirilmesi.
  • Önerilen bağıntının geçerlilik aralığının açık biçimde sınırlandırılması.

Temel yöntemsel sınırlılıklar

Deney yalnız 0,2 MPa basınçta gerçekleştirilmiştir. Bu nedenle basıncın ısı transferi üzerindeki bağımsız etkisi parametrik olarak test edilmemiştir.

Çalışma tek bir temel kanal geometrisine ve 200 mm kanal çapına dayanmaktadır. Kanal çapı, yüzey pürüzlülüğü veya farklı iç geometrilerin etkisi deneysel olarak karşılaştırılmamıştır.

Önerilen bağıntı mevcut deney veri setinden türetilmiş ve aynı deney aralığında değerlendirilmiştir. Başka bir laboratuvar, farklı dökme demir yapısı veya bağımsız bir veri setiyle dış doğrulama sunulmamıştır.

Isı kayıpları iyi yalıtım nedeniyle ihmal edilebilir kabul edilmiştir; ancak çalışma bütün çevrim için ayrıntılı bir enerji kaybı bütçesi veya uzun dönem dinamik şarj-deşarj verim haritası sunmamaktadır.

Çalışmanın hedefi yerel/ortalama taşınımla ısı transferinin mühendislik korelasyonunu geliştirmektir. Tam sistem seviyesi depolama verimi, pompa tüketimi, ekipman maliyeti, yaşam döngüsü maliyeti veya çevresel performans bu deneyin doğrudan çıktıları değildir.

Kaynak ve Yöntem Notu

Tam özgün çalışma adı: Experimental investigation of forced convective heat transfer of subcritical CO₂ in high-temperature cast-iron thermal storage channels

Yazarlar ve sıralama: Hongbo Wang; Cong Liu; Qimin Wang; Xiaohui Zhang.

Sorumlu yazar: Cong Liu.

Sorumlu yazar e-posta adresi: liucong@sie.edu.cn

Eş katkı/eş birinci yazarlık: Yüklenen sürümde eş katkı veya eş birinci yazarlık beyanı bulunmamaktadır.

Kurum: Shenyang Institute of Engineering, Shenyang 110136, China.

Kaynak türü: Yüksek sıcaklıklı katı duyulur ısı depolama sisteminde kritik altı CO₂ taşınımla ısı transferini inceleyen deneysel araştırma preprinti.

Platform: SSRN.

SSRN kayıt numarası: 6934134.

DOI: 10.2139/ssrn.6934134.

Resmî bağlantı:SSRN 6934134

Hakemlik durumu: İncelenen sürüm hakem değerlendirmesinden geçmemiş bir preprinttir. Bibliyografik kontrolde aynı başlık için kesin biçimde doğrulanmış ayrı bir hakemli dergi sürümü tespit edilmediğinden herhangi bir dergi veya hakemli DOI bilgisi eklenmemiştir.

Yayın yılı: 2026.

Deneysel yöntem: 0,2 MPa basınçta kapalı çevrim kritik altı CO₂ sistemi; elektrikle ısıtılmış dökme demir duyulur ısı depolama ortamı; 32 seri depolama ünitesi; değişken frekanslı pompa; sıcaklık, basınç ve debi ölçümleri; Reynolds, Prandtl ve Nusselt hesapları; deneysel belirsizlik analizi ve ampirik regresyon.

Termofiziksel veri kaynağı: Araştırmacılar CO₂'nin yoğunluk, viskozite, özgül ısı ve ısıl iletkenlik değerleri için NIST Chemistry WebBook SRD 69 kullanmıştır.

Önerilen korelasyon: \(Nu=0.202Re^{0.465}Pr^{-2.19}\).

Önerilen korelasyonun çalışma tarafından bildirilen performansı: R²=0,988; MAPE=%0,715; maksimum bağıl sapma=%2,404.

Geçerlilik sınırı: 70.000≤Re≤140.000; 9≤u≤15 m/s; 450≤T≤550 °C; p=0,2 MPa. Bu sınırların dışındaki kullanım çalışma tarafından deneysel olarak doğrulanmamıştır.

Kaynak içi sayısal tutarsızlık: Tablo 3'te Case B için Reynolds sayısı 8,7×105 olarak verilmiştir; bu değer çalışmanın tanımladığı 7,5×104–1,35×105 deney aralığıyla uyuşmamaktadır. Değer sessizce düzeltilmemiştir.

Kaynak içi numaralandırma tutarsızlığı: Dört yeni korelasyon ana metinde Denklem 14–17 olarak verilirken Tablo 6 ve ilgili tahmin grafiklerinde Eq.12–Eq.15 etiketleri kullanılmaktadır. Metin ve sonuç bölümü nihai öneriyi Denklem 16 olarak tanımladığından matematiksel bağıntı aynen korunmuş, numara uyuşmazlığı ayrıca belirtilmiştir.

Finansman: Yüklenen sürümde ayrı bir finansman beyanı tespit edilmemiştir.

Çıkar çatışması: Yüklenen sürümde ayrı bir çıkar çatışması beyanı tespit edilmemiştir.

Veri erişilebilirliği: Yüklenen sürümde ayrı bir veri erişilebilirliği beyanı bulunmamaktadır.

CRediT/yazar katkıları: Yüklenen sürümde ayrıntılı CRediT yazar katkı beyanı bulunmamaktadır.

Üretken yapay zekâ beyanı: Araştırmacılar çalışma metninin hazırlanması sırasında okunabilirlik ve mantıksal açıklığı geliştirmek amacıyla ChatGPT kullandıklarını, sonrasında içeriği gözden geçirip düzenlediklerini ve çalışmanın içeriğine ilişkin sorumluluğu üstlendiklerini belirtmektedir.

Lisans: İncelenen sürüm üzerinde Creative Commons veya eşdeğer açık yeniden kullanım lisansı tespit edilmemiştir. Bu nedenle kaynak çalışmanın özgün deney düzeneği fotoğrafları ve grafiklerinin bire bir kopyalanmasına dayalı bir görsel kullanım varsayılmamıştır.

Temel bilimsel sınır: Bulgular belirli bir 200 mm kanal geometrisi, 32 ünitelik dökme demir depolama sistemi ve 0,2 MPa kritik altı CO₂ koşullarına aittir. Korelasyonun farklı basınçlara, malzemelere, kanal geometrilerine veya süperkritik çalışma rejimine aktarılması ayrıca doğrulanmalıdır.

Bilimsel içerik ilkesi: Bu Verianla açıklamasındaki deney tasarımı, formüller, ölçüm değerleri, korelasyonlar, istatistikler ve sınırlılıklar incelenen çalışmaya dayanmaktadır. Dış bibliyografik kontrol yalnız SSRN kimliği ve DOI doğrulaması amacıyla kullanılmış; başka araştırmalardan yeni bilimsel sonuç ana anlatıya eklenmemiştir.


Paylaş:

Yorumlar incelendikten sonra yayımlanır.Gönderdiğiniz yorum onay sürecine alınır ve uygun bulunduğunda görünür hâle gelir.

Bir yorum bırakın

E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir

Your experience on this site will be improved by allowing cookies Cookie Policy