O que é um pré-condicionador escalável para Helmholtz de alta frequência?


15

Os métodos padrão de decomposição multigrid e de domínio não funcionam, mas eu tenho grandes problemas em 3D e os solucionadores diretos não são uma opção. Quais métodos devo tentar?

Como minhas escolhas são afetadas pelas seguintes considerações?

  • os coeficientes variam em várias ordens de magnitude, ou
  • elemento finito versus diferentes métodos finitos são usados

11
Em 3D, os solucionadores iterativos normalmente apresentam desempenho ruim, recomendo consultar alguns solucionadores diretos de reordenação da matriz HSS de Ming Gu, Xia e Chandrasekaran.
Shuhao Cao

Respostas:



5

Penso que, em geral, vale lembrar que os métodos mais eficientes que temos (multigrid geométricos e algébricos e, até certo ponto, decomposição de domínio) dependem do fato de que as soluções de EDPs geralmente são suaves e que a solução de um problema mais grosseiro pode gerar uma boa aproximação para o problema de escala fina. O problema com a equação de Helmholtz para altas frequências é que essa suposição não é verdadeira: você precisa de uma malha relativamente fina para representar a solução, e os solucionadores de malha grossa não conseguirão produzir nada que seja de grande utilidade. Consequentemente, as abordagens típicas para bons pré-condicionadores não funcionam nesse caso, e essa é a razão subjacente pela qual não há boas opções reais no seu caso, além de apenas jogar muitos processadores no problema;


4

O material da matriz H de Jack Poulson e Lexing Ying é o método mais eficiente que conheço. Isso deve ser divulgado na primavera, mas eles fizeram apresentações.


2
Provavelmente, eu deveria qualificar sua afirmação dizendo que resolvemos problemas de alta frequência em larga escala com uma abordagem PML em movimento eficientemente em milhares de núcleos para frequência muito alta, mas ainda não testamos a estratégia da matriz H nessa escala. O raciocínio é que ele não tem tanta justificativa teórica em 3d, apesar de ser mais escalável do ponto de vista da computação paralela.
precisa
Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.