Resumen de la Tesis:

El objetivo principal de esta tesis doctoral es el análisis matemático y numérico de dos problemas importantes en electromagnetismo; el primero de ellos está relacionado con diversas aplicaciones de la magnetohidrodinámica, mientras que el segundo tiene que ver con el estudio de corrientes inducidas. En lo que se refiere al primer problema, el estudio se centra en la aproximación numérica de los autovalores del operador rotacional, cuyas soluciones se denominan campos de Beltrami y que surgen en diversas áreas de la física. Para llevar a cabo este estudio, primero se analiza matemáticamente el problema espectral, para lo cual se propone una formulación variacional mixta mediante la que se obtiene una caracterización completa de las soluciones del problema de autovalores. Además, se considera una formulación primal la cual resulta ser "equivalente", bajo ciertas hipótesis, a dicho problema. Para la aproximación numérica del problema de autovalores se consideran esquemas de elementos finitos asociados a cada una de estas formulaciones. En ambos casos se obtienen aproximaciones espectrales con orden de convergencia óptimo, las cuales se corroboran mediante ejemplos numéricos. En la segunda parte de la tesis se aborda el análisis matemático y numérico de diversos problemas de corrientes inducidas en régimen transitorio, suponiendo dominios axisimétricos. Los materiales considerados son no lineales y pueden presentar o no histéresis magnética. Para ello, y motivados por las aplicaciones físicas, se consideran dos tipos datos: el primero de ellos se corresponde con una condición de Dirichlet no homogénea en la frontera del dominio (usualmente la intensidad de corriente), mientras que el segundo consiste en suponer conocido el flujo magnético que atraviesa una sección meridional del dominio. En ambos casos, se propone una formulación en términos del campo magnético, y se considera que la relación entre este campo y la inducción magnética está dada bien mediante una función no lineal, o bien mediante un operador de histéresis. Inicialmente se estudia el problema no lineal de corrientes inducidas considerando el flujo magnético como dato. Se demuestra la existencia y unicidad de solución de la formulación variacional correspondiente mediante un resultado abstracto. Para la aproximación numérica se propone una discretización espacial mediante elementos finitos para la cual se demuestran existencia de solución y una estimación de error. El esquema anterior se combina con un esquema de Euler implícito para la discretización temporal y se demuestran estimaciones óptimas de error. A continuación, se analiza el problema de corrientes inducidas con condición de Dirichlet no homogénea. En este caso, la existencia y unicidad de solución se basan en técnicas de discretización temporal, estimaciones a priori y paso al límite mediante compacidad. La aproximación numérica de este problema se estudia considerando un esquema de Euler implícito para la discretización temporal, que posteriormente se combina con un método de elementos finitos en espacio. Al igual que en el problema anterior, se demuestran estimaciones óptimas de error en las normas apropiadas, tanto para la semi-discretización temporal como para el problema completamente discreto. Para ambos problemas se muestran test de convergencia que confirman los resultados teóricos obtenidos. Finalmente, se estudia el problema axisimétrico de corrientes inducidas en el caso en que la relación entre el campo magnético y la indución magnética viene dada mediante un operador de histéresis. Se demuestra la existencia de solución del problema considerando un operador de histéresis general y las distintas condiciones de contorno. Al igual que en el problema sin histéresis, el estudio de la existencia de solución se basa en una discretización implícita del tiempo; este procedimiento de aproximación es utilizado con frecuencia en el análisis de ecuaciones que incluyen operadores con memoria. Para la aproximación numérica, se considera un esquema completamente discreto mediante elementos finitos y Euler implícito, con una elección particular del operador de histéresis dada por el operador de Preisach clásico. Al contrario que en los problemas sin histéresis, el análisis de convergencia del esquema utilizado se realiza únicamente mediante ejemplos numéricos.

Publicaciones Originadas de la Tesis (ISI)

Alfredo BERMúDEZ, Luc DUPRé, Dolores GóMEZ, Pablo VENEGAS: Electromagnetic computations with Preisach hysteresis model. Finite Elements in Analysis and Design, vol. 116, pp. 65-74 (2017).

Alfredo BERMúDEZ, Dolores GóMEZ, Rodolfo RODRíGUEZ, Pablo VENEGAS: Numerical analysis of a transient non-linear axisymmetric eddy current model. Computers & Mathematics with Applications, vol. 70, 8, pp. 1984-2005, (2015).

Rodolfo RODRíGUEZ, Pablo VENEGAS: Numerical approximation of the spectrum of the curl operator. Mathematics of Computation, vol. 83, 286, pp. 553-577, (2014).

Alfredo BERMúDEZ, Dolores GóMEZ, Rodolfo RODRíGUEZ, Pilar SALGADO, Pablo VENEGAS: Numerical solution of a transient non-linear axisymmetric eddy current model with non-local boundary conditions. Mathematical Models and Methods in Applied Sciences, vol. 23, 13, pp. 2495-2521, (2013).