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.

A Small Eddy Correction Method for Nonlinear Dissipative Evolutionary Equations

A two-grid method based on Newton iteration for the Navier-Stokes equations

Enhancing Low-Order Discontinuous Galerkin Methods with Neural Ordinary Differential Equations for Compressible Navier--Stokes Equations

An Eulerian finite element method for the linearized Navier--Stokes problem in an evolving domain

An Eulerian finite element method for the linearized Navier–Stokes problem in an evolving domain

Error estimation for non-linear singularly perturbed reaction–diffusion parabolic problems via element-free Galerkin method

A staggered discontinuous Galerkin method for a class of nonlinear elliptic equations

A Two-Level Method with Backtracking for the Navier--Stokes Equations

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

On a Full Discretisation for Nonlinear Second-Order Evolution Equations with Monotone Damping: Construction, Convergence, and Error Estimates

A narrow-stencil framework for convergent numerical approximations of fully nonlinear second order PDEs

Iterative methods for Navier-Stokes inverse problems

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

Two Robust, Efficient, and optimally Accurate Algorithms for parameterized stochastic navier-stokes Flow Problems

An Oscillation Free Local Discontinuous Galerkin Method for Nonlinear Degenerate Parabolic Equations

Iterative Methods for Navier--Stokes Inverse Problems

High-order exponential time differencing multi-resolution alternative finite difference WENO methods for nonlinear degenerate parabolic equations

An Adaptive Staggered Discontinuous Galerkin Method for the Steady State Convection–Diffusion Equation

A ghost perturbation scheme to solve ordinary differential equations

An Ensemble-Proper Orthogonal Decomposition Method for the Nonstationary Navier-Stokes Equations

A Robust High-Order Residual Distribution Type Scheme for Steady Euler Equations on Unstructured Grids