ترغب بنشر مسار تعليمي؟ اضغط هنا

Resolving velocity space dynamics in continuum gyrokinetics

57   0   0.0 ( 0 )
 نشر من قبل Michael Barnes
 تاريخ النشر 2009
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Many plasmas of interest to the astrophysical and fusion communities are weakly collisional. In such plasmas, small scales can develop in the distribution of particle velocities, potentially affecting observable quantities such as turbulent fluxes. Consequently, it is necessary to monitor velocity space resolution in gyrokinetic simulations. In this paper, we present a set of computationally efficient diagnostics for measuring velocity space resolution in gyrokinetic simulations and apply them to a range of plasma physics phenomena using the continuum gyrokinetic code GS2. For the cases considered here, it is found that the use of a collisionality at or below experimental values allows for the resolution of plasma dynamics with relatively few velocity space grid points. Additionally, we describe implementation of an adaptive collision frequency which can be used to improve velocity space resolution in the collisionless regime, where results are expected to be independent of collision frequency.

قيم البحث

اقرأ أيضاً

First-principles simulations of tokamak turbulence have proven to be of great value in recent decades. We develop a pseudo-spectral velocity formulation of the turbulence equations that smoothly interpolates between the highly efficient but lower res olution 3D gyrofluid representation and the conventional but more expensive 5D gyrokinetic representation. Our formulation is a projection of the nonlinear gyrokinetic equation onto a Laguerre-Hermite velocity-space basis. We discuss issues related to collisions, closures, and entropy. While any collision operator can be used in the formulation, we highlight a model operator that has a particularly sparse Laguerre-Hermite representation, while satisfying conservation laws and the H theorem. Free streaming, magnetic drifts, and nonlinear phase mixing each give rise to closure problems, which we discuss in relation to the instabilities of interest and to free energy conservation. We show that the model is capable of reproducing gyrokinetic results for linear instabilities and zonal flow dynamics. Thus the final model is appropriate for the study of instabilities, turbulence, and transport in a wide range of geometries, including tokamaks and stellarators.
Plasma turbulence is studied via direct numerical simulations in a two-dimensional spatial geometry. Using a hybrid Vlasov-Maxwell model, we investigate the possibility of a velocity-space cascade. A novel theory of space plasma turbulence has been r ecently proposed by Servidio {it et al.} [PRL, {bf 119}, 205101 (2017)], supported by a three-dimensional Hermite decomposition applied to spacecraft measurements, showing that velocity space fluctuations of the ion velocity distribution follow a broad-band, power-law Hermite spectrum $P(m)$, where $m$ is the Hermite index. We numerically explore these mechanisms in a more magnetized regime. We find that (1) the plasma reveals spectral anisotropy in velocity space, due to the presence of an external magnetic field (analogous to spatial anisotropy of fluid and plasma turbulence); (2) the distribution of energy follows the prediction $P(m)sim m^{-2}$, proposed in the above theoretical-observational work; and (3) the velocity-space activity is intermittent in space, being enhanced close to coherent structures such as the reconnecting current sheets produced by turbulence. These results may be relevant to the nonlinear dynamics weakly-collisional plasma in a wide variety of circumstances.
We present an energy-conserving discontinuous Galerkin scheme for the full-$f$ electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for $partial A_parallel/partial t$, the inductive com ponent of the parallel electric field, using a generalized Ohms law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Amp`ere cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.
The Fermi acceleration model was introduced to describe how cosmic ray particles are accelerated to great speeds by interacting with moving magnetic fields. We identify a new variation of the model where light ions interact with a moving wall while u ndergoing pitch angle scattering through Coulomb collisions due to the presence of a heavier ionic species. The collisions introduce a stochastic component which adds complexity to the particle acceleration profile and sets it apart from collisionless Fermi acceleration models. The unusual effect captured by this simplified variation of Fermi acceleration is the non-conservation of phase space, with the possibility for a distribution of particles initially monotonically decreasing in energy to exhibit an energy peak upon compression. A peaked energy distribution might have interesting applications, such as to optimize fusion reactivity or to characterize astrophysical phenomena that exhibit non-thermal features.
Magnetic confinement fusion reactors suffer severely from heat and particle losses through turbulent transport, which has inspired the construction of ever larger and more expensive reactors. Numerical simulations are vital to their design and operat ion, but particle collisions are too infrequent for fluid descriptions to be valid. Instead, strongly magnetised fusion plasmas are described by the gyrokinetic equations, a nonlinear integro-differential system for evolving the particle distribution functions in a five-dimensional position and velocity space, and the consequent electromagnetic field. Due to the high dimensionality, simulations of small reactor sections require hundreds of thousands of CPU hours on High Performance Computing platforms. We develop a Hankel-Hermite spectral representation for velocity space that exploits structural features of the gyrokinetic system. The representation exactly conserves discrete free energy in the absence of explicit dissipation, while our Hermite hypercollision operator captures Landau damping with few variables. Calculation of the electromagnetic fields becomes purely local, eliminating inter-processor communication in, and vastly accelerating, searches for linear instabilities. We implement these ideas in SpectroGK, an efficient parallel code. Turbulent fusion plasmas may dissipate free energy through linear phase mixing to fine scales in velocity space, as in Landau damping, or through a nonlinear cascade to fine scales in physical space, as in hydrodynamic turbulence. Using SpectroGK to study saturated electrostatic drift-kinetic turbulence, we find that the nonlinear cascade suppresses linear phase mixing at energetically-dominant scales, so the turbulence is fluid-like. We use this observation to derive Fourier-Hermite spectra for the electrostatic potential and distribution function, and confirm these spectra with simulations.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا