Abstract
In this paper, we present a method to improve the accuracy of the charge simulation method for numerical conformal mapping. The method constructs the constraint equation by using the charge simulation method. The charges and the conformal radius are computed by using the Runge-Kutta method based on the dynamic system of the constraint equation. By using this method, we can obtain new approximated conformal mapping function and improve the accuracy of numerical conformal mapping compared to the charge simulation method for numerical conformal mapping proposed by Amano. Furthermore, the corresponding numerical results are shown to illustrate the performance of the proposed method.
1. Introduction
Conformal mapping appears in complex function theory, which plays important roles for applications in fluid mechanics, image processing, plane elasticity theory, and so on [1, 2]. As is well known, solutions to conformal mapping can be divided into analytical and numerical calculation method. The theoretical basis of analytical method is the famous Riemann Theorem, which proves the existence of conformal mapping function but cannot give the expression for a specific function. This implies that we must use numerical calculation method to solve some related problems when the conformal mapping functions are complex and diverse. Symm [3, 4] proposed the integral equation method for the conformal mapping of a simply connected domain bounded by a closed Jordan curve onto the unit disk, of its exterior onto the exterior of the unit disk and of a doubly connected domain onto a concentric circular annulus. In [5], a numerical integral method is proposed to find the specific length-to-width ratio, which is able to improve the accuracy of the numerical conformal mapping on the boundaries and makes the mapping method more reliable when solving problems that are sensitive to accuracy. Nasser [6–9] proposed a method which can be used to compute the mapping function onto the five canonical regions in a unified way. This method is based on a boundary integral equation with the generalized Neumann kernel. In the year 1969, Steinberger firstly proposed the charge simulation method to study the electric field calculation [10]. Following the work of Steinberger, Amano [11–13] developed the charge simulation methods and obtained many nice and important results on the charge simulation method and numerical conformal mapping. In fact, compared to traditional methods, the charge simulation method has a higher accuracy evaluation.
The matrix of Nasser method is well-conditioned, but the coefficient matrix of the constraint equation of Amano method is ill-conditioned [8]. It is very important to solve the constraint equation in order to obtain the high accuracy charges and the conformal mapping radius. Wilkinson’s iteration refinement for linear systems [14] is well known for solution in which Cholesky decomposition is repeatedly used. It has been shown that if the linear equations are not too ill-conditioned, the iterative sequence of the approximate solution will be convergent to the solution. But this method may not be useful for the too ill-conditioned system of linear equations. Reference [15] established a new iteration improvement of the solution based on the dynamic system of ordinary differential equations (ODEs), and the new method is more effective than Wilkinson’s method. However, the number of iterations required for this method is often more. This paper applies the Runge-Kutta method for the constraint equations of doubly connected domains; the new approach reduced the number of iterations and improved the accuracy.
This paper is organized as follows. In Section 2, this paper describes the numerical conformal mapping of doubly connected domains based on the charge simulation method. In Section 3, a novel algorithm based on Runge-Kutta method for numerical conformal mapping of doubly connected domains is proposed. In Section 4, it gives some numerical examples and shows the superiority of developed method stated. Section 5 is devoted to the statement of conclusions.
2. Numerical Conformal Mapping of Doubly Connected Domains Based on the Charge Simulation Method
This section will describe the conformal mapping of a finite doubly connected domain bounded by two closed Jordan curves and given in the -plane onto a circular annulus in the -plane (see Figure 1) [13, 16–19]. The origin lies inside the domain bounded by . Then, the mapping function can be expressed asHere and are conjugate harmonic functions in and is the solution of the Dirichlet problem
Based on the charge simulation method, and can be approximated byrespectively. The charge points are arranged outside the given domain . More precisely, charge points are arranged outside the domain bounded by and inside the domain bounded by . The charges can be determined to satisfy the Dirichlet boundary condition of collection points arranged on the boundary and ; that is,According to the condition of existence of conjugate harmonic function, we have
From (5) and (6), we obtain the constraint equationsof doubly connected domains based on the charge simulation method. Here .
3. Runge-Kutta Method for Numerical Conformal Mapping of Doubly Connected Domain
It is important to find an effective way to solve the constraint equation (7). We may consider linear equations of the formwhere and and represents the number of simulation charge points.
Normally, the higher precision of conformal mapping, the more points of simulation charge points. In this case, the condition number of coefficient matrix is large. Hence (8) is ill-conditioned linear systems.
Consider the following preconditioned linear equations:Eq. (9) is on both sides plus , where . Thus, we can construct iterative formulaFor the simplicity of expression, we define , , and . Assume that , thenThis can be viewed as explicit Euler method [15, 20] with step for solving the system of ordinary differential equations
Runge-Kutta method is a high-precision single-step algorithm for numerical solutions of ordinary differential equations [21–24]. For integral differential equation (12), from to , we obtain
To improve the formula order, it is necessary to increase the accuracy of the quadrature formula of the right-hand. Obviously, it is necessary to increase the quadrature node, for which the right-hand quadrature formula of (14) can be expressed as
In general, the more the nodes , the higher the precision. In order to obtain the explicit method that facilitates calculating, formula (15) is expressed as in (16) were given bywhere , , , and are constants given by [22].
Here we prove the convergence of the algorithm described above, taking fourth-order Runge-Kutta method as an example, and it is easy to prove other Runge-Kutta methods converge.
Theorem 1. For the initial value problems for systems of ordinary differential equations When the function on its domain about satisfies Lipschitz condition, the increment function of the fourth-order Runge-Kutta method satisfies Lipschitz condition on , where , , and .
Proof. For any and we obtainSince on its domain about satisfies Lipschitz condition, we have From this we obtain where . Since and are arbitrary, therefore satisfies Lipschitz condition on .
Theorem 2. For the initial value problems for systems of ordinary differential equations When the function on its domain about satisfies Lipschitz condition, the fourth-order Runge-Kutta method convergence.
Proof. Note that ; then Theorem 1 implies fourth-order Runge-Kutta method convergence.
At this point, we can obtain the Runge-Kutta method for solving the constraint equations. Specific steps are given in Algorithm 1, where and represent the relative error and maximum number of iterations, respectively. Letting , , , and are determined by the constraint equations given in Section 4.
| ||||||||||||||||||||
According to the previous analysis, the Runge-Kutta method for numerical conformal mapping based on the charge simulation method is summarized as follows.(1)Given the number of charge points , charge points and constraint points .(2)By Algorithm 1, we calculate charges and the approximation of transformation radius.(3)After calculating and by (3) and (4) for each point in the closed region , including the boundary , substitute into the approximate conformal mapping function to calculate the corresponding points.
4. Numerical Examples
In this section, two numerical experiments are made to show the effectiveness of Runge-Kutta method in comparison with Amano’s method [13]. In the first example, we consider a doubly connected domain with exact solution given by [25]. In the second example, we consider the conformal mapping of a finite concentration ellipse onto a circular annulus. In the third example, we consider the conformal mapping of a finite doubly connected domain bounded by two squares onto a circular annulus. According to [13], the error formula of numerical conformal mapping of doubly connected domain isThe Runge-Kutta method for numerical conformal mapping of doubly connected domain has been implemented in MATLAB.
Example 1. Region is the doubly connected domain with the boundaries and , where (see Figure 2).
In this example, the analytic solution for domain is known, which is given by
The values of parameters , , and in this example are taken as , , and , respectively. Figure 2 shows the distribution of charge points for , half arranged outside the domain bounded by and half inside the domain bounded by . A comparison between the exact values of the conformal mapping and the approximate values, obtained by Amano’s method and by our method with the same number of collection points, at various points on the boundaries and , is given in Table 1. From Figure 3 and Table 1, we can see that our method can achieve higher numerical conformal mapping accuracy. The error of numerical conformal mapping is computed by (23). The approximate conformal mapping function maps Figures 4 and 5, from which we can see the effectiveness of our method. From Table 2, we can observe that the CPU time of our method is shorter than Amano’s method.
Example 2. Region is the doubly connected domain with the boundaries and (see Figure 6).
For this example, let ; the relative error and maximum number of iterations are and . Figure 6 shows the distribution of charge points for , half arranged outside the domain bounded by and half inside the domain bounded by . In Figure 3, we consider the numerical conformal mapping error behavior based on the Runge-Kutta methods when changing the number of charge points. Looking at Figure 7, we can observe that numerical conformal mapping error decreases by changing the number of charge points. The accuracy of numerical conformal mapping based on the Runge-Kutta is higher than Amano’s method when the number of charge points is more than . Some numerical results of CPU time are listed in Table 3; we can see that the CPU time of Amano’s method is longer than our method. When the number of charge points , the accuracy of numerical conformal mapping based on the Runge-Kutta method is , whereas for the Amano’s method the accuracy is .
By applying the Runge-Kutta method to the constrains equation, the approximate conformal mapping function is constructed by using new charge points and conformal mapping radius. In Figure 8 thick solid lines represent the contour boundaries; the thin solid line represents the contour lines. As can be seen from Figure 9, the approximate conformal mapping function maps the boundaries to concentric circles, and the inner region still corresponds to the interior region of concentric circles. This proves the effectiveness of the proposed method.
Example 3. Region is the doubly connected domain with the boundaries , and , (see Figure 10).
For this example, the values of parameters , , and are taken as , , and , respectively. Similar to Figure 6, Figure 10 shows the distribution of the charge points for . For Figure 11, we can find that the accuracy of numerical conformal mapping based on the Runge-Kutta method is better than that of Amano’s method. When the number of charge points is more than , the error of Amano’s method is obviously higher than that of our method. The error of numerical conformal mapping based on the Runge-Kutta method has reached when the number of charge simulation points . Figures 12 and 13 represent the contour lines and corresponding mapping, respectively. Figure 13 shows that the method used in this paper is effective. As can be seen from Table 4, the CPU time of our method is shorter than that of Amano’s method.
5. Conclusion
In this paper, a method has been proposed for improving the accuracy of the charge simulation method for numerical conformal mapping of doubly connected domain onto a circular annulus. The method constructs constraint equation by using the charge simulation method by Amano, and the constraint equations are solved by Runge-Kutta method. Numerical experiments showed that the proposed method can obtain higher numerical conformal mapping accuracy than Amano’s method. The results of numerical conformal mapping of doubly connected domain are simulated by using the contour lines. In a forthcoming work, this method is considered for computing conformal mapping of multiply connected domains.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (NSFC) (Grant no. 11461037).