14.03.2014 Views

A review of the dense Z-pinch

A review of the dense Z-pinch

A review of the dense Z-pinch

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!