next up previous
Next: V. Inboard Injection Up: MHD Effects on Previous: III. Computational Model

IV. Outboard Injection

All other things being equal, the 3D evolution depends on the initial location of the pellet perturbation in the poloidal plane, as well as the pellet density. In the case of outboard injection, the initial density tex2html_wrap_inline897 is shown in Fig.4(a), in the plane tex2html_wrap_inline899 where the density is a maximum. The pellet is located at the outboard midplane. The density enhancement tex2html_wrap_inline901 where tex2html_wrap_inline903 is the initial, uniform density, is 15. On a rapid time scale, during which there is approximate fast magnetosonic and shear Alfvén equilibration, the pellet expands outwards towards larger major radius.

Next, Fig.4(b) shows the density in the tex2html_wrap_inline899 plane at time tex2html_wrap_inline909 where tex2html_wrap_inline911 is the Alfvén time defined in terms of the average toroidal magnetic field and background density. The pellet both shifts outwards and expands along the magnetic field.

The displacement of the pellet can be measured by calculating the flux surface average tex2html_wrap_inline913 as in (8), using the toroidally averaged 2D magnetic field, its n=0 component. This is acceptable if the initial and asymptotic states have axisymmetric magnetic field. In fact even at intermediate times, tex2html_wrap_inline917 the n = 0 part of tex2html_wrap_inline921 changes relatively little in comparison with the changes in the average density profile. The averaged density tex2html_wrap_inline923 is shown in Fig.5 at times (a) t = 0, (b) t = 30, and (c) t = 96. The maximum shift occurs at about t = 30. The peak of the tex2html_wrap_inline923 profile moves to lower tex2html_wrap_inline917 which means that tex2html_wrap_inline937 has moved closer to the wall boundary at tex2html_wrap_inline939 Afterwards, the peak starts to bounce back towards its initial state. In principle it should return to its initial state, if the density in each magnetic flux tube is conserved. However, there is some loss caused by density diffusivity D and temperature diffusivity tex2html_wrap_inline943 with tex2html_wrap_inline945

The three dimensional evolution of the density can be seen in Fig.6(a,b) and Fig.7(a,b), which are isoplots of the density surface at tex2html_wrap_inline947 of its maximum. Initially the density is localized near toroidal angle tex2html_wrap_inline949 In Fig.6(b)at time t = 51, the pellet has spread more than once around the torus, which gives an idea of the sound wave transit time for this case. In Fig.7(a), at time t = 79, the pellet material has spread out twice around the torus, and forms a ``snake" like structure. The pellet material is not resonant with the q = 2 surface, since the pellet was injected on an initial tex2html_wrap_inline765 surface with q ;SPMgt; 2. Hence it smoothes out as time elapses, as seen in Fig.7(b) at t = 96.


next up previous
Next: V. Inboard Injection Up: MHD Effects on Previous: III. Computational Model

Hank Strauss
Wed Jan 7 14:23:34 EST 1998