English

A Non-compact Positivity-Preserving Numerical Scheme for Elliptic Differential Equations Based on Mathematical Expectation

Numerical Analysis 2026-04-06 v1 Numerical Analysis

Abstract

We propose a novel non-compact, positivity-preserving scheme for linear non-divergence form elliptic equations. Based on the Feynman--Kac formula, the solution is represented as a conditional expectation associated with a diffusion process.Instead of using compact Markov chain approximations, we construct a wide-stencil scheme by approximating the expectation with carefully designed transition probabilities, ensuring both consistency and positivity preservation. The method is effective for anisotropic diffusion problems with mixed derivatives, where classical schemes typically fail unless the covariance matrix is diagonally dominant. A key feature of the proposed framework is its robust treatment of boundary conditions. For Dirichlet boundaries, we introduce a quadtree-based non-uniform stopping-time strategy, achieving O(h)O(h) accuracy. For Neumann boundaries, a discrete specular reflection mechanism is employed, yielding O(h1/2)O(h^{1/2}) convergence. Periodic boundaries are handled through modular wrapping, also achieving O(h)O(h) accuracy. The resulting schemes are unconditionally stable and positivity-preserving due to their probabilistic structure. Numerical experiments confirm the theoretical convergence rates under all boundary conditions considered.

Keywords

Cite

@article{arxiv.2604.02797,
  title  = {A Non-compact Positivity-Preserving Numerical Scheme for Elliptic Differential Equations Based on Mathematical Expectation},
  author = {Haoran Xu and Kunyang Li and Xingye Yue},
  journal= {arXiv preprint arXiv:2604.02797},
  year   = {2026}
}

Comments

arXiv admin note: text overlap with arXiv:2601.10977

R2 v1 2026-07-01T11:52:28.322Z