Kısa bir süre önce ele alınan Darpa Suboff gövdesinin 2 Boyutlu direnç çözümlerinin devamı niteliğinde eşdeğer modelleme ve çözüm yaklaşımı kullanılmak suretiyle bu kez 3 Boyutlu bir değerlendirmeden bahsedilecek.

Suboff çıplak teknesinin yarım gövdesinden oluşan 3B modeli temelinde OpenFOAM ile elde edilen direnç değerlerinin 2B çözümler ile karşılaştırılması ve bir miktar renkli resimden müteşekkil bazı son işlem görüntülerinin sunulmasıyla da mevzu kapatılacak.

Önceki 2B inceleme ile yapılacak karşılaştırmanın tutarlı olabilmesi için söz konusu çalışma için tercih edilen eksensimetrik Suboff modeli örgüsünün boyuna eksen etrafında π kadar döndürülmesiyle elde edilen yarım gövde 3B model elde edilmiştir ve dönüş yönünde 48 düğüm noktası kullanılmıştır. Diğer bütün değişkenler ise önceki çalışma ile aynıdır.

Süre kısıtlaması sebebiyle bu kez yalnızca yArtı ~30 özelliğindeki sınır tabaka örgüsü ile daha önce en iyi sonucu veren türbülans modellerinden ikisi olan kOmegaTNT ve LienCubicKE ile üçüncü olarak da bu tür problemlerin çözümünde kullanılan en yaygın genel maksat modeli olan kOmegaSST tercih edilmiştir.

Bu durumda Reynolds sayısı olan 14e06 için gerçekleştirilecek hesaplamalarda kullanılan temel değişkenler [1] aşağıdaki gibidir:

• Ortam: Tatlı su
• Sıcaklık (t): 15°C
• Yoğunluk (ρ): 999,1026kg/m³
• Girdap luzûciyeti (ν): 1,1386e-06m²/sn
• Serbest türbülans yoğunluğu (I): %0,161
• Hız (U): 3,660m/sn
• Uzunluk (L): 4,3561m
• Islak Alan (S): 5,998m2 (tam model için!)

Resim.1) Suboff DTRC model 5470 için 3 Boyutlu çokyüzlü melez örgü ve akış hatları.
Yanlış hatırlamıyorsan burada hücum açısı 30derece.

Yukarıdaki resimde görülen söz konusu 3B hesaplama örgüsü, önceki bölümde verilen Suboff çıplak tekne kesit geometrisi kullanılarak kolayca imâl edilebilir. Eğer bu örgüyü Gmsh ile üretmek isterseniz bu geometri verisini doğrudan kullanarak Gmsh ile Eksensimetrik Melez Örgü İmâli adlı bölümde açıklandığı şekilde oluşturup OpenFOAM içine ithâl ederek çözebilirsiniz.

Çizelge.1) 3B çözüm ile elde edilen direnç sonuçlarının (yArtı ~30) havuz deneyi sonuçları [2] ve ITTC-1957 ile karşılaştırılması.

3B model kullanılmasıyla bu kez 0 dereceye ilave olarak 5, 10, 18 ve 30 derecelik açılar için de çözümler elde edilebilmiştir. Bununla birlikte konuya yeterli zaman ayrılamaması sebebiyle bu ilave hücum açıları ile elde edilen sonuçların deney verileriyle karşılaştırılması gerçekleştirilememiştir yine de aşağıdaki birkaç resim üzerinde 18derece için sonuçlara ait bir iki küçük karşılaştırma görülebilir...

Resim.2) Suboff DTRC model 5470 için 3 farklı türbülans modeli ile 18derece için elde edilen duvar kayma gerilmesi sonuçlarının rüzgâraltı tarafından görünüşü.

Resim.3) Suboff DTRC model 5470 için kuyruk bölgesindeki bir kesitte (x/L=0,95) hesaplanan akış hızları.

Resim.4) DTRC Model 5470 kestileri üzerinde hesaplanan basınç dağılımları.
Üstte x/L=0,5 ve altta x/L=0,95.
X ekseninin solu rüzgârüstü sağı ise rüzgâraltı tarafıdır.
Daha önce benzer bir çalışma ile gerçekleştirilen doğrulama denemesinde en iyi basınç dağılımını LienCubicKE sağlamıştı.

Resim.5) DTRC Model 5470 Kesitleri üzerinde hesaplanan bileşke kayma gerilmesi dağılımları.
Üstte x/L=0,5 ve altta x/L=0,95.
X ekseninin solu rüzgârüstü sağı ise rüzgâraltı tarafıdır.