skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A fully implicit numerical method for single-fluid resistive magnetohydrodynamics

Journal Article · · Journal of Computational Physics
 [1];  [2];  [3]
  1. Department of Mathematics, University of California at San Diego, 9500 Gilman Drive, Dept 0112, La Jolla, San Diego, CA 92093-0112 (United States)
  2. Princeton Plasma Physics Lab, MS 26, P.O. Box 451, Princeton, NJ 08543 (United States)
  3. Lawrence Livermore National Lab, P.O. Box 808, L-551, Livermore, CA 94551 (United States)

We present a nonlinearly implicit, conservative numerical method for integration of the single-fluid resistive MHD equations. The method uses a high-order spatial discretization that preserves the solenoidal property of the magnetic field. The fully coupled PDE system is solved implicitly in time, providing for increased interaction between physical processes as well as additional stability over explicit-time methods. A high-order adaptive time integration is employed, which in many cases enables time steps ranging from one to two orders of magnitude larger than those constrained by the explicit CFL condition. We apply the solution method to illustrative examples relevant to stiff magnetic fusion processes which challenge the efficiency of explicit methods. We provide computational evidence showing that for such problems the method is comparably accurate with explicit-time simulations, while providing a significant runtime improvement due to its increased temporal stability.

OSTI ID:
20840362
Journal Information:
Journal of Computational Physics, Vol. 219, Issue 1; Other Information: DOI: 10.1016/j.jcp.2006.03.022; PII: S0021-9991(06)00144-6; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
Country of Publication:
United States
Language:
English