Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics

This article serves as a summary outlining the mathematical entropy analysis of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD equations as they are particularly useful for mathematically modeling a wide variety of magnetized fluids. In order to be self-contained we first motivate the physical properties of a magnetic fluid and how it should behave under the laws of thermodynamics. Next, we introduce a mathematical model built from hyperbolic partial differential equations (PDEs) that translate physical laws into mathematical equations. After an overview of the continuous analysis, we thoroughly describe the derivation of a numerical approximation of the ideal MHD system that remains consistent to the continuous thermodynamic principles. The derivation of the method and the theorems contained within serve as the bulk of the review article. We demonstrate that the derived numerical approximation retains the correct entropic properties of the continuous model and show its applicability to a variety of standard numerical test cases for MHD schemes. We close with our conclusions and a brief discussion on future work in the area of entropy consistent numerical methods and the modeling of plasmas

Electromagnetic induction in non-uniform domains

Kinematic simulations of the induction equation are carried out for different setups suitable for the von-K\'arm\'an-Sodium (VKS) dynamo experiment. Material properties of the flow driving impellers are considered by means of high conducting and high permeability disks that are present in a cylindrical volume filled with a conducting fluid. Two entirely different numerical codes are mutually validated by showing quantitative agreement on Ohmic decay and kinematic dynamo problems using various configurations and physical parameters. Field geometry and growth rates are strongly modified by the material properties of the disks even if the high permeability/high conductivity material is localized within a quite thin region. In contrast the influence of external boundary conditions remains small. Utilizing a VKS like mean fluid flow and high permeability disks yields a reduction of the critical magnetic Reynolds number for the onset of dynamo action of the simplest non-axisymmetric field mode. However this decrease is not sufficient to become relevant in the VKS experiment. Furthermore, the reduction of Rm_c is essentially influenced by tiny changes in the flow configuration so that the result is not very robust against small modifications of setup and properties of turbulence

High-Resolution Mathematical and Numerical Analysis of Involution-Constrained PDEs

Partial differential equations constrained by involutions provide the highest fidelity mathematical models for a large number of complex physical systems of fundamental interest in critical scientific and technological disciplines. The applications described by these models include electromagnetics, continuum dynamics of solid media, and general relativity. This workshop brought together pure and applied mathematicians to discuss current research that cuts across these various disciplines’ boundaries. The presented material illuminated fundamental issues as well as evolving theoretical and algorithmic approaches for PDEs with involutions. The scope of the material covered was broad, and the discussions conducted during the workshop were lively and far-reaching

Robust Preconditioners for Incompressible MHD Models

In this paper, we develop two classes of robust preconditioners for the structure-preserving discretization of the incompressible magnetohydrodynamics (MHD) system. By studying the well-posedness of the discrete system, we design block preconditioners for them and carry out rigorous analysis on their performance. We prove that such preconditioners are robust with respect to most physical and discretization parameters. In our proof, we improve the existing estimates of the block triangular preconditioners for saddle point problems by removing the scaling parameters, which are usually difficult to choose in practice. This new technique is not only applicable to the MHD system, but also to other problems. Moreover, we prove that Krylov iterative methods with our preconditioners preserve the divergence-free condition exactly, which complements the structure-preserving discretization. Another feature is that we can directly generalize this technique to other discretizations of the MHD system. We also present preliminary numerical results to support the theoretical results and demonstrate the robustness of the proposed preconditioners