13 Fixed boundary tokamak equilibrium problem

The fixed boundary equilibrium problem (also called the “inverse equilibrium problem” by some authors) refers to the case where the shape of a boundary magnetic surface is given and one is asked to solve the equilibrium within this magnetic surface. To make it convenient to deal with the shape of the boundary, one usually uses a general coordinates system which has one coordinate surface coinciding with the given magnetic surface. This makes it trivial to deal with the irregular boundary. To obtain the equilibrium, one needs to solve the GS equation in the general coordinate system.

  13.1 Toroidal elliptic operator in general coordinates
  13.2 Finite difference form of toroidal elliptic operator in general coordinate system
  13.3 Special treatment at coordinate origin, wrong! to be deleted
  13.4 Pressure and toroidal field function profile
  13.5 Boundary magnetic surface and initial coordinates
  13.6 Fixed boundary equilibrium numerical code
  13.7 Benchmark of the code
  13.8 Low-beta equilibrium vs. high-beta equilibrium
  13.9 Analytical form of Jacobian (need cleaning up)
  13.10 Grad-Shafranov equation with prescribed safety factor profile (to be finished)