English

A Kernel Based High Order "Explicit" Unconditionally Stable Constrained Transport Method for Ideal Magnetohydrodynamics

Numerical Analysis 2019-08-06 v1 Numerical Analysis

Abstract

The ideal Magnetohydrodynamics (MHD) equations are challenging because one needs to maintain the divergence free condition, \Bv=0\nabla \cdot \Bv = 0. Many numerical methods have been developed to enforce this condition. In this work, we further our work on mesh aligned constrained transport by developing a new kernel based approach for the vector potential in 2D and 3D. The approach for solving the vector potential is based on the method of lines transpose and is A-stable, eliminating the need for diffusion limiters needed in our previous work in 3D. The work presented here is an improvement over the previous method in the context of problems with strong shocks due to the fact that we could eliminate the diffusion limiter that was needed in our previous version of constrained transport. The method is robust and has been tested on the 2D and 3D cloud shock, blast wave and field loop problems.

Keywords

Cite

@article{arxiv.1908.01023,
  title  = {A Kernel Based High Order "Explicit" Unconditionally Stable Constrained Transport Method for Ideal Magnetohydrodynamics},
  author = {Firat Cakir and Andrew Christlieb and Yan Jiang},
  journal= {arXiv preprint arXiv:1908.01023},
  year   = {2019}
}
R2 v1 2026-06-23T10:38:35.219Z