“JOURNAL OF RADIO ELECTRONICS” N 10, 2013

contents

SCATTERING OF DIPOLE FIELD BY PERFECTLY CONDUCTING DISK

Vadim A. Kaloshin 1, Kirill K. Klionovski 2

1 Kotelnikov Institute of Radioengineering and Electronics of RAS

2 JSC «Concern «Morinformsystem-Agat»»

The paper is received on December 2, 2013

 

Abstract. Asymptotic formulas for calculating a scattering field of an arbitrarily oriented magnetic and electric dipole located on axis of a perfectly conducting disk were obtained using the method of physical theory of diffraction and the method of uniform asymptotic theory of diffraction. Radiation patterns were calculated using the asymptotic formulas and a numerical solution of a singular integral equation. The obtained results were compared to ones available in literature. A dependence of the front-to-back ratio on the disk’s radius was investigated.

Keywords: perfectly conducting disk, spherical wave, dipole, physical theory of diffraction, uniform asymptotic theory of diffraction.

 

Introduction.

Scattering of electromagnetic waves by a perfectly conducting disk was researched in [1-33]. In particular, scattering of a plane electromagnetic wave by a perfectly conducting disk was considered in [1-15]. Scattering of field of an electric and magnetic dipole located on axis of a disk was researched in [16-32]. A paper [33] describes a scattering field of an electric dipole located arbitrarily relative to a disk. A special case of a dipole mounted on a surface of a disk was presented in [16-25]. Asymptotic expressions of a scattering field of an electric dipole with perpendicular to disc orientation are available in [19, 25, 29]. Asymptotic formulas of a scattering field of a magnetic dipole with parallel orientation are available in [20]. A rigorous solution is available in a case of a perpendicular electric [16, 17] and parallel magnetic [17] dipole. This solution is represented as a series of spheroidal functions.

  

The statement of problem.

The first purpose of present paper is obtaining asymptotic expressions for a scattering pattern of an electric and magnetic arbitrarily oriented dipole located on an axis of a perfectly conducting disk. The second purpose of the paper is obtaining a numerical solution of this problem. The third one is to compare numerical results obtained using the asymptotic expressions and the numerical solution with numerical results from the papers [19, 20, 29].

We considered two cases - the dipole oriented parallel and normally to the disk (Fig. 1). It is enough to determine the field of the arbitrarily oriented dipole.

Fig. 1. The dipoles on axis of the perfectly conducting disk

 

The asymptotic expressions for scattering pattern derived separately near and far from the axis Z. We use the uniform asymptotic theory of diffraction [34, 35] to determine the scattering field far from the axis, and results of paper [36] near the axis. The results of paper [36] were obtained in the physical theory of diffraction approximation [3].

 

Fields of the dipoles in free space.

Consider a perpendicular magnetic dipole

                                                                          (1)

and a parallel magnetic dipole

                                                                         (2)

in free space where, x, y, z are Cartesian coordinates; z0 and x0 are the unit vectors in direction of the Z-axis and the X-axis, respectively; δ(x) – is the Dirac delta function; mz, mx – are moments of the magnetic dipoles.

Azimuthal and meridional components of field of the magnetic dipoles in free space are 

                                                              (3)

in a case of the perpendicular magnetic dipole, and

                                                       (4)

in a case of the parallel magnetic dipole. Here, λ is wavelength; r, θ, φ are spherical coordinates. Bottom indexes θ and φ correspond to the meridional and azimuthal field component, respectively.

 Consider a perpendicular electric dipole

                                                                          (5)

and a parallel electric dipole

                                                                          (6)

Here, pz, px – are moments of the electric dipoles. Components of field of the perpendicular and parallel electric dipole in free space are respectively

                                               (7)

and

                                        (8)

 

The asymptotic expressions of field far from the axis.

Let’s define the asymptotic expressions for meridional  and azimuthal  components of field of the electric and magnetic dipoles. The dipoles locate on the axis of the disk of radius R on distance h above the disk (Fig. 1). These expressions for the angles far from the Z-axis in the uniform asymptotic theory of diffraction approximation are

 

           (9)

 

         (10)

The meridional  and azimuthal  component in (9), (10) corresponds to the meridional and azimuthal component in (3), (4), (7) or (8) of analyzed dipole , ,  or , respectively. The variable  for a perpendicular magnetic dipole or a parallel electric dipole. The variable  for a parallel magnetic dipole or a perpendicular electric dipole. The variable  for  perpendicular magnetic and electric dipoles, and  for parallel magnetic and electric dipoles. In (9), (10)  is the Heaviside function. Function  is first term of asymptotic representation of the Fresnel integral . The arguments of the functions  and  are

        (11)

Formulas (9) and (10) take into account the first and the second order of diffraction. The last term in (9) describes the second order of diffraction and it is determined using the uniform geometrical theory of diffraction.

 

The asymptotic expressions of field near the axis.

Let’s define asymptotic expressions near from the Z-axis. These expressions in the physical theory of diffraction approximation have the form

                          (12)

for the meridional component of field, and

                         (13)

for the azimuthal component of field. The scattering field in (12), (13) is determined as

                              (14)

for the perpendicular magnetic dipole,

      (15)               

for the parallel magnetic dipole,

                     (16)

for the perpendicular electric dipole, and

  (17)

for the parallel electric dipole.

The meridional component  and the azimuthal component  in (12)-(17) correspond to the meridional and azimuthal component in (3), (4), (7) or (8) of analyzed dipole , ,  or , respectively. In (14)-(17) Jn(u) is the Bessel function of order n and argument u, and

(18)

 

Let’s define the front-to-back ratio  for the parallel electric and magnetic dipoles. We define asymptotic expressions of the front-to-back ratio in the physical theory of diffraction using the asymptotic expressions (12), (15) and (17) near the Z-axis:

                                    (19)

for the parallel magnetic dipole, and

              (20)

for the parallel electric dipole, where

The expressions (19), (20) take the next form when the radius of the disk tends to infinity

                                                        (21)

for the parallel magnetic dipole, and

                                                                                                                (22)

for the parallel electric dipole. Expression (21) shows that the front-to-back ratio of the magnetic dipole oscillate relative a constant when the radius of the disk increase. This constant and amplitude of oscillation doesn’t depend on the disk’s radius. This dependence of  can be explained by the fact that an amplitude of field of the magnetic dipole on an edge of the disk decreases proportionally R. At the same time, amplitude of scattering field by the edge increases proportionally R. As a result, the front-to-back ratio oscillates relative to a constant value.

Expression (22) shows that  of the electric dipole is infinite when . In this case, amplitude of field of the dipole on the edge of the disk decreases proportionally , and  increases in proportion to the disk’s radius.

 

Numerical research of scattering by a disk.

The scattering by a disc was investigated by solving a singular integral equation of the first kind

                                                              (23)

numerically. Here,  is the Green’s tensor;  is the unknown electric current on the disk;  is electric field of a dipole in free space;  is surface of the disk. We used the method of moments and an algorithm described in [37] to solve the integral equation (23). The unknown electric current  shall be decomposed over basis of triangular elements:

                                                        (24)

Here,  are unknown current’s amplitudes; r0 is the radial unit vector; φ0 is the azimuthal unit vector. Basic functions with bases 2Tr and 2Tφ have the form

                                    (25)

for perpendicular magnetic dipole, and

                                    (26)

parallel electric and magnetic dipoles were considered. The method of moments reduces the problem of scattering to a system of linear algebraic equations. A matrix of the unknown current’s amplitudes can be determined solving this system:

                                                                                                           (27)

Elements of matrix of own and mutual impedances  are calculated by numerical integration of the Green's tensor components  and the basic functions:

           (28)

A representation of components  in spectral form was used during calculations of the own and "nearest" mutual impedances in the matrix . The basic function was decomposed into a two-dimensional Fourier integral over flat sheets of electric current. Integral over the surface in (28) was determined by integrating fields of the flat sheets with components . A representation  in source-wise form for a ring of radial and azimuthal current in the spherical coordinates was used for purpose of reducing computation time in calculating the "distant" mutual impedances in the matrix . Elements of a matrix  which describes an interaction of field of a dipole and the current of the disk have the following form:

                                                                    (29)

The scattering field is determined by summing the field of dipole (3), (4), (7) or (8), and the field of the electric current of the disk. The field of the disk’s electrical current is determined by numerical integration of fields of radial and azimuthal electrical current’s rings. The integration was performed over radius of the rings with amplitude distribution (24).

The numerical results.

A pattern of the perpendicular magnetic dipole located on distance h=0,4λ above the disk of radius R=2λ and R=0.5λ is shown in Fig. 2 and 3, respectively. Red curve plots were obtained using numerical solution of the integral equation. Blue curve plots were obtained using asymptotic equations of the field near the Z-axis (13), (14). Green curve plots were obtained using asymptotic equation of the field far from the Z-axis (10).

Figures 4-7 show patterns of the parallel magnetic dipole which mounts on the surface of the disk of radius R=2λ and R=0.5λ. Curves in Fig. 4, 6 show patterns in Å-plane . Curves in Fig. 5, 7 show patterns in Í-plane . Red curve plots were obtained using numerical solution of the integral equation. Blue curve plots were obtained using asymptotic equations of the field near the Z-axis (12), (15) in Fig. 4, 6, and (13), (15) in Fig. 5, 7. Green curve plots were obtained using asymptotic equation of the field far from the Z-axis (9) in Fig. 4, 6, and (10) in Fig. 5, 7. Violet and cyan curves plot were obtained using asymptotic equations [20] of the field far and near from the Z-axis, respectively.

 

Fig. 2. Gain of the perpendicular magnetic dipole above the disk of radius R=2λ

 

 

Fig. 3. Gain of the perpendicular magnetic dipole above the disk of radius R=0.5λ

 

 

Fig. 4. Gain of the parallel magnetic dipole which mounts on the surface of the disk of radius R=2λ in Å-plane

 

 

Fig. 5. Gain of the parallel magnetic dipole which mounts on the surface of the disk of radius R=2λ in Í-plane

 

 

Fig. 6. Gain of the parallel magnetic dipole which mounts on the surface of the disk of radius R=0.5λ in Å-plane

  

 

Fig. 7. Gain of the parallel magnetic dipole which mounts on the surface of the disk of radius R=0.5λ in Í-plane

 

Patterns for a case of the perpendicular electric dipole which mounts on the surface of the disk of radius R=2λ and R=0.5λ are shown in Fig. 8 and 9, respectively. Blue curve plots were obtained using asymptotic equations of the field near the Z-axis (12), (16). Green curve plots were obtained using asymptotic equation of the field far from the Z-axis (9). Violet and cyan curves plot were obtained using asymptotic equations [19] of the field far and near from the Z-axis, respectively.

 

Fig. 8. Gain of the perpendicular electric dipole which mounts on the surface of the disk of radius R=2λ

 

 

Fig. 9. Gain of the perpendicular electric dipole which mounts on the surface of the disk of radius R=0.5λ

 

Figures 10-13 show patterns for the case of the parallel electric dipole located on distance h=0,4λ above the disk of radius R=2λ and R=0.5λ. Curves in Fig. 10, 12 show patterns in Å-plane . Curves in Fig. 11, 13 show patterns in Í-plane . Red curve plots were obtained using numerical solution of the integral equation. Blue curve plots were obtained using asymptotic equations of the field near the Z-axis (12), (17) in Fig. 10, 12, and (13), (17) in Fig. 11, 13. Green curve plots were obtained using asymptotic equation of the field far from the Z-axis (9) in Fig. 10, 12, and (10) in Fig. 11, 13.

 


Fig. 10. Gain of the parallel electric dipole above the disk of radius R=2λ in Å-plane

 

 


Fig. 11. Gain of the parallel electric dipole above the disk of radius R=2λ in Í-plane

 

 

Fig. 12. Gain of the parallel electric dipole above the disk of radius R=0.5λ in Å-plane

 

 

Fig. 13. Gain of the parallel electric dipole above the disk of radius R=0.5λ in Í-plane

 

Plots 14, 15 show dependence of the front-to-back ration of the parallel magnetic and electric dipoles on radius of the disk. The magnetic and electric dipoles locate on distances h=0,5λ and h=0,75λ above the disk, respectively. The disk’s radius is presented in proportions of wavelength. Red curves in Fig. 14 and 15 were obtained using asymptotic equations (19) and (20), respectively. Blue curve in Fig. 14 plots was obtained using asymptotic equation (21). Green points in Fig. 14 and 15 t were obtained using numerical solution of the integral equation.

 

Fig. 14. The front-to-back ratio of the parallel magnetic dipole above the disk

 

 

Fig. 15. The front-to-back ratio of the parallel electric dipole above the disk

 

Let’s discuss the calculation results in a case when radius of the disk is 2λ. The asymptotic expressions of field near the axis describe radiation pattern with relative accuracy less than 1 dB for sector of angles  and . The asymptotic expressions of field far from the axis describe radiation pattern with relative accuracy less than 1 dB for sector of angles , except the case of the parallel magnetic dipole in H-plane for angles . Curves of radiation pattern obtained using the asymptotic expressions near and far from the axis coincide for sectors of angles  and . Sector of coinciding increases when the disk radius increases. The relative accuracy of obtained asymptotic expressions decreases to 4 dB for angles  in a case of the disk’s radius is 0.5λ.

The relative accuracy of obtained asymptotic expressions of field near the axis is better than relative accuracy of the asymptotic expressions [20]. Analysis shows that asymptotic formulas obtained in [29] for radiation pattern of perpendicular electric dipole located on the disk are incorrect. 

The asymptotic expressions of the front-to-back ratio of the parallel electric and magnetic dipoles allow determining this ratio with accuracy less than 0.3 dB for any radius of the disk.

 

Conclusions.

Asymptotic expressions for calculating a scattering field of an arbitrarily oriented magnetic and electric dipole located on axis of a perfectly conducting disk were obtained. The asymptotic formulas make possible to calculate scattering field in full space. A relative accuracy of the asymptotic formulas is less than 1 dB when radius of the disk more then 2λ. The asymptotic expressions for the front-to-back ratio have the accuracy less than 0.3 dB for any radius of the disk.

 

REFERENCES

1.     Ëóðüå Ê. À. Äèôðàêöèÿ ïëîñêîé ýëåêòðîìàãíèòíîé âîëíû íà èäåàëüíî ïðîâîäÿùåì òîíêîì äèñêå // ÆÒÔ. 1959. ¹ 12. Ñ. 1421–1433.

2.     Hansen E. B. Scalar diffraction by an infinite strip and a circular disc // Journ. of Math. and Physics. V. 41. ¹ 3. 1962. pp. 229–245.

3.     Óôèìöåâ Ï. ß. Ìåòîä êðàåâûõ âîëí â ôèçè÷åñêîé òåîðèè äèôðàêöèè. Ì.: Ñîâ. ðàäèî. 1962.

4.     Õ¸íë Õ., Ìàóý À., Âåñòïôàëü Ê. Òåîðèÿ äèôðàêöèè. Ì.: Ìèð. 1964.

5.    Diane P. Marsland, Constantine A. Balanis, Stephen A. Brumley Higher order diffractions from a circular disk // IEEE Trans. on Anten. and Prop. V. AP-35. ¹ 12. Dec. 1987. pp. 1436–1444.

6.  Akashi T., Ando M., Kinoshita T. Effect of multiple diffraction in PTD analysis of scattered field from a conducting disk // Trans. IEICE. V. E-72. ¹ 4. 1989. pp. 259–261.

7.     Ando M., Kinoshita T. Accuracy comparison of PTD and PO for plane wave diffraction from large circular disk // Trans. IEICE. V. E-72. ¹ 11. 1989. pp. 1212–1218.

8.     Ando M., Kinoshita T. PO and PTD analysis in polarization prediction for plane wave diffraction from large circular disk // Digests of 1989 IEEE AP/S Int. Symp. V. 3. 1989. pp. 1282–1285.

9.     Allen K. Dominek Transient scattering analysis for a circular disk // IEEE Trans. on Anten. and Prop. V. AP-39. ¹ 6. June 1991. pp. 815–819.

10. Li L. W., Kooi P. S., Qiu Y. L., Yeo T. S., and Leong M. S. Analysis of electromagnetic scattering of conducting circular disk using a hybrid method // PIER. V. 20. 1998. pp. 101–123.

11. Zalipaev V., Kostin A. Scattering of a plane electromagnetic wave by a perfectly conducting disk in the high-frequency approximation // in proc. of int. seminar «Day on Diffraction 2001». 2001. pp. 311–319.

12. Hongo K. and Naqvi Q. A. Diffraction of electromagnetic wave by disk and circular hole in a perfectly conducting plane // PIER. V. 68. 2007. pp. 113–150.

13. Balaban M. V., Sauleau R., Benson T. M., and Nosich A. I. Dual integral expressions technique in electromagnetic wave scattering by a thin disk // PIER. V. 16. 2009. pp. 107–126.

14. Ñóõàðåâñêèé È. Î., Çàëåâñêèé Ã. Ñ., Íå÷èòàéëî Ñ. Â., Ñóõàðåâñêèé Î.È. Ðàññåÿíèå ýëåêòðîìàãíèòíîé âîëíû êðóãëûì èäåàëüíî ïðîâîäÿùèì äèñêîì êîíå÷íîé òîëùèíû // Ýëåêòðîìàãíèòíûå âîëíû è ýëåêòðîííûå ñèñòåìû. 2010. ¹ 2. Ñ. 42–47.

15. Hongo K., Jafri A. D. U., and Naqvi Q. A. Scattering of electromagnetic plane wave by a circular disk with surface impedance // PIER. V. 127. 2012. pp. 501–522.

16. Leitner A., Spence R. D. Effect of a Circular Ground Plane on Antenna Radiation // Journ. of Appl. Physics. Oct. 1950. V. 21. ¹10.

17. Áåëêèíà Ì. Ã. Äèôðàêöèÿ ýëåêòðîìàãíèòíûõ âîëí íà äèñêå // Â ñá.: Äèôðàêöèÿ ýëåêòðîìàãíèòíûõ âîëí íà íåêîòîðûõ òåëàõ âðàùåíèÿ. Ì.: Ñîâ. ðàäèî. 1957. Ñ. 148–174.

18. Lopez A. R. The geometrical theory of diffraction applied to antenna pattern and impedance calculation // IEEE Trans. on Anten. and Prop. V. AP-14. ¹ 1. Jan. 1966. pp. 40–45.

19. Ïèìåíîâ Þ. Â. Èçëó÷åíèå ýëåìåíòàðíîãî ýëåêòðè÷åñêîãî âèáðàòîðà, ðàñïîëîæåííîãî â öåíòðå èäåàëüíî ïðîâîäÿùåãî êðóãëîãî äèñêà // Àíòåííû. 1968. âûï. 3. Ñ. 93–102.

20. Ïèìåíîâ Þ. Â., Áðàóäå Ë. Ã. Èçëó÷åíèå ýëåìåíòàðíîãî ùåëåâîãî âèáðàòîðà, ðàñïîëîæåííîãî â öåíòðå èäåàëüíî ïðîâîäÿùåãî äèñêà // Àíòåííû. 1969. âûï. 6. Ñ. 89–105.

21. Hahn K.F., Fikioris J.G. Impedance and Radiation Pattern of Antennas above Flat Discs // IEEE Trans. on Anten. and Prop. Jan. 1973. V. AP-21. ¹ 1. pp. 97–100.

22. Awadalla K. H., Maclean T. S. M. Monopole antenna at center of circular ground plane: input impedance and radiation pattern // IEEE Trans. on Anten. and Prop. V. AP-27. ¹ 2. Mar. 1979. pp. 151–153.

23. Jack H. Richmond Monopole antenna on circular disk // IEEE Trans. on Anten. and Prop. V. AP-32. ¹ 12. Dec. 1984. pp. 1282–1287.

24. McNamara D. A., Pistorius C. W. I., Malherbe J. A. G. Introduction to the uniform geometrical theory of diffraction. Norwood, MA.: Artech house. 1990.

25. Volakis J. L. Antenna Engineering Handbook Fourth Edition. McGraw Hill Professional. 2007.

26. Ãðèíáåðã Ã. À., Ïèìåíîâ Þ. Â. Ê âîïðîñó î äèôðàêöèè ýëåêòðîìàãíèòíûõ âîëí íà áåñêîíå÷íî òîíêèõ èäåàëüíî ïðîâîäÿùèõ ïëîñêèõ ýêðàíàõ // ÆÒÔ. 1957. ¹10. Ñ. 2326–2339.

27. Inawashiro S. Diffraction of electromagnetic waves from an electric dipole by a conducting circular disk // Journal of Physics Society. Japan. V. 18. 1963. pp. 273–287.

28. Ïèìåíîâ Þ. Â., Äàâûäîâ À. Ã. Î âëèÿíèè èäåàëüíî ïðîâîäÿùåé ïëîñêîñòè ñ êðóãëûì îòâåðñòèåì íà èçëó÷åíèå ýëåìåíòàðíîãî ìàãíèòíîãî âèáðàòîðà // Ðàäèîòåõíèêà. 1978. ò. 33. ¹ 7. Ñ. 91–94.

29. Ïèìåíîâ Þ. Â., Äàâûäîâ À. Ã. Èçëó÷åíèå âåðòèêàëüíîãî ýëåêòðè÷åñêîãî âèáðàòîðà, ðàñïîëîæåííîãî íà îñè èäåàëüíî ïðîâîäÿùåãî äèñêà // Èçâ. ÂÓÇîâ. Ðàäèîôèçèêà. 1980. ò. 23. ¹ 1. Ñ. 101–112.

30. Duan D. W., Rahmat-Samii Y. and Mahon J. P. Scattering from a circular disk: Comparative study of PTD and GTD techniques // Proc. of IEEE. V. 79. ¹ 10. 1991. pp. 1472–1480.

31.  Murasaki T., Sato M., Ando M., Goto N. Diffraction Analysis of a Flat Plate for Dipole Wave Incidence Using the Modified Edge Representation // Proc. of ISAP. 1992. pp. 569–572.

32. Êðàâ÷åíêî Â. Ô., Åðîôååíêî Â. Ò. Äèôðàêöèÿ ýëåêòðîìàãíèòíûõ âîëí, èçëó÷àåìûõ ìàãíèòíûì äèïîëåì, íà ñâåðõïðîâîäÿùåì òîíêîì äèñêå // Ðàäèîòåõíèêà. 1995. ¹10. Ñ. 65–70.

33. Hongo K., Jafri A. D. U., and Naqvi Q. A. Scattering of electromagnetic spherical wave by a perfectly conducting disk // PIER. V. 129. 2012. pp. 315–343.

34. Áîðîâèêîâ Â. À., Ýéäóñ À. Ã. Î ñðàâíåíèè äâóõ ìåòîäîâ âû÷èñëåíèÿ ïîëÿ â îáëàñòè ïîëóòåíè // Àíòåííû. Âûï. 29. 1981. Ñ. 49–55.

35. Kouyoumjian R. G., Pathak P. H.  A uniform geometrical theory of diffraction for an edge in a perfectly conducting surface // Proc. IEEE. V. 62. ¹ 11. 1974. pp. 14481461.

36. Êàëîøèí Â. À., Ïîïîâ À. Ï. Ðàññåÿíèå íà îñåñèììåòðè÷íîé êðîìêå // ÐÝ. 1984. ¹ 8. Ñ. 1502–1509.

37. Êëèîíîâñêè Ê. Ê. Õàðàêòåðèñòèêè íàïðàâëåííîñòè ïðîâîäÿùèõ ýêðàíîâ êðóãëîé ôîðìû // Àíòåííû. 2011. ¹ 12. Ñ. 31–37.