Numerical resolution of high-dimensional nonlinear PDEs remains a huge challenge due to the curse of dimensionality. Starting from the weak formulation of the Lawson-Euler scheme, this paper proposes a stochastic particle method (SPM) by tracking the deterministic motion, random jump, resampling and reweighting of particles. Real-valued weighted particles are adopted by SPM to approximate the high-dimensional solution, which automatically adjusts the point distribution to intimate the relevant feature of the solution. A piecewise constant reconstruction with virtual uniform grid is employed to evaluate the nonlinear terms, which fully exploits the intrinsic adaptive characteristic of SPM. Combining both, SPM can achieve the goal of adaptive sampling in time. Numerical experiments on the 6-D Allen-Cahn equation and the 7-D Hamiltonian-Jacobi-Bellman equation demonstrate the potential of SPM in solving high-dimensional nonlinear PDEs efficiently while maintaining an acceptable accuracy.
翻译:高维非线性偏微分方程的数值求解因维度灾难而仍是一个巨大挑战。本文从 Lawson-Euler 格式的弱形式出发,提出一种随机粒子方法,通过追踪粒子的确定性运动、随机跳跃、重采样与重加权过程来实现。该方法采用实值加权粒子来逼近高维解,其能自动调整点分布以捕捉解的相关特征。通过引入虚拟均匀网格的分段常数重构来评估非线性项,这充分利用了 SPM 固有的自适应特性。两者结合使 SPM 能够实现时间上的自适应采样。在六维 Allen-Cahn 方程和七维 Hamilton-Jacobi-Bellman 方程上的数值实验表明,SPM 在保持可接受精度的同时,具有高效求解高维非线性偏微分方程的潜力。