Endüstriyel Tesislerde CFD ile Soğutma Doğrulaması: CNC Atölyesi Örneği
Bir CNC atölyesinde soğutma ve havalandırmayı boyutlandırırken genellikle klasik ısı yükü hesabından çıkan bir toplam kapasite ve birkaç diki üfleme menfezi ile işe başlarız. Sorun şu ki, tezgâhların ısıl kararlılığı (dolayısıyla işleme hassasiyeti) ve operatör konforu, odanın ortalama sıcaklığından çok, iş hacmindeki yerel sıcaklık dağılımına ve hava akış paternine bağlıdır. Spindle motorları, sürücü panoları, hidrolik üniteleri ve kesme prosesi noktasal ısı kaynakları oluşturur; bu kaynaklar üzerinde oluşan yükselen sıcak hava kolonları (plume) ve menfez jetlerinin bunları süpürüp süpürmediği, tek başına yük hesabıyla görülmez. Bu yazıda, böyle bir atölye için yapılan HVAC tasarımının CFD ile nasıl doğrulandığını; konjuge ısı transferi (CHT), k-omega SST türbülans modeli, ısı kaynaklarının kaynak terimi (source term) olarak modellenmesi ve mesh bağımsızlığı açılarından ele alıyorum. Proje verileri anonimdir; amaç metot ve karar mantığını aktarmak.
Problemin tanımı ve modelleme stratejisi
Doğrulamak istediğimiz temel soru basitti: “Tasarlanan üfleme debisi ve menfez yerleşimi, iş bölgesinde hedef sıcaklık bandını sağlıyor ve tezgâh çevresinde durgun sıcak cepler bırakmıyor mu?” Bunu deneysel olarak ölçmek mümkün ama tesis daha kurulmadan karar vermek gerekiyorsa CFD pratik tek yol.
İlk stratejik karar, ısı kaynaklarını nasıl temsil edeceğimizdir. İki uç var:
- Tam CHT (konjuge ısı transferi): Her tezgâhın gövdesi, panosu, motoru katı bölge olarak modellenir; katıda iletim, havada taşınım eşzamanlı ve arayüzde birebir kuple çözülür. Fiziksel olarak en doğrusu, ama her makinenin gerçek geometrisini, malzeme iletkenliğini ve iç ısı üretimini bilmeyi gerektirir. Atölye ölçeğinde bu, hem veri hem hücre sayısı açısından pahalıdır.
- Hacimsel kaynak terimi (lumped source): Isı kaynağı, ilgili katıyı temsil eden bir hacme dağıtılmış enerji kaynağı (W/m³) olarak girilir; geometri kabaca bir kutu ile temsil edilir. Isının havaya geçiş dinamiği (plume, buoyancy) doğru yakalanır, ama katı içi sıcaklık dağılımı hedef değilse gereksiz detay elenmiş olur.
Bizim sorumuz oda hava akışıyla ilgili olduğundan, ana kaynakları enerji kaynak terimi olarak modelledik; sadece ısıl kütlesi büyük ve yüzey sıcaklığı kritik olan birkaç bileşen için lokal CHT bölgesi tanımladık. Menfez jetlerini gerçek geometriyle çözmek yerine, üfleme ağızlarını momentum + enerji kaynağı ile temsil ettik: verilen debi ve sıcaklıktaki havayı, jet yönünde momentum katarak besledik. Bu, difüzör iç geometrisini çözmeden doğru atım (throw) ve karışımı yakalamanın makul bir yoludur.
Kaldırma kuvveti (buoyancy) bu problemde belirleyicidir; sıcak plume’ler doğal taşınımla yükselir. Bunu ihmal edemeyiz, dolayısıyla yoğunluğu sıcaklığa bağlı (ideal gaz ya da sıcaklık farkı küçükse Boussinesq yaklaşımı) alıp yerçekimini aktif tuttuk. Grashof/Rayleigh mertebesi, zorlanmış ve doğal taşınımın birlikte etkin olduğu karışık taşınım (mixed convection) rejiminde olduğumuzu gösterir; bu yüzden ne sadece jetleri ne de sadece plume’leri baskın kabul edip basitleştiremeyiz.
Türbülans modeli ve duvar çözümü (y+)
Türbülans modeli olarak k-omega SST seçtik. Gerekçe standart: SST, duvar yakınında k-omega’nın iyi davranışını, serbest akışta ise k-omega’nın serbest-akım hassasiyetini k-epsilon benzeri davranışla harmanlar. Menfez jetlerinde ayrılma, tezgâh arkalarında ters basınç gradyanı ve resirkülasyon bölgeleri beklediğimiz için, bu tür ayrılmalı iç akışlarda SST’nin standart k-epsilon’a göre daha güvenilir olması bizim için belirleyiciydi.
Duvar çözümünde kritik parametre y+. SST’nin otomatik duvar muamelesi (automatic/blended wall treatment) hem düşük-Re çözümüne (y+ ≈ 1, viskoz alt tabaka çözülür) hem de duvar fonksiyonlarına (y+ ≈ 30–300) izin verir, ancak kaçınılması gereken bant y+ ≈ 5–30 tampon tabakasıdır — orada ne biri ne öteki doğru davranır. Atölye ölçeğinde her duvarda y+ ≈ 1 tutmak hücre bütçesini patlatır; bu yüzden duvarlarda duvar fonksiyonu bandını hedefledik ve ısı transferinin kritik olduğu yüzeylerde (kaynak yüzeyleri) meshi sıklaştırdık. Çözüm sonrası y+ dağılımını mutlaka kontrol ettik — bu, HVAC CFD’sinde en sık atlanan doğrulama adımıdır.
Sınır koşulları özetle şöyle kuruldu:
| Sınır | Tip | Tanım |
|---|---|---|
| Üfleme menfezleri | Hız girişi / momentum kaynağı | Tasarım debisi, üfleme sıcaklığı, türbülans yoğunluğu %5, hidrolik çap üzerinden uzunluk ölçeği |
| Emiş / egzoz | Basınç çıkışı | 0 Pa gauge, geri akışta oda sıcaklığı |
| Dış duvarlar / çatı | Duvar | Kabuk iletimi veya sabit ısı akısı (zarf yüküne göre) |
| Zemin | Duvar | Adyabatik kabul (konservatif) |
| Isı kaynakları (tezgâh, pano, hidrolik) | Enerji kaynak terimi | Bileşen başına ısı gücü, temsil hacmine W/m³ olarak dağıtılmış |
| Kritik bileşenler | CHT arayüzü | Katı-akışkan kuple, malzeme iletkenliği tanımlı |
Solver tarafında tipik kurulum parametreleri (illüstratif, projeye göre değişir):
# Çözüm kurulumu (özet)
solver : steady-state, incompressible (buoyant)
turbulence_model : k-omega SST
wall_treatment : automatic (low-Re + wall function blend)
gravity : (0, 0, -9.81) m/s^2 # buoyancy aktif
density : ideal_gas # veya Boussinesq (dT küçükse)
# Isı kaynağı (örnek bir tezgâh temsil hacmi)
heat_source:
type : energy (volumetric)
power : <bileşen ısı gücü> W
volume : <temsil kutusu hacmi> m^3 # -> W/m^3
# Ayrıklaştırma / yakınsama
schemes : 2nd order upwind (momentum, energy, turb)
residual_targets : U,p 1e-4 ; energy 1e-6
monitors : iş bölgesi ort. T, maks. yüzey T, çıkış T
Yakınsamayı yalnızca artıklara (residual) bakarak değil, fiziksel izleme büyüklüklerinin (monitor) sabitlenmesine bakarak değerlendirdik. Kaldırma kuvvetli akışlarda artıklar takılırken çözüm hâlâ salınabilir; iş bölgesi ortalama sıcaklığı düz bir çizgiye oturmadan yakınsadı demek yanıltıcı olur.
Mesh bağımsızlığı
Sonuçların meshten bağımsız olduğunu göstermeden hiçbir sayı savunulabilir değildir. Kaba, orta ve ince olmak üzere en az üç mesh seviyesi ürettik; her seviyede aynı izleme büyüklüğünü (iş bölgesi hacimsel ortalama sıcaklık ve maksimum yerel sıcaklık) karşılaştırdık. Menfez jetlerinin kesme bölgelerinde ve kaynak yüzeylerinde yerel sıklaştırma (refinement) uyguladık, çünkü gradyanların dik olduğu yer orasıdır.
Aşağıdaki tablo yöntemi göstermek için şematiktir (değerler illüstratif):
| Mesh | Hücre sayısı (mertebe) | İzlenen büyüklükteki değişim* |
|---|---|---|
| Kaba | ~1–2 M | referans |
| Orta | ~4–6 M | orta → kaba arası fark belirgin |
| İnce | ~10–15 M | ince → orta arası fark ihmal edilebilir |
*Kabul kriteri: ardışık iki seviye arasında izleme büyüklüğündeki değişim önceden belirlenen bir toleransın (ör. birkaç yüzde) altına inince mesh bağımsız kabul edilir. Mümkünse GCI (Grid Convergence Index) ile Richardson ekstrapolasyonu yaparak sayısal belirsizlik bandını raporlamak daha savunulabilir bir yaklaşımdır.
Pratikte “orta” mesh ile “ince” mesh arasındaki fark toleransın altına düştüğünde, mühendislik kararlarını orta mesh üzerinden veririz; ince mesh yalnızca doğrulama içindir.
Sonuçların yorumlanması
CFD çıktısını iki katmanda okuduk. Birincisi sıcaklık alanı: iş bölgesi kesitlerinde hedef sıcaklık bandının sağlanıp sağlanmadığı ve durgun sıcak ceplerin (özellikle tezgâh aralarında ve yüksek kaynakların üstünde) nerede oluştuğu. İkincisi akış paterni: menfez jetlerinin iş bölgesine ulaşıp ulaşmadığı, kısa devre yapıp yapmadığı (üflenen havanın işi görmeden doğrudan emişe kaçması) ve plume’lerin nasıl yükseldiği. Akım çizgileri ve hız vektörleri, “toplam debi yeterli ama dağılım kötü” durumunu ortaya çıkarır — ki klasik yük hesabının göremediği tam da budur. Örneğin bir kaynağın üstünde emiş menfezi yoksa, sıcak plume tavanda birikip resirkülasyonla iş bölgesine geri dönebilir; çözüm çoğu zaman kapasiteyi artırmak değil, menfez yönünü/yerini değiştirmektir.
Özet
- Atölye ölçeğinde ısı kaynaklarını çoğunlukla hacimsel enerji kaynak terimi, menfezleri ise momentum kaynağı olarak modellemek; yalnızca kritik bileşenlerde CHT kullanmak, doğruluk ile hesap maliyeti arasında sağlıklı dengeyi kurar.
- k-omega SST, ayrılmalı iç akışlar için uygundur; ama seçim, y+ dağılımının kontrolüyle tamamlanmalıdır — tampon tabaka bandından (y+ ≈ 5–30) kaçının.
- Kaldırma kuvveti ihmal edilemez; problem karışık taşınım rejimindedir, bu yüzden yerçekimi ve sıcaklığa bağlı yoğunluk aktif tutulmalıdır.
- Mesh bağımsızlığı ve fiziksel izleme büyüklüklerinin yakınsaması gösterilmeden hiçbir sayı raporlanmamalıdır; mümkünse GCI ile belirsizlik bandı verilmelidir.
- CFD’nin asıl katkısı toplam kapasiteyi teyit etmek değil, yerel dağılımı ve akış paternini görüp menfez yerleşimini optimize etmektir — kısa devre ve durgun sıcak cepler ancak böyle yakalanır.