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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Plasma Phys. Control. Fusion 53 (2011) 093001<br />

Topical Review<br />

Figure 29. Comparison <strong>of</strong> 2D MHD simulations with schlieren images (left <strong>of</strong> each pair on top<br />

line). Lower line images are simulations at 33 ns. Reprinted with permission from [197]. Copyright<br />

1997, American Institute <strong>of</strong> Physics.<br />

is despite <strong>the</strong> fact that linear stability <strong>the</strong>ory states that <strong>the</strong> growth rate increases with k,<br />

approximately as k 1/2 . One reason for this can be deduced by considering <strong>the</strong> nonlinear<br />

ρ(v ·∇)v term in <strong>the</strong> fluid equation <strong>of</strong> motion. The m = 0 mode leads to <strong>the</strong> development <strong>of</strong><br />

necking regimes where <strong>the</strong>re is a high J × B <strong>pinch</strong>ing force and bulges which expand almost<br />

freely with a velocity related to <strong>the</strong> sound speed. For <strong>the</strong> bulges if we assume that <strong>the</strong> nonlinear<br />

convective term balances <strong>the</strong> pressure gradient and that <strong>the</strong> J × B force is much weaker, i.e.<br />

ρ(v ·∇)v ∼ = −∇p, (3.42)<br />

it follows that v ∼ = cs , <strong>the</strong> sound speed (Ɣp/ρ) 1/2 , and that <strong>the</strong> radial displacement will be<br />

proportional to c s t. The mode evolves to this nonlinear state when ei<strong>the</strong>r ξ r is ∼βa or k −1 ,<br />

whichever is <strong>the</strong> smaller. (The value <strong>of</strong> β depends on <strong>the</strong> radial pr<strong>of</strong>ile <strong>of</strong> ξ r in <strong>the</strong> linear<br />

phase, e.g. in figure 27, and is approximately 0.5.) This means that <strong>the</strong> fastest growing short<br />

wavelength modes will become nonlinear earlier and will change from exponential growth to<br />

algebraic, in this case a linear dependence on t. If this occurs at time t 1 , and <strong>the</strong> linear growth<br />

rate γf is approximately (ka) 1/2 c s /a, we can write<br />

ξ k (t 1 ) = ξ 0k exp[(ka) 1/2 c s t 1 /a] = k −1 t 1 <strong>the</strong> displacement is given by<br />

ξ k = α 1 c s (t − t 1 ) + k −1 . (3.44)<br />

48

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

Saved successfully!

Ooh no, something went wrong!