Coupled Vlasov and two-fluid codes on GPUs

Abstract

We present a way to combine Vlasov and two-fluid codes for the simulation of a collisionless plasma in large domains while keeping full information of the velocity distribution in localized areas of interest. This is made possible by solving the full Vlasov equation in one region while the remaining area is treated by a 5-moment two-fluid code. In such a treatment, the main challenge of coupling kinetic and fluid descriptions is the interchange of physically correct boundary conditions between the different plasma models. In contrast to other treatments, we do not rely on any specific form of the distribution function, e.g. a Maxwellian type. Instead, we combine an extrapolation of the distribution function and a correction of the moments based on the fluid data. Thus, throughout the simulation both codes provide the necessary boundary conditions for each other. A speed-up factor of around 20 is achieved by using GPUs for the computationally expensive solution of the Vlasov equation and an overall factor of at least 60 using the coupling strategy combined with the GPU computation. The coupled codes were then tested on the GEM reconnection challenge.

References in corpus (6)

• A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects
• A Parallel Mesh-Adaptive Framework for Hyperbolic Conservation Laws
• Kinetic Vlasov Simulations of collisionless magnetic Reconnection
• Darwin-Vlasov Simulations of magnetized Plasmas
• Comparison of time splitting and backsubstitution methods for integrating Vlasov's equation with magnetic fields
• The Fluid-Kinetic Particle-in-Cell Solver for Plasma Simulations

Cited by in corpus (13)

• Vlasov methods in space physics and astrophysics
• PIC methods in astrophysics: Simulations of relativistic jets and kinetic physics in astrophysical systems
• Two-way coupling of magnetohydrodynamic simulations with embedded particle-in-cell simulations
• Temperature gradient driven heat flux closure in fluid simulations of collisionless reconnection
• Combining MHD and kinetic modelling of solar flares
• Magnetohydrodynamic with Adaptively Embedded Particle-in-Cell model: MHD-AEPIC
• A Parallel Low-Rank Solver for the Six-Dimensional Vlasov-Maxwell Equations
• An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
• The muphyII Code: Multiphysics Plasma Simulation on Large HPC Systems
• Enhanced conservation properties of Vlasov codes through coupling with conservative fluid models
• Adaptive Kinetic-Fluid Solvers for Heterogeneous Computing Architectures
• Multiphysics simulations of collisionless plasmas
• Verification of a hybrid gyrokinetic model using the advanced semi-Lagrange code ssV