Implement a Fusion Power Balance Formulary

[wdeshazer]

The plasma power balance in a basic form can be written as:

fusion heating + external heating >= radiation loss + transport loss. (Stacey, Weston "Fusion" 2010)

Relatively simple models for these various elements, many of which can be found in the aforementioned reference, can be collected and implemented in PlamaPy module resources so that plasma power balance studies could be done.

As an example, the fusion heating calculation is based on a fusion reactivity model that is fairly easy to implement (I have my own implementation in Mathcad that could be ported to python). The models exist for several fusion relevant reactions and can be extended for ones that I am not familiar with.

The radiation and transport losses have simple models as well that would not be too hard to implement.

External heating sources such as Neutral Beam Injection (NBI) Electron Cyclotron Resonance Heating (ECRH), and Ion Cyclotron Resonance Heating (ICRH) at a simple level just get modeled as a heat source. The community might have suggestions on how this part could be enriched.