A Flow-Based Hybrid Approach for Kinetic Plasma Simulations: Bridging Direct Vlasov and Particle Methods

TL;DR

This paper introduces a flow-based kinetic plasma simulation method that unifies Vlasov and particle approaches, offering higher accuracy, reduced noise, and better scalability for complex plasma phenomena.

Contribution

A novel flow-based approach for plasma simulation that directly computes distribution functions along characteristic curves, reducing reliance on Monte Carlo sampling and enhancing accuracy.

Findings

Achieves higher accuracy with fewer markers compared to PIC.

Reduces sampling noise and computational overhead.

Demonstrates robustness in simulating Landau damping, two-stream instability, and collisional relaxation.

Abstract

We present a novel flow-based kinetic approach, inspired by continuous normalizing flows, for plasma simulation that unifies the complementary strengths of direct Vlasov solvers and particle-based methods. By tracking the distribution function along the characteristic curves defined by the Newton--Lorentz equations, our method directly computes f(z(t)) at selected points in phase space without reliance on Monte Carlo sampling. We employ a scatter-point integration scheme using smoothing kernels reminiscent of Smoothed Particle Hydrodynamics (SPH), to calculate field quantities and moments, achieving higher accuracy with far fewer markers compared to Particle-in-Cell (PIC) methods. Unlike PIC, our approach supports strategic marker placement and dynamic refinement in regions of interest, thus reducing sampling noise and computational overhead. This capability is particularly advantageous in high-density plasmas, where PIC's particle requirements can be prohibitive. In addition, the method naturally accommodates collisional effects via an augmented phase-space flow description ensuring robust handling of both collisionless and collisional plasmas. Our simulations of Landau damping, two-stream instability, and collisional relaxation demonstrate reduced noise, accurate phase-space resolution with significantly fewer markers, and robust energy conservation. Moreover, the independent characteristic curves and local scatter integration are highly amenable to GPU acceleration, enabling efficient large-scale simulations. Overall, this flow-based framework offers a powerful, flexible, and computationally efficient alternative to traditional particle methods for kinetic plasma dynamics, with potential applications spanning inertial confinement, Zpinch, and other complex kinetic systems.