ukaea / PROCESS

PROCESS is a systems code at UKAEA that calculates in a self-consistent manner the parameters of a fusion power plant with a specified performance, ensuring that its operating limits are not violated, and with the option to optimise to a given function of these parameters.
https://ukaea.github.io/PROCESS/
MIT License
34 stars 11 forks source link

Self-inductance of TF coil #850

Open jonmaddock opened 5 years ago

jonmaddock commented 5 years ago

In GitLab by @mkovari on May 29, 2019, 11:29

When calculating the inductance of the TF coil, the contribution from the cross-sectional area of the coil itself is calculated by taking the field as B(r)/2. Although this is correct for the field, the contribution to the inductance is not correct, because the flux passing through the coil itself is only linked by the turns running outside that point in the coil.

(In the supercon coil the winding pack is fairly thin and the flux passing through the winding pack is unlikely to be a major contribution.)

See #848 for somewhat similar calculations of force.

Checklist

After implementing issue do the following

jonmaddock commented 3 years ago

In GitLab by @jlion on Jan 12, 2021, 12:26

One possible way to calculate the inductance using elliptical integrals was implemented in the previous stellarator routine, e.g. see commit: d6d318dac444c916159c3dd36b3bac1beb1610f3

I attach the respective function. It uses some values from the outside scope, like total current or major radius which should be easy to retrieve.

function_coil_energy.f90

mkovari commented 2 days ago

Just a comment that, as far as I know, the self-inductance of the TF coil set is only used to derive the stored energy, which is itself just an output and is not used in any other calculation.