alfkoehn/fusion_plots

Plots and figures relevant for a fusion research lecture

fusion_plots's Issues

Current sort isn't ideal

The table is currently sorted by date. I think it would be more intuitive to sort the table by company name, since you're much more likely to know the name of a company you want info about, rather than the date it was founded.

Ideally you could make each column sortable on the header, but it doesn't look like markdown supports that.

Evolution of triple product record vs year

Hi Alf,

I've seen a number of presentations over the years which show the triple-product record vs time.
These are then often compared to, e.g. Moore's law.

Do any of the triple-product datasets contain year information, or would you know where to find the data for this?

If I know where to get the data, I would volunteer a PR with the figure.

Cheers,
Thomas

pp reaction cross section

Hi! thanks for that great package. Very useful!

I'd like to plot cross sections with the pp reaction. Do you think there's a way of adding this to fusion_cross_sections.py ?

I've come up with that script so far but must've messed up something with the units...

import numpy as np
import matplotlib.pyplot as plt
from scipy.constants import e, atomic_mass, hbar


ev_to_joule = 1.60218e-19  # J/eV
barn_to_m2 = 1e-28  # m2/b

E = np.logspace(0, 6)  # eV


A1 = A2 = 1.0072765
mu = atomic_mass*A1*A2/(A1+A2)
vel = (2*E*ev_to_joule/mu)**0.5

eta = e**2/hbar/vel
S_0 = 4e-25*1e6  # eV b
S_prime_0 = 11.2  # b
S = S_0 + S_prime_0 * E
cross_section = S/E*np.exp(-2*np.pi*eta)
cross_section *= barn_to_m2

plt.plot(E, cross_section)
plt.xscale("log")
plt.yscale("log")
plt.show()

Referencing

Hi @alfkoehn and @jhillairet

I would like to use these lovely plots in my thesis, just wondering if you would like a reference to this repository?

Energy vs temperature in cross section

fusion_plots/fusion_cross_sections/fusion_cross_sections.py

Line 179 in aace104

# (actually energy, but in plasma physics we just call it temperature :-)

On the plot of cross section vs T, I always understood that the plots which show energy are "per-particle", whereas the plots which show temperature are integrated over a Maxwell-Boltzmann distribution -- this also leads to a reduction in the peak T/E for the DT distribution thanks to the Maxwell tail.

The comment on L179 of the script seems to suggest that T vs E is just taking temperature vs energy, presumably just considering E = E_thermal for a distribution?

Anyway, the long and short of my question is: Is it just sigma(energy), scaled by 1/1.1e4 to translate between energy (eV) and temperature (K), or are thermal distributions taken into account?