DeepAI AI Chat
Log In Sign Up

On the discretization of Laplace's equation with Neumann boundary conditions on polygonal domains

by   Jeremy Hoskins, et al.
Simons Foundation
Yale University

In the present paper we describe a class of algorithms for the solution of Laplace's equation on polygonal domains with Neumann boundary conditions. It is well known that in such cases the solutions have singularities near the corners which poses a challenge for many existing methods. If the boundary data is smooth on each edge of the polygon, then in the vicinity of each corner the solution to the corresponding boundary integral equation has an expansion in terms of certain (analytically available) singular powers. Using the known behavior of the solution, universal discretizations have been constructed for the solution of the Dirichlet problem. However, the leading order behavior of solutions to the Neumann problem is O(t^μ) for μ∈ (-1/2,0) depending on the angle at the corner (compared to O(C+t^μ) with μ>1/2 for the Dirichlet problem); this presents a significant challenge in the design of universal discretizations. Our approach is based on using the discretization for the Dirichlet problem in order to compute a solution in the "weak sense" by solving an adjoint linear system; namely, it can be used to compute inner products with smooth functions accurately, but it cannot be interpolated. Furthermore we present a procedure to obtain accurate solutions arbitrarily close to the corner, by solving a sequence of small local subproblems in the vicinity of that corner. The results are illustrated with several numerical examples.


An Integral Equation Method for the Cahn-Hilliard Equation in the Wetting Problem

We present an integral equation approach to solve the Cahn-Hilliard equa...

A Steklov-spectral approach for solutions of Dirichlet and Robin boundary value problems

In this paper we revisit an approach pioneered by Auchmuty to approximat...

Interaction with an obstacle in the 2d focusing nonlinear Schrödinger equation

We present a numerical study of solutions to the 2d focusing nonlinear S...

Laplace Green's Functions for Infinite Ground Planes with Local Roughness

The Green's functions for the Laplace equation respectively satisfying t...

Regularity of the solution of the scalar Signorini problem in polygonal domains

The Signorini problem for the Laplace operator is considered in a genera...

Mathematical aspects relative to the fluid statics of a self-gravitating perfect-gas isothermal sphere

In the present paper we analyze and discuss some mathematical aspects of...

Finding twin smooth integers by solving Pell equations

Any pair of consecutive B-smooth integers for a given smoothness bound B...