Radiation Integrals and Auxiliary Potential Functions Ranga Rodrigo June 23, 2010 Lecture notes are fully based on Balanis [?]. Some diagrams and text are directly from the books. Contents 1 The Vector Potential for an Electric Current Source J 4 2 Solution of the Inhomogeneous Vector Potential Equation 6 3 Far-Field Radiation 9 Introduction In the analysis of radiation problems, the usual procedure is to specify the sources and then require the fields radiated by the sources. This is in contrast to the synthesis problem where the radiated fields are specified, and we are required to determine the sources. It is a very common practice in the analysis procedure to introduce auxiliary functions, known as vector potentials, which will aid in the solution of the problems. The most common vector potential functions are the A (magnetic vector potential and F (electric vector potential. 1
Sources J, M Integration path 1 Radiated fields E, H Integration path 2 Differentiation path 2 Vector potentials A, F The Divergence The basic definition of divergence is [ A = lim v 0 S A ds v The expansion of divergence in Cartesian coordinates is ]. A = A x x + A y y + A z z. The Curl The basic definition of curl is A = lim S 0 [ C A dl ] i n. S max The expansion of curl in Cartesian coordinates is i x i y i z A = x y z. A x A y A z Vector Identities div curl = 0 A = 0 curl grad = 0 φ = 0 2
Maxwell s Equations in Integral Form E dl = d B ds. C dt S H dl = J ds + d C S dt D ds = ρdv. S V B ds = 0. S S D d S. (Faraday s law (Ampère s circuital law (Gauss law for the electric field (Gauss law for the magnetic field Auxiliary equation J ds = d dt S V ρdv. (Law of conservation of charges Other relations D = εe, H = B µ, Maxwell s Equations in Differential Form E = B t. H = J + D t. D = ρ. B = 0. Auxiliary equation J = ρ t. Maxwell s Equations in Free Space 3
E = B t. H = D t. D = 0. B = 0. 1 The Vector Potential for an Electric Current Source J The vector potential A is useful in solving for the EM field generated by a given harmonic electric current J. The magnetic flux B is always solenoidal; that is, B = 0. Therefore, it can be represented as the curl of another vector because it obeys the vector identity A = 0. where A is an arbitrary vector. Thus we define B A = µh A = A. (1 We can substitute this into Maxwell s curl equation E A = j ωµh A, This can also be written as E A = j ω A. [ E A + j ω A ] = 0. From the vector identity ( φ e = 0, it follows that E A + j ωa = φ e, or E A = φ e j ωa. 4
The scalar function φ e represents an arbitrary electric scalar potential which is a function of position. Taking the curl of both sides of Eq. 1 and using the vector identity A = ( A 2 A, reduces it to For homogeneous medium, (µh A = ( A 2 A µ H A = ( A 2 A Using Maxwell s equation, H A = J + j ωɛe A, µj + j ωµɛe A = ( A 2 A Substituting 11, 2 A + k 2 A = µj + ( A + j ωµɛφ e where k 2 = ω 2 µɛ. In 1 the curl of A was defined. Now we are at liberty to define the divergence of A, which is independent of its curl. Let A = j ωɛµφ e, which is known as the Lorentz condition. Substituting, In addition, 11 leads to 2 A + k 2 A = µj. 1 E A = φ e j ωa = j ωa j ( A. (2 ωµɛ Once A is known, H A can be found from Eq. 1 and E A from Eq. 2. E A can just as easily be found from Maxwell s equation, H A = J + j ωɛe A, with J = 0. It will be shown later how to find A in terms of the current density J. It will be a solution to the inhomogeneous Helmholtz equation of Eq. 11. Although magnetic currents appear to be physically unrealizable, equivalent magnetic currents arise when we use the volume or the surface equivalence theorems. The fields generated by a harmonic magnetic current in a homogeneous region, with J = 0 but M 0, must satisfy D = 0. 5
Therefore, E F can be expressed as the curl of the vector potential F by E F = 1 ɛ F. (3 By letting we can obtain F = j ωµɛφ m (4 2 F + k 2 F = ɛm. Once F is known, E F can be found from Eq. 3 and H F with M = 0. It will be shown later how to find F once M is known. It will be a solution to the inhomogeneous Helmholtz equation of 4. Total Fields Now, we have developed equations that can be used to find the electric and magnetic fields generated by an electric current source J and a magnetic current source M. The procedure requires that the auxiliary potential functions A and F generated, respectively, by J and M are found first. In turn, the corresponding electric and magnetic fields are then determined (E A,H A due to A and E F,H F due to F. The total fields are then obtained by the superposition of the individual fields due to A and F (J and M. E = E A + E F = j ωa j 1 ωµɛ ( A 1 ɛ F = 1 j ωɛ H A 1 ɛ F. H = H A + H F = 1 µ A j ωf j 1 ωµɛ ( F = 1 µ A 1 j ωµ E F. 2 Solution of the Inhomogeneous Vector Potential Equation The solution of the inhomogeneous vector wave equation is A = µ kr J 4π R dv. V To derive it, let us assume that a source with current density J z, which in the limit is an infinitesimal source, is placed at the origin of a x, y, z coordinate system. 6
z (x, y, z θ (x, y, z φ y x Since the current density is directed along the z-axis (J z, only an A z component will exist. Thus we can write 2 A + k 2 A = µj as 2 A z + k 2 A z = µj z. At points outside of the infinitesimal source, 2 A z + k 2 A z = 0. (5 Point source: A z is not a function of direction (θ and φ and A z = A z (r, where r is the radial distance. 2 A z (r + k 2 A z (r = 1 [ r 2 r 2 A ] z(r + k 2 A z (r = 0. r r d 2 A z (r dr 2 + 2 d A z (r + k 2 A z (r = 0. r dr This differential equation has two independent solutions j kr e A z1 = C 1 r +j kr e A z2 = C 2 r : radially outward traveling wave. : radially inward traveling wave. 7
Therefore, we choose j kr e A z = A z1 = C 1. r In the static case (ω = 0, k = 0, this simplifies to A z = C 1 r. which is the solution to the wave equation 6 when k = 0. Thus at points outside the source, the time-varying and static solutions differ only by e j kr, the phase retardation factor. In the presence of the source (J z 0, and k = 0, the wave equation reduces to 2 A z = µj z. (6 This equation is recognized to be Poisson s equation, the most familiar form being that relates the scalar electric potential φ to the electric charge density ρ. The solution is 2 φ = ρ ɛ. φ = 1 4πɛ V ρ r dv, where r is the distance from any point in the charge density to the observation point. Using a similar form A z = µ 4π V J z r dv. This represents the static solution. The time-varying solution is A z = µ j kr e J z dv. 4π V r If we consider the current densities in x- and y-directions (J x, and J y, A = µ kr J dv. 4π V r If the source is removed from the origin and placed at a position represented by (x, y, z, A(x, y, z = µ J(x, y, z kr 4π V R dv. 8
where (x, y, z is the observation point and R is the distance from any point on the source to the observation point. In a similar fashion we can obtain If J and M represent linear densities, F (x, y, z = ɛ M(x, y, z kr 4π V R dv. A(x, y, z = µ J s (x, y, z kr 4π S R ds. F (x, y, z = ɛ M s (x, y, z kr 4π S R ds. For electric and magnetic current I e and I m A(x, y, z = µ I e (x, y, z kr 4π C R dl. F (x, y, z = ɛ I m (x, y, z kr 4π C R dl. 3 Far-Field Radiation The far-fields radiated by antennas of finite dimensions are spherical waves. For these radiators, a general solution to the vector wave equation is spherical components is A = â r A r (r,θ,φ + â θ A θ (r,θ,φ + â φ A φ (r,θ,φ. The amplitude variations of r in each component are of the form 1/r n. Neglecting the higher order terms, for large r, Consider ] e A [â r A r (r,θ,φ + â θ A θ (r,θ,φ + â φa j kr φ (r,θ,φ. r 1 E A = φ e j ωa = j ωa j ( A. ωµɛ If we evaluate the second term in spherical coordinates, we will notice that higher powers of 1/r are produced. Neglecting them, E A j ωa for θ and φ components. 9
Neglecting the higher order terms of 1/r n, the radiated E- and H-fields have only θ and φ components. They can be expressed as H r 0, E r 0, E θ j ωa θ, E φ j ωa φ. H θ +j ω η A φ = E θ η, H φ j ω η A θ = + E φ η. Radial far-field terms exist only for higher order terms of 1/r n. In a similar manner, the far-zone fields due to magnetic source M (potential F can be written as E r 0, H r 0, H θ j ωf θ, H φ j ωf φ. E θ j ωηf φ = ηh φ, E φ +j ωηf θ = ηh θ The corresponding far-zone E- and H-field components are orthogonal to each other and form TEM (to r mode fields. 10