Calculates total cross section at a given deuterium energy, eb
,
for the Helium-3 branch of Deuterium-Deutrium nuclear reactions
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=Float64), | intent(in) | :: | eb | Deuterium energy [keV] |
Cross Section []
function d_d_fusion_he(eb) result(sigma)
!+Calculates total cross section at a given deuterium energy, `eb`,
!+for the Helium-3 branch of Deuterium-Deutrium nuclear reactions
!+
!+###Equation
!+$$ D + D \rightarrow He^3(0.82 MeV) + n(2.45 MeV) (50%)$$
!+
!+###References
!+* Equations 8-9
!+* Table IV in Ref. 7 [[atomic_tables(module)]]
real(Float64), intent(in) :: eb
!+ Deuterium energy [keV]
real(Float64) :: sigma
!+ Cross Section [\(cm^2\)]
real(Float64), dimension(5), parameter :: A = [ 5.3701d4, 3.3027d2, &
-1.2706d-1, 2.9327d-5, &
-2.5151d-9 ]
real(Float64), dimension(4), parameter :: B = [0.d0,0.d0,0.d0,0.d0]
real(Float64), parameter :: Bg = 31.3970
real(Float64) :: S, E
E = min(max(eb,0.5),4900.0)
S = (A(1) + E*(A(2) + E*(A(3) + E*(A(4) + E*A(5))))) / &
(1 + E*(B(1) + E*(B(2) + E*(B(3) + E*B(4)))))
sigma = (1.0d-27)*(S/(E*exp(Bg/sqrt(E))))
end function d_d_fusion_he