Modeling and meshing for tokamak edge plasma simulations

Abstract

The accurate representation of a problem domain in terms of a geometric model and its effective discretization into a mesh plays an important role in achieving higher-quality simulation results and better computational efficiency. In this work, recent developments to extend an automated tokamak modeling and meshing infrastructure to better support the near flux field following meshing requirements of the XGC Gyro-kinetic Code are presented. The first extension is a procedure to effectively represent a general set of O-point and X-points configurations. The second extension is a procedure to decompose the geometric model that includes the tokamak wall, selected flux curves and the separatrix curves in a manner such that an appropriate set of mesh generation procedures can be applied to each subregion. The third extension is the application of appropriate mesh generation procedures in each subregion to create the desired final mesh.