Abstract
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 RiemannHilbert (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.
Keywords:
mixed boundary value problem; RH problem; Fredholm integral equation; generalized Neumann kernel1 Introduction
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 [46] 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., [79]). 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 twodimensional 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.
2 Notations and auxiliary material
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).
Figure 1. An unbounded multiply connected regionGof connectivitym.
The curves are parametrized by 2πperiodic twice continuously differentiable complexvalued functions with nonvanishing 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 complexvalued 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
3 The mixed boundary value problem
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 singlevalued. The function F can be written as
where f is a singlevalued 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])
4 Integral equation
In this paper we assume that the function A is a continuously differentiable complexvalued 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
5 Reformulation of the mixed boundary value problem as an RH problem
The mixed boundary value problem can be reduced to an RH problem as follows. Let the boundary values of the multivalued analytic function F be given by
Although, the function is in general multivalued, its derivative is a singlevalued 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 CauchyRiemann 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 CauchyRiemann 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
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
6 The solution of the mixed boundary value problem
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
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).
7 Numerical implementations
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].
8 Numerical examples
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)
Figure 2. The region for Example 1.
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^{−15}. 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.
Figure 3. The error norm (54) for Example 1.
Figure 4. The absolute error for Example 1.
Figure 5. The graph of the approximate solution for Example 1.
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.
Figure 7. The error norm (56) for Example 2.
Figure 8. The absolute error for Example 2.
Figure 9. The graph of the approximate solution for Example 2.
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 starshape boundary with Dirichlet data , while on the ellipseshape boundary is the Neumann condition with Neumann data . The graph of the approximate solution is illustrated in Figure 11.
Figure 10. The region for Example 3.
Figure 11. The graph of the approximate solution for Example 3.
9 Conclusion
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.
Figure 12. The eigenvalues of the coefficient matrix of the linear systems obtained by discretizing integral equations (45) withfor the region in Example 1.
Figure 13. The eigenvalues of the coefficient matrix of the linear systems obtained by discretizing integral equations (45) withfor the region in Example 2.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally to the manuscript and typed, read and approved the final manuscript.
Acknowledgements
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.
References

Nasser, MMS, Murid, AHM, AlHatemi, 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)

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

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

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

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

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)

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

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

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

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

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

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

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

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

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

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

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