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 />

The first and second models are given in detail in [97] in which <strong>the</strong> transport coefficients in<br />

a magnetized plasma <strong>of</strong> Marshall [99] were employed. An asymmetry in heat flow is introduced<br />

by <strong>the</strong> inclusion <strong>of</strong> <strong>the</strong> enthalpy flow in <strong>the</strong> heat-flow equation and <strong>the</strong> <strong>the</strong>rmoelectric field in<br />

Ohm’s law, which are<br />

q z =−κ ∂T<br />

∂z + ψJ z (2.31)<br />

J z = σE z + φ ∂T<br />

(2.32)<br />

∂z<br />

when <strong>the</strong> coefficients are κ = λT 5/2 , ψ =−ζT, σ = αT 3/2 and φ = αγT 3/2 , λ, ζ , α<br />

and γ being constant apart from a weak ln dependence. The energy balance equation with<br />

bremsstrahlung loss included is<br />

J · E = ∂q z<br />

∂z + β bn 2 e T 1/2 + 1 ∂<br />

r ∂r (rq r), (2.33)<br />

where β b is <strong>the</strong> bremsstrahlung coefficient.<br />

On averaging equation (2.33) over radius, i.e. integrating ∫ a<br />

0<br />

2πr dr to <strong>the</strong> outer radius a<br />

at which <strong>the</strong> radial heat flux q r is assumed to be zero, it becomes<br />

N 2<br />

πa 2 ∂q z<br />

∂z + 4 3 β b<br />

πa T 1/2 = IE 2 z (2.34)<br />

On introducing equations (2.7), (2.31) and (2.32) to eliminate N(z), E z (z) and q z (z) a second<br />

order, nonlinear ordinary differential equation for T(z)is found<br />

πa 2 λ d ( )<br />

5/2<br />

dT<br />

T + I(ζ − γ) dT<br />

dz dz<br />

dz + I 2 (<br />

1 − I 2 )<br />

πa 2 aT 3/2 IPB<br />

2 = 0, (2.35)<br />

where <strong>the</strong> Pease–Braginskii current I PB [8, 9] is defined by<br />

( ) 192 1/2<br />

πk<br />

I PB =<br />

= 0.433(ln ) 1/2 × 10 6 A (2.36)<br />

αβ b µ 0<br />

and is derived later in section 2.8.<br />

The variables can be changed to ξ and x, defined by<br />

ξ = πa2 T 5/2<br />

, (2.37)<br />

z 0 I<br />

x = z (2.38)<br />

z 0<br />

so that equation (2.35) becomes<br />

2<br />

5 αλξ d2 ξ<br />

dx 2 + 4 25 αλ ( dξ<br />

dx<br />

) 2<br />

+ 2 α(ζ − γ)dξ<br />

5 dx +1− I 2<br />

IPB<br />

2 = 0. (2.39)<br />

When Marshall’s transport coefficients are employed it is found <strong>the</strong>re is a limiting negative<br />

value <strong>of</strong> dξ/dx in <strong>the</strong> domain <strong>of</strong> interest. This can be seen in a plot <strong>of</strong> <strong>the</strong> first term <strong>of</strong> equation<br />

(2.39) versus dξ/dx in figure 6 where (d 2 ξ/dx 2 ) → 0as(dξ/dx) →−(5/4λ)(ζ −γ)(1−s) =<br />

ξ ′ 1<br />

where s is given by<br />

[<br />

(<br />

4λ<br />

s = 1 −<br />

1 − I 2 )] 1/2<br />

α(ζ − γ) 2 IPB<br />

2 (2.40)<br />

<strong>the</strong> coefficient 4λ/α(ζ − γ) 2 having a value <strong>of</strong> 1/1.47. This limitation means that <strong>the</strong> heat flux<br />

to <strong>the</strong> cathode is very small and negative, while all <strong>the</strong> Joule heating apart from that radiated<br />

is transferred to <strong>the</strong> anode. Equation (2.39) can be analytically integrated once to give ξ as a<br />

17

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

Saved successfully!

Ooh no, something went wrong!