You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Plasma Phys. Control. Fusion 53 (2011) 093001<br />
Topical Review<br />
3.13. Micro-instabilities and anomalous resistivity<br />
Some form <strong>of</strong> micro- or velocity-space-instability will be triggered when <strong>the</strong> drift velocity<br />
v d (≈ J/n e e) <strong>of</strong> <strong>the</strong> current-carrying electron exceeds some critical velocity v crit , which is<br />
typically <strong>of</strong> order <strong>the</strong> sound speed, c s (= √ [(ZT e + T i )/m i ]). The instability will rapidly<br />
grow in time and saturate in amplitude. The resulting electrostatic or electromagnetic turbulent<br />
wave structure will interact with <strong>the</strong> drifting electrons, causing <strong>the</strong>m to be scattered. This<br />
scattering is usually much larger than that caused by Coulomb collisions with ions; hence an<br />
anomalous resistivity will be generated.<br />
The microinstability that is perhaps most relevant in <strong>the</strong> Z-<strong>pinch</strong> configuration is <strong>the</strong> lower<br />
hybrid drift instability, because <strong>the</strong> current and drift velocity are orthogonal to <strong>the</strong> magnetic<br />
field. It is related to <strong>the</strong> modified two-stream instability studied earlier by Buneman [221] and<br />
Ashby and Paton [222]. Krall and Liewer [223] derived a linear electrostatic kinetic model <strong>of</strong><br />
<strong>the</strong> lower hybrid instability. It does not have a critical velocity threshold, but nonlinearly its<br />
effective collision frequency varies as <strong>the</strong> square <strong>of</strong> <strong>the</strong> electron drift velocity.<br />
A quasi-linear model to find a value for <strong>the</strong> anomalous resistivity was developed by<br />
Davidson and Gladd [224]. Huba and Papdopoulos [225] considered electron resonance<br />
broadening as a saturation mechanism especially in finite β plasmas. Particle simulations<br />
by Winske and Liewer [226] in contrast find that ion trapping is a saturating mechanism in<br />
<strong>the</strong> high drift velocity regime, while Davidson [227] found current relaxation and plateau<br />
formation can cause saturation at low drift velocity.<br />
Probably <strong>the</strong> two most relevant papers to yield in 2D nonlinear modelling a value for <strong>the</strong><br />
anomalous collision frequency at saturation are Drake et al [228] and Brackbill et al [229].<br />
The former considered nonlinear mode–mode coupling <strong>of</strong> <strong>the</strong> electrostatic waves showing<br />
that energy is transferred from long wavelength modes to short wavelength modes which are<br />
Landau damped. The anomalous collision frequency ν anom so found from this model is<br />
ν anom = 2.4( e i ) 1/2 (v d /v i ) 2 . (3.54)<br />
Here ( e i ) 1/2 , where e and i are <strong>the</strong> cyclotron angular frequencies <strong>of</strong> <strong>the</strong> electrons and<br />
ions, respectively, is <strong>the</strong> lower hybrid resonance frequency, while v d and v i are <strong>the</strong> electron<br />
drift velocity and <strong>the</strong> ion <strong>the</strong>rmal speed, respectively.<br />
Brackbill et al [229] employed in contrast a 2D implicit e.m. particle-in-cell code, VENUS,<br />
initially with a Harris equilibrium. The saturation mechanism was considered to be electron<br />
kinetic dissipation with wavelengths greater than <strong>the</strong> electron Larmor radius. The resulting<br />
anomalous collision frequency found in this numerical simulation differs only slightly from<br />
equation (3.43), namely<br />
ν anom =<br />
C ( e i ) 1/2 (v d /v i ) 2 , (3.55)<br />
4πβ i<br />
where C is 0.13–0.38 and β i is <strong>the</strong> ratio <strong>of</strong> ion pressure to magnetic pressure.<br />
Before applying <strong>the</strong>se results to <strong>the</strong> Z-<strong>pinch</strong> it is perhaps worthwhile to recall <strong>the</strong> basic<br />
physics underlying <strong>the</strong> lower hybrid resonance. Although <strong>the</strong> ion cyclotron frequency is<br />
present in <strong>the</strong> formula, <strong>the</strong> ions are in fact to a good approximation unmagnetized. Figure 31<br />
illustrates how <strong>the</strong> lower hybrid resonance arises when <strong>the</strong> ion plasma frequency is much<br />
higher. Consider a plane e.s. wave propagating in <strong>the</strong> z-direction. The unmagnetized ions<br />
will oscillate to and fro in <strong>the</strong> ±z-direction. With a steady magnetic field in <strong>the</strong> y-direction,<br />
<strong>the</strong> magnetized electrons will have an oscillatory Ẽ/B drift in <strong>the</strong> ±x-direction. Superposed<br />
on this is <strong>the</strong> polarization drift, Ė/(B e ) in <strong>the</strong> ±z-direction, which will convert <strong>the</strong> electron<br />
guiding-centre trajectory into an ellipse. The resonance occurs, i.e. a tendency for zero charge<br />
separation when <strong>the</strong> spatial amplitude <strong>of</strong> <strong>the</strong> ion oscillation eẼ/m i ω 2 equals <strong>the</strong> width <strong>of</strong> <strong>the</strong><br />
ellipse Ẽ/(B e ); this gives ω = ( e i ) v 2<br />
52