An efficient stochastic particle method for high-dimensional nonlinear PDEs

Abstract

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.

References in corpus (14)

• Solving high-dimensional partial differential equations using deep learning
• Deep learning-based numerical methods for high-dimensional parabolic partial differential equations and backward stochastic differential equations
• Ten equivalent definitions of the fractional Laplace operator
• On multilevel Picard numerical approximations for high-dimensional nonlinear parabolic partial differential equations and high-dimensional nonlinear backward stochastic differential equations
• Machine Learning from a Continuous Viewpoint
• A Randomized Algorithm Based on Threshold Accepting to Approximate the Star Discrepancy
• High-Dimensional Stochastic Design Optimization by Adaptive-Sparse Polynomial Dimensional Decomposition
• Stable interpolation with isotropic and anisotropic Gaussians using Hermite generating function
• A massively parallel semi-Lagrangian solver for the six-dimensional Vlasov-Poisson equation
• Exponential methods for solving hyperbolic problems with application to kinetic equations
• The Wigner branching random walk: Efficient implementation and performance evaluation
• Overcoming the numerical sign problem in the Wigner dynamics via adaptive particle annihilation
• Branching random walk solutions to the Wigner equation
• From continuous-time formulations to discretization schemes: tensor trains and robust regression for BSDEs and parabolic PDEs