numerical analysis

An Adaptive Fourier Spectral Method for the Vlasov-Poisson system

arXiv:2607.12368

summary

The paper introduces an adaptive Fourier spectral method with high-order time splitting to dynamically expand the wavenumber domain and avoid artificial filtering when simulating the Vlasov‑Poisson system.

Abstract

Numerical simulation of the Vlasov-Poisson system faces fundamental challenges due to phase-space filamentation. Standard spectral methods rely on artificial filtering to suppress errors, which inadvertently degrades physical structures and accuracy over time. This paper proposes a dynamically adaptive Fourier spectral method combined with high-order time splitting to overcome these limitations. To prevent filamentation-induced aliasing, we dynamically expand the wavenumber domain via adaptive zero-padding whenever spectral tail monitoring detects emerging fine-scale structures. Acting as an exact trigonometric interpolation, it inherently avoids the artificial smoothing of traditional filters, preserving mass and effectively maintaining the -norm within the dynamically resolved scales. Numerical experiments on Landau damping and various instability problems validate the robustness of our scheme.

Topics & keywords

#vlasov-poisson#spectral methods#adaptive algorithms#phase-space filamentation#time splittingFourier spectral methodadaptive zero-paddingtrigonometric interpolationLandau dampingaliasing control
An Adaptive Fourier Spectral Method for the Vlasov-Poisson system · wovepaper