Accelerating the convergence of Newton's method for nonlinear elliptic PDEs using Fourier neural operators

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Aghili, Joubine, Franck, Emmanuel, Hild, Romain, Michel-Dansac, Victor, Vigon, Vincent
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866915054010499072
author Aghili, Joubine
Franck, Emmanuel
Hild, Romain
Michel-Dansac, Victor
Vigon, Vincent
author_facet Aghili, Joubine
Franck, Emmanuel
Hild, Romain
Michel-Dansac, Victor
Vigon, Vincent
contents It is well known that Newton's method can have trouble converging if the initial guess is too far from the solution. Such a problem particularly occurs when this method is used to solve nonlinear elliptic partial differential equations (PDEs) discretized via finite differences. This work focuses on accelerating Newton's method convergence in this context. We seek to construct a mapping from the parameters of the nonlinear PDE to an approximation of its discrete solution, independently of the mesh resolution. This approximation is then used as an initial guess for Newton's method. To achieve these objectives, we elect to use a Fourier neural operator (FNO). The loss function is the sum of a data term (i.e., the comparison between known solutions and outputs of the FNO) and a physical term (i.e., the residual of the PDE discretization). Numerical results, in one and two dimensions, show that the proposed initial guess accelerates the convergence of Newton's method by a large margin compared to a naive initial guess, especially for highly nonlinear and anisotropic problems, with larger gains on coarse grids.
format Preprint
id arxiv_https___arxiv_org_abs_2403_03021
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Accelerating the convergence of Newton's method for nonlinear elliptic PDEs using Fourier neural operators
Aghili, Joubine
Franck, Emmanuel
Hild, Romain
Michel-Dansac, Victor
Vigon, Vincent
Numerical Analysis
It is well known that Newton's method can have trouble converging if the initial guess is too far from the solution. Such a problem particularly occurs when this method is used to solve nonlinear elliptic partial differential equations (PDEs) discretized via finite differences. This work focuses on accelerating Newton's method convergence in this context. We seek to construct a mapping from the parameters of the nonlinear PDE to an approximation of its discrete solution, independently of the mesh resolution. This approximation is then used as an initial guess for Newton's method. To achieve these objectives, we elect to use a Fourier neural operator (FNO). The loss function is the sum of a data term (i.e., the comparison between known solutions and outputs of the FNO) and a physical term (i.e., the residual of the PDE discretization). Numerical results, in one and two dimensions, show that the proposed initial guess accelerates the convergence of Newton's method by a large margin compared to a naive initial guess, especially for highly nonlinear and anisotropic problems, with larger gains on coarse grids.
title Accelerating the convergence of Newton's method for nonlinear elliptic PDEs using Fourier neural operators
topic Numerical Analysis
url https://arxiv.org/abs/2403.03021