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

Title: Solving the 3D MHD equilibrium equations in toroidal geometry by Newton's method

Abstract

We describe a novel form of Newton's method for computing 3D MHD equilibria. The method has been implemented as an extension to the hybrid spectral/finite-difference Princeton Iterative Equilibrium Solver (PIES) which normally uses Picard iteration on the full nonlinear MHD equilibrium equations. Computing the Newton functional derivative numerically is not feasible in a code of this type but we are able to do the calculation analytically in magnetic coordinates by considering the response of the plasma's Pfirsch-Schlueter currents to small changes in the magnetic field. Results demonstrate a significant advantage over Picard iteration in many cases, including simple finite-{beta} stellarator equilibria. The method shows promise in cases that are difficult for Picard iteration, although it is sensitive to resolution and imperfections in the magnetic coordinates, and further work is required to adapt it to the presence of magnetic islands and stochastic regions.

Authors:
 [1];  [2];  [3]
  1. Princeton Plasma Physics Laboratory, Princeton University, NJ (United States). E-mail: h.oliver@niwa.co.nz
  2. Princeton Plasma Physics Laboratory, Princeton University, NJ (United States). E-mail: reiman@pppl.gov
  3. Princeton Plasma Physics Laboratory, Princeton University, NJ (United States). E-mail: monticello@pppl.gov
Publication Date:
OSTI Identifier:
20687269
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 211; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2005.05.007; PII: S0021-9991(05)00264-0; Copyright (c) 2005 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; COORDINATES; EQUATIONS; GEOMETRY; MAGNETIC FIELDS; MAGNETIC ISLANDS; MAGNETOHYDRODYNAMICS; MHD EQUILIBRIUM; NEWTON METHOD; NONLINEAR PROBLEMS; STOCHASTIC PROCESSES; THREE-DIMENSIONAL CALCULATIONS

Citation Formats

Oliver, H.J., Reiman, A.H., and Monticello, D.A. Solving the 3D MHD equilibrium equations in toroidal geometry by Newton's method. United States: N. p., 2006. Web. doi:10.1016/j.jcp.2005.05.007.
Oliver, H.J., Reiman, A.H., & Monticello, D.A. Solving the 3D MHD equilibrium equations in toroidal geometry by Newton's method. United States. doi:10.1016/j.jcp.2005.05.007.
Oliver, H.J., Reiman, A.H., and Monticello, D.A. Sun . "Solving the 3D MHD equilibrium equations in toroidal geometry by Newton's method". United States. doi:10.1016/j.jcp.2005.05.007.
@article{osti_20687269,
title = {Solving the 3D MHD equilibrium equations in toroidal geometry by Newton's method},
author = {Oliver, H.J. and Reiman, A.H. and Monticello, D.A.},
abstractNote = {We describe a novel form of Newton's method for computing 3D MHD equilibria. The method has been implemented as an extension to the hybrid spectral/finite-difference Princeton Iterative Equilibrium Solver (PIES) which normally uses Picard iteration on the full nonlinear MHD equilibrium equations. Computing the Newton functional derivative numerically is not feasible in a code of this type but we are able to do the calculation analytically in magnetic coordinates by considering the response of the plasma's Pfirsch-Schlueter currents to small changes in the magnetic field. Results demonstrate a significant advantage over Picard iteration in many cases, including simple finite-{beta} stellarator equilibria. The method shows promise in cases that are difficult for Picard iteration, although it is sensitive to resolution and imperfections in the magnetic coordinates, and further work is required to adapt it to the presence of magnetic islands and stochastic regions.},
doi = {10.1016/j.jcp.2005.05.007},
journal = {Journal of Computational Physics},
number = 1,
volume = 211,
place = {United States},
year = {Sun Jan 01 00:00:00 EST 2006},
month = {Sun Jan 01 00:00:00 EST 2006}
}
  • The Schur-decomposition for three-dimensional matrix equations is developed and used to directly solve the radiative discrete ordinates equations which are discretized by Chebyshev collocation spectral method. Three methods, say, the spectral methods based on 2D and 3D matrix equation solvers individually, and the standard discrete ordinates method, are presented. The numerical results show the good accuracy of spectral method based on direct solvers. The CPU time cost comparisons against the resolutions between these three methods are made using MATLAB and FORTRAN 95 computer languages separately. The results show that the CPU time cost of Chebyshev collocation spectral method with 3Dmore » Schur-decomposition solver is the least, and almost only one thirtieth to one fiftieth CPU time is needed when using the spectral method with 3D Schur-decomposition solver compared with the standard discrete ordinates method.« less
  • A computer code has been constructed that solves the full magnetohydrodynamic (MHD) equations in toroidal geometry. The code is applicable to toroidal devices, including tokamaks, stellarators, and reversed field pinches. A fully implicit numerical technique is used that allows linear eigenvalues and eigenfunctions to be found in a very few computational steps. Although the present work describes the solution of the linearized equations, generalization of the numerical method to the solution of the nonlinear problem is straightforward. Use of the code is illustrated by calculating the n = 1 instability for a tokamak configuration. The results show the structural changesmore » in the eigenfunctions as the plasma pressure is increased.« less
  • A solution to the Euler-Ostrogradskii equation is found for the toroidal current.(AIP)
  • Stable and efficient implicit and explicit fractional-step methods for solving three-dimensional, time-dependent Maxwell equations have been successfully developed. These numerical procedures are characteristic-based schemes with the intrinsically accurate noreflection wave condition on the boundaries of truncated computational domain. Excellent simulations for electromagnetic phenomena have been achieved for a three-dimensional wave guide and an oscillating electric dipole.
  • The few-group diffusion codes programmed for various machines are inadequate for classes of problems requiring higher approximations to the transport equation. It is possible, however, to do one-velocity P/sub 3/ or double-P/sub 1/ problems, in slab geometry, with the aid of a slightly modified diffusion code. This diffusion code is discussed. (A.C.)