gafusion / omas

Ordered Multidimensional Array Structure
http://gafusion.github.io/omas
MIT License
32 stars 15 forks source link

Include fusion power in neutrons #179

Closed bclyons12 closed 2 years ago

bclyons12 commented 2 years ago

Strangely, summary.fusion.power is just the power coupled to the plasma (according the IMAS schema). This updates it, assuming it's 100% D-T fusion.

I'm not sure if we should worry about a contribution from D-D fusion. I think there's an O(1%) correction from that?

TimSlendebroek commented 2 years ago

@bclyons12

Thermal TT and DD might not be super significant but beam-target is a very significant fraction (for example the DT JET campaigns).

What I find weird is that they do split up all the neutron_fluxes but not the neutron powers:

While there is only (for the powers):

Maybe something to bring up to the IMAS dataschema people @bclyons12 @orso82 ?

If that is added we could do summary.fusion.neutron_power_total = the sum of all the individual components and fill in DT in TGRYO_step

bclyons12 commented 2 years ago

@TimSlendebroek Maybe because each fusion reaction results in monoenergetic neutrons, so if you have the neutron flux, you can compute the power?

TimSlendebroek commented 2 years ago

@bclyons12 I agree but that is true for almost every derived quantity (that is almost the definition of a derived quantity) and there are many of those in the IMAS dataschema.

Practically on a now time-frame what do you think about adding a neutron_power_sum() function that sums up all contributions and put the dt neutron flux in it's specific location. Something tells me that this isn't super obvious for beam target though (conservation of energy?).