09.12.2012 Views

1 - Nuclear Sciences and Applications - IAEA

1 - Nuclear Sciences and Applications - IAEA

1 - Nuclear Sciences and Applications - IAEA

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

s.o<br />

• .5<br />

(a)<br />

<strong>IAEA</strong>-CN-50/A-II-l 91<br />

360<br />

0 60 120 180 2*0 300 360<br />

POLOIDAL ANGLE (DEGREES)<br />

FIG. 5. Resistive mode calculation for n = 1 beta collapse (55791): (a) calculated perturbation of<br />

poloidal magnetic flux; (b) phase <strong>and</strong> amplitude of Mirnov oscillations versus poloidal angle, comparing<br />

measured (•) <strong>and</strong> calculated (o) values.<br />

The low-n modes which cause beta saturation <strong>and</strong> collapse cannot<br />

be attributed to ideal modes. Their growth rates are slow (typically<br />

r ~ 10 4 TAifven) <strong>and</strong> they can persist unchanged for times approaching the<br />

resistive time of the plasma. Ideal stability calculations using GATO [13]<br />

(discussed further below) show that the point where the n = 4, 3, <strong>and</strong> 2<br />

modes are first observed is usually far from the kink mode stability limit.<br />

These low-n modes are most consistent with pressure driven resistive<br />

instabilities. The three-dimensional resistive code CART—II [14] was used<br />

to analyze a discharge similar to that of Fig. 4, at a time during the large<br />

n = 1 oscillation. The nonlinear calculation predicts a predominantly<br />

?n = 2/n = 1 instability, as shown in Fig. 5. The predicted mode structure<br />

agrees well with soft X-ray <strong>and</strong> magnetic fluctuation measurements,<br />

including the coupling to higher poloidal mode numbers which dominates<br />

the structure along the inboard wall. The nonlinear calculation was halted<br />

when the calculated value of Bg/Be at the outer wall was 1%, equal to<br />

the measured value at the time when the n = 1 mode in the discharge of<br />

Fig. 4 stopped rotating <strong>and</strong> the rate of energy loss accelerated. At this<br />

point, the width of the 2/1 isl<strong>and</strong> was calculated to be 11 cm <strong>and</strong> flux<br />

surfaces outside of the q = 2 surface had started to become stochastic.

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

Saved successfully!

Ooh no, something went wrong!