The beam fusion model is used to calculation fusion power from beam/plasma interactions.
beam_fusion() -> beamcalc() -> alpha_power_beam() -> bosch_hale_reactivity()
bosch_hale_reactivity should be passed a profile, currently it is passed a volume averaged temperature.
This results in a much higher value of beam_fusion being calculated than might be expected and is inconsistent with the way the function should work according to its own doc-string and reference..
This is also inconsistent with how the calculation is done elsewhere in the code, where fusion_rate_integral() is used.
This affects a correction ratio for the reactivity in alpha_power_beam() and has a significant impact on the beam fusion power calculated by PROCESS.
To amend this, we should pass in an ion temperature profile directly; and integrate this as is done elsewhere in the code, and use this integral to calculate the correction ratio.
The beam fusion model is used to calculation fusion power from beam/plasma interactions.
beam_fusion()->beamcalc()->alpha_power_beam()->bosch_hale_reactivity()bosch_hale_reactivityshould be passed a profile, currently it is passed a volume averaged temperature.This results in a much higher value of beam_fusion being calculated than might be expected and is inconsistent with the way the function should work according to its own doc-string and reference..
This is also inconsistent with how the calculation is done elsewhere in the code, where
fusion_rate_integral()is used.This affects a correction ratio for the reactivity in
alpha_power_beam()and has a significant impact on the beam fusion power calculated by PROCESS.To amend this, we should pass in an ion temperature profile directly; and integrate this as is done elsewhere in the code, and use this integral to calculate the correction ratio.