A structure-preserving finite element method for compressible ideal and resistive MHD

05/03/2021 ∙ by Evan S. Gawlik, et al. ∙ 0

We construct a structure-preserving finite element method and time-stepping scheme for compressible barotropic magnetohydrodynamics (MHD) both in the ideal and resistive cases, and in the presence of viscosity. The method is deduced from the geometric variational formulation of the equations. It preserves the balance laws governing the evolution of total energy and magnetic helicity, and preserves mass and the constraint divB = 0 to machine precision, both at the spatially and temporally discrete levels. In particular, conservation of energy and magnetic helicity hold at the discrete levels in the ideal case. It is observed that cross helicity is well conserved in our simulation in the ideal case.



There are no comments yet.


page 17

page 18

page 19

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.