The predictions of gyrokinetic and gyrofluid simulations of ion-temperature-gradient (ITG) instability and turbulence in tokamak plasmas as well as some tokamak plasma thermal transport models, which have been widely used for predicting the performance of the proposed ITER tokamak, are compared. These comparisons provide information on effects of differences in the physics content of the various models and on the fusion-relevant figures of merit of plasma performance predicted by the models. Many of the comparisons are undertaken for a simplified plasma model and geometry which is an idealization of the plasma conditions and geometry in a DIII-D H-mode experiment. Most of the models show good agreements in their predictions and assumptions for the linear growth rates and frequencies. There are some differences associated with different equilibria. However, there are significant differences in the transport levels between the models. The causes of some of the differences are examined in some detail, with particular attention to numerical convergence in the turbulence simulations (with respect to simulation mesh size, system size and, for particle-based simulations, the particle number). The implications for predictions of fusion plasma performance are also discussed.

The turbulent Rayleigh-Taylor instability is investigated in the limit of strong mode-coupling using a variety of high-resolution, multimode, three dimensional numerical simulations ͑NS͒. The perturbations are initialized with only short wavelength modes so that the self-similar evolution ͑i.e., bubble diameter D b ϰamplitude h b) occurs solely by the nonlinear coupling ͑merger͒ of saturated modes. After an initial transient, it is found that h b ϳ␣ b Agt 2 , where AϭAtwood number, gϭacceleration, and tϭtime. The NS yield D b ϳh b /3 in agreement with experiment but the simulation value ␣ b ϳ0.025Ϯ0.003 is smaller than the experimental value ␣ b ϳ0.057Ϯ0.008. By analyzing the dominant bubbles, it is found that the small value of ␣ b can be attributed to a density dilution due to fine-scale mixing in our NS without interface reconstruction ͑IR͒ or an equivalent entrainment in our NS with IR. This may be characteristic of the mode coupling limit studied here and the associated ␣ b may represent a lower bound that is insensitive to the initial amplitude. Larger values of ␣ b can be obtained in the presence of additional long wavelength perturbations and this may be more characteristic of experiments. Here, the simulation data are also analyzed in terms of bubble dynamics, energy balance and the density fluctuation spectra.

No abstract

Recent gyrokinetic simulations of electron temperature gradient (ETG) turbulence with flux-tube continuum codes vs. the global particle-in-cell (PIC) code GTC yielded different results despite similar plasma parameters. Differences between the simulations results were attributed to insufficient phase-space resolution and novel physics associated with toroidicity and/or global simulations. We have reproduced the results of the global PIC code using the flux-tube PIC code PG3EQ, thereby eliminating global effects as the cause of the discrepancy. We show that the late-time decay of ETG turbulence and the steady-state heat transport observed in these PIC simulations results from discrete particle noise. Discrete particle noise is a numerical artifact, so both these PG3EQ simulations and the previous GTC simulations have nothing to say about steady-state ETG turbulence and the associated anomalous heat transport. In the course of this work we develop three diagnostics which can help to determine if a particular PIC simulation has become dominated by discrete particle noise.-2-

Numerical simulations of electron temperature gradient (ETG) turbulence are presented which characterize the ETG fluctuation spectrum, establish limits to the validity of the adiabatic ion model often employed in studying ETG turbulence, and support the tentative conclusion that plasmaoperating regimes exist in which ETG turbulence produces sufficient electron heat transport to be experimentally relevant. We resolve prior

In this paper, particle simulation algorithms with time-varying weights for the gyrokinetic Vlasov-Poisson system have been developed. The primary purpose is to use them for the removal of the selected nonlinearities in the simulation of gradient-driven microturbulence so that the , relative importance of the various nonlinear effects can be assessed. It is hoped that the use of these procedures will result in a better understanding of the transport mechanisms and scaling in tokamaks. Another application of these algorithms is for the improvement of the numerical properties of the simulation plasma. For instance, implementations of such algorithms 1) enable us to suppress the intrinsic numerical noise in the simulation, and 2) also make it possible to regulate the weights of the fast-moving particles and, in mm, to eliminate the associated high frequency osciUations. Examples of their application to drift-type instabilities in slab geometry are given. We note that the work reported here represents the first successful use of the weighted algorithms in particle codes for the nonlinear simulation of plasmas.

In diverse areas of science and technology, including inertial confinement fusion (ICF), astrophysics, geophysics, and engineering processes, turbulent mixing induced by hydrodynamic instabilities is of scientific interest as well as practical significance. Because of the fundamental roles they often play in ICF and other applications, three classes of hydrodynamic instability-induced turbulent flows—those arising from the Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities—have attracted much attention. ICF implosions, supernova explosions, and other applications illustrate that these phases of instability growth do not occur in isolation, but instead are connected so that growth in one phase feeds through to initiate growth in a later phase. Essentially, a description of these flows must encompass both the temporal and spatial evolution of the flows from their inception. Hydrodynamic instability will usually start from potentially infinitesimal spatial perturbations, will eventually transition to a turbulent flow, and then will reach a final state of a true multiscale problem. Indeed, this change in the spatial scales can be vast, with hydrodynamic instability evolving from just a few microns to thousands of kilometers in geophysical or astrophysical problems. These instabilities will evolve through different stages before transitioning to turbulence, experiencing linear, weakly, and highly nonlinear states. The challenges confronted by researchers are enormous. The inherent difficulties include characterizing the initial conditions of such flows and accurately predicting the transitional flows. Of course, fully developed turbulence, a focus of many studies because of its major impact on the mixing process, is a notoriously difficult problem in its own right. In this pedagogical review, we will survey challenges and progress, and also discuss outstanding issues and future directions.

Three-dimensional high-resolution simulations (up to 8 billion zones) have been performed for a Richtmyer-Meshkov instability produced by passing a shock through a contact discontinuity with a two-scale initial perturbation. The setup approximates shock tube experiments with a membrane pushed through a wire mesh. The simulation produces mixing-layer widths similar to those observed experimentally. Comparison of runs at various resolutions suggests a transition from unstable to turbulent flow as the numerical Reynolds number is increased. At the highest resolutions, the spectrum exhibits a region of power-law decay, in which the spectral flux is approximately constant, suggestive of an inertial range, but with a weaker wavenumber dependence than Kolmogorov scaling, about k −5/4. Analysis of structure functions at the end of the simulation indicates the persistence of structures with velocities largest in the stream-wise direction. Comparison of three-dimensional and two-dimensional runs illustrates the tendency toward forward cascade in