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

Scrape-Off Layer Turbulence in Tokamaks Simulated with a Continuum Gyrokinetic Code

158   0   0.0 ( 0 )
 نشر من قبل Ammar Hakim
 تاريخ النشر 2016
  مجال البحث فيزياء
والبحث باللغة English




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

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).



قيم البحث

اقرأ أيضاً

An electrostatic gyrokinetic-based model is applied to simulate parallel plasma transport in the scrape-off layer to a divertor plate. The authors focus on a test problem that has been studied previously, using parameters chosen to model a heat pulse driven by an edge-localized mode (ELM) in JET. Previous work has used direct particle-in-cell equations with full dynamics, or Vlasov or fluid equations with only parallel dynamics. With the use of the gyrokinetic quasineutrality equation and logical sheath boundary conditions, spatial and temporal resolution requirements are no longer set by the electron Debye length and plasma frequency, respectively. This test problem also helps illustrate some of the physics contained in the Hamiltonian form of the gyrokinetic equations and some of the numerical challenges in developing an edge gyrokinetic code.
Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas (T. N. Bernard, et. al., Phys. of Plasmas 26, 042301 (2019)). Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs which can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, as well as improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.
A four-dimensional plasma model able to describe the scrape-off layer region of tokamak devices at arbitrary collisionality is derived in the drift-reduced limit. The basis of the model is provided by a drift-kinetic equation that retains the full no n-linear Coulomb collision operator and describes arbitrarily far from equilibrium distribution functions. By expanding the dependence of distribution function over the perpendicular velocity in a Laguerre polynomial basis and integrating over the perpendicular velocity, a set of four-dimensional moment equations for the expansion coefficients of the distribution function is obtained. The Coulomb collision operator, as well as Poissons equation, are evaluated explicitly in terms of perpendicular velocity moments of the distribution function.
3D2V continuum gyrokinetic simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using the full-$f$ discontinuous-Galerkin code Gkeyll. These simulations include the basic elements of a fusion-devi ce scrape-off layer: localized sources to model plasma outflow from the core, cross-field turbulent transport, parallel flow along magnetic field lines, and parallel losses at the limiter or divertor with sheath model boundary conditions. The set of sheath boundary conditions used in the model allows currents to flow through the walls. In addition to details of the numerical approach, results from numerical simulations of turbulence in the Large Plasma Device (LAPD), a linear device featuring straight magnetic field lines, are presented.
The first gyrokinetic simulations of plasma turbulence in the Texas Helimak device, a simple magnetized torus, are presented. The device has features similar to the scrape-off layer region of tokamaks, such as bad-curvature-driven instabilities and s heath boundary conditions on the end plates, which are included in these simulations. Comparisons between simulations and measurements from the experiment show similarities, including equilibrium profiles and fluctuation amplitudes that approach experimental values, but also some important quantitative differences. Both experimental and simulation results exhibit turbulence statistics that are characteristic of blob transport.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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