A new formulation for the numerical proof of the existence of solutions to elliptic problems

10/02/2019
by   Kouta Sekine, et al.
0

Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the infinite Newton-type fixed point equation w = - L^-1 F(û) + L^-1 G(w), where L is a linearized operator, F(û) is a residual, and G(w) is a local Lipschitz term. Therefore, the estimations of L^-1 F(û) and L^-1 G(w) play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement' for matrix problems, we represent the inverse operator L^-1 as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as L^-1 are presented in the appendix.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/16/2023

A score-based operator Newton method for measure transport

Transportation of probability measures underlies many core tasks in stat...
research
01/03/2021

Validated forward integration scheme for parabolic PDEs via Chebyshev series

In this paper we introduce a new approach to compute rigorously solution...
research
01/03/2022

Variational Inverting Network for Statistical Inverse Problems of Partial Differential Equations

To quantify uncertainties of the inverse problems governed by partial di...
research
11/14/2021

Numerical methods to evaluate Koopman matrix from system equations

A method that is employed to evaluate a Koopman matrix from a data set o...
research
07/19/2020

An Efficient Online-Offline Method for Elliptic Homogenization Problems

We present a new numerical method for solving the elliptic homogenizatio...
research
03/25/2021

A Matrix-free Multigrid Preconditioner for Jacobian-free Newton-Krylov Methods

In this work, we propose a multigrid preconditioner for Jacobian-free Ne...

Please sign up or login with your details

Forgot password? Click here to reset