Nonlinear magnetohydrodynamics simulation using high-order finite elements

  • Authors:
  • C. R. Sovinec;A. H. Glasser;T. A. Gianakon;D. C. Barnes;R. A. Nebel;S. E. Kruger;D. D. Schnack;S. J. Plimpton;A. Tarditi;M. S. Chu

  • Affiliations:
  • Department of Engineering Physics, University of Wisconsin-Madison, 1500 Engineering Drive, Madison, WI;Los Alamos National Laboratory, Los Alamos, NM;Los Alamos National Laboratory, Los Alamos, NM;Los Alamos National Laboratory, Los Alamos, NM;Los Alamos National Laboratory, Los Alamos, NM;Science Applications International Corporation, San Diego, CA;Science Applications International Corporation, San Diego, CA;Sandia National Laboratory, Albuquerque, NM;Advanced Space Propulsion Laboratory, National Aeronautics and Space Administration-Johnson Space Center, Houston, TX;General Atomics Corporation, San Diego, CA

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

Quantified Score

Hi-index 31.51

Visualization

Abstract

A conforming representation composed of 2D finite elements and finite Fourier series is applied to 3D nonlinear non-ideal magnetohydrodynamics using a semi-implicit time-advance. The self-adjoint semi-implicit operator and variational approach to spatial discretization are synergistic and enable simulation in the extremely stiff conditions found in high temperature plasmas without sacrificing the geometric flexibility needed for modeling laboratory experiments. Growth rates for resistive tearing modes with experimentally relevant Lundquist number are computed accurately with time-steps that are large with respect to the global Alfvén time and moderate spatial resolution when the finite elements have basis functions of polynomial degree (p) two or larger. An error diffusion method controls the generation of magnetic divergence error. Convergence studies show that this approach is effective for continuous basis functions with p ≥ 2, where the number of test functions for the divergence control terms is less than the number of degrees of freedom in the expansion for vector fields. Anisotropic thermal conduction at realistic ratios of parallel to perpendicular conductivity (χ||/χ⊥) is computed accurately with p ≥ 3 without mesh alignment. A simulation of tearingmode evolution for a shaped toroidal tokamak equilibrium demonstrates the effectiveness of the algorithm in nonlinear conditions, and its results are used to verify the accuracy of the numerical anisotropic thermal conduction in 3D magnetic topologies.