International Journal of Microwave Science and Technology
Volume 2015 (2015), Article ID 724702, 7 pages
http://dx.doi.org/10.1155/2015/724702

Research Article

An Effective Math Model for Eliminating Interior Resonance Problems of EM Scattering

Zhang Yun-feng,1,2 Zhou Zhong-shan,1,2 Su Zhi-guo,² Wang Rong-zhu,² and Chen Ze-huang²

¹Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China
²Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science, No. 219, Ningliu Road, Nanjing 210044, China

Received 22 May 2014; Revised 26 November 2014; Accepted 15 December 2014

Academic Editor: Giancarlo Bartolucci

Copyright © 2015 Zhang Yun-feng et al. 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.

It is well-known that if an -field integral equation or an -field integral equation is applied alone in analysis of EM scattering from a conducting body, the solution to the equation will be either nonunique or unstable at the vicinity of a certain interior frequency. An effective math model is presented here, providing an easy way to deal with this situation. At the interior resonant frequencies, the surface current density is divided into two parts: an induced surface current caused by the incident field and a resonance surface current associated with the interior resonance mode. In this paper, the presented model, based on electric field integral equation and orthogonal modal theory, is used here to filter out resonant mode; therefore, unique and stable solution will be obtained. The proposed method possesses the merits of clarity in concept and simplicity in computation. A good agreement is achieved between the calculated results and those obtained by other methods in both 2D and 3D EM scattering.

1. Introduction

Electric field integral equation (EFIE) and magnetic field integral equation (MFIE) have been widely employed to analyze electromagnetic scattering of conducting bodies [1–3]. However, interior resonance phenomena exist in solving electromagnetic scattering problems with surface integral equations. When working frequencies of conductors are near (or exactly at) to the frequencies associated with interior resonances, the single equation will become highly ill-conditioned (or singular) which makes the solution unstable or nonunique. Also, the interior resonance behavior has significant influence on the late time stability associated with time domain EFIE and MFIE [4, 5].

Several ways of dealing with this numerical problem have been proposed. Nowadays, the popular combined field integral equation (CFIE) technique to overcome this problem is a proper combination of the electric field integral equation and the magnetic field equation [6, 7]. The CFIE technique requires the calculation of both and impedance matrices and it is not suitable for aperture problems. The combined source integral equation (CSIE) [8, 9] technique makes up for aperture structures. Also a technique has been proposed by Mittra and Klein [10], involving application of the generalized boundary condition [11], and consists of additional points in the interior of the conductor and forces the field to be zero at those points. The problem with this technique is the fact that the chosen interior points must be carefully selected so as not to lie on nodal lines, which is not too practical for large bodies of simple shape, for which a slight change in frequency can take us from one resonance to the next, of different modal distribution. There are also some iterative methods reported in the works of Sarkar and Ergül [12, 13] to deal with this situationand they are used to compute the minimum norm solution (which produces the correct scattered fields but not the true tangential fields) and to calculate the LQSR. Another work related to the use of extended integral equations has been presented by Mautz and Harrington [8], which involves application of the boundary element method, with observation points lying on an internal closed surface near the boundary of the scatterer. Unfortunately, the resulting matrix equations are also ill-conditioned since the internal contour can resonate by itself. The authors’ proposed scheme of allowing the internal contour to vary with the wave number seems impractical and in reality does not solve the problem except in some isolated cases of electrically small simple shapes whose resonances are known. More recently, Canning [14] illustrated a matrix algebra technique known as the singular value decomposition (SVD) has been proposed for moment method calculations involving perfect conductors. Such a technique diagonalizes the matrix equation, isolating the resonant contribution, which is then omitted in the calculation. We also refer to the interesting work of Yaghjian, who originally presented his augmented electric or magnetic field equation [15], and more recently Tobin et al. [16] pointed out their drawbacks for an arbitrarily shaped, multiwavelength body and, most importantly, introduced a modification, the so-called dual-surface integral equation, which is applicable to the perfectly conducting bodies and supposedly eliminates all the spurious solutions. Finally, the modal orthogonal characteristics [17, 18] were applied to solve the same problem in two-dimensional EM scattering.

Here, we present an effective method to solve the scattering problem of conductor bodies at or in the neighborhood of the resonant frequencies in both 2D and 3D EM scattering. At the resonant frequencies, Inagaki modes, firstly applied to analyze the antenna array [19], are employed here to be validated away from resonance with a unique and stable solution to the equation. It can both stabilize the numerical calculation and yield to reliable results of the surface current density and exterior field for conductors at interior resonances.

2. Theory

The principle is elaborated as follows. The solution of anill-conditioned system of the equation will consist of (1) the correct physical solution to the problem and (2) the resonant solution, which is in conjunction with Green’s Theorem which produces the nonzero complementary resonant modes. When the orthogonal modes are used to solve the electric field integral equation, the solutions will be divided into the induced modes and the resonant modes corresponding to eigenvalues. We can easily obtain the current density and exterior field from the induced modes.

2.1. Why Is the Solution Unstable?

Consider a PEC scatterer excited by an incident wave and the scatterer defined by the surface , in an impressed electric field (see Figure 1).

Figure 1: PEC scatterer impressed by electromagnetic field.

According to the boundary condition, we get

An operator equation for the current on is

The operator is defined by [2]

where denotes a field point, denotes a source point, and , , and denote angular frequency, permeability, and wave number, respectively, of free space; is the free space Green’ function,

Equation (2) is discretized into a set of simultaneous linear algebraic equations. This set of equations can be compactly rewritten in a matrix form as

where is the general impedance matrix, denotes the general voltage vector, and stands for the coefficient of the current density to be determined. Usually (noninterior resonant condition) the current density of the conductor can be obtained from (5) with a unique and stable solution, but, at the frequencies associated with interior resonances, the situation will be on the contrary.

As a result of interior resonance, the current density on at discrete resonant frequencies will be composed of two parts: an induced surface caused by the incident field and a resonance surface current associated with the interior resonance mode:

where is the resonant current density which is determined by the conducting body alone and it is independent of the incident wave; is the induced current density which is jointly determined by the incident field and the conducting body itself.

If the incident is removed, the electric field produced by the resonant current on satisfies the following homogeneous -field integral equation:

In other words, the solution to the inhomogeneous EFIE is not unique since the non-zero- solution to its corresponding homogeneous equation exits at some discrete frequencies. Namely, in addition to the induced current determined jointly by the incident field and the conducting body, there exists also the resonant current at some discrete frequencies determined alone by the conducting body itself. It follows from the nonuniqueness of the solution to EFIE that the solution to its corresponding moment matrix equation (5) will also not be unique at the interior resonant condition. The nonuniqueness of the moment matrix equation infers that the moment matrix is singular.

Similar to the behavior of a conducting cavity, the resonant current density on a conducting body is also not able to produce scattered field in space external to . Therefore, the scattered field and the radar cross section (RCS) external to determined by the electric field integral equation should be unique theoretically. Although resonant current theoretically does not change the exterior scattered field, it does make the method of moments analysis unstable and unreliable since the homogeneous operator equation has non-zero-solution .

2.2. How to Stabilize the Solution?

For the sake of simplicity, only the electric field integral equation is involved here. To solve the EFIE by method of moments, first of all, we should choose a set of expansion functions and a set of weighting functions. Here, we choose basis functions aiming at the diagonalization of the moment matrix which was proposed by Inagaki and Garbacz [19], in which the expansion functions are chosen to be the eigenfunctions of a composite Hermitian operator , where is the adjoint operator of ,

the weighting functions to be the response of them

The orthogonal property of Inagaki modes [13] leads to the satisfaction of the moment matrix diagonalization condition; that is,

namely,

So the modal solution to the operator equation becomes

Here is the expansion function vector,

where is a column excitation vector given by

Then method of moments is employed here for computation of the modes. To do so, functions are taken as both expansion functions and weighting functions for a moment method analysis to the operator eigenvalue equation, giving

whose matrix form is

In this matrix eigenvalue equation, both and are matrices with their elements, respectively, as

And is a column vector composed of the unknown coefficients as

The matrix is simply the moment method matrix only if a least-squares method of moments is applied to the general operator equation in which are taken to be the expansion functions and their responses the weighting functions.

At the interior resonances, the surface current is composed of the resonant current and the induced current. The resonant current is the solution to the corresponding homogeneous -field integral equation. It follows from the nonuniqueness of the solution to -field integral equation that the solution to its corresponding moment matrix equation will also not be unique at the interior resonance condition. The nonuniqueness of the moment matrix equation infers that the moment matrix is singular.

Actually, the moment matrix will not be exactly singular, but rather, it will be highly ill-conditioned since the round-off and the truncated error exist during the computation of the matrix elements. According to the orthogonal property of Inagaki modes [13], the surface current of conductors is composed of serials of modal currents corresponding to the eigenvalue ,

At the interior resonances, when eigenvalue is much smaller than the other eigenvalue, the corresponding modal current is resonant current and the other is the induced current. These modal currents are mutually orthogonal, so to eliminate the interior resonant mode in the Inagaki mode solution, what we need to do is just to replace the number in []⁻¹ by zero. The incident current can be acquired from the following equation:

Theoretically, the scattered field and the radar cross section external to determined by the electric field integral equation should be unique. However, due to the ill-condition of the moment matrix, the exterior fields from conductors will be unstable and unreliable. After filtering out the resonant current, the exterior field, obtained from the induced current, is stable and reliable.

In fact, the resonance current does not really exist on the surface of the conductors. The interior resonance problem is just caused by the deficiency of the selected mathematical model. After getting rid of the virtual resonance current, we can obtain the stable and reliable current density and exterior field of conductors at interior resonances.

3. Numerical Results

Here we confine our attention to not only two-dimensional structures but also three-dimensional structures. The presented method is applied here to get the surface current density and exterior field of conductors when they are at (or near) the interior resonant frequencies.

3.1. An Infinitely Long Circular Cylinder

As the first case for testing the approach described above, scattering from a circular conducting cylinder was examined here. It was found that when the TM plane wave normally illuminates an infinite circular cylinder, the interior resonance takesplace if , where is the wave number and is the radius, which is the numerical resonant frequency point of E ₁₁ mode of the same surface circular cylinder cavity.

From Figures 2 and 3, we can see that the magnitude of induced current is much smaller than that of resonant current, so the interior resonant current has completely masked the true current responsible for the scattered field. The true surface current, computed by the presented method, agrees well with the exact analytical solution. The bistatic RCS of the circular cylinder computed from the induced current is compared with the exact solution. It is clear that there are excellent agreements between the calculated results and the analytic results, as it is shown in Figure 4.

Figure 2: Resonant current on a circular cylinder, = 3.8214.

Figure 3: Induced current on a circular cylinder, = 3.8214.

Figure 4: Bistatic RCS of a circular cylinder, = 3.8214.

3.2. An Infinitely Long Square Cylinder

Also, an infinitely long square conducting cylinder was considered here with a TM wave incident along the axis of the cylinder. The first and the second resonant frequency of the square cylinder are, respectively, at and at theoretically. We can get the first numerical resonant point at , using the presented method, as shown in Figure 5. Here, only the first resonant frequency is considered. At the first resonant frequency, the resonant current and induced current are shown, respectively, in Figures 6 and 7. From Figure 7, we can see that the induced current is compared well with that obtained by CFIE technique.

Figure 5: Minimal eigenvalue of square cylinder.

Figure 6: Resonant current on a square cylinder, = 0.70775.

Figure 7: Induced current on a square cylinder, = 0.70775.

In Figure 8, the bistatic RCS of the square cylinder obtained through three different methods is depicted. Due to the ill-conditioned equation, the RCS of square cylinder computed by single EFIE (dash line) obviates the true result. After we filtered out the resonant modal current, the RCS (dot-dash line) produced by the induced current coincides with that obtained by CFIE method (solid line).

Figure 8: Bistatic RCS on square cylinder, = 0.70775 obtained by presented method (labeled IM), CFIE method, and EFIE.

3.3. Two Conducting Spheres

In the end, two conducting spheres, respectively, at the and , are analyzed by the presented method, where denotes the radius of the sphere. At the interior resonant frequency, the bistatic RCS, obtained by the single EFIE method (dash line), is away from the right value. After being corrected by the presented method (dot-dash line), the bistatic RCS coincides with that calculated by CFIE technique (solid line), as shown in Figures 9 and 10.

(a)

(b)

Figure 9: Bistatic RCS of a resonant sphere, = 2.768 (a) -polarization, (b) -polarization obtained by presented method (labeled IM), CFIE method, and EFIE.

(a)

(b)

Figure 10: Bistatic RCS of a resonant sphere, = 4.518 (a) -polarization, (b) -polarization obtained by presented method (labeled IM), CFIE method, and EFIE.

4. Conclusion

A new scheme for eliminating interior resonance problems associated with surface integral equation is presented. The orthogonal property of Inagaki modes is used here to isolate the resonant mode, which is then omitted in the computation to obtain the right property of the conductors at the interior resonances. This simple technique has been proven to be effective in attenuating the resonant modes and getting the unique and stable solution to EFIE when conductors are at (or near) the frequencies associated with interior resonances. Excellent numerical results, away from resonance problems, have been obtained for some shapes not only in two dimensions but also in three dimensions. Compared with other techniques, the advantages of this approach are only involved EFIE equation and easy to get the right results, but, due to the determination of eigenvalues and eigenvectors, there is a bit time-consuming for property computation at or in the interior resonance.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

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

R. F. Harrington, Time Harmonic Electromagnetic Field, McGraw-Hill, New York, NY, USA, 1968.

W. C. Chew and J. M. Song, “Gedanken experiments to understand the internal resonance problems of electromagnetic scattering,” Electromagnetics, vol. 27, no. 8, pp. 457–471, 2007.

H. Lianrong, G. Tian, J. Fang, and G. Xiao, “The behavior of MFIE and EFIE at interior resonances and its impact in MOT late time stability,” American Journal of Electromagnetics and Applications, vol. 1, no. 2, pp. 30–37, 2013.

F. P. Andriulli, K. Cools, F. Olyslager, and E. Michielssen, “Time domain Calderón Identities and their application to the integral equation analysis of scattering by PEC objects part II: stability,” IEEE Transactions on Antennas and Propagation, vol. 57, no. 8, pp. 2365–2375, 2009.

J. R. Mautz and R. F. Harrington, “H-field, E-field and combined-field solutions for conducting bodies of revolution,” AEU, vol. 32, pp. 157–164, 1978.

F. P. Andriulli and E. Michielssen, “A regularized combined field integral equation for scattering from 2-D perfect electrically conducting objects,” IEEE Transactions on Antennas and Propagation, vol. 55, no. 9, pp. 2522–2529, 2007.

J. R. Mautz and R. F. Harrington, “A combined-source solution for radiation and scattering from a perfectly conducting body,” IEEE Transactions on Antennas and Propagation, vol. 27, no. 4, pp. 445–454, 1979.

K. F. A. Hussein, “Effect of internal resonance on the radar cross section and shield effectiveness of open spherical enclosures,” Progress in Electromagnetics Research, vol. 70, pp. 225–246, 2007.

R. Mittra and C. A. Klein, “Stability and convergence of moment method solutions,” in Numerical and Asymptotic Techniques in Electromagnetics, vol. 3 of Topics in Applied Physics, pp. 129–163, Springer, Berlin, Germany, 1975.

P. C. Waterman, “Matrix formulation of electromagnetic scattering,” Proceedings of the IEEE, vol. 53, no. 8, pp. 805–812, 1965.

T. K. Sarkar and S. M. Rao, “A simple technique for solving -Field integral equations for conducting bodies at internal resonances,” IEEE Transactions on Antennas and Propagation, vol. 30, no. 6, pp. 1250–1254, 1982.

Ö. Ergül and L. Gürel, “Efficient solution of the electric-field integral equation using the iterative LSQR algorithm,” IEEE Antennas and Wireless Propagation Letters, vol. 7, pp. 36–39, 2008.

F. X. Canning, “Singular value decomposition of integral equations of EM and applications to the cavity resonance problem,” IEEE Transactions on Antennas and Propagation, vol. 37, no. 9, pp. 1156–1163, 1989.

A. D. Yaghjian, “Augmented electric and magnetic field equations,” Radio Science, vol. 16, no. 6, pp. 987–1001, 1981.

A. R. Tobin, A. D. Yaghjian, and M. M. Bell, “Surface integral equations for multiwavelength, arbitrarily shaped, perfectly conducting bodies,” in Proceedings of the URSI National Radio Science Meeting Digest, pp. 12–15, Boulder, Colo, USA, January 1987.

W. Cao and J. Chen, “The application of bi-orthogonal mode analysis approach to in electromagnetic scattering of conducting bodies at interior resonances,” Chinese Journal of Radio Science, vol. 10, pp. 16–22, 1995.

Y.-F. Zhang, C.-G. Jiang, and W. Cao, “Inagaki mode approach to electromagnetic scattering of conducting bodies at interior resonances,” Journal of Electronics and Information Technology, vol. 28, no. 9, pp. 1735–1739, 2006.

N. Inagaki and R. J. Garbacz, “Eigenfunctions of composite Hermitian operator with application to discrete and continuous radiating systems,” IEEE Transactions on Antennas and Propagation, vol. 30, no. 4, pp. 571–575, 1982.