pith. sign in

arxiv: 2604.18871 · v1 · submitted 2026-04-20 · 🧮 math.PR · math.AP

Quantitative approximation of the Vlasov(-Fokker-Planck)-Navier-Stokes system by stochastic particle systems

Ludovic Gouden\`ege , Christian Olivera , Gabriela Planas , Alexandre Richard This is my paper

Pith reviewed 2026-05-10 03:11 UTC · model grok-4.3

classification 🧮 math.PR math.AP
keywords Vlasov-Navier-Stokesstochastic particlesmean-field limitconvergence ratesfluid-particle systemNavier-StokesVlasov equationdrag force
0
0 comments X

The pith

A system of N stochastic particles approximates the Vlasov-Navier-Stokes equations with an explicit convergence rate in dimensions 2 and 3.

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

This paper proves that in two and three dimensions, a system of N interacting stochastic particles coupled with a fluid converges at a positive rate to the Vlasov(-Fokker-Planck)-Navier-Stokes system as N tends to infinity. The empirical measure of the particles approaches the kinetic density in the Vlasov or Vlasov-Fokker-Planck equation, while the fluid velocity field approaches the solution of the incompressible Navier-Stokes equations. This quantitative approximation is relevant for modeling two-phase flows such as aerosols and sprays, where particle methods are used in practice. The analysis combines stochastic calculus for the particles with PDE techniques to obtain energy estimates and commutator estimates that control the difference between discrete and continuous systems.

Core claim

In dimensions d in {2,3}, the empirical measure associated with the N-particle system converges to the Vlasov(-Fokker-Planck) component, while the fluid velocity converges to the Navier-Stokes component of the coupled system, with an explicit rate as N to infinity. This holds both for the case with noise and for the case where noise vanishes as N increases, leading to the Vlasov-Navier-Stokes system.

What carries the argument

The empirical measure of the N interacting stochastic particles, whose difference from the continuous Vlasov density is controlled by energy estimates and commutator estimates obtained via stochastic calculus and PDE methods.

If this is right

  • The approximation error between the particle system and the target PDE system goes to zero at a controlled rate as N increases.
  • Both the Vlasov-Navier-Stokes and Vlasov-Fokker-Planck-Navier-Stokes cases are covered by the convergence result.
  • The proofs establish parallel energy and commutator estimates for the discrete particle system and the continuous PDE system.
  • Convergence of the fluid velocity field accompanies the convergence of the particle distribution.

Where Pith is reading between the lines

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

  • This result supports the use of stochastic particle methods for accurate simulation of aerosol and spray dynamics with predictable error bounds.
  • The approach of vanishing noise may be adaptable to derive other deterministic limits from stochastic particle systems in fluid mechanics.
  • Similar quantitative rates could potentially be obtained for related coupled systems, such as those with different interaction forces or in bounded domains.
  • The combination of stochastic and deterministic estimates might inform numerical schemes that preserve the convergence properties.

Load-bearing premise

The initial data and the solutions of the systems are assumed to have sufficient regularity to allow the energy estimates and commutator estimates to close.

What would settle it

A counterexample where the convergence rate does not hold for a specific choice of initial data satisfying the paper's assumptions, or a numerical experiment showing no convergence for large N, would falsify the claim.

read the original abstract

This paper is concerned with a fluid-particle system given by the incompressible Navier-Stokes equations coupled with the Vlasov(-Fokker-Planck) equation through a drag force. Such a model arises naturally in the study of aerosols, sprays, and more generally two-phase flows. In dimensions $d\in \{2,3\}$, we establish a rate of convergence for a system of $N$ interacting stochastic particles coupled with a fluid, towards the Vlasov(-Fokker-Planck)-Navier-Stokes system, as $N\to \infty$. The case of particles with a noise that vanishes as $N\to \infty$ is considered and leads specifically to the Vlasov-Navier-Stokes system. More precisely, we prove that the empirical measure associated with the particle system converges to the Vlasov(-Fokker-Planck) component, while the fluid velocity converges to the Navier-Stokes component of the coupled system. The proofs combine stochastic calculus and PDE techniques to establish energy estimates and commutator estimates for both the discrete and continuous systems.

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

2 major / 2 minor

Summary. The paper establishes quantitative convergence rates as N→∞ for a system of N interacting stochastic particles coupled to an incompressible fluid, approximating the Vlasov(-Fokker-Planck)-Navier-Stokes system in dimensions d=2,3. The empirical measure of the particles converges to the Vlasov(-Fokker-Planck) component while the fluid velocity converges to the Navier-Stokes component. Proofs combine stochastic calculus with PDE energy estimates and commutator estimates; the vanishing-noise case yields the Vlasov-Navier-Stokes limit.

Significance. If the estimates close under the stated assumptions, the work supplies explicit error rates for stochastic particle approximations of two-phase fluid models arising in aerosols and sprays. This is valuable for both the theory of mean-field limits for coupled kinetic-fluid systems and for numerical error control. The direct use of stochastic calculus to derive the rates, rather than purely qualitative propagation of chaos, is a clear strength.

major comments (2)
  1. [§2, Theorem 2.1] §2, Theorem 2.1 (main convergence statement): the rate for d=3 is stated without explicit restriction on the time interval or smallness/regularity assumptions on the data. The commutator estimates in §4 (controlling the difference between the empirical measure and the Vlasov component) close only when ||∇u||_∞ remains bounded, but global smooth solutions to the 3D Navier-Stokes equations are not known to exist for arbitrary initial data. The theorem must therefore be qualified as local-in-time or conditional on small data to support the claimed global rate.
  2. [§3] §3 (model assumptions and initial data): the precise scaling of the noise term (vanishing as N→∞) and the form of the drag force are load-bearing for the energy estimates to close, yet the precise Sobolev or moment assumptions on the initial particle distribution and fluid velocity are not listed explicitly in the theorem statement. These must be stated verbatim so that the reader can verify the estimates apply.
minor comments (2)
  1. [§1] The introduction would benefit from a short paragraph contrasting the quantitative rates obtained here with prior qualitative mean-field results for Vlasov-Navier-Stokes systems.
  2. [Notation section] Notation for the empirical measure μ^N_t and the limiting density f_t should be introduced once and used consistently; occasional switches between μ^N and the particle positions X^i_t create minor ambiguity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and will incorporate the suggested clarifications to strengthen the presentation.

read point-by-point responses

  1. Referee: [§2, Theorem 2.1] §2, Theorem 2.1 (main convergence statement): the rate for d=3 is stated without explicit restriction on the time interval or smallness/regularity assumptions on the data. The commutator estimates in §4 (controlling the difference between the empirical measure and the Vlasov component) close only when ||∇u||_∞ remains bounded, but global smooth solutions to the 3D Navier-Stokes equations are not known to exist for arbitrary initial data. The theorem must therefore be qualified as local-in-time or conditional on small data to support the claimed global rate.

    Authors: We agree that the commutator estimates require ||∇u||_∞ to remain bounded and that global smooth solutions to the 3D Navier-Stokes equations are not known to exist for arbitrary data. The manuscript derives the rates under the standing assumption that a sufficiently regular solution to the coupled Vlasov-Navier-Stokes system exists on the time interval under consideration. To eliminate any ambiguity, we will revise Theorem 2.1 to state explicitly that the result holds locally in time or conditionally on the existence of a smooth solution with bounded velocity gradient. This qualification will also be noted in the introduction. revision: yes

  2. Referee: [§3] §3 (model assumptions and initial data): the precise scaling of the noise term (vanishing as N→∞) and the form of the drag force are load-bearing for the energy estimates to close, yet the precise Sobolev or moment assumptions on the initial particle distribution and fluid velocity are not listed explicitly in the theorem statement. These must be stated verbatim so that the reader can verify the estimates apply.

    Authors: We thank the referee for this observation. The vanishing scaling of the noise with N and the precise form of the drag force are indeed central to closing the estimates and are fully specified in Section 3. To improve immediate verifiability, we will add a concise, verbatim summary of the required Sobolev regularity, moment bounds on the initial particle distribution, and fluid velocity assumptions directly into the statement of Theorem 2.1, with an explicit reference to the drag force as given in the model equations. revision: yes

Circularity Check

0 steps flagged

No circularity; standard energy/commutator estimates for particle-to-PDE convergence

full rationale

The paper derives quantitative convergence rates for the empirical measure and fluid velocity using stochastic calculus combined with energy estimates and commutator estimates on both the discrete particle system and the continuous Vlasov(-Fokker-Planck)-Navier-Stokes limit. These steps are constructed directly from the model equations and standard PDE techniques without any self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations that presuppose the target result. The derivation chain remains independent of the claimed rates and does not reduce to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on the abstract only; no free parameters, axioms, or invented entities are mentioned. The result rests on standard frameworks from stochastic analysis and PDE theory.

pith-pipeline@v0.9.0 · 5498 in / 1186 out tokens · 50064 ms · 2026-05-10T03:11:19.280583+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

34 extracted references · 34 canonical work pages

  1. [1]

    Grégoire Allaire,Homogenization of the Navier-Stokes equations in open sets perforated with tiny holes. I. Abstract framework, a volume distribution of holes, Arch. Rational Mech. Anal.113 (1990), no. 3, 209–259,doi. MR-1079189

    work page 1990

  2. [2]

    343, Springer, Heidelberg, 2011,doi

    Hajer Bahouri, Jean-Yves Chemin, and Raphaël Danchin,Fourier analysis and nonlinear partial differential equations, Grundlehren der mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences], vol. 343, Springer, Heidelberg, 2011,doi. MR-2768550

    work page 2011

  3. [3]

    Baranger, L

    C. Baranger, L. Boudin, P.-E. Jabin, and S. Mancini,A modeling of biospray for the upper airways, CEMRACS 2004—mathematics and applications to biology and medicine, ESAIM Proc., vol. 14, EDP Sci., Les Ulis, 2005, pp. 41–47. MR-2226800

    work page 2004

  4. [4]

    Etienne Bernard, Laurent Desvillettes, François Golse, and Valeria Ricci,A derivation of the Vlasov-Navier-Stokes model for aerosol flows from kinetic theory, Commun. Math. Sci.15(2017), no. 6, 1703–1741,doi. MR-3668954

    work page 2017

  5. [5]

    Integral Equ.22(2009), no

    Laurent Boudin, Laurent Desvillettes, Céline Grandmont, and Ayman Moussa,Global existence of solutions for the coupled Vlasov and Navier-Stokes equations, Differ. Integral Equ.22(2009), no. 11-12, 1247–1271. MR-2555647

    work page 2009

  6. [6]

    Surveys, vol

    Laurent Boudin, Céline Grandmont, Bérénice Grec, Sébastien Martin, Amina Mecherbet, and FrédériqueNoël,Fluid-kinetic modelling for respiratory aerosols with variable size and temperature, CEMRACS 2018—numerical and mathematical modeling for biological and medical applications: deterministic, probabilistic and statistical descriptions, ESAIM Proc. Surveys,...

    work page 2018

  7. [7]

    Laurent Boudin, Céline Grandmont, Alexander Lorz, and Ayman Moussa,Modelling and numerics for respiratory aerosols, Commun. Comput. Phys.18(2015), no. 3, 723–756,doi. MR-3396231

    work page 2015

  8. [8]

    Kleber Carrapatoso and Matthieu Hillairet,On the derivation of a Stokes-Brinkman problem from Stokes equations around a random array of moving spheres, Comm. Math. Phys.373(2020), no. 1, 265–325,doi. MR-4050096

    work page 2020

  9. [9]

    Differential Equations251(2011), no

    MyeongjuChae, KyungkeunKang, andJihoonLee,Global existence of weak and classical solutions for the Navier-Stokes-Vlasov-Fokker-Planck equations, J. Differential Equations251(2011), no. 9, 2431–2465,doi. MR-2825335

    work page 2011

  10. [10]

    9, 3309–3336,doi

    Young-Pil Choi and Bongsuk Kwon,Global well-posedness and large-time behavior for the inho- mogeneous Vlasov-Navier-Stokes equations, Nonlinearity28(2015), no. 9, 3309–3336,doi. MR- 3403400

    work page 2015

  11. [11]

    Laurent Desvillettes, François Golse, and Valeria Ricci,The mean-field limit for solid particles in a Navier-Stokes flow, J. Stat. Phys.131(2008), no. 5, 941–967,doi. MR-2398959

    work page 2008

  12. [12]

    E. B. Fabes, B. F. Jones, and N. M. Rivière,The initial value problem for the Navier-Stokes equations with data inLp, Arch. Rational Mech. Anal.45(1972), 222–240,doi. MR-316915

    work page 1972

  13. [13]

    Franco Flandoli, Marta Leocata, and Cristiano Ricci,The Navier–Stokes–Vlasov–Fokker–Planck system as a scaling limit of particles in a fluid, J. Math. Fluid Mech.23(2021), no. 2, Paper No. 40, 39,doi. MR-4234291 56

    work page 2021

  14. [14]

    MR-3157048

    Isabelle Gallagher, Laure Saint-Raymond, and Benjamin Texier,From Newton to Boltzmann: hard spheres and short-range potentials, Also available as arXiv:1208.5753, Zurich Lectures in Advanced Mathematics, European Mathematical Society (EMS), Zürich, 2013. MR-3157048

    work page arXiv 2013

  15. [15]

    Corcoran, and Norman Chigier,A numerical and experimental study of spray dynamics in a simple throat model, Aerosol Science & Technology36(2002), no

    Tevfik Gemci, Timothy E. Corcoran, and Norman Chigier,A numerical and experimental study of spray dynamics in a simple throat model, Aerosol Science & Technology36(2002), no. 1, 18–38

    work page 2002

  16. [16]

    Thierry Goudon, Lingbing He, Ayman Moussa, and Ping Zhang,The Navier-Stokes-Vlasov- Fokker-Planck system near equilibrium, SIAM J. Math. Anal.42(2010), no. 5, 2177–2202,doi. MR-2729436

    work page 2010

  17. [17]

    Thierry Goudon, Pierre-Emmanuel Jabin, and Alexis Vasseur,Hydrodynamic limit for the Vlasov- Navier-Stokes equations. I. Light particles regime, Indiana Univ. Math. J.53(2004), no. 6, 1495– 1515,doi. MR-2106333

    work page 2004

  18. [18]

    ,Hydrodynamic limit for the Vlasov-Navier-Stokes equations. II. Fine particles regime, Indiana Univ. Math. J.53(2004), no. 6, 1517–1536,doi. MR-2106334

    work page 2004

  19. [19]

    Daniel Han-Kwan and David Michel,On hydrodynamic limits of the Vlasov-Navier-Stokes system, Mem. Amer. Math. Soc.302(2024), no. 1516, v+115,doi. MR-4813036

    work page 2024

  20. [20]

    Daniel Han-Kwan, Évelyne Miot, Ayman Moussa, and Iván Moyano,Uniqueness of the solution to the 2D Vlasov-Navier-Stokes system, Rev. Mat. Iberoam.36(2020), no. 1, 37–60,doi. MR- 4061981

    work page 2020

  21. [21]

    840, Springer-Verlag, Berlin-New York, 1981

    Daniel Henry,Geometric theory of semilinear parabolic equations, Lecture Notes in Mathematics, vol. 840, Springer-Verlag, Berlin-New York, 1981. MR-610244

    work page 1981

  22. [22]

    Math.214(2018), no

    Pierre-Emmanuel Jabin and Zhenfu Wang,Quantitative estimates of propagation of chaos for stochastic systems withW −1,∞ kernels, Invent. Math.214(2018), no. 1, 523–591,doi. MR- 3858403

    work page 2018

  23. [23]

    Lanford, III,Time evolution of large classical systems, Dynamical systems, theory and applications (Rencontres, Battelle Res

    Oscar E. Lanford, III,Time evolution of large classical systems, Dynamical systems, theory and applications (Rencontres, Battelle Res. Inst., Seattle, Wash., 1974), Lecture Notes in Phys., Vol. 38, Springer, Berlin-New York, 1975, pp. 1–111,doi. MR-479206

    work page 1974

  24. [24]

    Theory Related Fields74(1987), no

    Karl Oelschläger,A fluctuation theorem for moderately interacting diffusion processes, Probab. Theory Related Fields74(1987), no. 4, 591–616. MR-876258

    work page 1987

  25. [25]

    Rational Mech

    ,On the connection between Hamiltonian many-particle systems and the hydrodynamical equations, Arch. Rational Mech. Anal.115(1991), no. 4, 297–310,doi. MR-1120850

    work page 1991

  26. [26]

    Christian Olivera, Alexandre Richard, and Milica Tomašević,Quantitative particle approximation of nonlinear Fokker-Planck equations with singular kernel, Ann. Sc. Norm. Super. Pisa Cl. Sci. (5)24(2023), no. 2, 691–749,doi. MR-4630042

    work page 2023

  27. [27]

    1, Paper No

    ,Quantitative approximation of the Burgers and Keller-Segel equations by moderately in- teracting particles, Potential Anal.64(2026), no. 1, Paper No. 22, 31,doi. MR-5000960

    work page 2026

  28. [28]

    Cristyan Pinheiro and Gabriela Planas,On theα-Navier-Stokes-Vlasov and theα-Navier-Stokes- Vlasov-Fokker-Planck equations, J. Math. Phys.62(2021), no. 3, Paper No. 031507, 22,doi. MR-4224710 57

    work page 2021

  29. [29]

    25, Springer, Cham, 2021, pp

    Mario Pulvirenti and Sergio Simonella,A brief introduction to the scaling limits and effective equations in kinetic theory, Trails in kinetic theory–foundational aspects and numerical methods, SEMA SIMAI Springer Ser., vol. 25, Springer, Cham, 2021, pp. 183–208,doi. MR-4293996

    work page 2021

  30. [30]

    Robinson, José L

    James C. Robinson, José L. Rodrigo, and Witold Sadowski,The three-dimensional Navier-Stokes equations, Classical theory, Cambridge Studies in Advanced Mathematics, vol. 157, Cambridge University Press, Cambridge, 2016,doi. MR-3616490

    work page 2016

  31. [31]

    1971, Springer-Verlag, Berlin, 2009,doi

    Laure Saint-Raymond,Hydrodynamic limits of the Boltzmann equation, Lecture Notes in Mathe- matics, vol. 1971, Springer-Verlag, Berlin, 2009,doi. MR-2683475

    work page 1971

  32. [32]

    MR-891189

    Hans-Jürgen Schmeisser and Hans Triebel,Topics in Fourier analysis and function spaces, A Wiley-Interscience Publication, John Wiley & Sons, Ltd., Chichester, 1987. MR-891189

    work page 1987

  33. [33]

    J.169(2020), no

    Sylvia Serfaty,Mean field limit for Coulomb-type flows, Duke Math. J.169(2020), no. 15, 2887– 2935,doi. MR-4158670

    work page 2020

  34. [34]

    Cheng Yu,Global weak solutions to the incompressible Navier-Stokes-Vlasov equations, J. Math. Pures Appl. (9)100(2013), no. 2, 275–293,doi. MR-3073216 58

    work page 2013