34 Search Results

Abstract The magnetic separatrix surface is designed to provide the final and critical confinement to the hot stationary-operation core plasma in modern tokamak reactors in the absence of an external magnetic perturbation or transient MHD perturbation, while diverting the exhaust heat to divertor plates. All the stationary operational boundary plasma studies and reactor designs have been performed under this assumption. However, there has been a long-standing suspicion that a stationary-operation tokamak plasma even without an external magnetic perturbations (MPs) or edge localized modes (ELMs) activities may not have a stable closed separatrix surface, especially near the magnetic X-point. Here, themore » first gyrokinetic numerical observation is reported that the divertor separatrix surface, due to homoclinic tangles caused by intrinsic electromagnetic turbulence, is not a stable closed surface in a stationary operation phase even without MPs or ELMs. Unlike the MP- or ELM-driven homoclinic tangles that could cause deleterious effects to core confinement or divertor plates, it is found that the micro-turbulence driven homoclinic tangles could connect the divertor plasma to the pedestal plasma in a constructive way by broadening the divertor heat-exhaust footprint and weakening the pedestal slope to the ELM-safe direction. Micro-turbulent homoclinic tangles can open a new research direction in understanding and controlling these two most troublesome and non-locally connected edge-plasma issues in a tokamak fusion reactor.« less

Here, the application of a bundling technique to model the diverse charge states of tungsten impurity species in total-f gyrokinetic simulations is demonstrated. The gyrokinetic bundling method strategically groups tungsten ions of similar charge, optimizing computational efficiency. The initial radial configuration of these bundles and their respective charges are derived from a coronal approximation and the quasi-neutrality of the plasma. A low-density JET H-mode like plasma is simulated using the neoclassical version of XGC across the entire plasma volume, spanning from the magnetic axis to the divertor. An accumulation of tungsten is observed at the pedestal top, as a resultmore » of low-Z tungsten ions moving inward from the scrape-off-layer into the core region and high-Z tungsten ions moving outward from the core into the pedestal. This organization of the fluxes cannot be captured by a single tungsten-ion simulation. Large up-down poloidal asymmetries of tungsten form in the pedestal and strongly influence the direction of neoclassical fluxes. The temperature screening effect and its correlation with asymmetries are analyzed.« less

Abstract Ion orbit loss has been used to model the formation of a strong negative radial electric field E _r in the tokamak edge, as well as edge momentum transport and toroidal rotation. To quantitatively measure ion orbit loss, an orbit-flux formulation has been developed and numerically applied to the gyrokinetic particle-in-cell code XGC. We study collisional ion orbit loss in an axisymmetric DIII-D L-mode plasma using gyrokinetic ions and drift-kinetic electrons. Numerical simulations, where the plasma density and temperature profiles are maintained through neutral ionization and heating, show the formation of a quasisteady negative E _rmore » in the edge. We have measured a radially outgoing ion gyrocenter flux due to collisional scattering of ions into the loss orbits, which is balanced by the radially incoming ion gyrocenter flux from confined orbits on the collisional time scale. This suggests that collisional ion orbit loss can shift E _r in the negative direction compared to that in plasmas without orbit loss. It is also found that collisional ion orbit loss can contribute to a radially outgoing (counter-current) toroidal-angular-momentum flux, which is not balanced by the toroidal-angular-momentum flux carried by ions on the confined orbits. Therefore, the edge toroidal rotation shifts in the co-current direction on the collisional time scale.« less

We report the simplified δf mixed-variable/pullback electromagnetic simulation algorithm implemented in XGC for core plasma simulations by Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total- f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence, and neoclassical physics. The δf mixed-variable/pullback electromagnetic implementation is based on the pioneering work by Kleiber and Mischenko et al. [Kleiber et al., Phys. Plasmas 23, 032501 (2016); Mishchenko et al., Comput. Phys. Commun. 238, 194 (2019)]. An electromagnetic demonstration simulation is performed in a DIII-D-like, H-modemore » boundary plasma, including a corresponding comparative electrostatic simulation, which confirms that the electromagnetic simulation is necessary for a higher fidelity understanding of the electron particle and heat transport even at the low- β pedestal foot in the vicinity of the magnetic separatrix.« less

Isotope effects under the influence of a radial electric field are examined in a helical magnetic field configuration. We perform global gyrokinetic simulations with additional poloidal rotations to estimate quasi-linear heat flux due to ion temperature gradient mode under the mixing length model. In single-ion-species plasmas, the mass number dependency of heat flux agrees with gyro-Bohm scaling in the absence of a radial electric field. Favorable mass number dependencies violating gyro-Bohm scaling are observed in the presence of a global radial electric field or a heavy hydrogen component in multi-ion-species plasmas. The radial electric field and the heavy hydrogen componentmore » affect the heat flux through an increase of wavelength as well as mode stabilization. Poloidal Mach number characterizes the transition from unfavorable to favorable mass number dependency under radial electric fields. While the heat flux is independent of mass number for a given poloidal Mach number, the heat flux decreases for higher mass numbers in a given radial electric field. The heat flux is also independent of average mass number in multi-ion-species plasmas because the heavy hydrogen component effectively enhances the light hydrogen heat flux. The present results are potentially relevant to the violation of gyro-Bohm scaling observed in the recent deuterium experiments in the Large Helical Device.« less

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effectsmore » of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high-β are reduced to the level of that of collisionless trapped electron modes.« less

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirablymore » large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.« less

A pellet source model has been implemented into the gyrokinetic code XGC and applied to neoclassical geodesic acoustic mode (GAM) study as the first step to investigate potentially rapid kinetic spread of pellet-born particles following their injection into plasma. Here in this study, GAM oscillations of the radial electric field, as well as of their effect on the radial particle and energy fluxes, are studied and the observed frequencies agree with theoretical expectations. It is verified that (i) GAM oscillations are driven both within the pellet-fuelled region and for locations at a larger radius than this region and not formore » locations radially inside the pellet-fuelled region, and that, (ii) even though the poloidal spread along the magnetic field lines is fast, the time-averaged radial plasma transport under GAM oscillations is kept to roughly the same level as the neoclassical transport in the absence of plasma turbulence.« less

A fully imore » mplicit particle-in-cell method for handling the $v_{\parallel }$-formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the $v_{\parallel }$-formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the $p_{\parallel }$-formalism. The $v_{\parallel }$-formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.« less

An encoder–decoder neural network has been used to examine the possibility for acceleration of a partial integro-differential equation, the Fokker–Planck–Landau collision operator. This is part of the governing equation in the massively parallel particle-in-cell code XGC, which is used to study turbulence in fusion energy devices. The neural network emphasizes physics-inspired learning, where it is taught to respect physical conservation constraints of the collision operator by including them in the training loss, along with the $$\ell _2$$ loss. In particular, network architectures used for the computer vision task of semantic segmentation have been used for training. A penalization method ismore » used to enforce the ‘soft’ constraints of the system and integrate error in the conservation properties into the loss function. During training, quantities representing the particle density, momentum and energy for all species of the system are calculated at each configuration vertex, mirroring the procedure in XGC. This simple training has produced a median relative loss, across configuration space, of the order of $$10^{-4}$$ , which is low enough if the error is of random nature, but not if it is of drift nature in time steps. The run time for the current Picard iterative solver of the operator is $O(n^2)$$ , where $$n$$ is the number of plasma species. As the XGC1 code begins to attack problems including a larger number of species, the collision operator will become expensive computationally, making the neural network solver even more important, especially since its training only scales as $$O(n)$ . A wide enough range of collisionality has been considered in the training data to ensure the full domain of collision physics is captured. An advanced technique to decrease the losses further will be subject of a subsequent report. Eventual work will include expansion of the network to include multiple plasma species.« less