bluemira.plasma_physics.reactions
=================================

.. py:module:: bluemira.plasma_physics.reactions

.. autoapi-nested-parse::

   Fusion reactions



Functions
---------

.. autoapisummary::

   bluemira.plasma_physics.reactions.E_DT_fusion
   bluemira.plasma_physics.reactions.E_DD_fusion
   bluemira.plasma_physics.reactions.n_DT_reactions
   bluemira.plasma_physics.reactions.n_DD_reactions
   bluemira.plasma_physics.reactions.r_T_burn
   bluemira.plasma_physics.reactions.r_D_burn_DT
   bluemira.plasma_physics.reactions.reactivity


Module Contents
---------------

.. py:function:: E_DT_fusion() -> float

   Calculates the total energy released from the D-T fusion reaction

   :rtype: The energy released from a single D-T fusion reaction [J]

   .. rubric:: Notes

   .. math::
       {^{2}_{1}H}+{^{3}_{1}H}~\rightarrow~{^{4}_{2}He}~
       (3.5~\text{MeV})+\text{n}^{0} (14.1 ~\text{MeV})

       \Delta E = \Delta m c^2


.. py:function:: E_DD_fusion() -> float

   Calculates the total energy released from the D-D fusion reaction

   :rtype: The energy released from a single D-D fusion reaction [J]

   .. rubric:: Notes

   .. math::
       {^{2}_{1}H}+{^{2}_{1}H}~\rightarrow~{^{3}_{1}H}
       (1.01 ~\text{MeV})+\text{p} (3.02~\text{MeV})~~[50 \textrm{\%}]
       ~~~~~~~~~~\rightarrow~{^{3}_{2}He} (0.82~\text{MeV})+\text{n}^{0} (2.45~\text{MeV})~~[50 \text{\%}]

       \Delta E = \Delta m c^2


.. py:function:: n_DT_reactions(p_fus: float) -> float

   Calculates the number of D-T fusion reactions per s for a given D-T fusion
   power

   :math:`n_{reactions} = \frac{P_{fus}[W]}{17.58 [MeV]eV[J]} [1/s]`

   :param p_fus: D-T fusion power [W]

   :rtype: Number of D-T reactions per second [1/s]


.. py:function:: n_DD_reactions(p_fus: float) -> float

   Calculates the number of D-D fusion reactions per s for a given D-D fusion
   power

   :math:`n_{reactions} = \frac{P_{fus}[W]}{E_{DD} [MeV] eV[J]} [1/s]`

   :param p_fus: D-D fusion power [W]

   :rtype: Number of D-D reactions per second [1/s]


.. py:function:: r_T_burn(p_fus: float) -> float

   Calculates the tritium burn rate for a given fusion power

   :math:`\dot{m_{b}} = \frac{P_{fus}[MW]M_{T}[g/mol]}{17.58 [MeV]eV[J]N_{A}[1/mol]} [g/s]`

   :param p_fus: D-T fusion power [W]

   :rtype: T burn rate in the plasma [kg/s]


.. py:function:: r_D_burn_DT(p_fus: float) -> float

   Calculates the deuterium burn rate for a given fusion power in D-T

   :param p_fus: D-T fusion power [W]

   :rtype: D burn rate in the plasma [g/s]

   .. rubric:: Notes

   .. math::
       \dot{m_{b}} = \frac{P_{fus}[MW]M_{D}[g/mol]}
       {17.58 [MeV]eV[J]N_{A}[1/mol]} [kg/s]


.. py:function:: reactivity(temp_k: float | numpy.ndarray, reaction: str | Reactions = Reactions.D_T, method: str | ReactivityMethod = ReactivityMethod.BOSCH_HALE) -> float | numpy.ndarray

   Calculate the thermal reactivity of a fusion reaction in Maxwellian plasmas,
   \t:math:`<\sigma v>`

   :param temp_k: Temperature [K]
   :param reaction: The fusion reaction
   :param method: The parameterisation to use when calculating the reactivity

   :returns: Reactivity of the reaction at the specified temperature(s) [m^3/s]