Existem abordagens de divisão de operadores para PDEs multifísicos que alcançam convergência de alta ordem?

16

Dada uma evolução PDE

vocêt=UMAvocê+Bvocê

onde são operadores diferenciais (possivelmente não lineares) que não se deslocam, uma abordagem numérica comum é alternar entre resolverUMA,B

vocêt=UMAvocê

e

vocêt=Bvocê.

A implementação mais simples disso é conhecida como divisão de Godunov e é precisa de 1ª ordem. Outra abordagem bem conhecida, conhecida como divisão de Strang, é de segunda ordem precisa. Existem métodos de divisão de operadores de ordem superior (ou abordagens alternativas de discretização multifísica)?

David Ketcheson
fonte
11
Os termos são rígidos ou não rígidos? Você tem uma função que aplica A e B ou apenas um algoritmo que avança o estado de para t n + 1 ? No caso em que um é rígido e outro não, existem muitos métodos interessantes. tntn+1 1
Jed Brown

Respostas:

7

Entendi que a fórmula BCH era uma maneira sistemática de aproximar o exponencial da matriz de duas matrizes não comutativas.

Matt Knepley
fonte
Mas isso não leva a termos complexos, mesmo quando o PDE é real? As pessoas o usam para discretização superior à 2ª ordem?
David Ketcheson
11
Não da minha memória (ou da página da web). Isso leva a muitos comutadores. No corpo quântico de muitos corpos, há boas maneiras de simplificar essas expressões.
Matt Knepley
7

Se você considera os operadores gerais A e B e se deseja apenas realizar etapas positivas no tempo (que normalmente é necessário para resolver problemas parabólicos), existe uma barreira de ordem de 2, ou seja, usando qualquer tipo de divisão, não é possível obter uma taxa de convergência superior a dois. Uma prova elementar é dada em um artigo recente de S. Blanes e F. Casas, http://www.gicas.uji.es/Fernando/MyPapers/2005APNUM.pdf .

No entanto, existem várias maneiras, se você souber um pouco mais sobre o seu problema:

  • Suponha que você possa resolver suas equações ao contrário no tempo (o que é comum em, por exemplo, equações de Schrödinger), então existem muitas splings disponíveis, consulte o livro "Integração Numérica Geométrica" ​​de Hairer, Lubich e Wanner.
  • Se seus operadores geram semigrupos analíticos, ou seja, você pode inserir valores complexos para t (típico para equações parabólicas), observou-se recentemente que é possível obter divisões de ordem superior entrando no plano complexo. Os primeiros artigos nessa direção são de E. Hansen e A. Ostermann, http://www.maths.lth.se/na/staff/eskil/dataEskil/articles/Complex.pdf , e F. Castella, P. Chartier , S. Descombes e G. Vilmart. A escolha de separações complexas que são "ideais", em certo sentido, é um tópico da pesquisa atual; você pode encontrar vários artigos sobre o tópico no arxiv.

Resumindo: se você colocar algumas suposições sobre o seu problema, poderá obter algo, mas se não, o pedido 2 será o máximo.

PS: Eu tive que retirar o link do Castella et al-paper devido à prevenção de spam, mas você pode encontrá-lo facilmente no google.

Philipp Dörsek
fonte
5

O grupo CCSE da LBNL usou recentemente métodos de Correção Adiada Espectral (SDC) em um fluxo de baixo número de mach com química complexa. Eles comparam os resultados do SDC com a divisão Strang e os resultados são muito promissores.

Aqui está um rascunho do documento com os detalhes: Uma estratégia de acoplamento de correção diferida para baixo fluxo de número Mach com química complexa

Observe que o esquema SDC é um esquema iterativo que converge para uma solução de disposição precisa de alta ordem, mas é construído a partir de métodos de primeira ordem.

Matthew Emmett
fonte
2

O erro de divisão pode, pelo menos em princípio, ser reduzido por métodos de correção diferida espectral. No entanto, essa parece ser uma área de pesquisa ativa e não é realmente algo pronto para uso geral.

Brian
fonte