Targeting realistic geometry in tokamak code Gysela
From MaRDI portal
Publication:4615491
Abstract: In magnetically confined plasmas used in Tokamak, turbulence is responsible for specific transport that limits the performance of this kind of reactors. Gyrokinetic simulations are able to capture ion and electron turbulence that give rise to heat losses, but require also state-of-the-art HPC techniques to handle computation costs. Such simulations are a major tool to establish good operating regime in Tokamak such as ITER, which is currently being built. Some of the key issues to address more realistic gyrokinetic simulations are: efficient and robust numerical schemes, accurate geometric description, good parallelization algorithms. The framework of this work is the Semi-Lagrangian setting for solving the gyrokinetic Vlasov equation and the Gyseka code. In this paper, a new variant for the interpolation method is proposed that can handle the mesh singularity in the poloidal plane at r=0 (polar system is used for the moment in Gysela). A non-uniform meshing of the poloidal plane is proposed instead of uniform one in order to save memory and computations. The interpolation method, the gyroaverage operator, and the Poisson solver are revised in order to cope with non-uniform meshes. A mapping that establish a bijection from polar coordinates to more realistic plasma shape is used to improve realism. Convergence studies are provided to establish the validity and robustness of our new approach.
Recommendations
- Development of semi-Lagrangian gyrokinetic code for full-f turbulence simulation in general tokamak geometry
- A 5D gyrokinetic full-\(f\) global semi-Lagrangian code for flux-driven ion turbulence simulations
- Guiding-center simulations on curvilinear meshes
- Field-aligned interpolation for semi-Lagrangian gyrokinetic simulations
- Optimization of the gyroaverage operator based on Hermite interpolation
Cites work
- scientific article; zbMATH DE number 2152342 (Why is no real title available?)
- A 5D gyrokinetic full-\(f\) global semi-Lagrangian code for flux-driven ion turbulence simulations
- A drift-kinetic semi-Lagrangian 4D code for ion turbulence simulation
- A note on finite difference discretizations for Poisson equation on a disk
- A parallel Vlasov solver based on local cubic spline interpolation on patches
- A simple compact fourth-order Poisson solver on polar geometry
- Elliptic Differential Equations
- Gradual time reversal in thermo- and photo-acoustic tomography within a resonant cavity
- Guiding center simulations on curvilinear grids
- Numerical treatment of polar coordinate singularities
- On the numerical solution of \(u_ t = u_{rr} +{2\over r} u_ r\)
- Solving the Vlasov equation in complex geometries
- The semi-Lagrangian method for the numerical resolution of the Vlasov equation
- The semi-Lagrangian method on curvilinear grids
Cited in
(10)- Solver comparison for Poisson-like equations on tokamak geometries
- Implicitly extrapolated geometric multigrid on disk-like domains for the gyrokinetic Poisson equation from fusion plasma applications
- Guiding-center simulations on curvilinear meshes
- New system of coordinates for tokamaks
- Energy-minimizing, symmetric discretizations for anisotropic meshes and energy functional extrapolation
- Development of a new zonal flow equation solver by diagonalisation and its application in non-circular cross-section tokamak plasmas
- Development of semi-Lagrangian gyrokinetic code for full-f turbulence simulation in general tokamak geometry
- Solving hyperbolic-elliptic problems on singular mapped disk-like domains with the method of characteristics and spline finite elements
- Field-aligned interpolation for semi-Lagrangian gyrokinetic simulations
- Optimization of the gyroaverage operator based on Hermite interpolation
This page was built for publication: Targeting realistic geometry in tokamak code Gysela
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4615491)
