A Stochastic Genetic Interacting Particle Method for Reaction-Diffusion-Advection Equations (2511.12275v1)

Abstract: We develop and analyze a stochastic genetic interacting particle method (SGIP) for reaction-diffusion-advection (RDA) equations. The method employs operator splitting to approximate the advection-diffusion and reaction processes, treating the former by particle drift-diffusion and the latter by exact or implicit integration of reaction dynamics over bins where particle density is estimated by a histogram. A key innovation is the incorporation of adaptive resampling to close the loop of particle and density field description of solutions, mimicking mutation in genetics. Resampling is also critical for maintaining long-term stability by redistributing particles with the evolving density field. We provide a comprehensive error analysis, and establish rigorous convergence bounds under appropriate regularity assumptions. Numerical experiments from one to three space dimensions demonstrate the method's effectiveness across various reaction types (FKPP, cubic, Arrhenius) and flow configurations (shear, cellular, cat's eye, ABC flows), showing excellent agreement with finite difference method (FDM) while offering computational advantages for complex flow geometries and higher-dimensional problems.