The generalized Dirichlet-Neumann map for linear elliptic PDEs and its numerical implementation

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:219

Article

Sifalakis, A. G.¹  Fokas, A. S.²  Fulton, S. R.³  Saridakis, Y. G.¹ 

[1] Tech Univ Crete, Dept Sci, Appl Math & Comp Lab, Khania 73100, Greece

[2] Univ Cambridge, Dept Appl Math & Theoret Phys, Cambridge CB3 0WA, England

[3] Clarkson Univ, Dept Math & Comp Sci, Potsdam, NY 13699 USA

关键词: elliptic PDEs;    Dirichlet-Neumann map;    global relation;    collocation;    iterative methods;   

DOI  :  10.1016/j.cam.2007.07.012

来源: Elsevier

【 摘 要 】

A new approach for analyzing boundary value problems for linear and for integrable nonlinear PDEs was introduced in Fokas [A unified transform method for solving linear and certain nonlinear PDEs, Proc. Roy. Soc. London Ser. A 53 (1997) 1411-1443]. For linear elliptic PDEs, an important aspect of this approach is the characterization of a generalized Dirichlet to Neumann map: given the derivative of the solution along a direction of an arbitrary angle to the boundary, the derivative of the Solution perpendicularly to this direction is computed without solving on the interior of the domain. This is based on the analysis of the so-called global relation, an equation which couples known and unknown components of the derivative on the boundary and which is valid for all values of a complex parameter k. A collocation-type numerical method for solving the global relation for the Laplace equation in an arbitrary bounded convex polygon was introduced in Fulton et al. [An analytical method for linear elliptic PDEs and its numerical implementation, J. Comput. Appl. Math. 167 (2004) 465-483]. Here, by choosing a different set of the collocation points (values for k), we present a significant improvement of the results in Fulton et al. [An analytical method for linear elliptic PDEs and its numerical implementation, J. Comput. Appl. Math. 167 (2004) 465-483]. The new collocation points lead to well-conditioned collocation methods. Their combination with sine basis functions leads to a collocation matrix whose diagonal blocks are point diagonal matrices yielding efficient implementation of iterative methods; numerical experimentation suggests quadratic convergence. The choice of Chebyshev basis functions leads to higher order convergence. which for regular polygons appear to be exponential. (C) 2007 Elsevier B.V. All rights reserved.

【 授权许可】

Free