Search

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality. 2010 MSC: 00-01, 99-00, Kinetic theory, Vlasov equation, Magnetic fusion energy, Tokamak, Edge localized mode

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. In order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

We present an approach for constructing finite-volume methods for flux-divergence forms to any order of accuracy defined as the image of a smooth mapping from a rectangular discretization of an abstract coordinate space. Our approach is based on two ideas. The first is that of using higher-order quadrature rules to compute the flux averages over faces that generalize a method developed for Cartesian grids to the case of mapped grids. The second is a method for computing the averages of the metric terms on faces such that freestream preservation is automatically satisfied. We derive detailed formulas for the cases of fourth-order accurate discretizations of linear elliptic and hyperbolic partial differential equations. For the latter case, we combine the method so derived with Runge-Kutta time discretization and demonstrate how to incorporate a high-order accurate limiter with the goal of obtaining a method that is robust in the presence of discontinuities and underresolved gradients. For both elliptic and hyperbolic problems, we demonstrate that the resulting methods are fourth-order accurate for smooth solutions.

With the next generation of high-power laser facilities for inertial fusion coming online1,2, ensuring laser beam propagation through centimetre-scale plasmas is a key physics issue for reaching ignition. Existing experimental results3,4,5 including the most recent one6 are limited to small laser spots, low-interaction laser beam energies and small plasma volumes of 1–2 mm. Here, we demonstrate the propagation of an intense, high-energy, ignition-size laser beam through fusion-size plasmas on the National Ignition Facility (NIF) and find the experimental measurements to agree with full-scale modelling. Previous attempts to apply computer modelling as a predictive capability have been limited by the inherently multiscale description of the full laser–plasma interaction processes7,8,9,10,11. The findings of this study validate supercomputer modelling as an essential tool for the design of future ignition experiments.

The numerical simulation of plasmas is a critical tool for inertial confinement fusion (ICF). We have been working to improve the predictive capability of a continuum laser plasma interaction code pF3d, which couples a continuum hydrodynamic model of an unmagnetized plasma to paraxial wave equations modeling the laser light. Advanced numerical techniques such as local mesh refinement, multigrid, and multifluid Godunov methods have been adapted and applied to nonlinear heat conduction and to multifluid plasma models. We describe these algorithms and briefly demonstrate their capabilities.

We investigate the use of adaptive mesh refinement in the simulation of laser plasma filamentation. A numerical algorithm is constructed to solve model equations consisting of a fluid approximation of a quasineutral plasma combined with a paraxial light propagation model. The algorithm involves high-resolution plasma and light model discretizations on a block-structured, locally refined grid hierarchy, which is dynamically modified during the time integration to follow evolving fine-scale solution features. Comparisons of the efficiency of this approach to that of uniform grid calculations are presented.

A strategy for implementing source iteration on massively parallel computers for use in solving multigroup discrete ordinates neutron transport equations on three-dimensional Cartesian grids is proposed and analyzed. Based on an analysis of the memory requirement and floating-point complexity of the formal matrix-vector multiplication effected by a single source iteration, a data decomposition and communication strategy is presented that is designed to achieve good scalability with respect to all phase-space variables, i.e., neutron position, energy, and direction. A performance model is developed to analyze the scalability properties of the algorithm and to provide computational and heuristic strategies for determining a data decomposition that minimizes wall clock execution time. Numerical results are presented to demonstrate the performance of a specific implementation of this approach on a 1,024-node nCUBE/2.

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.

Summary form only given. Experiments on the National Ignition Facility (NIF) have employed the first four beams to measure propagation and laser backscattering losses in large ignition-size plasmas. Gas-filled targets between 2 mm and 7 mm length have been heated from one side by overlapping the focal spots of the four beams from one quad operated at 351 nm with a total intensity of 2times1015 W cm-2. The targets were filled with 1 atm of CO2 producing of up to 7 mm long homogeneously heated plasmas with densities of ne = 5times1020 cm-3 and temperatures of Te = 2 keV. The high energy in a NIF quad of beams of 16 kJ, illuminating the target from one direction, creates unique conditions for the study of laser plasma interactions at scale lengths not previously accessible. The propagation through the large-scale plasma was measured with a gated X-ray imager that was filtered for 3.5 keV X-rays. These data indicate that the beams interact with the full length of this ignition-scale plasma during the last ~1 ns of the experiment when applying full laser beam smoothing consisting of phase plates, smoothing by spectral dispersion and polarization smoothing. Measurements that only apply phase plates show, laser beam filamentation and reduced propagation speed. These results demonstrate the NIF experimental capabilities and further provide a benchmark for three-dimensional modeling of the laser-plasma interactions at ignition-size scale lengths

Geodesic acoustic modes (GAMs) are an important phenomenon in a tokamak edge plasma. They regulate turbulence in a low confinement (L-mode) regime and can play an important role in the low to high (L–H) mode transition. It is therefore of considerable importance to develop a detailed theoretical understanding of their dynamics and relaxation processes. The present work reports on the numerical modelling of collisionless GAM relaxation, including the effects of a strong radial electric field characteristic of a tokamak pedestal in a high confinement (H-mode) regime. The simulations demonstrate that the presence of a strong radial electric field enhances the GAM decay rate, and heuristic arguments elucidating this finding are provided. The numerical modelling is performed by making use of the continuum gyrokinetic code COGENT.

We are developing a new class of finite-volume methods on locally-refined and mapped grids, which are at least fourth-order accurate in regions where the solution is smooth. This paper discusses the implementation of such methods for time-dependent problems on both Cartesian and mapped grids with adaptive mesh refinement. We show 2D results with the Berger-Colella shock-ramp problem in Cartesian coordinates, and fourth-order accuracy of the solution of a Gaussian pulse problem in a polytropic gas in mapped coordinates.

We present edge gyrokinetic simulations of tokamak plasmas using the fully non-linear (full-f) continuum code TEMPEST. A non-linear Boltzmann model is used for the electrons. The electric field is obtained by solving the 2D gyrokinetic Poisson equation. We demonstrate the following. (1) High harmonic resonances (n > 2) significantly enhance geodesic-acoustic mode (GAM) damping at high q (tokamak safety factor), and are necessary to explain the damping observed in our TEMPEST q-scans and consistent with the experimental measurements of the scaling of the GAM amplitude with edge q 95 in the absence of obvious evidence that there is a strong q-dependence of the turbulent drive and damping of the GAM. (2) The kinetic GAM exists in the edge for steep density and temperature gradients in the form of outgoing waves, its radial scale is set by the ion temperature profile, and ion temperature inhomogeneity is necessary for GAM radial propagation. (3) The development of the neoclassical electric field evolves through different phases of relaxation, including GAMs, their radial propagation and their long-time collisional decay. (4) Natural consequences of orbits in the pedestal and scrape-off layer region in divertor geometry are substantial non-Maxwellian ion distributions and parallel flow characteristics qualitatively like those observed in experiments.

Experiments show that power is transferred between two copropagating 351nm laser beams crossing in an Al plasma when the frequency of the driven ion wave is shifted by a Mach 1 flow. The resonant amplification of a low-intensity (⩽2.5×1014W∕cm2) beam intersected by a high-intensity (7.0×1014W∕cm2) pump beam is determined by comparing the transmitted beam power to that measured in experiments where the plasma flow direction is reversed and the ion wave is evidently detuned. The polarization of the amplified light is also observed to align to the pump polarization consistent with ion-wave scattering. The amplification is found to reduce with probe-beam intensity demonstrating a nonlinear saturation mechanism that is effective when the ion-wave damping is weak, which is modeled with a calculation including both the nonlinear ion-wave frequency shifts due to ion trapping and whole-beam pump depletion.