Numerical solution of three-dimensional magnetic differential equations


Journal Article

A computer code is described that solves differential equations of the form B · {down triangle, open}f = 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-Schlüter currents for a stellarator. © 1988.

Full Text

Duke Authors

Cited Authors

  • Reiman, AH; Greenside, HS

Published Date

  • January 1, 1988

Published In

Volume / Issue

  • 75 / 2

Start / End Page

  • 423 - 443

Electronic International Standard Serial Number (EISSN)

  • 1090-2716

International Standard Serial Number (ISSN)

  • 0021-9991

Digital Object Identifier (DOI)

  • 10.1016/0021-9991(88)90121-0

Citation Source

  • Scopus