ECE Spring Prof. David R. Jackson ECE Dept. Notes 5

Similar documents
ECE Spring Prof. David R. Jackson ECE Dept. Notes 2

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 22

EKT 345 MICROWAVE ENGINEERING CHAPTER 2: PLANAR TRANSMISSION LINES

Lecture 2 Date:

EKT 356 MICROWAVE COMMUNICATIONS CHAPTER 2: PLANAR TRANSMISSION LINES

Antennas & Propagation

ANTENNAS. Vector and Scalar Potentials. Maxwell's Equations. D = εe. For a linear, homogeneous, isotropic medium µ and ε are contant.

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

12th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 29-31,

Contribution to the cavity model for analysis of microstrip patch antennas

ECE Spring Prof. David R. Jackson ECE Dept. Notes 1

ECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson ECE Dept. Notes 13

4. Electrodynamic fields

Solutions. V in = ρ 0. r 2 + a r 2 + b, where a and b are constants. The potential at the center of the atom has to be finite, so a = 0. r 2 + b.

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

ECE Spring Prof. David R. Jackson ECE Dept. Notes 7

Review: Electrostatics and Magnetostatics

ECE Spring Prof. David R. Jackson ECE Dept. Notes 5

Mutual impedance between linear elements: The induced EMF method. Persa Kyritsi September 29th, 2005 FRB7, A1-104

Introduction to Arrays

University of Illinois at Chicago Department of Physics. Electricity & Magnetism Qualifying Examination

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

5. Plane Electromagnetic Waves

1 Spherical multipole moments

Multipole Radiation. February 29, The electromagnetic field of an isolated, oscillating source

rect_patch_cavity.doc Page 1 of 12 Microstrip Antennas- Rectangular Patch Chapter 14 in Antenna Theory, Analysis and Design (4th Edition) by Balanis

Magnetic field due to a current loop.

Experiment I Voltage Variation and Control

Notes 11 Gauss s Law II

Collaborative ASSIGNMENT Assignment 3: Sources of magnetic fields Solution

Fields and Waves I Spring 2005 Homework 8. Due: 3 May 2005

11) A thin, uniform rod of mass M is supported by two vertical strings, as shown below.

ECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 4

Analytical calculation of the power dissipated in the LHC liner. Stefano De Santis - LBNL and Andrea Mostacci - CERN

PHYS 1444 Lecture #5

Force and Work: Reminder

Magnetic Dipoles Challenge Problem Solutions

Electromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology

1 Fundamental Solutions to the Wave Equation

Antennas produce fields which add in phase at certain points of space. Consider a loop of wire that carries a current. R 1 R 2. R 2 -R 1 0.1λ. D 0.

3. Electromagnetic Waves II

Appendix B The Relativistic Transformation of Forces

Physics 506 Winter 2006 Homework Assignment #9 Solutions

Question 1: The dipole

SAMPLE PAPER I. Time Allowed : 3 hours Maximum Marks : 70

Qualifying Examination Electricity and Magnetism Solutions January 12, 2006

ECE 3318 Applied Electricity and Magnetism. Spring Prof. David R. Jackson Dept. Of ECE. Notes 20 Dielectrics

Power efficiency and optimum load formulas on RF rectifiers featuring flow-angle equations

Physics 313 Practice Test Page 1. University Physics III Practice Test II

CHAPTER 25 ELECTRIC POTENTIAL

15 Solving the Laplace equation by Fourier method

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

Chapter 3 Optical Systems with Annular Pupils

SEE LAST PAGE FOR SOME POTENTIALLY USEFUL FORMULAE AND CONSTANTS

Conventional Paper-I (a) Explain the concept of gradient. Determine the gradient of the given field: ( )

Chapter 2: Basic Physics and Math Supplements

B. Spherical Wave Propagation

Physics 2212 GH Quiz #2 Solutions Spring 2016

PY208 Matter & Interactions Final Exam S2005

LC transfer of energy between the driving source and the circuit will be a maximum.

Faraday s Law. Faraday s Law. Faraday s Experiments. Faraday s Experiments. Magnetic Flux. Chapter 31. Law of Induction (emf( emf) Faraday s Law

Gauss s Law: Circuits

[ ] = jω µ [ + jω ε E r

ELECTROSTATICS::BHSEC MCQ 1. A. B. C. D.


EM Boundary Value Problems

Coupled Electromagnetic and Heat Transfer Simulations for RF Applicator Design for Efficient Heating of Materials

Swissmetro: design methods for ironless linear transformer

1D2G - Numerical solution of the neutron diffusion equation

PHYS 110B - HW #7 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased

The condition for maximum intensity by the transmitted light in a plane parallel air film is. For an air film, μ = 1. (2-1)

#7A&B Magnetic force between two conductors

Today in Physics 218: radiation from moving charges

Electrostatics (Electric Charges and Field) #2 2010

Fields and Waves I Spring 2005 Homework 4. Due 8 March 2005

Calculate the electric potential at B d2=4 m Calculate the electric potential at A d1=3 m 3 m 3 m

Physics 122, Fall October 2012

OSCILLATIONS AND GRAVITATION

17.1 Electric Potential Energy. Equipotential Lines. PE = energy associated with an arrangement of objects that exert forces on each other

COLLISIONLESS PLASMA PHYSICS TAKE-HOME EXAM

RESONANCE SERIES RESONANT CIRCUITS. 5/2007 Enzo Paterno 1

Force between two parallel current wires and Newton s. third law

Physics 2112 Unit 14

Current, Resistance and

f(k) e p 2 (k) e iax 2 (k a) r 2 e a x a a 2 + k 2 e a2 x 1 2 H(x) ik p (k) 4 r 3 cos Y 2 = 4

Phys-272 Lecture 17. Motional Electromotive Force (emf) Induced Electric Fields Displacement Currents Maxwell s Equations

3-Axis Vector Magnet: Construction and Characterisation of Split Coils at RT. Semester Project Petar Jurcevic

School of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007

Physics 122, Fall October 2012

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

Supplementary Figure 1. Circular parallel lamellae grain size as a function of annealing time at 250 C. Error bars represent the 2σ uncertainty in

Gauss Law. Physics 231 Lecture 2-1

Graphs of Sine and Cosine Functions

PHYS 301 HOMEWORK #10 (Optional HW)

[Griffiths Ch.1-3] 2008/11/18, 10:10am 12:00am, 1. (6%, 7%, 7%) Suppose the potential at the surface of a hollow hemisphere is specified, as shown

3. Magnetostatic fields

Fields. Coulomb s Law

Homework # 3 Solution Key

cos kd kd 2 cosθ = π 2 ± nπ d λ cosθ = 1 2 ± n N db

Welcome to Physics 272

Transcription:

ECE 6345 Sping 15 Pof. David R. Jackson ECE Dept. Notes 5 1

Oveview This set of notes discusses impoved models of the pobe inductance of a coaxially-fed patch (accuate fo thicke substates). A paallel-plate waveguide model is initially assumed (at the end of the notes we will also look at the actual finite patch). a h ε, µ b x

Oveview (cont.) The following models ae investigated: Cosine-cuent model Gap-souce model Fill model Deivations ae given in the Appendix. Even moe details may be found in the efeence below. Refeence: H. Xu, D. R. Jackson, and J. T. Williams, Compaison of Models fo the Pobe Inductance fo a Paallel Plate Waveguide and a Micostip Patch, IEEE Tans. Antennas and Popagation, vol. 53, pp. 39-335, Oct. 5. 3

Cosine Cuent Model I () a h ε, µ x We assume a tube of cuent (as in Notes 4) but with a vaiation. cos I k h Note: The deivative of the cuent is eo at the top conducto (PEC). Z in I P c P c = complex powe adiated by pobe cuent P 1 = E J ds * c s S p 4

Cosine Cuent Model (cont.) Final esult: 1 1 ( ) () Zin = ( kh ) η sec kh ε Im k m(1 m) H ( k ma) J( k ma) 8 ρ + δ ε ρ ρ m= whee I m ( kh) sin( kh) = 1 + δ mo ( kh) ( mπ) k ρm = k1 mπ h 1/ Note: The wavenumbe k ρm is chosen to be a positive eal numbe o a negative imaginay numbe. kρm = k / k ρm δ m 1, m = =, m 5

Gap Souce Model a h ε, µ + - 1V x An ideal gap voltage souce of height is assumed at the bottom of the pobe. Z in I 1 6

Gap Souce Model (cont.) Final esult: Y in () 1 a H ( k ma) ρ mπ = j4π ε sinc () η h m= (1 + δm) kρmh ( kρma) h whee k ρm = k1 mπ h 1/ Note: The wavenumbe k ρm is chosen to be a positive eal numbe o a negative imaginay numbe. kρm = k / k ρm δ m 1, m = =, m 7

Fill Model a h ε, µ b x A magnetic fill of adius b is assumed on the mouth of the coax. M ˆ Eˆ ˆρE s ρ M sφ = E ρ Choose: E ρ 1 1 ρ lnb/a Z in I 1 (TEM mode of coax, assuming 1 V) 8

Fill Model (cont.) Final esult: Y in ( ) ( ) 1 1 1 H ( kρmb) H ( kρma) = j 4πε ( ) kh ln ( b/ a) η m= ( kρm) (1 + δm) H ( kρma) whee k ρm = k1 mπ h 1/ Note: The wavenumbe k ρm is chosen to be a positive eal numbe o a negative imaginay numbe. kρm = k / k ρm δ m 1, m = =, m 9

Compaison of Models Next, we show esults that compae the vaious models, especially as the substate thickness inceases. a h ε, µ b x 1

Xin (Ω) Compaison of Models Models ae compaed fo changing substate thickness. a h ε, µ b x 3 1-1 - Unifom cuent model Fill model Cosine cuent model HFFS simulation -3.5.1.15..5 h (m) ε =. a =.635 mm f = GH Z = 5 Ω (b =.19 mm) (λ = 15 cm) 11

Rin (Ω) Compaison of Models (cont.) Models ae compaed fo vaying substate thickness. a h ε, µ b x 5 45 4 35 3 5 Unifom cuent model Fill model Cosine cuent model HFFS simulation ε =. a =.635 mm f = GH Z = 5 Ω (b =.19 mm) 15 (λ = 15 cm) 1 5.5.1.15..5 h (m) λ / 4 =.5 [m] d Note: 1

Xin (Ω) Compaison of Models (cont.) Fo the gap-souce model, the esults depend on. a h ε, µ + - 1V x 3 1-1 - Delta=.1mm Delta=.1mm Delta=1mm ε =. a =.635 mm f = GH Z = 5 Ω (b =.19 mm) (λ = 15 cm) -3.5.1.15..5 h (m) 13

Rin o Xin (Ω) 5 4 3 1 Compaison of Models (cont.) The gap-souce model is compaed with the fill model, fo vaying, fo a fixed h. R h = mm Gap souce model eactance Fill model eactance Gap souce model esistance Fill model esistance a h ε, µ + - 1V ε =. a =.635 mm f = GH Z = 5 Ω (b =.19 mm) (λ = 15 cm) x -1 - X -3.1..3.4.5.6.7.8.9.1 (m) b a =.5 mm 3 14

Compaison of Models (cont.) These esults suggest the 1/3 ule: The best is chosen as = b a 3 This ule applies fo a coax feed that has a 5 Ω impedance. a h ε, µ b x a h ε, µ + - 1V x 15

Xin (Ω) Compaison of Models (cont.) The gap-souce model is compaed with the fill model, using the optimum gap height (1/3 ule). Reactance a h ε, µ b x 1-1 Gap souce model Fill model ε =. a =.635 mm f = GH Z = 5 Ω (b =.19 mm) - (λ = 15 cm) -3.5.1.15..5 h (m) 16

Rin (Ω) Compaison of Models (cont.) The gap-souce model is compaed with the fill model, using the optimum gap height (1/3 ule). 5 Resistance a h ε, µ b x 4 3 Gap souce model Fill model ε =. a =.635 mm f = GH Z = 5 Ω (b =.19 mm) 1 (λ = 15 cm).5.1.15..5 h (m) 17

Pobe in Patch A pobe in a patch does not see an infinite paallel-plate waveguide. y (x, y ) W L x Exact calculation of pobe eactance: Z jx Z cavity in p in X Im Z p in f Z in may be calculated by HFSS o any othe softwae, o it may be measued. f = fequency at which R in is maximum X cavity in 18

Pobe in Patch (cont.) Cavity Model Using the cavity model, we can deive an expession fo the pobe eactance (deivation given late) y PMC (x, y ) W e x L e This fomula assumes that thee is no vaiation of the pobe cuent o cavity fields (thinsubstate appoximation), but it does accuately account fo the actual patch dimensions. 19

Pobe in Patch (cont.) Final esult: X p ωμ μ h mπx nπy nπw p cos cos sinc L W e e m,n 1, WL e e 1 m 1 mπ nπ n δ δ 4 k L W e e 3/ wp e a a = pobe adius (x, y ) = pobe location

Pobe in Patch (cont.) Image Theoy Image theoy can be used to impove the simple paallel-plate waveguide model when the pobe gets close to the patch edge. Using image theoy, we have an infinite set of image pobes. 1

Pobe in Patch (cont.) A simple appoximate fomula is obtained by using two tems: the oiginal pobe cuent in a paallel-plate waveguide and one image. This should be an impovement when the pobe is close to an edge. η X =X + X two pobe image in in in ( ) ( ) η ( s) ( a) two 1 1 in 4 4 X = khy ka J ka khy k J k Oiginal k = k ε η = η ε / Image

Pobe in Patch (cont.) As shown on the next plot, the best oveall appoximation in obtained by using the following fomula: ( pobe two ) X = max X, X in in in modified CAD fomula 3

Xf (Ω) Pobe in Patch (cont.) Results show that the simple fomula ( modified CAD fomula ) woks faily well. 35 3 5 ( pobe two ) X = max X, X in in in Unifom cuent model Cavity model Unifom cuent model (with one image) Modified CAD fomula HFFS simulation y (x, y ) W 15 L x 1 5 x = x L/ L /.1..3.4.5.6.7.8.9 1 cente X edge 4

Appendix Next, we investigate each of the impoved pobe models in moe detail: Cosine-cuent model Gap-souce model Fill model 5

Cosine Cuent Model I () a h ε, µ x Assume that I ( ) = cos k ( h) Note: I() = cos ( kh) 6

Cosine Cuent Model (cont.) Cicuit Model: I() coax feed Z in 1 Pc = Zin I() Z in I P c 7

Cosine Cuent Model (cont.) I () a h ε, µ x Z in I P c Repesent the pobe cuent as: mπ I( ) = Im cos m= h This will allow us to find the fields and hence the powe adiated by the pobe cuent. 8

Cosine Cuent Model (cont.) Using Fouie-seies theoy: h h m π mπ m π I( )cos d = I cos cos d h h h m m= The integal is eo unless m = m. Hence h mπ h I( )cos d = Im ( 1+ δm ) h h mπ Im = I( )cos d h h ( 1 δ ) + m 9

Cosine Cuent Model (cont.) o h mπ Im = cos k( h)cos d h + h ( 1 δ ) mo Result: I m ( kh) sin( kh) = 1 + δ mo ( kh) ( mπ) (deivation omitted) 3

Cosine Cuent Model (cont.) Note: We have both E and E ρ To see this: E = (Time-Hamonic Fields) 1 1 Eφ E ( ρe ) ρ + + = ρ ρ ρ φ so E ρ 31

Cosine Cuent Model (cont.) Fo E, we epesent the field as follows: ρ < a mπ E A J k = mcos m= h ρm ( ρ ) mπ = ρ > a + + () E A cos H ( k ρ ) m ρm m= h whee k ρm k = mπ h ( k km ) = 1/ 1/ 3

Cosine Cuent Model (cont.) At ρ = a E E = (BC 1) + so ( ) = ( ) A H k a A J k a + () m ρm m ρm 33

Cosine Cuent Model (cont.) Also we have H H = J φ φ1 s (BC ) whee H φ 1 Eρ E = jωµ ρ To solve fo E ρ, use H = jωε E 34

Cosine Cuent Model (cont.) so jωε E E ρ ρ 1 H H = ρ φ = 1 H jωε φ φ Hence we have H φ 1 1 Hφ E = jωµ jωε ρ Fo the m th Fouie tem: 1 1 ( ) E H km H jωµ jωε ρ ( m) ( m) ( m) φ = φ 35

Cosine Cuent Model (cont.) so that E ρ ( m) ( m) ( m) φ m φ = ωε kh k H j whee k k = k ρ m m Hence H ( m) φ = jωε E k ρ ρm ( m) 36

Cosine Cuent Model (cont.) H H = J = φ φ1 s I π a Fo the m th Fouie tem: ( ) ( ) ( ) H H = J = m m m m φ φ1 s I π a whee ( m) H ( m) φ = jωε E k ρ ρm 37

Cosine Cuent Model (cont.) Hence jωε + () ( ) ( ) k ( ) m AmH ρ kρma AmJ kρma = k ρm I m π a ( ) = ( ) A H k a A J k a + () m ρm m ρm Using (BC 1) we have () H () ( k ma) k + ρ Im ρm AH m ( k ma) A + ρ m J ( kpma) = J( kρma) π a jωε 38

Cosine Cuent Model (cont.) o k + () () Im m A ρ m J( k ma) H ρ ( kρma) H ( k ma) J ( k ma) ρ ρ = J( kρma) π a jωε o I k m m A m j + ρ = J( kρma) π kρma π a jωε Hence (using the Wonskian identity) 1 k + ρm Am = Im J kρma 4 ωε ( ) 39

Cosine Cuent Model (cont.) We now find the complex powe adiated by the pobe: P c 1 * = E J sds = 1 s π h h E ( a) J * s = π a E ( a) J d * s a d dφ π a h * = E ( ) ( ) a I d π a 1 h + () mπ * m π = A mh ( kρma)cos Im cos d m= h m' = h 4

Cosine Cuent Model (cont.) Integating in and using othogonality, we have: 1 P AI H k a h + * () c= m m ( ρm ) 1+ m m= ( δ ) = + k A m + coefficient h () 1 ρm * ( 1 δm) H ( kρma) Im J( kρma) Im 4 m= 4 ωε Hence, we have: h 1 () Pc =+ Im kpm(1 + δm) H ( kρma) J( kρma) 16 ωε m= 41

Cosine Cuent Model (cont.) Z in = P c cos ( kh) Theefoe, h 1 () Zin = sec ( kh) Im kρm ( 1 + δm) H ( kρma) J( kρma) 8 ωε m= Define: k ρm = k ρm k mπ = εµ kh 4

Cosine Cuent Model (cont.) Also, use k = ωε ω µ ωε ε = η ε We then have 1 1 ( ) () Zin = ( kh ) η sec kh ε Im kρm(1 + δm) H ( kρma) J( kρma) 8 ε m= The pobe eactance is: X = Im( Z ) p in 43

Cosine Cuent Model (cont.) Thin substate appoximation 1 1 ( ) () Zin = ( kh ) η sec kh ε Im k m(1 m) H ( k ma) J( k ma) 8 ρ + δ ε ρ ρ m= kh<< 1: I m Keep only the m = tem ( kh) sin( kh) = 1 + δ mo ( kh) ( mπ) The esult is 1 () Zin η ( kh) µ J( ka) H ( ka) 4 (same as pevious esult using unifom model) 44

Gap Model a h ε, µ + - 1V x ( ( ) ) ( ) mπ E, ρ = Bm H kρmρ cos m= h E ( a, ) 1/, < < =, othewise. Note: It is not clea how best to choose, but this will be e-visited late. 45

Gap Model (cont.) a h ε, µ + - 1V x At ρ = a: ( ( ) ) mπ 1/, < < E a, = Bm H ( kρma) cos = m= h, othewise. Fom Fouie seies analysis (details omitted): B m mπ = sinc h + H k a h ( 1 δ ) m ( ) ( ρm ) 46

Gap Model (cont.) a h ε, µ + - 1V x Y in = π ajs ( ) whee J ( ) = H ( ) s φ The magnetic field is found fom E, with the help of the magnetic vecto potential A (the field is TM ): 1 A H φ = µ ρ Use: ( ( ) ) mπ A, ρ = Am H ( kρmρ) cos m= h whee 1 jωµε E = + k A Setting ρ = a allows us to solve fo the coefficients A m. 47

Gap Model (cont.) a h ε, µ + - 1V x Final esult: () 1 a H ( k ma) ρ mπ Yin = j4π k sinc η () h m= (1 + δ m) kρmh ( kρma) h 48

a Fill Model h ε, µ To find the cuent I (), use ecipocity. b x M s φ 1 1 = ρ ln ( b/ a) Z in 1V fill I 1 Intoduce a ing of magnetic cuent K = 1 in the φ diection at (the testing cuent B ). A B S F a b I H M dv A,B B, A V V S F H H b b a M dv M ds s 49

Fill Model (cont.) A B S F ( ) b I = H M s ds S = F π b b Hφ (,) Ms φ d d a b a φ ( b a) ( ) b 1 1 b = π Hφ ( ρ,) ρdρ ρ ln ( b/ a a ) π = ln / ρ ρ ρ φ b = π H ρ, M ρdρ b a H b φ sφ ( ρ,) dρ 5

Fill Model (cont.) The magnetic cuent ing B may be eplaced by a 1V gap souce of eo height (by the equivalence pinciple). b π gap I( ) = Hφ ( ρ,) dρ ln / ( b a) b a Let : b The field of the gap souce is then calculated as was done in the gapsouce model, using =. 51

Fill Model (cont.) a h ε, µ b x Final esult: ( ) ( ) 1 1 H ( kρmb) H ( kρma) Yin = j 4π k ( ) hη ln ( b/ a) m= ( kρm) (1 + δm) H ( kρma) 5