Search

We present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

An analysis of the effects of ion trapping on ion acoustic waves excited by the stimulated Brillouin scattering of crossing intense laser beams is presented. Ion trapping alters the dispersion of ion acoustic waves by nonlinearly shifting the normal mode frequency and by reducing the ion Landau damping. This in turn can influence the energy transfer between two crossing laser beams in the presence of plasma flows such that stimulated Brillouin scattering (SBS) occurs. The same ion trapping physics can influence the saturation of SBS in other circumstances. A one-dimensional analytical model is presented along with reasonably successful comparisons of the theory to results from particle simulations and laboratory experiments. An analysis of the vulnerability of the National Ignition Facility Inertial Confinement Fusion point design [S. W. Haan et al., Fusion Sci. Technol. 41, 164 (2002)] is also presented.

The first 4D (axisymmetric) high-order continuum gyrokinetic transport simulations that span the magnetic separatrix of a tokamak are presented. The modeling is performed with the COGENT code, which is distinguished by fourth-order finite-volume discretization combined with mapped multiblock grid technology to handle the strong anisotropy of plasma transport and the complex X-point divertor geometry with high accuracy. The calculations take into account the effects of fully nonlinear Fokker-Plank collisions, electrostatic potential variations, and anomalous radial transport. Topics discussed include: (a) ion orbit loss and the associated toroidal rotation and (b) edge plasma relaxation in the presence of anomalous radial transport.

The development of the continuum gyrokinetic code COGENT for edge plasma simulations is reported. The present version of the code models a nonlinear axisymmetric 4D (R, v∥, μ) gyrokinetic equation coupled to the long-wavelength limit of the gyro-Poisson equation. Here, R is the particle gyrocenter coordinate in the poloidal plane, and v∥ and μ are the guiding center velocity parallel to the magnetic field and the magnetic moment, respectively. The COGENT code utilizes a fourth-order finite-volume (conservative) discretization combined with arbitrary mapped multiblock grid technology (nearly field-aligned on blocks) to handle the complexity of tokamak divertor geometry with high accuracy. Topics presented are the implementation of increasingly detailed model collision operators, and the results of neoclassical transport simulations including the effects of a strong radial electric field characteristic of a tokamak pedestal under H-mode conditions.

Two-dimensional Vlasov simulations of nonlinear electron plasma waves are presented, in which the interplay of linear and nonlinear kinetic effects is evident. The plasma wave is created with an external traveling wave potential with a transverse envelope of width Δy such that thermal electrons transit the wave in a “sideloss” time, tsl~Δy/ve. Here, ve is the electron thermal velocity. The quasisteady distribution of trapped electrons and its self-consistent plasma wave are studied after the external field is turned off. In cases of particular interest, the bounce frequency, ωbe=keϕ/me, satisfies the trapping condition ωbetsl>2π such that the wave frequency is nonlinearly downshifted by an amount proportional to the number of trapped electrons. Here, k is the wavenumber of the plasma wave and ϕ is its electric potential. For sufficiently short times, the magnitude of the negative frequency shift is a local function of ϕ. Because the trapping frequency shift is negative, the phase of the wave on axis lags ...