2025/07/20 by Liet Vo, Vo, Liet
Biochemistry, Genetics and Molecular Biology · Mathematics · #FOS: Mathematics #Gene Regulatory Network Analysis #Mathematical Biology Tumor Growth #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.2507.15103
openalex publication_date 2025/07/20 · openalex created_date 2025/10/16 · openalex updated_date 2026/07/28
We develop and analyze numerical methods for a stochastic Keller-Segel system perturbed by Stratonovich noise, which models chemotactic behavior under randomly fluctuating environmental conditions. The proposed fully discrete scheme couples a Crank-Nicolson time discretization with a splitting mixed finite element method in space. We rigorously prove the stability of the numerical scheme and establish strong convergence rates of order O(k1/2 + k-1/2h2), where k and h denote the time and spatial step sizes, respectively. Notably, the presence of stochastic forcing leads to an inverse dependence on k in the error estimates, distinguishing the convergence behavior from that of the deterministic case. Numerical experiments are presented to validate the theoretical results and demonstrate the effectiveness and accuracy of the proposed methods.