Electrothermal instability

From Wikipedia, the free encyclopedia
Jump to: navigation, search

The electrothermal instability (also known as the ionization instability or Velikhov instability in the literature) is a magnetohydrodynamic (MHD) instability appearing in magnetized non-thermal plasmas used in MHD converters. It was first theoretically discovered in 1962 and experimentally measured into a MHD generator in 1964 by Evgeny Velikhov.[1][2][3]

Physical explanation and characteristics[edit]

Evolution of the electrothermal instability in a Faraday MHD converter. Electric current lines.

This instability is a turbulence of the electron gas in a non-equilibrium plasma (i.e. where the electron temperature Te is greatly higher than the overall gas temperature Tg). It arises when a magnetic field powerful enough is applied in such a plasma, reaching a critical Hall parameter βcr.

Locally, the number of electrons and their temperature fluctuate (electron density and thermal velocity) as the electric current and the electric field.

The Velikhov instability is a kind of ionization wave system, almost frozen in the two temperature gas. The reader can evidence such a stationary wave phenomenon just applying a transverse magnetic field with a permanent magnet on the low-pressure control gauge (Geissler tube) provided on vacuum pumps. In this little gas-discharge bulb a high voltage electric potential is applied between two electrodes which generates an electric glow discharge (pinkish for air) when the pressure has become low enough. When the transverse magnetic field is applied on the bulb, some oblique grooves appear in the plasma, typical of the electrothermal instability.

The electrothermal instability occurs extremely quickly, in a few microseconds. The plasma becomes non-homogeneous, transformed into alternating layers of high free electron and poor free electron densities. Visually the plasma appears stratified, as a "pile of plates".

Hall effect in plasmas[edit]

The Hall effect in ionized gases has nothing to do with the Hall effect in solids (where the Hall parameter is always very inferior to unity). In a plasma, the Hall parameter can take any value.

The Hall parameter β in a plasma is the ratio between the electron gyrofrequency Ωe and the electron-heavy particles collision frequency ν:


\beta \, = \, \frac {\Omega_e}{\nu} \, = \, \frac {e\ B}{m_e\ \nu}

where

e is the electron charge (1.6 × 10−19 coulomb)
B is the magnetic field (in teslas)
me is the electron mass (0.9 × 10−30 kg)

The Hall parameter value increases with the magnetic field strength.

Physically, when the Hall parameter is low, the trajectories of electrons between two encounters with heavy particles (neutral or ion) are almost linear. But if the Hall parameter is high, the electron movements are highly curved. The current density vector J is no more colinear with the electric field vector E. The two vectors J and E make the Hall angle θ which also gives the Hall parameter:


\ \beta \, = \, \tan \theta

Plasma conductivity and magnetic fields[edit]

In a non-equilibrium ionized gas with high Hall parameter, Ohm's law,


\mathbf{J} = \sigma\mathbf{E}

where σ is the electrical conductivity (in siemens per metre),

is a matrix, because the electrical conductivity σ is a matrix:


\sigma = \sigma_s \begin{Vmatrix} \dfrac{1}{1+\beta^2} & \dfrac{-\beta}{1+\beta^2} \\ \dfrac{\beta}{1+\beta^2} & \dfrac{1}{1+\beta^2} \end{Vmatrix}

σS is the scalar electrical conductivity:


\sigma_s = \frac {n_e\ e^2}{m_e\ \nu}

where ne is the electron density (number of electrons per cubic meter).

The current density J has two components:


J_{\parallel} = \frac {n_e\ e^2}{m_e\ \nu} \ \frac {E}{1+\beta^2} \qquad \text{and} \qquad J_{\perp} = \frac {-n_e\ e^2}{m_e\ \nu} \ \frac {\beta\ E}{1+\beta^2}

Therefore


J_{\perp} = J_{\parallel}\ \beta

The Hall effect makes electrons "crabwalk".

When the magnetic field B is high, the Hall parameter β is also high, and \frac {1}{1+\beta^2} \ll 1

Thus both conductivities


\sigma_{\parallel} \approx \frac {\sigma_s}{\beta^2} \qquad \text{and} \qquad \sigma_{\perp} \approx \frac {\sigma_s}{\beta}

become weak, therefore the electric current cannot flow in these areas. This explains why the electron current density is weak where the magnetic field is the strongest.

Critical Hall parameter[edit]

The electrothermal instability occurs in a plasma at a (Te > Tg) regime when the Hall parameter is higher that a critical value βcr.

We have


f = \frac{\left( \frac{\delta \mu}{\mu} \right)}{\left( \frac{\delta n_e}{n_e} \right)}

where μ is the electron mobility (in m2/(V·s))

and


s = \frac{2\ k\ T_e^2}{E_i\; (T_e - T_g)} \times \frac {1}{1 + \dfrac{3}{2}\ \dfrac{k\; T_e}{E_i}}

where Ei is the ionization energy (in electron volts) and k the Boltzmann constant.

The growth rate of the instability is


g = \frac{\sigma\ E^2}{n_e\; \left( E_i + \frac{3}{2} k\; T_e \right)\; \left( 1 + \beta^2 \right)}\; (\beta - \beta_{cr})

And the critical Hall parameter is


\beta_{cr} = 1.935 f + 0.065 + s
~

The critical Hall parameter βcr greatly varies according to the degree of ionization α :


\alpha = \frac{n_i}{n_n}

where ni is the ion density and nn the neutral density (in particles per cubic metre).

The electron-ion collision frequency νei is much greater than the electron-neutral collision frequency νen.

Therefore with a weak energy degree of ionization α, the electron-ion collision frequency νei can equal the electron-neutral collision frequency νen.

  • For a weakly ionized gas (non-Coulombian plasma, when νei < νen ):

\beta_{cr} \approx (s^2 + 2s)^{\frac{1}{2}}
  • For a fully ionized gas (Coulombian plasma, when νei > νen ):

\beta_{cr} \approx (2 + s)

NB: The term "fully ionized gas", introduced by Lyman Spitzer, does not mean the degree of ionization is unity, but only that the plasma is Coulomb-collision dominated, which can correspond to a degree of ionization as low as 0.01%.

Technical problems and solutions[edit]

A two-temperature gas, globally cool but with hot electrons (Te >> Tg) is a key feature for practical MHD converters, because it allows the gas to reach sufficient electrical conductivity while protecting materials from thermal ablation. This idea was first introduced for MHD generators in the early 1960s by Jack L. Kerrebrock[4] and Alexander E. Sheindlin.[5]

But the unexpected large and quick drop of current density due to the electrothermal instability ruined many MHD projects worldwide, while previous calculation envisaged energy conversion efficiencies over 60% with these devices. Whereas some studies were made about the instability by various researchers,[6][7] no real solution was found at that time. This prevented further developments of non-equilibrium MHD generators and compelled most engaged countries to cancel their MHD power plants programs and to retire completely from this research field in the early 1970s, because this technical problem was considered as an impassable stumbling block in these days.

Nevertheless experimental studies about the growth rate of the electrothermal instability and the critical conditions showed that a stability region still exists for high electron temperatures.[8][9] The stability is given by a quick transition to "fully ionized" conditions (fast enough to overtake the growth rate of the electrothermal instability) where the Hall parameter decreases cause of the collision frequency rising, below its critical value which is then about 2. Stable operation with several megawatts in power output had been experimentally achieved as from 1967 with high electron temperature.[10][11][12] But this electrothermal control does not allow to decrease Tg low enough for long duration conditions (thermal ablation) so such a solution is not practical for any industrial energy conversion.

Another idea to control the instability would be to increase non-thermal ionisation rate thanks to a laser which would act like a guidance system for streamers between electrodes, increasing the electron density and the conductivity, therefore lowering the Hall parameter under its critical value along these paths. But this concept has never been tested experimentally.

In the 1970s and more recently, some researchers tried to master the instability thanks to oscillating fields. Oscillations of the electric field or of an additional RF electromagnetic field locally modify the Hall parameter.[13][14]

Finally, a solution has been found in the early 1980s to annihilate completely the electrothermal instability within MHD converters, thanks to non-homogeneous magnetic fields. A strong magnetic field implies a high Hall parameter, therefore a low electrical conductivity in the medium. So the idea is to make some "paths" linking an electrode to the other, where the magnetic field is locally attenuated. Then the electric current tends to flow in these low B-field paths as thin plasma cords or streamers, where the electron density and temperature increase. The plasma becomes locally Coulombian, and the local Hall parameter value falls, while its critical threshold is risen. Experiments where streamers do not present any inhomogeneity has been obtained with this method.[15][16] This effect, strongly nonlinear, was unexpected but led to a very effective system for streamer guidance.

But this last working solution was discovered too late, 10 years after all the international effort about MHD power generation had been abandoned in most nations. Vladimir S. Golubev, coworker of Evgeny Velikhov, who met Jean-Pierre Petit in 1983 at the 9th MHD International conference in Moscow, made the following comment to the inventor of the magnetic stabilization method:

You bring the cure, but the patient already died...

However it should be noted that this electrothermal stabilization by magnetic confinement, if found too late for the development of MHD power plants, might be of interest for future applications of MHD to aerodynamics (magnetoplasma-aerodynamics for hypersonic flight).

See also[edit]

External links[edit]

References[edit]

  1. ^ E.P. Velikhov (1962). "Hall instability of current-carrying slightly-ionized plasmas". Newcastle-upon-Tyne, England. 1st International Conference on MHD Electrical Power Generation, Paper 47. 
  2. ^ E.P. Velikhov; A.M. Dykhne (13–18 July 1963). "Plasma turbulence due to the ionization instability in a strong magnetic field". In Paris, France. Proceedings. 6th International Conference on Ionization Phenomena in Gases 4. 
  3. ^ E.P. Velikhov; A.M. Dykhne, I.Ya Shipuk (1965). "Ionization instability of a plasma with hot electrons". Belgrade, Yugoslavia. 7th International Conference on Ionization Phenomena in Gases. 
  4. ^ J.L. Kerrebrock (1 November 1960). "Non-equilibrium effects on conductivity and electrode heat transfer in ionized gases". Technical Note No. 4 (Guggenheim Jet Propulsion Center, Caltech, Pasadena, California.: AFOSR-165). OSTI 4843920
  5. ^ A.E. Sheindlin; V.A. Batenin, E.I. Asinovsky (6 July 1964). "investigation of non-equilibrium ionization in a mixture of argon and potassium". In Paris, France. CONF-640701-102. International symposium on magnetohydrodynamic electric power generation.  OSTI 5024025
  6. ^ A. Solbes (24–30 July 1968). "A quasi linear plane wave study of electrothermal instabilities". In International Atomic Energy Agency, Warsaw, Poland. SM/107/26. 8th International Conference on MHD Electrical Power Generation. 
  7. ^ A.H. Nelson; M.G. Haines (26–28 March 1969). "Analysis of the nature and growth of electrothermal waves". In MIT, Cambridge, MA, USA. Proceedings. 10th Symposium in Engineering Aspects of MHD. doi:10.1088/0032-1028/11/10/003. 
  8. ^ J.P. Petit; J. Valensi, J.P. Caressa (24–30 July 1968). "Theoretical and experimental study in shock tube of non-equilibrium phenomena in a closed-cycle MHD generator". In International Atomic Energy Agency, Warsaw, Poland. Proceedings. 8th International Conference on MHD Electrical Power Generation 2. pp. 745–750. 
  9. ^ J.P. Petit; J. Valensi (1 September 1969). "Growth rate of electrothermal instability and critical Hall parameter in closed-cycle MHD generators when the electron mobility is variable". Comptes rendus de l'Académie des sciences (Paris: French Academy of Sciences) (269): 365–367. 
  10. ^ J.P. Petit; J. Valensi, J.P. Caressa (24–30 July 1968). "Electrical characteristics of a converter using as a conversion fluid a binary mix of rare gases with non-equilibrium ionization". In International Atomic Energy Agency, Warsaw, Poland. Proceedings. 8th International Conference on MHD Electrical Power Generation 3. 
  11. ^ J.P. Petit; J. Valensi, D. Dufresnes, J.P. Caressa (27 January 1969). "Electrical characteristics of a Faraday linear generator using a binary mix of rare gases, with non-equilibrium ionization". Comptes rendus de l'Académie des sciences (Paris: French Academy of Sciences) 268 (A): 245–247. 
  12. ^ S. Hatori; S. Shioda (March 1974). "Stabilization of Ionization Instability in an MHD Generator". Journal of the Physical Society of Japan (Tokyo Institute of Technology, Yokohama, Japan) 36 (3): 920–920. doi:10.1143/JPSJ.36.920. 
  13. ^ G.I. Shapiro; A.H. Nelsone (12 April 1978). "Stabilization of ionization instability in a variable electric field". Pis'ma v Zhurnal Tekhnicheskoi Fiziki (Akademiia Nauk SSSR, Institut Problem Mekhaniki, Moscow, USSR) 4 (12): 393–396. 
  14. ^ T. Murakami; Y. Okuno; H. Yamasaki (December 2005). "Dynamic stabilization of the electrothermal instability". Applied Physics Letters (Tokyo Institute of Technology, Yokohama, Japan) 86 (19): 191502–191502.3. Bibcode:2005ApPhL..86s1502M. doi:10.1063/1.1926410. 
  15. ^ J.P. Petit; M. Billiotte (27 April 1981). "Method for eliminating the Velikhov instability". Comptes-rendus de l'Académie des Sciences (Paris: French Academy of Sciences): 158–161. 
  16. ^ J.P. Petit (September 1983). "Cancellation of the Velikhov instability by magnetic confinment" (PDF). Moscow, Russia. 8th International Conference on MHD Electrical Power Generation.