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.

Variational formulation based on duality to solve partial differential equations: Use of B-splines and machine learning approximants

Hidden convexity in the heat, linear transport, and Euler's rigid body equations: A computational approach

Robust SDE-Based Variational Formulations for Solving Linear PDEs via Deep Learning

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

Gauss Newton Method for Solving Variational Problems of PDEs with Neural Network Discretizaitons

New Generalized Derivatives for Solving Variational Inequalities Using the Nonsmooth Newton Methods

Variationally Correct Neural Residual Regression for Parametric PDEs: On the Viability of Controlled Accuracy

Variational Monte Carlo - Bridging Concepts of Machine Learning and High Dimensional Partial Differential Equations

A Hidden Convexity of Nonlinear Elasticity

Variational Iteration Method for Partial Differential Equations with Piecewise Constant Arguments

A Variational Perspective on High-Resolution ODEs

Variational framework for distance-minimizing method in data-driven computational mechanics

Algorithms based on DQM with new sets of base functions for solving parabolic partial differential equations in $(2+1)$ dimension

Learning solutions to partial differential equations using LS-SVM

Variational Principles for Some Nonlinear Partial Differential Equations with Variable Coefficients

Meshless Method Based on Radial Basis Functions for Solving Parabolic Partial Differential Equations with Variable Coefficients

Variational problems and partial differential equations on implicit surfaces

Numerical Solution and Bifurcation Analysis of Nonlinear Partial Differential Equations with Extreme Learning Machines

Machine Learning Approximation Algorithms for High-Dimensional Fully Nonlinear Partial Differential Equations and Second-order Backward Stochastic Differential Equations

Biquadratic element discrete duality finite volume method for solving elliptic equations on quadrilateral mesh