Waveguides Free Space. Modal Excitation. Daniel S. Weile. Department of Electrical and Computer Engineering University of Delaware

Similar documents
Summary: Method of Separation of Variables

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)

Patch Antennas. Chapter Resonant Cavity Analysis

Waveguide Guide: A and V. Ross L. Spencer

Consequently, the temperature must be the same at each point in the cross section at x. Let:

The Dirac distribution

Cartesian Coordinates

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

3 Mathematics of the Poisson Equation

1 E3102: a study guide and review, Version 1.0

Phys 6321 Final Exam - Solutions May 3, 2013

Chapter 5 Waveguides and Resonators

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Orthogonal functions

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

This final is a three hour open book, open notes exam. Do all four problems.

c n φ n (x), 0 < x < L, (1) n=1

1 1D heat and wave equations on a finite interval

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6

Candidates must show on each answer book the type of calculator used.

( ) 2. ( ) is the Fourier transform of! ( x). ( ) ( ) ( ) = Ae i kx"#t ( ) = 1 2" ( )"( x,t) PC 3101 Quantum Mechanics Section 1

Phys. 506 Electricity and Magnetism Winter 2004 Prof. G. Raithel Problem Set 1 Total 30 Points. 1. Jackson Points

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Aike ikx Bike ikx. = 2k. solving for. A = k iκ

Conservation Law. Chapter Goal. 5.2 Theory

Lecture 13 - Linking E, ϕ, and ρ

Green function and Eigenfunctions

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

Partial Differential Equations

The Dirichlet Problem in a Two Dimensional Rectangle. Section 13.5

Lecture Outline. Dispersion Relation Electromagnetic Wave Polarization 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3c

DEFINITION The inner product of two functions f 1 and f 2 on an interval [a, b] is the number. ( f 1, f 2 ) b DEFINITION 11.1.

Physics 3323, Fall 2016 Problem Set 7 due Oct 14, 2016

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015

Problem Set 3 Solutions

New Expansion and Infinite Series

Pressure Wave Analysis of a Cylindrical Drum

SYDE 112, LECTURES 3 & 4: The Fundamental Theorem of Calculus

Jackson 2.7 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

Introduction to Group Theory

Theoretical foundations of Gaussian quadrature

g i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f

Bernoulli Numbers Jeff Morton

Wave Phenomena Physics 15c. Lecture 12 Multi-Dimensional Waves

ODE: Existence and Uniqueness of a Solution

Convergence of Fourier Series and Fejer s Theorem. Lee Ricketson

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Quantum Physics III (8.06) Spring 2005 Solution Set 5

PDE Notes. Paul Carnig. January ODE s vs PDE s 1

Point Lattices: Bravais Lattices

1 E3102: A study guide and review, Version 1.2

ENGI 9420 Lecture Notes 7 - Fourier Series Page 7.01

Riemann Sums and Riemann Integrals

Physics 215 Quantum Mechanics 1 Assignment 2

Read section 3.3, 3.4 Announcements:

Problems for HW X. C. Gwinn. November 30, 2009

The Velocity Factor of an Insulated Two-Wire Transmission Line

Homework Assignment 5 Solution Set

Main topics for the First Midterm

Riemann Sums and Riemann Integrals

Continuous Quantum Systems

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

The solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr

Best Approximation in the 2-norm

Section 17.2 Line Integrals

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Homework Assignment #1 Solutions

Chapter E - Problems

Theoretische Physik 2: Elektrodynamik (Prof. A.-S. Smith) Home assignment 4

Numerical integration

221B Lecture Notes WKB Method

13: Diffusion in 2 Energy Groups

Improper Integrals, and Differential Equations

Chapter 28. Fourier Series An Eigenvalue Problem.

Lesson 1: Quadratic Equations

Math 426: Probability Final Exam Practice

x = b a n x 2 e x dx. cdx = c(b a), where c is any constant. a b

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Green s functions. f(t) =

Green s function. Green s function. Green s function. Green s function. Green s function. Green s functions. Classical case (recall)

Review SOLUTIONS: Exam 2

Quantum Physics I (8.04) Spring 2016 Assignment 8

) 4n+2 sin[(4n + 2)φ] n=0. a n ρ n sin(nφ + α n ) + b n ρ n sin(nφ + β n ) n=1. n=1. [A k ρ k cos(kφ) + B k ρ k sin(kφ)] (1) 2 + k=1

Cylindrical Radiation and Scattering

10 Vector Integral Calculus

Lecture 1. Functional series. Pointwise and uniform convergence.

Review of Calculus, cont d

Section 6.1 INTRO to LAPLACE TRANSFORMS

Chapter 0. What is the Lebesgue integral about?

Math 8 Winter 2015 Applications of Integration

Homework Assignment 9 Solution Set

Lecture 24: Laplace s Equation

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term Solutions to Problem Set #1

Phys 4321 Final Exam December 14, 2009

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

7.3 Problem 7.3. ~B(~x) = ~ k ~ E(~x)=! but we also have a reected wave. ~E(~x) = ~ E 2 e i~ k 2 ~x i!t. ~B R (~x) = ~ k R ~ E R (~x)=!

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

Physics 712 Electricity and Magnetism Solutions to Final Exam, Spring 2016

Transcription:

Modl Excittion Dniel S. Weile Deprtment of Electricl nd Computer Engineering University of Delwre ELEG 648 Modl Excittion in Crtesin Coordintes

Outline 1 Aperture Excittion Current Excittion

Outline 1 Aperture Excittion Current Excittion 2 Apertures in Ground Plne Plne Current Sheets

Outline Aperture Excittion Current Excittion 1 Aperture Excittion Current Excittion 2 Apertures in Ground Plne Plne Current Sheets

Excittion of Aperture Excittion Current Excittion Our discussion so fr hs ignored how energy might be coupled into guides. After ll mode is by definition solution tht cn exist without source. How re modes excited? Let us consider our stndrd guide. Over the z = plne, let us ssume E x = E y = f (x, y) We further ssume the wveguide is mtched somewhere for z > nd seek the solution for the modes lunched.

Modl Expnsion Aperture Excittion Current Excittion Since E x = we look for TE x solutions. We cn thus write Modl Expnsion for ψ ψ = m=1 n= A mn sin mπx To find the A mn we must impose f (x, y) = E y z= = m=1 n= How do we find the coefficients? nπy cos b e γmnz γ mn A mn sin mπx nπy cos b

Aside: Useful Nottion Aperture Excittion Current Excittion Two functions tht often rise in mthemticl physics re the Kronecker Delt δ mn = { 1 if m = n otherwise nd the ever present nd strngely nmed Neumnn s Number ɛ n = { 1 if n = 2 if n

Aside: Orthogonl Functions Aperture Excittion Current Excittion Note the Orthogonlity Reltions sin mπx cos mπx sin mπx Know these! nπx sin dx = 2 δ mn for m, n = 1, 2,... nπx cos dx = δ mn for m, n =, 1, 2,... ɛ m nπx cos dx = for m, n =, 1, 2,...

Finding the Coefficients Aperture Excittion Current Excittion To find the coefficients, we multiply nd integrte: = b Thus m=1 n= f (x, y) sin m πx γ mn A mn = b m=1 n= A mn = 2ɛ n bγ mn b cos n πy b dxdy sin mπx γ mn A mn b 2ɛ n nπy cos b sin m πx f (x, y) sin mπx cos n πy b dxdy δ mm δ nn = γ m n A b m n 2ɛ n nπy cos b dxdy

Observtions Aperture Excittion Current Excittion A similr procedure cn be followed if the mgnetic field is given. Excittion by field with both x- nd y-components cn be solved by superposition. Note tht A finite number of modes propgtes. An infinite number of modes evnesces. How cn we excite given mode? Cn we void exciting some modes entirely? How?

Power Flow Aperture Excittion Current Excittion Becuse our modes re TE x, we cn compute the power flow in the guide using Plugging in P = dx b dy P = [ m=1 n= b p=1 q= [E y H x ] z= dxdy. E mn sin mπx ] nπy cos b k 2 ( pπ ) 2 jωµγpq Epq sin pπx cos qπy b Why do we use different vribles in the two summtions?

Power Flow Aperture Excittion Current Excittion Thus P = m=1 n= p=1 q= = k 2 ( pπ ) 2 jωµγpq b m=1 n= dy cos nπy b E pqe mn qπy cos b b k 2 ( mπ ) 2 2ɛ n jωµγmn = where (Y ) mn is the modl dmittnce. Thus, ech mode crries its own power. E mn 2 m=1 n= dx sin mπx pπx sin b 2ɛ n (Y ) mn E mn 2

Wveguide Junction Aperture Excittion Current Excittion Consider the junction between the wveguides below. Their width is. y c b z We ssume the field in the junction is { sin πx E y z= y < c y > c

Wveguide Junction Aperture Excittion Current Excittion Now, we hve seen tht the field in the junction is given by f (x, y) = E y z= = m=1 n= γ mn A mn sin mπx nd tht the coefficients cn be computed from γ mn A mn = E mn = 2ɛ n b b f (x, y) sin mπx Immeditely, we see the only excited modes hve by orthogonlity. m = 1 nπy cos b nπy cos b dxdy

Computing the Coefficients Aperture Excittion Current Excittion Plugging in we find nd E 1 = 2 c sin 2 πx b dxdy = c b, E 1n = 4 c sin 2 πx nπy cos b b dxdy = 2 nπc sin nπ b. Most of these modes re of course evnescent.

Power Flow nd Admittnce Aperture Excittion Current Excittion Applying our power formul P = c2 [ ( 2b (Y ) 1 + 2 sin nπ )] 2 (Y ) b 1n nπc n=1 We cn lso define n perture dmittnce by defining the perture voltge s the line integrl of E cross the center of the perture. Then Y = 2P V 2 b

Outline Aperture Excittion Current Excittion 1 Aperture Excittion Current Excittion 2 Apertures in Ground Plne Plne Current Sheets

Aperture Excittion Current Excittion Current Excittion of In perture excittion, we mtch mode series to given field in the junction. In current excittion, we re given current. The electric field is continuous t the current. The mgnetic field is discontinuous t the current. Note tht in this wy, perture excittion cn be thought of s mgnetic current excittion. Suppose for simplicity we ssume current sheet in the z = plne with J s = u x f (x, y) nd tht the guide is mtched in both directions. We cn ssume the wves generted re TM x. (Why?)

Aperture Excittion Current Excittion Current Excittion of The most generl solution for TM wves cn be written s ψ + = ψ = m= n=1 m= n=1 B + mn cos mπx B mn cos mπx nπy sin b e γmnz z > nπy sin b e+γmnz z < Continuity of the tngentil E-field (which is proportionl to ψ nd its trnsverse derivtives) requires B mn = B + mn = B mn

Aperture Excittion Current Excittion Current Excittion of The mgnetic field hs to be discontinuous by the current, i.e. Thus, u z ( H + H ) = J s f (x, y) = [ H y H y + ]z= = m= n=1 2γ mn B mn cos mπx nπy sin b nd we hve the Field Coefficients 2γ mn B mn = J mn = 2ɛ m b dx b dy f (x, y) cos mπx nπy sin b

Power Flow Aperture Excittion Current Excittion The power supplied cn be computed from wht we know of the Poynting Theorem: P = dx b dy E x Jx z= After substituting the Fourier series for E x nd J x nd ppeling to orthogonlity, we find n expression for Power Flow P = m= n=1 b 4ɛ m (Z ) mn J mn 2

An Exmple Aperture Excittion Current Excittion Consider coxilly fed wveguide with the probe ffixed to the opposite wll: y c b x The probe cts s n ntenn; by imge theory, we expect the current to be pproximtely Probe Current J x = I cos[k( x)]δ(y c)

An Exmple Aperture Excittion Current Excittion J mn = 2ɛ m b dx b = 2I ɛ m b dy f (x, y) cos mπx nπc sin b nπy sin b dx cos[k( x)] cos mπx = 2I ɛ m k sin k sin nπc b b[(k) 2 (mπ) 2 ] The impednce seen by the probe cn be computed in the usul wy: where I i is the feed current Z i = 2P I i 2 I i = I cos k

Outline Apertures in Ground Plne Plne Current Sheets 1 Aperture Excittion Current Excittion 2 Apertures in Ground Plne Plne Current Sheets

Hve We Seen This Before? Apertures in Ground Plne Plne Current Sheets We lredy know one wy to compute the field rdited from n perture, but consider: All methods must give the sme result becuse of uniqueness. This method is spectrl method. Spectrl methods converge fstest when sptil methods converge slowest. Why?

Hve We Seen This Before? Apertures in Ground Plne Plne Current Sheets We lredy know one wy to compute the field rdited from n perture, but consider: All methods must give the sme result becuse of uniqueness. This method is spectrl method. Spectrl methods converge fstest when sptil methods converge slowest. Why? This method elucidtes different physics. Thus, suppose we re given the perture field E x (x,, z) nd we wnt the field for y >. How should we proceed?

Modl Expnsion Apertures in Ground Plne Plne Current Sheets As before, we cn expnd the field in modes. Now, however, we hve continuum of modes, rther thn discrete set. Therefore, we choose the Wve Function ψ(x, y, z) = 1 4π 2 dk x dk z f (k x, k z )e jkx x e jky y e jkzz Note tht this is subject s usul to the usul seprtion or dispersion reltion k 2 x + k 2 y + k 2 z = k 2 Our unknown function is the Fourier Trnsform of ψ: ψ(k x, k z ) = f (k x, k z )e jky y

Modl Expnsion Apertures in Ground Plne Plne Current Sheets Given E x in constnt y plne, we choose TE z modes. From before, these modes re relted to the wve function with E x = ψ y H x = 1 E y = ψ x H y = 1 E z = H z = 1 jωµ We cn esily compute the 2 ψ jωµ x z 2 ψ jωµ y z ( 2 z 2 + k 2 ) ψ Spectrl Domin TE z Equtions Ẽ x = jk y ψ Ẽ y = jk x ψ Ẽ z = Hx = Hy = kx kz jωµ ky kz jωµ Hz = k 2 k 2 z jωµ ψ ψ ψ

The Solution Apertures in Ground Plne Plne Current Sheets Since we know the electric field in the perture, we cn write Ẽ x (k x, k z ) = jk y f (k x, k z ) y= or f (k x, k z ) = 1 Ẽ x (k x, k z ) jk y y= Of course, cre must be tken in the computtion of k y : { j kx k y = 2 + kz 2 k 2 for k < kx 2 + kz 2 k 2 kx 2 kz 2 for k > kx 2 + kz 2 How re the signs chosen? Wht does this sy bout determining the current fr wy from it?

An Exmple: The Infinite Slot Apertures in Ground Plne Plne Current Sheets To simplify the exmple, we use 2D problem independent of z. Consider the infinite slot shown below: x y If we consider the TEM mode, H is z directed everywhere. We therefore use H z s our wve function: H z = 1 2π f (k x )e jkx x e jky y dk x

The Infinite Slot Apertures in Ground Plne Plne Current Sheets From this expression, we see immeditely H z = f (k x )e jky y To get the electric field, we use Ampere s Lw: H = jωɛe Now, in the spectrl domin, this becomes k H = jωɛẽ From symmetry, this eqution becomes k x H z = jωɛe y k y H z = jωɛe x

The Infinite Slot Apertures in Ground Plne Plne Current Sheets Now, ssume we know the electric field in the slot. We hve just shown Ẽ x y= = k y ωɛ f (k x) Now, given the incident wve nd the symmetry in the z-direction, we cn pproximte Therefore Ẽ x y= = E E x y= { E x < 2 x > 2 2 2 e jkx x dx = 2E k x sin k x 2

The Infinite Slot Apertures in Ground Plne Plne Current Sheets From this we find tht f (k x ) = ωɛ 2E sin k x k y k x 2 The remining field cn be found from this vi the formuls on the previous pges. As usul, the vlue of k y is dictted by ensure tht the wve Evnesces (nd does not explode) in the direction of propgtion, or, Trvels wy from (nd not towrds) the perture. This implies k y = { j k 2 x k 2 for k < k x k 2 k 2 x for k > k x

Outline Apertures in Ground Plne Plne Current Sheets 1 Aperture Excittion Current Excittion 2 Apertures in Ground Plne Plne Current Sheets

Plne Current Sheets Apertures in Ground Plne Plne Current Sheets Plne current sheets excittion of free spce cn be computed in the sme wy. Assuming the sheet is in the y = plne, nd A = µψu z, we cn write the Fourier trnsformed (with respect to x nd z) equtions for TM z propgtion: H x = jk y ψ H y = jk x ψ H z = E x = kx kz jωɛ ky kz Ey = jωɛ E z = k 2 k 2 z jωɛ ψ ψ ψ

Plne Current Sheets Apertures in Ground Plne Plne Current Sheets We need to expnd ψ on either side of the sheet: ψ + = f + (k x, k z )e jk + y y ψ = f (k x, k z )e jk y y for y > for y < Of course, k + y = k y = { j kx 2 + kz 2 k 2 k 2 kx 2 kz 2 for k < kx 2 + kz 2 for k > kx 2 + kz 2

Boundry Condition Apertures in Ground Plne Plne Current Sheets At the sheet E x nd E z hve to be continuous, so f + (k x, k z ) = f (k x, k z ) = f (k x, k z ) We lso need u y [ H + H ] = J jk y + ψ jky ψ = J x (k x, k z ) y= where J x (k x, k z ) = 2jk + y f (k x, k z ) = J x (k x, k z ) f (k x, k z ) = 1 2jk + y J x (k x, k z ) J x (k x, k z )e jkx x e jkzz dxdz

An Importnt Exmple Apertures in Ground Plne Plne Current Sheets Consider current sheet with J z = Ilδ(x)δ(z). We lredy know this cretes mgnetic vector potentil with A z = µψ nd e jk x 2 +y 2 +z 2 ψ = Il 4π x 2 + y 2 + z 2 On the other hnd, we cn tret it s current sheet. The trnsform is J z = 1 4π 2 J x (k x, k z )e jkx x e jkzz dxdz = Il 4π 2

An Importnt Exmple Apertures in Ground Plne Plne Current Sheets Given this so tht, for y >, f (k x, k z ) = 1 2jk + y Il 4π 2 ψ(k x, y, k z ) = jil 8π 2 k y e jky y where k y = k + y = { j kx 2 + kz 2 k 2 k 2 kx 2 kz 2 for k < kx 2 + kz 2 for k > kx 2 + kz 2

An Importnt Exmple Apertures in Ground Plne Plne Current Sheets Therefore, we cn find ψ from n inverse trnsform: For y > ψ(x, y, z) = jil 8π 2 1 k y e jkx x e jky y e jkzz dxdz For y <, both y nd k y chnge sign, so we cn write ψ(x, y, z) = jil 1 8π 2 e jkx x e jky y e jkzz dxdz k y without cvet, with the understnding tht k y chnges sign t the y = plne.

The Weyl Identity Apertures in Ground Plne Plne Current Sheets Invoking equivlence (how?) we cn equte the two different expressions for ψ. This gives rise to the Weyl Identity e jkr r = 1 1 e jkx x e jky y e jkzz dxdz 2πj k y Wht does this formul men?