The Approximate Solution of 2D Dirichlet Problem in Doubly Connected Domains

We propose a newmethod for constructing an approximate solution of the two-dimensional Laplace equation in an arbitrary doubly connected domain with smooth boundaries for Dirichlet boundary conditions. Using the fact that the solution of the Dirichlet problem in a doubly connected domain is represented as the sum of a solution of the Schwarz problem and a logarithmic function, we reduce the solution of the Schwartz problem to the Fredholm integral equation with respect to the boundary value of the conjugate harmonic function. The solution of the integral equation in its turn is reduced to solving a linear system with respect to the Fourier coefficients of the truncated expansion of the boundary value of the conjugate harmonic function. The unknown coefficient of the logarithmic component of the solution of the Dirichlet problem is determined from the following fact.TheCauchy integral over the boundary of the domainwith a density that is the boundary value of the analytical in this domain function vanishes at points outside the domain. The resulting solution of the Dirichlet problem is the sum of the real part of the Cauchy integral in the given domain and the logarithmic function. In order to avoid singularities of the Cauchy integral at points near the boundary, the solution at these points is replaced by a linear function.The resulting numerical solution is continuous in the domain up to the boundaries.Three examples of the solution of theDirichlet problem are given: one example demonstrates the solutionwith constant boundary conditions in the domain with a complicated boundary; the other examples provide a comparison of the approximate solution with the known exact solution in a noncircular domain.


Introduction
Here we introduce the Cauchy integral method for the solution of the Dirichlet problem in doubly connected domains.The proposed method gives an analytical approximate solution to this problem.This analytical solution is differentiable at the interior points and can be calculated at any point in the domain of solution and this is the main advantage of the Cauchy integral method.The method can be applied to arbitrary simply and multiply connected domains with smooth boundaries.Furthermore, it is applicable for the domains with piecewise smooth boundary curves which can be approximated by Fourier polynomial.This technic can be extended to apply the method to Poisson equation, biharmonic equation, and some other types of problems with different types of boundary conditions (Neumann and mixed).It can be applied for solving Riemann-Hilbert problem for analytical functions.The base of these problems is the Dirichlet problem.Here we introduce the solution of 2D Laplace equation with Dirichlet conditions in order to declare the basic idea of the method.In Section 3 we have to improve the behavior of the solution in the form of the Cauchy integral with the help of the linearization technic due to the singularity of Cauchy integral at the points closed to the boundaries.
The mathematical theory regarding Laplace's equation is often referred to as the potential theory, given the significance the equation holds for describing physical phenomena such as gravitational and electrical potentials.Laplace equation with Dirichlet boundary conditions arises in different areas such as electrostatics (where it describes the electrostatic potential in a charge-free region), gravitation (where it describes the gravitational potential in free space), steady state flow of inviscid fluids, and steady state heat conduction.Many authors (e.g., [1][2][3][4]) introduced integral equation methods for the two-dimensional Laplace equation solution 2 Advances in Mathematical Physics in order to calculate the potential field.The boundary value problem was reduced to an integral equation.The singular kernel in the obtained boundary integral equation was manipulated so that the numerical integration methods or any other numerical methods can be applied directly to evaluate the integrals without any difficulty.
The numerical solution of two-dimensional Laplace equation with Dirichlet boundary conditions in doubly connected domain has been introduced by many authors; for example, the complex variable boundary element methods has been presented in [5].A highly accurate collocation Trefftz method appeared in [6].
The proposed method constructs the approximate Cauchy integral solution of 2D Dirichlet problem in doubly connected domains.The method is based on the reduction of the problem to the Fredholm integral equation of the second kind for the boundary values of the conjugate harmonic function.The singularity of the obtained integral equation is overcome by using the Hilbert formula.The solution of the resulting integral equation is reduced to the solution of a truncated linear system by using the truncated Fourier series.Finally, the solution of Dirichlet problem has the form of the real part of the Cauchy integral.The solution at the points near the boundaries is approximated by linear functions.

2D Dirichlet Problem for the Laplace Equation
Let Ω be a doubly connected domain and Ω =  0 ∪  1 be the smooth boundary of Ω, where  0 is the outer smooth curve and  1 is the inner smooth curve.The corresponding Dirichlet problem for the Laplace equation is to find the doubly differentiable in Ω function (, ), which is continuous in Ω ∪ Ω and satisfies the Laplace equation and the boundary conditions in the form (2)

The Numerical Solution of the 2D Dirichlet
Problem in Doubly Connected Domain Assume that () is analytic in Ω; the solution of the Dirichlet problem according to [7] has the form (, ) = Re(()) +  log ||,  =  + , where  is a real-valued parameter which will be calculated later.So the boundary values of the analytic in Ω function () can be written in the parametric form as follows: Due to [7] the function (, ) is analytic in Ω if and only if where the corresponding singular integral at the right-hand side of the relation is a principle value integral and   (),  = {0,1}, are written in the form of truncated Fourier polynomials: After separating the imaginary part in (4) we obtain the following system of equations: We consider the factor (  −   ) in the expression (  () −   ()),  = 0, 1 in order to separate the improper PV integral in the Fredholm equation of the second kind as follows [8,9]: log so ( 7)-( 8) take the form Here Let us search for the solution of the system of Fredholm integral equations (11) in the form of Fourier series as follows: According to [8,10], (11) has a unique solution if we set the summands The solvability of ( 11) is proved in [8] where the Fourier series solution form of the functions  s (), s = 0, 1, leads us to an infinite linear system of equations which can be reduced to a finite one according to the following lemma.
According to the previous lemma and [10,11] the unique solution of the second kind Fredholm integral equation system (11) can be expressed in the form of truncated Fourier series with arbitrary summands as follows: where  s = (1/2) ∫ 2 0  s () = 0,  = 0, 1.Put relation (15) into (11) and write down the finite linear system in the matrix form as follows [12,13]: Here is the vector of unknown coefficients, are the right-hand side vectors,  is the 4 × 4 identity matrix, and  is a 4 × 4 coefficient block matrix.The elements of these matrices can be easily calculated by writing formulas (11) into the matrix form.
The solution of the linear system of ( 16) is the vector with the elements that are linear functions of the parameter :   () = ĝ () + g  ().Since the point  = 0 is an exterior point of the domain Ω and   () −  log |  ()| + (ĝ  () + g  ()),  = {0, 1}, are the boundary values of the analytic in Ω function, the value of  can be calculated from the following formula: We can satisfy (18) for  = 0, 1, . . ., N, if we apply least square method.Finally, the solution (, ) of Dirichlet problem (1)-( 2) in the doubly connected domain Ω can be expressed as

Improving the Approximate Solution at the Points Near
Boundaries.Here we introduce a new technic to improve the approximate solution at the points near boundaries.
The technic depends on the approximation of the harmonic function solution of Dirichlet problem at these points by linear functions.Let us define the reference curves: where   ,  = 0, 1 are the radii of curvature of the boundaries   (),  = 0, 1,   < 1.Let the approximate solution at the points of reference curves calculated by the Cauchy integral method be   (),  = 0, 1.
Define Ω is the domain with boundaries   ().The solution of Dirichlet problem at the points of Ω\ Ω is approximated with linear functions as follows.
Firstly, for every point , we construct a straight line which is perpendicular to the curve   () and passes through this point.The straight line provides the projections of the point on the curves   (),   () which will be defined as   ( 0 ),   ( 0 ),  = 0, 1.
Secondly, the approximate solution of 2D Dirichlet problem is calculated at this point by substituting in the following linear equation:

Numerical Examples
The Cauchy integral method was applied to the 2D Dirichlet problem and highly accurate results for regular and irregular doubly connected domains with smooth boundaries are obtained.Numerical examples are presented to verify the accuracy of the proposed method in the earlier sections.
Example 1.Let us define the doubly connected domain with nonstarlike boundaries as in [10].The boundaries are composed of the outer and inner parametric curves, respectively, as follows: The method was applied and we use a constant boundary values which equal  0 = 1,  1 = 2 and the contour plot of the solution is shown in Figure 1.The figure shows that the extreme values of the solution are on the boundaries.
The boundary conditions are derived from the closed form of the exact solution (, ) =   cos .The method was applied for a = 4; b = 1 and the value of constant  was calculated and equal to 7.503 − 08.The absolute error along the points of circle with radius  = max ∈[0,2]  1 () is shown in Figure 3.The figure demonstrates the accuracy of the method with absolute error less than 1.2325×10 −5 which is better than the Trefftz method introduced in [6].This example shows the effectivity of our method at the points near the boundaries.
To test the method we consider the exact solution (, ) =  2 −  2 .
Figure 5 shows the absolute error which is calculated along the points of circle with radius equal to 0.6.The maximum error is 1.315E-08 which shows that the proposed method gives highly accurate results for the interior domain when compared with the exact solution.

Conclusions
The Cauchy integral method gives highly accurate results for the solution of 2D Dirichlet problem for irregular doubly connected domains.The proposed algorithm improves the approximate solution at the points near the boundaries.The method is applicable for domains bounded by any smooth curve approximated by Fourier polynomial.Numerical experiments are given to verify the efficiency of the method.

Figure 1 :Figure 2 :
Figure 1: The contour plot of the approximate solution with constant boundary conditions in Example 1.

Figure 3 :
Figure 3: The absolute error along circle with radius max  1 in Example 2.