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. Here, we investigate a high order accurate, fully conservative, finite element discretization of the nonlinear multispecies Landau integral with adaptive mesh refinement using the P...

Full description

Saved in:
Bibliographic Details
Published in:SIAM journal on scientific computing 2017-01, Vol.39 (6), p.C452-C465
Main Authors: Adams, Mark F., Hirvijoki, Eero, Knepley, Matthew G., Brown, Jed, Isaac, Tobin, Mills, Richard
Format: Article
Language:English
Subjects:
fusion plasma physics
Landau collision integral
MATHEMATICS AND COMPUTING
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The Landau collision integral is an accurate model for the small-angle dominated Coulomb collisions in fusion plasmas. Here, we investigate a high order accurate, fully conservative, finite element discretization of the nonlinear multispecies Landau integral with adaptive mesh refinement using the PETSc library (ulwww.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:1064-8275
1095-7197
DOI:10.1137/17M1118828