Bu çalışma, türbülanslı bir kanal içindeki suyun duvardan üstel biçimde artan ısı akısına verdiği tek fazlı tepkiyi, çekirdekli kaynamanın başlangıcına kadar yüksek çözünürlüklü büyük girdap benzetimleriyle incelemektedir. Araştırmanın temel amacı, hızlı bir güç artışı sırasında termal sınır tabakanın yalnızca sıcaklık bakımından değil, hız alanı, girdap yapıları, viskoz gerilmeler ve yönsel ısı taşınımı bakımından nasıl yeniden örgütlendiğini ortaya koymaktır.

Simülasyonlarda sıcaklık pasif bir izleyici olarak değil, suyun dinamik viskozitesini ve ısıl iletkenliğini değiştiren aktif bir skaler olarak ele alınmıştır. Böylece sıcaklık alanı akış özelliklerini, akış alanı da sıcaklık dağılımını eşzamanlı olarak etkilemektedir. Hesaplanan duvar sıcaklıkları, kaynama başlamadan önceki deneysel ölçümlerle genel olarak yüzde 5’in altında farkla uyuşmuştur.

Çalışmanın sonuçlarına göre hızlı ısıtma, sıcak duvar yakınındaki viskoziteyi düşürerek akış yönündeki hızı artırmakta, açıklık doğrultusundaki vortisiteyi güçlendirmekte ve tampon tabakada yaklaşık [ y^+=60\text{–}140 ] arasında Couette akışını andıran, neredeyse doğrusal bir hız profili oluşturmaktadır. Isıtılmış akış yapıları duvara dik yönde sıkışırken kanal genişliği doğrultusunda yayılmaktadır. Bu değişiklikler, ısının duvara yakın bölgede tutulmasına katkıda bulunmaktadır.

Hızlı geçişte Stanton sayısı yaklaşık yüzde 30 yükselmekte ve termal olarak tam gelişmiş bölge kanalın daha erken bir konumunda oluşmaktadır. Buna karşılık yavaş, yarı kararlı ısıtmada Stanton sayısı geçiş boyunca yaklaşık yüzde 8 içinde kalmakta ve termal giriş bölgesi daha uzun sürmektedir.

Kaynama başlangıcına hemen önce hızlı ısıtma, duvar boyunca geniş ölçüde homojen bir aşırı ısıtılmış sıvı tabakası meydana getirmektedir. Yarı kararlı ısıtmada ise aşırı ısıtılmış sıvı, türbülanslı düşük hızlı çizgilerle ilişkili akış yönündeki şeritler içinde kalmaktadır. Böylece ısıtma süresi yalnızca toplam ısı aktarımını değil, ilk kabarcıkların nerede ve hangi uzamsal düzen içinde ortaya çıkabileceğini de belirlemektedir.

Araştırma, kabarcık oluşumunu veya iki fazlı kaynamayı doğrudan simüle etmemektedir. Bütün hesaplamalar, deneylerde ilk çekirdekli kaynama belirtisinin görüldüğü [ t^*=0 ] anına kadar olan tek fazlı sıvı rejimiyle sınırlıdır.

Araştırmanın ele aldığı problem

Duvara aktarılan ısının milisaniye ölçeğinde hızla artması; nükleer reaktörler, yüksek güçlü elektronik sistemler, roket motorları ve gelecekteki füzyon reaktörlerinin plazmaya bakan bileşenleri açısından kritik bir ısıl güvenlik problemidir. Isı yükü yeterince hızlı yükseldiğinde duvarın yakınındaki sıvı, kanalın geri kalanı hâlâ güçlü biçimde soğutulmuş olsa bile doyma sıcaklığının üzerine çıkabilir.

Bu koşullarda sistem üç farklı yola ilerleyebilir:

• Soğutma tek fazlı zorlanmış taşınım rejiminde kalabilir.
• Çekirdekli kaynama başlayarak ısı aktarımı geçici olarak güçlenebilir.
• Daha ağır koşullarda buhar filmi oluşabilir ve duvarın sıvıyla teması azalarak ısı uzaklaştırma kapasitesi tehlikeli biçimde düşebilir.

Özellikle araştırma reaktörlerinde bir reaktivite eklenmesi kazası, yakıt gücünde yaklaşık üstel bir artışa yol açabilir. Soğutucunun tek fazlı sınır tabakasının bu hızlı artışa nasıl tepki verdiği, yakıt kılıfının ulaştığı sıcaklığı ve kaynamanın ne zaman başlayacağını etkiler.

Önceki çalışmalar çoğunlukla ortalama duvar sıcaklığı, ortalama ısı aktarım katsayısı veya toplam ısı akısı gibi bütünleşik ölçümlere dayanmıştır. Bu ölçümler duvarın iç tarafında gelişen üç boyutlu hız, vortisite, sıcaklık ve gerilme alanlarını doğrudan göstermez. Araştırma, bu eksikliği yüksek çözünürlüklü büyük girdap benzetimleriyle gidermeyi amaçlamaktadır.

Çalışmanın araştırma soruları

Çalışmada aşağıdaki temel sorular ele alınmaktadır:

• Üstel biçimde artan duvar ısı akısı, türbülanslı kanal akışının yakın duvar yapısını nasıl değiştirmektedir?
• Sıcaklığa bağlı viskozite ve ısıl iletkenlik, momentum ve ısı aktarımı arasında nasıl bir çift yönlü bağ oluşturmaktadır?
• Hızlı bir geçiş ile yarı kararlı bir geçiş arasında termal giriş uzunluğu, Stanton sayısı ve sıcaklık dağılımı bakımından hangi farklar vardır?
• Kaynama başlamadan hemen önce aşırı ısıtılmış sıvı homojen bir film mi, yoksa türbülanslı şeritler hâlinde mi oluşmaktadır?
• Duvara dik, akış yönündeki ve kanal genişliği yönündeki ısı akıları hangi büyüklük sırasını izlemektedir?
• Farklı zamanlarda ve kanal geometrilerinde sıcaklık profilleri ortak bir matematiksel bağıntıyla temsil edilebilir mi?

Deneysel sistem ve işletme koşulları

Sayısal model, deiyonize su kullanılan azaltılmış ölçekli bir akış çevriminin deneysel koşullarına dayanmaktadır. İki ince dikdörtgen kanal geometrisi incelenmiştir:

İşletme basıncı 10 bar, giriş sıcaklığı 353 K yani 80 °C olarak tutulmuştur. Bu koşullarda doyma sıcaklığı yaklaşık 453 K olduğundan giriş soğutma derecesi:

\[ \Delta T_{\mathrm{sub}}=T_{\mathrm{sat}}-T_{\mathrm{in}}=100\ \mathrm{K} \]

olarak tanımlanmıştır.

Isıtıcı eleman 20 mm uzunluğunda ve 4 mm genişliğindedir. Safir bir altlık üzerine buhar biriktirme yöntemiyle uygulanmış, nanometrik pürüzlülüğe sahip titanyum film kullanılmıştır. Safirin kızılötesi ışınıma kısmen geçirgen olması, ısıtıcı yüzey sıcaklığının alt taraftan kızılötesi termografiyle ölçülmesini sağlamıştır.

Deneysel Joule gücü ile sıvıya geçen net duvar ısı akısı aynı değildir. Safir altlık belirli miktarda enerji depoladığı için özellikle [ \tau=5\ \mathrm{ms} ] geçişinde sıvıya aktarılan net ısı akısı kusursuz bir üstel eğri izlememektedir. Sayısal sınır koşulu olarak idealize edilmiş elektrik gücü değil, deneyden elde edilen zaman bağımlı net ısı akısı kullanılmıştır.

Üstel ısı akısı ve zaman değişkenleri

Duvara uygulanan ısı akısının genel davranışı:

\[ q''_w(t)\propto \exp\left(\frac{t}{\tau}\right) \]

biçimindedir. Burada [ \tau ], güç artışının karakteristik zamanıdır. Küçük [ \tau ] değeri daha hızlı, büyük [ \tau ] değeri daha yavaş bir güç artışı anlamına gelir.

Karşılaştırmalar için duvar sıcaklığı ve zaman boyutsuzlaştırılmıştır:

\[ T_w^*=\frac{T_w-T_{\mathrm{sat}}}{\Delta T_{\mathrm{sub}}} \]

\[ t^*=\frac{t-t_{\mathrm{ONB}}}{\tau} \]

\[ t_{\mathrm{ONB}} ], deneylerde çekirdekli kaynamanın başladığı andır. Buna göre:

• [ t^*<0 ] tek fazlı kaynama öncesi dönemi,
• [ t^*=0 ] deneysel kaynama başlangıcını,
• [ t^*>0 ] ise simülasyonun fiziksel kapsamı dışında kalan iki fazlı dönemi temsil etmektedir.

Sayısal düzenek

Şekil 2.1’de akış yönünde iki ardışık hesap bölgesi gösterilmektedir. İlk bölüm, tam gelişmiş türbülans üretmek için kullanılan hidrolik geri dolaşım kutusudur. Bu kutunun çıkışı her zaman adımında ısıtılmış test kesitinin girişine aktarılmaktadır.

Test kesitinde:

• Alt duvara zamana bağlı ve uzayda eşit ısı akısı uygulanmaktadır.
• Üst duvar adyabatiktir.
• Her iki duvarda kaymama koşulu uygulanmaktadır.
• Kanal genişliği doğrultusunda periyodik sınır koşulu kullanılmaktadır.
• Çıkış 10 bar işletme basıncında tutulmaktadır.
• Giriş akışkanı 353 K sıcaklıktadır.

Öncül simülasyonlar, birinci dereceden türbülans istatistikleri kararlı hâle gelene kadar 40–100 yıkama zamanı boyunca çalıştırılmıştır:

[ \tau_{\mathrm{wash}}=\frac{L_x}{U} \]

Burada [ L_x ] hesap alanının akış yönündeki uzunluğu, [ U ] ortalama akış hızıdır.

Çözülen korunum denklemleri

Akış sıkıştırılamaz kabul edilmiştir. Süreklilik denklemi:

\[ \nabla\cdot\mathbf{u}=0 \]

Momentum denklemi:

\[ \frac{\partial\mathbf{u}}{\partial t} +\nabla\cdot(\mathbf{u}\mathbf{u}) = -\frac{1}{\rho}\nabla P +\nabla\cdot\left[\nu(T)\nabla\mathbf{u}\right] -\nabla\cdot\boldsymbol{\tau}^{\mathrm{SGS}} \]

Sıcaklık denklemi ise:

\[ \frac{\partial T}{\partial t} +\mathbf{u}\cdot\nabla T = \frac{1}{\rho c_p} \nabla\cdot\left[k(T)\nabla T\right] + \nabla\cdot\left(\alpha_t\nabla T\right) \]

biçiminde çözülmüştür.

\[ \nu(T) ] sıcaklığa bağlı kinematik viskoziteyi, [ k(T) ] sıcaklığa bağlı ısıl iletkenliği, [ \boldsymbol{\tau}^{\mathrm{SGS}} ] çözülemeyen küçük ölçeklerin gerilmesini ve [ \alpha_t ] alt ağ ölçeği ısıl yayınımını göstermektedir.

Yoğunluk [ \rho ] ve özgül ısı [ c_p ] sabit kabul edilmiştir. Yoğunluk giriş sıcaklığında, özgül ısı ise doyma sıcaklığında değerlendirilmiştir. Bu yaklaşım, sıcaklığa bağlı bütün termofiziksel özellikleri içeren tam sıkıştırılabilir bir model değildir.

Yerçekimi etkisi ihmal edilmiştir. Çalışmada hesaplanan toplu Richardson sayısı:

[ Ri=\frac{(\Delta\rho/\rho)gL}{U^2} \]

için [ Ri<0{,}01 ] olduğundan zorlanmış taşınımın kaldırma kuvvetine baskın olduğu kabul edilmiştir.

Sıcaklığa bağlı su özellikleri

Su özellikleri CATHARE termal-hidrolik özellik kütüphanesinden alınmıştır. Şekil 2.2, 10 bar basınçta viskozite, ısıl iletkenlik ve Prandtl sayısının sıcaklıkla değişimini göstermektedir.

Dinamik viskozite

Dinamik viskozite sıcaklık arttıkça yaklaşık üstel biçimde azalmaktadır. Simülasyonda kapsanan sıcaklık aralığında viskozite yaklaşık üç kat düşmektedir. Bu düşüş duvar yakınındaki momentum yayınımını ve viskoz gerilmeleri önemli ölçüde değiştirmektedir.

Isıl iletkenlik

Isıl iletkenlik monoton değildir. Yaklaşık yüzde 10’luk bir aralıkta değişmekte ve 416 K civarında en yüksek değerine ulaşmaktadır. Duvara en yakın aşırı ısıtılmış filmde iletkenlik yaklaşık yüzde 5 azalmaktadır. Isı duvardan sıvıya ilk olarak iletimle geçtiği için bu görece küçük azalma bile geçişin son aşamasında önem taşımaktadır.

Prandtl sayısı

Moleküler Prandtl sayısı:

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

olarak tanımlanmıştır. [ T>443\ \mathrm{K} ] olduğunda [ Pr<1 ] düzeyine düşmektedir. Böylece kanalın soğuk çekirdeğinde en küçük termal ölçek Batchelor ölçeğiyle, sıcak duvar yakınında ise Kolmogorov ölçeğiyle ilişkilendirilmektedir.

Şekil 3.4’te ısıl iletkenlik, viskozite ve Prandtl sayısının duvara dik–kanal genişliği kesitleri gösterilmektedir. Viskozitedeki değişiklikler [ y^+>200 ] bölgesine kadar uzanırken, iletkenlikteki en belirgin değişiklikler daha ince bir duvar yakın bölgesinde kalmaktadır.

Büyük girdap benzetimi ve hesaplama çözünürlüğü

Hesaplamalar TRUST/TrioCFD yazılımıyla ve Fransız TGCC süper bilgisayar kaynaklarında gerçekleştirilmiştir. En küçük durumlarda 10–40 milyon, en büyük durumda ise iki ayrı hesap bölümünde yaklaşık 650’şer milyon hücre kullanılmıştır.

Başlıca sayısal yöntemler şunlardır:

Duvara yakın ağ aralıkları yaklaşık [ \Delta y_w^+=0{,}5 ], akış yönündeki aralık [ \Delta x^+=8 ] ve kanal genişliği yönündeki aralık [ \Delta z^+=4 ] olarak verilmiştir. Fiziksel duvar-normal hücre yüksekliği koşula bağlı olarak yaklaşık 0,9–2,13 µm arasındadır.

Küçük ölçekli momentum aktarımı için WALE modeli ve [ C_w=0{,}5 ] kullanılmıştır. Termal alt ağ modeli için türbülanslı Prandtl sayısı:

\[ Pr_t=\frac{\nu_t}{\alpha_t}=0{,}9 \]

olarak seçilmiştir.

Ek B’deki kalite değerlendirmesinde alt ağ viskozitesinin toplam etkideki payı yüzde 10’un, modellenmiş termal akının toplam termal aktarım içindeki payı ise yüzde 1,2’nin altında kalmıştır. Bu sonuçlar, incelenen yakın duvar mekanizmalarının büyük bölümünün doğrudan çözüldüğünü göstermektedir. Bununla birlikte hesaplamalar doğrudan sayısal simülasyon değil, LES’tir.

Duvar sıcaklığı doğrulaması

Şekil 3.1’de [ \tau=5 ], 20 ve 100 ms için boyutsuz duvar sıcaklıkları deneylerle karşılaştırılmıştır. Kaynama başlamadan önce simülasyon ve deney eğrileri genel olarak yakın seyretmektedir.

\[ t^*<-1 ] döneminde en büyük fark hızlı geçişlerde görülmektedir. Araştırmacılar bunu, standart girdap-viskozitesi LES modellerinin türbülans üretimi ile yutulması arasında yerel denge varsaymasına bağlamaktadır. Çok hızlı ve denge dışı akışlarda bu varsayım zayıflayabilir.

[ -1<t^*<0 ] aralığında fark azalmaktadır. Kaynama başlangıcının hemen öncesindeki bazı zaman adımlarında deneysel ve sayısal sıcaklıklar neredeyse çakışmaktadır.

Tek fazlı dönemde hesaplanan duvar sıcaklıkları çoğunlukla deneysel değerlerin biraz üzerindedir. Bu durum simülasyonun ısı aktarım katsayısını hafifçe düşük hesapladığını göstermektedir. Fark çoğunlukla yüzde 5’in altında kalmıştır.

[ t^*>0 ] olduğunda deneyde kabarcık oluşurken model tek fazlı kaldığı için hata hızla yükselmektedir. Bu nedenle çalışma [ t^*=0 ] sonrasına ilişkin fiziksel sonuç üretmemektedir.

Kaynama başlangıcındaki sıcaklık aşımı

Hızlı ısı akısı artışları, çekirdekli kaynama başlamadan önce duvarın doyma sıcaklığının daha fazla üzerine çıkmasına izin vermektedir. Doğrulama eğrilerinde yarı kararlı geçiş için yaklaşık 20 K, hızlı geçiş için yaklaşık 65 K duvar aşırı sıcaklığı bildirilmektedir.

Akışın termal olarak gelişmiş aşağı akış bölgesindeki ortalama duvar-normal profillerde ise [ t^*=0 ] anında yaklaşık:

[ T^*\approx0{,}75\quad(\tau=5\ \mathrm{ms}) \]

\[ T^*\approx0{,}30\quad(\tau=100\ \mathrm{ms}) \]

değerleri görülmektedir. [ \Delta T_{\mathrm{sub}}=100\ \mathrm{K} ] olduğundan bunlar sırasıyla yaklaşık 75 K ve 30 K yerel/alt bölge ortalamalı aşırı sıcaklıklara karşılık gelir. Bu değerler, yüzeyin tamamına ait doğrulama ortalamalarıyla aynı uzaysal ortalama değildir; dolayısıyla doğrudan birbirinin yerine kullanılmamalıdır.

Anlık sıcaklık alanları ve türbülanslı çizgiler

Şekil 3.2, [ y^+=13 ] düzlemindeki sıcaklık ve yerel sıcaklık değişim hızı alanlarını göstermektedir. Sıcak bölgeler akış yönünde uzayan ince şeritler hâlindedir. Şeritlerin kanal genişliği yönündeki aralığı yaklaşık 60–120 viskoz birimdir. Bu aralık, yakın duvar düşük hızlı türbülans çizgilerinin bilinen karakteristik aralığıyla uyumludur.

Girişe yakın ilk yaklaşık [ 1000\delta_\nu ] mesafede ısı büyük ölçüde viskoz alt tabakada kalmaktadır. Daha aşağı akışta ısı tampon tabakaya ulaşmakta, sıcak akışkan duvardan dışarı fırlatılmakta ve soğuk çekirdek akışından gelen büyük girdaplarla karışmaktadır.

Sıcaklık değişim hızı alanında yan yana bulunan ısınma ve soğuma bölgeleri görülmektedir. Türbülanslı girdapların sıcak ve soğuk sıvı arasındaki arayüzü katlayıp uzatması, termal gradyanların ve yerel enerji aktarımının artmasına yol açmaktadır.

Girdap yapılarının değişimi

Şekil 3.3’te girdap çekirdekleri [ \lambda_2 ] ölçütüyle tanımlanmıştır. Geçişin erken döneminde farklı büyüklüklerde yoğun ve duvardan kanal çekirdeğine kadar uzanan vorteks yapıları görülmektedir.

\[ t^*=0 ] anında hem hızlı hem yarı kararlı ısıtmada:

• Duvara bağlı vortekslerin yüksekliği azalmakta,
• Vortekslerin yoğunluğu ve uzunluğu düşmekte,
• Yapılar daha ince bir yakın duvar bölgesine sıkışmaktadır.

Çalışma bunu sıcaklık arttıkça viskozitenin azalması ve yakın duvar türbülans kinetik enerjisinin yeniden düzenlenmesiyle ilişkilendirmektedir.

Hızlı geçişte [ t^*=0 ] anında ısıtılan yüzeyin hemen üzerindeki sıvının neredeyse tamamı doyma sıcaklığını aşmaktadır. Yarı kararlı geçişte ise aşırı ısıtılmış alanlar, düşük hızlı çizgiler boyunca uzanan akış yönlü bantlar hâlindedir ve aralarında daha soğuk sıvı bulunmaktadır.

Bu görüntüler iki farklı kaynama başlangıcı yolu önermektedir:

• Hızlı geçiş: Yüzeye yayılmış, uzamsal olarak daha eş zamanlı kaynama başlangıcı.
• Yarı kararlı geçiş: Türbülanslı çizgilere bağlı, şeritler hâlinde ve önce aşağı akışta başlayan kaynama.

Yakın duvar hızının artması

Boyutsuz hız:

[ U^+=\frac{U_x}{u_\tau} \]

ve duvar koordinatı:

\[ y^+=\frac{y}{\delta_\nu} \]

üzerinden incelenmiştir. Klasik izotermal viskoz alt tabakada yaklaşık [ U^+=y^+ ] doğrusal ilişkisi beklenir.

Şekil 3.5’e göre [ t^*=0 ] anında yakın duvar hız profilinin eğimi bu kanonik değere göre yaklaşık yüzde 30 artmaktadır. Aynı bölgede boyutsuz açıklık vortisitesinin büyüklüğü yaklaşık iki katına çıkmaktadır. Etki [ \tau=5\ \mathrm{ms} ] geçişinde daha belirgindir.

Vortisite ile hız gradyanı arasındaki ilişki nedeniyle bu bulgu, sıcak duvar yakınında daha güçlü bir kesme tabakası oluştuğunu göstermektedir.

Couette benzeri ara tabaka

\[ y^+=60\text{–}140 ] aralığında hız profili neredeyse kusursuz doğrusal hâle gelmektedir. Doğrusal uyum için bildirilen belirleme katsayısı:

[ R^2>0{,}99 \]

düzeyindedir. Aynı aralıkta açıklık vortisitesi yaklaşık sabit bir plato oluşturmaktadır.

Sabit hız gradyanı ve yaklaşık sabit vortisite, sabit kesmeli Couette akışına benzeyen bir ara tabakayı göstermektedir. Bu tabaka klasik iki plaka arasındaki bağımsız bir Couette akışı değildir; ısıtma ve viskozite değişimi nedeniyle türbülanslı kanalın tampon bölgesinde ortaya çıkan yerel bir benzerliktir.

Viskoz kayma gerilmesi

LES çerçevesindeki toplam viskoz gerilme:

\[ \tau_{ij} = \left[\mu(T)+\mu_t\right] \left( \frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i} \right) \]

olarak tanımlanmıştır.

Duvar yakınındaki ana kayma gerilmesi [ \tau_{xy} ], viskozitenin düşmesi nedeniyle geçiş ilerledikçe azalmaktadır. Kanal merkezinde işaret değiştirmeye devam etmekte ve sıfır noktası merkez çizgisinde kalmaktadır.

Hızlı geçişte [ y^+\approx60\text{–}140 ] aralığında [ \tau_{xy} ] profili belirgin bir plato oluşturmaktadır. Bu plato, hız ve vortisite profillerinden belirlenen Couette benzeri tabakayla aynı bölgede ortaya çıkmaktadır. Isıtılmayan karşı duvarda benzer bir plato görülmemesi, etkinin genel bir anlık türbülans dalgalanması değil, ısıl özellik değişimiyle ilişkili olduğunu desteklemektedir.

Normal gerilmeler: duvara dik sıkışma ve yanal yayılma

Şekil 3.7’de [ \tau_{yy} ] ve [ \tau_{zz} ] normal gerilme bileşenleri incelenmektedir. Bu bileşenlerin büyüklüğü ana [ \tau_{xy} ] kayma gerilmesinden yaklaşık iki mertebe küçüktür.

Geçiş ilerledikçe:

• [ \tau_{yy} ] daha negatif hâle gelerek duvara dik sıkışmayı,
• [ \tau_{zz} ] daha pozitif hâle gelerek kanal genişliği doğrultusunda yayılmayı göstermektedir.

Başlangıçta [ \tau_{yy} ] ekstremumu yaklaşık [ y^+=20 ], [ \tau_{zz} ] ekstremumu yaklaşık [ y^+=5 ] konumundadır. [ t^*=0 ] anında bu iki karakteristik konum yaklaşık [ y^+=10 ] ve [ y^+=9 ] düzeyinde birbirine yaklaşmaktadır.

Bu gerilme düzeni, düşük hızlı çizgilerde oluşan sıcak lekelerin duvara dik yönde incelirken kanal genişliği yönünde uzamasını desteklemektedir. Araştırmacılar bunun ısıyı yakın duvar bölgesinde tutarak çekirdeğe doğru etkin uzaklaştırmayı sınırlandırdığını belirtmektedir.

Duvar-normal sıcaklık profilleri

Akışkan sıcaklığı:

\[ T^*=\frac{T-T_{\mathrm{sat}}}{\Delta T_{\mathrm{sub}}} \]

ile boyutsuzlaştırılmıştır. Kanal boyunca giriş, orta ve çıkış olmak üzere üç alt bölge karşılaştırılmıştır.

Bütün zamanlarda sıcaklık artışı büyük ölçüde [ y^+\lesssim100 ] içinde kalmakta, kanal çekirdeği giriş sıcaklığına yakın kalmaktadır. Bu nedenle kısa geçiş süresi boyunca toplam sıvı kütlesi değil, esas olarak ince yakın duvar tabaka ısınmaktadır.

Hızlı geçişte giriş, orta ve çıkış profilleri birbirine daha yakındır. Yarı kararlı geçişte aşağı akış yönündeki sıcaklık gelişimi daha belirgin ve termal giriş bölgesi daha uzundur.

\[ t^*=0 ] anında hızlı geçiş tüm ısıtılmış yüzey boyunca doyma sıcaklığının üzerinde bir yakın duvar tabaka oluştururken, yarı kararlı geçişte giriş bölgesinin bir kısmı hâlâ doyma sıcaklığının altında kalmaktadır.

Stanton sayısı ve ısı aktarım şiddeti

Yerel Stanton sayısı:

[ St(x^*,t^*) = \frac{q''_w(t^*)} {\rho c_pU\left[T_w(x^*,t^*)-T_{\mathrm{in}}\right]} \]

olarak tanımlanmıştır.

Şekil 3.9’da Stanton sayısı girişte en yüksek değerine ulaşmakta, aşağı akışta azalarak bir platoya yaklaşmaktadır. Plato, termal sınır tabakanın akış yönünde tam gelişmiş olduğu bölgeyi tanımlamaktadır.

Yarı kararlı geçiş

\[ \tau=100\ \mathrm{ms} ] durumunda Stanton sayısı geçiş boyunca yalnızca yaklaşık yüzde 8 değişmektedir. Bu, yerel ısı aktarım katsayısının duvar yükündeki artışa büyük ölçüde uyum sağladığını ve akışın yarı kararlı bir taşınım rejimine yakın kaldığını göstermektedir.

Hızlı geçiş

[ \tau=5\ \mathrm{ms} ] durumunda Stanton sayısı [ t^*=-2 ] ile [ t^*=-1 ] arasında bütün ısıtılmış uzunluk boyunca yaklaşık yüzde 30 artmakta, ardından kaynama başlangıcına kadar bu yüksek düzeye yakın kalmaktadır.

Bu artış, hızlı geçişte ısı aktarımının mutlak olarak zayıf olduğu anlamına gelmez. Tam tersine toplam aktarım güçlenmektedir; ancak ısı artış hızı türbülanslı hız alanının bütünüyle yeniden ayarlanabileceğinden daha yüksek olduğundan geçici iletim ve yerel enerji depolama daha önemli hâle gelmektedir.

Termal giriş uzunluğu

Akış yönü kanal yarı yüksekliğiyle:

[ x^*=\frac{x}{e} \]

biçiminde boyutsuzlaştırılmıştır.

Hızlı geçişte Stanton platosuna yaklaşık:

\[ x^*\approx2\pi \]

konumunda ulaşılmaktadır. Yarı kararlı geçişte termal olarak tam gelişmiş durum ancak yaklaşık:

\[ x^*\approx5\pi \]

civarında, test kesiti çıkışına yakın oluşmaktadır.

Hızlı ısıtmada daha kısa termal giriş uzunluğu, yüzeyin ani yüklemeye yerel ısı depolama ve geçici iletimle tepki vermesiyle ilişkilendirilmektedir. Yavaş ısıtmada türbülanslı taşınımın akış yönünde gelişmesi için daha fazla zaman bulunmakta ve termal yapı daha uzağa yayılmaktadır.

İki parametreli sigmoid sıcaklık yasası

Termal olarak tam gelişmiş bölgede sıcaklık profilleri şu değişkenle yeniden ölçeklenmiştir:

\[ \Theta=\frac{T-T_w}{T_{\mathrm{in}}-T_w} \]

\[ \Theta=0 ] duvar sıcaklığına, [ \Theta=1 ] giriş sıcaklığına karşılık gelmektedir.

Profillerin logaritmik duvar koordinatında iki parametreli bir sigmoid bağıntıya uyduğu gösterilmiştir:

[ \Theta(y^+) = \frac{(y^+)^k} {(y_0^+)^k+(y^+)^k} \]

Eşdeğer lojistik gösterim:

\[ \Theta(y^+) = \frac{1} {1+\exp\left[-k\left(\ln y^+-\ln y_0^+\right)\right]} \]

şeklindedir.

\[ y_0^+ ], [ \Theta=1/2 ] olan geçiş orta noktasını; [ k ] ise logaritmik koordinattaki geçiş keskinliğini göstermektedir. Büyük [ k ] daha keskin ve duvara daha sıkışmış, küçük [ k ] daha yaygın bir termal tabaka anlamına gelir.

TS1 uyum parametreleri

TS2 uyum parametreleri

Bütün uyumlarda [ R^2>0{,}99 ] elde edilmiştir. Hızlı geçişler genel olarak daha küçük [ y_0^+ ] ve daha büyük [ k ] üretmektedir. Bu sonuç, termal geçişin hızlı ısıtmada duvara daha yakın ve daha keskin olduğunu nicel olarak doğrulamaktadır.

[ y_0^+ ] değerinin önce dışa, ardından kaynama başlangıcına yakın yeniden duvara doğru hareket etmesi, normal gerilme analizindeki duvara dik sıkışmayla tutarlıdır.

Yerel Péclet göstergesi

Yerel taşınım ve iletim büyüklükleri:

[ Pe= \frac{\left|\rho c_p\mathbf{u}\cdot\nabla T\right|} {\left|\nabla\cdot(k\nabla T)\right|} \]

ile karşılaştırılmıştır. Çok büyük [ Pe ] değerlerini sınırlı bir aralığa dönüştürmek için:

\[ \xi=\frac{Pe}{1+Pe} \]

tanımlanmıştır.

• [ \xi=0 ] yalnızca iletim,
• [ \xi=1/2 ] eşit iletim ve taşınım,
• [ \xi=1 ] yalnızca taşınım baskınlığı anlamına gelmektedir.

Şekil 3.11’de [ t^*=-1 ] anındaki aşırı ısıtılmış bölgeler çoğunlukla akış yönündeki, iletim ağırlıklı şeritler hâlindedir. [ t^*=0 ] anında termal giriş bölgesinin aşağısında sıvının büyük bölümü doyma sıcaklığının üzerine çıkmakta; iletim baskın, karma ve taşınım baskın küçük bölgeler aynı film içinde birlikte bulunmaktadır.

Bu aşırı ısıtılmış film metastabildir. Çalışma, uygun bir bozucu etki veya çekirdeklenme bölgesi oluştuğunda geniş bir yüzey alanında kabarcık oluşumuna yatkın olduğunu ileri sürmektedir. Kabarcık çekirdeklenmesi hesap modelinde yer almamaktadır.

Isı denkleminin duvar-normal bölgeleri

Depolama, iletim ve taşınım terimleri anlık duvar ısı akısının kanal yarı yüksekliğine oranıyla ölçeklendirilmiştir:

\[ \frac{q''_w}{e} \]

Şekil 3.12’ye göre termal sınır tabaka yaklaşık üç bölgeye ayrılmaktadır:

Hızlı geçişte operatörlerin ekstremumları daha küçük [ y^+ ] konumlarında bulunmaktadır. Kaynama başlangıcına yaklaşıldıkça enerji aktarım etkinliği giderek duvara daha yakın bir bölgede yoğunlaşmaktadır.

Yönsel ısı akısı ayrışımı

Isı denklemi akı biçiminde:

\[ \rho c_p\frac{\partial T}{\partial t} = \nabla\cdot\left(\mathbf{j}_q+\mathbf{j}_{\mathrm{adv}}\right) \]

olarak yazılmıştır. İletim ve taşınım akıları:

\[ \mathbf{j}_q=-k\nabla T \]

\[ \mathbf{j}_{\mathrm{adv}} =-\rho c_p\mathbf{u}(T-T_{\mathrm{in}}) \]

biçimindedir.

Şekil 3.13’te mutlak akı büyüklükleri karşılaştırıldığında dört mertebeye yayılan şu sıra elde edilmiştir:

\[ |j_{\mathrm{adv}}^x| > |j_{\mathrm{adv}}^z| > |j_{\mathrm{adv}}^y| \approx |j_q^y| > |j_q^z| > |j_q^x| \]

Akış yönündeki taşınım

En büyük bileşendir; çünkü ortalama akış hızı diğer hız bileşenlerinden çok daha büyüktür. Isıyı ısıtıcıdan aşağı akışa taşır fakat duvardan kanal çekirdeğine doğrudan uzaklaştırmaz.

Kanal genişliği yönündeki taşınım

İkinci büyük bileşendir. Yarı akış yönlü vortekslerin sıcaklığı düşük hızlı çizgilerden yanlara taşımasını ve sıcaklık alanının kanal genişliği boyunca homojenleşmesini destekler.

Duvara dik taşınım ve iletim

Bu iki bileşen benzer büyüklüktedir. Yarı kararlı geçişte duvara dik taşınım, uygulanan ısı akısına göre daha güçlüdür. Akış alanı yavaş yüklemeye uyum sağlamak için yeterli zaman bulduğundan ısı sıcak duvardan daha etkin biçimde uzaklaştırılmaktadır.

Hızlı geçişte duvara dik taşınımın göreli gelişimi daha sınırlı kalmakta, iletim ve yerel enerji depolama daha büyük rol oynamaktadır. Bu durum daha yüksek duvar aşırı sıcaklığına rağmen sıcak tabakanın duvara yakın kalmasına katkıda bulunmaktadır.

İletim bileşenleri

Kanal genişliği ve akış yönündeki iletim, en küçük akı bileşenleridir. Kanal genişliği iletimi yaklaşık [ y^+\approx7 ], akış yönündeki iletim ise yaklaşık [ y^+\approx9 ] konumunda en yüksek değerine ulaşmaktadır.

Hızlı ve yarı kararlı ısıtmanın karşılaştırılması

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

• Sayısal model doğrudan deneysel duvar ısı akısı geçmişiyle çalıştırılmıştır.
• Duvar sıcaklıkları deneysel ölçümlerle karşılaştırılmıştır.
• Sıcaklığa bağlı viskozite ve ısıl iletkenlik hesaba katılarak iki yönlü termal–hidrodinamik bağ kurulmuştur.
• Akışın yalnızca ortalama değerleri değil, anlık üç boyutlu sıcaklık ve girdap yapıları incelenmiştir.
• Çok yüksek hücre sayıları ve ince yakın duvar ağı kullanılmıştır.
• Alt ağ modeli katkıları nicel olarak değerlendirilmiştir.
• Hız, vortisite, kayma gerilmesi, normal gerilmeler, Stanton sayısı ve yönsel ısı akıları birbirini tamamlayan biçimde analiz edilmiştir.
• Sıcaklık profilleri için iki parametreli ve yüksek uyumlu bir sigmoid bağıntı önerilmiştir.
• Hızlı ve yarı kararlı geçişlerin kaynama başlangıcına farklı tek fazlı yollardan ilerlediği gösterilmiştir.

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

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Simülasyonlar tek fazlıdır ve ilk kabarcık oluştuktan sonraki dönemi kapsamaz.
• Kabarcık çekirdeklenmesi, kabarcık büyümesi, yüzey boşlukları ve iki fazlı arayüz dinamiği modellenmemiştir.
• Yoğunluk ve özgül ısı sabit tutulmuştur; yalnızca viskozite ve ısıl iletkenlik sıcaklığa bağlıdır.
• Akış sıkıştırılamaz kabul edilmiştir.
• Yerçekimi, düşük toplu Richardson sayısına dayanılarak ihmal edilmiştir.
• LES modeli hızlı denge dışı geçişlerde türbülans üretimi ve yutulması arasındaki yerel dengeyi tam temsil etmeyebilir.
• Deneysel ısı akısı sınır koşulu uzaysal ortalamadır; ısıtıcı kenarlarındaki yerel kayıplar modele dâhil edilmemiştir.
• Sonuçlar iki kanal geometrisi, belirli Reynolds sayıları, 10 bar basınç ve 100 K giriş soğutmasıyla sınırlıdır.
• Çalışma ağırlıklı olarak 5 ve 100 ms uç durumlarını karşılaştırmakta; 20 ms ara durumunun ayrıntılı mekanizma analizi daha sınırlıdır.
• Zamansal ortalama ve Reynolds ayrışımı uygulanmamış, anlık alanların uzaysal ortalamalarına odaklanılmıştır.
• Ham üç boyutlu LES alanları veri büyüklüğü nedeniyle açık bir arşive yüklenmemiştir.
• Sigmoid bağıntının farklı akışkanlar, basınçlar ve daha geniş Reynolds sayısı aralığında genellenebilirliği gösterilmemiştir.

Çalışma neyi göstermektedir?

• Hızlı duvar ısıtması, kaynama başlamadan önce yakın duvar hız alanını değiştirmektedir.
• Sıcaklığa bağlı viskozite azalması, hız artışı ve daha güçlü açıklık vortisitesiyle ilişkilidir.
• Tampon tabakada Couette benzeri sabit kesmeli bir bölge oluşabilmektedir.
• Isıtılmış yapılar duvara dik yönde sıkışmakta ve kanal genişliği yönünde yayılmaktadır.
• Hızlı geçişte Stanton sayısı yaklaşık yüzde 30 yükselmektedir.
• Hızlı geçişin termal giriş uzunluğu yarı kararlı geçişten daha kısadır.
• Termal olarak tam gelişmiş bölgedeki sıcaklık profilleri iki parametreli sigmoid yasayla çok yüksek doğrulukta temsil edilmektedir.
• Hızlı ve yarı kararlı geçişler farklı uzamsal aşırı ısıtma desenleri oluşturmaktadır.

Çalışma neyi göstermemektedir?

• Kabarcığın tam olarak hangi yüzey kusurunda ve hangi mikroskobik mekanizmayla çekirdeklendiğini göstermemektedir.
• Çekirdekli kaynama başladıktan sonra ısı aktarım katsayısının nasıl değiştiğini hesaplamamaktadır.
• Kritik ısı akısını veya film kaynamasına geçişi öngörmemektedir.
• Sigmoid sıcaklık yasasının evrensel bir türbülans yasası olduğunu kanıtlamamaktadır.
• Bütün nükleer reaktör geometrilerinin aynı yüzde 30 Stanton artışını göstereceğini ortaya koymamaktadır.
• Değişken viskozitenin tek başına bütün türbülans değişikliklerinin nedeni olduğunu bağımsız kontrol simülasyonlarıyla kesin olarak ayırmamaktadır.
• Hızlı geçişte kaynamanın kesin olarak bütün yüzeyde aynı anda başlayacağını doğrudan simüle etmemektedir; yalnızca tek fazlı sıcaklık alanı bu yönde bir başlangıç durumu göstermektedir.

Bilimsel ve teknolojik önemi

Çalışma, hızlı ısı yüklerinin güvenli biçimde yönetilmesinde yalnızca ortalama duvar sıcaklığının yeterli olmadığını göstermektedir. Aynı duvar sıcaklığına yaklaşan iki geçiş, yakın duvar türbülansı ve aşırı ısıtılmış sıvının uzamsal yapısı bakımından farklı olabilir.

Nükleer reaktör güvenliği açısından bulgular, reaktivite artışı sırasında kaynama başlangıcını öngören modellerde ısıtma zaman ölçeğinin, değişken akışkan özelliklerinin ve termal giriş uzunluğunun ayrıca dikkate alınması gerektiğini göstermektedir.

Elektronik soğutma ve roket motorlarında hızlı güç değişimleri, duvarın yakınında benzer biçimde termal depolama ve akışın gecikmeli yeniden düzenlenmesine yol açabilir. Bununla birlikte çalışma bu sistemlere ait özgül geometrileri doğrudan modellememiştir; bağlantı yöntemsel düzeydedir.

Gelecekteki çalışmaların tek fazlı LES sonuçlarını kabarcık çekirdeklenme modelleri, yüzey mikrogeometrisi ve arayüz izleme yöntemleriyle birleştirmesi; böylece hesaplanan aşırı ısıtılmış filmin gerçek kaynama başlangıcına nasıl dönüştüğünü incelemesi gerekmektedir.

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

Teknik yöntem özeti

Başlıca nicel bulgular

Şekil ve tablo bulgularının özeti

Veri işleme ve istatistiksel değerlendirme

Çalışma zamansal ortalama kullanmamaktadır. Anlık alanlar kanal genişliği doğrultusunda ve termal olarak tam gelişmiş aşağı akış bölgelerinde uzaysal olarak ortalanmıştır.

Deneysel doğrulamada göreli duvar sıcaklığı farkı kullanılmıştır. Sigmoid sıcaklık profilleri için [ R^2 ] belirleme katsayısı hesaplanmıştır. Hipotez testi, p değeri, güven aralığı veya olasılıksal belirsizlik yayılımı uygulanmamıştır.

Çalışma, en büyük türbülanslı hareketlerin yavaş geçişlerde test kesitinden geçmesi nedeniyle doğrulama farkında salınımlar oluşabildiğini belirtmektedir. Bu salınımlar ayrı bir spektral veya istatistiksel modelle nicel olarak ayrıştırılmamıştır.

Kaynak ve Yöntem Notu

Çalışmanın tam özgün adı:Exponential Wall Heat-Flux Transients in Turbulent Channel Flow

Yazarlar ve sıraları: Elie Roumet; Raksmy Nop; Nicolas Dorville; Christophe Bourcier; Elie Saikali; Marie-Christine Duluc.

Yazar–kurum eşleşmeleri:

• Elie Roumet: Université Paris-Saclay, CEA, Service de Thermo-hydraulique et de Mécanique des Fluides; ayrıca Conservatoire National des Arts et Métiers.
• Raksmy Nop: Université Paris-Saclay, CEA, Service de Thermo-hydraulique et de Mécanique des Fluides.
• Nicolas Dorville: Université Paris-Saclay, CEA, Service de Thermo-hydraulique et de Mécanique des Fluides.
• Christophe Bourcier: Université Paris-Saclay, CEA, Service de Génie Logiciel pour la Simulation.
• Elie Saikali: Université Paris-Saclay, CEA, Service de Génie Logiciel pour la Simulation.
• Marie-Christine Duluc: Conservatoire National des Arts et Métiers.

Kurum adresleri:

• Université Paris-Saclay, CEA, 91191 Gif-sur-Yvette, Fransa.
• Conservatoire National des Arts et Métiers, 75003 Paris, Fransa.

Sorumlu yazar: Elie Roumet.

İletişim adresi: elie.roumet@cea.fr

Eş birinci yazar: Eş katkı veya eş birinci yazarlık beyanı bulunmamaktadır.

DOI:10.2139/ssrn.6947712

Resmî kayıt:SSRN çalışma sayfası

Yayın platformu: SSRN.

SSRN yayımlanma tarihi: 15 Haziran 2026.

Yayın yılı: 2026.

Sayfa sayısı: 34.

Kaynak türü: Termal-hidrolik, hesaplamalı akışkanlar dinamiği ve türbülanslı ısı aktarımı alanında araştırma preprinti.

Dergi: Hakemli bir dergi adı, cilt veya sayı bilgisi bulunmamaktadır.

Özgün yayınevi: Ayrı bir hakemli akademik yayınevi doğrulanmamıştır.

Hakemlik durumu: Çalışma hakem değerlendirmesinden geçmemiştir. PDF’nin her sayfasında “This preprint research paper has not been peer reviewed” uyarısı bulunmaktadır.

Yazar katkıları: Elie Roumet kavramsallaştırma, yöntem, doğrulama, biçimsel analiz, araştırma, görselleştirme ve ilk taslak yazımını; Raksmy Nop kavramsallaştırma, kaynaklar, finansman temini, proje yönetimi ve gözden geçirmeyi; Nicolas Dorville kavramsallaştırma, proje yönetimi ve gözden geçirmeyi; Christophe Bourcier ile Elie Saikali yazılım, veri düzenleme ve gözden geçirmeyi; Marie-Christine Duluc kavramsallaştırma, danışmanlık ve gözden geçirmeyi üstlenmiştir.

Finansman ve hesaplama kaynağı: Çalışma, GENCI tarafından sağlanan A0152A13031 ve A0172A13031 tahsisleri kapsamında TGCC yüksek başarımlı hesaplama kaynaklarını kullanmıştır.

Çıkar çatışması: Yazarlar, çalışmayı etkileyebilecek bilinen mali çıkar veya kişisel ilişki bulunmadığını beyan etmiştir.

Veri erişimi: Ham üç boyutlu LES alanları yüksek veri hacmi nedeniyle açık bir depoya yüklenmemiştir. Şekilleri ve temel nicel sonuçları destekleyen işlenmiş verilerin sorumlu yazardan makul talep üzerine alınabileceği belirtilmiştir.

Yapay zekâ kullanım beyanı: Yazarlar, makalenin okunabilirliğini geliştirmek amacıyla OpenAI ChatGPT kullandıklarını; metni daha sonra gözden geçirip düzenlediklerini ve içerikten tam sorumluluk aldıklarını bildirmiştir.

Bu Türkçe makale, yüklenen çalışmanın tamamı; formülleri, sınır koşulları, sayısal yöntemleri, grafik ve alan görselleri, tablo değerleri, ekleri ve yazar beyanları incelenerek hazırlanmıştır. Bilimsel açıklamalar yalnızca çalışmada sunulan yöntem, veriler ve yorumlarla sınırlandırılmış; dış kaynaklardan ek bir bilimsel bulgu eklenmemiştir.

Çalışmanın ana yöntemsel sınırı, yalnızca çekirdekli kaynama başlangıcından önceki tek fazlı dönemi modellemesidir. Bu nedenle sonuçlar kabarcık çekirdeklenmesi, iki fazlı ısı aktarımı, kritik ısı akısı veya film kaynaması için doğrudan öngörü olarak kullanılmamalıdır. Ayrıca çalışma henüz bağımsız hakem değerlendirmesinden geçmemiştir.