Equações de Euler em 2d

10

Como trabalho na faculdade, fiz uma simulação em 1d. A declaração do problema era resolver o problema do tubo de choque 1d envolvendo gás ideal compressível como fluido de trabalho. Para esse problema, resolvi o sistema de equações de Eulers usando o solucionador Riemann de Roe. Quero saber, para resolver as equações de Euler em 2 ou 3 dimensões, por onde devo começar? Qual é o problema do teste, devo considerar primeiro? (Por favor, não sugira solucionadores comerciais. Quero escrever meu próprio código) só preciso de ajuda para escrever meu próprio código.

Quais são os bons recursos que introduzem o problema 2d da maneira mais prática?

Subodh
fonte
Recentemente eu aprendi a propriedade de invariância rotacional das equações de Euler, isso pode ser usado para a codificação 2D? Como o primeiro, alinhamos o sistema de coordenadas perpendicularmente a uma face e resolvemos o problema 1d, depois giramos o quadro perpendicular à face ortogonal e repetimos?
Subodh
2
Eu recomendo os capítulos 18-21 do livro de FV de LeVeque: depts.washington.edu/clawpack/book.html A única desvantagem disso é que a abordagem é baseada em flutuações e não em fluxos (os últimos são mais comuns).
David Ketcheson

Respostas:

5

Como David Ketcheson sugere, o livro de LeVeque é um grande recurso; no entanto, concentra-se naturalmente nos métodos de solução de Riemann usados ​​no CLAWPACK. Uma alternativa a esses métodos são os esquemas centrais não-oscilatórios (por exemplo, o de Kurganov & Tadmor ), que são fáceis de implementar e se estendem naturalmente a 2-D ou 3-D (em vez de usar a divisão dimensional). Há código-fonte aberto e grande variedade de documentos em http://www.cscamm.umd.edu/centpack/ .

Um artigo de Liska & Wendroff (SIAM Journal on Scientific Computing 25 (3), 2003, 995-1017) compara vários esquemas (tanto de resolução de Riemann quanto centrais) para as equações de Euler em 1-D e 2-D, e possui um número de problemas de teste que podem ser úteis.

Chris Johnson
fonte
4

Sim, o problema 2D mais ou menos se torna resolver um problema 1D na direção X seguido por um problema 1D na direção Y. As estruturas de dados são um pouco mais complicadas, a matriz Roe precisa ser ajustada para a projeção característica na segunda direção, e sua condição de CFL também precisa ser modificada.

O Clawpack e, em particular, os exemplos 2D do Clawpack (http://depts.washington.edu/clawpack/users-4.6/claw/doc/gallery/gallery_2d.html) podem ser úteis para você, se você quiser espiar o código de outra pessoa .

Rhys Ulerich
fonte