Free-Boundary MHD Equilibrium Solver for Axisymmetric Tokamak
This code solves Grad-Shafranov Equation for axisymmetric Tokamak plasma with given poloidal field coil and limiter configuration. The main algorithm is based on this article.
Currently, a basic version of the code is available here. Other versions like plasma position constrained, betap constrained and inverse MHD Equilibrium solver will also be uploaded here soon.
To use the code, open the ipynb file using the Jupyter Notebook/Lab. The required Python-3 libraries are numpy, scipy, matplotlib and scikit-image. I expect the code to work with the latest version of these libraries, so if you encounter an issue please let me know.
If you used this code for your research, please cite this repository and the following article: https://doi.org/10.1063/5.0290301. This article is also available on arXiv.
Critical updates coming soon:
- A more robust algorithm to find psi surfaces, the current one might fail in some rare cases.
- A faster computation of Green's functions part using Lackner and Von Hagenow formulation.
If output G-EQDSK file does not work with GENRAY code, reduce the length of the header line to 61 characters by removing any non-essential characters (like in-between spaces).
For any feedback, query or discussion, please contact me on email or Telegram. Developement of this code has been supported by Institute for Plasma Research (India).