A Finite Element Method for MHD that Preserves Energy, Cross-Helicity, Magnetic Helicity, Incompressibility, and div B = 0
We construct a structure-preserving finite element method and time-stepping scheme for inhomogeneous, incompressible magnetohydrodynamics (MHD). The method preserves energy, cross-helicity (when the fluid density is constant), magnetic helicity, mass, total squared density, pointwise incompressibility, and the constraint div B = 0 to machine precision, both at the spatially and temporally discrete levels.
READ FULL TEXT