Search

A new approach to high pressure magnetically-confined plasmas is necessary to design efficient fusion devices. This paper presents an equilibrium combining two solutions of the Grad-Shafranov equation, which describes the magnetohydrodynamic equilibrium in toroidal geometry. The outer equilibrium is paramagnetic and confines the inner equilibrium, whose strong diamagnetism permits to balance large pressure gradients. The existence of both equilibria in the same volume yields a dual equilibrium structure. Their combination also improves free-boundary mode stability.

Because experiment/model comparisons in magnetic confinement fusion have not yet satisfied the requirements for validation as understood broadly, a set of approaches to validating mathematical models and numerical algorithms are recommended as good practices. Previously identified procedures, such as verification, qualification, and analysis of error and uncertainty, remain important. However, particular challenges intrinsic to fusion plasmas and physical measurement therein lead to identification of new or less familiar concepts that are also critical in validation. These include the primacy hierarchy, which tracks the integration of measurable quantities, and sensitivity analysis, which assesses how model output is apportioned to different sources of variation. The use of validation metrics for individual measurements is extended to multiple measurements, with provisions for the primacy hierarchy and sensitivity. This composite validation metric is essential for quantitatively evaluating comparisons with experiments. To mount successful and credible validation in magnetic fusion, a new culture of validation is envisaged., Comment: 27 pages, 1 table, 6 figures

One approach to improving the real-time efficiency of plasma turbulence calculations is to use a parallel algorithm. A parallel algorithm for plasma turbulence calculations was tested on the Intel iPSC/860 hypercube and the Touchtone Delta machine. Using the 128 processors of the Intel iPSC/860 hypercube, a factor of 5 improvement over a single-processor CRAY-2 is obtained. For the Touchtone Delta machine, the corresponding improvement factor is 16. For plasma edge turbulence calculations, an extrapolation of the present results to the Intel (sigma) machine gives an improvement factor close to 64 over the single-processor CRAY-2.

Electron instabilities of magnetized spherical shell distributions in velocity space with a colder Maxwellian background are investigated analytically with simulations using electrostatic particle codes. The resonant and nonresonant instabilities observed in the particle simulations are in agreement with zeros of the dielectric function, as found from the resonant approximations for waves with an electric field component along the magnetic field or by computation with a root solver code in the case of perpendicular propagation. Saturation of the instabilities is by nonlinear cyclotron resonance with the cold background in the resonant case or by nonstochastic cyclotron harmonic damping by the cold background in the nonresonant case. Instabilities invariably lead to perpendicular acceleration and heating of the cold background to velocities sometimes exceeding the shell velocity.

The physical processes of a linearly polarized electromagnetic pulse of highly relativistic amplitude in an underdense plasma accelerating particles to very high energies are studied through computer simulation. An electron-positron plasma is considered first. The maximum momenta achieved scale as the square of the wave amplitude. This acceleration stops when the bulk of the wave energy is converted to particle energy. The pulse leaves behind as a wake a vacuum region whose length scales as the amplitude of the wave. The results can be explained in terms of a snow plow or piston-like action of the radiation on the plasma. When a mass ratio other than unity is chosen and electrostatic effects begin to play a role, first the ion energy increases faster than the electron energy and then the electron energy catches up later, eventually reaching the same value.

Global three-dimensional computer simulations of the magnetosphere using a particle MHD code, reproduce the steady-state Dungey magnetospheric topology in three dimensions. The formation of a compression zone downstream of the tail neutral line that is probably bounded by wake shocks is observed. This compression zone changes its cross-section with distance downstream.

Results are presented of a computer simulation study of electron cyclotron harmonic waves; an electrostatic finite size particle simulation code is used. The initial electron velocity distribution function is modeled as a hot ring perpendicular to the ambient magnetic field and a cold Maxwellian component. Nonlinear cyclotron resonance is found to be the saturation mechanism of the unstable waves and the heating mechanism of the cold electrons. The results can be related to certain ionospheric phenomena.

The time-dependent magnetohydrodynamic interaction of the solar wind with a two-dimensional dipole magnetic field has been simulated using a novel Lagrangian particle type of MHD code that can treat local low density or vacuum regions without numerical instability. This enables one to simulate the time-dependent magnetic tail. When the solar wind field is southward, a magnetic field line topology consistent with Dungey's model emerges in steady state. The tail, however, is short, and the x-points are only slightly shifted from their vacuum locations, because of strong numerical resistivity. Different configurations resulting from different relative orientations of the solar wind magnetic field and dipole axis are also presented. While the magnetic field is relatively steady, the density and flow in the magnetosheath are turbulent, as are the bow shock and magnetopause; the Kelvin-Helmholtz instability may account for these phenomena. We also model a 'substorm' as the passage of a rotational discontinuity in the solar wind over the dipole. Both 90 and 180 deg shifts to a southward solar wind field cause a violent readjustment of the magnetic tail which eventually settles down to the Dungey configuration.

The paper presents preliminary results from time-dependent two-dimensional numerical modelling of the magnetohydrodynamic interaction of the solar wind with the magnetosphere. A southward solar wind-field produces a magnetospheric topology consistent with Dungey's (1961) model. The interaction appears to be fundamentally unsteady; the shock, magnetosheath, and magnetopause are highly turbulent. A 'substorm' is modelled as the passage of a rotational discontinuity over the magnetosphere; the onset of enhanced reconnection in the magnetospheric tail produces a closed magnetic island which convects downstream.