Research Article  Open Access
Boundary Galerkin Method of a SkewDerivative Problem in the Exterior of an Open Arc Based on Chebyshev Polynomials
Abstract
A problem modeling Hall effect in a semiconductor film from an electrode of arbitrary shape is considered, which is a skewderivative problem. Boundary Galerkin method for solving the problem in Sobolev spaces is developed firstly. The solution is represented in the form of the combined angular potential and singlelayer potential. The final integral equations do not contain hypersingular integrals. Uniqueness and existence of the solution to the equations are proved. The weakly singular and Cauchy singular integral arising in these equations can be computed directly by truncated series of Chebyshev polynomials with their weighting function without approximation. The numerical simulation showing the high accuracy of the scheme is presented.
1. Introduction
Hall effect is the production of a voltage difference across a semiconductor, transverse to an electric current in the conductor and a magnetic field perpendicular to the current. It is widely used in industrial automation, detection, and information processing, which is the basic method for the study of semiconductor material properties. In this paper, we concern an electric current in a semiconductor coming from an electrode, which can be modeled by an open arc. A magnetic is orthogonal to the film. The directions of electric current and electric intensity are not the same because of Hall effect. From a mathematical standpoint, the problem is a skewderivative problem in the exterior of an open arc in a plane.
The boundary integral has been a means of solving problems in the exterior of an open arc. The difficulty of the problem is that the solution is continuous at the tips, but its gradient may have weak singularity. Numerical methods of problems are considered in Hölder space in a series of papers [1–3]. The analysis in [4] is performed within a Sobolev space setting and a weak solution concept, but no numerical method is presented. These authors construct the solution, via a combination of a single and a doublelayer potential, which can not describe the singularity at the tips directly. Thus, another more difficult problem is introduced, that is, hypersingularity.
In [5], the author constructs the solution by singlelayer potential and nonclassical angular potential [6–10] and proves the existence and uniqueness of the solution. An angular potential can be seen as a double layer potential under certain conditions. It has the same order of singularity as a singlelayer potential, which is different from a doublelayer potential. This makes the final integral system much simpler. The solution which is a combination of a singlelayer and angular potential may have weak singularity at the tips of arc if the density function has the form .
In [11], a numerical method for solving this problem is first considered. It numerically solved a Fredholm integral equation of the second kind in Hölder space. The method can avoid the hypersingularity caused from the doublelayer potential. But due to Cauchy singular integral of the system and the special form of density function, approximate calculation of singular integral can not be avoided. In recent years, the inverse problems for arc have been considered by some researchers. In [12–14], the authors obtain asymptotic analysis in terms of the size of the solutions to scattering problems by cracks. Such asymptotics are useful for solving inverse and imaging problems.
In the present paper, we consider the problem in Sobolev spaces firstly. The equations which consist of a Cauchy integral equation of the first kind with additional condition are derived from summing up a singlelayer potential and an angular potential. We discuss the uniqueness and existence of the solution to boundary Galerkin method. The purpose of this paper is to develop an efficient solution method. Due to the properties of Chebyshev polynomials, several researchers [1, 3] use them for solving integral equations. In this paper, we approximate the unknown function by means of truncated series of Chebyshev polynomials with their weighting function. Thus, the numerical solution may have weak singularity at the tips of the arc and the singular integral can be computed directly. The convergence is proved. Some numerical examples are presented in order to test the method by comparison with explicit solutions and the numerical solutions of [11].
The plan of the paper is as follows: in Section 2, model problem is listed. We introduce our formulation of skewderivative problem and proceed by establishing the existence and uniqueness of a solution to integral equations. Then in Section 3, we develop a spectral boundary element method for solving the integral equations numerically and discuss some numerical techniques. Section 4 presents several numerical experiments to illustrate the accuracy.
2. Uniqueness and Existence
In this section, we introduce the model problem firstly. Thus the problem is reduced to an integral equation with some additional conditions. Then we establish the uniqueness and existence of the solution to the integral equations in Sobolev spaces, by inversion of the integral operator.
Consider a constant magnetic field which is perpendicular to a plane semiconductor film . An electrode is placed in the unbounded film. The electrode is modeled by an open curve , . Set By we denote the two end points of . Denote the lefthand of by when increases, while the opposite side is called . The outward normal of is denoted by . The tangent direction of is directed into the direction of increasing. Obviously, , where .
The electricfield potential is the solution of the skewderivative problem as follows: where , is the mobility of the carriers, and is the projection of the magnetic induction onto axis. Suppose is a real constant. The electric current satisfies the equations as follows: where and is a constant. It is easy to see that the direction of current density and the intensity of the electric field do not coincide if . These equations lead to the first formula of (2). The boundary condition on in (2) is deduced by
Problem (2) is considered by [5]. The author establishes the existence of a solution by a combined singlelayer potential and angular potential approach in Hölder space and proves uniqueness up to arbitrary constant.
Consider the angular potential introduced (see [6]): where is defined by From the definition, it is easy to see that an angular potential is a multivalued function. In order to make it singlevalued, it is necessary to require the following additional condition (see [6]): Integrating it by parts, becomes a doublelayer potential as follows: where and .
In [5], the solution of (2) is constructed in the following form: where is an arbitrary constant. is the singlelayer potential If in (11) behaves as , where is a continuous function, the solution has singularity at the endpoints of , that is, with and .
We put (11) in the boundary condition of (2) and get the following integral equation: where is the angle between the vector and the direction of the normal . If in (11) satisfies (14), it is the solution of problem (2).
First, we extend the arc to a piecewise smooth, simple connected, closed curve and define the following spaces (see [15]): where denote the usual Sobolev spaces. That is to say, contains functions such that their extension by zero to is in . In [16], the author characterizes by using the norm where . It is shown that the norms are equivalent.
Now define the dual space of by and the dual space of by . Hence we have the chain Note that also coincides with .
In order to prove the existence and uniqueness of solutions for (8) and (14), we introduce two boundary integral operators which are defined by with . It is easy to see is the limit of normal derivative of angular potential on .
Lemma 1. Under the condition of (8) One rewrites (14) with (8) in operator form
Lemma 2. has a bounded inverse operator.
Proof. Let . It is easy to see and . Let . Obviously, and . There is a homeomorphism from to . Define an operator . Thus is the bounded inverse operator of .
It is known that is positive and bounded below up to a compact perturbation, that is, there exists a compact operator
such that
where denotes the duality pairing between and . In other words, is a Fredholm operator with zero index. It can be easily proved that the homogeneous equation (14) has only a trivial solution. As noted above, we arrive at the following assertion.
Theorem 3. Equation (22) has a unique solution for any .
3. Spectral Boundary Galerkin Method
We approximate the unknown function in (22) by truncated series of Chebyshev polynomials with weighting function . As mentioned above, the solution may have weak singularity. Let and denote the Chebyshev polynomials of first and second kinds with real values over . Then Chebyshev polynomials are orthogonal with respect to the weighting functions
Next we construct the approximation spaces: Let us show approximation property over the associated Sobolev space.
Lemma 4.
(1) The space is a closed subspace of .
(2) is dense in .
Proof. Obviously, belongs to . For , , s.t.
According to the Weierstrass approximation theorem, there is an integer such that
where is the Bernstein polynomial of order . On the basis of the definition of and , we have
Take the transform
Obviously, . So there is a polynomial such that
Thus
The proof is complete.
Obviously, the next lemma follows.
Lemma 5.
(1) The space is a closed subspace of .
(2) is dense in .
The idea of constructing the finite element space is to describe through the truncated expansion:
The spectral Galerkin variational formulation of problem (22) is to find , and , Thanking [17], the following estimation holds:
By using of the relation (22) and (14) and the finite dimensional subspace introduced above, we can rewrite the spectral Galerkin method as to find , such that It can be shown that Rewrite (37) in the form It is convenient to evaluated integrals in (40) by using the GaussChebyshev integration rule of the first kind with , , and the second kind where , .
For the first integral in (40), it is Cauchy intergal and the following formula is useful:
If , , the integral can be evaluated by (42): In addition, there is an important condition (8), which can also be computed by (41):
If , , (8) can be satisfied naturally according to the definition of . We have which is obtained by directly using of the orthogonal relation of Chebyshev polynomials .
Assume that the stiffness matrix constructed by the first term in (40) is denoted by . In the second case, is a diagonal matrix
Consider the computation of the second term of (40). The contribution for stiffness matrix from the second term in (40) is denoted by . According to (41) and (42), the entries of matrix
Therefore, the stiffness matrix for (40) is . The third term can also be evaluated by (42):
Finally, we get a linear algebraic system of equations with respect to .
Next, we recover singlelayer potential and angular potential by (41):
It is easy to see that we can get approximate electricfield potential by
Although (6) can also be expressed by a doublelayer potential, such an expression has more strong singularity near compared with (6). So (51) is deduced by (6).
4. Numerical Examples
In Section 3 we develop a method for finding the numerical solution of the integral equation (22). More precisely, a spectral Galerkin method is used. We reduce it to a linear algebraic system. In this section, we describe the results of the numerical tests mentioned above. We examine the feasibility of the numerical method for solving (22) and compare the numerical solution with these solutions of [11] for certain function on the righthand side of the integral equation under the assumption that is a straight line or an arc of the unit circle.
Example 6. Let be a segment of the straight line: and the parameters , . The corresponding exact solution of integral equation (22) can be written as The solution of problem (2) can be expressed by (11). The exact analytical expressions for potentials under the assumption are where , , , .
Now, let us consider the numerical solution at this time. Obviously, The spectral Galerkin approximate solution of is equal to the exact solution.
This situation allows us to test the accuracy of the numerical method of potentials. We choose as the points where we compute the potentials. errors are listed in Table 1, where the numerical solutions are compared with analytical solutions . Table 1 also contains the errors of the potentials which are computed by the methods in [11].

It can be seen that the computation errors of the potentials is on the good control. The errors increase as decreases, because logarithmic function in (50) and in (51) have singularities when is very small.
As shown in Table 1, the numbers of singlepotential are superior to those in [11] when , but inferior to those when . That is because (50) is not as accurate as the approximation formula in [11] when is very small. But it is adequate enough for our purpose and our formula is easier to implement.
According to Table 1, our results of angular potential have higher accuracy than results of [11]. In [11], the author subdivides each interval of length into some subintervals uniformly in order to improve the accuracy. The values of on subintervals are computed by interpolation. In this paper, we need to increase so that to decrease the errors. The potential can be evaluated easily without interpolation since the numerical result of is a function. The interpolation procedure may deduced the related accuracy.
The straightline case is too special to show the accuracy of the method for arbitrary arcs. Therefor we consider further numerical tests on a curved arc.
Example 7. Suppose is an arc of the unit circle If , then the density function can be given by If , the exact density function is The potentials under these two assumptions cannot be given by explicit expressions. The comparison between numerical solution and the analytical solution is not convenient because of the space they belong to. Let , . errors between and are listed in Tables 2 and 4. We observe excellent agreement between them. The accuracy of increases if increases. But no significant change can be seen whatever increases or not. The numerical solution of [11] is compared with the analytical solution in Tables 3 and 5. Results in Tables 2 and 4 are smaller than those in Tables 3 and 5. The computational schemes in this paper are stable and efficient.




Let , be the angular potential and singlepotential computed by (50) and (51). We use the exact density function to get them. Let , be the angular potential and singlepotential used by numerical approximation of the density function.
We compute the potentials in a parallel surface of defined by where . Results of the tests are presented in Tables 6 and 7. Obviously, errors decrease as increases. When , it can be seen that the approximation is satisfactory even when . The simulation of is not as good as the simulation of . The errors when with are almost the same as the errors when under the condition of . But the accuracy for both cases is high, when .


Example 8. In this example, we take the arc and the parameters We compute the singlepotential and angular potential on a semicircular arc: see Figures 1 and 2. We choose and to compute the electricfield potential which is shown in Figure 3.
5. Conclusions
In this paper, we have demonstrated an efficient and highly accurate boundary Galerkin method for a skewderivative problem. The problem is considered in Sobolev spaces firstly. In our method, approximate calculation of singular integral can be avoided. The density function can be evaluated at any point without interpolation once the series coefficients have been determined. It is direct, simple, and fast.
Acknowledgment
The work of the first author is supported by the Youth Natural Science Foundation of Heilongjiang (Grant no. QC2011c103).
References
 S. A. Yang, “A numerical method for scattering from acoustically soft and hard thin bodies in two dimensions,” Journal of Sound and Vibration, vol. 250, no. 5, pp. 773–793, 2002. View at: Publisher Site  Google Scholar
 R. Kress and K.M. Lee, “Integral equation methods for scattering from an impedance crack,” Journal of Computational and Applied Mathematics, vol. 161, no. 1, pp. 161–177, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. JerezHanckes and J. C. Nedelec, “Hybrid Galerkin method for elliptic and wave propagation problems in ${\mathbb{R}}^{3}$ over planar structures,” in Integral Methods in Science and Engineering, vol. 2, chapter 19, Computational Methods, 2009. View at: Google Scholar
 F. Cakoni and D. Colton, “The linear sampling method for cracks,” Inverse Problems, vol. 19, no. 2, pp. 279–295, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. A. Krutitskii, N. C. Krutitskaya, and G. Y. Malysheva, “A problem related to the hall effect in a semiconductor with an electrode of an arbitrary shape,” Mathematical Problems in Engineering, vol. 5, no. 1, pp. 83–95, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. A. Krutitskiĭ, “The Dirichlet problem for the Helmholtz equation in the exterior of cuts in the plane,” Computational Mathematics and Mathematical Physics, vol. 34, no. 89, pp. 1237–1258, 1994. View at: Google Scholar  MathSciNet
 P. A. Krutitskiĭ and V. V. Kolybasova, “A generalization of the Neumann problem for the Helmholtz equation outside cuts on the plane,” Differential Equations, vol. 41, no. 9, pp. 1213–1224, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. A. Krutitskiĭ, “On the properties of solutions of the Dirichlet problem for the Helmholtz equation in a twodimensional domain with cuts,” Differential Equations, vol. 43, no. 9, pp. 1200–1212, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 P. A. Krutitskiĭ, “A boundary value problem for the Laplace equation outside cuts on the plane with different conditions of the third kind on opposite sides of the cuts,” Differential Equations, vol. 45, no. 1, pp. 86–100, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 J. Liu, P. A. Krutitskii, and M. Sini, “Numerical solution of the scattering problem for acoustic waves by a twosided crack in 2dimensional space,” Journal of Computational Mathematics, vol. 29, no. 2, pp. 141–166, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. A. Krutitskii, D. Y. Kwak, and Y. K. Hyon, “Numerical treatment of a skewderivative problem for the Laplace equation in the exterior of an open arc,” Journal of Engineering Mathematics, vol. 59, no. 1, pp. 25–60, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Ammari, J. Garnier, H. Kang, W.K. Park, and K. Sølna, “Imaging schemes for perfectly conducting cracks,” SIAM Journal on Applied Mathematics, vol. 71, no. 1, pp. 68–91, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 H. Ammari, H. Kang, and H. Lee, Layer Potential Techniques in Spectral Analysis, vol. 153, American Mathematical Society, Providence, RI, USA, 2009. View at: MathSciNet
 H. Ammari, H. Kang, H. Lee, and W.K. Park, “Asymptotic imaging of perfectly conducting cracks,” SIAM Journal on Scientific Computing, vol. 32, no. 2, pp. 894–922, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge University Press, Cambridge, UK, 2000. View at: MathSciNet
 G. C. Hsiao and W. L. Wendland, Boundary Integral Equations, vol. 164 of Applied Mathematical Sciences, Springer, Berlin, Germany, 2008. View at: Publisher Site  MathSciNet
 R. Li, The Galerkin Method of Boundary Value Problem, Shanghai Scitific and Technical Publishers Press, 1988.
Copyright
Copyright © 2013 Wei Sun and Fuming Ma. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.