اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدLaguerre-Hermite Pseudo-Spectral Velocity Formulation of Gyrokinetics
48
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
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 resolution 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.
قيم البحث
اقرأ أيضاً
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. C
onsequently, 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.
Plasma turbulence is investigated using high-resolution ion velocity distributions measured by the Magnetospheric Multiscale Mission (MMS) in the Earths magnetosheath. The particle distribution is highly structured, suggesting a cascade-like process
in velocity space. This complex velocity space structure is investigated using a three-dimensional Hermite transform that reveals a power law distribution of moments. In analogy to hydrodynamics, a Kolmogorov approach leads directly to a range of predictions for this phase-space cascade. The scaling theory is in agreement with observations, suggesting a new path for the study of plasma turbulence in weakly collisional space and astrophysical plasmas.
We study the stability of spatially periodic, nonlinear Vlasov-Poisson equilibria as an eigenproblem in a Fourier-Hermite basis (in the space and velocity variables, respectively) of finite dimension, $N$. When the advection term in Vlasov equation i
s dominant, the convergence with $N$ of the eigenvalues is rather slow, limiting the applicability of the method. We use the method of spectral deformation introduced in [J. D. Crawford and P. D. Hislop, Ann. Phys. 189, 265 (1989)] to selectively damp the continuum of neutral modes associated with the advection term, thus accelerating convergence. We validate and benchmark the performance of our method by reproducing the kinetic dispersion relation results for linear (spatially homogeneous) equilibria. Finally, we study the stability of a periodic Bernstein-Greene-Kruskal mode with multiple phase space vortices, compare our results with numerical simulations of the Vlasov-Poisson system and show that the initial unstable equilibrium may evolve to different asymptotic states depending on the way it was perturbed.
Turbulence at kinetic scales is an unresolved and ubiquitous phenomenon that characterizes both space and laboratory plasmas. Recently, new theories, {it in-situ} spacecraft observations and numerical simulations suggest a novel scenario for turbulen
ce, characterized by a so-called phase space cascade -- the formation of fine structures, both in physical and velocity space. This new concept is here extended by directly taking into account the role of inter-particle collisions, modeled through the nonlinear Landau operator or the simplified Dougherty operator. The characteristic times, associated with inter-particle correlations, are derived in the above cases. The implications of introducing collisions on the phase space cascade are finally discussed.
We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309
(2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
