Reaction rate calculation

[kiselevnightkeo]

Can you explain how you calculate the reaction rate in the "bt_maxwellian_eb" function?

I thought it was calculating the beam-target reaction rate, but using the formula in the attached picture did not give the desired result.

Cross section from atomic_tables.h5 (!!)

I used the following article. fusion_cross_sections_and_reactivities_by_G_Miley_H_Towner_N_Ivich.pdf

[lstagner]

What result did your calculation give?

[kiselevnightkeo]

beam energy eb = 100 keV Ion temperature Ti = 10 keV

D-D reaction, neutron branch.

Reaction rate by fidasim <sigmav>= 6.95e-22 m^3/s Reaction rate by formulae <sigmav>= 1.39e-23 m^3/s

Could this be related to the 1/2 coefficient from Kronecker’s delta?

But for other lower energy and temperature difference is bigger.

beam energy eb = 10 keV Ion temperature Ti = 1 keV Reaction rate by fidasim <sigmav>= 7.97e-27 m^3/s Reaction rate by formulae <sigmav>= 9.4e-26 m^3/s

[lstagner]

So a factor of 50. Below is how we calculate

reaction_rates.pdf