(196b) Fully-Implicit Finite Element Formulations for Stability and Bifurcation Analysis of Large-Scale Resistive Magneto-Hydrodynamic Systems

Ionized fluids with strong electromagnetic effects occur frequently in nature and are critical for many important technological applications. Examples include stellar interiors, gaseous nebula, the earth's magnetoshpere, and Tokamak and Z-pinch physics. These systems are described by a set of partial differential equations that conserve momentum, mass, charge, and energy along with magnetic flux for the electric and magnetic fields (Maxwell's equations). The resulting magnetohydrodynamics (MHD) equations are strongly coupled, highly nonlinear, and span a large range of time and length scales, making the scalable, robust, and accurate solution of such systems extremely challenging.

In this presentation, we will discuss the development of multiple MHD formulations based on unstructured stabilized finite element methods. The formulations are designed for weak enforcement of the solenoidal constraint (divergence?free magnetic field). The formulations include a 2D vector potential formulation, a 3D projection formulation and a 3D Lagrange multiplier formulation using variational multi-scale stabilization. The resulting set of nonlinear equations are solved using a fully-coupled Newton-Krylov solver with nonlinear globalization techniques. Linear systems are solved using a multi-level preconditioned GMRES iterative technique. We will present numerical performance, accuracy, and scalability studies of the formulations. Additionally we will present an application of our solvers to perform a stability and bifurcation analysis of the hydromagnetic Rayleigh-Bernard problem. The bifurcation analysis uses the minimally augmented bifurcation formulation that has been extended to large-scale systems.

*This work was partially funded by the DOE Office of Science AMR Program, and was carried out at Sandia National Laboratories operated for the U.S. Department of Energy under contract no. DE-ACO4-94AL85000