Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Numerical solution of three-dimensional magnetic differential equations

Journal Article · · J. Comput. Phys.; (United States)
A computer code is described that solves differential equations of the form B x delf = h for a single-valued solution f, given a toroidal three-dimensional divergence-free field B and a single-valued function h. The code uses a new algorithm that Fourier decomposes a given function in a set of flux coordinates in which the field lines are straight. The algorithm automatically adjusts the required integration lengths to compensate for proximity to low order rational surfaces. Applying this algorithm to the Cartesian coordinates defines a transformation to magnetic coordinates, in which the magnetic differential equation can be accurately solved. Our method is illustrated by calculating the Pfirsch--Schlueter currents for a stellarator. copyright 1988 Academic Press, Inc.
Research Organization:
Plasma Physics Laboratory, Princeton University, Princeton, New Jersey 08544
OSTI ID:
5267605
Journal Information:
J. Comput. Phys.; (United States), Journal Name: J. Comput. Phys.; (United States) Vol. 75:2; ISSN JCTPA
Country of Publication:
United States
Language:
English