bluemira.plasma_physics.reactions

Fusion reactions

Functions

`E_DT_fusion`(→ float)	Calculates the total energy released from the D-T fusion reaction
`E_DD_fusion`(→ float)	Calculates the total energy released from the D-D fusion reaction
`n_DT_reactions`(→ float)	Calculates the number of D-T fusion reactions per s for a given D-T fusion
`n_DD_reactions`(→ float)	Calculates the number of D-D fusion reactions per s for a given D-D fusion
`r_T_burn`(→ float)	Calculates the tritium burn rate for a given fusion power
`r_D_burn_DT`(→ float)	Calculates the deuterium burn rate for a given fusion power in D-T
`reactivity`(→ float \| numpy.ndarray)	Calculate the thermal reactivity of a fusion reaction in Maxwellian plasmas,

Module Contents

bluemira.plasma_physics.reactions.E_DT_fusion() → float

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

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

Notes

\[ \begin{align}\begin{aligned}{^{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\end{aligned}\end{align} \]

bluemira.plasma_physics.reactions.E_DD_fusion() → float

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

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

Notes

\[ \begin{align}\begin{aligned}{^{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\end{aligned}\end{align} \]

bluemira.plasma_physics.reactions.n_DT_reactions(p_fus: float) → float

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

\(n_{reactions} = \frac{P_{fus}[W]}{17.58 [MeV]eV[J]} [1/s]\)

Parameters:: p_fus (float) – D-T fusion power [W]
Return type:: Number of D-T reactions per second [1/s]

bluemira.plasma_physics.reactions.n_DD_reactions(p_fus: float) → float

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

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

Parameters:: p_fus (float) – D-D fusion power [W]
Return type:: Number of D-D reactions per second [1/s]

bluemira.plasma_physics.reactions.r_T_burn(p_fus: float) → float

Calculates the tritium burn rate for a given fusion power

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

Parameters:: p_fus (float) – D-T fusion power [W]
Return type:: T burn rate in the plasma [kg/s]

bluemira.plasma_physics.reactions.r_D_burn_DT(p_fus: float) → float

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

Parameters:: p_fus (float) – D-T fusion power [W]
Return type:: D burn rate in the plasma [g/s]

Notes

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

bluemira.plasma_physics.reactions.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>

Parameters:

temp_k (float | numpy.ndarray) – Temperature [K]
reaction (str | Reactions) – The fusion reaction
method (str | ReactivityMethod) – The parameterisation to use when calculating the reactivity

Returns:

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

Return type:

float | numpy.ndarray