Bir çözümün bittiğine karar vermek, çözümü başlatmaktan daha çok deneyim ister. En sık görülen kabul şudur: artıklar belirli bir eşiğin altına indi, öyleyse iş tamam. Bu tek başına hiçbir şey söylemez [1].
Artık (residual) gerçekte neyi ölçer
Artık, korunum denklemlerinin ayrıklaştırılmış hâlinin ne kadar sağlanmadığını gösterir. Ölçeklendirilmiş bir büyüklüktür ve ölçekleme başlangıç koşuluna göre yapılır. Kötü bir başlangıçtan başlarsanız artıklar hızla "düşer"; bu düşüş çözümün doğruluğu hakkında bilgi vermez.
Aynı problemi iki farklı başlangıç koşulundan başlatın: artık eğrileri bambaşka görünür, çözüm aynıdır. Ölçüt olarak kullanılan şeyin, ölçmek istediğiniz şeyle ilgisi bu kadar zayıftır.
Yakınsamanın üç ölçütü
1. Artıklar düşmüş ve yatay seyrediyor olmalı. Değerinden çok davranışı önemlidir: hâlâ düşüyorsa çözüm oturmamıştır, salınıyorsa fizik zamana bağlı olabilir.
2. İzlenen büyüklükler sabitlenmiş olmalı. Kararı verecek büyüklükleri (basma yüksekliği, sürükleme kuvveti, çıkış sıcaklığı, basınç kaybı) iterasyon boyunca kaydedin. Artıklar yatayken bu büyüklükler hâlâ sürükleniyorsa çözüm yakınsamamıştır.
3. Kütle ve enerji dengesi kapanmalı. Girenle çıkan arasındaki fark, giren debinin binde birinden küçük olmalı; enerji için de aynısı. Bu denge kapanmıyorsa geri kalan her sayı şüphelidir.
Salınım her zaman kusur değildir
İzlenen büyüklük düzenli bir genlikle salınıyorsa bu sayısal bir kusur değil, fiziğin kendisi olabilir: girdap kopması, kanat geçiş etkisi, çalkalanma. Bu durumda kalıcı rejim çözücüsüyle ısrar etmek hatadır; problem zamana bağlıdır ve öyle çözülmelidir.
Ayırt etmenin pratik yolu: ağı sıklaştırın ve gevşetme katsayılarını düşürün. Salınım genliği küçülüyorsa sayısaldır, değişmiyorsa fizikseldir.
Rapora ne yazılmalı
Artık eğrileri tek başına kanıt değildir. İzlenen büyüklüklerin iterasyon grafiği, kütle ve enerji denge tablosu, yakınsama ölçütünün ne olduğu yazılmalıdır.
Kaynaklar
- Ferziger, J. H., Perić, M. & Street, R. L. (2020). Computational Methods for Fluid Dynamics, 4th ed. Springer. DOI: 10.1007/978-3-319-99693-6
- Celik, I. B., Ghia, U., Roache, P. J., Freitas, C. J., Coleman, H. & Raad, P. E. (2008). Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. Journal of Fluids Engineering, 130(7), 078001. DOI: 10.1115/1.2960953