Parallel filtering in global gyrokinetic simulations

In this work, a Fourier solver [B.F. McMillan, S. Jolliet, A. Bottino, P. Angelino, T.M. Tran, L. Villard, Comp. Phys. Commun. 181 (2010) 715] is implemented in the global Eulerian gyrokinetic code GT5D [Y. Idomura, H. Urano, N. Aiba, S. Tokuda, Nucl. Fusion 49 (2009) 065029] and in the global Particle-In-Cell code ORB5 [S. Jolliet, A. Bottino, P. Angelino, R. Hatzky, T.M. Tran, B.F. McMillan, O. Sauter, K. Appert, Y. Idomura, L. Villard, Comp. Phys. Commun. 177 (2007) 409] in order to reduce the memory of the matrix associated with the field equation. This scheme is verified with linear and nonlinear simulations of turbulence. It is demonstrated that the straight-field-line angle is the coordinate that optimizes the Fourier solver, that both linear and nonlinear turbulent states are unaffected by the parallel filtering, and that the k_∥ spectrum is independent of plasma size at fixed normalized poloidal wave number.     

Dimits shift in realistic gyrokinetic plasma-turbulence simulations

In simulations of turbulent plasma transport due to long wavelength (k perpendicular rhoi < or = 1) electrostatic drift-type instabilities, we find a persistent nonlinear up-shift of the effective threshold. Next-generation tokamaks will likely benefit from the higher effective threshold for turbulent transport, and transport models should incorporate suitable corrections to linear thresholds. The gyrokinetic simulations reported here are more realistic than previous reports of a Dimits shift because they include nonadiabatic electron dynamics, strong collisional damping of zonal flows, and finite electron and ion collisionality together with realistic shaped magnetic geometry. Reversing previously reported results based on idealized adiabatic electrons, we find that increasing collisionality reduces the heat flux because collisionality reduces the nonadiabatic electron microinstability drive.     

Magnetic stochasticity in gyrokinetic simulations of plasma microturbulence

Analysis of the magnetic field structure from electromagnetic simulations of tokamak ion temperature gradient turbulence demonstrates that the magnetic field can be stochastic even at very low plasma pressure. The degree of magnetic stochasticity is quantified by evaluating the magnetic diffusion coefficient. We find that the magnetic stochasticity fails to produce a dramatic increase in the electron heat conductivity because the magnetic diffusion coefficient remains small.     

Unlike-particle collision operator for gyrokinetic particle simulations

Plasmas in modern tokamak experiments contain a significant fraction of impurity ion species in addition to main deuterium background. A new unlike-particle collision operator for δf particle simulation has been developed to self-consistently study the non-local effects of impurities on neoclassical transport in toroidal plasmas. A new algorithm for simulation of cross-collisions between different ion species includes test-particle and conserving field-particle operators. The field-particle operator is designed to enforce conservation of number, momentum and energy. It was shown that the new operator correctly simulates the thermal equilibration of different plasma components. It was verified that the ambipolar radial electric field reaches steady state when the total radial guiding center particle current vanishes.     

Particle dynamics simulations of triboelectric charging in granular insulator systems

Particle dynamics simulations are carried out to study triboelectric charging in granular systems composed of a single insulating material. The simulations implement a model in which electrons trapped in localized high energy states can be transferred during collisions to low energy states in the other particle. It is shown that this effect alone can generate electrostatic charging in the system, and cause net electron transfer from larger particles to smaller particles. The magnitude of charging is small for systems of a single particle size but becomes much greater for a system with polydispersal particle sizes, due to the net electron transfer from larger to smaller particles. The negative charge of smaller particles, and positive charge of larger particles has been observed in field studies and laboratory experiments of granular systems.     

Collisional gyrokinetic full-f particle-in-cell simulations on open field lines with PICLS

Applying gyrokinetic simulations in theoretical turbulence and transport studies for the plasma edge and scrape-off layer (SOL) presents significant challenges. To particularly account for steep density and temperature gradients in the SOL, the “full-f” code PICLS was developed. PICLS is a gyrokinetic particle-in-cell (PIC) code, is based on an electrostatic model with a linearized field equation, and uses kinetic electrons. In previously published results, we applied PICLS to the well-studied 1D parallel transport problem during an edge-localized mode (ELM) in the SOL without collisions. As an extension to this collision-less case and in preparation for 3D simulations, in this work, a collisional model will be introduced. The implemented Lenard–Bernstein collision operator and its Langevin discretization will be shown. Conservation properties of the collision operator, as well as a comparison of the collisional and non-collisional case, will be discussed.     

Nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence for the optimized Wendelstein 7-X stellarator

Ion-temperature-gradient turbulence constitutes a possibly dominant transport mechanism for optimized stellarators, in view of the effective suppression of neoclassical losses characterizing these devices. Nonlinear gyrokinetic simulation results for the Wendelstein 7-X stellarator [G. Grieger, in (IAEA, Vienna, 1991) Vol. 3, p. 525]-assuming an adiabatic electron response-are presented. Several fundamental features are discussed, including the role of zonal flows for turbulence saturation, the resulting flux-gradient relationship, and the coexistence of ion-temperature-gradient modes with trapped ion modes in the saturated state.     

Effects of anisotropic thermal conductivity in magnetohydrodynamics simulations of a reversed-field pinch

A compressible magnetohydrodynamics simulation of the reversed-field pinch is performed including anisotropic thermal conductivity. When the thermal conductivity is much larger in the direction parallel to the magnetic field than in the perpendicular direction, magnetic field lines become isothermal. As a consequence, as long as magnetic surfaces exist, a temperature distribution is observed displaying a hotter confined region, while an almost uniform temperature is produced when the magnetic field lines become chaotic. To include this effect in the numerical simulation, we use a multiple-time-scale analysis, which allows us to reproduce the effect of a large parallel thermal conductivity. The resulting temperature distribution is related to the existence of closed magnetic surfaces, as observed in experiments. The magnetic field is also affected by the presence of an anisotropic thermal conductivity.     

Particle simulations of space weather

We review the application of particle simulation techniques to the full kinetic study of space weather events. We focus especially on the methods designed to overcome the difficulties created by the tremendous range of time and space scales present in the physical systems. We review the aspects of the derivation of the particle in cell (PIC) method relevant to the discussion. We consider first the explicit formulation highlighting its severe limitations due to the presence of stability constraints. Next we introduce implicit methods designed to remove such constraints. We describe both fully implicit methods based on the use of non-linear iteration solvers and semi-implicit methods based on the linearization of the coupling and on simpler linear solvers. We focus the discussion on the implicit moment method but remark its differences from the direct implicit method. The application of adaptive methods within PIC is discussed. Finally practical considerations about the implementation of the implicit PIC method on massively parallel computers to conduct studies of space weather events are given.     

Particle pinch with fully noninductive lower hybrid current drive in Tore Supra

Recently, plasmas exceeding 4 min have been obtained with lower hybrid current drive (LHCD) in Tore Supra. These LHCD plasmas extend for over 80 times the resistive current diffusion time with zero loop voltage. Under such unique conditions the neoclassical particle pinch driven by the toroidal electric field vanishes. Nevertheless, the density profile remains peaked for more than 4 min. For the first time, the existence of an inward particle pinch in steady-state plasma without toroidal electric field, much larger than the value predicted by the collisional neoclassical theory, is experimentally demonstrated.     

Particle simulations of radio-frequency glow discharges

Particle-in-cell (PIC) simulations are used to study the structure of radio-frequency (RF) glow discharges in helium between parallel-plate electrodes. The authors have examined a range of conditions and report on a variety of observed phenomena. Comparisons to experiment and analytical models are made, when possible. The differences between discharges in which secondary electrons play a key role in sustaining the discharge and those in which secondary electrons are unimportant are examined in three cases which illustrate the importance of the discharge-sustaining mechanisms. Electron-energy distributions are found to be, in general, non-Maxwellian, with shapes that depend in complex ways on discharge conditions. In the absence of secondary electron emission, electron heating in the sheath regions of the discharge is enhanced at higher voltages compared to ohmic heating in the bulk of the plasma. Fast electrons accelerated by the advancing sheath can carry a substantial fraction of the conduction current in the bulk of the discharge, reducing the effective bulk ohmic heating of electrons. Ion-energy distributions at electrode surfaces have been predicted and are compared to experimental measurements     

Particle imaging sizing: GLMT simulations

Defocused off-axis images of particles are computed in the framework of the generalized Lorenz-Mie theory. Two exemplifying cases are studied: interferometric sizing at large off-axis angles and imaging in near forward directions (shadow Doppler velocimetry)     

Free energy cascade in gyrokinetic turbulence

In gyrokinetic theory, the quadratic nonlinearity is known to play an important role in the dynamics by redistributing (in a conservative fashion) the free energy between the various active scales. In the present study, the free energy transfer is analyzed for the case of ion temperature gradient driven turbulence. It is shown that it shares many properties with the energy transfer in fluid turbulence. In particular, one finds a (strongly) local, forward (from large to small scales) cascade of free energy in the plane perpendicular to the background magnetic field. These findings shed light on some fundamental properties of plasma turbulence, and encourage the development of large-eddy-simulation techniques for gyrokinetics.     

The gyrokinetic water-bag modeling in toroidal geometry

The present paper addresses the gyrokinetic water-bag model in toroidal geometry. The previous works were focused on the water-bag concept in magnetized cylindrical plasmas. Here we report on the possibility to improve the water-bag model by taking into account the curvature and gradient drifts. After a presentation of the model, a local linear analysis with some approximations is performed. Interchange and ion temperature gradient instabilities are examined with this new gyro-water-bag model in order to show its ability and its theoretical interest in describing kinetic instabilities in toroidal geometry.     

Density fluctuations in lattice-Boltzmann simulations of multiphase fluids in a closed system

A two-dimensional single component two-phase lattice Boltzmann model was used to simulate the Rayleigh–Taylor instability in a closed system. Spatiotemporally variable densities were generated by gravity acting on the fluid density. The density fluctuations were triggered by rapid changes in the fluid velocity induced by changes in the interface geometry and impact of the dense fluid on the rigid lower boundary of the computational domain. The ratio of the maximum density fluctuations to the maximum fluid velocity increased more rapidly at low velocities than at high velocities. The ratio of the maximum density fluctuations in the dense phase to its maximum velocity was on the order of the inverse of the sound speed. The solution became unstable when the density-based maximum local Knudsen number exceeded 0.13.     

Particle subgrid scale modelling in large-eddy simulations of particle-laden turbulence

This study is concerned with particle subgrid scale (SGS) modelling in large-eddy simulations (LESs) of particle-laden turbulence. Although many particle-laden LES studies have neglected the effect of the SGS on the particles, several particle SGS models have been proposed in the literature. In this research, the approximate deconvolution method (ADM) and the stochastic models of Fukagata et al. (Dynamics of Brownian particles in a turbulent channel flow, Heat Mass Transf. 40 (2004), 715–726) Shotorban and Mashayek (A stochastic model for particle motion in large-eddy simulation, J. Turbul. 7 (2006), 1–13) and Berrouk et al. (Stochastic modelling of inertial particle dispersion by subgrid motion for LES of high Reynolds number pipe flow, J. Turbul. 8 (2007), pp. 1–20) are analysed. The particle SGS models are assessed using both a priori and a posteriori simulations of inertial particles in a periodic box of decaying, homogeneous and isotropic turbulence with an initial Reynolds number of Re_λ = 74. The model results are compared with particle statistics from a direct numerical simulation (DNS). Particles with a large range of Stokes numbers are tested using various filter sizes and stochastic model constant values. Simulations with and without gravity are performed to evaluate the ability of the models to account for the crossing trajectory and continuity effects. The results show that ADM improves results but is only capable of recovering a portion of the SGS turbulent kinetic energy. Conversely, the stochastic models are able to recover sufficient SGS energy, but show a large range of results dependent on the Stokes number and filter size. The stochastic models generally perform best at small Stokes numbers, but are unable to predict preferential concentration.     

Dynamic weighting in simulations of spin systems

We establish the invariance property of a dynamically weighted Monte Carlo process and apply the method to the simulation of spin glasses and Ising models. For the two-dimensional Edwards-Anderson model, we obtain an ergodicity time of O(L^2.44), where L is the linear dimension of the lattice. The results suggest that dynamic weighting is a highly promising new tool for Monte Carlo simulation.