Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures

The Landau collision integral is an accurate model for the small-angle dominated Coulomb collisions in fusion plasmas. We investigate a high order accurate, fully conservative, finite element discretization of the nonlinear multi-species Landau integral with adaptive mesh refinement using the PETSc...

Full description

Saved in:
Bibliographic Details
Published inarXiv.org
Main Authors Adams, M F, Hirvijoki, E, Knepley, M G, Brown, J, Isaac, T, Mills, R
Format Paper Journal Article
LanguageEnglish
Published Ithaca Cornell University Library, arXiv.org 01.03.2017
Subjects
Algorithms
Computer Science - Computational Engineering, Finance, and Science
Computer Science - Mathematical Software
Coulomb collisions
Finite element method
Grid refinement (mathematics)
Integrals
Microprocessors
Plasmas (physics)
Vector processing (computers)
Online AccessGet full text

Cover

More Information
Summary:The Landau collision integral is an accurate model for the small-angle dominated Coulomb collisions in fusion plasmas. We investigate a high order accurate, fully conservative, finite element discretization of the nonlinear multi-species Landau integral with adaptive mesh refinement using the PETSc library (www.mcs.anl.gov/petsc). We develop algorithms and techniques to efficiently utilize emerging architectures with an approach that minimizes memory usage and movement and is suitable for vector processing. The Landau collision integral is vectorized with Intel AVX-512 intrinsics and the solver sustains as much as 22% of the theoretical peak flop rate of the Second Generation Intel Xeon Phi, Knights Landing, processor.
ISSN:2331-8422
DOI:10.48550/arxiv.1702.08880