Abstract:Nonlinear partial differential equations (PDEs) are used to model dynamical processes in a large number of scientific fields, ranging from finance to biology. In many applications standard local models are not sufficient to accurately account for certain non-local phenomena such as, e.g., interactions at a distance. Non-local nonlinear PDE models can accurately capture these phenomena, but traditional numerical approximation methods are infeasible when the considered non-local PDE is high-dimensional. In this article we propose two numerical methods based on machine learning and on Picard iterations, respectively, to approximately solve non-local nonlinear PDEs. The proposed machine learning-based method is an extended variant of a deep learning-based splitting-up type approximation method previously introduced in the literature and utilizes neural networks to provide approximate solutions on a subset of the spatial domain of the solution. The Picard iterations-based method is an extended variant of the so-called full history recursive multilevel Picard approximation scheme previously introduced in the literature and provides an approximate solution for a single point of the domain. Both methods are mesh-free and allow non-local nonlinear PDEs with Neumann boundary conditions to be solved in high dimensions. In the two methods, the numerical difficulties arising due to the dimensionality of the PDEs are avoided by (i) using the correspondence between the expected trajectory of reflected stochastic processes and the solution of PDEs (given by the Feynman–Kac formula) and by (ii) using a plain vanilla Monte Carlo integration to handle the non-local term. We evaluate the performance of the two methods on five different PDEs arising in physics and biology. In all cases, the methods yield good results in up to 10 dimensions with short run times. Our work extends recently developed methods to overcome the curse of dimensionality in solving PDEs.

An Extreme Learning Machine-Based Method for Computational PDEs in Higher Dimensions

Local Extreme Learning Machines and Domain Decomposition for Solving Linear and Nonlinear Partial Differential Equations

On Computing the Hyperparameter of Extreme Learning Machines: Algorithm and Application to Computational PDEs, and Comparison with Classical and High-Order Finite Elements

Algorithms for Solving High Dimensional PDEs: from Nonlinear Monte Carlo to Machine Learning

Solving PDEs on unknown manifolds with machine learning

A deep learning solution approach for high-dimensional random differential equations

Solving High-Dimensional Partial Differential Equations Using Deep Learning.

Numerical solution for high-dimensional partial differential equations based on deep learning with residual learning and data-driven learning

A Nonoverlapping Domain Decomposition Method for Extreme Learning Machines: Elliptic Problems

Machine Learning For Elliptic PDEs: Fast Rate Generalization Bound, Neural Scaling Law and Minimax Optimality

Deep learning approximations for non-local nonlinear PDEs with Neumann boundary conditions

Solving High Dimensional Partial Differential Equations Using Tensor Neural Network and A Posteriori Error Estimators

Solving High-Dimensional Partial Differential Equations Using Tensor Neural Network and A Posteriori Error Estimators

Global Convergence of Deep Galerkin and PINNs Methods for Solving Partial Differential Equations

Solving Partial Differential Equations Using Deep Learning and Physical Constraints

Physics Informed Extreme Learning Machine (PIELM) -- A rapid method for the numerical solution of partial differential equations

Tackling the curse of dimensionality with physics-informed neural networks

Extremization to Fine Tune Physics Informed Neural Networks for Solving Boundary Value Problems

Solving two-dimensional linear partial differential equations based on Chebyshev neural network with extreme learning machine algorithm

Neural networks catching up with finite differences in solving partial differential equations in higher dimensions