THE UNIVERSITY OF MICHIGAN STUDIES IN RADAR CROSS-SECTIONS XXII ELEMENTARY SLOT RADIATORS R. F. Gpodrich, A. L. Maffett, N Reitlinger C E. Schensted, and K. M. Siegel C-. ontrect- AF.33(03)0-2 863L. —-".CHon'-PO0 LA-2F'5~o8,-r"l ~,^ 3.,.,......~.~~ 1. tti 5Noveber 6s,.,.''. -'.'2 7,24,. t 2l'72-13 -T THE UNIVERSITY OF MICHIGAN ENGINEERING RESEARCH INSTITUTE ANN ARBOR, MICHIGAN

1s d,',3, i o 111`1 K t~),1 1

THE UNIVERSITY OF MICH IGAN 2172-13-T STUDIES IN RADAR CROSS-SECTIONS I Scattering by a Prolate Spheroid, by F. V. Schultz (UMM-42, March 1950), W-33(038)-ac-14222, UNCLASSIFIED. II The Zeros of the Associated Legendre Functions Pmn (/) of NonIntegral Degree, by K. M. Siegel, D. M. Brown, H. E. Hunter, H. A. Alperin, and C. W. Quillen (Ubw-82, April 1951), W-33(038) -ac-14222, UNCLASSIFIED. III Scattering by a Cone, by K. M. Siegel and H. A. Alperin (11I-87, January 1952), AF-30(602)-9, UNCLASSIFIED. IV Comparison Between Theory and Experiment of the Cross-Section of a Cone, by K. M, Siegel, H. A. Alperin, J. W. Crispin, Jr., H. E. Hunter, R. E. Kleinman, W. C. Orthwein, and C. E. Schensted (Wbl-92, February 1953), AF-30(602)-9, UNCLASSIFIED. V An Examination of Bistatic Early Warning Radars, by K. M. Siegel (UIM-98, August 1952), W-33(038)-ac-14222, SECRET. VI Cross-Sections of Corner Reflectors and Other Multiple Scatterers at Microwave Frequencies, by R. R. Bonkowski, C. R. Lubitz, and C. E. Schensted (UMM-106, October 1953), AF-30(602)-9, SECRET - UNCLASSIFIED when Appendix is removed. VII Summary of Radar Cross-Section Studies Under Project Wizard, by K. M. Siegel, J. W. Crispin, Jr., and R. E. Kleinman (Ubaiv-108, November 1952), W-33(038)-ac-14222, SECRET. VIII Theoretical Cross-Section as a Function of Separation Angle Between Transmitter and Receiver at Small Wavelengths, by K. 1. Siegel, H. A. Alperin, R. R. Bonkowski, J. W. Crispin, Jr., A. L. Maffett, C. E. Schensted, and I. V. Schensted (IMM-115, October 1953), W-33(038)-ac-14222, UNCLASSIFIED. IX Electromagnetic Scattering by an Oblate Spheroid, by L. M. Rauch (UMM-116, October 1953), AF-30(602)-9, UNCLASSIF.IED. X Scattering of Electromagnetic Waves by Spheres, by H. Weil, M. L. Barasch, and T. A. Kaplan (2255-20-T, July 1956), AF-30(602)-1070, UNCLASSIFIED. ii

THE UNIVERSITY OF MICHIGAN 2472-13-T XI The Numerical Determination of the Radar Cross-Section of a Prolate Spheroid, by K. M. Siegel, B. H. Gere, I. Marx, and F. B. Sleator (UNIr-126, December 1953), AF-30(602)-9, UIJCLASSIFIED. XII Summary of Radar Cross-Section Studies Under Project MIRO, by K. M. Siegel, M. E. Anderson, R. R. Bonkowski, and W. C. Orthwein (UliIv-127, December 1953), AF-30(602)-9, SECRET. XIII Description of a Dynamic Measurement Program, by K. M. Siegel and J. M. Wolf (UIM-128, May 1954), W-33(038)-ac-14222, CONFIDENTIAL. XIV Radar Cross-Section of a Ballistic ilissile, by K. M. Siegel, i. L. Barasch, J. W. Crispin, Jr., W. C. Orthwein, I. V. Schensted, md H. Weil (UM32-134, September 1954), W-33(038)-ac-14222, SECRET. xV Radar Cross-Sections of B-47 and B-52 Aircraft, by C. E. Schensted, J. W. Crispin, Jr., and K. i. Siegel (2260-1-T, August 1954), AF-33(616)-2531, CONFIDENTIAL. XVI Microwave Reflection Characteristics of Buildings, by H. Weil, H. R. Bonkowski, T. A. Kaplan, and M. Leichter (2255-12-T, May 1955), AF-30(602)-1070, SECRET. XVII Complete Scattering Matrices and Circular Polarization CrossSections for the B-47 Aircraft at S-bana, by A. L. i, iaffett, V1. L. Barasch, w. E. Burdick, R. F. Goodrich, w. C. Orthwein, C. E. Schensted, and K. 1. Siegel (2260-6-T, June 1955), AF-33(61O)-2531, CON5FIDDENT IAL. XVIII Airborne Passive Measures and Countermeasures, by K. M. Siegel, 1. L. Barasch, J. W. Crispin, Jr., R. F. Goodrich, A. H. Halpin, A. L. i affett, W. C. Orthwein, C. E. Schensted, and C. J. Titus (2260-29-F, January 1956), AF-33(616)-253-, SECRET. XIX Radar Cross-Section of a Ballistic lMissile - II, by K. M, Siegel, M. L. Barasch, H. Brysk, J. w. Crispin, Jr., T. B. Curtz, and T. A. Kaplan (2428-3-T, January 1956), AF-04(645)-33, SECRET. XX Radar Cross-Section of Aircraft and Missiles, by K. M. Siegel, W. E. Burdick, J. W. Crispin, Jr., and S. Chapman (VJR-31-J, 1 March 1956), SEC ~ET. XXI Radar Cross-Section of a Ballistic Missile - III, by K. M. Siegel, H. Brysk, J. w. Crispin, Jr., and R. E. Kleinman (2428-19-T, October 1956) AF-04(645)-33, SECRET. XXII Elementary Slot Radiators, R.F.Goodrich, A.L.Maffett, N.Rei.tlinger, C.E.Schensted, and K.M.Siegel, (2472-13-T, November 1956), AF 33(038)-28631; HAC-PO L-265165-F31. UNCLASSIFIED......iii

THE UNIVERSITY OF MICHIGAN 2472-13-T TABLE OF CONTENTS Pa~e List of Studies in Radar Cross-Sections ii Preface v Chapter 1. Introduction and Statement of Problem. 1 Chapter 2. General Discussion of Radiation and Scattering Problems. 3 2.1 Reciprocity. 3 2.2 Geometric Optics. 4 2.3 Refinements of Geometric Optics. 7 Chapter 3. Cone Radiation Patterns by Optical Techniques. 15 Chapter 4. An Optimization Technique. 22 Chapter 5. Wiener Integral Methods. 32 Chapter 6. Conclusion and Prospectus. 43 Appendix A. Wedge 46 A-1 Electromagnetic Field for an Infinite Perfectly Conducting Wedge with Uniformly Excited Slots Parallel to its Edge. 46 A-l.l Exact Solution. 47 A-1.2 Asymptotic Solutions and Estimate of Error of Asymptotic Solutions. 50 A-1.3 Geometric Optics Solution. 68 A-l.4 Comparison Between the Exact and Asymptotic Solution and Calculation of the Error of the Asymptotic Expression. 70 Appendix B. Cone. 74 Appendix C. Graphical Presentations. 80 References. 99 ivs

THE UNIVE RS ITY OF MI CHIGAN 2472-13-T PREFACE This paper is the twenty-second in a series growing out of studies of radar cross-sections at the Engineering Research Institute of The University of Michigan. The primary aims of this program are: 1. To show that radar cross-sections can be determined analytically. 2, A. To determine means for computing the radiation patterns from antennas by approximate techniques which determine the pattern to the accuracy required in military problems but which do not require the unique determination of exact solutions. B. To determine means for computing the radar cross-sections of various objects of military interest, (Since 2A and 2B are inter-related by the reciprocity theorem it is necessary to solve only one of these problems) 3. To demonstrate that these theoretical cross-sections and theoretically determined radiation patterns are in agreement with experimentally determined ones. Intermediate objectives are: 1. A. To compute the exact theoretical cross-sections of various simple bodies by solution of the approximate boundary-value problems arising from electromagnetic theory. o.

THE UNIVERSITY OF MICH I GAN 2472-13-T B. Compute the exact radiation patterns from infinitesimal solutions on the surface of simple shapes by the solution of appropriate boundary-value problems arising from electromagnetic theory. (Since 1A and lB are inter-related by the reciprocity theorem it is necessary to solve only one of these problems) 2. To examine the various approximations possible in this problem and to determine the limits of their validity and utility. 3. To find means of combining the simple-body solutions in order to determine the cross-sections of composite bodies. h. To tabulate various formulas and functions necessary to enable such computations to be done quickly for arbitrary objects. 5. To collect, summarize, and evaluate existing experimental data. Titles of the papers already published or presently in process of publication are listed on the preceding page. The major portion of the effort in this report was performed for the Hughes Aircraft Company under purchase order L-265165-F31 under Air Force Contract AF33(038)-28634. K, M. Siegel vi

THE UNIVERSITY OF M I CHIGAN 2472-13-T CHAPTER I INTRODUCTION AND STATEMENT OF PROBLEM This report concerns methods of readily computing radiation patterns, to the accuracy required in practical problems, from slot sources on various perfectly conducting shapes, with particular emphasis on the semi-infinite cone. The methods used are based on the equivalence between a slot source (voltage impressed across a slot) and a magnetic dipole source. The approximations of significance in the applications with which we are concerned are those used for a wavelength limit which is small in comparison to all dimensions of the body in question. The purpose of this study is to determine the behavior of the elementary slot radiators. Such a study is a necessary prelude to the development of arrays of such elements, arrays which will serve as useful flush-mounted antennas. In Chapter 2 radiation problems are discussed from the point of view of known solutions of the reciprocal scattering problems. Chapter 3 then lays the foundation for the calculation of radiation patterns arising from various excitations of a circumferential slot on a perfectly conducting semi-infinite cone. A representative set of patterns, computed on the basis of this chapter, is presented in Appendix C. In Chapter 4 a particular pattern optimization technique, with application to a circumferential slot on a cone, is given. The _______________________i_______ I ii i i ii

THE UNIVE RSITY OF MICHIGAN 2472-13-T material of Chapters 3 and 4 is based on optical techniques; more precise treatments for wedges and cones are given in Appendices A and B, respectively. An alternative approach to electromagnetic boundary value problems: is suggested in Chapter 5 by a discussion of Wiener integral methods. Finally, Chapter 6 presents the conclusions we have drawn to date and offers a prospectus for future endeavor. 2

THE UNIVERS ITY OF MICHIGAN 2472-13-T CHAPTER 2 GENERAL DISCUSSION OF RADIATION AND SCATTERING PROBLEMS 2.1. Reciprocity In order to be able to draw on the large body of knowledge of scattering problems in electromagnetic theory we make use of the reciprocity properties of the electromagnetic field. For our purpose we state the Lorentz reciprocity theorem in the form H1 M2 dv =JH M1 dv (2.1) where H1(H2) is the field due to the magnetization M1(M2) and the integration is over all space. If Ml are'of the form m 2(r-rl 2)' i.e., point sources, then we find Hl(r2)'m2 = H2(rl) ml (2.2) Consider now a magnetic dipole on the surface of a perfect conductor S having a position vector r and a magnetic dipole of the same strength located at some position rl in space (see the following figure). Under these conditions it follows that by varying the orientation of ml, the field induced at r2 on the surface of S by the source at l determines the fieldat r due to the source on the surface S at r. ^~. 2 * )~~~~~~~~

THE UNIVERS ITY OF MI CHIGAN 2472-13-T A r2 r1 If we let Irlincrease without limit, we see that the field induced en S at'2 due to the incoming plane waves determines the radiation pattern due to a source at r2 in the direction of the incoming plane wave. The exact solution of the scattering problem is no simpler to obtain than the exact solution of the reciprocal radiation problem. In contrast, we have readily available a large number of approximate solutions of various scattering problems. We now propose to consider a number of these approximate solutions, most of which are based on some assumption about the field induced on the scatterer, i.e., the radiation field of the reciprocal problem. 2.2. Geometric Optics All of the approximations we will examine are for short wavelengths. We start with the simplest, gometric optics, the exact limit of vanishingly small wavelength. For a finite wavelength the geometric optics approximation is equivalent to replacing the body S on which the source is located at r2 by an infinite perfectly conducting plane tangent to S at r2 with the source at the point of tangency. It is apparent that the pertinent parameters are the radii of curvature at r2 in wavelengths. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ii:- i i -i

THE UNIVERSITY OF MICHIGAN 2472-13-T With these assumptions, the electric field E produced by a radiating slot on a body with voltage Vo across the slot takes the form E Vocurl ikRdl (2.3) 2ir R visible portion of slot where R is the distance between field point and integration point along the slot and dl is an infinitesimal of length in the direction of the magnetic dipole. In particular, when the body containing the slot becomes an infinite perfectly conducting plane, the expression (2.3) becomes the exact solution for the boundary value prbblem of a radiating slot on such a body. The expression (2.3) for the electric field E will depend on the body under consideration only in the sense that it will depend on the position and orientation of the slot on that body; that is, the optics method typified by Equation 2.3 will not yield any information concerning possible diffraction effects due to the body itself. Thus the form of Equation 2.3 will be similar for all bodies whose shapes are in a certain sense similar indeed, when an approximation of the form (2.3) is applied to an arbitrary convex body of revolution having a circumferential slot in a plane normal to the axis of revolution, the problem reduces to that for a cone, tangent to the body of revolution at the slot, with a slot at the circle of tangency.

THE UNIVERSITY OF MICHIGAN 2472-13-T Schensted (Ref. 1) has shown that the geometric optics answer is the exact scattering answer for a plane wave incident along the axis of a paraboloid. Thus, applying the reciprocity theorem, the far field produced by a magnetic dipole located arbitrarily on the surface of a paraboloid is given exactly in the direction of the axis of the paraboloid by the geometric optics field. The expression for the magnetic far field in the direction of the axis of the paraboloid is just twice the free ik(z-zo) space expression for the magnetic dipole (i.e., H 2 eik ( zto A for a magnetic dipole oriented in the i - direction, where z is the distance to field point parallel to the axis of the paraboloid, and z is the corresponding distance to the dipole.) In this same vein, it is worth noting that similar reasoning can be applied to the cone. We recall that the bistatic radar cross-section is very closely approximated by the physical optics formulation for illumination along the axis of the cone. This suggests that the use of the physical optics field in the reciprocity theorem will give a higher order approximation for the field in the direction of the cone axis due to a slot excitation of the cone. In the case of a circumferential slot on either the paraboloid or the cone we note that, of cos n$ excitations, only the cos 0 excitation will produce a non-vanishing contribution along the axis of symmetry. Using the above method for the paraboloid will, as indicated, give the ~.-6 _____6__________

THE UNIVERSITY OF MICHIGAN 2472-13-T exact field along the axis of symmetry for the cosine excited circumferential slot. Although the reciprocity theorem is not obeyed by the physical optics field, Felsen (Ref. 6) has demonstrated that for plane wave illumination along the cone axis the physical optics field agrees with the first order small cone angle approximation to the exact scattered field outside the region of specular reflection 9 <(2o -a). Recently Felsen (Ref. 2) has indicated that this result can be continued past the specular reflection region, and hence, to the surface of the cone, where the agreement with physical optics also obtains. Thus we may use the physical optics field on the surface of the cone. We expect to have reciprocity obeyed by physical optics to first order in the small cone angle approximation and to first order in (ka)'l where a is the distance from the slot to the tip and k = 2r/. In this way we obtain the field on the axis of symmetry produced by the slot on the cone to the same order of approximation. 2.3 Refinements of Geometric optics. In looking for refinements of geometric optics we turn to the reciprocal scattering proolems - in particular to the method of Fock based on an approximate formulation of the scattering problem and to the use of approximations of the exact solution for separable surfaces. The first approach, that of Fock (Ref. 3), is based on an approximation to Maxwell's equations depending upon the physical assumption of a 7 ~

THE UNIVE RSITY OF MI CHIGAN 2472-13-T sufficiently short wavelength. This formulation we can consider as a modification of geometric optics based on a local analysis of the field in the region of the shadow boundary. In particular, Fock defines two universal scalar functions whose argument is a reduced distance measured from the shadow boundary. Depending upon the polarization of the incident radiation with respect to the tangent to the shadow boundary, one of the Fock functions is approximately proportional to the field induced on the scatterer. Following is an account of the general procedure. Let a convex closed surface S, f(x,y,z) = O, be illuminated by a plane wave incident in the direction of the x-axis. The geometrical shadow is then given by the two equations f(x,y,z) =, 0 f 0. Let the origin be located at a point on the shadow boundary with the z-axis the outward directed normal to S and the y-axis chosen to form a right-handed system. Using the geometric assumption that the surface can be approximated by a paraboloid at any point, i.e. z + 1/2 (ax2 + 2bxy + cy2) 0, so that ax + by and the physical assumption that the variation of this field in the z-direction is much smaller than that in either the xor y-direction for sufficiently small A, Fock obtains an approximation to Maxwell's equations which lead to the solutions, _______________________________ 8 _______________________________

THE UNIVERSITY OF MICHIGAN 24L72-13-T H (E) i H eikx F( ) Hz = O on the surface. The incident field is given by o0 0 Ho = (0, Hy, Hz) while the functions G and F have the asymptotic behavior lim G(l)lwhere 5 is a reduced distance from the shadow boundary given by 1/3 =('a) (ax ) by). We have a modification of geometric optics field induced on S which gives a smooth transition through the shadow boundary. Hence, we have an approximate solution to the reciprocal radiation problems as noted above. ______________________________ 9 ~

THE UNIVERS ITY OF MICH IGAN 2472-13-T By way of illustration, the application of Fock's method to the radiation pattern from an axial slot in an infinite perfectly conducting cylinder has been compared with Bailin's evaluation of the exact series (Ref. 4) for this problem. If a is the radius of the cylinder, and k = 2'/A ( the wavelength), ka = 12, so that ka = ) x The corna parison along with the geometric optics result appears in Fig. 2-1. Following the method of N. Logan (Ref. 18) the oscillations in the Fock expression were found by including the contribution going around the rear of the cylinder. Logan has made the same comparison for ka = 8. In detail we substitute for x, which is the correct variable near its shadow boundary, the path length S along the surface of the cylinder. The magnitude of the field is then proportional to 1E - ~G( ) + G( ) where \(k)/V3 S where S is the path length from one shadow boundar, that from the other. 10where S is the path length from one shadow boundary, S that from the other

Geometric Optics 0~ 20 |0 60 80 100 120 = 160 - L9o, ^ 1 1 ~1 ^ ^1 I I I I Cz ~~~~^FIG. 2. I F S1 AC ~~~~~4~ ~~~~~APOXMTO ~~~W~ Fock L'I.4,2.1 0 20 40 60 80 100 120 140 160'18 % IN DEGREES 0 FIG. 2.1 COMPARISON OF SCATTERING AMPLITUDES FROM THE EXACT SERIES, FOCKtS CURRENT DISTRIBUTION, AND THE GEOMETRICAL OPTICS Z APPROXIMATION FOR AN AXIAL HALF-WAVELENGTH SLOT ON AN INFINITE CYLINDER WITH ka = 12

THE UNIVERSITY OF MICHIGAN 2472-13-T In terms of polar coordinates S = a( - ) S-a( - s's a(31r- 0 hence, + aG ((kal/3 () In applying Fock s method to general shapes we must exercise some caution. First, since the method is in fact a modification of geometric optics and therefore a short wavelength approximation, all dimensions of the surface must be large with respect to a wavelength. Second, the radii of' curvature must be continuous. To illustrate this, consider a modified spindle shape. 12

THE UNIVERSITY OF M ICH I GAN 2472-13-T Let a source be located at S. Then, so long as the dimensions a and b are large with respect to the wavelength, the Fock approach is useful in determining the radiation pattern in region 1 but must not be used in region 2 since the radius of curvature is discontinuous at the tip. The second approach depends upon the decomposition of the exact solution into the geometric optics contribution plus the remainder which we will call the diffraction term. This method has been applied in various separable cases for which the exact solution is available. The use of this method to refine the geometric optics contribution depends upon a ready approximation of the diffracted term. This has been achieved by Franz et al. (Ref. 5) for the infinite circular cylinder and the sphere and by Oberhettinger (Ref. 16) and others for the infinite wedge. The diffraction term appearing in the exact solution of the infinite right circular cone, however, is not so easily approximated except in the limits of a large and small cone angle (Ref. 6). Because of our particular interest we consider the cone and its twodimensional counterpart, the wedge, at length. The details of the decomposition into the geometrical optics and diffraction terms appear in Appendices A and B for the wedge and cone, respectively. Except as noted above we do not have much of a hold on the diffracted term for the cone but since we can form an estimate of its size we can discover the range of usefulness of the optics result for sources located at various places on the cone. We find that if the source is a radial distance a from the 13 ____

THE UNI VE RS ITY OF M I CH I GAN 2472-13-T tip of the cone, the diffracted term is of order l/ka as compared with the optics term as long as the source is in view. As we move the observer out of view of the source on the cone, the optics solution is discontinuous, and we must make use of a more sensitive technique. Since we eventually intend to discuss a cylindrically symmetric distribution of sources about the cone axis, we make the point that the effect of the discontinuity in the optics result becomes negligible for a sufficiently dense distribution of sources. In the following chapter a method of using the optics solution in obtaining the entire radiation pattern from variously excited circumferential slots on a cone will be given. ~ ~~~ ~ ~~~~~~~~~~~~~III ~~

THE UNI VERS ITY OF MICHIGAN 2472-13-T CHAPTER 3 CONE RADIATION PATTERNS BY OPTICAL TECHNIQUES In determining the radiation from slots in a perfectly conducting semi-infinite cone by means of the geometric optics approximation we note that the case of radial slots can be subsumed underthe general theory of linear slots in an infinite plane and hence presents no new difficulty. Contrariwise, the case of the circumferential slot introduces a more difficult problem of characterizing the radiation pattern arising from an arbitrary excitation. The case of a cone, 9 = 6o, where 9 is the usual spherical polar variable, with an arbitrarily excited circumferential slot at a distance "at from its tip will be illustrated in some detail. Let (r, 9, 0) and (a, Go0,) designate field and integration points (spherical coordinate system), respectively. Then Equation 2.3 takes the form - Vo eikR A A E c- curl J e. (cos,- i sin ) a sin 90 f( ) d3 (3.1) 27r R ^21 visible portion of slot where i and j are unit vectors in the x and y directions, respectively, Vo is a voltage across the slot, and =r2 + a2 2ar cos os cs 0o + sin sin 90 cos (0-3)] _________________________1^_______________

THE UNIVERSITY OF MICHIGAN 2472-13-T The function f( ) is an excitation, arbitrary to the extent that it is representable by a Fourier decomposition oo f(3) =2 C eing _ i n=-oo where Cm = f() e-3 d/ 7r1 Then for Je9<'r- 90, the region where the entire slot can be seen, and for r>>a, Equation 3.1 becomes 0oo _ ik VO eikr -ika cos 9 cos 90 a sin o0 ein 2 r n=o00 n [ +iTncos (3.2) where 9 and ~ are unit vectors in the usual spherical coordinate system and Sn n [Jn (ka sin E sin 80)- J (ka sin sin 0o) n in1 n- n+ Tn = [n (ka sin e sin o)+ Jnl (ka sin e sin 90) (3.3) When 9 >T - 80, the integration of (3.1) is accomplished by the method of stationary phase; the error thus incurred is comparable to the error of the _________.....___ 16

THE UNIVERSITY OF MICHIGAN 2472-13-T integral itself. Then for large ka, since there is just one stationary phase point at / = 0, E Vo eikr ka sin 90 -ika cos (0 - e) ()). ) ^^o~-, V ~sinO —— ~ ei'f(0)e 8< 9. (3.4) Using the excitation f (1) i, Bailin and Silver (Ref. 7) have calculated from an exact series the expression ika VO ) E-BE (3.5) for four points, when 80 = 165~ and ka = b1T. A comparison between their results and the results from Equation 3.4 is given in Table 3.1. TABLE 3.1 60 = 165~ ka = 5 0' 9 600 75~ 90~ 105~ BE from Ref.7.182 +.087i -.102 -.1101.088 +.lli.108 +.171i BE from Eq.3.5.144 +.098i -.117 - 117i.092 +.134i +.117 17i Figure 3.1 gives a graphical picture, IBEI versus 9, of Table 3.1. We felt that it would be valuable to present additional patterns for various excitations and several different cones with circumferential slots at various locations on the cone. In particular, patterns have been computed for cost (Fig. 3-2) and cos 2/ (Fig. 3-3) excitations with the values of 9o and ka used above. These are presented on the succeeding pages. __________________________ 17 ~~1

THE UNIVERSITY OF MICH I GAN 2472-13-T 90 60 /~~ o ~1525 slot Approximate Calculation x x x Bailin and Silver Calculation FIG 3.1 NORMALIZED FIELD INTENSITY, I BE0I, VS 9 FOR A UNIFORMLY EXCITED CIRCUMFERENTIAL SLOT ON A 30~ CONE WITH ka - 50n1 18

THE UNI VERSITY OF MICHIGAN 2472-13-T 1800 210~ ~ 10.5~0 150~ 00. FIG. 3.2 CONTOUR PLOT OF RADIATION PATTERN FROM A 30~ CONE WITH CIRCUMFERENTIAL SLOT HAVING COS 0 EXCITATION (ka = SOTr) ~19 0.5 /,5 330~ 0.5 300 00 FIG. 3.2 CONTOUR PLOT OF RADIATION PATTERN FROI4 A 300 CONE WITH CIRCtUMFERENTIAL SLOT HAVING COS 0 EXCITATION

THE UNIVERSITY OF M ICH IGAN 2472-13-T 180~ 004 FIG. 3.3 CONTOUR PIljT OF RADIATION PATTERN FROM A 300 CONE WITH CIRCUMFERENTIAL SLOT HAVING COS 2# EXCITATION (ka 507~) ~20 lOo 150 330~ 00 FIG. 3-3 CONTOUR PLOT OF RADIATION PATTERN FROM A 30- CONE WITH CIRCUMFERENTIAL SLOT HAVIG COS 2V EXCITATION (ka = 50r ),, ~~~20....

THE UNIVERSITY OF MICHIGAN 2472-13-T In addition, patterns have been computed for all combinations of these parameters: cone angle 90 = 1600, 165~, 170~; ka = 50r, 75n; and excitations ein, n = 1, 2, 3. These patterns are presented in Appendix C. In the following chapter we present a beam optimumization technique using the above results. 21

TH E UNI VERS ITY OF MICHIGAN 2472-13-T CHAPTER 4 AN OPTIMIZATION TECHNIQUE One method of designing an array of radiators to produce a narrow beam is to fix the power radiated in some particular direction and then to minimize the total power radiated in all other directions. As an example, the case of a single circumferential slot, which is especially simple due to the orthogonality of eina, is studied below. Here the excitation coefficients are obtained directly in terms of the Lagrange multipliers (Equations 4.18). In the case of M circumferential slots we would have to solve M simultaneous equations for the magnitudes of the ein0 excitation in the M slots. ^.e will concern ourselves here only with the forward directions contained in the extension of the cone. In this case, the field produced by a circumferential slot with an e+ino excitation has, to a good approximation, the following simple form (see Equations 3.2, 3.3): _- ikVo eikr -ika cos 0 cos 90 +ino ^A E X r _ a sin 80 e e S + iT cos e 2ff r 0 ^n n where Sn ~i J nl(ka sin sin 0,) + Jn+l(ka sin G sin 00) n in-l. n-1 nl and T -= T J (ka sin 9 sin 9) + Jn+l(ka sin 9 sin ) n n-l - L 22

THE UNIVERSITY OF MICHIGAN 2472-13-T Here the symbols have the same meaning they have had in Chapter 3. What we want to do is to find an excitation 7 ^ (An + iBn) i e n=-N which will give a pencil beam in the direction 8 =, $- =. We proceed by first specifying the field in the direction e = P, 0 = 0, to be ikV eikr -ika cos cos 8 + (^2) E(&,0) ~ ~ _ a sin 0 e (+ ) (4.2) 2rr r This done, we minimize the energy radiated in the forward cone. In order to minimize the radiated energy while simultaneously keeping the field fixed in the specified direction, it is necessary to have a narrow beam pointed in the specified direction. We now proceed to evaluate the excitation coefficients An and Bn in n-1 n-1 accord with the above program. We let Sn = Sn i and Tn Tn i cos 8 so as to be able to deal with real quantities. Now (4.2) gives us the following requirements N T (An + iBn) Sn( ) -l n=-N N (An iBn) Tn(P) = 1 (4.3) n —N 23 ~.

THE UNIVERSITY OF MICHIGAN 2472-13-T The minimum radiated energy requirement means that we minimize f J sin O f Z (An iBn) Sn (0) e c f i s( T(ej do n-N n-N n-N Cn Sin ( )Sn ()+ T n () d (45) 0 ~fidta msN in minimize + A 92 (.7) ~n=~~-N sub t, of curs, to (3 If we n(let n=-N n=-N An =... C p sin (^ (9+ Tn ( 9)) Bn -n= — n-.24 Cn ~~-~ 2 ~* i. _ i..,, _ i! i. ~.

THE UNIVERS ITY OF M I C H I GAN 2472-13-T The Lagrange multipliers A,,/,, v are determined by substituting (4.7) into (4.3). We find /OV = 0 and N S ( ) + N Sn (1) Tn() ( C2 2 E2 n=-N n n =- n (4.8) N S (S ) Tn (') N Tn ( ) or S () Tn () 0 - ~_______n __ L.2 C2 n n 2 N Sk ) Sk (T ) Tk ( ) Zk.... C~ 2 k=-N k k=N C'- 2, J1 E - C2 C2 k~N S k=TN T A ~ n N S2 () N () Tk () (k ) k k=-N k k=-N k 2 k=-N k =k 25 -

THE UNIVERSITY OF MICHIGAN 2472-13-T The derivation given above was for the case where circular polarization is desired in the beam direction (see Equation 4.2). A similar derivation could be just as readily carried out for any other polarization. We have presented above a method for determining the excitation of a circumferential slot so as to obtain a narrow beam. The presentation of the method was based on a particular approximate formula for the field which holds only in the forward regions. In the remainder of the chapter we show that the method holds in general and give general formulas which can use more exact representations of the field (either theoretical or experimental). First of all we need to observe that if the slot has an ein$ excitation then the only 0 dependence of the 0 and 0 components of the electric field will be a factor ein$. We can prove this as follows: Suppose the excitation ein$ gives the following field for 0 0: ein -En(8, ) 4 [fn() + ign(9)] 9 [h + in() +i()] (4.10) Then the excitation ein(0 -') will give the same field for 0 =O since the problems are identical: ein($ - )En(,f) fn(O) + ign(8)] + [hn() + iin() r (h.ll) The field for the excitation e in can be obtained from that for the excitation ein( -c) simply by multiplying by einc. Then replacingoc by 0 we get the desired relation... 26

THE UNIVERSITY OF MICHIGAN 2472-13-T i. ( 0){[ (r 0 ^ [ ( )A A} in0 eikr ein-.PEn(O, ) fn(G) + ign()j + h() + ( in(9)j }e (4.12) Now let us suppose that we require that the field in the direction (a, O) be eikr E (, o)a( + i$) e (h.13) where the excitation is N 7 (An +iBn) ei n=-N We now require that the quantity @o 2r 2 f J sin e E(e), 0) do d9 0 - be minimized. This gives us the condition that 0 21T 2 o o (4.14) 0JJ0 sin e inii (en + iBt) [tn() cons]trn d0 d + 7si J4 (An + iBn) [hn(Q)+ in() ein do d 0 0 be minimized subject to the constraints 27

THE UNIVE RS ITY OF MICHIGAN 2472-13-T N r ZI (An + iBn) [fn(~) + ign( ) ] n —N (4..5) (An + iBr) n(t) + iJn() n=-N Now we let 0 r 2 cl-s / sin 8 f2(e) + gn(e) + hn(e) jn(e de (4.16) o Using Lagrange multipliers,r,4 /, V we minimize Cn (A2 + Bn) -2A [Anfn( ) -Bn gn( -2 )]2 [n gn( ) +Bnfn( (u..17) The result is Xfn(5a) +/ gn( 1 ) +p hn(5) + Ujn( ) A = n C~ n~ (L~. ~(418),Ufn(t) -Xgn( ) + Vhn($) -Pin( ~) Bn 2 28

THE UNIVE RSITY OF MI CHIGAN 2472-13-T or (+ iA 4 )[fn() - ign(~)+ (+ i~ )[hn( ) - in( S An + iBn = (4.19) It is convenient to deal with the following complex numbers Dn An + iBn L - h + i> (4.20) P = p+ iv Fn(9) fn() + ign(9) Hn(e)- hn(9) + ijn(e) Then (4.19) becomes LFn"') + PHn(') Dn C (4.21) n C2 while (.15) taes the form N L Fn( e)1 + PFn(~) Hn() C2....... ~~.. ) n = i. (4.22) n=-N Cn 29

THE UNIVERSITY OF MICHIGAN 2472-13-T Solving (4.21) and (4.22) for Dn we find N2 % N f Fn()l N Fn( F ) Hn() 1Cn -.C2 n=-N n n-N n i 21(~Z I( 02 2 N N F() Hn() N () 2 n-=N n n-N n __ (4.23) Dn =N |Fn(*) 2 Fn(H) HT( r) n=-N n n=-N n n=-N n-.N n This gives the required excitation when circular polarization is desired in the beam direction. If it is desired that the field in the direction (t, O) have the form E (,0)(or,- +/) - (4h.24) where oc andP are arbitrary complex numbers, then the only change in (4.23) is to replace the column....... 30

THE UNIVERSITY OF MI CH I GAN 2472-13-T 11 by ( in the upper determinant. One of the interesting conclusions is that for a beam in the forward direction ('r= 0) it is not possible to reduce the total radiated power to an arbitrarily small value (while keeping the radiated power in the forward direction fixed) when using a single circumferential slot, since for?= 0 we have Dn = 0 unless n = + 1.,i~ ~ P

THE UNIVERSITY OF MICHIGAN 2472-13-T CHAPTER 5 WIEdER INTEGRAL METHODS Because standard techniques of solving scattering and radiation problems have proved inadequate for some problems of interest, there has been a certain amount of interest in developing new approaches. One such new method, the Wiener integral method, is discussed here. It is not yet clear whether or not this method will prove practicable; but the possibility that it might, may justify this discussion. The Wiener integral method actually is a method of attack for boundary value problems involving the diffusion equation. Since our interest centers in the wave equation rather than the diffusion equation, we will start the discussion by pointing out a connection between the two problems. Such a connection is mentioned, among other places, in Reference 8. Let, (x,y,z,t) be the solution of the equation v2s 4 U -4f(x,y,z) S(t) (5.1) which is zero for t< and which satisfies a linear, time-independent boundary condition on some surface. Here S(t) is the Dirac delta function. Now consider the function G(x,y,z,p) e P2t(x,y,z,t) dt (5.2) -oo 32

THE UNIVERSITY OF MICHIGAN 2472-13-T If we differentiate under the integral sign we find that VG - p2G -4f(x,y,z), (5.3) and that G satisfies the same boundary condition asZ/. The equation for which we would actually like a solution is V2G + k2G = -4f(x,y,z). (5.4) If we assume that G(x,y,z,p) is an analytic function of p, then we can replace p by ik or -ik to get solutions of (5.3). In order to proceed in this way it will be very desirable to obtain G as a function of p rather than obtaining G numerically for various values of p. The reason that we do not take p = I ik directly in (5.2) is that then the integral would not, in general, converge, Although we have mentioned only a scalar problem, it is evident that the above reasoning applies equally well to vector problems. Having made the connection, we will now restrict our attention to diffusion problems. In order to see the relationship between Wiener integrals and diffusion problems it is simplest to start with the freespace problem. In this case it is well-known (and easily verified) that the solution of (5.1) is 8~(x a z ~t) e=J t ~ f(x-;,y-,z- g) d5 d d. ^(x1y~zt) r _ fqx )/t)+d7d ~-OD -00 (5.)3 33

THE UNIVERSITY OF MICHIGAN 2472-13-T We can write this concisely in terms of Wiener integrals by using the following relationship w O x2 J [x(tSj dwx= J (x)dx. (5.6) c -0o This equation, as well as a brief discussion of what a Wiener integral is, will be found in Reference 9. The result of using (5.6) in (5.5) is that we can write w w w WWW 8(x Y zt) 2 ji f[x^- (t) (t )(t) z-gt)j d c d 5 d d d c c c (5.7) In order to see how to proceed in the case when a scattering surface is present it is convenient to give a physical interpretation to (5.7). The interpretation we use is the following. At time t = O we release particles with a density f(x,y,z). A typical particle follows a path which at time s has coordinates x- (t-s), y- (t-s), z- (t-s). The probability of this particular path is dw; dw4 dwS. Then bl(x,y,z,t) is the density of particles at time t. Using this type of interpretation, we can conceive problems whose solutions in terms of Wiener integrals are obvious, but which are at the same time solutions of (5.1) subject to certain boundary conditions. For 3h

THE UNIVERSITY OF M I CH IGAN 2472-13-T example, suppose we want bY (x,y,z,t) to vanish on some given surface. We can accomplish this in the following way. We again let Ab be a density of particles. We again start the particles out at a density f(x,y,z) and let them follow paths x - (t-s),.... However, now whenever a particle hits the given surface we assume that it is absorbed. The density of particles will satisfy (5.1). On the surface the density of particles will be zero since all particle s on the surface are absorbed. Thus the particle density is a solution of the desired problem. But, from the description of what happens to the particles, we see that the density of particles can be written as w w w (x.y.z.t) X( i, ) fJx- (t), y-(t), z- t(t)ddwdw c cc (5.8) where X is zero if the path x-. (t-s),... touches the given surface and is one otherwise. This is a well known way of satisfying this problem and is mentioned, for example, in Reference 8. The above model is not the only one which will enable us to solve this problem in terms of Wiener integrals. The following is another of the possible models. We have two kinds of particles which we call positive particles and negative particles. We start with positive particles having a density f(x,y,z). Now whenever a particle hits the given surface it is reflected off specularly and simultaneously is changed into a particle of the other kind. If we look at the particles near the surface we find that 35

THE UNIVERSITY OF MICHIGAN 2472-13-T half of the particles have just hit the surface and thus changed sign so that the numbers of positive and negative particles are equal. If we take~ to be the difference in the density of the positive and negative particles we see that 2 is the solution to the desired problem. In terms of Wiener integrals the solution to this problem can be written immediately in the form w w w s= J f(xiy*,z*) dS 4d d. (5.9) c c c Here the path x-; (t-s),... is modified by reflection as indicated above (see Figure 5-1), n is the number of reflections the path undergoes and *, y*, z are the points from which the modified path starts. We see x-9 (t),y-7 (t),z-r (t) x,yz Scattering x*,y,z Surface Unmodified Path Modified Path FIG 5-1 PATH OF A PARTICLE IN FREE SPACE (UNMODIFIED PATH) AND IN THE PRESENCE OF A BODY (MODIFIED PATH) 36

THE UNIVERSITY OF MICHIGAN 2472-13-T that (5.8) and (5.9) give us two quite different expressions for the same thing. We could also construct other models which would give us still more expressions. For example we could think of a surface which absorbs some particles while reflecting and changing the sign of others. Also the reflection would not have to be specular. Other models might involve the creation of a number of particles whenever a particle hits the surface. All of these different models give us many Wiener integrals all having the same value. Now only in very rare cases will we be able to evaluate the integrals in a simple closed form. Generally we will have to resort to a series expansion or something similar. Thus the many different models give us some freedom in trying to choose a model for which the expansion will be rapidly convergent. Also it is possible to use the results for a second model as a check on the results of the first model. Equations 5.8 and 5.9 have been applied to scattering by an infinite plane. A sequence of approximate evaluations of the Wiener integrals was used. (Cf. Equation 5.13). The approximate evaluations of (5.9) could all be carried out in closed form and all agree with the exact answer. The first two approximations to (5.8) were obtained. Figure 5-2 illustrates the relation between the approximations and the exact answer as a function of distance from the plane at a particular instant of time. 37

THE UNIVERSITY OF MICHIGAN 2472-13-T O..06- / \~ 1 Approximation to Eq. 5.8 // 2d Approximation to Eq. 5.8 0~~.4 t / ~ ~~~ Eq. 5.9 0.22 0 0o.5 1.5 1 2.0 2.5. 3.0 3.5 Distance From Plane FIG 5-2 COMPARISON OF EXACT AND APPROXIMATE INFINITE PLANE SOLUTIONS The above discussion has been for a scalar problem. We will now mention a model which may be used for electromagnetic problems where the bodies are perfect conductors. As a preliminary let us reconsider the model discussed above which involved two kinds of particles, positive and negative. A slightly different language to describe the situation is found to be convenient. We now assume that we have only one kind of particle, but that the particle carries along a label which can be either positive or negative. Whenever the particle hits a scatterer it is specularly reflected and its label changes sign. Now to compute,/ at a point we take a small volume, v, about the point and let =ZL,n, (5.10) V ^~f n

THE UNIVERSITY OF M I CH I GA. N 2472-13-T where Ln is the label ( +1 or -1) of the nth particle in the volume. In the vector case we proceed similarly, except that now the label is a unit vector. If we are trying to compute the electric field, we require that on reflection of a particle the tangential component of its label change sign while the normal component remains unchanged. The formula for lis now - v i E n* (5.11) We could equivalently use a model involving six kinds of particles corresponding to positive and negative components in three directions. The above considerations allow us to write down (Wiener) integral representations for the solutions to scattering and radiation problems involving arbitrarily shaped bodies. However, in order for these representations to be useful we must have practical methods for evaluating the Wiener integrals involved. This is the problem which still awaits a completely satisfactory solution. We will discuss below an approach which can be used but which is quite laborious. When one is faced with a Riemann integral which he cannot evaluate in a simple closed form, one frequently resorts to numerical integration. Certain general formulas have been developed for this purpose such as the 39..

THE UNIVERSITY OF MICHIGAN 2472-13-T trapezoidal rule and Simpsonts rule. For Wiener integrals Cameron (Ref. 10) has given some analogous approximation formulae. The simplest formula he gives is w CO CD 2 2 n FF[X.)]dxn^ f...fe n F J'1sin(j.l)7r() -o n/2 =l 2j-1 d -co -aO (5.12) Thus the Wiener integral is approximated by an n-fold Riemann integral. Cameron has shown under certain conditions that the right side of (5.12) approaches the left side of (5.12) when n —ao. Let us take a look at the application of (5.12) to (5.8). In order to use (5.12) we find that we must evaluate the functional for a path which is given by a trigonometric series with arbitrary coefficients (we integrate with respect to the coefficients from- oo to co). This means that, for (5.8), we must take a path whose x, y, and z components are trigonometric series in time with arbitrary coefficients, and find out if the path hits the body or not. We then integrate the coefficients over the region in which the path does not hit the body. Determining whether the path hits the body involves the solution of a complicated transcendental equation. The situation would be much simpler if the paths we hnad to consider were of some simpler form such as broken line segments. We can readily obtain a modification of (5.12) 40

THE UNIVERSITY OF MICHIGAN 2472-13-T in which the paths considered are broken line segments. We quote the formula without proof: /F[x(.i)]d^y F[L Z (J7)d (5.13) e d M^j=l (' /2 df''T dfn n where E.j t.(t) is the function shown in Figure 5-3. j=l r ~ ~ t rn 0 1/n 2/n 1 t FIG 5-3 TYPICAL PATH When we use (5.13) in conjunction with the model for the vector problem discussed in connection with (5.11), we find that for n = 1 we get just geometric optics. For larger n we get corrections to geometric optics including contributions from the shadow regions. 41

THE UNIVE RS ITY OF MICH IGAN 2472-13-T The sequence of approximations obtained is not an asympototic sequence although the greatest accuracy is probably obtained (for a given n) for short wavelengths. If n is taken sufficiently large, we can get as accurate an answer as desired, but the amount of labor involved in evaluating the multiple integrals is very large. Thus it would be desirable to get improved integration formulas. 42

THE U N I VER S I TY OF M I C H I GAN 2172-13-T CHAPTER 6 CONCLUSIONS AND PROSPECTUS In conclusion we find, on the basis of the above analysis (see Chapters 2,3,4),that radiation patterns of single circumferential or radial slots with arbitrary excitation can be given by optics techniques with sufficient accuracy that we may now turn to the problem of using such elements to produce a useful antenna. The solution of this problem requires that we determine for a given surface the distribution and excitation of such elements necessary to produce a radiation pattern of given beamwidth, side-lobe level, gain, and scan capabilities. Hence, during the coming year, The University of Michigan plans to investigate the following problems: (1) the problem of determining pattern features, corresponding to scan angles 00 - 80~, resulting from (a) a linear array along a cone generator, and (b) a system of such linear arrays spaced at certain azimuthal intervals and excited with a constant, or cosine, azimuthal distribution. In connection with this problem an attempt will be made to ascertain the degree to which the control of beam shape can be specified in terms of phase distributions (possibly non-linear) along an array. (2) the problem similar to 1 and for (a) an annular array, and (b) a set of such arrays. h3

THE UNIVERSITY OF MI CHIGAN 2472-13-T Problems 1 and 2 will undoubtedly necessitate refinement and extension of the approximation methods of this report so as to obtain beam scan information to within prescribed limits of beam deterioration and cross-polarization effects. A study will be made of (3) the problem of determining what alterations in pattern features and scanning capabilities of a dipole-activated prolate spheroid result from alterations of the surface of the spheroid. Surface alteration would be defined to mean placement of dipoles on the surface in various ways in possible conjunction with bumps and indentations of the surface itself. The Polytechnic Institute of Brooklyn and The University of Illinois are engaged in efforts complementary to those of The University of Michigan. These efforts are described below. The Electrophysics group of the Microwave Research Institute, Polytechnic Institute of Brooklyn will continue to investigate the problems of scattering by and radiation from infinite and finite cones. In particular they will continue their study of the vector field problem of radiation from circumferential and infinitesimal slots on a cone via Green's function and modal techniques; in this consideration they will give particular attention to alternative representations which would permit casting the results into a rapidly convergent form applicable to the numerical description of fields in illuminated, shadow, ~ u~ ~~~~~~~W

THE UNIVERSITY OF MICHIGAN 2472-13-T and transition regions surrounding various conical geometries. This group will investigate fields set up on coupled radiating slots and slot arrays on cones via integral equation and variational techniques; from this information they hope to apply the Green's function results alluded to above and thereby obtain radiation patterns from slot arrays of arbitrary elements. In connection with arrays of slots, they will investigate the ttsurfacett and "leaky" waves whose propagation along the surface of a cone, in consequence of slot arrays thereon, and whose reflections by the cone boundaries,can modify the radiation pattern. The Electrophysics group will attempt to apply the above results to numerical calculation of radiation patterns for specific array geometrie s. The University of Illinois will analyze a method of replacing an antenna system which can be enclosed by a given surface in space, by an array of slots in a conductor which coincides with the surface. In particular, they will attempt to determine the minimum density of slots on the surface of a sphere (and ultimately on a finite cone) necessary to reproduce the pattern arising from a continuous tangential field distribution (or a dish in free space) according to variously prescribed criteria.

THE UNIVE RS ITY OF MICHIGAN 2472-13-T APPENDIX A WEDGE The radiation from two slots parallel to and equidistant from the edge of an infinite perfectly conducting wedge is treated below: the exact boundary value problem is solved, and then, since the solution is not amenable to rapid computation, an asymptotic approximation as well as the geometrical optics solution is developed, and the results obtained from the exact and approximate methods are compared. A more general solution for the radiation from a wedge, involving slots with various orientations relative to the edge, has been studied by Felsen and is given in Reference 11. A-1 Electromagnetic Field for an Infinite Perfectly Conducting Wedge with Uniformly Excited Slots Parallel to its Edge Let Vo — voltage across the slots a - distance of slots from the edge e i time dependence of fields. The electric field for the slotted wedge of Figure A-1 is given by TM modes, with the electric field perpendicular to the edge of the wedge. 46

THE UNIVERSITY OF MICHIGAN 2472-13-T z _ I x Direction Direction to of Incidence Field point FIG A-1 A-l.l Exact Solution The exact solution is obtained from the boundary problem (V2+k2)y= 0 (A.1) where E 1Vx ~ L laY And (AA2) E = --- x k -- ik ik r rwith boundary conditions 1k 1r dY |_ Vo8 () A.3) since the tangential fields must vanish on the wedge, except at the slots. Solutions of Equation A.1 are: 47

THE UNIVERSITY OF MICHIGAN 2472-13-T = Jv (kr) Avcos2/ + Bsin2j (A.4) for r small, so that the fields remain bounded near the edge, and Y'r= (kr) Acos V+Bsin (A 5) for r large; this corresponds to an outgoing wave. From Equation A.3 d s o at ~ V andV, Bn= O for n even. An 0 for all n since the V component of the electric fields as the edge is approached from above the 0 O axis must be the negative of the ~ component as the edge is approached from below the 0 = O axis. Then oo Y-I B2n-1 H(1) (kr) sin 2-) 0 for r large, (A.6) n-l n and oo Y B2n1 J (kr) sin ( 2n- ) for r small. (A.7) n=l 2n-1 2n- 2 The coefficients may be found by use of the Lorentz Reciprocity theorem in the form (Ref. 7): O x ] ~ E n x HdS 4~ __

THE UNIVERSITY OF M I CH I GAN 2L72-13-T where E, tare the fields obtained from or given by Equation (A.6), and EM, HM, are the modal solutions corresponding to a given index n. S is the sum of SO, S, S2 and Soo as shown in Figure A-2. / FIG A-2 REGIONS OF INTEGRATION FOR THE WEDGE Integration over the surface yields: 2S =0 = II + 2 Vo(-1l) J (ka) m = 2n-1 (A.8) Using Bm from Equation A.8 in Equation A.6 and applying Equation A.2: iV0 ar E n (ka) H(1) (kr) sin 2n- (A.) __________________________ )49 A*/O LY —

THE UNIVERSITY OF MICHIGAN 2472-13-T for r large. For r small the Bessel and Hankel functions are reversed. A-1.2 Asymptotic Solutions and Estimate of Error of Asymptotic Solutions Approximate expressions for Ep for r>a and ka>l are obtained through three approaches. The first approximation may be obtained by considering a related scattering problem. The second method consists in expressing Equation A.9 in terms of integrals obtained by Oberhettinger (Ref. 12) and by obtaining asymptotic expansions of these integrals. A comparison between the exact solution (in terms of the integrals of Reference 12) and the approximations used provides an estimate of the error of the asymptotic approximation. A third method expresses the integrals encountered in the second method in terms of continued fractions at large distances. The same expression is obtained from the first and second methods; a different and more accurate one from the third method. Method 1: Asymptotic Solution Obtained by Considering a Scattering Problem a) Development of Asymptotic Solution For r>>a, Equation A.9 may be expressed as:.2 eikr+ ~ 0 -i?(n+),2 — E- i, e ~ ik + - e X J (ka)sin (n+)1 2'r 22 e2r2 JrJ (A.10) 50

THE UNIVERSITY OF MICHIGAN 2472-13-T The expression in brackets may be obtained from a two-dimensional scattering problem. For this purpose consider a perfectly conducting, non-slotted wedge in the polar coordinate system, a, 0, as shown in Figure A-3: Y x A-, FIG A-3 For this wedge geometry the geometric optics fields are: -ika cos(0-o0) YrG.o.(o) = 2e u(7+* 4o) (A.ll) and -ika cos( )+ 00) Y/G.,(0Ro)= 2e u((-0o) (A.l2) where u is the Heaviside unit function. The scattered fields corresponding to the above geometric optics fields are (Ref. 13): oo~ -in 2 37r -ie J (ka)(1-4 6 )cos nf4o) (A.13) e nr no5 ~~ 51 ______________

THE UNIVERS ITY OF MICHIGAN 2472-13-T and 2 OD nIT n7(n 2 X e i' n^ (ka) ()n(l ) cosA.) 7o n=O From this it appears that the bracketed expression in Equation A.10 is (0) - (-$o). The total geometrical optics solution G.O. given by the sum of Equation A.ll and Equation A.12, may be expanded into the set of functions sin(n + -):; (n +) 1 5'G =.On fn(ka) sin (n 0 (A.15) where, for ka large, n 2 Tr ika 1T 2 Mika +~] fn(ka)( (_l)n 2 [e ika I +2 e. (A.16) a~i7~ ~?'Le'"' i) Tcos(n+. + )e By separating out the geometrical optics solution expressed in Equation A.15 and Equation A.16 from Equation A.10, and by using the asymptotic form of fn(ka) and J( T (ka), and then summing the terms not involving G.O. ikr + ika + (A iVo e nH e 2IT......... + (A.17) ^"2~" ~. ^ 0. 2~. ~ ^2 ~

THE UNIVERSITY OF MICHIGAN 2472-13-T 1 - Tr(Tr*~) -a (Tr-0) { I cos cos 1 1 This asymptotic expression may be obtained through other methods (see Section A-1.2, Method 2a). It is not applicable to the transition regions $= + (~T- o) for which an asymptotic formula must be derived separately. b) Comparison of Scattering Problem with Results Obtained from the Pulse Solution of Keller and Blank The bracketed expression in Equation A.17 may be obtained by considering the solution obtained by Keller and Blank (Ref. 11, pp. 75-94) for the two dimensional scattering problem of a plane pulse incident on a perfectly conducting infinite wedge. Again considering the polar coordinate system a, 0 and letting V(cC, 0, t) be the pulse solution, the time harmonic solution is given by Duhamel's theorem (Ref. 14, p. 90): oo i t iwt ) i(t) e dt (A.18) -oo and oo V(t) -y e dt (A*19) -oo00 ~ ^ ~^~~~~~~~5

THE UNIVERSITY OF MICHIGAN 2472-13-T The pulse solution corresponding to the geometric optics field -ika cos((-0o) -ika cos(0+ )) e.o. 2 u(2+V )-0 e u(- -o) (A.20), V GO u t+ cos (-o)j u(\+ Vo )-u [(\tSAC cos(V+001 u(r 0) j (A.21) The diffraction field obtained from Keller and Blank (Ref, 14, p.81) is: j Ds r[(1- 2) sin Xw V'tan-'.. ~.~ (1 + p )cos X + -2p sin X (l + 2cos T+ 2 f sin), where a = ~- n and the tan-1 lies between ct + c2t2 -a2'' 0 and T. The diffraction field may be obtained by applying Equation A.18 to Equation A.22 00 r: ibWS e -ftanl Ja )2 si (l1+g )cos>Xr-2X sin I (A.23) -tan-l (1- 2A )sinXir dt i(l+ ~ 2)cos r 2'sinX ) J 54

THE UNIVE RS ITY OF MICHIGAN 2h72-13-T The integrals may be evaluated through rapidly convergent expansions (Ref. 15), of which only the first term is used here. By examining the value of the integrand for t a + * ( e small,) it is observed that, C except for $ + (r- 0o), the integrand varies as. Constant terms, however, are brought about by values of tan'l close to;hence,to avoid this, the tan"1 will be chosen between - and. This will be indicated by t Equation A.23 is then in such a form that the rapidly convergent expansions mentioned above may be used. 2e [u ( 00+) -u(-o00-0)] 4/ 2 2 (1 2) sin A,rr 2~I~i 2 ~,2) a/ l ^ (I + e2^) cos Xt - 2 se sin- (1 - p2 ) sin A IT ~~-t^r^^~^~~~ ~dt. (1p2N) CcosXi + 2e sinX j The field of Equation A.21 persists only until t = a; the terms c introduced by the shutting off of Equation A.21 cancel the first term of Equation A.24 and the scattered field becomes (by using Equation 1 of Reference 15): ~ ^ ~~~~~~~~~~~~~~s

THE UNIVERSITY OF MICHIGAN 2472-13-T e ika )1 (1 ^ ) sinAir G.. + 31 iL \ +,0 2A )cosAIT -2/ A sinA0 (A.25);1 (1-p2A ) sinT A1I (l+p )cosAqT +2p' sinA$ As ka -oo, this may be reduced to the expression between the brackets of Equation A.17. Method 2: Asymptotic Expressions Obtained by Use of Oberhettinger's Formulas for Diffraction by a Wedge. Estimate of Error. a) Asymptotic Expression Excluding the Transition Regions, and Estimate of Error. Oberhettinger obtained the Dirichlet Green's function for the two dimensional case of an incident cylindrical wave with direction of incidence L and axis parallel to the edge. Y Direction to Field Point 00,/ / ^ Direction of,, ~. Incidence x FIG A-4 56

THE UNIVERS ITY OF MICHIGAN 2472-13-T The expression is: G( (1) 0o0o (1) ( G(^ k) = i[JO(ka) Ho (kr) + Jn(ka) Hnl(kr) (A.26) L nl 2j Ij 5Co ($ -) + Cos o (0 + S+ 200) 9, from expression (A.6) with the proper coefficients, may be expressed in terms of this Green function. Then kV0 [G(o0)- G(-00)] (A.27) By use of Equations 29 - 35 of Reference 12, this may be written as: kV (1) 2H 2H() [k(r2+ a2 - 2ar cos(0- 0o))J u(T+ - 00) 2 2 1/21 - 2H(1) k(r2 + a 2ar cos(% + 00)) u( r- - 0 ) (A.28) H ) k(r2 + a2 2ar cosh x) 1/ dx + 1 cos Ir (it+ r) coM...._,J ~ 2) cosh X - sin I (Jr+ $) (This equation is continued on next page). - ~ -~ ~5?

THE UNIVE RS ITY OF MI CHIGAN2h72-13-T H) [k(r2 + a2 + 2ar cosh x)/ ] dx - 1 cos h ( _ /+ H) k(r + a2 + 2ar csh x) ] dx''; / cosh 7x - sin ) (- _ ( ) ~o 220Zo ^o H(1) [k(r2 + a2 + 2ar cosh x) 1/2] dx o o Estimate of Error An asymptotic expression as well as an estimate of the error was found for integrals of the type 0 4) (kR) ____________ (k)dx (A.29) coshlx + sin/ x o wheref = - f+ )0) and R2= r2 +a2 +2ar cosh x.'or r>> a HO () (kR) eik - = A1) (kR). (A.30) 58

THE U N IVE RS I T Y OF M I CHIGAN 2472-13-T Equation A.29 may be written as 00 (I) o H(o) A0 (k) d 1 (kR) - A o R) dx + D. (A.31) J 7lr + 7 r cosh + sin cosh + sin 0'0o o 0n I may be evaluated by using a method which is elaborated in the following section:00 co I=3 A4l)(kR) dxT i -i ekRd.......... _ e -' coshlrx sin J, f [cosh sin3 ] e - -i ~ 1 e e......................... (A.32) (1 sin/5 ) 3 2kar' I P1dR-a r+a ikR. ikJ ix' ar sinh x(cosh 1 + sin3) o L _I 1 dx j2kar(R-a-r) (1 sin6) J I 1 J eikR cB dx sp ik I do = Isp - J59

THE UNIVERSITY OF MICHIGAN 2472-13-T 00 (1) T H (kR) dx Then I = IS p J + D J cosh rx + sin p ik (a + r) SP k fk'~ (1 I sin/) H(l(kR) dx By using S co = Isp in Equation A.28, (A.33) cosh 7+ sinp r ak V0 H) [ k r2 + a2 _2ar cos (000 1u(f+ 0 00) (A.34) 1 k Vo H1) [k Jr2 + a2 - 2ar cos (00o) u(i - 00) V eik (a+r) I 1 0o0 cos I ^ cos r (7,- The asymptotic form of E obtained from Equation A.34 corresponds to Equation A.17. 60

THE UNIVE RS ITY OF MI CHIGAN 2472-13-T In order to estimate the error of the approximation, bounds will be found for J and D. From the results of the following section jI I< 2(n+l)A (A.3) -.i~ ~i~ (A.35) where n is the number of zeroes of - and A is a bound of B. ax It was found that IBL| 1 (A.36) 4a TWP (1+ sin6) which implies 1 I1k 31~~ (A.37) 2k32 a Ar (1+ sin /3) Only one zero was used for B in Equation A.35, since the other zeroes ax occur at the end points where B = 0. A bound for D is obtained as follows: 0o (A.38) 77 rF(v3/2) for z,3/2 ~~Then Hi~~~~o (kR) A kR ok (1Then () A(1)( 1) 1 TX (A.39) o o (kR) (k (_61 61

THE UNIVERSITY OF MICHIGAN 2472-13-T and using this in D gives a bound:.00 (1) (1) o os 1 (kR) aJ r - A) 3 () + )csc (+ ) cosh - sin/ General Discussion of Error in the Stationary Phase Method The object of our attention here will be the integral I = g(t)eikf(t)dt 0 where f'(t)> 0 except for t = 0 and f' (0) = 0. By a change of variable the integral can be rewritten in the form f(oo) I = | g(t) eikf df. ft(t) f(O) Near t = 0 we have +t3 f(t) = f(O) + f"(0) t + f'(O) +... 2 6 f'(t) = f, (O)t +... =2f'T(0) o[f - f +.. Making use of this behavior we split the integral into two parts, the first of which is the stationary phase contribution, as follows: 62

THE UNIVE RS ITY OF M I CH IGAN 2472-13-T f(oo) ikff(oo) ^I [S g(0) eikf df g( t) g() ikfd f(0) f (}f -f(0] f I(t) efd(0)Lf-f(0 On integrating the second integral by parts we get I = I + g(t) g(0) eikf(t) Sp ft. i Lf'(t) 42f"(0) [f(t)-f(0) ]' ik 0 00 I | eikf(t) TTT J2f (t) () f (O) g(0) g0(o)' ft'(o) - eif(0) 0 +ft() g(O) 3f"(0) ik - 1 ekt)d g(t) g(O) 0P f fI,(t) 42f" 0 (0) [f(t)-f(0)]J Now by obtaining a bound for the last integral, J, in the above equation we can obtain a bound for the error incurred by using stationary phase. Such a bound is 63

THE UNI VE RS ITY OF MI CH I GAN 2472-13-T 1 a 0k L dt 42f,(O) [f(t)-f(O )] 0 The integral just above can be rewritten in the following form tfl II.17 i'g(t) g(0) j+1 IJI t) 42f1 (0i) f (t)_-f(O)]- where to = 0, tn+1 = oo and g(tj) g(O) (t - 32f,(0) [ f(t ) -f(0)] Thus, if we can get a bound of the form g(t) _ g(0 A, If (T J,,'2f (O0) [f(t)-f(O)] then we have IJI - 2(n+l)A k64 64

THE UNIVERSITY OF MICHIGAN 2h72-13-T b) Asymptotic Expression for the Transition Region For D = (r- 00) the asymptotic expression (A.17) is not applicable. For these values of 0, the integrals of Equation A.28 are of the form: co ~1 Ang 4 H(l) (kR) dx (A.41) I ($,2 o =- sin 0Io (A 41. cosh-I. X- cos V where ^ has the values: 81 =+ - 0o; 2=- 2 - O + 0o; 3 - - o0 0; 6' + 0+ o Using the method or the preceding section, these integrals are given by: ik(a+r) t r 2I (,200) 2. - (A.h2) 2o k JZ For 6 and 8 there are no singularities. For and S3 there are singularities at 0 = + (W- o0). For such integrals, Oberhettinger (Ref. 16) obtained a series representation of which the first term is used here. Oberhettinger's integrand numerator was e-ika cosh x; for this reason i was changed into -i and the factor eikr " was added to Oberhettinger's results in order to obtain the integration of the asymptotic form of Equation A.41. 65

THE UNIVERSITY OF MICHIGAN 2472-13-T Then k ikr- i -i ka cos( 0o)- uo) ^-^ JL e ^ e - eU(ir + 0 -0o)' 2 kr -ika cos($ + 0) -e U- u- (o) (A.43) + T*( (1,2o0) + I( 62,200) - * ( 3,2o0) - I( 8,20) where the I terms are given by Equation A.42 and I* by iika cos- -", I* (8,2$00) ~ i ei{ eika C058 mS* (2J 2ka [sin 8|)sgn 8 tl'' ir i (A.44) + e e Ao( &,200) Ao (&,2r) - Ao (21-& 2, ) where 2D 2 Z it2 S *(z) J e dt and Ao( 2,20) = cot z This leads to an expression which is bounded for 0 = ( (- -o), and which for +- (i- 0o) and ka large, gives a field E which reduces to Equation A.17. 66

THE UNIVERSITY OF MICHIGAN 2472-13-T Method 3: Asymptotic Expansions by Means of Continued Fractions. A different asymptotic form of Equation A.28 may be obtained by expressing integrals of the type of Equation A.29 in the form of asymptotic series, and then converting them into continued fractions. Consider JS 00 sin Tr- e-ika cosh x I(S,200) = sin.fcoshfX - cost V (A.15) These integrals are of the type of Equation A.29 for r>> a if i-.-i and the factor 2 eikr- ir/4 is added. irkr The asymptotic form of A-45 is: -i(ka + Ar(S,2$0)p~n + I( s,200) ef' n< nze0(ika) n.(A.46) where sin 2~0 sine= An(,200o ) tn (I t,< ) 2(200) t+2 cosh(f cosh1 (l+t)-cos I nO 67

THE UNIVERSITY OF MICHIGAN 2472-13-T This asymptotic series may be converted into a continued fraction (Ref.17) cot i~ -i(ka + 3) _7o ~ka e 4 I( 8,2 ) 2 1~ (A.47) 1 12 + 2w2 (sin w )-2 + ika +. 1 1 -02o J.. There are indications that the continued fraction may converge and represent I ( 8200), and further that the continued fraction, when cut off at a point corresponding to that at which the asymptotic series is cut off, will always give a better answer than the asymptotic series. The fields obtained from this method were calculated and compared with those obtained from t he exact solution; this method was found to give a better approximation than Equation A.17. A-1.3 Geometric Optics Solution By adapting the field of a slot in an infinite perfectly conducting plane to an infinite perfectly conducting wedge, one obtains, for the case of short wavelengths, a first approximation to the radiation problem: this is a geometric optics approximation. The field of a slot of length L and width W on an infinite perfectly conducting plane (see Fig. A-5) is derived as follows. Let L <<X and 1 <<X. The tangential component of the electric field must be zero outside the slot and equal to V/WT across the slot. Such a field may be produced by a magnetic dipole placed in the plane of the slot along the slot: its magnetic moment is derived below. The field of a magnetic dipole at x=O, y=0 oriented along the x-axis is ~' K 1~ (1 - i;) rx n where K must be determined. I~ — ~p - ~ 68

THE UNI VERS ITY OF MICH I GAN 2h72-13-T Y K ~ L w V FIG A-5 SLOT CONFIGURATION At a distance A above the slot, for kr << 1: iK(xi. + y + A k)xi iK( A j- yk) k( + x2 y2 +2 (.28) A As A-O0, the j component should vanish except at x y 0 O. In order to evaluate K, the tangential component is integrated over x and y. For the slot the integral has the value - L Vo j, where the direction of the field has been taken to be given by the right hand rule. Then 00 00 - V = l ikA fdx dy = 2iiK ^ ^ ^ ^ [ f'-~ — I a.49) A O0 k J T = 2 + A2 k 00-oo -00 and ikr ikr ikVoL e ^ VoL e EI __L __ (1I- 1) r x i V x ~ ~ (A.50) 2w r 2 r This gives the exact field in the case of an infinite plane. when Equation A.50 is used for an arbitrary body, integration takes place over the portion of the slot visible from the field point (Fig. A-6). 69

THE UNIVERSITY OF MICHIGAN 2472-13-T Yx FIG A-6 For a slot parallel to and placed at a distance a from the edge of a wedge oo iJr2 - 2ar cos (-)+ a2 + E v xf k dz u( +0-0o) 2m J -r2 2ar cos(-00) +a2 +z2' which becomes ~ ^ ik[r - a cos(-00)] - E Q V0 2 w - e - u(*r+$-o) for r>> a. (A.l) This is in agreement with the geometric optics portion of Equation A.17 for f (0o) G.O. A-1,4 Comparison betweenthe Exact and Asymptotic Solution and Calculation of the Error of the Asymptotic Expression. Fields for a0 = and 0= I 5 and - were calculated from the 2 and ~o =, - exact solution (A.9), the asymptotic solution (A.17) and the asymptotic expression obtained through the use of continued fractions. The fields from Equation A.9 and Equation A.17 are plotted in Figure A-7, and the 70

THE UNIVERSITY OF MICHIGAN 2472-13-T fields from A-9 and those obtained through the use of continued fractions in Figure A-8; for the latter calculation two terms of Aw48 were used. From the graphs it is apparent that the agreement is much better for the approximation obtained through continued fractions, The approximation is also best for thin wedges and deteriorates as $ decreases. This may be expected from the estimated error of Equation A.17: the calculation was made through the use of formulae given in Section A-1.2 Method (2a) and results show an increase in the radius of error from 0o = to 00o = The equation for the approximate radius of error is: Q OkaI Cos c os Ico l where From ( gr t n = to From these graphs it is apparent that the error increases from O I r to o 3=. For Oo = ~ the approximate solution is exact and the error will increase from o= - to Bo = 71

THE UNIVERSITY OF MICHIGAN 1.4 1.2 0.2.4...8 1.0 1.4 1.6 i.6 rl / 61 -42, D- i~2~ ~ ~I~ ~ - 2 I I } JI O.2.4.,.8 1.0 o 1i(kr I -ka cos ( o. 1.2 For ka = 1 thru 7 00=1T 0 ~ ~ ~ ~^ ~ ~ ~^ ~ ~ For = L,,.1.2 2 \ 4 -i.4 -i.6 72~Exac o.2.4.6 8 1.0 1.2 I. l.6 FIG. A.7 E FROM EXACT SOLUTION (9) AND ASYMPTOTIC SOLUTION (17) i.~ ~~~72

THE UNIVERSITY OF MICH I GAN i,r, i.4.I~ ". E i(kr _-ka cos (~4 ) ^~ 0 5 1.0 (-i)Voe 4 /k i-4 For ka = 1 thru 7 For 57 3,WI 0-0_ 6 4 and ~2 e/ -., Approximate -i.4 _I -— * — Exact OO/ 0.5 1.0 1.5 i.4 21L 0 0.~ 1.O -~.5 FIG. A.8 Ed FROM EXACT SOLUTION (9) AND ASYMPTOTIC SOLUTION OBTAINED BY CONTINUED FRACTIONS 73 73

THE UNIVERSITY OF MICHIGAN 2h72-13-T APPENDIX B. COME We determine the field due to a magnetic dipole source in the presence of a perfectly conducting semi-infinite cone. With a time harmonic dependence of the form e-it Maxwell's equations are Vx E = ik H + M, k = _ (B.I) c VxH =-ikE, where we assume a point source M m 6 (r'-r) Since the cone is a surface separable in spherical coordinates, we can make use of the modal representations Emode: E =Tx IT (B.2) H -ik LT H mode: E = ik L (B.3) H = vx LB where L =-i r x7is the angular momentum operator. The scalar Hertz potentials satisfy the equation (V2 +k2)To(T) 0 (B.4) away from the source. 74

THE UNI VE RS ITY OF M I CH I GAN 2472-13-T If we eliminate the electric field from Maxwell's equations, we find that the E-mode Hertz functionIT satisfies (V + k2) LT= -M. (B.5) Since the operator L commutes with the Laplacian, we form (v2 + k2) L2T= -L* (B.6) or (V2 + k2) _ 1 L'*. M (r-.') (B.7) L2t' where the primes indicate an operation with respect to the variable (xt, yt, z'). Inverting the differential operator of Equation B.6 T= i L'* I m GD (r, r' ) (B.8),2 L' where GD is the Dirichlet Green's function for the cone. Similarly, eliminating the magnetic field the H-mode Hertz function satisfies (v2 + k2) L= - Vx M. (B.9) So by the above reasoning w~e find ( + k2)" - 1 _ L* xV &( ) (B.10) ik L7 75

THE UNIVE RS ITY OF MICHIGAN 2472-13-T and TT 1 1 Li x7 ~ m GN (r,r' ) (B.11) ik L where GN (rri') is the Neumann Green's function for the cone. From the above we see that the decomposition of the exact solution into essentially a geometric optics part plus a part asymptotic in the wave number can result only from such a decomposition of the scalar Green's functions. That this is indeed the case has been shown by Felsen (see Ref. 2). Consider a ring source of magnetic dipoles centered about the axis of a perfectly conducting semi-infinite cone of angle 90>. Performing the operation indicated in (B.8) and(B.2) above, we have for a uniformly excited source that H(rr ) = (-i) dv (2v + (B12) H im> (B.12) 2T G (O,+ 1)?'.J (kr ) h ) (kr>).. G (,9V) where the contour includes the positive real axis in the complex plane and (cos 9 )P (-cos o>) sinvTr 1 P (-cos Qo) n+ FT (cos 0l PV (cos 9 ) Pv (cos o) ~ sin-7 P7 (cos 6 o) 76

THE UNIVERSITY OF MI CHIGAN 2472-13-T is the decomposition of the scalar Green's function mentioned above. Provided 9 + O9 < 2G - n, the contour can be deformed in the case of the second term to run from-+ i oo to - - i oo. This gives the diffracted term:' m -ikr x H dX H diff 8 X2 + (B13) H(1) (kr) H(1) (kr') K (-cos 90) ix ix _ x. K' (cos e)K' (cos o') cosh T x K (cos 9) x m where Kx (cos ) = P 1 + ix (cos 9). For sufficiently large kr and kr' the asymptotic form of H(l)may be employed in Equation B.13 to give ix Hdiff m eik(r+r) dx x e (B.14 rrr 2 (x2 + 1)cosh iTX -oo Kx (cos 9)Kx (cos 91). K (X os ). Kx (cos Go) For values of 9onear T(small cones), the employment of the asymptotic form Kx(-cos 0) n -1, Q_ (B.15) Kx(cos Go) 2 (I-0o0) cosh Tx 77

THE UNIVERSITY OF MICH I GAN 2472-13-T in Equation B.14 permits the integration of (B.14) and yields, to the first order, for a small cone approximation H ~ ^1eik(r+r') 21 tan ta (B.16) Hdiff rr - )O, 8 \ cos 2 +cos 9 If the source is permitted to go to the surface of the cone, the same form for Hdiff of Equation B.16 results; hence, an evaluation is obtained for the higher order terms in the case of a circumferential slot on a small cone. If the ring source is on the surface ot the cone, the decomposition, as was pointed out by Felsen, depends upon breaking up the Legendre functions into positive and negative exponentials. That is, let (1) (2) P (-cos 9) - e +V Pp (cos 9) Dv (G)~ ~___________(B.I7) 2 i sin e so that (2) expvli- i- L.) - o2 r v sin ie for sin >>~l. Substitution of Equation B317 in Equation B.12 yields integrals of the form 78

THE UNIVE RS ITY OF MI CHIGAN 2472-13-T 1 ) iT P (cos 9) I dv j (krt) ei(V+ 7 )P (cos ) C1 ~~ v (I(+ 1) (1) (e) r _+ j_ (, e -i((,))i P2 (cos ) D(2) (e0) (v+ 1) D() (9 ) (9 (cos 8) +f vd - j (kr) e-i (v+ (cos V(+ 1) D(1) (eo) 2 d 7v+ -i (v + P (cos 9) D(2) (60) - I ~dv j (kr) e J v( + ) D(1) (90) PV (cos 9o) C2 where the contours are indicated below: v -plane Cl Cl C2 2 I The integrals over C1' and C2 can then be evaluated, using the asymptotic form of j. for large kr'. These lead to the diffracted term, while a saddle point evaluation of the integral over Cl leads to the geometric optics term. 79 1

THE UNIVERSITY OF MICHIGAN 2472-13-T APPENDIX C. GRAPHICAL PRESENITATIONS The expression (3.2) of Chapter 3 has been employed to obtain radiation patterns in the region 0 e6< -e80 for the following values of the parameters ka, g0, and n: ka = 50w, 75i; 9 -= 1600, 165~, 1700; n 1, 2, 3. These patterns are shown in Figures C-1 through C-18. The expression actually calculated is sin ~ ka 2 cos2 where B is given by Equation 3.5 and St and T' are given by Equation 3.3. 80

THE UNIVE R S ITY OF M I CH I GAN 2472-13-T 2.8 0 2.4 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1DEGREES 1.6 ~ ~ IBEI 01 DEGREES ~ 281

THE UNIVE RS ITY OF MICHIGAN 2472-13-T 2.8 2.4 Oo = 160" ka = 50 - n 2.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 IBE.4 0 2 4 6 8 10 12 14 16 6 - DEGREES FIG. C.2 40~ CONE WITH e2i EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 82

THE UNI VE RS ITY OF M I CHIGAN 2472-13-T 2.8 2.4 00= 160' ka = 507r n = 3.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 E - lEI 1.2.8.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG C.3 40~ CONE WITH e3i0 EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 83

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 00= 165 ka = 507r n= 1.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BEI 1.2.8.4 0 2 4 6 8 10 12 14 16 0 - DEGREES FIG. C.4 300 CONE WITH eiV EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 84

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 Oo 1650 ka = 507r n = 2.0 2.0....... 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 | BE| 1.2...8.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.5 30~ CONE WITH e2iV EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 85

THE UNIVERSITY OF MICH I GAN 2472-13-T 2.8 2.4 0 = 165" ka = 507r n =3.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 IBEI 1.2 0 2 4 6 8 10 12 14 16 - DEGREES FIG C.6 30~ CONE WITH e3i EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP ~ ~~ ~ ~~~~~~~86 ~

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 800 1700 ka= 507r n= 1.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 IBE 1.2.8 0 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.7 200 CONE WITH e EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 87

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 _~~___ 2.4 60 = 1700 ka 50rn = 2.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BE 1.2.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.8 20~ CONE WITH e2 EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 88

THE UNIVE RSITY OF MICHIGAN 2472-13-T 2.8. 2.4 60 = 170' ka= 50Or n = 3.0 0 2 4 6 8 10 12 14 16 2.0 DEGREES 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BE 1.2 FIG. C.9 200 CONE WITH e3i EXCITED CIRCUMFERENTIAL SLOT 25 WAVELENGTHS FROM TIP 89

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 _ 2.4 00 = 160~ ka= 75 7r n= 1.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 IBE 1.2.8.4 0 2 4 6 8 10 12 14 16 0 - DEGREES FIG C.10 400 CONE WITH ei EXCITED CIRCUNFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP 90

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 0o = 160' ka = 75 r n= 2.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6BE.8.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.ll 400 CONE WITH e2i~ EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELNGTHS FROM TIP 91

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 - 0 = 1600 ka= 75 r n = 3.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 IBE 1.2.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.12 40~ CONE WITH e3i EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP 92

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 oo 1650 ka = 75 r n= 1.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BE 1.2.8.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.13 300 CONE WITH ei- EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELEGTHS FROM TIP 93

THE UNI VERSITY OF M I CH I GAN 2472-13-T 2.81.8 2.0 2.4 0= 165 ka = 75 7r "=2.0~ 2.0 0.2.4.6,.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BE 1.2.4 0 2 4 6 8 10 12 14 16 0 - DEGREES FIG. C.14 300 CONE WITH e2i EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP __________________________ 9~-4

THE UNIVERSITY OF MICH I GAN,2472-13-T 2.8 2.4 00 = 1650 ka = 75 7r n = 3.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 | B E| I~1 1.2 0 4~ ~ I~~ 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.15 30~ CONE WITH e3iV EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP 95

THE UNIVERSITY OF MICHIGAN 2k72-13-T 2.8 2.4 00 = 170' ka = 75 r n = 1.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 181 1.2.8.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.16 200 CONE WITH eiV EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP 96

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 00 = 170~ ka = 75 nr n = 2.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BE 1.2.8 __.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.17 200 CONE WITH e2i~ EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP 97

THE UNIVERSITY OF MICHIGAN 2472-13-T 2.8 2.4 60 = 1700 ka = 75 r 2 = 3.0 2.0 0.2.4.6.8 1.0 1.2 1.4 1.6 1.8 2.0 1.6 BEI 1.2.8 ~.4 0 2 4 6 8 10 12 14 16 - DEGREES FIG. C.18 20~ CONE WITH e EXCITED CIRCUMFERENTIAL SLOT 37.5 WAVELENGTHS FROM TIP 98

THE UNIVERSITY OF MICHIGAN 2472-13-T REFEREN ICES 1. C. E. Schensted, "Electromagnetic and Acoustic Scattering by a Semi-infinite Body of Revolution", Jour. Appl. Phys., 26, No. 3, 306-309 (1955). 2. L. B. Fclsen, Private communications to the authors, 3. V. A. Fock, Journal of Physics X, 399 (1946). 4.L. L. ailin, "Field Produced by a Slot in a Large Circular Cylinder", IRE Trans. Antennas and Propagation, AP-3, No. 3, 128-138 (1955). 5. W. Franz and K. Depperman, Ann. Phys. 10, 361 (1952). I. Imai, Z. Phys. 137, 31 (IT9) W. Franz, Z. fur Naturforshung 9a, 705 (1954). 6. L. B. Felsen, "Back Scattering from Wide-Angle and Narrow-Angle Cones", J. A. P. 26, 138-151 (1955). 7. L. L. Bailin and S. Silver, "Exterior Electromagnetic Boundary Value Problems for Spheres and Cones"t, IRE Trans. Antennas and Propagation, AP-h, No. 1, 5-16 (1956). 8. Daniel Ray, "On Spectra of Second-Order Differential Operators", NBS Report No. 2075. 9. R. H. Cameron and v. T. Martin, "An Expression for the Solution of a Class of Non-Linear Integral Equations", American J. Math., 66, 281-298 (1944). 10. R. H. Cameron, "A'Simpson's Rule' for the Numerical Evaluation of Wiener Wiener's Integrals in -unction Space." Duke Math. Journ. 18, 111-130 (1951). 11. L. B. Felsen, "Radiation from Slots on a Perfectly Conducting Wedge", R-L92-56 Memorandum No. 13, PIB-422, Polytechnic Institute of Brooklyn. 12. F. Oberhettinger, "Diffraction of Waves by a Wedge." Communications on Pure and Applied Mathematics 7, 551-563, (1954). 99

THE UNIVERSITY OF MICHIGAN 2472-13-T 13. W. Hansen and L. I. Schiff, "Theoretical Study of Electronagnetic Waves Scattered from Shaped Metal Surfaces." Microwave Laboratory, Stanford University, Quarterly Report (Feb. 1948). 14. J. Keller and A. Blank, "Diffraction and Reflection of Pulses by Wedges and Corners, N.Y.U. Symposium on the Theory of Electromagnetic Waves (1950). 15. J. Franklin and B. Friedman, "A Convergent Asymptotic Representation for Integrals", New York University, Institute of Mathematical Sciences, Division of Electromagnetic Research, Research Report No. BR-9. 16. F. Oberhettinger, "On Asymptotic Series for Functions Occurring in the Theory of Diffraction of Waves by Wedges", Journal of Mathematics and Physics 34, No. l, January 1956. See Eq. 32. 17. H. S.'Wall, "Continued Fractions", D.Van Nostrand Company Inc., (1948) 18. N. Logan, Private communication to the authors. ______________________________ 100