Lawson criterion

Last updated
Lawson criterion of important magnetic confinement fusion experiments Fustion triple-product diagram Horvath.webp
Lawson criterion of important magnetic confinement fusion experiments

The Lawson criterion is a figure of merit used in nuclear fusion research. It compares the rate of energy being generated by fusion reactions within the fusion fuel to the rate of energy losses to the environment. When the rate of production is higher than the rate of loss, the system will produce net energy. If enough of that energy is captured by the fuel, the system will become self-sustaining and is said to be ignited.

Contents

The concept was first developed by John D. Lawson in a classified 1955 paper [1] that was declassified and published in 1957. [2] As originally formulated, the Lawson criterion gives a minimum required value for the product of the plasma (electron) density ne and the "energy confinement time" that leads to net energy output.

Later analysis suggested that a more useful figure of merit is the triple product of density, confinement time, and plasma temperature T. The triple product also has a minimum required value, and the name "Lawson criterion" may refer to this value.

On August 8, 2021, researchers at Lawrence Livermore National Laboratory's National Ignition Facility in California confirmed to have produced the first-ever successful ignition of a nuclear fusion reaction surpassing the Lawson's criteria in the experiment. [3] [4]

Energy balance

The central concept of the Lawson criterion is an examination of the energy balance for any fusion power plant using a hot plasma. This is shown below:

Net power = Efficiency × (Fusion − Radiation loss − Conduction loss)

  1. Net power is the excess power beyond that needed internally for the process to proceed in any fusion power plant.
  2. Efficiency is how much energy is needed to drive the device and how well it collects energy from the reactions.
  3. Fusion is rate of energy generated by the fusion reactions.
  4. Radiation loss is the energy lost as light (including X-rays) leaving the plasma.
  5. Conduction loss is the energy lost as particles leave the plasma, carrying away energy.

Lawson calculated the fusion rate by assuming that the fusion reactor contains a hot plasma cloud which has a Gaussian curve of individual particle energies, a Maxwell–Boltzmann distribution characterized by the plasma's temperature. Based on that assumption, he estimated the first term, the fusion energy being produced, using the volumetric fusion equation. [5]

Fusion = Number density of fuel A × Number density of fuel B × Cross section(Temperature) × Energy per reaction

  1. Fusion is the rate of fusion energy produced by the plasma
  2. Number density is the density in particles per unit volume of the respective fuels (or just one fuel, in some cases)
  3. Cross section is a measure of the probability of a fusion event, which is based on the plasma temperature
  4. Energy per reaction is the energy released in each fusion reaction

This equation is typically averaged over a population of ions which has a normal distribution. The result is the amount of energy being created by the plasma at any instant in time.

Lawson then estimated [5] the radiation losses using the following equation:

where N is the number density of the cloud and T is the temperature. For his analysis, Lawson ignores conduction losses. In reality this is nearly impossible; practically all systems lose energy through mass leaving the plasma and carrying away its energy.

By equating radiation losses and the volumetric fusion rates, Lawson estimated the minimum temperature for the fusion for the deuterium tritium (D-T) reaction

to be 30 million degrees (2.6 keV), and for the deuteriumdeuterium (D-D) reaction

to be 150 million degrees (12.9 keV). [2] [6]

Extensions into nτE

The confinement time measures the rate at which a system loses energy to its environment. The faster the rate of loss of energy, , the shorter the energy confinement time. It is the energy density (energy content per unit volume) divided by the power loss density (rate of energy loss per unit volume):

For a fusion reactor to operate in steady state, the fusion plasma must be maintained at a constant temperature. Thermal energy must therefore be added at the same rate the plasma loses energy in order to maintain the fusion conditions. This energy can be supplied by the fusion reactions themselves, depending on the reaction type, or by supplying additional heating through a variety of methods.

For illustration, the Lawson criterion for the D-T reaction will be derived here, but the same principle can be applied to other fusion fuels. It will also be assumed that all species have the same temperature, that there are no ions present other than fuel ions (no impurities and no helium ash), and that D and T are present in the optimal 50-50 mixture. a Ion density then equals electron density and the energy density of both electrons and ions together is given, according to the ideal gas law, by

where is the temperature in electronvolt (eV) and is the particle density.

The volume rate (reactions per volume per time) of fusion reactions is

where is the fusion cross section, is the relative velocity, and denotes an average over the Maxwellian velocity distribution at the temperature .

The volume rate of heating by fusion is times , the energy of the charged fusion products (the neutrons cannot help to heat the plasma). In the case of the D-T reaction, .

The Lawson criterion, or minimum value of (electron density * energy confinement time) required for self-heating, for three fusion reactions. For DT, ntE minimizes near the temperature 25 keV (300 million kelvins). Fusion ntau.svg
The Lawson criterion, or minimum value of (electron density * energy confinement time) required for self-heating, for three fusion reactions. For DT, nτE minimizes near the temperature 25 keV (300 million kelvins).

The Lawson criterion requires that fusion heating exceeds the losses:

Substituting in known quantities yields:

Rearranging the equation produces:

The quantity is a function of temperature with an absolute minimum. Replacing the function with its minimum value provides an absolute lower limit for the product . This is the Lawson criterion.

For the deuterium tritium reaction, the physical value is at least

The minimum of the product occurs near .

Extension into the "triple product"

A still more useful figure of merit is the "triple product" of density, temperature, and confinement time, nTτE. For most confinement concepts, whether inertial, mirror, or toroidal confinement, the density and temperature can be varied over a fairly wide range, but the maximum attainable pressure p is a constant. When such is the case, the fusion power density is proportional to p2v>/T 2. The maximum fusion power available from a given machine is therefore reached at the temperature T where <σv>/T 2 is a maximum. By continuation of the above derivation, the following inequality is readily obtained:

The fusion triple product condition for three fusion reactions Fusion tripleprod.svg
The fusion triple product condition for three fusion reactions

The quantity is also a function of temperature with an absolute minimum at a slightly lower temperature than .

For the D-T reaction, the minimum occurs at T = 14 keV. The average <σv> in this temperature region can be approximated as [7]

so the minimum value of the triple product value at T = 14 keV is about

This number has not yet been achieved in any reactor, although the latest generations of machines have come close. JT-60 reported 1.53x1021 keV.s.m−3. [8] For instance, the TFTR has achieved the densities and energy lifetimes needed to achieve Lawson at the temperatures it can create, but it cannot create those temperatures at the same time. ITER aims to do both.

As for tokamaks, there is a special motivation for using the triple product. Empirically, the energy confinement time τE is found to be nearly proportional to n1/3/P 2/3[ citation needed ]. In an ignited plasma near the optimum temperature, the heating power P equals fusion power and therefore is proportional to n2T 2. The triple product scales as

The triple product is only weakly dependent on temperature as T -1/3. This makes the triple product an adequate measure of the efficiency of the confinement scheme.

Inertial confinement

The Lawson criterion applies to inertial confinement fusion (ICF) as well as to magnetic confinement fusion (MCF) but in the inertial case it is more usefully expressed in a different form. A good approximation for the inertial confinement time is the time that it takes an ion to travel over a distance R at its thermal speed

where mi denotes mean ionic mass. The inertial confinement time can thus be approximated as

By substitution of the above expression into relationship ( 1 ), we obtain

This product must be greater than a value related to the minimum of T 3/2/<σv>. The same requirement is traditionally expressed in terms of mass density ρ = <nmi>:

Satisfaction of this criterion at the density of solid D-T (0.2 g/cm3) would require a laser pulse of implausibly large energy. Assuming the energy required scales with the mass of the fusion plasma (Elaser ~ ρR3 ~ ρ−2), compressing the fuel to 103 or 104 times solid density would reduce the energy required by a factor of 106 or 108, bringing it into a realistic range. With a compression by 103, the compressed density will be 200 g/cm3, and the compressed radius can be as small as 0.05 mm. The radius of the fuel before compression would be 0.5 mm. The initial pellet will be perhaps twice as large since most of the mass will be ablated during the compression.

The fusion power times density is a good figure of merit to determine the optimum temperature for magnetic confinement, but for inertial confinement the fractional burn-up of the fuel is probably more useful. The burn-up should be proportional to the specific reaction rate (n2<σv>) times the confinement time (which scales as T -1/2) divided by the particle density n:

Thus the optimum temperature for inertial confinement fusion maximises <σv>/T3/2, which is slightly higher than the optimum temperature for magnetic confinement.

Non-thermal systems

Lawson's analysis is based on the rate of fusion and loss of energy in a thermalized plasma. There is a class of fusion machines that do not use thermalized plasmas but instead directly accelerate individual ions to the required energies. The best-known examples are the migma, fusor and polywell.

When applied to the fusor, Lawson's analysis is used as an argument that conduction and radiation losses are the key impediments to reaching net power. Fusors use a voltage drop to accelerate and collide ions, resulting in fusion. [9] The voltage drop is generated by wire cages, and these cages conduct away particles.

Polywells are improvements on this design, designed to reduce conduction losses by removing the wire cages which cause them. [10] Regardless, it is argued that radiation is still a major impediment. [11]

See also

Notes

^a It is straightforward to relax these assumptions. The most difficult question is how to define when the ion and electrons differ in density and temperature. Considering that this is a calculation of energy production and loss by ions, and that any plasma confinement concept must contain the pressure forces of the plasma, it seems appropriate to define the effective (electron) density through the (total) pressure as . The factor of is included because usually refers to the density of the electrons alone, but here refers to the total pressure. Given two species with ion densities , atomic numbers , ion temperature , and electron temperature , it is easy to show that the fusion power is maximized by a fuel mix given by . The values for , , and the power density must be multiplied by the factor . For example, with protons and boron () as fuel, another factor of must be included in the formulas. On the other hand, for cold electrons, the formulas must all be divided by (with no additional factor for ).

Related Research Articles

<span class="mw-page-title-main">Uncertainty principle</span> Foundational principle in quantum physics

The uncertainty principle, also known as Heisenberg's indeterminacy principle, is a fundamental concept in quantum mechanics. It states that there is a limit to the precision with which certain pairs of physical properties, such as position and momentum, can be simultaneously known. In other words, the more accurately one property is measured, the less accurately the other property can be known.

In statistical mechanics, the virial theorem provides a general equation that relates the average over time of the total kinetic energy of a stable system of discrete particles, bound by a conservative force with that of the total potential energy of the system. Mathematically, the theorem states that where T is the total kinetic energy of the N particles, Fk represents the force on the kth particle, which is located at position rk, and angle brackets represent the average over time of the enclosed quantity. The word virial for the right-hand side of the equation derives from vis, the Latin word for "force" or "energy", and was given its technical definition by Rudolf Clausius in 1870.

In physics, a Langevin equation is a stochastic differential equation describing how a system evolves when subjected to a combination of deterministic and fluctuating ("random") forces. The dependent variables in a Langevin equation typically are collective (macroscopic) variables changing only slowly in comparison to the other (microscopic) variables of the system. The fast (microscopic) variables are responsible for the stochastic nature of the Langevin equation. One application is to Brownian motion, which models the fluctuating motion of a small particle in a fluid.

<span class="mw-page-title-main">Fokker–Planck equation</span> Partial differential equation

In statistical mechanics and information theory, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. The Fokker-Planck equation has multiple applications in information theory, graph theory, data science, finance, economics etc.

The Ising model, named after the physicists Ernst Ising and Wilhelm Lenz, is a mathematical model of ferromagnetism in statistical mechanics. The model consists of discrete variables that represent magnetic dipole moments of atomic "spins" that can be in one of two states. The spins are arranged in a graph, usually a lattice, allowing each spin to interact with its neighbors. Neighboring spins that agree have a lower energy than those that disagree; the system tends to the lowest energy but heat disturbs this tendency, thus creating the possibility of different structural phases. The model allows the identification of phase transitions as a simplified model of reality. The two-dimensional square-lattice Ising model is one of the simplest statistical models to show a phase transition.

<span class="mw-page-title-main">Equipartition theorem</span> Theorem in classical statistical mechanics

In classical statistical mechanics, the equipartition theorem relates the temperature of a system to its average energies. The equipartition theorem is also known as the law of equipartition, equipartition of energy, or simply equipartition. The original idea of equipartition was that, in thermal equilibrium, energy is shared equally among all of its various forms; for example, the average kinetic energy per degree of freedom in translational motion of a molecule should equal that in rotational motion.

<span class="mw-page-title-main">Drude model</span> Model of electrical conduction

The Drude model of electrical conduction was proposed in 1900 by Paul Drude to explain the transport properties of electrons in materials. Basically, Ohm's law was well established and stated that the current J and voltage V driving the current are related to the resistance R of the material. The inverse of the resistance is known as the conductance. When we consider a metal of unit length and unit cross sectional area, the conductance is known as the conductivity, which is the inverse of resistivity. The Drude model attempts to explain the resistivity of a conductor in terms of the scattering of electrons by the relatively immobile ions in the metal that act like obstructions to the flow of electrons.

In physics, the S-matrix or scattering matrix is a matrix which relates the initial state and the final state of a physical system undergoing a scattering process. It is used in quantum mechanics, scattering theory and quantum field theory (QFT).

In solid-state physics, the free electron model is a quantum mechanical model for the behaviour of charge carriers in a metallic solid. It was developed in 1927, principally by Arnold Sommerfeld, who combined the classical Drude model with quantum mechanical Fermi–Dirac statistics and hence it is also known as the Drude–Sommerfeld model.

<span class="mw-page-title-main">Wiedemann–Franz law</span> Law of physics

In physics, the Wiedemann–Franz law states that the ratio of the electronic contribution of the thermal conductivity (κ) to the electrical conductivity (σ) of a metal is proportional to the temperature (T).

The diffusion of plasma across a magnetic field was conjectured to follow the Bohm diffusion scaling as indicated from the early plasma experiments of very lossy machines. This predicted that the rate of diffusion was linear with temperature and inversely linear with the strength of the confining magnetic field.

<span class="mw-page-title-main">Toroidal coordinates</span>

Toroidal coordinates are a three-dimensional orthogonal coordinate system that results from rotating the two-dimensional bipolar coordinate system about the axis that separates its two foci. Thus, the two foci and in bipolar coordinates become a ring of radius in the plane of the toroidal coordinate system; the -axis is the axis of rotation. The focal ring is also known as the reference circle.

Neutral-beam injection (NBI) is one method used to heat plasma inside a fusion device consisting in a beam of high-energy neutral particles that can enter the magnetic confinement field. When these neutral particles are ionized by collision with the plasma particles, they are kept in the plasma by the confining magnetic field and can transfer most of their energy by further collisions with the plasma. By tangential injection in the torus, neutral beams also provide momentum to the plasma and current drive, one essential feature for long pulses of burning plasmas. Neutral-beam injection is a flexible and reliable technique, which has been the main heating system on a large variety of fusion devices. To date, all NBI systems were based on positive precursor ion beams. In the 1990s there has been impressive progress in negative ion sources and accelerators with the construction of multi-megawatt negative-ion-based NBI systems at LHD (H0, 180 keV) and JT-60U (D0, 500 keV). The NBI designed for ITER is a substantial challenge (D0, 1 MeV, 40 A) and a prototype is being constructed to optimize its performance in view of the ITER future operations. Other ways to heat plasma for nuclear fusion include RF heating, electron cyclotron resonance heating (ECRH), ion cyclotron resonance heating (ICRH), and lower hybrid resonance heating (LH).

In 3-dimensional topology, a part of the mathematical field of geometric topology, the Casson invariant is an integer-valued invariant of oriented integral homology 3-spheres, introduced by Andrew Casson.

In many-body theory, the term Green's function is sometimes used interchangeably with correlation function, but refers specifically to correlators of field operators or creation and annihilation operators.

Resonance fluorescence is the process in which a two-level atom system interacts with the quantum electromagnetic field if the field is driven at a frequency near to the natural frequency of the atom.

<span class="mw-page-title-main">Viscoplasticity</span> Theory in continuum mechanics

Viscoplasticity is a theory in continuum mechanics that describes the rate-dependent inelastic behavior of solids. Rate-dependence in this context means that the deformation of the material depends on the rate at which loads are applied. The inelastic behavior that is the subject of viscoplasticity is plastic deformation which means that the material undergoes unrecoverable deformations when a load level is reached. Rate-dependent plasticity is important for transient plasticity calculations. The main difference between rate-independent plastic and viscoplastic material models is that the latter exhibit not only permanent deformations after the application of loads but continue to undergo a creep flow as a function of time under the influence of the applied load.

Dynamical mean-field theory (DMFT) is a method to determine the electronic structure of strongly correlated materials. In such materials, the approximation of independent electrons, which is used in density functional theory and usual band structure calculations, breaks down. Dynamical mean-field theory, a non-perturbative treatment of local interactions between electrons, bridges the gap between the nearly free electron gas limit and the atomic limit of condensed-matter physics.

In statistical mechanics, the Griffiths inequality, sometimes also called Griffiths–Kelly–Sherman inequality or GKS inequality, named after Robert B. Griffiths, is a correlation inequality for ferromagnetic spin systems. Informally, it says that in ferromagnetic spin systems, if the 'a-priori distribution' of the spin is invariant under spin flipping, the correlation of any monomial of the spins is non-negative; and the two point correlation of two monomial of the spins is non-negative.

In quantum mechanics, a quantum speed limit (QSL) is a limitation on the minimum time for a quantum system to evolve between two distinguishable (orthogonal) states. QSL theorems are closely related to time-energy uncertainty relations. In 1945, Leonid Mandelstam and Igor Tamm derived a time-energy uncertainty relation that bounds the speed of evolution in terms of the energy dispersion. Over half a century later, Norman Margolus and Lev Levitin showed that the speed of evolution cannot exceed the mean energy, a result known as the Margolus–Levitin theorem. Realistic physical systems in contact with an environment are known as open quantum systems and their evolution is also subject to QSL. Quite remarkably it was shown that environmental effects, such as non-Markovian dynamics can speed up quantum processes, which was verified in a cavity QED experiment.

References

  1. Lawson, J. D. (December 1955). Some criteria for a useful thermonuclear reactor (PDF) (Technical report). Atomic Energy Research Establishment, Harwell, Berkshire, U. K.[ dead link ]
  2. 1 2 Lawson, J. D. (December 1955). "Some Criteria for a Power Producing Thermonuclear Reactor". Proceedings of the Physical Society, Section B. 70 (1): 6–10. Bibcode:1957PPSB...70....6L. doi:10.1088/0370-1301/70/1/303.
  3. "Scientists Achieved Self-Sustaining Nuclear Fusion… But Now They Can't Replicate It". Sciencealert. August 16, 2022.
  4. Abu-Shawareb, H.; Acree, R.; Adams, P.; Adams, J.; Addis, B.; Aden, R.; Adrian, P.; Afeyan, B. B.; Aggleton, M.; Aghaian, L.; Aguirre, A.; Aikens, D.; Akre, J.; Albert, F.; Albrecht, M. (2022-08-08). "Lawson Criterion for Ignition Exceeded in an Inertial Fusion Experiment". Physical Review Letters. 129 (7): 075001. Bibcode:2022PhRvL.129g5001A. doi: 10.1103/PhysRevLett.129.075001 . hdl: 10044/1/99300 . ISSN   0031-9007. PMID   36018710.
  5. 1 2 Spitzer, Lyman; Seeger, Raymond J. (November 1963). "Physics of Fully Ionized Gases". American Journal of Physics. 31 (11): 890–891. Bibcode:1963AmJPh..31..890S. doi:10.1119/1.1969155. ISSN   0002-9505.
  6. "Energy Converter". www.phys.ksu.edu. Kansas State University. Retrieved 2023-02-17.
  7. Wesson, J. (2004). "Tokamaks". Oxford Engineering Science Series (48) (3 ed.). Oxford: Clarendon Press.
  8. "World Highest Fusion Triple Product Marked in High-βp H-mode Plasmas". Archived from the original on 2013-01-06.
  9. Hirsch, Robert L. (October 1967). "Inertial-Electrostatic Confinement of Ionized Fusion Gases". Journal of Applied Physics. 38 (11): 4522–4534. Bibcode:1967JAP....38.4522H. doi:10.1063/1.1709162. ISSN   0021-8979.
  10. Bussard, Robert W (2006-10-02). "The Advent of Clean Nuclear Fusion: Superperformance Space Power and Propulsion". 57th International Astronautical Congress. Reston, Virigina: American Institute of Aeronautics and Astronautics. doi:10.2514/6.iac-06-d2.8.05. ISBN   978-1-62410-042-0.
  11. Rider, Todd H. (1997-04-01). "Fundamental limitations on plasma fusion systems not in thermodynamic equilibrium". Physics of Plasmas. 4 (4): 1039–1046. Bibcode:1997PhPl....4.1039R. doi:10.1063/1.872556. hdl: 1721.1/11412 . ISSN   1070-664X.