Solucionador FEM 3D de código abierto para electromagnetismo (Time-Harmonic Maxwell)

8

Me preguntaba si existe algún solucionador de FEM de código abierto bueno (preciso / rápido / fácil de usar) para las ecuaciones de Maxwell 3D armónicas en el tiempo. Estoy buscando simular sistemas de unas pocas longitudes de onda grandes en las dimensiones X / Y y decenas de longitudes de onda en la dimensión Z, por lo que problemas de tamaño relativamente pequeño. Todos mis dieléctricos son homogéneos y perfectamente rectangulares (cajas).

Hice algunas búsquedas en línea y encontré varios paquetes FEM de código abierto maduros, pero no parecen admitir Maxwell. ¡Cualquier sugerencia sería muy apreciada!

Costis
fuente
¿Probaste FEniCS?
facetus

Respuestas:

6

deal.II (ver http://www.dealii.org/ ) admite elementos Nedelec y, como consecuencia, puede resolver los problemas que le interesan. (Descargo de responsabilidad completo: soy uno de los principales desarrolladores de trato .II.)

Wolfgang Bangerth
fuente
¡Gracias, profesor Bangerth! ¿Deal.II tiene algún soporte incorporado para materiales PML o necesitaría codificarse por separado?
Costis
deal.II no implementa ninguna ecuación o formulación en particular, solo proporciona las herramientas para hacerlo. Por lo tanto, en su caso, deberá proporcionar la forma bilineal que desea resolver. Si eso incluye la PML o cualquier otra condición límite absorbente, depende de usted.
Wolfgang Bangerth
¿Deal.ii admite números complejos hoy? Creo que es útil porque resolver un sistema lineal complejo es mucho más natural que resolver el sistema real equivalente de doble tamaño.
Hui Zhang
No, todavía usa dos sistemas de valor real. Es demasiado difícil asegurarse de que uno distinga entre transposición y Hermitian, los productos de punto regulares y complejos, etc.
Wolfgang Bangerth
4

Hypre tiene varios preacondicionadores integrados para resolver las ecuaciones de Maxwell. Hay varios paquetes que lo conectan (puede usar hypre de PETSc) como solucionador para sistemas algebraicos lineales, pero también tiene una cuadrícula estructurada y una interfaz de elementos finitos.

Daniel Shapero
fuente