
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 | İşlem | Kaynak |
|---|---|---|
| Isıl şarj | Dö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ı çekimi | CO₂, dökme demir içindeki dairesel kanallardan geçerken zorlanmış taşınımla ısı alır. | Bölüm 2.1–2.3 |
| Ölçüm | Giriş/çı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 parametre | Sembol | Değer |
|---|---|---|
| Isı transfer kanalı çapı | d1 | 200 mm |
| Isıtma kanalı çapı | d2 | 253 mm |
| Isıtma kanalı uzunluğu | L1 | 820 mm |
| Isıl depolama ünitesi yüksekliği | h | 1300 mm |
| Kısa simetrik delik aralığı | δ1 | 600 mm |
| Uzun simetrik delik aralığı | δ2 | 990 mm |
| Delik–dış duvar mesafesi | δ3 | 125 mm |
| Delik merkezleri arası mesafe | S1 | 250 mm |
| Isıl depolama ünitesi sayısı | N | 32 |
| Toplam etkin ısı transfer uzunluğu | L | 26.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?
| Parametre | Araştırmada verilen aralık/değer |
|---|---|
| CO₂ basıncı | 0,2 MPa |
| Kütlesel debi | 33.994–53.764 kg/h |
| Reynolds sayısı | 7,5×104–1,35×105 |
| Akış rejimi | Yüksek Reynolds sayılı zorlanmış taşınım |
| Isıl depolama ünitesi | 32 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ük | Belirsizlik |
|---|---|
| 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ı Re | Nusselt sayısı Nu | Kaynak |
|---|---|---|---|
| İ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 denklem | Ampirik 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şken | Geçerlilik aralığı |
|---|---|
| Reynolds sayısı | 70.000 ≤ Re ≤ 140.000 |
| CO₂ akış hızı | 9 ≤ u ≤ 15 m/s |
| Sıcaklık | 450 ≤ 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 |
| Çevrim | Kapalı döngü |
| Basınç | 0,2 MPa |
| Isıl şarj | Elektrikli ı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ğu | 26,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:
| Vaka | Re | Kütlesel debi | Basınç | Giriş sıcaklığı |
|---|---|---|---|---|
| Case A | 7,6×104 | 3,3×104 kg/h | 0,2 MPa | 443,1 °C |
| Case B | 8,7×105 — kaynakta bu şekilde yazılmıştır | 3,6×104 kg/h | 0,2 MPa | 430,8 °C |
| Case C | 1,3×105 | 5,3×104 kg/h | 0,2 MPa | 425,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österge | Değ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.

Bir yorum bırakın
E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir