Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userVerification of long wavelength electromagnetic modes with a gyrokinetic-fluid hybrid model in the XGC code
110
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
As an alternative option to kinetic electrons, the gyrokinetic total-f particle-in-cell (PIC) code XGC1 has been extended to the MHD/fluid type electromagnetic regime by combining gyrokinetic PIC ions with massless drift-fluid electrons analogous to Chen and Parker, Physics of Plasmas 8, 441 (2001). Two representative long wavelength modes, shear Alfven waves and resistive tearing modes, are verified in cylindrical and toroidal magnetic field geometries.
rate research
Read More
A fully implicit particle-in-cell method for handling the $v_parallel$-formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the $v_parallel$-formalism, we avoid introducing the non-physical skin terms in Amp`{e}res 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 $delta 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{e}n wave propagation in addition to a case demonstrating the ITG-KBM transition. The ITG-KBM transition case is compared to results obtained from other $delta f$ gyrokinetic codes/schemes, whose verification has already been archived in the literature.
This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a ``hybrid model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.
The linear properties of the geodesic acoustic modes (GAM) in tokamaks are investigated by means of the comparison of analytical theory and gyrokinetic numerical simulations. The dependence on the value of the safety factor, finite-orbit-width of the ions in relation to the radial mode width, magnetic-flux-surface shaping, and electron/ion mass ratio are considered. Nonuniformities in the plasma profiles (such as density, temperature, or safety factor), electro-magnetic effects, collisions and presence of minority species are neglected. Also, only linear simulations are considered, focusing on the local dynamics. We use three different gyrokinetic codes: the lagrangian (particle-in-cell) code ORB5, the eulerian code GENE and semi-lagrangian code GYSELA. One of the main aims of this paper is to provide a detailed comparison of the numerical results and analytical theory, in the regimes where this is possible. This helps understanding better the behavior of the linear GAM dynamics in these different regimes, the behavior of the codes, which is crucial in the view of a future work where more physics is present, and the regimes of validity of each specific analytical dispersion relation.
We are developing a new continuum gyrokinetic code, Gkeyll, for use in edge plasma simulations, and here present initial simulations of turbulence on open field lines with model sheath boundary conditions. The code implements an energy conserving discontinuous Galerkin scheme, applicable to a general class of Hamiltonian equations. Several applications to test problems have been done, including a calculation of the parallel heat-flux on divertor plates resulting from an ELM crash in JET, for a 1x/1v SOL scenario explored previously, where the ELM is modeled as a time-dependent intense upstream source. Here we present initial simulations of turbulence on open field lines in the LAPD linear plasma device. We have also done simulations in a helical open-field-line geometry. While various simplifications have been made at present, this still includes some of the key physics of SOL turbulence, such as bad-curvature drive for instabilities and rapid parallel losses with sheath boundary conditions. This is useful for demonstrating the overall feasibility of this approach and for initial physics studies of SOL turbulence. We developed a novel version of DG that uses Maxwellian-weighted basis functions while still preserving exact particle and energy conservation. The Maxwellian-weighted DG method achieves the same error with 4 times less computational cost in 1v, or 16 times lower cost in the 2 velocity dimensions of gyrokinetics (assuming memory bandwidth is the limiting factor).
Vlasiator is a new hybrid-Vlasov plasma simulation code aimed at simulating the entire magnetosphere of the Earth. The code treats ions (protons) kinetically through Vlasovs equation in the six-dimensional phase space while electrons are a massless charge-neutralizing fluid [M. Palmroth et al., Journal of Atmospheric and Solar-Terrestrial Physics 99, 41 (2013); A. Sandroos et al., Parallel Computing 39, 306 (2013)]. For first global simulations of the magnetosphere, it is critical to verify and validate the model by established methods. Here, as part of the verification of Vlasiator, we characterize the low-$beta$ plasma wave modes described by this model and compare with the solution computed by the Waves in Homogeneous, Anisotropic Multicomponent Plasmas (WHAMP) code [K. Ronnmark, Kiruna Geophysical Institute Reports 179 (1982)], using dispersion curves and surfaces produced with both programs. The match between the two fundamentally different approaches is excellent in the low-frequency, long wavelength range which is of interest in global magnetospheric simulations. The left-hand and right-hand polarized wave modes as well as the Bernstein modes in the Vlasiator simulations agree well with the WHAMP solutions. Vlasiator allows a direct investigation of the importance of the Hall term by including it in or excluding it from Ohms law in simulations. This is illustrated showing examples of waves obtained using the ideal Ohms law and Ohms law including the Hall term. Our analysis emphasizes the role of the Hall term in Ohms law in obtaining wave modes departing from ideal magnetohydrodynamics in the hybrid-Vlasov model.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
