Regularized Maxwell equations and nodal finite elements for electromagnetic field computations in frequency domain
Authors:   , Otín R. (Rubén) 
Published: October 2023

En esta tesis se presenta una alternativa a la formulación clásica de elementos finitos basada en elementos de arista (edge elements) y las ecuaciones de Maxwell curl-curl. Este enfoque alternativo se basa en elementos nodales y las ecuaciones de Maxwell regularizadas. Las ventajas son que, sin añadir incógnitas adicionales (tales como los multiplicadores de Lagrange), ofrece soluciones sin modos espurios y matrices bien acondicionadas. Además, su representación integral tiene en un orden de singularidad inferior (orden 1 en vez de 3), lo que la hace más adecuada para la hibridación con técnicas numéricas basadas en métodos integrales. Por otro lado, surgen nuevas dificultades que no estaban presentes en la formulación clásica. El principal inconveniente es que se obtiene una solución globalmente errónea si el campo electromagnético tiene una singularidad en el dominio del problema. Además, las condiciones de contorno y las discontinuidades de los campos al cambiar de medio son más laboriosas de implementar. Este trabajo explica la manera de superar estas dificultades y demuestra que se pueden obtener soluciones precisas con elementos nodales y la formulación regularizada.

También se presenta ERMES, un código en C++ que es la implementación de la formulación de elementos finitos mencionada en el párrafo anterior. Calcularemos con ERMES los parámetros de scattering de filtros de microondas y la tasa de absorción específica inducida en un cuerpo cuando está expuesto a campos electromagnéticos. ERMES es también la herramienta de cálculo utilizado en dos nuevos modelos numéricos presentados por primera vez en esta tesis. El primer modelo sirve para caracterizar procesos de conformado electromagnético y el segundo para obtener la impedancia de transferencia del apantallamiento de cables coaxiales perforados.

El modelo de conformado electromagnético calcula la fuerza de Lorentz impulsora del proceso y estima a su vez la frecuencia óptima a la que puede obtenerse la deformación máxima de la pieza. La principal ventaja del método es que proporciona una relación explícita entre la capacidad del banco de condensadores y la frecuencia de la descarga, que es un parámetro clave en el diseño de un sistema de conformado electromagnético. El uso de la formulación regularizada en esta aplicación pone de manifiesto su excelente comportamiento en problemas de baja frecuencia (cuasi-estáticos).

El segundo modelo presentado en este trabajo calcula la impedancia de transferencia de los apantallamientos perforados de cables coaxiales. Este modelo reproduce el comportamiento en alta frecuencia de la impedancia de transferencia con mayor precisión que los métodos encontrados en la literatura y, además, es capaz de analizar una mayor variedad de geometrías y materiales.

In this thesis, we present an alternative approach to the usual finite element formulation based on edge elements and double-curl Maxwell equations. This alternative approach is based on nodal elements and regularized Maxwell equations. The advantages are that, without adding extra unknowns (such as Lagrange multipliers), it provides spurious-free solutions and well-conditioned matrices. Besides, its integral representation involves a less singular kernel (order 1 instead of 3), which makes this approach best suited to hybridization with integral numerical techniques. On the other hand, new difficulties arise that were not present in the classical formulation. The main drawback is that a globally wrong solution is obtained when the electromagnetic field has a singularity in the problem domain. Also, boundary conditions and field discontinuities are more laborious to implement. This work explains how to overcome these difficulties and demonstrates that accurate solutions can be obtained with nodal elements and the regularized formulation.

We also present ERMES, the C++ implementation of the finite element approach depicted above and the main deliverable of this work. We compute with ERMES the scattering parameters of microwave filters and the specific absorption rate induced in a body when exposed to electromagnetic fields. ERMES is also the computational tool used in two novel numerical models introduced in this thesis. The first one characterizes electromagnetic metal forming processes and the second one the transfer impedance of cable shields.

The electromagnetic metal forming model calculates the driving Lorentz force and estimates the optimum frequency at which it is attained the maximum workpiece deformation. The main advantage of the approach is that it provides an explicit relation between the capacitance of the capacitor bank and the frequency of the discharge, which is a key parameter in the design of an electromagnetic forming system. The use of the regularized formulation in this application reveals its excellent behavior in the low-frequency (quasi-static) regime.

The second numerical model introduced in this work computes the transfer impedance of cable shields. The model reproduces the high frequency behavior of the transfer impedance more accurately than the approaches found in the literature and, moreover, it is able to analyze a wider variety of geometries and materials.