Bu çalışma, düşük sıcaklıklı güneş ısısının veya atık ısının uzun süre saklanması için kullanılabilecek, stronsiyum bromür hidrat esaslı açık çevrimli bir termokimyasal enerji depolama reaktörünü üç boyutlu sayısal modelle incelemektedir. Reaktörde katı SrBr₂·H₂O parçacıkları yerçekimi etkisiyle aşağı doğru hareket ederken nemli hava yatağın içinden yatay yönde geçmektedir. Su buharı ile tuz arasındaki ekzotermik hidratasyon reaksiyonu, depolanmış kimyasal enerjiyi yeniden ısıya dönüştürmekte ve reaktörden çıkan havanın sıcaklığını yükseltmektedir.

Araştırmacılar COMSOL Multiphysics ortamında gözenekli ortam akışını, su buharı taşınımını, reaksiyon kinetiğini, katı ve gaz fazlarındaki ısı transferini ve hareketli tuz yatağını birleştiren üç boyutlu bir model geliştirmiştir. Katı parçacıkların aşağı doğru ilerlemesi, Arbitrary Lagrangian–Eulerian (ALE) adı verilen hareketli ağ yöntemiyle temsil edilmiştir. Böylece reaktör içindeki reaksiyon cephesi, sıcaklık alanı, su buharı dağılımı, basınç kaybı ve yerel tuz dönüşümü zamana ve konuma bağlı olarak izlenebilmiştir.

Modelin temel sonucu, reaktör performansını belirleyen en güçlü işletme değişkeninin giriş havasındaki nem miktarı olduğudur. Giriş nemi 4 g/kg’dan 12 g/kg’a yükseltildiğinde yatak sıcaklığı yaklaşık 26°C’den 39°C’ye, havaya aktarılan yararlı hacimsel ısı ise yaklaşık 2 kW/m³’ten 5,5 kW/m³’e çıkmıştır. Daha yüksek nem, su buharı kısmi basıncını ve hidratasyonun termodinamik itici kuvvetini artırmış; daha sıcak, daha geniş ve daha belirgin bir reaksiyon bölgesi oluşturmuştur.

Hava debisinde en dengeli sonuç 110 m³/saat koşulunda elde edilmiştir. 80 m³/saatte hava tuzla daha uzun süre temas etse de toplam ısı taşıma kapasitesi düşük kalmıştır. 140 m³/saatte ise hava daha hızlı hareket etmiş, ancak kısa temas süresi ve güçlü konvektif soğutma nedeniyle reaksiyon ve yatak sıcaklığı sınırlanmıştır. Katı yatağın 0,5–2 cm/saat arasındaki hareket hızı sıcaklık ve anlık ısı çıkışını yalnızca sınırlı ölçüde değiştirmiş; buna karşılık daha hızlı katı hareketi tuzun reaktörde kalma süresini kısaltarak tamamlanmamış hidratasyon riskini artırmıştır.

Model daha önce yayımlanmış deneysel bir hareketli yatak reaktörü çalışmasının sıcaklık, nem ve ısı çıkışı verileriyle karşılaştırılmıştır. Eğilimler ve genel büyüklükler deneylerle uyumludur. Bununla birlikte PDF’deki doğrulama tablosunun doğrudan hesaplanması, bazı koşullardaki göreli farkların özette bildirilen yüzde 1–4 aralığını aşarak yararlı ısıda yüzde 8’e, yatak sıcaklığında yaklaşık yüzde 7’ye ulaşabildiğini göstermektedir.

Çalışma yeni bir fiziksel reaktör veya yeni deneysel veri sunmamaktadır. Bulgular, literatürdeki deneysel düzeneğin üç boyutlu modellenmesine ve parametrik sayısal simülasyonlara dayanmaktadır. Malzemenin uzun dönemli çevrim kararlılığı, aglomerasyon, aşınma, delikuesans, parçacık kırılması ve gerçek tesis kontrolü incelenmemiştir. Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.

Termokimyasal enerji depolama neden önemlidir?

Güneş enerjisi, endüstriyel atık ısı ve diğer yenilenebilir ısı kaynakları üretildikleri anda her zaman kullanılamaz. Özellikle bina ısıtmasında yazın veya gündüz elde edilen ısının haftalar ya da aylar sonra kullanılabilmesi, yüksek enerji yoğunluklu ve düşük bekleme kayıplı depolama yöntemleri gerektirir.

Isıl enerji depolama teknolojileri genel olarak üç gruba ayrılmaktadır:

• Duyulur ısı depolama: Su, taş veya beton gibi bir malzemenin sıcaklığının yükseltilmesine dayanır.
• Gizli ısı depolama: Faz değiştiren malzemelerin erime veya katılaşma enerjisini kullanır.
• Termokimyasal enerji depolama: Isıyı tersinir bir kimyasal reaksiyonun ürünlerinde depolar.

Termokimyasal depolamada şarj edilmiş reaksiyon ürünleri ortam sıcaklığında birbirinden ayrı tutulabildiği için uzun bekleme süresinde sürekli sıcak bir tankın korunması gerekmez. Enerji, sıcaklık farkı yerine kimyasal bileşimde saklanır. Bu nedenle uzun dönemli veya mevsimsel depolamada ısı kayıpları oldukça düşük olabilir.

Tuz hidratları nasıl ısı depolamaktadır?

Çalışmada kullanılan reaksiyon çifti stronsiyum bromür monohidrat ile stronsiyum bromür hekzahidrat arasındadır. PDF’de reaksiyon aşağıdaki biçimde verilmiştir:

\[ SrBr_2\cdot H_2O(s) + 5H_2O(g) \rightleftharpoons SrBr_2\cdot 6H_2O(s) + 5\Delta H \]

Şarj aşamasında hidratlı tuza sıcak ve kuru hava verilerek bağlı suyun bir bölümü uzaklaştırılır. Bu endotermik dehidratasyon için dışarıdan ısı sağlanır ve enerji susuzlaştırılmış tuz yapısında kimyasal potansiyel olarak depolanır.

Deşarj aşamasında ise soğuk ve nemli hava susuzlaştırılmış tuzla temas eder. Tuz, havadaki su buharını yeniden bağlayarak hidratlanır. Hidratasyon ekzotermik olduğu için açığa çıkan ısı havayı ısıtır. Isınmış çıkış havası ayrı bir ısı değiştirici üzerinden bina veya bölgesel ısıtma sistemine enerji aktarabilir.

Reaksiyonun gerçekleşebilmesi yalnızca havada su bulunmasına bağlı değildir. Havadaki gerçek su buharı basıncının, tuzun bulunduğu sıcaklıktaki denge buharı basıncını aşması gerekir. Bu nedenle nem, sıcaklık ve denge basıncı birbirine bağlıdır.

Açık ve kapalı termokimyasal sistemler arasındaki fark

Kapalı bir termokimyasal sistemde tuz ile saf su buharı çoğunlukla vakum altında reaksiyona girer. Bu yapı daha kontrollü olabilir; ancak buharlaştırıcı, yoğuşturucu, sıvı su deposu, vakum ekipmanı ve sızdırmazlık gerektirir.

İncelenen açık sistemde gaz reaktifi olarak atmosferik basınçtaki nemli hava kullanılmaktadır. Bu sayede ayrı bir su buharı çevrimine ihtiyaç azalır. Buna karşılık sistem dış ortam nemine çok daha duyarlı hâle gelir. Fazla nem tuzun yapışmasına, aglomerasyona veya delikuesansa; düşük nem ise yetersiz reaksiyon hızına neden olabilir.

Hareketli yataklı reaktör nasıl çalışmaktadır?

Reaktörün üst bölümünden kuru veya düşük hidratlı SrBr₂ parçacıkları beslenmektedir. Katı yatak yerçekimi etkisiyle dikey olarak aşağı doğru hareket eder. Alt taraftaki döner vana, katının reaktörden çıkış hızını ve dolayısıyla katı hızını kontrol eder.

Nemli hava reaktörün bir yanından girerek katı yatağa çapraz yönde ilerler ve karşı taraftan çıkar. Bu çapraz akış düzeninde:

• Katı faz yukarıdan aşağıya,
• Gaz faz soldan sağa,
• Su buharı havadan katıya,
• Reaksiyon ısısı katıdan ve reaksiyon bölgesinden havaya

taşınmaktadır.

PDF’deki reaktör şemasında, katı giriş hunisi üstte, hareketli reaktif yatak merkezde, katı çıkış mekanizması altta ve nemli hava bağlantıları reaktör gövdesinin iki yanında yer almaktadır. Merkezdeki dar dikey yatak, genişleyen iki yan hava dağıtım hacmiyle çevrilmektedir.

Reaktör geometrisinin temel ölçüleri

Reaktif yatağın üzerindeki iki paslanmaz çelik duvar, havanın üstteki boş bölgeden doğrudan çıkışa kaçmasını engellemek ve akışı tuz yatağına zorlamak amacıyla kullanılmıştır. Bu duvarlar aynı zamanda ek hidrolik direnç oluşturmuştur.

Araştırmanın hedeflediği bilimsel boşluk

Önceki hareketli yatak çalışmaları çoğunlukla iki boyutlu, yığılmış parametreli veya yüksek sıcaklıklı ve tersinmez reaksiyonlara yönelik modeller kullanmıştır. Bu yaklaşımlar aşağıdaki etkileri aynı anda ayrıntılı biçimde yakalamakta sınırlı kalmıştır:

• Katı yatağın sürekli hareketi,
• Üç boyutlu hava dağılımı ve yan kaçaklar,
• Su buharı kısmi basıncına bağlı reaksiyon kinetiği,
• Reaksiyon cephesinin uzay ve zaman içindeki ilerlemesi,
• Yerel gözeneklilik ve geçirgenlik etkileri,
• Katı kalış süresi ile dönüşüm arasındaki ilişki,
• Isı çıkışı ile toplam enerji depolama verimi arasındaki ödünleşme.

Çalışmanın amacı bu süreçleri tek bir üç boyutlu, sürekli işletim modelinde birleştirmek ve model sonuçlarını daha önce yayımlanmış deneysel reaktör verileriyle karşılaştırmaktır.

Modelde kullanılan çalışma koşulları

Modelin temel varsayımları

• Isı ve kütle taşınımını esas olarak hava akışı sürüklemektedir.
• Katı yatak çok yavaş hareket ettiği için hava dinamiği üzerindeki geri etkisi ihmal edilmiştir.
• Gaz akışı kararlı ve giriş kesitinde düzgün kabul edilmiştir.
• Katı parçacıklar sabit hızla aşağı doğru hareket etmektedir.
• Tuzun reaksiyon sırasında kimyasal reaktivitesini ve mekanik yapısını koruduğu varsayılmıştır.
• Yatak gözenekliliği, bypass etkisinin ayrıca tanımlandığı durumlar dışında sabit kabul edilmiştir.
• Nemli hava ideal gaz karışımı olarak modellenmiştir.
• Uzun süreli çevrimsel bozunma, topaklanma, parçacık kırılması ve aşınma modele dâhil edilmemiştir.

Su buharı ve tuz dönüşümünün korunumu

Gözenekli yataktaki su buharı dengesi, birikim, reaksiyonla tüketim ve gaz akışıyla taşınım terimlerinden oluşmaktadır. PDF’deki gösterimin özlü biçimi şöyledir:

\[ \varepsilon n_h\frac{\partial y_v}{\partial t} = \phi_v(1-y_v) - n_h\mathbf{u}_h\cdot\nabla y_v \]

• ε: Yatak gözenekliliği,
• n_h: Nemli havanın molar yoğunluğunu temsil eden terim,
• y_v: Su buharı mol veya kütle fraksiyonu,
• φ_v: Reaksiyon nedeniyle oluşan veya tüketilen buhar kaynak terimi,
• u_h: Nemli havanın hız vektörüdür.

Tuzun yerel hidratasyon dönüşümü için verilen denge:

\[ (1-\varepsilon)\frac{\partial X}{\partial t} = -\frac{\phi_v}{\nu} - \mathbf{u}_s\cdot\nabla X \]

• X: Yerel hidratasyon dönüşümü,
• ν: Stokiyometrik katsayı,
• u_s: Katı yatağın aşağı yönlü hızıdır.

İkinci denklemde ilk terim kimyasal reaksiyonla dönüşümü, ikinci terim ise dönüşmüş veya dönüşmemiş katının yatak boyunca taşınmasını temsil etmektedir.

Gözenekli ortam akışı: Brinkman denklemi

Reaktördeki hava akışı, Darcy yasasına viskoz kayma etkilerini ekleyen Brinkman denklemiyle modellenmiştir:

\[ \rho(\mathbf{u}\cdot\nabla)\mathbf{u} = \nabla\cdot\left[-p\mathbf{I} +\mu\left(\nabla\mathbf{u}+(\nabla\mathbf{u})^T\right)\right] -\mu\alpha\mathbf{u} +\rho\mathbf{F} \]

• ρ: Gaz karışımının yoğunluğu, kg/m³,
• u: Hız vektörü, m/s,
• p: Basınç, Pa,
• μ: Dinamik viskozite, Pa·s,
• α: Gözenekli ortamın Darcy direnciyle ilişkili terim, m⁻²,
• F: Birim kütleye etkiyen cisim kuvvetidir.

Denklem, yan hava girişinde oluşan hızlı jetleri, genişleme bölgelerindeki dolaşımı, tuz yatağındaki basınç kaybını ve çıkış kanalındaki viskoz direnci hesaplamak için kullanılmıştır.

Isı transferi ve reaksiyon ısısı

Gaz fazındaki enerji denklemi, konveksiyon, iletim ve hacimsel ısı kaynaklarını birlikte içermektedir:

\[ \rho C_p\mathbf{u}\cdot\nabla T+\nabla\cdot\mathbf{q} = Q+Q_p+Q_{vd} \]

\[ \mathbf{q}=-k\nabla T \]

• C_p: Özgül ısı kapasitesi, J/(kg·K),
• T: Sıcaklık, K veya °C farkı,
• q: İletimle ısı akısı, W/m²,
• Q: Kimyasal reaksiyonun hacimsel ısı kaynağı, W/m³,
• Q_p: Basınç işi,
• Q_vd: Viskoz yayılım terimidir.

Katı fazında hareket ve konveksiyon ihmal edildiğinde enerji denklemi:

\[ \rho C_p\frac{\partial T}{\partial t}+\nabla\cdot\mathbf{q} = Q+Q_{ted} \]

biçiminde verilmiştir. Termal genleşme terimi Q_ted çalışmada ihmal edilmiştir.

Hidratasyonun ekzotermik ısı kaynağı:

\[ Q=R_{react}\Delta H_{hydration} \]

şeklinde uygulanmıştır. Reaksiyon hızı yükseldikçe aynı bölgede açığa çıkan hacimsel ısı da artmaktadır.

Nem ve sıcaklığa bağlı reaksiyon kinetiği

Yerel reaksiyon hızı için Arrhenius biçimli ve buhar basıncı itici kuvveti içeren denklem kullanılmıştır:

\[ X_a=A_r\exp\left(-\frac{E_a}{RT}\right)(1-X) \left(1-\frac{p_v}{p_{eq}}\right) \]

• X_a: Yerel reaksiyon hızı, 1/s,
• A_r: Ön üstel katsayı; PDF’de 1,63×10⁴ olarak verilmiştir,
• E_a: Aktivasyon enerjisi, 56,6 kJ/mol,
• R: Evrensel gaz sabiti,
• T: Mutlak sıcaklık, K,
• X: Yerel dönüşüm,
• p_v: Su buharı kısmi basıncı,
• p_eq: Tuz-hidrat dengesi için buhar basıncıdır.

Denklemdeki işaretin hidratasyon yönündeki kullanımının uygulama ayrıntıları PDF’de tam olarak açıklanmamıştır. Metin, ısı kaynağının yalnızca su buharı basıncının denge eşiğini aştığı ve hava-tuz temasının bulunduğu bölgelerde etkin olduğunu belirtmektedir.

Denge buharı basıncı

Denge basıncı Clausius–Clapeyron türü bağıntıyla hesaplanmıştır:

\[ P_{eq}=P_{ref}\exp\left(-\frac{\Delta H^0}{RT} +\frac{\Delta S^0}{R}\right) \]

• P_ref: Referans basınç,
• ΔH⁰: Standart reaksiyon entalpisi; model tablosunda 67.400 J/mol,
• ΔS⁰: Standart reaksiyon entropisi; model tablosunda 175 J/(mol·K).

Bu ilişki, sıcaklık yükseldikçe tuzun su buharıyla denge koşulunun nasıl değiştiğini belirlemektedir. Reaktördeki gerçek buhar basıncı ile P_eq arasındaki fark, hidratasyonun sürmesini sağlayan temel termodinamik etkidir.

Gözenekli yatağın etkin ısıl özellikleri

Gaz ile katı fazın hacimsel ortalama ısıl özellikleri aşağıdaki şekilde birleştirilmiştir:

\[ (\rho C_p)_{eff} = (1-\varepsilon)(\rho C_p)_s +\varepsilon(\rho C_p)_v \]

\[ k_{eff} = (1-\varepsilon)k_s+\varepsilon k_v \]

Tuzun özellikleri dönüşüme göre monohidrat ve hekzahidrat arasında doğrusal olarak değiştirilmiştir:

\[ (\rho C_p)_s = (1-X)(\rho C_p)_{s1} +X(\rho C_p)_{s0} \]

\[ k_s=(1-X)k_{s1}+Xk_{s0} \]

PDF’de s1 ve s0 alt indislerinin sırasıyla hekzahidrat ve monohidratı gösterdiği belirtilmiştir. Ancak denklemlerde X’in hangi uç durumda hangi hidratı temsil ettiği yeterince açık değildir; bu tanım bağımsız yeniden üretimde netleştirilmelidir.

Havaya aktarılan yararlı ısı

Çıkış havasının taşıdığı yararlı ısı aşağıdaki duyulur ve su buharına bağlı entalpi farkından hesaplanmıştır:

\[ q_{heat} = \dot{m}_{a,in} \left[ (C_{pa}+w_{out}C_{pv})T_{out} - (C_{pa}+w_{in}C_{pv})T_{in} \right] \]

• \(\dot{m}_{a,in}\): Giriş kuru hava kütlesel debisi, kg/s,
• C_pa: Kuru havanın özgül ısısı,
• C_pv: Su buharının özgül ısısı,
• w_in, w_out: Giriş ve çıkış mutlak nemleri,
• T_in, T_out: Giriş ve çıkış sıcaklıklarıdır.

Sonuç grafiklerinde bu değer reaktör hacmine göre kW/m³ olarak raporlanmıştır.

Sınır koşulları

Giriş su buharı derişimi ve kısmi basıncı:

\[ c_{in}=\frac{P_v}{RT_{in}} \]

\[ P_v=P_s\cdot RH \]

bağıntılarıyla oluşturulmuştur. Burada P_s, giriş sıcaklığındaki doymuş su buharı basıncı; RH ise bağıl nemdir.

ALE hareketli ağ yöntemi

Hareketli yatağın aşağı yönlü ilerlemesi, geometrinin sabit bir katı akış alanı olarak bırakılması yerine ağ düğümlerinin katının yer değiştirmesini izlemesini sağlayan Arbitrary Lagrangian–Eulerian yöntemiyle modellenmiştir.

ALE uygulamasında:

• Katı çıkışında aşağı yönlü normal ağ hızı tanımlanmış,
• Yan duvarlarda teğetsel kısıtlaması olmayan kayma koşulları kullanılmış,
• Ağın istenmeyen bozulmasını önlemek için simetri veya makaralı kısıtlar uygulanmış,
• Ağ deformasyonunda Yeoh yöntemi ve 10 değerinde sertleştirme kullanılmıştır.

Bu hareketli ağ; akış, su buharı taşınımı, reaksiyon ve enerji denklemleriyle eş zamanlı çözülmüştür. Böylece yeni tuzun reaktöre girmesi ve kısmen hidratlanmış tuzun aşağı taşınması modellenebilmiştir.

Ağ bağımsızlığı

Araştırmacılar hesaplama maliyeti ile doğruluk arasındaki denge nedeniyle 644.783 elemanlı ağı kullanmıştır. Bununla birlikte tabloda sapmaların hangi ağ çiftine veya hangi nihai referansa göre hesaplandığı açık biçimde tanımlanmamıştır.

Metin üç boyutlu hesaplama bölgesinde “serbest üçgensel elemanlar” kullanıldığını söylemektedir. Üç boyutlu hacim ağlarında tetrahedral eleman ifadesi beklenir; PDF’deki ağ görseli de tetrahedral hacim ağına benzemektedir. Bu terminoloji yeniden üretim açısından açıklığa kavuşturulmalıdır.

Modelin deneysel verilerle karşılaştırılması

Model sıcaklık, giriş-çıkış nemi, yararlı ısı ve basınç kaybı açısından daha önce yayımlanmış hareketli yatak deneyleriyle karşılaştırılmıştır. Giriş neminin deneyde ilk yaklaşık 0,22 saat içinde 0’dan 9 g/kg’a yükselmesi, deney grafiğinden sayısallaştırılmış zamana bağlı bir interpolasyon fonksiyonuyla modele aktarılmıştır.

Sıcaklık karşılaştırma grafiğinde yatak içindeki iki ölçüm hattı ve çıkış sıcaklığı yaklaşık sekiz saat boyunca gösterilmektedir. Simülasyon sıcaklık eğrileri deneysel eğilimleri takip etmektedir. Çıkış sıcaklığı simülasyonda yaklaşık 32,7°C, deneyde ise yaklaşık 31°C çevresindedir. İç yatak sıcaklıkları yaklaşık 35–37°C aralığındadır.

Su buharı grafiğinde giriş nemi yaklaşık 0,009 kg/kg düzeyine, çıkış nemi ise yaklaşık 0,005 kg/kg düzeyine ulaşmaktadır. Çıkış neminin girişten düşük olması, su buharının tuz tarafından tüketildiğini göstermektedir.

Doğrulama tablosunun ayrıntılı değerlendirmesi

Tablodaki sayılardan hesaplanan göreli farklar yararlı ısı için yaklaşık yüzde 1,3–8,0, yatak sıcaklığı için yaklaşık yüzde 0,46–7,07 aralığındadır. Dolayısıyla model genel eğilim ve büyüklükleri yakalamakla birlikte, özet bölümündeki “± yüzde 1–4” ifadesi bütün doğrulama koşullarını kapsamamaktadır.

Basınç dağılımı ve yan akış yapısı

Akım çizgileri ve basınç dağılımı, hava girişinin dar bağlantı kanalında hızlandığını ve geniş yan hacme ulaştığında büyük dolaşım bölgeleri oluşturduğunu göstermektedir. Reaktif yatağın giriş tarafında basınç yaklaşık 1.600 Pa düzeyine çıkmakta ve gözenekli ortam boyunca sürekli azalmaktadır.

Model ile deney arasındaki toplam basınç kaybı yaklaşık 1.500–1.700 Pa aralığında uyumludur. Bu sonuç, yatak direncinin yanında hava yönlendirme duvarlarının ve genişleme-daralma geometrisinin de hidrolik performansta önemli olduğunu göstermektedir.

Akış çizgileri, yan hacimlerin tamamının eşit biçimde kullanılmadığını; bazı bölgelerde dolaşım ve kısa devre akışı oluşabildiğini göstermektedir. Bu nedenle çalışmada bypass faktörü ayrı bir performans değişkeni olarak ele alınmıştır.

Katı hızının etkisi

Katı hızı 0,5 cm/saatten 2 cm/saate yükseltildiğinde yatak sıcaklığı yaklaşık 35°C’den 35,5°C’ye çıkmıştır. Soğuk ve düşük hidratlı tuzun daha hızlı beslenmesi yerel reaksiyonu artırabilmekte, ancak aynı zamanda her bir katı kütlesinin yatakta kalma süresini azaltmaktadır.

Çıkış sıcaklığı ve havaya aktarılan yararlı ısı eğrileri dört hız için birbirine oldukça yakındır. Yararlı ısı başlangıçtaki hızlı yükselişten sonra yaklaşık 3 kW/m³ çevresinde kararlı hâle gelmektedir.

Düşük katı hızlarında tuz daha uzun süre reaktörde kaldığı için birim katı kütlesi başına daha ileri hidratasyon mümkündür. Yüksek hızlarda daha fazla katı geçirilebilse de her parçacığın dönüşümü daha düşük kalabilir. Bu durum, güç çıkışı ile toplam depolanan kimyasal enerjinin geri kazanılması arasında bir ödünleşme oluşturur.

Bununla birlikte PDF’de katı hızına karşı doğrudan toplam dönüşüm veya enerji verimi grafiği verilmemiştir. Bu nedenle daha yüksek hızın dönüşüm verimini ne ölçüde azalttığı nicel olarak gösterilmemektedir.

Reaksiyon cephesinin zamansal gelişimi

Başlangıçta bütün yatak yaklaşık 22°C’dir. Nemli hava girişi başladıktan sonra reaksiyon önce girişe yakın bölgede ortaya çıkmaktadır:

• t = 0 saat: Reaktör soğuk ve reaksiyon öncesi durumdadır.
• t = 0,2 saat: Giriş tarafında sarı-yeşil bir sıcaklık bölgesi oluşmaktadır.
• t = 0,3 saat: Reaksiyon bölgesi yatak boyunca genişlemekte ve sıcak çekirdek belirginleşmektedir.
• t = 1 saat: Reaktif yatağın önemli bölümü yaklaşık 34–38°C aralığına ulaşmaktadır.

Kontur görselleri sıcaklığın hava girişinden çıkış yönüne doğru yükseldiğini ve merkezdeki dikey reaktif yatak boyunca ilerlediğini göstermektedir. Reaksiyonun başladığı bölge, su buharı kısmi basıncı ile denge basıncı arasındaki farkın yeterince büyük olduğu bölgedir.

Metinde bu başlangıç durumu yanlışlıkla “dehidratasyon başlamadan önce” biçiminde tanımlanmıştır. İncelenen aşama deşarj ve hidratasyon olduğundan burada “hidratasyon başlamadan önce” ifadesi daha tutarlıdır.

Hava debisinin etkisi

Hava debisi 80, 110 ve 140 m³/saat için karşılaştırılmıştır. Sonuçlar tekdüze bir “daha fazla debi daha iyi” ilişkisi göstermemektedir.

80 m³/saat

Düşük debide hava tuzla daha uzun süre temas etmektedir. Hidratasyon için yeterli temas süresi sağlanmasına rağmen toplam hava kütlesel debisi düşük olduğu için reaksiyon ısısının reaktörden taşınması sınırlı kalmaktadır. Yararlı ısı yaklaşık 3 kW/m³ düzeyindedir.

110 m³/saat

Orta debi, yeterli temas süresi ile yeterli ısı taşıma kapasitesi arasında en iyi dengeyi oluşturmuştur. Yatak ve çıkış sıcaklığı en yüksek, yararlı ısı yaklaşık 4–4,5 kW/m³ ile üç koşulun en yükseğidir.

140 m³/saat

Yüksek debi reaktöre daha fazla hava taşımakta, ancak gazın yatakta kalma süresini azaltmaktadır. Güçlü konvektif soğutma yatağın daha soğuk ve sıcaklık dağılımının daha düzensiz olmasına neden olmaktadır. Çıkış sıcaklığı 110 m³/saat durumundan daha düşüktür.

Hız konturlarında hava giriş jetinin debiyle birlikte güçlendiği; 140 m³/saatte giriş bağlantısında yaklaşık 2,5 m/s’ye ulaşan yerel hızlar oluştuğu görülmektedir. Yüksek yerel hız her zaman daha yüksek dönüşüm anlamına gelmemektedir; çünkü reaksiyon için gaz-katı temas süresi de gereklidir.

Giriş neminin etkisi

Nem, çalışmanın en güçlü sonuç değişkenidir. Hava debisi 110 m³/saat ve katı hızı 0,5 cm/saat tutulduğunda:

Nem arttığında su buharı kısmi basıncı yükselmekte, denge basıncına göre hidratasyon itici kuvveti büyümekte ve daha fazla tuz aynı zaman aralığında reaksiyona girmektedir. Bunun sonucu olarak:

• Reaksiyon cephesi daha keskin hâle gelmekte,
• Yüksek sıcaklıklı bölge genişlemekte,
• Su buharı tüketimi artmakta,
• Çıkış havası daha sıcak hâle gelmekte,
• Havaya aktarılan ısı yükselmektedir.

Su buharı derişim alanı

Su buharı konturlarında 4 g/kg giriş neminde gaz derişimi yaklaşık 0,2–0,3 mol/m³ düzeyindedir. 9 g/kg durumunda yaklaşık 0,5–0,7 mol/m³, 12 g/kg durumunda ise yaklaşık 0,8 mol/m³’e ulaşmaktadır.

Giriş tarafındaki yüksek buhar derişimi, reaktif yatağa girdikçe tuz tarafından tüketilmektedir. Çıkışa doğru derişimin azalması, hidratasyon reaksiyonunun konumsal izini göstermektedir.

Yerel hidratasyon dönüşümü

Dönüşüm konturları X değerini 0–0,09 aralığında göstermektedir. En yüksek yerel dönüşüm yaklaşık yüzde 9’dur. Bu nedenle simülasyon süresinde tam hidratasyon elde edilmemiştir.

• 4 g/kg: Reaksiyon zayıf, dar ve dağınıktır.
• 9 g/kg: Reaksiyon cephesi belirginleşmekte ve yatak boyunca uzamaktadır.
• 12 g/kg: Daha fazla konumda yüksek yerel dönüşüm oluşmakta, ancak X yine 1’e yaklaşmamaktadır.

PDF, 9 g/kg nemde deneysel toplam dönüşümün yaklaşık yüzde 1,7–8,7 aralığında olduğunu ve sayısal sonuçların aynı büyüklük mertebesinde kaldığını belirtmektedir. Bu oldukça düşük dönüşüm, kararlı bir ısı gücü elde edilse bile katının kimyasal depolama kapasitesinin büyük bölümünün tek geçişte kullanılmadığını göstermektedir.

Buhar basıncı kontrollü rejim ne anlama gelmektedir?

Yatak sıcaklığı ve yararlı ısı katı hızına kıyasla giriş nemine çok daha güçlü tepki vermiştir. Araştırmacılar bu durumu “buhar basıncı kontrollü rejim” olarak tanımlamaktadır.

Bu rejimde reaksiyon hızını belirleyen ana etken, tuzun yatakta bulunma süresinden önce havanın gerçek su buharı basıncı ile denge buharı basıncı arasındaki farktır. Katı daha uzun süre yatakta kalsa bile su buharı itici kuvveti düşükse reaksiyon yavaş kalmaktadır.

Bununla birlikte tam dönüşüm için yalnızca yüksek nem yeterli değildir. Reaksiyon bölgesinin bütün katı hacmine ulaşabilmesi için yeterli zaman, uygun gaz dağılımı ve düşük bypass gerekmektedir.

Bypass faktörünün etkisi

Bypass faktörü, giriş havasının reaktif tuzla etkili temas kurmadan reaktörden geçebilen bölümünü temsil etmektedir. 110 m³/saat hava debisi, 22°C giriş sıcaklığı, 9 g/kg nem ve 0,5 cm/saat katı hızında:

Bypass arttıkça daha az hava tuzla reaksiyona girmekte, su buharı tüketimi ve reaksiyon ısısı azalmaktadır. Bu sonuç, yalnızca tuzun kimyasal özelliklerinin değil, giriş dağıtıcısı, duvar geometrisi ve yatak homojenliğinin de reaktör performansı için kritik olduğunu göstermektedir.

Tuzun ısıl iletkenliğinin etkisi

Tuzun ısıl iletkenliği 0,7–8 W/(m·K) arasında değiştirildiğinde yatak ve çıkış sıcaklığı eğrileri neredeyse üst üste gelmiştir. Araştırmacılar bunu reaktörün ince yatak yapısında ısı iletiminin, konvektif hava taşınımı ve reaksiyon ısısına göre ikincil kalmasıyla açıklamaktadır.

Bu sonuç incelenen geometri ve çalışma koşulları için geçerlidir. Daha kalın yataklarda, daha düşük hava hızlarında veya yerel sıcak noktaların bulunduğu farklı tasarımlarda tuz iletkenliğinin etkisi daha büyük olabilir.

Geçirgenliğin etkisi

Geçirgenlik 10⁻⁸–10⁻¹⁰ m² arasında değiştirildiğinde çıkış sıcaklığı üzerindeki etki sınırlıdır. Daha düşük geçirgenlikte hava akışı zorlaşmakta, fakat gazın reaktif katıyla temasının artması ve iç bypassın azalması yararlı ısıyı yükseltebilmektedir.

PDF’de 10⁻¹² m² koşulunda yararlı ısının belirgin biçimde yükseldiği görülmektedir. Ancak bu kadar düşük geçirgenlik aynı zamanda çok daha yüksek fan gücü ve basınç kaybı gerektirebilir. Parametrik bölümde basınç kaybı veya fan enerjisi hedef fonksiyona dâhil edilmediği için bu sonucun sistem düzeyindeki net yararı gösterilmemiştir.

Kinetik katsayının etkisi

Kinetik faktör yaklaşık 0,8×10⁴–2×10⁴ s⁻¹ aralığında artırıldığında çıkış sıcaklığı ve yararlı ısı yükselmiştir. Daha hızlı reaksiyon:

• Su buharının daha kısa mesafede tüketilmesine,
• Aktif reaksiyon bölgesinin genişlemesine,
• Daha hızlı ısı salımına,
• Çıkış havasının daha fazla ısınmasına

neden olmaktadır.

Ancak kinetik katsayının hangi malzeme modifikasyonu veya parçacık yapısıyla pratik olarak değiştirilebileceği modellenmemiştir. Sonuç, kinetiği geliştirilmiş kompozit veya yapılandırılmış tuzların potansiyelini göstermektedir.

Çalışmanın en önemli mühendislik bulguları

• Reaktör performansı katı hızından çok giriş havasının nemine duyarlıdır.
• En yüksek incelenen ısı çıkışı 12 g/kg giriş neminde elde edilmiştir.
• Pratik kış koşulu olarak önerilen 9 g/kg nem, ısı çıkışı ile ulaşılabilir nem arasında bir uzlaşma olarak sunulmuştur.
• Hava debisi için 110 m³/saat, temas süresi ile ısı taşıma kapasitesini en iyi dengeleyen koşuldur.
• Yüksek katı hızı sürekli güç çıkışını koruyabilse de tek geçiş dönüşümünü sınırlayabilir.
• Bypassın azaltılması, kimyasal kinetiği değiştirmeden reaktör performansını artırabilecek önemli bir tasarım stratejisidir.
• Reaksiyon cephesi hava giriş tarafında oluşmakta ve nem yükseldikçe yatak boyunca genişlemektedir.
• Simülasyon süresi içinde tam hidratasyon gerçekleşmemiş; en yüksek yerel dönüşüm yaklaşık yüzde 9’da kalmıştır.

Bina ısıtması açısından sonuçların anlamı

Çıkış havasının yaklaşık 33–35°C düzeyine ulaşması, düşük sıcaklıklı ısıtma sistemleri veya bir ısı pompasının kaynak sıcaklığını yükseltme gibi uygulamalar açısından yararlı olabilir. Ancak bu çalışma doğrudan bir bina ısıtma yükü, radyatör sistemi, fan gücü, ısı değiştirici performansı veya yıllık iklim profili simülasyonu gerçekleştirmemiştir.

Raporlanan kW/m³ değeri reaktör hacmi başına yararlı ısıdır. Gerçek bir binada sağlanabilecek toplam kW değeri, reaktif hacim, çalışma süresi, tuz miktarı, fan gücü, ısı değiştirici verimi ve çevrim stratejisine bağlı olacaktır.

Çalışmanın güçlü yönleri

• Hareketli yataklı açık termokimyasal reaktörü üç boyutlu olarak modellemektedir.
• Katı hareketini basit bir sabit kaynak terimi yerine ALE hareketli ağla temsil etmektedir.
• Gözenekli akış, reaksiyon, nem taşınımı ve ısı transferini eş zamanlı çözmektedir.
• Reaksiyon cephesini, su buharını, sıcaklığı ve dönüşümü uzaysal olarak göstermektedir.
• Model sonuçlarını sıcaklık, nem, yararlı ısı ve basınç kaybı açısından deneysel literatürle karşılaştırmaktadır.
• Nem, hava debisi, katı hızı, bypass, geçirgenlik, iletkenlik ve kinetik faktörü ayrı ayrı incelemektedir.
• Sürekli ısı gücü ile katı dönüşümü arasındaki farkı görünür hâle getirmektedir.
• Reaktör geometrisindeki dolaşım ve kısa devre akışlarını üç boyutlu akım çizgileriyle göstermektedir.

Çalışmanın temel sınırlılıkları

• Araştırma yeni deneysel veri veya fiziksel prototip sunmamaktadır.
• Model literatürdeki tek bir temel reaktör geometrisine dayanmaktadır.
• Tuzun uzun dönemli çevrim kararlılığı modellenmemiştir.
• Aglomerasyon, delikuesans, parçacık kırılması, aşınma ve mekanik sıkışma dikkate alınmamıştır.
• Gözeneklilik çoğunlukla sabit kabul edilmiş; reaksiyonla hacim değişimi ayrıntılı olarak çözülmemiştir.
• Katı hareketinin hava akışına geri etkisi ihmal edilmiştir.
• Parçacık boyutu dağılımı ve yatak çökelmesi raporlanmamıştır.
• Fan gücü ve pompalama enerjisi optimizasyona dâhil edilmemiştir.
• Çıkış ısısının haricî ısı değiştiricide ne kadarının kullanılabileceği hesaplanmamıştır.
• İklim ve mevsimsel nem değişimleri zaman serisi olarak modellenmemiştir.
• Yüksek nemde delikuesans veya akışkanlık kaybı incelenmemiştir.
• Katı hızına bağlı dönüşüm kaybı doğrudan bir dönüşüm grafiğiyle nicelendirilmemiştir.
• Simülasyon dosyaları, ağ, kod veya COMSOL modeli paylaşılmamıştır.
• Çözücü ayarları, zaman adımı, yakınsama ölçütleri ve hesaplama donanımı ayrıntılı olarak verilmemiştir.
• Çalışma hakem değerlendirmesinden geçmemiştir.

Çalışma neyi göstermektedir?

• Nemli hava ile çalışan SrBr₂ esaslı hareketli yatak reaktörünün üç boyutlu ve sürekli işletim altında sayısal olarak modellenebildiğini göstermektedir.
• Giriş neminin sıcaklık ve yararlı ısı üzerinde baskın bir etkiye sahip olduğunu göstermektedir.
• Orta hava debisinin düşük ve yüksek debiden daha yüksek yararlı ısı sağlayabileceğini göstermektedir.
• Reaksiyon cephesinin giriş tarafında oluştuğunu ve nemle genişlediğini göstermektedir.
• Bypassın azaltılmasının hava-tuz temasını ve ısı çıkışını artırabileceğini göstermektedir.
• Kararlı bir ısı çıkışının tam tuz dönüşümü gerçekleşmeden de elde edilebildiğini göstermektedir.
• ALE yönteminin hareketli katı yatağı, reaksiyon ve taşınım süreçleriyle birleştirebildiğini göstermektedir.

Çalışma neyi kanıtlamamaktadır?

• Reaktörün gerçek bir binada mevsimsel depolama ihtiyacını karşılayacağını kanıtlamamaktadır.
• 12 g/kg nem koşulunun bütün iklimlerde güvenli veya ekonomik olduğunu göstermemektedir.
• 9 g/kg değerinin bütün kış iklimleri için evrensel optimum olduğunu göstermemektedir.
• Daha düşük geçirgenliğin fan enerjisi hesaba katıldığında daha verimli olduğunu kanıtlamamaktadır.
• Stronsiyum bromürün binlerce çevrim boyunca yapısını koruyacağını göstermemektedir.
• Reaktörde topaklanma, delikuesans veya tıkanma oluşmayacağını kanıtlamamaktadır.
• Modelin farklı tuz hidratlarına veya farklı geometrilere doğrudan aktarılabileceğini göstermemektedir.
• Yararlı ısı değerlerinin net sistem verimi veya depolama çevrim verimi olduğunu göstermemektedir.
• Hakemli bilimsel uzlaşı veya ticari ürün doğrulaması sunmamaktadır.

Gelecekte yapılması gereken çalışmalar

Modelin bağımsız deneylerle doğrulanması için aynı reaktörde üç boyutlu sıcaklık sensörleri, yerel nem probları, basınç ölçümleri ve çıkış katısının hidratasyon analizi birlikte kullanılmalıdır. Reaksiyon cephesinin konumu, örneğin nötron görüntüleme, X-ışını yöntemleri veya farklı yüksekliklerden alınan katı numuneleriyle doğrulanabilir.

Uzun dönemli çevrim çalışmaları tuzun parçalanması, hacim değişimi, gözeneklilik kaybı, topaklanma ve delikuesans davranışını belirlemelidir. Bu etkiler modele geçirgenlik ve reaksiyon kinetiğinin çevrime bağlı değişimi olarak eklenebilir.

Reaktör optimizasyonunda yalnızca yararlı ısı değil; fan elektrik tüketimi, basınç kaybı, toplam dönüşüm, tuz verimi, ısı değiştirici performansı ve sermaye maliyeti de hedef olarak kullanılmalıdır.

Bypassı azaltmak için hava dağıtıcıları, delikli plakalar, yapılandırılmış reaktif malzemeler, kompozit tuzlar veya çok bölgeli yatak tasarımları incelenebilir. Aktif öğrenme veya vekil modellerle üç boyutlu tasarım uzayı daha düşük hesaplama maliyetiyle taranabilir.

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

Teknik yöntem özeti

Temel parametrik sonuçlar

Nem koşullarının karşılaştırılması

Hava debisi ödünleşmesi

Doğrulama verilerinden hesaplanan göreli farklar

Bu yüzdeler PDF’deki deney ve simülasyon değerlerinden hesaplanmıştır. Model genel eğilimleri yakalasa da doğruluk bütün koşullarda yüzde 1–4 aralığında değildir.

Reaktör tasarımı için çıkarılabilecek kararlar

1. İşletme nemi, reaktörün sıcaklık ve güç kontrolünde birincil ayar değişkenidir.
2. Hava debisi yalnızca maksimuma çıkarılmamalı; gaz-katı temas süresine göre optimize edilmelidir.
3. Katı hızı güç talebine göre ayarlanırken tek geçiş dönüşümü ayrıca izlenmelidir.
4. Hava dağıtım geometrisi bypassı ve basınç kaybını azaltacak şekilde geliştirilmelidir.
5. Düşük geçirgenlikle elde edilen termal kazanç, fan elektrik tüketimiyle birlikte değerlendirilmelidir.
6. Tam enerji depolama verimi için yalnızca çıkış ısısı değil, çıkış katısının hidratasyon derecesi de ölçülmelidir.

Kaynak ve Yöntem Notu

Çalışmanın özgün adı: 3D Multiphysics Modeling of a Moving-Bed Thermochemical Reactor: Performance Analysis and Optimization under Humid Air Conditions

Yazarlar ve sıraları: Sanaz Akbarzadeh; Abhishek K. Singh; Haoshui Yu; Ahmad Arabkoohsar.

Eş katkı bilgisi: PDF’de eş katkı veya eş birinci yazarlık beyanı bulunmamaktadır.

Sorumlu yazar: Ahmad Arabkoohsar.

Sorumlu yazar e-postası: ahmar@dtu.dk

Kurumlar:

1. Department of Civil and Mechanical Engineering, Technical University of Denmark, Kgs. Lyngby, Denmark.
2. Department of Thermal and Fluid Engineering, University of Twente, the Netherlands.
3. Department of Energy, Aalborg University, DK-9220 Aalborg East, Denmark.

Yazar-kurum eşleştirmesi: Sanaz Akbarzadeh ve Ahmad Arabkoohsar Technical University of Denmark; Abhishek K. Singh University of Twente; Haoshui Yu Aalborg University.

ORCID bilgileri: PDF’de yazarların ORCID numaraları verilmemiştir.

DOI: 10.2139/ssrn.6945146

Kaynak türü: Üç boyutlu sayısal modelleme ve parametrik reaktör performans analizi preprinti.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir.

Yayın platformu: SSRN.

SSRN gönderim tarihi: 15 Haziran 2026.

Sayfa sayısı: 46. Ana metin ve kaynakça ilk 24 sayfada, büyütülmüş şekiller 25–46. sayfalarda yer almaktadır.

Dergi: Hakemli bir dergi adı veya kabul edilmiş nihai yayın bilgisi bulunmamaktadır.

Özgün dergi yayınevi: Bu sürüm için doğrulanamamıştır. SSRN, Elsevier tarafından işletilen açık erişimli bir erken araştırma ve preprint platformudur; çalışma hakemli bir Elsevier dergi makalesi olarak sunulmamalıdır.

Resmî bağlantılar:SSRN kayıt sayfası ve DOI bağlantısı.

Finansman: Çalışma Avrupa Birliği tarafından desteklenmiştir. PDF’de “A Game-Changing Solar-Thermochemical Heat Storage and Intensification Technology for the Process Heat Sector” başlıklı, proje numarası 101147098 olan Marie Skłodowska-Curie Action / Horizon Europe projesinin parçası olduğu belirtilmiştir.

Çıkar çatışması: PDF’de ayrı bir çıkar çatışması beyanı bulunmamaktadır.

Yazar katkıları: PDF’de CRediT veya benzeri ayrıntılı yazar katkısı beyanı bulunmamaktadır.

Veri ve kod erişimi: COMSOL model dosyası, geometrik model, ağ dosyası, çözücü ayarları veya analiz kodu paylaşılmamıştır. Doğrulamada kullanılan deneysel veriler daha önce yayımlanmış kaynaklardan ve bir grafikten sayısallaştırılan giriş nemi verisinden alınmıştır.

Yapay zekâ kullanım beyanı: Yazarlar metnin dil bilgisi ve yazımını iyileştirmek için Microsoft Copilot kullandıklarını; içeriği daha sonra gözden geçirip düzenlediklerini ve bütün sorumluluğu üstlendiklerini belirtmiştir.

Bu Türkçe makale, yüklenen 46 sayfalık PDF’nin metni, formülleri, tabloları, reaktör şemaları, ağ görseli, sıcaklık-nem-dönüşüm konturları, basınç dağılımı, akım çizgileri ve parametrik grafiklerinin tamamı incelenerek hazırlanmıştır. PDF’de bulunmayan deneysel başarı, ekonomik fizibilite, bina enerji tasarrufu, uzun dönemli malzeme ömrü veya ticari uygulanabilirlik iddiası eklenmemiştir. Dış kaynaklar yalnızca başlık, yazarlar, DOI, sorumlu yazar, SSRN kayıt tarihi ve platform statüsünün bibliyografik doğrulanması için kullanılmıştır.

PDF’de belirlenen iç tutarsızlıklar ve yeniden üretilebilirlik sorunları

• Hava debisi aralığı: Yöntem bölümünde giriş debisinin 25–100 m³/saat olduğu yazılmasına rağmen sonuçlarda 110 ve 140 m³/saat koşulları kullanılmıştır. İncelenen gerçek parametrik aralık açıkça 80–140 m³/saattir.
• Başlangıç hidrat derecesi: Tuz “monohidrat” olarak tanımlanırken SrBr₂·XH₂O, X≈0 ifadesi kullanılmıştır. Monohidrat için X=1 olması beklenir; X≈0 susuz veya susuza yakın durumu ifade eder.
• Hidrat alt indisleri: Etkin özellik denklemlerinde s1 ve s0 sırasıyla hekzahidrat ve monohidrat olarak açıklanmış, fakat X’in uç değerlerinin hangi faza karşılık geldiği net değildir.
• Doğrulama kaynağı: Metin bir yerde deneysel doğrulamayı Kaynak [15] ile ilişkilendirmektedir. Kaynak [15] sayısal çalışmadır; kullanılan deneysel reaktör verileri şekil ve tablolara göre Kaynak [28]’den gelmektedir.
• Doğrulama hata aralığı: Özet ± yüzde 1–4 sapma bildirmektedir. Tablo 4 değerlerinden hesaplanan göreli farklar yararlı ısıda yaklaşık yüzde 8’e, yatak sıcaklığında yüzde 7,07’ye kadar çıkmaktadır.
• Özgül ısı birimi: Tablo 3’te tuzun özgül ısısı 970 J/(g·K) olarak yazılmıştır. Bu değer fiziksel olarak olağan dışıdır ve büyük olasılıkla 970 J/(kg·K) olmalıdır.
• Tablo numaralandırması: Hem reaktör ölçü tablosu hem ağ bağımsızlığı tablosu “Tablo 2” olarak adlandırılmıştır.
• Ağ elemanı terminolojisi: Üç boyutlu hacim modeli için “serbest üçgensel eleman” ifadesi kullanılmıştır. Görsel ve COMSOL terminolojisi tetrahedral hacim ağını düşündürmektedir.
• Ağ sapmalarının referansı: Ağ bağımsızlığı tablosundaki yüzdelerin hangi ağ çözümüne göre hesaplandığı açıklanmamıştır.
• Reaksiyon yönü: Hidratasyon denklemi ve reaksiyon hızı ifadesindeki işaretlerin COMSOL uygulamasında nasıl sınırlandırıldığı yeterince gösterilmemiştir.
• Bypass faktörünün denklemi: F_b tanımlanmakta ve parametrik olarak değiştirilmektedir; fakat gösterilen reaksiyon hızı veya taşınım denklemlerinde nasıl uygulandığı açıkça yazılmamıştır.
• Dehidratasyon/hidratasyon ifadesi: Deşarj aşamasındaki Şekil 7 açıklaması başlangıcı yanlışlıkla “dehidratasyon başlamadan önce” olarak adlandırmaktadır. İncelenen süreç hidratasyondur.
• Hacimsel ısı birimleri: Geçirgenlik ve kinetik faktör grafiklerinin sağ eksenlerinde yaklaşık 1.000–4.000 değerleri “kW/m³” olarak etiketlenmiştir. Diğer bütün sonuçlar 2–5,5 kW/m³ aralığındadır. Bu eksenlerin büyük olasılıkla W/m³ olması veya ölçek katsayısı içermesi gerekir.
• Katı hızı ve dönüşüm: Sonuç bölümü yüksek katı hızının dönüşüm verimini azalttığını belirtmektedir; ancak katı hızına karşı doğrudan dönüşüm verisi veya enerji verimi grafiği sunulmamıştır.
• Pratik nem optimumu: 9 g/kg değeri pratik kış koşulu olarak önerilmiştir. Farklı iklimlerin nem profili, delikuesans sınırı veya iklimsel erişilebilirlik analizi yapılmadığından bu değer evrensel optimum olarak yorumlanmamalıdır.
• Eksik hesaplama ayrıntıları: Zaman adımı, doğrusal olmayan çözücü, yakınsama eşikleri, stabilizasyon yöntemleri, donanım ve toplam hesaplama süresi belirtilmemiştir.

Bu tutarsızlıklar, giriş neminin ve hava dağılımının reaktör performansı üzerindeki güçlü etkisini gösteren genel sonuçları ortadan kaldırmamaktadır. Ancak modelin bağımsız olarak yeniden üretilmesi, hata sınırlarının doğru değerlendirilmesi ve tasarım önerilerinin gerçek bir sisteme aktarılması için bu noktaların düzeltilmesi gerekir.