-
High-order stochastic integration schemes for the Rosenbluth-Trubnikov collision operator in particle simulations
Authors:
Zhixin Lu,
Guo Meng,
Tomasz Tyranowski,
Alex Chankin
Abstract:
In this study, we consider a numerical implementation of the nonlinear Rosenbluth-Trubnikov collision operator for particle simulations in plasma physics in the framework of the finite element method (FEM). The relevant particle evolution equations are formulated as stochastic differential equations, both in the Stratonovich and Itô forms, and are then solved with advanced high-order stochastic nu…
▽ More
In this study, we consider a numerical implementation of the nonlinear Rosenbluth-Trubnikov collision operator for particle simulations in plasma physics in the framework of the finite element method (FEM). The relevant particle evolution equations are formulated as stochastic differential equations, both in the Stratonovich and Itô forms, and are then solved with advanced high-order stochastic numerical schemes. Due to its formulation as a stochastic differential equation, both the drift and diffusion components of the collision operator are treated on an equal footing. Our investigation focuses on assessing the accuracy of these schemes. Previous studies on this subject have used the Euler-Maruyama scheme, which, although popular, is of low order, and requires small time steps to achieve satisfactory accuracy. In this work, we compare the performance of the Euler-Maruyama method to other high-order stochastic methods known in the stochastic differential equations literature. Our study reveals advantageous features of these high-order schemes, such as better accuracy and improved conservation properties of the numerical solution. The main test case used in the numerical experiments is the thermalization of isotropic and anisotropic particle distributions.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Gyrokinetic analysis and simulation of pedestals, to identify the culprits for energy losses using fingerprints
Authors:
M. Kotschenreuther,
X. Liu,
D. R. Hatch,
S. Mahajan,
L. Zheng,
A. Diallo,
R. Groebner,
the DIII-D TEAM,
J. C. Hillesheim,
C. F. Maggi,
C. Giroud,
F. Koechl,
V. Parail,
S. Saarelma,
E. Solano,
JET Contributors,
A. Chankin
Abstract:
Fusion performance in tokamaks hinges critically on the efficacy of the Edge Transport Barrier (ETB) at suppressing energy losses. The new concept of fingerprints is introduced to identify the instabilities that cause the transport losses in the ETB of many of today's experiments, from widely posited candidates. Analysis of the Gyrokinetic-Maxwell equations, and gyrokinetic simulations of experime…
▽ More
Fusion performance in tokamaks hinges critically on the efficacy of the Edge Transport Barrier (ETB) at suppressing energy losses. The new concept of fingerprints is introduced to identify the instabilities that cause the transport losses in the ETB of many of today's experiments, from widely posited candidates. Analysis of the Gyrokinetic-Maxwell equations, and gyrokinetic simulations of experiments, find that each mode type produces characteristic ratios of transport in the various channels: density, heat and impurities. This, together with experimental observations of transport in some channel, or, of the relative size of the driving sources of channels, can identify or determine the dominant modes causing energy transport. In multiple ELMy H-mode cases that are examined, these fingerprints indicate that MHD-like modes are apparently not the dominant agent of energy transport; rather, this role is played by Micro-Tearing Modes (MTM) and Electron Temperature Gradient (ETG) modes, and in addition, possibly Ion Temperature Gradient (ITG)/Trapped Electron Modes (ITG/TEM) on JET. MHD-like modes may dominate the electron particle losses. Fluctuation frequency can also be an important means of identification, and is often closely related to the transport fingerprint. The analytical arguments unify and explain previously disparate experimental observations on multiple devices, including DIII-D, JET and ASDEX-U, and detailed simulations of two DIII-D ETBs also demonstrate and corroborate this.
△ Less
Submitted 24 March, 2019;
originally announced March 2019.
-
Testing nonlocal models of electron thermal conduction for magnetic and inertial confinement fusion applications
Authors:
Jonathan Peter Brodrick,
Robert J. Kingham,
Michael M. Marinak,
Mehul V. Patel,
Alex V. Chankin,
John Omotani,
Maxim Umansky,
Dario Del Sorbo,
Ben Dudson,
Joseph Thomas Parker,
Gary D. Kerbel,
Mark Sherlock,
Christopher P Ridgers
Abstract:
Three models for nonlocal electron thermal transport are here compared against Vlasov-Fokker-Planck (VFP) codes to assess their accuracy in situations relevant to both inertial fusion hohlraums and tokamak scrape-off layers. The models tested are (i) a moment-based approach using an eigenvector integral closure (EIC) originally developed by Ji, Held and Sovinec; (ii) the non-Fourier Landau-fluid (…
▽ More
Three models for nonlocal electron thermal transport are here compared against Vlasov-Fokker-Planck (VFP) codes to assess their accuracy in situations relevant to both inertial fusion hohlraums and tokamak scrape-off layers. The models tested are (i) a moment-based approach using an eigenvector integral closure (EIC) originally developed by Ji, Held and Sovinec; (ii) the non-Fourier Landau-fluid (NFLF) model of Dimits, Joseph and Umansky; and (iii) Schurtz, Nicolaï and Busquet's multigroup diffusion model (SNB). We find that while the EIC and NFLF models accurately predict the damping rate of a small-amplitude temperature perturbation (within 10% at moderate collisionalities), they overestimate the peak heat flow by as much as 35% and do not predict preheat in the more relevant case where there is a large temperature difference. The SNB model, however, agrees better with VFP results for the latter problem if care is taken with the definition of the mean free path. Additionally, we present for the first time a comparison of the SNB model against a VFP code for a hohlraum-relevant problem with inhomogeneous ionisation and show that the model overestimates the heat flow in the helium gas-fill by a factor of ~2 despite predicting the peak heat flux to within 16%.
△ Less
Submitted 6 September, 2017; v1 submitted 28 April, 2017;
originally announced April 2017.