Improved flux-surface parameterization through constrained nonlinear optimization

Parameterization of magnetic flux-surfaces is often used for magnetohydrodynamic stability analysis and microturbulence modeling in tokamaks. Shape parameters for such local parameterization of a (numerical) equilibrium are traditionally computed analytically using geometrically derived quantities. However, often the shape is approximated by the average of values for different sections of the flux-surface contour or a truncated series, which does not guarantee an optimal fit. Here, instead nonlinear least squares optimization is used to compute these parameters, with a weighted sum of squared error cost function that is robust to outliers. This method results in a lower total absolute error for both the parameterization of the flux-surface contour and the poloidal magnetic field density than current methods for several parameterizations based on the well-known "Miller geometry." Furthermore, rapid convergence of shape parameters is achieved, no approximate geometric measurements of the contour are needed, and the method is applicable to any analytical shape parameterization. Validation with local, linear gyrokinetic simulations using these optimized shape parameters showed reduced root mean square errors in both the growth rate and frequency spectra when compared with simulations based on numerical equilibria. In particular, the popular Turnbull–Miller parameterization benefits from this approach, extending its usability closer toward the last-closed flux-surface for cases with minor up-down asymmetry.