Solucionador 3D FEM de código aberto para eletromagnetismo (Maxwell harmônico no tempo)

8

Eu queria saber se existe algum bom (preciso / rápido / fácil de usar) resolvedores FEM de código aberto para as equações de Maxwell com harmônicas no tempo 3D. Estou procurando simular sistemas com alguns comprimentos de onda grandes nas dimensões X / Y e dezenas de comprimentos de onda na dimensão Z, portanto, problemas de tamanho relativamente pequeno. Todos os meus dielétricos são homogêneos e perfeitamente retangulares (caixas).

Pesquisei on-line e encontrei vários pacotes FEM maduros de código aberto, mas eles não parecem suportar o Maxwell. Qualquer sugestão será muito bem-vinda!

Costis
fonte
Você experimentou o FEniCS?
Facetus 28/05

Respostas:

6

deal.II (consulte http://www.dealii.org/ ) suporta elementos Nedelec e, como conseqüência, pode resolver os problemas nos quais você está interessado. (Isenção de responsabilidade completa: sou um dos principais desenvolvedores de negócios .II.)

Wolfgang Bangerth
fonte
Obrigado, professor Bangerth! O deal.II tem algum suporte interno para materiais de PML ou precisaria ser codificado separadamente?
Costis 29/05
deal.II não implementa nenhuma equação ou formulação específica - apenas fornece as ferramentas para isso. Portanto, no seu caso, você precisará fornecer a forma bilinear que deseja resolver. Depende de você incluir a PML ou qualquer outra condição de contorno de absorção.
Wolfgang Bangerth
O deal.ii suporta números complexos hoje? Eu acho que é útil porque resolver um sistema linear complexo é muito mais natural do que resolver o sistema real equivalente de tamanho duplo.
Hui Zhang
Não, ele ainda usa dois sistemas com valor real. É muito difícil garantir uma distinção entre transposição e Hermitian, produtos regulares e complexos, etc.
Wolfgang Bangerth
4

A Hypre possui vários pré-condicionadores embutidos para resolver as equações de Maxwell. Existem vários pacotes que fazem interface com ele (você pode usar o hypre do PETSc) como um solucionador de sistemas algébricos lineares, mas também possui uma grade estruturada e uma interface de elementos finitos também.

Daniel Shapero
fonte