Como construir volume finito bem equilibrado e métodos descontínuos de Galerkin para EDPs hiperbólicas com termos de origem?

13

Termos de origem, como os devidos à batimetria nas equações de águas rasas, precisam ser integrados de uma maneira especial para preservar os estados físicos estáveis. Existe uma maneira geral de construir métodos bem equilibrados ou requer técnicas especiais para cada equação?

Jed Brown
fonte

Respostas:

7

A resposta curta é: requer trabalho específico para diferentes equações, mas existem algumas técnicas gerais que sugerem como fazê-lo. Essencialmente, dada uma evolução de primeira ordem PDE

vocêt=UMAvocê+Bvocê

Onde UMA,B como alguns operadores (possivelmente diferenciais), os estados estacionários são aqueles para os quais

UMAvocê+Bvocê=0

É comum usar uma abordagem de divisão em que UMA e Bsão discretizados de maneiras diferentes. Em seguida, haverá erros de truncamento associados a cada uma dessas discretizações, e os erros de truncamento geralmente não serão cancelados, mesmo no caso de um estado estacionário. Um exemplo clássico (como mencionado na pergunta) são as equações de águas rasas com batimetria, nas quaisUMA representa termos convectivos e Brepresenta a força do momento devido à variável altura do fundo. Existem muitos artigos publicados nos últimos anos que oferecem maneiras diferentes de manter exatamente as soluções em estado estacionário.

Uma abordagem que eu gosto é o uso de solucionadores de onda F de Riemann, conforme proposto por Bale et. al. . A ideia é discretizar os termos convectivos com um método Godunov-tipo, mas para subtrair ao largo da contribuição dos outros termos dentro do solucionador de Riemann. Então, no caso do estado estacionário, nenhuma onda é gerada. No entanto, isso exige que os termos convectivos e de origem sejam calculados exatamente (para cancelar exatamente). Isso é possível para as equações de águas rasas, mas mais difícil para muitos outros sistemas.

David Ketcheson
fonte