Stabilized finite element approximation of the incompressible MHD equations

Abstract

No es frecuente encontrar un campo donde dos ramas principales de la Física estén involucradas. La Magnetohidrodinámica es uno de tales campos debido a que involucra a la Mecánica de Fluidos y al Electromagnetismo. Aun cuando puede parecer que esas dos ramas de la Física tienen poco en común, comparten similitudes en las ecuaciones que gobiernan los fenómenos involucrados en ellas. Las ecuaciones de Navier-Stokes y las ecuaciones de Maxwell, ambas en la raíz de la Magnetohidrodinámica, tienen una condición de divergencia nula y es esta condición de divergencia nula sobre la velocidad del fluido y el campo magnético lo que origina algunos de los problemas numéricos que surgen en la modelación de los fenómenos donde el flujo de fluidos y los campos magnéticos están acoplados.<br/>El principal objetivo de este trabajo es desarrollar un algoritmo eficiente para la resolución mediante elementos finitos de las ecuaciones de la Magnetohidrodinámica de fluidos incompresibles.<br/>Para lograr esta meta, los conceptos básicos y las características de la Magnetohidrodinámica se presentan en una breve introducción informal.<br/>A continuación, se da una revisión completa de las ecuaciones de gobierno de la Magnetohidrodinámica, comenzando con las ecuaciones de Navier-Stokes y las ecuaciones de Maxwell. Se discute la aproximación que da origen a las ecuaciones de la Magnetohidrodinámica y finalmente se presentan las ecuaciones de la Magnetohidrodinámica.<br/>Una vez que las ecuaciones de gobierno de la Magnetohidrodinámica han sido definidas, se presentan los esquemas numéricos desarrollados, empezando con la linealización de las ecuaciones originales, la formulación estabilizada y finalmente el esquema numérico propuesto. En esta etapa se presenta una prueba de convergencia.<br/>Finalmente, se presentan los ejemplos numéricos desarrollados durante este trabajo.<br/>Estos ejemplos pueden dividirse en dos grupos: ejemplos numéricos de comparación y ejemplos de internes tecnológico. Dentro del primer grupo están incluidas simulaciones del flujo de Hartmann y del flujo sobre un escalón. El segundo grupo incluye simulaciones del flujo en una tobera de inyección de colada continua y el proceso Czochralski de crecimiento de cristales. It is not frequent to find a field where two major branches of Physics are involved. Magnetohydrodynamics is one of such fields because it involves Fluid Mechanics and Electromagnetism. Although those two branches of Physics can seem to have little in common, they share similarities in the equations that govern the phenomena involved. The Navier-Stokes equations and the Maxwell equations, both at the root of Magnetohydrodynamics, have a divergence free condition and it is this divergence free condition over the velocity of the fluid and the magnetic field what gives origin to some of the numerical problems that appear when approximating the equations that model the phenomena where fluids flow and magnetic fields are coupled.<br/>The main objective of this work is to develop an efficient finite element algorithm for the incompressible Magnetohydrodynamics equations.<br/>In order to achieve this goal the basic concepts and characteristics of Magnetohydrodynamics are presented in a brief and informal introduction.<br/>Next, a full review of the governing equations of Magnetohydrodynamics is given, staring from the Navier-Stokes equations and the Maxwell equations. The MHD approximation is discussed at this stage and the proper Magnetohydrodynamics equations for incompressible fluid are reviewed.<br/>Once the governing equations have been defined, the numerical schemes developed are presented, starting with the linearization of the original equations, the stabilization formulations and finally the numerical scheme proposed. A convergence test is shown at this stage.<br/>Finally, the numerical examples performed while this work was developed are presented. These examples can be divided in two groups: numerical benchmarks and numerical examples of technological interest. In the first group, the numerical simulations for the Hartmann flow and the flow over a step are included. The second group includes the simulation of the clogging in a continuous casting nozzle and Czochralski crystal growth process.