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

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:

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

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