FACIT -- Generalize ICRH effects

[cjperks7]

I can do this one too

Issue

The formula used in FACIT for the ICRH effects on poloidal asymmetry assumes just D(H) minority heating. Future fusion devices (SPARC) will use other schemes such as D(He3), D/T(He3).

Plan of action

In the original derivation by M. Reinke, they obtained a general expression and then specified it to D(H) which was grabbed for use it FACIT. Pretty straight forward to change that to account for other schemes.

[fsciortino]

Of interest for @DanielFajardoJ

[DanielFajardoJ]

This would be a good improvement for ICRH applications, thanks @cjperks7

[odstrcilt]

More accurate formulas for poloidal asymmetry by ICRH are here:
https://doi.org/10.1088/1741-4326/aa5fd6
but still, you will need Tper/Tpar on LFS from TORIC.

[cjperks7]

Oops, seems like I forgot to do this... ICRF folks are (finally) starting to make a database of Tperp/Tpar for SPARC and the broadband SXR diagnostic is coming up on Final Design Review, so it would be of interest soon.

Thanks for the reference and the reminder @odstrcilt !

[cjperks7]

@odstrcilt I was looking through this new paper and see you were second author. I'm curious what in your opinion makes it more accurate v. say Matt Reinke's C-Mod work? I particularly noticed this paper seems to be all theory, so have you tried anything comparison with experiment, i.e. DIII-D? Maybe I could get Matt's C-Mod shots if you want to look at that to quantify improvement.

[odstrcilt]

It is, in principle, the same idea, the similar formula, just more carefully derived with less approximations.

Most important is the eq. 15 which provides you the poloidal distribution of the high-Z impurity

it depends on a Nm/NmLFS - which is the poloidal profile of density of fast ions.

I have compared it with experiment on DIII-D.
https://iopscience.iop.org/article/10.1088/1361-6587/aa8690
the same formula for asymmetry is derived there at well, just I went directly to the approximation.

You can either use fast ion profile from NUBEAM for NBI or calculate fast ion profile from toric.
The second part depends on formulas in the previous part of the paper.
Depending on the version of toric, it calculates the LFS Tpero and Tpar or flux surface averaged p_perp and p_par. The toric in TRANSP does the second.

Here is an example how I did it for NEO:
https://github.com/gafusion/OMFIT-source/blob/auto_merge/modules/NEO_GACODE/SCRIPTS/add_fast_ions.py
This script was working with some older version of NEO module, it will require some fixes to make it working again.

The asymmetry by beam ions is not included. If I'm lucky, I can find the scripts somewhere...