Related papers: gyaradax: Local Gyrokinetics JAX Code
The first nonlinear gyrokinetic simulations obtained using a moment approach based on the Hermite-Laguerre decomposition of the distribution function are presented, implementing advanced models for the collision operator. Turbulence in a…
In robot control, planning, and learning, there is a need for rigid-body dynamics libraries that are highly performant, easy to use, and compatible with CPUs and accelerators. While existing libraries often excel at either low-latency CPU…
Predicting the dynamics of a thermonuclear plasma during a magnetic confinement experiment is fundamental in order to make nuclear fusion a reliable source of energy. The development of a set of equations describing the plasma evolution on…
Simulation-based plasma scenario development, optimization and control are crucial elements towards the successful deployment of next-generation experimental tokamaks and Fusion power plants. Current simulation codes require extremely…
One of the most well-established codes for modeling non-linear Magnetohydrodynamics (MHD) for tokamak reactors is JOREK, which solves these equations with a B\'ezier surface based finite element method. This code produces a highly sparse…
Integrating computational fluid dynamics (CFD) software into optimization and machine-learning frameworks is hampered by the rigidity of classic computational languages and the slow performance of more flexible high-level languages.…
A new global two-fluid electromagnetic turbulence code, CENTORI, has been developed for the purpose of studying magnetically-confined fusion plasmas on energy confinement timescales. This code is used to evolve the combined system of…
The field of plasma physics heavily relies on simulations to model various phenomena, such as instabilities, turbulence, and nonlinear behaviors that would otherwise be difficult to study from a purely theoretical approach. Simulations are…
We present new stellarator equilibria that have been optimized for reduced turbulent transport using nonlinear gyrokinetic simulations within the optimization loop. The optimization routine involves coupling the pseudo-spectral GPU-native…
We have used the local-$\delta{f}$ gyrokinetic code GS2 to perform studies of the effect of flux-surface shaping on two highly-shaped, low- and high-$\beta$ JT-60SA-relevant equilibria, including a successful benchmark with the GKV code. We…
Optimizing the performance of magnetic confinement fusion devices is critical to achieving an attractive fusion reactor design. Negative triangularity (NT) scenarios have been shown to achieve excellent levels of energy confinement, while…
We developed a portable code for dissipative particle dynamics (DPD) simulations. This Fortran program named CAMUS has a couple of notable features. One is the omission of constructing the so-called neighboring particles list, providing a…
Recent development of real-time equilibrium code Equinox [1] using a fixed-point algorithm [2] allow major plasma magnetic parameters to be identified in real-time, using rigorous analytical method. The code relies on the boundary flux code…
A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure…
In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube, full-flux-surface, and radially global domains. Two problems are studied: the linear relaxation of zonal flows and…
The main idea of Gyrokinetic dynamical reduction consists in systematical removing of fastest scale of motion (the gyro motion) from plasma's dynamics, resulting in a considerable model simplification and gain of computing time. Gyrokinetic…
The nonlinear gyrokinetic equations describe plasma turbulence in laboratory and astrophysical plasmas. To solve these equations, massively parallel codes have been developed and run on present-day supercomputers. This paper describes…
Previous limiter-biasing experiments on the Texas Helimak, a simple magnetized torus, have been inconclusive on the effect of flow shear on turbulence levels. To investigate this, the first gyrokinetic simulations of limiter biasing in the…
Full-f global gyrokinetic simulations of the plasma boundary have until now required heroic computational efforts and case-by-case expert intervention, precluding systematic parameter scans. Here we demonstrate a paradigm shift: hundreds of…
We present linrax, the first simplex based linear program (LP) solver compatible with the JAX ecosystem. In many control algorithms, LPs are often automatically generated and frequently solved either offline or online in the control loop.…