Plasma scaling

Range of plasmas. Density increases upwards, temperature increases towards the right. The free electrons in a metal may be considered an electron plasma

The parameters of plasmas, including their spatial and temporal extent, vary by many orders of magnitude. Nevertheless, there are significant similarities in the behaviors of apparently disparate plasmas. Understanding the scaling of plasma behavior is of more than theoretical value. It allows the results of laboratory experiments to be applied to larger natural or artificial plasmas of interest. The situation is similar to testing aircraft or studying natural turbulent flow in wind tunnels with smaller-scale models.

Similarity transformations (also called similarity laws) help us work out how plasma properties changes in order to retain the same characteristics. A necessary first step is to express the laws governing the system in a nondimensional form. The choice of nondimensional parameters is never unique, and it is usually only possible to achieve by choosing to ignore certain aspects of the system.

One dimensionless parameter characterizing a plasma is the ratio of ion to electron mass. Since this number is large, at least 1836, it is commonly taken to be infinite in theoretical analyses, that is, either the electrons are assumed to be massless or the ions are assumed to be infinitely massive. In numerical studies the opposite problem often appears. The computation time would be intractably large if a realistic mass ratio were used, so an artificially small but still rather large value, for example 100, is substituted. To analyze some phenomena, such as lower hybrid oscillations, it is essential to use the proper value.

A commonly used similarity transformation

One commonly used similarity transformation was derived for gas discharges by James Dillon Cobine (1941), Alfred Hans von Engel and Max Steenbeck (1934), and further applied by Hannes Alfvén and Carl-Gunne Fälthammar to plasmas. They can be summarised as follows:

Similarity Transformations Applied to Gaseous Discharges and some Plasmas

Property	Scale Factor
length, time, inductance, capacitance	x1	x
particle energy, velocity, potential, current, resistance	x0=1	Unchanged
electric and magnetic fields, conductivity, neutral gas density, ionization fraction	x−1	1/x
current density, electron and ion densities	x−2	1/x2

This scaling applies best to plasmas with a relatively low degree of ionization. In such plasmas, the ionization energy of the neutral atoms is an important parameter and establishes an absolute energy scale, which explains many of the scalings in the table:

• Since the masses of electrons and ions cannot be varied, the velocities of the particles are also fixed, as is the speed of sound.
• If velocities are constant, then time scales must be directly proportional to distance scales.
• In order that charged particles falling through an electric potential gain the same energy, the potentials must be invariant, implying that the electric field scales inversely with the distance.
• Assuming that the magnitude of the E-cross-B drift is important and should be invariant, the magnetic field must scale like the electric field, namely inversely with the size. This is also the scaling required by Faraday's law of induction and Ampère's law.
• Assuming that the speed of the Alfvén wave is important and must remain invariant, the ion density (and with it the electron density) must scale with B², that is, inversely with the square of the size. Considering that the temperature is fixed, this also ensures that the ratio of thermal to magnetic energy, known as beta, remains constant. Furthermore, in regions where quasineutrality is violated, this scaling is required by Gauss's law.
• Ampère's law also requires that current density scales inversely with the square of the size, and therefore that current itself is invariant.
• The electrical conductivity is current density divided by electric field and thus scales inversely with the length.
• In a partially ionized plasma, the electrical conductivity is proportional to the electron density and inversely proportional to the neutral gas density, implying that the neutral density must scale inversely with the length, and ionization fraction scales inversely with the length.

Limitations

While these similarity transformations capture some basic properties of plasmas, not all plasma phenomena scale in this way. Consider, for example, the degree of ionization, which is dimensionless and thus would ideally remain unchanged when the system is scaled. The number of charged particles per unit volume is proportional to the current density, which scales as x ^-2, whereas the number of neutral particles per unit volume scales as x ^-1 in this transformation, so the degree of ionization does not remain unchanged but scales as x ^-1.

Astrophysical application

As an example, take an auroral sheet with a thickness of 1 km. A laboratory simulation might have a thickness of 10 cm, a factor of 10⁴ smaller. To satisfy the condition of this similarity transformation, the gaseous density would have to be increased by a factor of 10⁴ from 10⁴ m⁻³ to 10⁸ m⁻³ (10¹⁰ cm⁻³ to 10¹⁴ cm⁻³), and the magnetic field would have to be increased by the same factor from 50 microteslas to 500 milliteslas (0.5 gauss to 5 kilogauss). These values are large but within the range of technology. If the experiment captures the essential features of the aurora, the processes will be 10⁴ times faster so that a pulse that takes 100 s in nature would take only 10 ms in the laboratory.

Similarity transformations applied to some astrophysical plasmas
Actual plasma properties compared to a laboratory plasma if the scale length is reduced to 10 cm.

Region	Characteristic dimension (cm)			Density (particles/cm3)		Magnetic field (gauss)		Characteristic time
	Actual	Scaled	Scale Factor	Actual	Scaled	Actual	Scaled	Actual	Scaled
Ionosphere	106 - 107	10	10−5 - 10−6	1010	1015 - 1016	0.5	5x104 - 5x105	Period of Giant pulsation
								100 s	0.1 - 1 ms
Exosphere	109	10	10−8	105 - 10	1013 - 109	0.5 - 5x10−4	5x107 - 5x104	One Day
								105 s	1 ms
Interplanetary space	1013	10	10−12	1 - 10	1012 - 1013	10−4	108	One Solar Rotation
								2x106 s	2 μs
Interstellar medium	3x1022	10	3x10−22	1	3x1021	10−6 - 10−5	3x1015 - 3x1016	Period of galactic rotation
								1x1016 s	3 μs
Intergalactic space	>3x1027	10	<3x10−27	10−4?	>3x1022	10−7?	>3x1019	Age of the Universe
								4x1017s	1x10−9s
Solar chromosphere	108	10	10−7	1011 - 1014	1018 - 1021	103 - 1	1010 - 107	Life of Solar Flare
								103 s	100 μs
								Life of Solar Prominence
								105 s	10 ms
Solar corona	1010 - 1011	10	10−9 - 10−10	108 - 106	1017 - 1016	102 - 10−1	1011 - 109	Life of Coronal Arc
								103 s	10−1 to 1 µs
								Solar Cycle
								22 years	70 to 700 ms

Particle density of the Earth's atmosphere at sea level is 10¹⁹ per cm³
Small bar magnet = 100 milliteslas. Big electromagnet = 2 teslas
10⁹ cm = 10,000 km

The table shows the properties of some actual space plasma (see the columns labelled Actual). It also shows how other plasma properties would need to be changed, if (a) the characteristic length of a plasma were reduced to just 10 cm, and (b) the characteristics of the plasma were to remain unchanged.

The first thing to notice is that many cosmic phenomena cannot be reproduced in the laboratory because the necessary magnetic field strength is beyond the technological limits. Of the phenomena listed, only the ionosphere and the exosphere can be scaled to laboratory size. Another problem is the ionization fraction. When the size is varied over many orders of magnitude, the assumption of a partially ionized plasma may be violated in the simulation. A final observation is that the plasma densities needed in the laboratory are sizeable, up to 10¹⁶ cm⁻³ for the ionosphere, compared to the atmospheric density of about 10¹⁹ particles per cm³. In other words, the laboratory analogy of a low density space plasma is not a "vacuum chamber", but laboratory plasma with a pressure, when the higher temperature is taken into consideration, which can approach atmospheric pressure.

Dimensionless parameters in tokamaks

One of the central questions in fusion power research is to predict the energy confinement time in machines that are larger than any that have ever been built. A widely accepted approach to doing this is to express the scaling in terms of nondimensional parameters. Geometrical parameters, such as the ratio of the major to the minor radius, the shape of the plasma cross section, and the angle of the magnetic field, can be chosen in current experiments to equal the value desired for a full scale reactor. The remaining (dimensional) parameters can be taken to be the particle density n, the temperature T, the magnetic field B, and the size (major radius) R. These can be combined into the three dimensionless parameters β (the ratio of plasma pressure to magnetic pressure), ν^* (the product of the collision frequency and the thermal transit time), and ρ^* (the ratio of the Larmor radius to the torus radius). These have the following scalings:

β ~ nTB ^-2

ν^* ~ nT ^-2R

ρ^* ~ T ^1/2B ^-1R ^-1

The radius R can be varied while keeping these three parameters constant if n, T, and B are scaled in this way:

n ~ R ^-2

T ~ R ^-1/2

B ~ R ^-5/4

Note that this similarity transformation is distinct from that considered above, which would yield n ~ R ^-1, T ~ R ⁰, and B ~ R ^-1. This is because the physical effects to be studied are different.

The scaling of the magnetic field with the minus 5/4 power of the size implies that a 1:3 scale model of a power-producing tokamak with a magnetic field of 10 T at the coils would require a field of 30 T, which is technologically infeasible.

The next best alternative is to allow ρ^* to vary and to extrapolate according to the dependence found. ρ^* is the parameter considered least likely to harbor surprises, partly for theoretical considerations, but also simply because it is, in contrast to β and ν^*, already much larger than unity. This can be done in a single machine (constant R) by varying the magnetic field and scaling density and temperature as:

n ~ B ^4/3

T ~ B ^2/3

It should be kept in mind that the assumption has been made that the important turbulent transport processes depend only on the parameters chosen. It is only physical reasoning, not mathematical necessity, that concludes that the ratio of the torus radius to the Larmor radius is important, and not, for example, the ratio to the Debye length. In the same way, it has been assumed that the absolute energy levels of atomic physics do not dictate an absolute temperature dependence, or equivalently, that the boundary layer where atomic physics is important, is small enough not to determine the overall energy confinement.

See also

• Similitude (model)