A fast spectral particle method for the Landau equation
Abstract
We propose a fast deterministic particle method for the spatially homogeneous Landau equation. The method combines a particle representation of the solution with a Fourier approximation of the nonlinear collision flux. Nonuniform fast Fourier transforms are used to reconstruct the density from the particles and to evaluate the flux and density at the particle locations, while the convolutional structure of the flux enables its efficient computation by FFTs. For a fixed transform tolerance, the cost per time step is 𝒪(N+Mdlog M), where N is the number of particles and M the number of Fourier modes per velocity dimension. We establish a consistency estimate for the reconstructed velocity field on regions where the reference density is bounded away from zero. The estimate separates the spectral truncation error from the particle-density reconstruction error, showing how the latter can dominate for sufficiently smooth densities. Two-dimensional numerical experiments with Maxwellian interactions illustrate the accuracy and efficiency of the method, its sensitivity to particle and spectral resolution, and the effects of filtering. Comparisons with a direct blob implementation demonstrate improved accuracy at a strongly reduced computational cost.
Related
- Weak-strong uniqueness for the Landau equation by a relative entropy method
- A fast Laguerre--Fourier spectral method for the polyatomic Boltzmann collision operator with continuous internal energy
- A Fully Discrete, Mass-Conservative Adaptive Discontinuous Spectral-Element Method for the One-Dimensional Richards Equation in Mixed Form
- Landau Levels of Landau-Raychaudhuri Equation and Wheeler-DeWitt Crystals.
- Topological defect trains excited by a boundary discontinuity in the two-dimensional complex Ginzburg-Landau equation
BibTeX
@misc{borghi2026fast,
title = {{A fast spectral particle method for the Landau equation}},
author = {Giacomo Borghi and Lorenzo Pareschi},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.23086},
eprint = {2609.23086},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.23086},
}
RIS
TY - UNPB
TI - A fast spectral particle method for the Landau equation
AU - Giacomo Borghi
AU - Lorenzo Pareschi
PY - 2026
DA - 2026/09/19
PB - arXiv
DO - 10.48550/arxiv.2609.23086
UR - https://arxiv.org/abs/2609.23086
AB - We propose a fast deterministic particle method for the spatially homogeneous Landau equation. The method combines a particle representation of the solution with a Fourier approximation of the nonlinear collision flux. Nonuniform fast Fourier transforms are used to reconstruct the density from the particles and to evaluate the flux and density at the particle locations, while the convolutional structure of the flux enables its efficient computation by FFTs. For a fixed transform tolerance, the cost per time step is 𝒪(N+M^dlog M), where N is the number of particles and M the number of Fourier modes per velocity dimension. We establish a consistency estimate for the reconstructed velocity field on regions where the reference density is bounded away from zero. The estimate separates the spectral truncation error from the particle-density reconstruction error, showing how the latter can dominate for sufficiently smooth densities. Two-dimensional numerical experiments with Maxwellian interactions illustrate the accuracy and efficiency of the method, its sensitivity to particle and spectral resolution, and the effects of filtering. Comparisons with a direct blob implementation demonstrate improved accuracy at a strongly reduced computational cost.
ER -