A boundary integral equation with the generalized Neumann kernel for a mixed boundary value problem in unbounded multiply connected regions

In this paper we propose a new method for solving the mixed boundary value problem for the Laplace equation in unbounded multiply connected regions. All simple closed curves making up the boundary are divided into two sets. The Dirichlet condition is given for one set and the Neumann condition is given for the other set. The mixed problem is reformulated in the form of a Riemann-Hilbert (RH) problem which leads to a uniquely solvable Fredholm integral equation of the second kind. Three numerical examples are presented to show the effectiveness of the proposed method.

In the present paper, we continue the research concerned with the study of mixed boundary value problems in the plane started in [1]. We consider a mixed boundary value problem for the Laplace equation in an unbounded multiply connected regions. Recently, the interplay of the RH boundary value problem and integral equations with the generalized Neumann kernel on unbounded multiply connected regions has been investigated in [2]. Based on the reformulations of the Dirichlet problem, the Neumann problem and conformal mappings as RH problems, boundary integral equations with the generalized Neumann kernel have been implemented successfully in [3] to solve the Dirichlet problem and the Neumann problem and in [4-6] to compute the conformal mappings of unbounded multiply connected regions onto the classical canonical slit domains.

The mixed boundary value problem also can be reformulated as an RH problem (see, e.g., [7-9]). Recently, Nasser et al.[1] have presented a uniquely solvable boundary integral equation with the generalized Neumann kernel for solving the mixed boundary value problem in bounded multiply connected regions. The idea of this paper is to reformulate the mixed boundary value problem to the form of the RH problem in unbounded multiply connected regions. Based on this reformulation, we present a new boundary integral equation method for two-dimensional Laplace’s equation with the mixed boundary condition in unbounded multiply connected regions. The method is based on a uniquely solvable boundary integral equation with the generalized Neumann kernel.

This paper is organized as follows. After presenting some auxiliary materials in Section 2, we present in Section 3 the mixed boundary value problem in unbounded multiply connected regions. In Section 4, we give an explanation of an integral equation with the generalized Neumann kernel and its solvability. The reduction of the mixed boundary value problem to the form of the RH problem is given in Section 5. In Section 6, we present the solution of the mixed boundary problem via an integral equation method. In Section 7, we explain briefly the numerical implementation of the method. In Section 8, we illustrate the method by presenting two numerical examples with exact solutions and also one example without an exact solution.

In this section, we review some properties of the generalized Neumann kernel from [2,3,5,10].

We consider an unbounded multiply connected region G of connectivity with boundary consisting of m clockwise oriented smooth closed Jordan curves , . The complement consists of m bounded simply connected components interior to , . We assume , (see Figure 1).

The curves are parametrized by 2π-periodic twice continuously differentiable complex-valued functions with non-vanishing first derivatives, i.e.,

The total parameter domain J is the disjoint union of the intervals , . We define a parametrization of the whole boundary Γ as the complex-valued function η defined on J by

Let H be the space of all real Hölder continuous functions on the boundary Γ. In view of the smoothness of η, a function can be interpreted via , , as a real Hölder continuous 2π-periodic function of the parameter , i.e.,

with real Hölder continuous 2π-periodic functions defined on . So, here and in what follows, we do not distinguish between functions of the form and .

The subspace of H which consists of all piecewise constant functions defined on Γ is denoted by S, i.e., a function has the representation

where are real constants. For simplicity, the function h is denoted by

Let and be two subsets of the set such that

Let n be the exterior unit normal to Γ and let be a given function. We consider the mixed boundary value problem

for a real function u in G. We call (6b) and (6c) Dirichlet conditions and Neumann conditions, respectively.

Problem (6a)-(6c) reduces to the Dirichlet problem for and to the Neumann problem for . Both problems have been considered in [3]. So, we assume in this paper that and .

The mixed boundary value problem (6a)-(6c) is uniquely solvable [11]. Its unique solution u can be regarded as a real part of an analytic function F in G which is not necessary single-valued. The function F can be written as

where f is a single-valued analytic function in G, are fixed points in , ; and are real constants uniquely determined by ϕ (see [12]). Without lost of generality, we assume that . The constants are chosen to ensure that (see [[12], p.149] and [3])

i.e., are given by (see [3])

The constants satisfy

In this paper we assume that the function A is a continuously differentiable complex-valued function given by

where θ is the real piecewise constant function

with either or , . Here the function is different from the ones used in [1,3]. The generalized Neumann kernel formed with A and is defined by

We also define a real kernel M by

The kernel N is continuous and the kernel M has a cotangent singularity type (see [2] for more details). Hence, the operator

is a Fredholm integral operator and the operator

is a singular integral operator.

The solvability of boundary integral equations with the generalized Neumann kernel is determined by the index (the change of the argument of A on the curves divided by 2π) of the function A (see [2]). For the function A given by (10), the indices of A on the curves and the index of A on the whole boundary curve Γ are given by

The generalized Neumann kernel for an integral equation associated with the mixed boundary value problem which will be presented in this paper is different from the generalized Neumann kernel for the integral equation considered in [1,3]. Thus, not all of the properties of the generalized Neumann kernel which have been studied in [3] are valid for the generalized Neumann kernel which will be studied in this paper. For example, it is still true that is not an eigenvalue of the generalized Neumann kernel which means that the presented integral equation is uniquely solvable.

By using the same approach used in [3] for unbounded multiply connected regions, we can prove that the properties of the generalized Neumann kernel proved in [3], except Theorem 8, Theorem 10 and Corollary 2, are still valid for the generalized Neumann kernel formed with the function in (10) above (see [5]).

Thus, we have from [5] the following theorem (see also [2,10]).

Theorem 1For a given function, there exist unique functionsandsuch that

are boundary values of a unique analytic functioninGwith. The functionμis the unique solution of the integral equation

and the functionhis given by

The mixed boundary value problem can be reduced to an RH problem as follows. Let the boundary values of the multi-valued analytic function F be given by

Although, the function is in general multi-valued, its derivative is a single-valued analytic function on G. The boundary values of the function are given by

For the Dirichlet conditions, i.e., and , the functions are equal to the known functions (see (6b)). Thus, the function satisfies the RH problem

The Neumann conditions can also be reduced to an RH problem by using the Cauchy-Riemann equations and integrating along the boundaries , . Let be the unit tangent vector and be the unit external normal vector to Γ at . Let also be the angle between the normal vector and the positive real axis, i.e., . Then

Thus,

Since , then by the Cauchy-Riemann equations, we have

Thus, the function satisfies the RH problem

If we define the real piecewise constant function

the boundary values of the function satisfy on the boundary Γ the condition

where

is known and

The functions for are given by (6b) and (6c). The functions can be then computed for and by integrating the functions . Then it follows from (7), (26) and (27) that the function is a solution of the RH problem

or briefly,

where

for . In view of (8) and (28), the real constants are known for and are given by

However, for , the real constants are unknown. Thus, the boundary condition (29) can be written as

where the function is known and is given by

Obviously, the functions are known explicitly for with . However, for with , it is required to integrate to obtain .

The functions are not necessary 2π-periodic. In order to keep dealing with periodic functions numerically, we do not compute directly by integrating the functions . Instead, we integrate the functions

According to the definitions of the constants and the functions , we have

which implies that the functions are always 2π-periodic. By using the Fourier series for with , the functions can be written as

Then the functions are given for with by

where are undetermined real constants and the functions are given by

Hence, the boundary condition (33) can then be written as

where is the real piecewise constant function

and the function is given by

Let (unknown real constant) and be the analytic function defined on G by

Then is analytic on G with . The function is a solution of the RH problem

where the function is given by (10) and the function is defined by

Let . Then the boundary values of the function are given on the boundary Γ by

where , are knowns and h, μ are unknowns. The real constants are known for and unknown for .

For , let the functions and be the unique solutions of the integral equations

respectively, and be given by

Then it follows from Theorem 1 that

are boundary values of an analytic function . By the uniqueness of the functions h and μ in (44), it follows from (44) and (47) that

and

Equation (49a) with the following equation (from (9)),

represents a linear system of m equations. Since from (43) the function is given by

only the constants c, for and for are unknowns. Thus, linear equations (49a) and (49b) represent a linear system of equations in unknowns for and for .

By obtaining the values of the constants , we obtain the functions μ from (48) and h from (49a). Consequently, the boundary values of the function g are given by

where

The function can be computed for by the Cauchy integral formula. Then the function is computed from

Finally, the solution of the mixed boundary value problem can be computed from , where is given by (7).

Since the functions and are 2π-periodic, the integrals in the operators N and M in integral equations (45) are best discretized on an equidistant grid by the trapezoidal rule [13]. The computational details are similar to previous works in [4,5,10,14]. For analytic integrands, the rate of convergence is better than for any positive integer k (see, e.g., [[15], p.83]). The obtained approximate solutions of the integral equations converge to the exact solutions with a similarly rapid rate of convergence (see, e.g., [[13], p.322]). Since the smoothness of the integrands in (45) depends on the smoothness of the function , results of high accuracy can be obtained for very smooth boundaries.

By using the trapezoidal rule with n (an even positive integer) equidistant collocation points on each boundary component, solving integral equations (45) reduces to solving mn by mn linear systems. Since integral equations (45) are uniquely solvable, then for sufficiently large values of n, the obtained linear systems are also uniquely solvable [13].

In this paper, the linear systems are solved using the Gauss elimination method. By solving the linear systems, we obtain approximations to and for , which give approximations to and for from (46). By solving (49a) and (49b), we get approximations to the constants c, for and for . These give approximations to the boundary values of the function from (50). Then the values of for are calculated by the Cauchy integral formula. For points z which are not close to the boundary Γ, the integrals in the Cauchy integral formula are approximated by the trapezoidal rule. However, for points z near the boundary Γ, the integrand is nearly singular. For the latter case, the integral in the Cauchy integral formula can be calculated accurately using the method suggested in [[16], Eq. (23)]. Then approximate values of the function are computed from (52). Finally, in view of (7), the approximate solution of the mixed boundary value problem can be computed from

In this paper, we have considered regions with smooth boundaries. For some ideas on how to solve numerically boundary integral equations with the generalized Neumann kernel on regions with piecewise smooth boundaries, see [14].

In this section, the proposed method is used to solve three mixed boundary value problems in unbounded multiply connected regions with smooth boundaries.

Example 1 In this example, we consider an unbounded multiply connected region of connectivity 4 bounded by the four circles (see Figure 2)

where .

We assume that the conditions on the boundaries , are the Neumann conditions and the conditions on the boundaries , are the Dirichlet conditions. The functions in (6b)-(6c) are obtained based on choosing an exact solution of the form

We use the error norm

where is the exact solution of the mixed boundary value problem and is the approximate solution obtained with n collocation points. The error norm vs. the total number of calculation points n by using the presented method is shown in Figure 3, where the integral in (54) is discretized by the trapezoidal rule. By using only (256 calculation points on the whole boundary), we obtain error norm less that 10⁻¹⁵. The absolute errors at selected points in the entire domain are plotted in Figure 4. The graph of the approximate solution is illustrated in Figure 5.

Example 2 In this example, we consider an unbounded multiply connected region of connectivity 6 (see Figure 6). The boundary is parametrized by

where the values of the complex constants and the real constants , , are as in Table 1.

The region in this example has been considered in [3,17,18] for the Dirichlet problem and the Neumann problem. In this example, we consider a mixed boundary value problem where we assume that the conditions on the boundaries , , are the Dirichlet conditions and the conditions on the boundaries , , are the Neumann conditions. The functions in (6b)-(6c) are obtained based on choosing an exact solution of the form

where the values of the complex constants are as in Table 1. For the error, we use the error norm (see Figure 7)

The absolute errors at selected points in the entire domain are plotted in Figure 8. The graph of the approximate solution is shown in Figure 9.

Example 3 This example aims to give an impression how the method works for a problem with an unknown exact solution. We assume that the boundaries of an unbounded doubly connected region are represented as follows (see Figure 10):

We assume the Dirichlet condition on the star-shape boundary with Dirichlet data , while on the ellipse-shape boundary is the Neumann condition with Neumann data . The graph of the approximate solution is illustrated in Figure 11.

We have constructed a new boundary integral equation with the generalized Neumann kernel for solving a mixed boundary value problem in unbounded multiply connected regions. The generalized Neumann kernel used in this paper is formed with which is different from the ones used in [1,3]. Numerical evidences show that Theorem 8 in [3], which claims that the eigenvalues of the generalized Neumann kernel lie in , is no longer true for the function of this paper (see Figures 12 and 13). Three numerical examples are presented to illustrate the accuracy of the presented method. The numerical examples illustrate that the proposed method yields approximations of high accuracy. This paper only applies to the explicitly mixed Dirichlet and Neumann boundary conditions in unbounded multiply connected regions, but the method can be extended to a boundary with both mixed boundary conditions in a boundary component . For this case, the function is discontinuous on , where on the part of corresponding to the Dirichlet condition and on the part of corresponding to the Neumann condition. Hence, the RH problem (42) will be a problem with discontinuous coefficient . Thus, further investigations are required. This will be considered in future work.

The authors declare that they have no competing interests.

All authors contributed equally to the manuscript and typed, read and approved the final manuscript.

The authors would like to thank the editor and referees for their helpful comments and suggestions which improved the presentation of the paper. The authors acknowledge the financial support for this research by the Malaysian Ministry of Higher Education (MOHE) through UTM GUP Vote Q.J130000.7126.01H75.

1. Nasser, MMS, Murid, AHM, Al-Hatemi, SAA: A boundary integral equation with the generalized Neumann kernel for a certain class of mixed boundary value problem. J. Appl. Math.. 2012, Article ID 254123. doi:10.1155/2012/254123 (2012)

2. Wegmann, R, Nasser, MMS: The Riemann-Hilbert problem and the generalized Neumann kernel on multiply connected regions. J. Comput. Appl. Math.. 214, 36–57 (2008). Publisher Full Text

3. Nasser, MMS, Murid, AHM, Ismail, M, Alejaily, EMA: Boundary integral equations with the generalized Neumann kernel for Laplace’s equation in multiply connected regions. Appl. Math. Comput.. 217, 4710–4727 (2011). Publisher Full Text

4. Nasser, MMS: Numerical conformal mapping via a boundary integral equation with the generalized Neumann kernel. SIAM J. Sci. Comput.. 31, 1695–1715 (2009). Publisher Full Text

5. Nasser, MMS: Numerical conformal mapping of multiply connected regions onto the second, third and fourth categories of Koebe’s canonical slit domains. J. Math. Anal. Appl.. 382, 47–56 (2011). Publisher Full Text

6. Yunus, AAM, Murid, AHM, Nasser, MMS: Conformal mapping of unbounded multiply connected region onto canonical slit regions. Abstr. Appl. Anal.. 2012, Article ID 293765. doi:10.1155/2012/293765 (2012)

7. Gakhov, FD: Boundary Value Problem, Pergamon, Oxford (1966)

8. Haas, R, Brauchli, H: Fast solver for plane potential problems with mixed boundary conditions. Comput. Methods Appl. Mech. Eng.. 89, 543–556 (1991). Publisher Full Text

9. Haas, R, Brauchli, H: Extracting singularities of Cauchy integrals - a key point of a fast solver for plane potential problems with mixed boundary conditions. J. Comput. Appl. Math.. 44, 167–185 (1992). Publisher Full Text

10. Nasser, MMS: A boundary integral equation for conformal mapping of bounded multiply connected regions. Comput. Methods Funct. Theory. 9, 127–143 (2009)

11. Gonzlez, R, Kress, R: On the Treatment of a Dirichlet-Neumann mixed boundary value problem for harmonic functions by an integral equation method. SIAM. J. Math. Anal.. 8, 504–517 (1977). Publisher Full Text

12. Mikhlin, SG: Integral Equations and Their Applications to Certain Problems in Mechanics Mathematical Physics and Technology, Pergamon, New York (1957)

13. Atkinson, KE: The Numerical Solution of Integral Equations of the Second Kind, Cambridge University Press, Cambridge (1997)

14. Nasser, MMS, Murid, AHM, Zamzamir, Z: A boundary integral method for the Riemann-Hilbert problem in domains with corners. Complex Var. Elliptic Equ.. 53(11), 989–1008 (2008). Publisher Full Text

15. Krommer, AR, Ueberhuber, CW: Numerical Integration on Advanced Computer Systems, Springer, Berlin (1994)

16. Helsing, J, Ojala, R: On the evaluation of layer potentials close to their sources. J. Comput. Phys.. 227, 2899–2921 (2008). Publisher Full Text

17. Greenbaum, A, Greengard, L, McFadden, GB: Laplace’s equation and the Dirichlet-Neumann map in multiply connected domains. J. Comput. Phys.. 105(2), 267–278 (1993). Publisher Full Text

18. Helsing, J, Wadbro, E: Laplace’s equation and the Dirichlet-Neumann map: a new mode for Mikhlin’s method. J. Comput. Phys.. 202, 391–410 (2005). PubMed Abstract | Publisher Full Text