Viscous regularization of the MHD equations

Tuan Anh Dao,Lukas Lundgren,Murtazo Nazarov
2024-02-06
Abstract:Nonlinear conservation laws such as the system of ideal magnetohydrodynamics (MHD) equations may develop singularities over time. In these situations, viscous regularization is a common approach to regain regularity of the solution. In this paper, we present a new viscous flux to regularize the MHD equations which holds many attractive properties. In particular, we prove that the proposed viscous flux preserves positivity of density and internal energy, satisfies the minimum entropy principle, is consistent with all generalized entropies, and is Galilean and rotationally invariant. We also provide a variation of the viscous flux that conserves angular momentum. To make the analysis more useful for numerical schemes, the divergence of the magnetic field is not assumed to be zero. Using continuous finite elements, we show several numerical experiments including contact waves and magnetic reconnection.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the singularity problem that may occur in the numerical simulation of the ideal magnetohydrodynamics (MHD) equations. Specifically, the ideal MHD equations may develop singularities in certain cases, resulting in unstable or inaccurate numerical solutions. To overcome these problems, the author proposes a new viscous regularization method to regularize the MHD equations. ### Specific Problems and Solutions 1. **Singularity Problem**: - The ideal MHD equations are a system of nonlinear conservation laws, and these equations may produce singularities in certain cases, causing the numerical solutions to lose smoothness and physical meaning. - For example, the fluid density or internal energy becoming negative will lead to the failure of physical laws, which is unacceptable. 2. **Limitations of Existing Methods**: - **Resistive - type MHD Flux**: Although the resistive term can be introduced to regularize the equations, this method has some drawbacks, such as not satisfying the minimum entropy principle. - **Global Regularization**: The global regularization method also has limitations. For example, it is not compatible with the divergence - free property of the magnetic field, cannot preserve angular momentum conservation, and adds equal viscosity to all conservative variables. 3. **Newly Proposed Viscous Regularization Method**: - The author proposes a new method that combines the Guermond - Popov (GP) viscous flux and the resistive - type electromagnetic flux. - This new method has the following advantages: - Maintains the positivity of density and internal energy. - Satisfies the minimum entropy principle. - Is compatible with all generalized entropies. - Has Galilean invariance and rotational invariance. - Provides a variant to preserve angular momentum conservation. ### Numerical Experiment Verification To verify the effectiveness of the new method, the author carried out a number of numerical experiments using the continuous finite - element method, including the simulation of contact waves and magnetic reconnection phenomena. The results of these experiments show that the newly proposed viscous regularization method can improve the stability and accuracy of the numerical solutions while maintaining the physical characteristics. ### Conclusion By introducing the new viscous regularization method, the author has successfully solved the singularity problem that may occur in the numerical simulation of the ideal MHD equations and ensured the physical rationality and stability of the numerical solutions. This provides a more reliable basis for future research.