pith. sign in

arxiv: 2306.14605 · v2 · submitted 2023-06-26 · 🧮 math.NA · cs.NA

A structure and asymptotic preserving scheme for the Vlasov-Poisson-Fokker-Planck model

Alain Blaustein (Penn State) , Francis Filbet (IMT) This is my paper

Pith reviewed 2026-05-24 08:29 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Vlasov-Poisson-Fokker-Planckstructure preserving schemeasymptotic preservinghypocoercivityHermite spectral methodfinite volume discretizationexponential relaxationkinetic plasma models
0
0 comments X

The pith

A Hermite spectral finite-volume scheme for the Vlasov-Poisson-Fokker-Planck model preserves stationary solutions and yields uniform exponential relaxation to equilibrium.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a numerical method for the Vlasov-Poisson-Fokker-Planck model that decomposes the velocity variable in a Hermite function basis and applies a finite volume discretization in space. This combination is built so that stationary solutions and a linearized free-energy estimate are preserved exactly by the discrete system. Hypocoercivity techniques are then transferred to the discrete level to obtain explicit exponential decay rates toward equilibrium for the linearized problem; these rates stay positive and uniform when the collision scaling or the mesh size changes. Numerical tests on the nonlinear equations confirm that the same discretization handles plasma echoes at low collision rates and relaxation at high collision rates while remaining stable for any parameter choice.

Core claim

The scheme naturally preserves both stationary solutions and linearized free-energy estimate. Adaptation of hypocoercivity methods provides quantitative estimates ensuring the exponential relaxation to equilibrium of the discrete solution for the linearized Vlasov-Poisson-Fokker-Planck system, uniformly with respect to both scaling and discretization parameters.

What carries the argument

Hermite function spectral decomposition in velocity combined with a structure-preserving finite volume scheme in space, which carries the preservation properties and permits direct application of hypocoercivity estimates at the discrete level.

If this is right

  • The discrete solution relaxes exponentially to equilibrium for the linearized system at a rate independent of collision scaling and mesh parameters.
  • The scheme remains unconditionally stable and asymptotically preserving across all collisional regimes.
  • It reproduces plasma echo behavior in weakly collisional regimes and monotonic relaxation in strongly collisional regimes.
  • Long-time simulations of the nonlinear system remain consistent without artificial numerical dissipation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The uniform decay estimates imply that the discrete model can be used for arbitrarily long integration times without retuning parameters when collision strength varies.
  • The structure-preserving construction may allow similar Hermite-finite-volume combinations for related kinetic equations that possess a free-energy dissipation structure.
  • Because the estimates hold independently of discretization, the method supplies a candidate for multi-scale plasma simulations where local collision frequency changes by orders of magnitude.

Load-bearing premise

Hypocoercivity arguments transfer to the discrete Hermite-finite-volume system without extra conditions on grid size or collision scaling.

What would settle it

A direct computation of the decay rate of the linearized discrete solution on a sequence of successively refined meshes showing that the rate drops to zero or becomes mesh-dependent for some fixed collision scaling.

Figures

Figures reproduced from arXiv: 2306.14605 by Alain Blaustein (Penn State), Francis Filbet (IMT).

Figure 4.1
Figure 4.1. Figure 4.1: Asymptotic-preserving properties: time development of (a) the potential energy (b) ∥f − f∞∥ L2  f −1 ∞  (in log scale). As predicted by our analysis of the linearized model, we observe exponential relaxation towards equi￾librium at a rate which is proportional to 1/ε. The scheme is uniformly stable with respect to ε and the solution converges to the discrete equilibrium f∞ when ε → 0 [PITH_FULL_IMAGE:… view at source ↗
Figure 4.2
Figure 4.2. Figure 4.2: Asymptotic-preserving properties: re-scaled time development (s ← ε t) of the potential energy (in log scale). The left chart in [PITH_FULL_IMAGE:figures/full_fig_p014_4_2.png] view at source ↗
Figure 4.3
Figure 4.3. Figure 4.3: Perturbation of non uniform density for τ0 = 104 (weakly collisional regime): time development of the potential energy in log scale (for (a) the linearized Vlasov-Poisson-Fokker-Planck system and (b) the nonlinear Vlasov-Poisson-Fokker-Planck system. In [PITH_FULL_IMAGE:figures/full_fig_p015_4_3.png] view at source ↗
Figure 4.4
Figure 4.4. Figure 4.4: Perturbation of non uniform density (weakly collisional regime, τ0 = 104 ): snapshots of the difference between the solution f and the equilibrium f∞ at time t = 4, 8, 16, 30, 40 and 70. On the other hand, we study the influence of the collision frequency τ0 and perform several numerical simulations for the nonlinear Vlasov-Poisson-Fokker-Planck system (1.2) with the same initial data for τ0 = 10k , with… view at source ↗
Figure 4.5
Figure 4.5. Figure 4.5: Perturbation of non uniform density: time development of (a) the potential energy (b) ∥f −f∞∥ L2  f −1 ∞  for various τ0 = 1, . . . , 104 (in log scale). 4.3. Plasma echo. We now investigate a much more intricate problem where the non-linearity plays the main role. Following the work [24, 33] or more recently [1, 25], we will consider a perturbation of an homogeneous Maxwellian distribution f∞(x, v) :=… view at source ↗
Figure 4.6
Figure 4.6. Figure 4.6: Perturbation of non uniform density (moderate collisional regime, τ0 = 102 ): snapshots of the difference between the solution f and the equilibrium f∞ at time t = 4, 8, 16, 30, 40 and 70. nonlinear system. The numerical solution corresponding to the linearized system exhibits a simple Landau damping, when t ≥ 5, with a decay rate corresponding to the predicted value γL = 0.355, whereas the numerical sol… view at source ↗
Figure 4.7
Figure 4.7. Figure 4.7: Plasma echo for τ0 = 106 (weakly collisional regime): time development of the potential energy (top) and square of the k-th mode of the electric field for k = k1, ..., 4k1 in log scale (bottom) for (a) the linearized Vlasov￾Poisson-Fokker-Planck system and (b) the nonlinear Vlasov-Poisson-Fokker-Planck system. For this weakly collisional regime (τ0 = 106 ), we also report (on [PITH_FULL_IMAGE:figures/fu… view at source ↗
Figure 4.8
Figure 4.8. Figure 4.8: Plasma echo for τ0 = 106 (weakly collisional regime): time development of ∥f − f∞∥L2(f −1 ∞ ) , ∥f − ρM∥L2(f −1 ∞ ) and ∥ρ − ρ∞∥L2(ρ −1 ∞ ) in log scale for (a) the linearized Vlasov-Poisson-Fokker-Planck system and (b) the nonlinear Vlasov-Poisson-Fokker-Planck system. This can be also viewed in [PITH_FULL_IMAGE:figures/full_fig_p020_4_8.png] view at source ↗
Figure 4.9
Figure 4.9. Figure 4.9: Plasma echo for τ0 = 106 (weakly collisional regime): snapshots of the difference between the solution f and the equilibrium f∞ at time t = 0, 20, 30, 40 and 50. with 2 ≤ k ≤ 6, and also on the bottom charts of [PITH_FULL_IMAGE:figures/full_fig_p021_4_9.png] view at source ↗
Figure 4.10
Figure 4.10. Figure 4.10: Plasma echo for τ0 = 1, . . . , 106 (various regimes): time development of the potential energy (top) and ∥f − f∞∥L2(f −1 ∞ ) (bottom) in log scale for (a) the linearized Vlasov-Poisson-Fokker-Planck system and (b) the nonlinear Vlasov-Poisson-Fokker-Planck system. also be observed on [PITH_FULL_IMAGE:figures/full_fig_p022_4_10.png] view at source ↗
Figure 4.11
Figure 4.11. Figure 4.11: Two-stream : time development of (a) ∥E − E∞∥L2(T) (b) and ∥f − f∞∥L2(f −1 ∞ ) for various τ0 = 102 , . . . , 106 (in log scale). (a) τ0 = 1 (b) τ0 = 10 (c) τ0 = 100 (d) τ0 = 1000 [PITH_FULL_IMAGE:figures/full_fig_p023_4_11.png] view at source ↗
Figure 4.12
Figure 4.12. Figure 4.12: Two-stream : time development of ∥f − f∞∥L2(f −1 ∞ ) , ∥f − ρM∥L2(f −1 ∞ ) and ∥ρ − ρ∞∥L2(ρ −1 ∞ ) for various τ0 (in log scale). equation, to the nonlinear scheme (2.11)-(2.16) by proving its asymptotic preserving properties and expo￾nential trend towards equilibrium of discrete solutions. This might be doable in a perturbative setting by controlling the nonlinear contribution using discrete Sobolev in… view at source ↗
Figure 4.13
Figure 4.13. Figure 4.13: Two-stream : snapshots of the distribution function f at time t = 8, 16, 30 and 60 for various τ0. result has been laid in [6], where we propagated discrete H1 norms in the linear setting. Regarding simulations, the study of echoes also raises interesting perspectives. In [25] were constructed theoretical solutions to the Vlasov-Poisson equation which display infinite cascades of echoes and for which La… view at source ↗
read the original abstract

We propose a numerical method for the Vlasov-Poisson-Fokker-Planck model written as an hyperbolic system thanks to a spectral decomposition in the basis of Hermite functions with respect to the velocity variable and a structure preserving finite volume scheme for the space variable. On the one hand, we show that this scheme naturally preserves both stationary solutions and linearized free-energy estimate. On the other hand, we adapt previous arguments based on hypocoercivity methods to get quantitative estimates ensuring the exponential relaxation to equilibrium of the discrete solution for the linearized Vlasov-Poisson-Fokker-Planck system, uniformly with respect to both scaling and discretization parameters. Finally, we perform substantial numerical simulations for the nonlinear system to illustrate the efficiency of this approach for a large variety of collisional regimes (plasma echos for weakly collisional regimes and trend to equilibrium for collisional plasmas) and to highlight its robustness (unconditional stability, asymptotic preserving properties).

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript proposes a Hermite spectral decomposition in velocity combined with a structure-preserving finite-volume scheme in space for the Vlasov-Poisson-Fokker-Planck system. It establishes direct preservation of stationary solutions and the linearized free-energy dissipation from the scheme stencil, adapts hypocoercivity arguments to derive quantitative exponential relaxation rates for the linearized discrete system that are uniform with respect to the collision scaling and discretization parameters (mesh size and Hermite truncation), and presents numerical experiments on the nonlinear system demonstrating asymptotic-preserving behavior and unconditional stability across weakly to strongly collisional regimes.

Significance. If the uniform hypocoercivity estimates hold as stated, the work supplies a theoretically grounded asymptotic-preserving method for kinetic plasma models that bridges collisional regimes without parameter-dependent restrictions. The explicit verification of structure preservation directly from the stencil and the tracking of constants in the discrete hypocoercivity analysis are notable strengths that support robustness claims.

minor comments (3)
  1. [Section 4] The statement of the main hypocoercivity theorem would benefit from an explicit remark on the precise range of admissible Hermite truncation degrees for which the constant uniformity holds without further restrictions.
  2. [Section 5] In the numerical section, the choice of time-step size relative to the spatial mesh and Hermite degree is not tabulated; adding a short discussion or table entry would clarify the unconditional stability claim.
  3. A few figure captions (e.g., those showing plasma echo tests) omit the precise values of the collision frequency and mesh parameters used, which would aid reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The referee's summary correctly identifies the main contributions: the structure-preserving properties derived directly from the stencil, the uniform hypocoercivity estimates for the linearized discrete system, and the numerical demonstration of asymptotic-preserving behavior across collisional regimes.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The claimed preservation of stationary solutions and linearized free-energy dissipation is verified directly from the finite-volume stencil and Hermite spectral discretization without reduction to fitted inputs. The hypocoercivity adaptation for uniform exponential relaxation is carried out with explicit constant tracking that remains independent of collision scaling and discretization parameters (mesh size, truncation degree). No self-citation chain, ansatz smuggling, or renaming of known results is load-bearing for the central claims; the derivation chain is self-contained against the stated assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review yields no identifiable free parameters, invented entities, or non-standard axioms; the approach relies on established spectral and finite-volume techniques whose details are not supplied.

axioms (1)
  • standard math Standard properties of Hermite orthogonal polynomials and finite-volume conservation
    Invoked for the velocity decomposition and spatial discretization.

pith-pipeline@v0.9.0 · 5704 in / 1133 out tokens · 41158 ms · 2026-05-24T08:29:42.389454+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Bedrossian

    J. Bedrossian. Suppression of Plasma Echoes and Landau Damping in Sobolev Spaces by Weak Collisions in a Vlasov- Fokker-Planck Equation. Annals of PDE , 3(19), 2017. 24

    work page 2017

  2. [2]

    Bessemoulin-Chatard and F

    M. Bessemoulin-Chatard and F. Filbet. A finite volume scheme for nonlinear degenerate parabolic equations. SIAM Journal on Scientific Computing , 34(5):B559–B583, 2012

    work page 2012

  3. [3]

    Bessemoulin-Chatard and F

    M. Bessemoulin-Chatard and F. Filbet. On the stability of conservative discontinuous Galerkin/Hermite spectral methods for the Vlasov-Poisson system. J. Comput. Phys. , 451:Paper No. 110881, 28, 2022

    work page 2022

  4. [4]

    Bessemoulin-Chatard and F

    M. Bessemoulin-Chatard and F. Filbet. On the convergence of discontinuous Galerkin/Hermite spectral methods for the Vlasov-Poisson system. SIAM J. Num. Analysis , 2023

    work page 2023

  5. [5]

    Bessemoulin-Chatard, M

    M. Bessemoulin-Chatard, M. Herda, and T. Rey. Hypocoercivity and diffusion limit of a finite volume scheme for linear kinetic equations. Math. Comp., 89(323):1093–1133, 2020

    work page 2020

  6. [6]

    Blaustein and F

    A. Blaustein and F. Filbet. On a discrete framework of hypocoercivity for kinetic equations. to appear in AMS Math. Computation, 2023

    work page 2023

  7. [7]

    Burger, J

    M. Burger, J. A. Carrillo, and M.-T. Wolfram. A mixed finite element method for nonlinear diffusion equations. Kinetic & Related Models, 3(1):59, 2010

    work page 2010

  8. [8]

    J. A. Carrillo, L. Wang, W. Xu, and M. Yan. Variational Asymptotic Preserving Scheme for the Vlasov–Poisson–Fokker– Planck System. Multiscale Modeling & Simulation , 19(1):478–505, 2021

    work page 2021

  9. [9]

    Chainais-Hillairet and F

    C. Chainais-Hillairet and F. Filbet. Asymptotic behaviour of a finite-volume scheme for the transient drift-diffusion model. IMA journal of numerical analysis , 27(4):689–716, 2007

    work page 2007

  10. [10]

    Chang and G

    J. Chang and G. Cooper. A practical difference scheme for Fokker-Planck equations. Journal of Computational Physics , 6(1):1–16, 1970

    work page 1970

  11. [11]

    Chaturvedi, J

    S. Chaturvedi, J. Luk, and T. T. Nguyen. The Vlasov–Poisson–Landau system in the weakly collisional regime, 2022

    work page 2022

  12. [12]

    Coughlin and J

    J. Coughlin and J. Hu. Efficient dynamical low-rank approximation for the Vlasov-Amp` ere-Fokker-Planck system. J. Comput. Phys., 470:Paper No. 111590, 20, 2022

    work page 2022

  13. [13]

    Crestetto, N

    A. Crestetto, N. Crouseilles, and M. Lemou. Kinetic/fluid micro-macro numerical schemes for Vlasov-Poisson-BGK equa- tion using particles. Kinet. Relat. Models, 5(4):787–816, 2012

    work page 2012

  14. [14]

    Degond and F

    P. Degond and F. Deluzet. Asymptotic-preserving methods and multiscale models for plasma physics. J. Comput. Phys. , 336:429–457, 2017

    work page 2017

  15. [15]

    Desvillettes and C

    L. Desvillettes and C. Villani. On the trend to global equilibrium in spatially inhomogeneous entropy-dissipating systems: the linear Fokker-Planck equation. Comm. Pure Appl. Math. , 54(1):1–42, 2001

    work page 2001

  16. [16]

    Dolbeault, C

    J. Dolbeault, C. Mouhot, and C. Schmeiser. Hypocoercivity for linear kinetic equations conserving mass. Transactions of the American Mathematical Society , 367(6):3807–3828, 2015

    work page 2015

  17. [17]

    Dujardin, F

    G. Dujardin, F. H´ erau, and P. Lafitte. Coercivity, hypocoercivity, exponential time decay and simulations for discrete Fokker–Planck equations. Numerische Mathematik, 144(3):615–697, 2020

    work page 2020

  18. [18]

    Einkemmer, J

    L. Einkemmer, J. Hu, and Y. Wang. An asymptotic-preserving dynamical low-rank method for the multi-scale multi- dimensional linear transport equation. J. Comput. Phys. , 439:Paper No. 110353, 21, 2021

    work page 2021

  19. [19]

    F. Filbet. Convergence of a finite volume scheme for the Vlasov-Poisson system. SIAM J. Numer. Anal. , 39(4):1146–1169, 2001

    work page 2001

  20. [20]

    Filbet and C

    F. Filbet and C. Negulescu. Fokker-Planck multi-species equations in the adiabatic asymptotics. J. Comput. Phys. , 471:Paper No. 111642, 2022

    work page 2022

  21. [21]

    Filbet and E

    F. Filbet and E. Sonnendr¨ ucker. Comparison of Eulerian Vlasov solvers. Comput. Phys. Comm. , 150(3):247–266, 2003

    work page 2003

  22. [22]

    Filbet and T

    F. Filbet and T. Xiong. Conservative Discontinuous Galerkin/Hermite Spectral Method for the Vlasov–Poisson System. Commun. Appl. Math. Comput. , 2020

    work page 2020

  23. [23]

    E. H. Georgoulis. Hypocoercivity-compatible finite element methods for the long-time computation of kolmogorov’s equa- tion. SIAM Journal on Numerical Analysis , 59(1):173–194, 2021

    work page 2021

  24. [24]

    Plasma wave echo

    Gould, R W and O’Neil, TM and Malmberg, JH. Plasma wave echo. Physical Review Letters, 19(5):219, 1967

    work page 1967

  25. [25]

    Grenier, T

    E. Grenier, T. T. Nguyen, and I. Rodnianski. Plasma Echoes Near Stable Penrose Data. SIAM Journal on Mathematical Analysis, 54(1):940–953, 2022

    work page 2022

  26. [26]

    K. J. Havlak and H. D. Victory Jr. On deterministic particle methods for solving Vlasov–Poisson–Fokker–Planck systems. SIAM journal on numerical analysis , 35(4):1473–1519, 1998

    work page 1998

  27. [27]

    H´ erau and L

    F. H´ erau and L. Thomann. On global existence and trend to the equilibrium for the Vlasov-Poisson-Fokker-Planck system with exterior confining potential. J. Funct. Anal., 271(5):1301–1340, 2016

    work page 2016

  28. [28]

    Herda and M

    M. Herda and M. Rodrigues. Large-time behavior of solutions to Vlasov-Poisson-Fokker-Planck equations: from evanescent collisions to diffusive limit. J. Stat. Phys. , 170(5):895–931, 2018

    work page 2018

  29. [29]

    H. J. Hwang and J. Jang. On the Vlasov-Poisson-Fokker-Planck equation near Maxwellian. Discrete Contin. Dyn. Syst. Ser. B, 18(3):681–691, 2013

    work page 2013

  30. [30]

    J. Jang, F. Li, J.-M. Qiu, and T. Xiong. High order asymptotic preserving DG-IMEX schemes for discrete-velocity kinetic equations in a diffusive scaling. J. Comput. Phys. , 281:199–224, 2015

    work page 2015

  31. [31]

    X. S. Li. An overview of SuperLU: Algorithms, implementation, and user interface. ACM Trans. Math. Softw., 31(3):302– 325, September 2005

    work page 2005

  32. [32]

    Liu and S

    L. Liu and S. Jin. Hypocoercivity based sensitivity analysis and spectral convergence of the stochastic galerkin ap- proximation to collisional kinetic equations with multiple scales and random inputs. Multiscale Modeling & Simulation , 16(3):1085–1114, 2018

    work page 2018

  33. [33]

    Malmberg, C

    J. Malmberg, C. Wharton, R. Gould, and T. O’Neil. Plasma wave echo experiment. Physical Review Letters , 20(3):95, 1968

    work page 1968

  34. [34]

    Pareschi and T

    L. Pareschi and T. Rey. Residual equilibrium schemes for time dependent partial differential equations. Computers & Fluids, 156:329–342, 2017

    work page 2017

  35. [35]

    Sonnendr¨ ucker, F

    E. Sonnendr¨ ucker, F. Filbet, A. Friedman, E. Oudet, and J.-L. Vay. Vlasov simulations of beams with a moving grid. Computer Physics Communications , 164(1-3):390–395, 2004. 25

    work page 2004

  36. [36]

    C. Villani. Hypocoercivity. Memoirs Amer. Math. Soc. , 2009

    work page 2009

  37. [37]

    Wollman and E

    S. Wollman and E. Ozizmir. Numerical approximation of the Vlasov–Poisson–Fokker–Planck system in one dimension. Journal of Computational Physics , 202(2):602–644, 2005. 26

    work page 2005