Green s functions for planarly layered media
|
|
- Joleen Anderson
- 6 years ago
- Views:
Transcription
1 Green s functions for planarly layered media Massachusetts Institute of Technology lecture notes Introduction: Green s functions The Green s functions is the solution of the wave equation for a point source (dipole). For scalar problems, the wave equation is written as (k 0 = ω ɛµ): and the solution for an unbounded medium is: ( 2 + k 2 0) g( r, r ) = δ( r r ), () g( r, r ) = eik 0 r r 4π r r. (2) From Maxwell s equations in frequency domain with an e iωt dependency, the wave equation for electric field Ē( r) is: Ē( r) k2 0Ē( r) = iω J( r), Therefore, the free-space dyadic Green s function satisfies = 0 for source free case. (3) which solution is G( r, r ) k 2 0G( r, r ) = I δ( r r ), (4) G( r, r ) = ( I + ) k0 2 g( r, r ), (5) (check using (Ig) = g ( g)i and Eq. ()). For the use of Green s functions in scattering problems, it is useful to express the Green s function in the same coordinates as the problem, which can be rectangular, cylindrical, spherical, etc. Here we shall concentrate on the rectangular representation (Cartesian).
2 2 Section 2. Cartesian coordinates 2 Cartesian coordinates 2. Scalar Green s function The formulae are derived from Eq. () and the Fourier transform of the quantities: g( r, r ) = (2π) 3 δ( r r ) = (2π) 3 d k e i k ( r r ) g( k), d k e i k ( r r ), (6a) (6b) where k = k xˆx + k y ŷ + k z ẑ and d k = dk x dk y dk z. Upon using Eq. (), we write: ( 2 + k0) 2 d k e i k ( r r ) g( k) = d k e i k ( r r ), (7) Introducing the differential operator ( 2 = 2 x y z 2 ) we write: ( 2 + k0) 2 d k e i k ( r r ) g( k) = d k ( 2 + k0)e 2 i k ( r r ) g( k) from which we conclude that = d k ( kx 2 ky 2 kz 2 + k0)e 2 i k ( r r ) g( k) = d k e i k ( r r ), (8) g( k) = kx 2 + ky 2 + kz 2 k0 2. (9) Using Eq. (6a), we therefore need to evaluate the following integral: g( r, r ) = (2π) 3 dk x dk y dk z k 2 x + k 2 y + k 2 z k 2 0 e i k ( r r ). (0) Note that Eq. (0) can be integrated along one of the three axis. In remote sensing application, the vertical axis is usually taken to be the z axis, (xy) being the transverse plane (planar components). We therefore choose to evaluate Eq. (0) along k z, and we split: k = k xˆx + k y ŷ + k z ẑ = k + k z ẑ, r = r + zẑ, r = r + z ẑ. (a) (b) (c)
3 3 We will perform the integral of Eq. (0) in the complex plane, using Cauchy s theorem and the Residue theorem. Before doing this, we have to be careful not to have divergent integrals. Since we integrate in k z, the condition is: If we write k z as k z = k z + ik z (k z R, k z R), we see that lim k eikzz < +. (2) z if z > 0, we have to choose k z > 0, which means that for complex plane integration, we need to deform the contour into the upper plane. if z < 0, we have to choose k z < 0, which corresponds to a deformation into the lower plane. In addition, we see from Eq. (0) that the integrand has a pole at k 2 0 z = k 2 0 k 2 x k 2 y = k 2 0 k 2. (3) We therefore need to evaluate Eq. (0) via the Residue theorem. Calculus: let us just write the integral in dk z for z > 0: + kz 2 k0 2 e i k ( r r ) = 2iπ Res z kz 2 k0 2 e i k ( r r ) z = 2iπ lim k z k 0z k z k 0z (k z k 0z )(k z + k 0z ) ei k ( r r ) e ikz(z z ) = 2iπ 2k 0z e i k ( r r ) e ik 0z (z z ), (4) so that g( r, r ) = i (2π) 2 d k 2k 0z e i k ( r r ) e ik 0z (z z ), for z z > 0. (5) The treatment for z < 0 follows the same reasoning so that we write for all (z z ): (z z ) R, g( r, r ) = i (2π) Dyadic Green s function From Eq. (6) and Eq. (26), we can get the dyadic Green s functions. d k 2k 0z e i k ( r r ) e ik 0z z z. (6) Note that is a dyadic operator (give a dyad when applied to a scalar) and can be exchanged with the integral sign. In addition, it only applies to the exponential terms so that we actually need to evaluate: e i k ( r r ) e ik 0z z z (7)
4 4 2.2 Dyadic Green s function or, by a simple change of variables: e i k r e ik 0z z (8) Calculus: Let us first consider z > 0 and Various derivatives will be: x x... x f(x, y, z) = k2 x f(x, y, z), y f(x, y, z) = k xk y f(x, y, z), and identically for z < 0. At z = 0 however, 2 z 2 eik 0z z = {ik 0z e ik 0z z z } z z (e i k r e ik 0z z ) = f(x, y, z). (9) { ( ) = ik 0z ik 0z e ik 0z 2 z z z + ik 0z e ik 0z z 2 z 2 z = 2ik 0z δ(z) k 2 0 z e ik 0z z. (20) } Using these results, we write: 2 z 2 g( r) = i (2π) 2 = δ(z) (2π) 2 = δ( r) i d k e i k r 2ik 0z δ(z) k 20z e ik 0z z 2k 0z d k e i k r i d k k 0z e i k r e ik 0z z d k k 0z e i k r e ik 0z z. (2) Again, all the other terms of applied to the integrand give k k so that the Green s function becomes: where G( r, r ) = ẑẑ δ( r) k i d k I k k k 0z k0 2 K K d k I k 0z k 2 0 e i k r e i K r for z > 0, for z < 0, (22) k = k xˆx + k y ŷ + k 0z ẑ, K = k xˆx + k y ŷ k 0z ẑ. (23a) (23b) Some notes:
5 5. The Dirac delta function is known as the singularity of the Green s function and is important in calculating the fields in the source region. 2. The different signs ensure that the integral converges for evanescent waves, i.e. when k 2 x + k 2 y > k The square bracket in the expression of the Green s functions can be expressed in terms of superposition of TE and TM waves, as we shall see. 2.3 Superposition of TE and TM waves Based on k, we can form an orthonormal system for TE/TM polarized waves: TE: ê(k 0z ) = k ẑ k ẑ = ˆxk y ŷk x = (ˆxk y ŷk x ), kx 2 + ky 2 k ρ (24a) TM: ĥ(k 0z ) = k 0 ê(k 0z ) k = k 0 z k o k ρ (ˆxk y + ŷk x ) + k ρ k o ẑ. (24b) The three vectors ˆk, ĥ and ê form an orthonormal system, in which: I = ˆkˆk + ê(k 0z )ê(k 0z ) + ĥ(k 0 z )ĥ(k 0 z ). (25) After translating to the origin, we get for the dyadic Green s functions: G( r, r ) = ẑẑ δ( r) k0 2 + i where d k ê(k 0z )ê(k 0z ) + k ĥ(k 0 )ĥ(k z 0 z ) 0z d k ê( k 0z )ê( k 0z ) + k ĥ( k 0 )ĥ( k z 0 z ) 0z e i k ( r r ) e i K ( r r ) for z > z, for z < z, (26) ê( k 0z ) = ê(k 0z ), (27a) ĥ( k 0z ) = ê( k 0 z ) K k 0. (27b) (Note that K, ê( k 0z ) and ĥ( k 0 z ) form another orthonormal set of vectors about K). 2.4 Treatment of layered media Depending upon the medium under study and the location of the source, the kernel of Eq. (26) will have to be modified. To make it more clear, we can gather the terms in the Green s function relative to the source (primed coordinates) and those relative to the source.
6 6 2.4 Treatment of layered media G( r, r )= ẑẑ δ( r) k i d k d k k 0z {ê(k 0z ) e i k r ê(k 0z ) e i k r +ĥ(k 0z ) ei k r ĥ(k 0z ) e i k r } for z > z, {ê( k k 0z ) e i K r ê( k 0z ) e i K r +ĥ( k 0z 0z ) ei K r ĥ( k 0z ) e i K r } for z < z, (28) If we now consider a layered medium problem, with an arbitrary number of layers and a source in the top region (incident wave), we write: with G( r, r ) i0 =. For z < z, i = 0: i { } dk x dk y K e ê( k 0z ) e i K r + k K h ĥ( k 0z ) e i K r 0z (29) K e =ê( k 0z ) e i K r + R T E ê(k 0z ) e k 0z e i k r K h =ĥ( k 0 z ) e i K r + R T M ĥ(k 0z ) e k 0z e i k r (30a) (30b) 2. For region l, i = l: K e =A l ê(k lz ) e i k l r + B l ê( k lz ) e i K l r K h =C l ĥ(k lz ) e i k l r + D l ĥ( k lz ) e i K l r (3a) (3b) (3c) 3. For region t, i = t: where K e =T T E ê( k tz ) e i K t r K h =T T M ĥ( k tz ) e i K t r k lz = kl 2 k2 x ky 2, k l =ˆxk x + ŷk y + ẑk lz, K l =ˆxk x + ŷk y ẑk lz, (32a) (32b) (32c) (33a) (33b) (33c) and the coefficients A l, B l, C l and D l are determined from the boundary conditions. The boundary conditions apply to the tangential electric and magnetic fields. Thus, in terms of Green s functions, we need to satisfy the continuity of ẑ G( r, r ) and ẑ
7 7 G( r, r ). Let us write this at the interface between media (l) and (l + ), by separating the TE and TM components: A l e ik lz z + B l e ik lz z = A l+ e ik lz z + B l+ e ik lz z k lz A l e ik lz z B l e ik lz z = k zl+ A l+ e ik lz z B l+ e ik lz z k lz A l e ik lz z B l e ik lz z = k zl+ A l+ e ik lz z B l+ e ik lz z k l k l+ k lz C l e ik lz z + D l e ik lz z = k zl+ C l+ e ik lz z + D l+ e ik lz z (34a) (34b) (34c) (34d) With the conditions in the first and last layer as: A 0 = R T E, B 0 =, C 0 = R T M, D 0 =, (35a) A t = 0, B t = T T E, C t = 0, D t = T T M. (35b) Before evaluating these coefficients, we can build a recursive scheme to calculate the amplitudes from region l to region l +. For example, it is straightforward to build a propagation matrix for TE modes from Eq. (34a) and (34b): ( ) ( ) A l+ e ikzl+ d l+ = V T E A l e ikzl d l (36) B l+ e ikzl+ d l+ B l e ikzl d l A similar procedure of course applied to the TM modes: ( ) ( ) C l+ e ikzl+ d l+ = V T M C l e ikzl d l D l+ e ikzl+ d l+ D l e ikzl d l In order to end up the recursive method, we have to express the reflection and transmission coefficient in the first and last regions, respectively. We shall only illustrated this point here, as it has been developed in previous classes. Let us consider a plane wave incident from region 0, with its plane of incidence parallel to the (xy) plane. All fields vectors are independent on y, so that y = 0 in Maxwell s equations. Thus, we can decompose the fields into their TE and TM components. We get in region l: TE modes: (37) z E ly, H lx = iωµ l H lz = iωµ l x E ly, ( 2 x 2 + y + ω2 ɛ l µ l ) E ly = 0. (38a) (38b) (38c)
8 8 2.4 Treatment of layered media TM modes: E lx = iωɛ l z H ly, E lz = ( 2 iωɛ l x H ly, x 2 + y + ω2 ɛ l µ l ) H ly = 0. (39a) (39b) (39c) For a TE wave inside the stratified medium: E ly =(A l e ik lz z + B l e ik lz z ) e ikxx, (40a) H lx = k l z ωµ l (A l e ik lz z B l e ik lz z ) e ikxx, (40b) H lz = k x ωµ l (A l e ik lz z + B l e ik lz z ) e ikxx. (40c) By matching the boundary conditions, and upon using the already known notation we get the recursive relation: p l(l+) = µ lk z(l+) µ l+ k lz, (4a) R l(l+) = p l(l+) + p l(l+). (4b) A l = e2iklz dl + /R2 l(l+) e 2ik (l+)z +k lz d l B l R l(l+) /R l(l+) e 2ik (l+) z d l + A, (42) l+ B l+ with the limiting condition: A t B t = 0,, A 0 B 0 = R. (43) Example: for a two-layer medium (t = 2): R = R 0 + R 2 e 2ik z (d d 0 ) + R 0 R 2 e 2ik z (d d 0 ) e2ikzd 0. (44)
FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES
Progress In Electromagnetics Research, PIER 5, 3 38, 000 FRACTIONAL DUAL SOLUTIONS AND CORRESPONDING SOURCES Q. A. Naqvi and A. A. Rizvi Communications Lab. Department of Electronics Quaidi-i-Azam University
More informationELECTROMAGNETIC WAVE INTERACTION WITH STRATIFIED NEGATIVE ISOTROPIC MEDIA
Progress In Electromagnetics Research, PIER 35, 1 52, 2002 ELECTROMAGNETIC WAVE INTERACTION WITH STRATIFIED NEGATIVE ISOTROPIC MEDIA J. A. Kong Department of Electrical Engineering and Computer Science
More informationCharacterization of Left-Handed Materials
Characterization of Left-Handed Materials Massachusetts Institute of Technology 6.635 lecture notes 1 Introduction 1. How are they realized? 2. Why the denomination Left-Handed? 3. What are their properties?
More informationTheory of Electromagnetic Nondestructive Evaluation
EE 6XX Theory of Electromagnetic NDE: Theoretical Methods for Electromagnetic Nondestructive Evaluation 1915 Scholl Road CNDE Ames IA 50011 Graduate Tutorial Notes 2004 Theory of Electromagnetic Nondestructive
More informationIntegral Equations in Electromagnetics
Integral Equations in Electromagnetics Massachusetts Institute of Technology 6.635 lecture notes Most integral equations do not have a closed form solution. However, they can often be discretized and solved
More informationEM waves: energy, resonators. Scalar wave equation Maxwell equations to the EM wave equation A simple linear resonator Energy in EM waves 3D waves
EM waves: energy, resonators Scalar wave equation Maxwell equations to the EM wave equation A simple linear resonator Energy in EM waves 3D waves Simple scalar wave equation 2 nd order PDE 2 z 2 ψ (z,t)
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationElectromagnetism HW 1 math review
Electromagnetism HW math review Problems -5 due Mon 7th Sep, 6- due Mon 4th Sep Exercise. The Levi-Civita symbol, ɛ ijk, also known as the completely antisymmetric rank-3 tensor, has the following properties:
More informationEITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity
EITN90 Radar and Remote Sensing Lecture 5: Target Reflectivity Daniel Sjöberg Department of Electrical and Information Technology Spring 2018 Outline 1 Basic reflection physics 2 Radar cross section definition
More informationDYADIC GREEN S FUNCTIONS FOR UNBOUNDED AND TWO-LAYERED GENERAL ANISOTROPIC MEDIA
Progress n Electromagnetics Research B, Vol. 30, 27 46, 2011 DYADC GREEN S FUNCTONS FOR UNBOUNDED AND TWO-LAYERED GENERAL ANSOTROPC MEDA Y. Huang and J. K. Lee Department of Electrical Engineering and
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationFiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers:
Waves can be guided not only by conductors, but by dielectrics. Fiber optics cable of silica has nr varying with radius. Simplest: core radius a with n = n 1, surrounded radius b with n = n 0 < n 1. Total
More informationGeneralized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields
PHYSICAL REVIEW E 71, 5661 5 Generalized optical theorem for reflection, transmission, and extinction of power for electromagnetic fields D. R. Lytle II Department of Electrical and Computer Engineering,
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationElectrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Electro Dynamic
Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Name Electro Dynamic Instructions: Use SI units. Short answers! No derivations here, just state your responses clearly. 1. (2) Write an
More information8.03 Lecture 12. Systems we have learned: Wave equation: (1) String with constant tension and mass per unit length ρ L T v p = ρ L
8.03 Lecture 1 Systems we have learned: Wave equation: ψ = ψ v p x There are three different kinds of systems discussed in the lecture: (1) String with constant tension and mass per unit length ρ L T v
More informationElectromagnetic Theory (Hecht Ch. 3)
Phys 531 Lecture 2 30 August 2005 Electromagnetic Theory (Hecht Ch. 3) Last time, talked about waves in general wave equation: 2 ψ(r, t) = 1 v 2 2 ψ t 2 ψ = amplitude of disturbance of medium For light,
More informationIntroduction to electromagnetic theory
Chapter 1 Introduction to electromagnetic theory 1.1 Introduction Electromagnetism is a fundamental physical phenomena that is basic to many areas science and technology. This phenomenon is due to the
More informationBound and Scattering Solutions for a Delta Potential
Physics 342 Lecture 11 Bound and Scattering Solutions for a Delta Potential Lecture 11 Physics 342 Quantum Mechanics I Wednesday, February 20th, 2008 We understand that free particle solutions are meant
More informationUniform Plane Waves Page 1. Uniform Plane Waves. 1 The Helmholtz Wave Equation
Uniform Plane Waves Page 1 Uniform Plane Waves 1 The Helmholtz Wave Equation Let s rewrite Maxwell s equations in terms of E and H exclusively. Let s assume the medium is lossless (σ = 0). Let s also assume
More informationPHY 396 K. Solutions for problems 1 and 2 of set #5.
PHY 396 K. Solutions for problems 1 and of set #5. Problem 1a: The conjugacy relations  k,  k,, Ê k, Ê k, follow from hermiticity of the Âx and Êx quantum fields and from the third eq. 6 for the polarization
More informationFourier Approach to Wave Propagation
Phys 531 Lecture 15 13 October 005 Fourier Approach to Wave Propagation Last time, reviewed Fourier transform Write any function of space/time = sum of harmonic functions e i(k r ωt) Actual waves: harmonic
More informationRECTANGULAR CONDUCTING WAVEGUIDE FILLED WITH UNIAXIAL ANISOTROPIC MEDIA: A MODAL ANALYSIS AND DYADIC GREEN S FUNCTION
Progress In Electromagnetics Research, PIER 5, 111 19, 000 RECTANGULAR CONDUCTING WAVEGUIDE FILLED WITH UNIAXIAL ANISOTROPIC MEDIA: A MODAL ANALYSIS AND DYADIC GREEN S FUNCTION S. Liu, L. W. Li, M. S.
More informationElectromagnetic Waves Across Interfaces
Lecture 1: Foundations of Optics Outline 1 Electromagnetic Waves 2 Material Properties 3 Electromagnetic Waves Across Interfaces 4 Fresnel Equations 5 Brewster Angle 6 Total Internal Reflection Christoph
More informationElectromagnetic Wave Propagation Lecture 3: Plane waves in isotropic and bianisotropic media
Electromagnetic Wave Propagation Lecture 3: Plane waves in isotropic and bianisotropic media Daniel Sjöberg Department of Electrical and Information Technology September 2016 Outline 1 Plane waves in lossless
More informationWaves in Linear Optical Media
1/53 Waves in Linear Optical Media Sergey A. Ponomarenko Dalhousie University c 2009 S. A. Ponomarenko Outline Plane waves in free space. Polarization. Plane waves in linear lossy media. Dispersion relations
More informationOptimal Electromagnetic Measurements
Optimal Electromagnetic Measurements Margaret Cheney and Gerhard Kristensson Department of Mathematical Sciences Rensselaer Polytechnic Institute Troy, NY 12180 USA chenem@rpi.edu Department of Electromagnetic
More informationEECS 117. Lecture 25: Field Theory of T-Lines and Waveguides. Prof. Niknejad. University of California, Berkeley
EECS 117 Lecture 25: Field Theory of T-Lines and Waveguides Prof. Niknejad University of California, Berkeley University of California, Berkeley EECS 117 Lecture 25 p. 1/2 Waveguides and Transmission Lines
More informationd 3 k In the same non-relativistic normalization x k = exp(ikk),
PHY 396 K. Solutions for homework set #3. Problem 1a: The Hamiltonian 7.1 of a free relativistic particle and hence the evolution operator exp itĥ are functions of the momentum operator ˆp, so they diagonalize
More informationGuided waves - Lecture 11
Guided waves - Lecture 11 1 Wave equations in a rectangular wave guide Suppose EM waves are contained within the cavity of a long conducting pipe. To simplify the geometry, consider a pipe of rectangular
More informationJoel A. Shapiro January 21, 2010
Joel A. shapiro@physics.rutgers.edu January 21, 20 rmation Instructor: Joel Serin 325 5-5500 X 3886, shapiro@physics Book: Jackson: Classical Electrodynamics (3rd Ed.) Web home page: www.physics.rutgers.edu/grad/504
More informationReview of Fundamental Equations Supplementary notes on Section 1.2 and 1.3
Review of Fundamental Equations Supplementary notes on Section. and.3 Introduction of the velocity potential: irrotational motion: ω = u = identity in the vector analysis: ϕ u = ϕ Basic conservation principles:
More informationCausality. but that does not mean it is local in time, for = 1. Let us write ɛ(ω) = ɛ 0 [1 + χ e (ω)] in terms of the electric susceptibility.
We have seen that the issue of how ɛ, µ n depend on ω raises questions about causality: Can signals travel faster than c, or even backwards in time? It is very often useful to assume that polarization
More informationELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER
Progress In Electromagnetics Research Letters, Vol. 31, 55 64, 2012 ELECTROMAGNETIC FIELD OF A HORIZONTAL IN- FINITELY LONG MAGNETIC LINE SOURCE OVER THE EARTH COATED WITH A DIELECTRIC LAYER Y.-J. Zhi
More informationPropagation of Radio Frequency Waves Through Density Filaments
PSFC/JA-15-13 Propagation of Radio Frequency Waves Through Density Filaments A. K. Ram and K. Hizanidis a May 015 a National Technical University of Athens (part of HELLAS) School of Electrical and Computer
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report February 2011
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 11-11 February 211 Notes on solving Maxwell equations Part 2: Green s function for stratified media by R. Rook
More informationEECS 117. Lecture 23: Oblique Incidence and Reflection. Prof. Niknejad. University of California, Berkeley
University of California, Berkeley EECS 117 Lecture 23 p. 1/2 EECS 117 Lecture 23: Oblique Incidence and Reflection Prof. Niknejad University of California, Berkeley University of California, Berkeley
More information5 Electromagnetic Waves
5 Electromagnetic Waves 5.1 General Form for Electromagnetic Waves. In free space, Maxwell s equations are: E ρ ɛ 0 (5.1.1) E + B 0 (5.1.) B 0 (5.1.3) B µ 0 ɛ 0 E µ 0 J (5.1.4) In section 4.3 we derived
More informationMultilayer Reflectivity
Multilayer Reflectivity John E. Davis jed@jedsoft.org January 5, 2014 1 Introduction The purpose of this document is to present an ab initio derivation of the reflectivity for a plane electromagnetic wave
More informationFields of a Dipole Near a Layered Substrate
Appendix C Fields of a Dipole Near a Layered Substrate µ z θ µ 1 ε 1 µ 2 ε 2 µ 3 ε 3 d z o x,y Figure C.1: An electric dipole with moment µ is located at r o = (,, z o ) near a layered substrate. The fields
More informationWaves. Daniel S. Weile. ELEG 648 Waves. Department of Electrical and Computer Engineering University of Delaware. Plane Waves Reflection of Waves
Waves Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Waves Outline Outline Introduction Let s start by introducing simple solutions to Maxwell s equations
More informationThe Helmholtz theorem at last!
Problem. The Helmholtz theorem at last! Recall in class the Helmholtz theorem that says that if if E =0 then E can be written as E = φ () if B =0 then B can be written as B = A (2) (a) Let n be a unit
More informationReflection/Refraction
Reflection/Refraction Page Reflection/Refraction Boundary Conditions Interfaces between different media imposed special boundary conditions on Maxwell s equations. It is important to understand what restrictions
More informationChapter 2 Basic Optics
Chapter Basic Optics.1 Introduction In this chapter we will discuss the basic concepts associated with polarization, diffraction, and interference of a light wave. The concepts developed in this chapter
More information2nd Year Electromagnetism 2012:.Exam Practice
2nd Year Electromagnetism 2012:.Exam Practice These are sample questions of the type of question that will be set in the exam. They haven t been checked the way exam questions are checked so there may
More informationH ( E) E ( H) = H B t
Chapter 5 Energy and Momentum The equations established so far describe the behavior of electric and magnetic fields. They are a direct consequence of Maxwell s equations and the properties of matter.
More informationCHAPTER 9 ELECTROMAGNETIC WAVES
CHAPTER 9 ELECTROMAGNETIC WAVES Outlines 1. Waves in one dimension 2. Electromagnetic Waves in Vacuum 3. Electromagnetic waves in Matter 4. Absorption and Dispersion 5. Guided Waves 2 Skip 9.1.1 and 9.1.2
More informationA. Bouzidi * and T. Aguili Syscom Laboratory, Engineer National School, B.P. 37, Le Belvedere 1002, Tunis, Tunisia
Progress In Electromagnetics Research M, Vol., 41 55, 01 RCS PREDICTION FROM PLANAR NEAR-FIELD MEASUREMENTS A. Bouzidi * and T. Aguili Syscom Laboratory, Engineer National School, B.P. 37, Le Belvedere
More informationPhysics 322 Midterm 2
Physics 3 Midterm Nov 30, 015 name: Box your final answer. 1 (15 pt) (50 pt) 3 (0 pt) 4 (15 pt) total (100 pt) 1 1. (15 pt) An infinitely long cylinder of radius R whose axis is parallel to the ẑ axis
More informationElectromagnetic Waves
Electromagnetic Waves Maxwell s equations predict the propagation of electromagnetic energy away from time-varying sources (current and charge) in the form of waves. Consider a linear, homogeneous, isotropic
More informationContents. 1 Basic Equations 1. Acknowledgment. 1.1 The Maxwell Equations Constitutive Relations 11
Preface Foreword Acknowledgment xvi xviii xix 1 Basic Equations 1 1.1 The Maxwell Equations 1 1.1.1 Boundary Conditions at Interfaces 4 1.1.2 Energy Conservation and Poynting s Theorem 9 1.2 Constitutive
More informationSimple medium: D = ɛe Dispersive medium: D = ɛ(ω)e Anisotropic medium: Permittivity as a tensor
Plane Waves 1 Review dielectrics 2 Plane waves in the time domain 3 Plane waves in the frequency domain 4 Plane waves in lossy and dispersive media 5 Phase and group velocity 6 Wave polarization Levis,
More informationWaveguides and Cavities
Waveguides and Cavities John William Strutt also known as Lord Rayleigh (1842-1919) September 17, 2001 Contents 1 Reflection and Transmission at a Conducting Wall 2 1.1 Boundary Conditions...........................
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationβi β r medium 1 θ i θ r y θ t β t
W.C.Chew ECE 350 Lecture Notes Date:November 7, 997 0. Reections and Refractions of Plane Waves. Hr Ei Hi βi β r Er medium θ i θ r μ, ε y θ t μ, ε medium x z Ht β t Et Perpendicular Case (Transverse Electric
More informationFundamental Limits to the Coupling between Light and 2D Polaritons by Small Scatterers SUPPORTING INFORMATION
Fundamental Limits to the Coupling between Light and D Polaritons by Small Scatterers SUPPORTING INFORMATION Eduardo J. C. Dias and F. Javier García de Abajo ICFO-Institut de Ciéncies Fotòniques The Barcelona
More information1 Fundamentals of laser energy absorption
1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms
More informationPHOTON PROPAGATOR (A REVIEW) The inhomogeneous wave equation for the four-dimensional vector
PHOTON PROPAGATOR A REVIEW I. THE VECTOR POTENTIAL The inhomogeneous wave equation for the four-dimensional vector potential A µ = A 0 = Φ, A, A 2, A 3 is 2 A µ x x 0 2 2 A µ x = J µ x where Φ is the scalar
More information1. Reminder: E-Dynamics in homogenous media and at interfaces
0. Introduction 1. Reminder: E-Dynamics in homogenous media and at interfaces 2. Photonic Crystals 2.1 Introduction 2.2 1D Photonic Crystals 2.3 2D and 3D Photonic Crystals 2.4 Numerical Methods 2.4.1
More information- HH Why Can Light Propagate in Vacuum? Hsiu-Hau Lin (Apr 1, 2014)
- HH0124 - Why Can Light Propagate in Vacuum? Hsiu-Hau Lin hsiuhau@phys.nthu.edu.tw (Apr 1, 2014) Sunlight is vital for most creatures in this warm and humid planet. From interference experiments, scientists
More informationGreen s Tensors for the Electromagnetic Field in Cylindrical Coordinates
Chapter 25 Green s Tensors for the Electromagnetic Field in Cylindrical Coordinates The solutions of the vector Helmholtz equation in three dimensions can be expressed by a complete set of vector fields
More informationClassical field theory 2012 (NS-364B) Feynman propagator
Classical field theory 212 (NS-364B Feynman propagator 1. Introduction States in quantum mechanics in Schrödinger picture evolve as ( Ψt = Û(t,t Ψt, Û(t,t = T exp ı t dt Ĥ(t, (1 t where Û(t,t denotes the
More informationCartesian Coordinates
Cartesian Coordinates Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Cartesian Coordinates Outline Outline Separation of Variables Away from sources,
More informationCalculation of Reflection and Transmission Coefficients in scuff-transmission
Calculation of Reflection and Transmission Coefficients in scuff-transmission Homer Reid May 9, 2015 Contents 1 The Setup 2 2 Scattering coefficients from surface currents 4 2.1 Computation of b(q).........................
More informationCHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS
CHAPTER 4 ELECTROMAGNETIC WAVES IN CYLINDRICAL SYSTEMS The vector Helmholtz equations satisfied by the phasor) electric and magnetic fields are where. In low-loss media and for a high frequency, i.e.,
More informationProblem Set 10 Solutions
Massachusetts Institute of Technology Department of Physics Physics 87 Fall 25 Problem Set 1 Solutions Problem 1: EM Waves in a Plasma a Transverse electromagnetic waves have, by definition, E = Taking
More informationIntroduction to Seismology
1.510 Introduction to Seismology Lecture 5 Feb., 005 1 Introduction At previous lectures, we derived the equation of motion (λ + µ) ( u(x, t)) µ ( u(x, t)) = ρ u(x, t) (1) t This equation of motion can
More informationPHYS 408, Optics. Problem Set 1 - Spring Posted: Fri, January 8, 2015 Due: Thu, January 21, 2015.
PHYS 408, Optics Problem Set 1 - Spring 2016 Posted: Fri, January 8, 2015 Due: Thu, January 21, 2015. 1. An electric field in vacuum has the wave equation, Let us consider the solution, 2 E 1 c 2 2 E =
More informationRadio Propagation Channels Exercise 2 with solutions. Polarization / Wave Vector
/8 Polarization / Wave Vector Assume the following three magnetic fields of homogeneous, plane waves H (t) H A cos (ωt kz) e x H A sin (ωt kz) e y () H 2 (t) H A cos (ωt kz) e x + H A sin (ωt kz) e y (2)
More informationRECTANGULAR MICROSTRIP RESONATOR ILLUMI- NATED BY NORMAL-INCIDENT PLANE WAVE
Progress In Electromagnetics Research, Vol. 120, 83 97, 2011 RECTANGULAR MICROSTRIP RESONATOR ILLUMI- NATED BY NORMAL-INCIDENT PLANE WAVE M. Pergol * and W. Zieniutycz Gdansk University of Technology,
More informationGreen s Function for Problems with Layers
Chapter 9 Green s Function for Problems with Layers 9. Continuity conditions The medium comprises two or more straight or curved layers. The properties of the medium are constant within each layer; they
More informationTheory of Ship Waves (Wave-Body Interaction Theory) Quiz No. 2, April 25, 2018
Quiz No. 2, April 25, 2018 (1) viscous effects (2) shear stress (3) normal pressure (4) pursue (5) bear in mind (6) be denoted by (7) variation (8) adjacent surfaces (9) be composed of (10) integrand (11)
More informationELECTROMAGNETIC WAVES
Physics 4D ELECTROMAGNETIC WAVE Hans P. Paar 26 January 2006 i Chapter 1 Vector Calculus 1.1 Introduction Vector calculus is a branch of mathematics that allows differentiation and integration of (scalar)
More informationarxiv: v1 [physics.optics] 11 Mar 2009
Brewster cross-polarization A. Aiello, 1,2 M. Merano, 1, J. P. Woerdman 1 arxiv:0903.1950v1 [physics.optics] 11 Mar 2009 1 Huygens Laboratory, Leiden University, P.O. Box 9504, 2300 RA Leiden, The Netherlands
More informationA Brief Revision of Vector Calculus and Maxwell s Equations
A Brief Revision of Vector Calculus and Maxwell s Equations Debapratim Ghosh Electronic Systems Group Department of Electrical Engineering Indian Institute of Technology Bombay e-mail: dghosh@ee.iitb.ac.in
More informationRadiation Integrals and Auxiliary Potential Functions
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
More informationChapter 1 Mathematical Foundations
Computational Electromagnetics; Chapter 1 1 Chapter 1 Mathematical Foundations 1.1 Maxwell s Equations Electromagnetic phenomena can be described by the electric field E, the electric induction D, the
More informationMultipole Expansion for Radiation;Vector Spherical Harmonics
Multipole Expansion for Radiation;Vector Spherical Harmonics Michael Dine Department of Physics University of California, Santa Cruz February 2013 We seek a more systematic treatment of the multipole expansion
More informationTypical anisotropies introduced by geometry (not everything is spherically symmetric) temperature gradients magnetic fields electrical fields
Lecture 6: Polarimetry 1 Outline 1 Polarized Light in the Universe 2 Fundamentals of Polarized Light 3 Descriptions of Polarized Light Polarized Light in the Universe Polarization indicates anisotropy
More informationPhysics 342 Lecture 30. Solids. Lecture 30. Physics 342 Quantum Mechanics I
Physics 342 Lecture 30 Solids Lecture 30 Physics 342 Quantum Mechanics I Friday, April 18th, 2008 We can consider simple models of solids these highlight some special techniques. 30.1 An Electron in a
More informationCHAPTER 2. COULOMB S LAW AND ELECTRONIC FIELD INTENSITY. 2.3 Field Due to a Continuous Volume Charge Distribution
CONTENTS CHAPTER 1. VECTOR ANALYSIS 1. Scalars and Vectors 2. Vector Algebra 3. The Cartesian Coordinate System 4. Vector Cartesian Coordinate System 5. The Vector Field 6. The Dot Product 7. The Cross
More informationPHYS4210 Electromagnetic Theory Quiz 1 Feb 2010
PHYS4210 Electromagnetic Theory Quiz 1 Feb 2010 1. An electric dipole is formed from two charges ±q separated by a distance b. For large distances r b from the dipole, the electric potential falls like
More informationThe Kubo formula of the electric conductivity
The Kubo formula of the electric conductivity We consider an electron under the influence of an electric field E(t). This system may be described by a Hamiltonian of the form Ĥ(t) = Ĥ0 e i E i (t)x i,
More informationNotes on Fourier Series and Integrals Fourier Series
Notes on Fourier Series and Integrals Fourier Series et f(x) be a piecewise linear function on [, ] (This means that f(x) may possess a finite number of finite discontinuities on the interval). Then f(x)
More informationDoes negative refraction make a perfect lens?
Does negative refraction make a perfect lens? arxiv:0809.080v1 [cond-mat.mtrl-sci] 1 Sep 008 A. G. Ramm Mathematics Department, Kansas State University, Manhattan, KS 66506-60, USA email: ramm@math.ksu.edu
More informationLecture 11 March 1, 2010
Physics 54, Spring 21 Lecture 11 March 1, 21 Last time we mentioned scattering removes power from beam. Today we treat this more generally, to find the optical theorem: the relationship of the index of
More informationQBX and the DPIE for the Maxwell Equations
University of Illinois @ Urbana-Champaign Fall 2017 - CS 598 APK The Goal For this project, the goals were to implement a model for the Python package pytential for tackling the Maxwell Equations for perfect
More informationGreen s Function for Tenuous Random Media
EECS 730, Winter 2009 c K. Sarabandi Green s Function for Tenuous Random Media In the previous section we demonstrated the application of wave theory to a tenuous medium excited by a plane wave for evaluating
More informationElectromagnetic Wave Propagation Lecture 1: Maxwell s equations
Electromagnetic Wave Propagation Lecture 1: Maxwell s equations Daniel Sjöberg Department of Electrical and Information Technology September 3, 2013 Outline 1 Maxwell s equations 2 Vector analysis 3 Boundary
More informationLecture 9. Transmission and Reflection. Reflection at a Boundary. Specific Boundary. Reflection at a Boundary
Lecture 9 Reflection at a Boundary Transmission and Reflection A boundary is defined as a place where something is discontinuous Half the work is sorting out what is continuous and what is discontinuous
More informationPhysics 442. Electro-Magneto-Dynamics. M. Berrondo. Physics BYU
Physics 44 Electro-Magneto-Dynamics M. Berrondo Physics BYU 1 Paravectors Φ= V + cα Φ= V cα 1 = t c 1 = + t c J = c + ρ J J ρ = c J S = cu + em S S = cu em S Physics BYU EM Wave Equation Apply to Maxwell
More informationWavepackets. Outline. - Review: Reflection & Refraction - Superposition of Plane Waves - Wavepackets - ΔΔk Δx Relations
Wavepackets Outline - Review: Reflection & Refraction - Superposition of Plane Waves - Wavepackets - ΔΔk Δx Relations 1 Sample Midterm (one of these would be Student X s Problem) Q1: Midterm 1 re-mix (Ex:
More informationIt is obvious that any boundary is closed. That is, the boundary of a boundary is equal to zero. For example, [kb ij] =2 [k i A j] =0.
VIOLATION OF THE GAUGE EQUIVALENCE R. I. Khrapko 1 Abstract F. V. Gubarev et al. [4] ( On the significance of the vector potential squared ) have argued that the minimum value of the volume integral of
More informationFourier Series and Integrals
Fourier Series and Integrals Fourier Series et f(x) beapiece-wiselinearfunctionon[, ] (Thismeansthatf(x) maypossessa finite number of finite discontinuities on the interval). Then f(x) canbeexpandedina
More informationRadiation by a dielectric wedge
Radiation by a dielectric wedge A D Rawlins Department of Mathematical Sciences, Brunel University, United Kingdom Joe Keller,Cambridge,2-3 March, 2017. We shall consider the problem of determining the
More informationElectromagnetism. Christopher R Prior. ASTeC Intense Beams Group Rutherford Appleton Laboratory
lectromagnetism Christopher R Prior Fellow and Tutor in Mathematics Trinity College, Oxford ASTeC Intense Beams Group Rutherford Appleton Laboratory Contents Review of Maxwell s equations and Lorentz Force
More informationScattering. 1 Classical scattering of a charged particle (Rutherford Scattering)
Scattering 1 Classical scattering of a charged particle (Rutherford Scattering) Begin by considering radiation when charged particles collide. The classical scattering equation for this process is called
More informationChapter 9. Electromagnetic waves
Chapter 9. lectromagnetic waves 9.1.1 The (classical or Mechanical) waves equation Given the initial shape of the string, what is the subsequent form, The displacement at point z, at the later time t,
More informationTHE WAVE EQUATION. The Wave Equation in the Time Domain
THE WAVE EQUATION Disturbances of materials can be described as a waveform propagating through a medium (the material itself) as particles either push or shear neighboring particles from their equilibrium
More information1 Infinite-Dimensional Vector Spaces
Theoretical Physics Notes 4: Linear Operators In this installment of the notes, we move from linear operators in a finitedimensional vector space (which can be represented as matrices) to linear operators
More information