Çok Fazlı Akışların Modellenmesinde CFD'nin Rolü
Hesaplamalı Akışkanlar Dinamiği (CFD), motorlardaki yanma odalarından meteorolojik modellere, uçak kanatlarındaki buzlanmadan kan damarlarındaki ilaç taşınımına kadar sayısız alanda mühendislerin en güçlü araçlarından biridir. Akışkanın içine katı parçacıkların veya sıvı damlacıkların karıştığı sistemler "çok fazlı akışlar" (multiphase flows) olarak adlandırılır.
Bir CFD simülasyonu kurarken en kritik karar, yazılımın parçacıkları nasıl izleyip hesaplayacağını belirlemektir. Milyonlarca parçacığın tek tek izlenmesi mi (Lagrangian yaklaşımı), yoksa parçacıkların ayrı bir akışkan gibi mi kabul edilmesi (Eulerian yaklaşımı) gerektiğine karar verilirken matematiksel pusula Stokes sayısı ($Stk$)'dır. Modelinizi kurmadan önce sisteminizin parametrelerini Stokes Sayısı Hesaplama aracımızı kullanarak analiz edebilirsiniz.
Stokes Sayısının CFD Yöntem Seçimine Etkisi
Stokes sayısı, parçacığın gevşeme süresinin akışın zaman ölçeğine oranıdır. CFD simülasyonlarında karakteristik akış zamanı genellikle bir engel boyutu üzerinden veya türbülans girdaplarının ömrü üzerinden (Kolmogorov zaman ölçeği) tanımlanır.
Stokes sayısının büyüklüğüne göre uygun matematiksel modeller şunlardır:
- $Stk \ll 1$ (Çok Düşük Stokes Sayısı): Parçacıklar akışkanın her hareketine anında uyum sağlar. Bu durumda Eulerian-Eulerian modeller kullanılabilir. Parçacıkların kendi ataleti olmadığı için ayrı ayrı izlenmelerine gerek yoktur; akışkanla karışmış yoğun bir gaz (veya karışım) gibi modellenerek hesaplama maliyetinden (CPU süresinden) devasa bir tasarruf sağlanır.
- $Stk \approx 1$ (Orta Stokes Sayısı): Parçacıklar akıştan kısmen bağımsız hareket eder, yörüngeleri karmaşıktır. Euler-Lagrange yaklaşımı olan Kesikli Faz Modeli (Discrete Phase Model - DPM) zorunludur. Burada akışkan Eulerian bir ağ üzerinde (grid), her bir parçacık veya parçacık paketi ise Lagrangian yöntemiyle kendi yörüngesinde tek tek hesaplanarak izlenir.
- $Stk \gg 1$ (Yüksek Stokes Sayısı): Parçacıkların ataleti çok yüksektir, aerodinamik sürükleme kuvvetinden ziyade yerçekimi, kendi kinetik enerjileri ve çarpışmalar (parçacık-duvar, parçacık-parçacık) hareketi dikte eder. DPM kullanılır ancak çözücü ayarlarında sürükleme (drag) denklemlerinden çok temas denklemlerine odaklanılır.
CFD Türbülans Modülasyonu ve Parçacık Etkileşimi
Stokes sayısı sadece yörüngeyi değil, parçacıkların akışkanı nasıl etkilediğini de (coupling effect) belirler:
- Tek Yönlü Etkileşim (One-way coupling): Parçacıkların kütle oranı çok düşükse akışkan parçacıkları sürükler, ancak parçacıklar akışkanı etkilemez.
- İki Yönlü Etkileşim (Two-way coupling): Parçacık yoğunluğu arttıkça (ve özellikle $Stk \sim 1$ civarında iken), parçacıklar türbülansı etkilemeye başlar (Türbülans Modülasyonu). Yüksek $Stk$ değerine sahip büyük parçacıklar geçtikleri yerde yeni girdaplar oluşturarak (wake) türbülansı artırırken, düşük $Stk$ değerine sahip ince parçacıklar türbülans enerjisini emerek (viskoz sönümleme) akışı sakinleştirir. CFD analizlerinde bu ayrımı yapmak için parçacık Stokes sayısının türbülans girdap boyutuna göre hesaplanması hayati önem taşır.
Gerçekçi Örnek: Egzoz Borusundaki Kurum Simülasyonu
Senaryo: Bir dizel motorun egzoz borusundaki kurum (is) parçacıklarının akışını ANSYS Fluent veya OpenFOAM'da modelleyeceksiniz. Egzoz gazı hızı yüksek, kurum parçacıkları ise oldukça küçüktür. Hangi modeli kullanmanız gerektiğine (DPM mi, Yoksa Mixture/Eulerian model mi) karar vermek istiyorsunuz.
Veriler:
- Kurum parçacık çapı ($d_p$): $2 \ \mu m \ (2 \times 10^{-6} \ m)$
- Kurum yoğunluğu ($\rho_p$): $1800 \ kg/m^3$
- Egzoz gazı hızı ($U$): $25 \ m/s$
- Egzoz borusu çapı ($L$ - karakteristik uzunluk): $0.1 \ m$
- Egzoz gazı viskozitesi ($\mu$): $3 \times 10^{-5} \ Pa\cdot s$ (sıcak gaz olduğu için havadan yüksektir)
- Cunningham düzeltmesi ($C_c$): $\approx 1$
Hesaplama Adımları:
- Parçacık Gevşeme Süresi ($\tau_p$):
$\tau_p = \frac{\rho_p \cdot d_p^2 \cdot C_c}{18 \cdot \mu} = \frac{1800 \cdot (2 \times 10^{-6})^2 \cdot 1}{18 \cdot 3 \times 10^{-5}} = \frac{7.2 \times 10^{-9}}{5.4 \times 10^{-4}} \approx 1.33 \times 10^{-5} \ s$ - Akış Zaman Ölçeği ($\tau_f$):
$\tau_f = \frac{L}{U} = \frac{0.1}{25} = 0.004 \ s$ - Stokes Sayısı ($Stk$):
$Stk = \frac{\tau_p}{\tau_f} = \frac{1.33 \times 10^{-5}}{0.004} \approx 0.0033$
Yorum: Hesaplanan Stokes sayısı ($Stk = 0.0033$), $0.1$'den çok küçüktür. Bu parçacıklar, içinden geçtikleri gazın ataletinden veya yörüngesinden neredeyse hiç sapmazlar.
CFD Stratejisi: Bu değer, CFD analistine çok net bir mesaj verir: "Parçacıkların her birini Lagrangian DPM modeliyle tek tek izleyerek milyonlarca ekstra denklem çözmeye gerek yoktur." Eğriler boyunca akışkanla tam bir uyum içinde olduklarından, Eulerian skaler taşıma denklemi (veya Mixture modeli) kullanmak yeterlidir. Bu karar, bilgisayarın çözüm süresini günlerden saatlere indirebilir. Farklı partikül boyutlarında $Stk$'nın kritik seviye olan $1$'e ne zaman ulaşacağını görmek için Stokes Sayısı Hesaplama aracımızı kullanabilirsiniz.
İlgili Uyarılar ve Sınır Durumları
CFD analizlerinde Stokes sayısını yorumlarken şu kısıtlamalar dikkate alınmalıdır:
- Zaman Ölçeği Seçimindeki Belirsizlik: Yukarıdaki örnekte makroskopik sistem boyutunu ($0.1 \ m$) baz aldık. Ancak konu "parçacık-türbülans etkileşimi" olduğunda, makroskopik hız yerine türbülansın en küçük ölçeği olan Kolmogorov girdap hızları ve boyutları kullanılmalıdır. Bu durumda elde edilen sayıya "Parçacık Stokes Sayısı" veya "Türbülans Stokes Sayısı ($St_{\eta}$)" denir. Aynı sistemde makro $Stk \ll 1$ çıkarken, mikro (Kolmogorov) $Stk > 1$ çıkabilir.
- Ağ (Mesh) Boyutu Etkisi: CFD çözümlerinde, DPM (Lagrangian) paketlerinin boyutu hesaplama ağının (mesh) hücresinden büyük olmamalıdır. Düşük Stokes sayısında parçacıklar spesifik hücrelerde birikme eğilimi göstereceğinden sayısal kararsızlıklara (numerical instability) yol açabilirler.
- Faz Değişimi: Bulaşıcı veya buharlaşan damlacıklar (örneğin sprey yanmasında veya buzlanmada) modellendiğinde, parçacığın kütlesi ($\rho_p$) ve çapı ($d_p$) her iterasyonda değişir. Bu nedenle Stokes sayısı simülasyon süresince sabit kalmaz, dinamik olarak yerel şartlara göre güncellenmek zorundadır.
