Initial NSF attempt

[nelsonand]

Hi @matt-pharr - I'm trying here to recreate your NSF POPCON on the new refactored branch...

Despite setting the profiles as indicated in the scripts you sent me on Slack, they come out incredibly hot:

Any idea why this might be? Could also be related to the hot profiles we saw with MANTA.

[matt-pharr]

@nelsonand I think it is due to a difference in averaging. I used volume-averaged temperature as "Average Temperature" when setting up the POPCONs -- since outer flux surfaces have much more volume than inner, this will dramatically lower the volume-averaged temperature compared to the line-averaged temperature, which weights interior and exterior flux surfaces equally.

[nelsonand]

@matt-pharr got it! Thanks! ...but I'm not sure how to change the script. I can drop the tmax here:

OpenPOPCON/resources/examples/NSF/NSF_POPCON.yaml

Line 59 in 59e797d

Tmax_keV: 5. # maximum volume-averaged temperature in keV

but then it struggles to converge

[matt-pharr]

@nelsonand I finally fixed the issue in #17 that was causing non convergence at low-temperatures. I also implemented D-D fuel cycles. This is what I have for a basic example now:

I have notebooks going for MANTA and NSF and will add those to the PR shortly. I will use your setup here to replace my NSF example, once we have an agreement that we have reached something physically satisfactory!

The grayed-out area is where impurity radiation exceeds W/tau_E (and therefore the operating point is not achievable with the requested level of impurity). I used an impfrac of 1e-8 n_i tungsten.