Categories:
Physics,
Plasma physics.
Lawson criterion
For sustained nuclear fusion to be possible,
the Lawson criterion must be met,
from which some required properties
of the plasma and the reactor chamber can be deduced.
Suppose that a reactor generates a given power Pfus by nuclear fusion,
but that it leaks energy at a rate Ploss in an unusable way.
If an auxiliary input power Paux sustains the fusion reaction,
then the following inequality must be satisfied
in order to have harvestable energy:
Ploss≤Pfus+Paux
We can rewrite Paux using the definition
of the energy gain factor Q,
which is the ratio of the output and input powers of the fusion reaction:
Q≡PauxPfus⟹Paux=QPfus
Returning to the inequality, we can thus rearrange its right-hand side as follows:
Ploss≤Pfus+QPfus=Pfus(1+Q1)=Pfus(QQ+1)
We assume that the plasma has equal species densities ni=ne,
so its total density n=2ni.
Then Pfus is as follows,
where fii is the frequency
with which a given ion collides with other ions,
and Efus is the energy released by a single fusion reaction:
Pfus=fiiniEfus=(ni⟨σv⟩)niEfus=4n2⟨σv⟩Efus
Where ⟨σv⟩ is the mean product
of the velocity v and the collision cross-section σ.
Furthermore, assuming that both species have the same temperature Ti=Te=T,
the total energy density W of the plasma is given by:
W=23kBTini+23kBTene=3kBTn
Where kB is Boltzmann’s constant.
From this, we can define the confinement time τE
as the characteristic lifetime of energy in the reactor, before leakage.
Therefore:
τE≡PlossW⟹Ploss=τE3nkBT
Inserting these new expressions for Pfus and Ploss
into the inequality, we arrive at:
τE3nkBT≤4n2⟨σv⟩Efus(QQ+1)
This can be rearranged to the form below,
which is the original Lawson criterion:
nτE≥Q+1Q⟨σv⟩Efus12kBT
However, it turns out that the highest fusion power density
is reached when T is at the minimum of T2/⟨σv⟩.
Therefore, we multiply by T to get the Lawson triple product:
nTτE≥Q+1Q⟨σv⟩Efus12kBT2
For some reason,
it is often assumed that the fusion is infinitely profitable Q→∞,
in which case the criterion reduces to:
nTτE≥⟨σv⟩Efus12kBT2
References
- M. Salewski, A.H. Nielsen,
Plasma physics: lecture notes,
2021, unpublished.