Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures
From MaRDI portal
Abstract: 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.
Recommendations
- Numerical Solver for the Boltzmann Equation with Self-Adaptive Collision Operators
- Collisional N-body numerical integrator with applications to charged particle dynamics
- Block-structured adaptive mesh refinement algorithms for Vlasov simulation
- On architecture and performance of adaptive mesh refinement in an electrostatics particle-in-cell code
- Boundary layer mesh generation with fast collision detection
- Leap frog integrator modifications in highly collisional particle-in-cell codes
- Adaptive mesh refinement on graphics processing units for applications in gas dynamics
- A scalable parallel method for large collision detection problems
- Adaptive mesh refinement and its application to magneto-hydrodynamics
Cites work
- A fully non-linear multi-species Fokker-Planck-Landau collision operator for simulation of fusion plasma
- A mass, momentum, and energy conserving, fully implicit, scalable algorithm for the multi-dimensional, multi-species Rosenbluth-Fokker-Planck equation
- An overview of SuperLU
- Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator
- scientific article; zbMATH DE number 3024704 (Why is no real title available?)
- p4est: scalable algorithms for parallel adaptive mesh refinement on forests of octrees
- Recursive algorithms for distributed forests of octrees
Cited in
(6)- Scalable implicit solvers with dynamic mesh adaptation for a relativistic drift-kinetic Fokker-Planck-Boltzmann model
- A multiscale hybrid Maxwellian-Monte-Carlo Coulomb collision algorithm for particle simulations
- A performance portable, fully implicit Landau collision operator with batched linear solvers
- A projection method for particle resampling
- Structure preservation using discrete gradients in the Vlasov-Poisson-Landau system
- Metriplectic relaxation to equilibria
This page was built for publication: Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4597610)