The collocation method for mixed boundary value problems on domains with curved polygonal boundaries (original) (raw)
Related papers
A numerical technique for linear elliptic partial differential equations in polygonal domains
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2015
Integral representations for the solution of linear elliptic partial differential equations (PDEs) can be obtained using Green's theorem. However, these representations involve both the Dirichlet and the Neumann values on the boundary, and for a well-posed boundary-value problem (BVPs) one of these functions is unknown. A new transform method for solving BVPs for linear and integrable nonlinear PDEs usually referred to as the unified transform ( or the Fokas transform ) was introduced by the second author in the late Nineties. For linear elliptic PDEs, this method can be considered as the analogue of Green's function approach but now it is formulated in the complex Fourier plane instead of the physical plane. It employs two global relations also formulated in the Fourier plane which couple the Dirichlet and the Neumann boundary values. These relations can be used to characterize the unknown boundary values in terms of the given boundary data, yielding an elegant approach for...
Efficient discretization of Laplace boundary integral equations on polygonal domains
Journal of Computational Physics, 2010
We describe a numerical procedure for the construction of quadrature formulae suitable for the efficient discretization of boundary integral equations over very general curve segments. While the procedure has applications to the solution of boundary value problems on a wide class of complicated domains, we concentrate in this paper on a particularly simple case: the rapid solution of boundary value problems for Laplace's equation on two-dimensional polygonal domains. We view this work as the first step toward the efficient solution of boundary value problems on very general singular domains in both two and three dimensions. The performance of the method is illustrated with several numerical examples.
Journal of Computational and Applied Mathematics, 2012
A new and novel approach for analyzing boundary value problems for linear and for integrable nonlinear PDEs was recently introduced. For linear elliptic PDEs, an important aspect of this approach is the characterization of a generalized Dirichlet-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. For this computation, a collocation-type numerical method has been recently developed. Here, we study the collocation's coefficient matrix properties. We prove that, for the Laplace's equation on regular polygon domains with the same type of boundary conditions on each side, the collocation matrix is block circulant, independently of the choice of basis functions. This leads to the deployment of the FFT for the solution of the associated collocation linear system, yielding significant computational savings. Numerical experiments are included to demonstrate the efficiency of the whole computation.
Piecewise Polynomial Collocation for Boundary Integral Equations
SIAM Journal on Scientific Computing, 1995
This paper considers the numerical solution of boundary integral equations of the second kind for Laplace's equation Au 0 on connected regions D in R with boundary S. The bounda S is allowed to be smooth or piecewisc smooth, and we let {AK _< K _< N} be a triangulation of S. The numerical method is collocation with approximations which arc pieccwise quadratic in the parametrization variables, leading to a numerical solution UN. Superconvergence results for UN are given for S a smooth surface and for a special type of refinement strategy for the triangulation. We show that u us is 0(84 log 8) at the collocation node points, with 8 being the mesh size for A K }. Error analyses are given are given for other quantities, and an important error analysis is given for the approximation of S by piecewise quadratic interpolation on each triangular element, with S either smooth or piecewise smooth. The convergence result we prove is only 0(82) but the numerical experiments suggest the result is 0(84) for the error at the collocation points, especially when S is a smooth surface. The numerical integration of the collocation integrals is discussed, and extended numerical examples are given for problems involving both smooth and piecewise smooth surfaces.
Zeitschrift für angewandte Mathematik und Physik ZAMP, 1973
A method is given for solving two dimensional harmonic mixed boundary value problems in simply-connected polygonal domains with re-entrant boundaries. The method consists of a numerical conformal mapping together with three other conformal transformations. The numerical mapping transforms the original domain onto the unit circle, which in turn is mapped onto a rectangle by means of two bilinear and one Schwarz-Christoffel transformations. The transformed problem in the rectangle is solved by inspection.
A domain-type boundary-integral-equation method for two-dimensional biharmonic Dirichlet problem
Engineering Analysis with Boundary Elements, 2006
This paper reports a new boundary-integral-equation method (BIEM) for numerically solving biharmonic problems with Dirichlet boundary conditions. For the solution of these problems in convex polygons, it was found that the accuracy of the conventional BIEM is significantly reduced, and spurious oscillatory behaviour is often observed in the boundary solutions especially for areas near corners (Mai-Duy N, Tanner RI. An effective high order interpolation scheme in BIEM for biharmonic boundary value problems. Eng Anal Bound Elem 2005; 29:210-23). In this study, a new treatment for these difficulties is proposed. The unknown functions in boundary integrals are approximated using a domain-type interpolation scheme rather than traditional boundary-type interpolation schemes. Two test problems are considered to validate the formulation and to demonstrate the attractiveness of the proposed method.
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2011
A recently discovered transform approach allows a large class of PDEs to be solved in terms of boundary and/or contour integrals. We introduce here a spectrally accurate numerical discretization of this approach for the case of Laplace's equation on a polygonal domain, and compare it against an also spectrally accurate implementation of the traditional boundary integral formulation.
Communications in Numerical Methods in Engineering, 2002
The singular function boundary integral method is applied for the solution of a Laplace equation problem over an L-shaped domain. The solution is approximated by the leading terms of the local asymptotic solution expansion, while the Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers. Estimates of great accuracy are obtained for the leading singular coe cients, as well as for the Lagrange multipliers. Comparisons are made with recent numerical results in the literature.