MatthewGrim
commented
5 years ago The Bennett equilibrium can be calculated easily and the global beta is easy to calculate. I noticed a numerical issue close to the axis where my local beta will blow up, but I have not investigated this further as the global beta is what is considered significant.
Now I can:
- Implement the same functions for the slightly more complex screw pinch
- Implement the other plasma parameters: kink safety factor and q.
After that I can experiment with the optimised profiles, but I am not sure if this will be useful as none of this analysis takes into account stability constraints.
Based off the analysis in Hartmut Zohm and Friedberg's books, this issue documents the implementation of a Bennett equilibrium solver. This is a basic integration to solve the pressure balance, but I am using it to get an better feel for plasma parameters and different equilibrium types.