2d'de Euler denklemleri


10

Üniversitede ödev olarak 1d simülasyon yaptım. Sorun ifadesi, çalışma sıvısı olarak sıkıştırılabilir ideal gazı içeren 1d şok tüpü problemini çözmekti. Bu problem için, Roe Riemann çözücüsünü kullanarak Eulers denklemleri sistemini çözdüm. Bilmek istiyorum, Euler denklemlerini 2 veya 3 boyutta çözmek için nereden başlamalıyım? Test problemi hangisi, önce düşünmeliyim? (Lütfen ticari çözücüler önermeyin. Kendi kodumu yazmak istiyorum) Sadece kendi kodumu yazmak için biraz yardıma ihtiyacım var.

2D problemini en pratik şekilde tanıtan iyi kaynaklar nelerdir?


Son zamanlarda Euler denklemlerinin dönme değişmezlik özelliğini öğrendim, bu 2d kodlama için kullanılabilir mi? Böyle bir ilk koordinat sistemini bir yüze dik olarak hizalayıp 1d problemini çözdükten sonra çerçeveyi dik yüze dik olarak döndürüp tekrarlıyoruz?
Subodh

2
LeVeque'nin FV kitabının 18-21 bölümlerini öneriyorum: depts.washington.edu/clawpack/book.html Bunun tek dezavantajı, yaklaşımın akılardan ziyade dalgalanmalara dayanmasıdır (ikincisi daha yaygındır).
David Ketcheson

Yanıtlar:


5

David Ketcheson'un da belirttiği gibi, LeVeque'in kitabı harika bir kaynaktır; ancak, doğal olarak CLAWPACK'te kullanılan Riemann çözme yöntemlerine odaklanır. Bu yöntemlere bir alternatif, doğal olarak uygulanması ve 2-D veya 3-D'ye (boyutsal bölme kullanmak yerine) uzatılması kolay olan salınımlı olmayan merkezi şemalardır (örneğin Kurganov ve Tadmor'unki ). Http://www.cscamm.umd.edu/centpack/ adresinde açık kaynak kodu ve büyük kağıt kaynağı var .

Bir kağıt (Bilimsel Bilgi İşlem 25 (3), 2003, 995-1017 ile SIAM Journal) Liska ve Wendroff 1-D, 2-D Euler denklemleri için çeşitli programları (Riemann çözme ve merkezi hem de) karşılaştırır, ve var faydalı olabilecek test problemlerinin sayısı.


4

Evet, 2D problemi aşağı yukarı X yönünde 1D problemini ve ardından Y yönünde 1D problemini çözer. Veri yapıları biraz daha karmaşıktır, Karaca matrisinin ikinci yönde karakteristik projeksiyon için ayarlanması gerekir ve CFL durumunuz da değiştirilmelidir.

Clawpack ve özellikle Clawpack'ın 2D örnekleri (http://depts.washington.edu/clawpack/users-4.6/claw/doc/gallery/gallery_2d.html) başka birinin koduna bakmak istiyorsanız sizin için yararlı olabilir .

Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.