y+ nedir?
Boyutsuz duvar mesafesi (y+ veya y-artı olarak yazılır), hesaplamalı akışkanlar dinamiğinde (CFD) önemli bir parametredir ve türbülanslı bir sınır tabakasındaki ilk ağ hücresinin yüksekliğini viskoz uzunluk ölçeğine göre tanımlar. Duvar yakınındaki ağınızın viskoz alt tabakayı mı, geçiş (buffer) tabakasını mı çözdüğünü, yoksa logaritmik yasa bölgesinde duvar fonksiyonlarına mı dayandığını gösterir.
Bu hesaplama aracı nasıl kullanılır?
Akışkan yoğunluğu \(\rho\), dinamik viskozite \(\mu\), duvar kayma gerilmesi \(\tau_w\) ve duvardan ilk hücre merkezine olan \(y\) mesafesini girin. Araç önce sürtünme hızı \(u_\tau\) değerini, ardından y+ değerini hesaplar. Duvarda çözüm yapan (düşük Reynolds) modeller için \(y^+ \approx 1\), duvar fonksiyonları kullanırken ise yaklaşık 30–300 aralığını hedefleyin.
Formülün açıklaması
Önce sürtünme hızı, duvar kayma gerilmesinden bulunur:
$$u_\tau = \sqrt{\frac{\tau_w}{\rho}}$$Bu hız ölçeği, duvar yakınındaki türbülansı karakterize eder. Ardından
$$y^+ = \frac{\rho \, u_\tau \, y}{\mu}$$ifadesi, fiziksel \(y\) mesafesini viskoz uzunluk \(\nu/u_\tau\) ile ölçeklendirir; burada \(\nu = \mu/\rho\) kinematik viskozitedir. Eşdeğer olarak \(y^+ = u_\tau \, y / \nu\) şeklinde de yazılabilir.
Örnek hesaplama
Hava için \(\rho = 1{,}225 \ \text{kg/m}^3\), \(\mu = 1{,}81\times10^{-5} \ \text{Pa}\cdot\text{s}\), \(\tau_w = 0{,}1 \ \text{Pa}\) ve \(y = 5\times10^{-5} \ \text{m}\) alalım. Bu durumda
$$u_\tau = \sqrt{\frac{0{,}1}{1{,}225}} = 0{,}28571 \ \text{m/s}$$$$y^+ = \frac{1{,}225 \times 0{,}28571 \times 5\times10^{-5}}{1{,}81\times10^{-5}} \approx 0{,}967$$olur. \(y^+ \approx 1\) hedefine ulaşmak için buradaki ilk hücre zaten oldukça iyi konumlanmıştır.
Sıkça Sorulan Sorular
Hangi y+ değerini hedeflemeliyim? Duvarda çözüm yapan modeller için (örneğin k-ω SST düşük Reynolds) yaklaşık 1, standart duvar fonksiyonları kullanırken ise 30–300 arası.
y+ neden boyutsuzdur? Duvar mesafesinin \(\nu/u_\tau\) viskoz uzunluk ölçeğine bölünmesiyle elde edilir, böylece birimler birbirini götürür.
\(\tau_w\) değerini bilmiyorsam ne yapmalıyım? Cilt sürtünme katsayısından tahmin edebilirsiniz: \(\tau_w = \tfrac{1}{2} \rho U^2 C_f\); burada \(C_f\), \(0{,}058 \cdot Re^{-0{,}2}\) gibi korelasyonlardan gelir.
