A Computational Method in Plasma Physics

In this book, we report on research in methods of computational magneto- hydrodynamics supported by the United States Department of Energy under Contract EY-76-C-02-3077 with New York University. The work has re- sulted in a computer code for mathematical analysis of the equilibrium and stability of a plasma in three dimensions with toroidal geometry but no sym- metry. The code is listed in the final chapter. Versions of it have been used for the design of experiments at the Los Alamos Scientific Laboratory and the Max Planck Institute for Plasma Physics in Garching. We are grateful to Daniel Barnes, Jeremiah Brackbill, Harold Grad, William Grossmann, Abraham Kadish, Peter Lax, Guthrie Miller, Arnulf Schliiter, and Harold Weitzner for many useful discussions of the theory. We are especially indebted to Franz Herrnegger for theoretical and pedagogical comments. Constance Engle has provided outstanding assistance with the typescript. We take pleasure in acknowledging the help of the staff of the Courant Mathematics and Com- puting Laboratory at New York University. In particular we should like to express our thanks to Max Goldstein, Kevin McAuliffe, Terry Moore, Toshi Nagano and Tsun Tam. Frances Bauer New York Octavio Betancourt September 1978 Paul Garabedian v Contents Chapter 1. Introduction 1 1. 1 Formulation of the Problem 1 1. 2 Discussion of Results 2 Chapter 2. The Variational Principle 4 4 2. 1 The Magnetostatic Equations 6 2. 2 Flux Constraints in the Plasma . 7 2. 3 The Ergodic Constraint .