Sandpiles with bistable automata rules: towards a minimal model of pedestal formation and structure.

A sandpile with two stable and two unstable ranges of slopes is presented as a minimal model for the study of H-mode pedestal formation and dynamics. Pedestals are observed to form and expand inward with increasing deposition. Transport bifurcation is not critical to pedestal formation, though the pedestal structure obtained with a second, hard stability boundary is qualitatively different from that found in standard sandpiles. Nonperiodic particle ejection events are observed, but do not initiate pedestal collapse. Pedestal formation alters the spectrum of transport avalanches.

Mitigation of tokamak disruptions using high-pressure gas injection.

High-pressure gas-jet injection of neon and argon is shown to be a simple and robust method to mitigate the deleterious effects of disruptions on the DIII-D tokamak. The gas jet penetrates to the central plasma at its sonic velocity. The deposited species dissipates >95% of the plasma by radiation and substantially reduces mechanical stresses on the vessel caused by poloidal halo currents. The gas-jet species-charge distribution can include >50% fraction neutral species which inhibits runaway electrons. The favorable scaling of this technique to burning fusion plasmas is discussed.

Modified lattice Boltzmann method for compressible fluid simulations.

A modified lattice Boltzmann algorithm is shown to have much better stability to growing temperature perturbations, when compared with the standard lattice Boltzmann algorithm. The damping rates of long-wavelength waves, which determine stability, are derived using a collisional equilibrium distribution function which has the property that the Euler equations are obtained exactly in the limit of zero time step. Using this equilibrium distribution function, we show that our algorithm has inherent positive hyperviscosity and hyperdiffusivity, for very small values of viscosity and thermal diffusivity, which are lacking in the standard algorithm. Short-wavelength modes are shown to be stable for temperatures greater than a lower limit. Results from a computer code are used to compare these algorithms, and to confirm the damping rate predictions made analytically. Finite amplitude sound waves in the simulated fluid steepen, as expected from gas dynamic theory.