Yüksek frekanslı Helmholtz için ölçeklenebilir bir ön koşul nedir?


15

Standart multigrid ve etki alanı ayrıştırma yöntemleri işe yaramaz, ancak büyük 3D sorunlarım var ve doğrudan çözücüler bir seçenek değil. Hangi yöntemleri denemeliyim?

Seçimlerim aşağıdaki hususlardan nasıl etkileniyor?

  • katsayılar birkaç büyüklük sırasına göre değişir veya
  • Sonlu elemanlara karşı sonlu farklı yöntemler kullanılır

1
3D'de, yinelemeli çözücüler normal olarak kötü performans gösterir, Ming Gu, Xia ve Chandrasekaran'ın doğrudan çözücülerini yeniden sıralayan bazı HSS matrisine başvurmanızı öneririm.
Shuhao Cao

Yanıtlar:



5

Genel olarak, sahip olduğumuz en verimli yöntemlerin (geometrik ve cebirsel multigridin yanı sıra bir dereceye kadar alan ayrışması) PDE'lerin çözümlerinin genellikle pürüzsüz olduğuna ve daha kaba bir problemin çözülmesinin ince ölçek problemi için iyi yaklaşım. Yüksek frekanslar için Helmholtz denklemi ile ilgili sorun, bu varsayımın doğru olmamasıdır: çözümü temsil etmek için nispeten ince bir ağa ihtiyacınız vardır ve kaba ağ çözücüler, çok kullanışlı bir şey üretemeyecektir. Sonuç olarak, iyi önkoşullara yönelik tipik yaklaşımlar bu durumda işe yaramamaktadır ve durumunuzda çok fazla işlemciyi soruna atmaktan kısa bir süre içinde gerçek iyi seçeneklerin bulunmamasının altında yatan neden budur;


4

Jack Poulson ve Lexing Ying'den H-matrix malzemeleri bildiğim en etkili yöntem. Bu baharda serbest bırakılmalıdır, ancak üzerinde sunum yaptılar.


2
Muhtemelen çok yüksek frekans için binlerce çekirdek üzerinde hareketli bir PML yaklaşımı ile büyük ölçekli yüksek frekanslı problemleri verimli bir şekilde çözdüğümüzü söyleyerek ifade etmeliyim , ancak bu ölçekte H-matris stratejisini henüz test etmedik. Akıl yürütme, paralel bilgi işlem perspektifinden daha ölçeklenebilir olmasına rağmen, 3D'de çok fazla teorik gerekçeye sahip olmamasıdır.
Jack Poulson
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.