International Journal of Differential EquationsVolume 2011 (2011), Article ID 169707, 9 pagesdoi:10.1155/2011/169707

Research Article

A Comparison of Two Methods for Solving Electromagnetic Field Integral Equation

L. Jebli

Groupe Canal, Radio et Propagation, LPHE-Modélisation et Simulation, Avenue Ibn Battouta, Faculté des Sciences, Rabat, Morocco

Received 19 December 2010; Revised 11 March 2011; Accepted 2 May 2011

Academic Editor: Sabri Arik

Copyright © 2011 L. Jebli. 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.

Abstract.

The present paper aims to compare Harrington's direct method of moment (MoM) with the conjugate gradient method (CGM) by evaluating the total current solving the electric field integral equation (EFIE). Based on their performances, the number of iterations needed for convergence, storage, and the level of precision, it is found that the direct MoM is more efficient than other iterative CGM.

1. Introduction

In this paper, we consider the integral electromagnetic field equation (EFIE) for computing the total current distribution on antenna surfaces, and we study the two standard methods by comparing the solutions. These methods are the direct method of moment (MoM for short) developed by Harrington [1] and the iterative conjugate gradient method (CGM) [2]. This comparative study, which is done at the level of numerical solutions, is achieved in terms of computational advantages and disadvantages as well as their ability to determine the current distributions on metal surfaces. Based on their performance, the number of iterations needed for convergence, the storage, and the level of precision, it comes out from our analysis, developed in the next sections, that the direct method MoM is more efficient than other iterative CGM.

Recall that the electromagnetic field has been successfully used in the computation of the current distribution over the surface with few geometries yielding analytic solutions. To deal with the study of the above-mentioned integral electromagnetic field equation, the most popular method for getting numerical solutions is given by the method of moment with advantages and disadvantages. Regarding the disadvantages of this approach, we mention the one related to the convergence problem and the second connected with the restriction of the explicit storage of a the dense matrix. Furthermore, the monotonic convergence in the MoM solution can be achieved under only stringent condition which are difficult to realize in case of complex problems [2]. In attempt to address the above disadvantages of the direct MoM, another alternative approach has been suggested; this concerns the CGM for solving iteratively the electromagnetic problem [3, 4]. Since then, CGM has been successfully applied to a number of radiation and scattering problems and has been quite popular because of the O(N) memory requirement. Using discretization of EFIE in the framework of the MoM to get the impedance matrix [5], Sarkar and Rao applied to the matrix the CGM of steepest descent to calculate the charge distribution over conducting surfaces [6]. However, few attempts have been done to apply CGM to electrostatic problems involving conducting surfaces. The presentation of this paper is as follows: in Section 2, we give some generalities on EFIE. In Section 3, we review some useful aspects on the method of moment. In Section 4, we describe the algorithm of the CGM. In Section 5, we work out the numerical solutions of EFIE using both MoM and CGM. Last section is devoted to the conclusion.

2. Generalities on EFIE

Electromagnetic field has been successfully used in the computation of finding the current distribution over the surface. Solving current distributions on metal surface usually incorporates formulation of EFIE (electric field integral equation) [7]. Few geometries yields solutions (analytical, numerical approximation) have been successfully used. But most of the popular method is the method of moment.

By the method the initial operator equation can be formulated by employing boundary condition for tangential incident () and impressed () electric field; Scattering or radiation problems are essentially identical, the only difference is that the “incident” field for the driven antenna is the applied electric field in the feed. The total electric filed is a combination of the incident field and the scattered field, The incident electric field is either the incoming signal or the excitation electric field ion the antenna feed. The scattered electric field is due to surface currents and free charges on the metal surface , The magnetic vector potential describes surface current radiation whereas the electric potential describes radiation of surface. Free charges are respectively, and the three-dimensional Green's function is where . The electric field integrodifferential equation for the induced current is obtained by requiring that the total tangential electric field, , vanishes on the conductor surface , thus, giving the electric field integral equations by (2.1), The corresponding weak form of the equation is obtained by testing it with a vector-valued weighting function defined on and tangent to . Using the identity and the divergence theorem, the term involving the scalar potential is integrated by parts to obtain the weak form Concerning the surface current density, is expanded into the basis functions in the form , and using   as a test function, the moment equations are obtained below: where , The square impedance matrix determines electromagnetic interaction between different segments elements.

3. Useful Tools on MoM

The method of moments (MoM) is a well-known technique for solving linear equations. In antenna analysis, the MoM is used to convert the electric field integral equation into a matrix equation or system of linear equations (Harrington, 1968) [1]. Let us consider the inhomogeneous equation: where is a linear operator, is known, and is to be determined (is the unknown function). We will now perform the two essential steps we have highlighted above. In order to create the matrix equation, the unknown function is defined to be the sum of a set of known independent functions, , called basis or expansion functions with unknown amplitudes, .(i)Firstly, let the basis functions be . The unknown function is expanded in terms of a linear combination of these basis functions: The linear combination of should represent the unknown in the domain. Substitute (3.1) into (2.9), we have that, upon using the linearity of the operator , the unknown amplitudes can be brought out of the operator giving where ,  are unknown coefficients that have to be determined.(ii)Secondly, we define a set of weighting functions (testing functions), , which are integrated with (3.3) to give different linear equations, and take the inner product of the previous equation with : The system can now be written in matrix form as where If the matrix is not singular, the unknowns are simply given by and the original function can be reconstructed using (3.3). We can now generalize the following definitions.

Basis FunctionsBasis functions can use full domain functions such as special functions, polynomials, and so forth. A set that is useful for practical problem is the subsectional basis function. This means that each is only nonzero over a subsection of the domain of . A common choice is the pulse function:

Weighting Functions (Testing Functions) Two common choices are as follow.(i)Point matching: taking Dirac delt functions as testing functions ( where .(ii)Galerkin Method: if we used the same function for both the basis and weighting , that is. If we used the same function for both where .

4. Conjugate Gradient Method

The conjugate gradient method is the most prominent iterative method for solving sparse systems of linear equations. It is applied to the analysis of radiation from thin wire antennas. With this iterative technique, it is possible to solve electrically large arbitrarily oriented wire structures without storing any matrices as is conventionally done in the method of moments. The basic difference between the proposed method and Galerkin's method, for the same expansion functions, is that for the iterative technique we are solving a least squares problem. The conjugate gradient method (CGM) is an algorithm for the iterative solution of an operator equation of the form set in a Hilbert space. If is positive definite, the classical CGM converges to the exact solution in at most steps in an -dimensional space, and it proceeds as follows [2, 8]: It has been shown [8] that sequence of solution generated by this algorithm minimizes the error function where is the exact solution.

5. Numerical Results

The results are given through two examples concerning straight wire scatterers.

5.1. First Example

As a first example of a straight-wire scatterer, we consider the simple straight-wire scatterer of length and radius . The scatterer is excited by a plane wave with a 45° angle of incidence; this example has been taken from a well-known book [7]. Figure 1 show the computed currents (direct solution MoM) according to Pocklington equations compared with conjugate gradient method. In these figures, we present the result for the current distribution along the antenna that we obtain through our own developments using MATLAB. This comparison shows that the movement is retained with some minor differences in middle of the curve. Notice that our result (under MATLAB) agrees with the one published in [9]. Notice also that one of the essential tasks in our way of doing is the use of the moment method and the choice of test and basis functions.

(a)

(b)

Figure 1: Current distributions: comparison of results for the MoM and CGM currents. On the left current distribution for a length dipole antenna and on the right a current distribution for a length dipole antenna.

5.2. Second Example

In this example we consider the scattering from the straight-wire scatterer of length and radius shown in Figure 2; the current has a delta-gap compared with the sinusoidal current and with the current found by CGM.

Figure 2: Current distribution for a length dipole antenna, comparison of results for the MoM, CGM current, and sinusoidal current.

Figure 3: Solution for total current : correspondence between direct MoM solution and iterative CGM.

6. Conclusion

In this paper, we have presented a comparative study of the direct method of moment and the iterative conjugate gradient method. These methods are used to solve the differential linear equations; in particular, the integral electromagnetic field equation (EFIE) for computing the total current distribution. Based on their performances, the number of iterations needed for convergence, storage, and the level of precision, we have found that the direct MoM is more efficient than other iterative CGM.

Acknowledgments

The author would like to thank professor E. H Saidi for helpful discussions and URAC 09 for support.

References

R. F. Harrington, Field Computation by Moment Methods, Macmillan, New York, NY, USA, 1968.

V. V. Sanjaynath, N. Balakrishnan, and G. R. Nagabhushana, “Application of conjugate gradient method for static problems involving conductors of arbitrary shape,” IEEE Transactions on Antennas and Propagation, vol. 42, no. 7, pp. 1028–1033, 1994.

T. K. Sarkar and S. M. Rao, “The application of the conjugate gradient method for the solution of electromagnetic scattering from arbitrarily oriented wire antennas,” IEEE Transactions on Antennas and Propagation, vol. 32, no. 4, pp. 398–403, 1984.

T. K. Sarkar and E. Arvas, “On a class of finite step iterative methods (conjugate directions) for the solution of an operator equation arising in electromagnetics,” IEEE Transactions on Antennas and Propagation, vol. 33, pp. 1058–1066, 1985.

S. M. Rao, A. W. Glisson, D. R. Wilton, and B. S. Vidula, “A simple numerical solution procedure for statics problems involving arbitrary shaped conductors,” IEEE Transactions on Antennas and Propagation, vol. 27, pp. 604–607, 1979.

T. K. Sarkar and S. M. Rao, “An iterative method for solving electrostatic problems,” IEEE Transactions on Antennas and Propagation, vol. 30, no. 4, pp. 611–616, 1982.

R. F. Harrington, Field Computation by Moment Methods, Wiley-IEEE Press, New York, NY, USA, 1993.

R. M. Hayes, “Iterative methods of solving linear problems on Hilbert space,” in Contributions to the solution of systems of linear equations and the determination of eigenvalues, National Bureau of Standards Applied Mathematics Series 39, pp. 71–103, U. S. Deptartment of Commerce, Washington, DC, USA, 1954.

V. Barrera-Figueroa, J. Sosa-Pedroza, and J. López-Bonilla, “Simplification of Pocklington's integral equation for arbitrary bent thin wires,” in Boundary Elements XXVII, vol. 39 of WIT Transactions on Modelling and Simulation, pp. 563–574, WIT Press, Southampton, UK, 2005.