An adaptive finite element method for magnetohydrodynamics

  • Authors:
  • S. Lankalapalli;J. E. Flaherty;M. S. Shephard;H. Strauss

  • Affiliations:
  • Third Wave Systems, Minneapolis, MN, United States;Scientific Computation Research Center, Rensselaer Polytechnic Institute, Troy, NY 12180-3590, United States;Scientific Computation Research Center, Rensselaer Polytechnic Institute, Troy, NY 12180-3590, United States;Courant Institute of Mathematical Sciences, New York University, NY, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.46

Visualization

Abstract

We describe a procedure for the adaptive h-refinement solution of the incompressible MHD equations in stream function form using a stabilized finite element formulation. The mesh is adapted based on a posteriori spatial error estimates of the magnetic field using both recovery and order extrapolation techniques. The step size for time integration is chosen so that temporal discretization errors are small relative to spatial errors. The adaptive procedure is applied to study singular current sheets in the tilt instability problem of ideal magnetohydrodynamics. Numerical results indicate a more accurate resolution of current sheets with higher-order methods than with piecewise-linear approximations.