d_d_fusion_he Function

  1. atomic_tables.f90
  2. atomic_tables
  3. d_d_fusion_he

public 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

References

Arguments

TypeIntentOptionalAttributesName
real(kind=Float64), intent(in) :: eb

Deuterium energy [keV]

Return Value real(kind=Float64)

Cross Section []


Called by

proc~~d_d_fusion_he~~CalledByGraph proc~d_d_fusion_he d_d_fusion_he proc~write_bb_d_d write_bb_D_D proc~write_bb_d_d->proc~d_d_fusion_he program~generate_tables generate_tables program~generate_tables->proc~write_bb_d_d

Contents

Source Code

d_d_fusion_he

Source Code

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