Saya mensimulasikan aliran yang tidak dapat dimampatkan di atas silinder pada bilangan Reynold 500. Saya memecahkan persamaan nokes stokes menggunakan metode koreksi tekanan. Solusi saya menjadi tidak stabil setelah waktu tertentu (sekitar 5 detik).
Saya telah mencoba memperbaiki mesh saya, stepsize (0,05) (memastikan CFL saya <1, meskipun saya menggunakan metode implisit)
Kondisi batas saya, mesh dan hasil yang tidak stabil ditunjukkan pada gambar terlampir. Domain ini sekitar 25 kali lebih besar dari diameter silinder.
Saya telah mencoba mensimulasikan masalah ini O grid (yang menjadi tidak stabil segera).
Tautan berikut berisi gambar-gambar kondisi batas dan hasil.
Saya akan berterima kasih jika ada yang bisa berbagi pemikiran / pengalaman mereka tentang masalah ini. Terimakasih banyak.
diedit:
Permintaan maaf untuk kesalahan pengetikan:
Saya menggunakan kondisi batas berikut: Batas Neumann
pada Dirichlet Boundary
diedit:
saya telah menerapkan kondisi batas kecepatan pada node di sekitar batas dirichlet. Node sudut kanan-atas dan kanan-bawah adalah batas dirichlet dengan kecepatan 1.
Setelah itu, saya melihat lebih dalam pada hasil simulasi, saya perhatikan bahwa ketidakstabilan mulai merayap di persimpangan inflow / outflow.
sumber
Jawaban:
Saya sudah menemukan masalahnya. Saya harus meningkatkan ukuran domain lebih jauh untuk menghilangkan efek batas. Selain itu, saya harus mengurangi angka CFL menjadi sekitar 0,5-1,0
Saya pikir nomor CFL perlu dikurangi lebih lanjut untuk nomor reynolds yang lebih tinggi.
Awalnya, saya berpikir bahwa saya telah mengurangi ukuran langkah yang cukup, tetapi bukan itu masalahnya.
sumber