MHD equilibrium is an important topic for fusion (and other MHD applications). A tokamak, in principle, is a toroidally symmetric fusion device and so MHD equilibrium can be reduced to solving the time independent MHD equations in axisymmetry. This produces the Grad-Shafranov equation (a two dimensional, nonlinear PDE) which has been solved using various techniques in the fusion community including finite difference, finite elements and spectral methods. A similar PDE exists if there is a plasma column with helical symmetry. Non-axisymmetric plasmas do occur in tokamaks as a result of instabilities and applied fields. However, if there is no symmetry angle there is no PDE to be solved. The current workhorse for finding non-axisymmetric equilibria uses energy minimization to find the equilibrium. New approaches to this problem that can use state of the art techniques are desirable. The speaker has formulated a coupled set of PDEs for the non-axisymmetric MHD equilibrium problem assuming that flux surfaces are nested (i.e. there are no magnetic islands) and has written this in weak form to use finite element method to solve the equations. The questions are around whether there is an optimal way to try to formulate the problem for FEM and to couple the equations, what sort of elements to use, if other solution techniques would be better suited and so on.
