EE 5337 Computational Electromagnetics
|
|
- Kathleen Holland
- 6 years ago
- Views:
Transcription
1 Instucto D. Ramond Rumpf (95) EE 5337 Computational Electomagnetics Lectue # WEM Etas Lectue These notes ma contain copighted mateial obtained unde fai use ules. Distibution of these mateials is stictl pohibited Slide Outline Using 3D WEM fo D and D Analsis, and using D WEM fo D Analsis Visualiing the fields Fomulation of efficient 3D WEM Fomulation of efficient D WEM WEM with anisotopic mateials Incopoating dispesion Reduced Bloch mode epansion (RBME) technique Genealied smmet Lectue Slide
2 Using 3D WEM fo D and D Analsis, and Using D WEM fo D Analsis Lectue Slide 3 Using 3D WEM fo 3D, D & D Analsis 3D Analsis with 3D WEM D Analsis with 3D WEM D Analsis with 3D WEM along Q along R along along Q along R along along Q along R along Conventional 3D WEM The device is unifom along the ais so R can be set to. The device is unifom along the and aes so and Q can both be set to. Lectue Slide 4
3 Using D WEM fo D & D Analsis D Analsis with D WEM D Analsis with D WEM along Q along along Q along Conventional D WEM The device is unifom along the ais so Q can be set to. Lectue Slide 5 Visualiing the Fields Lectue Slide 6 3
4 Eigen Vectos Contain Field Infomation Fo D WEM, we solved the following eigen value poblem. k V K K K K s s Recall the field epansion into a plane wave basis jk p q eigen vectos λ eigen values,,, p ˆ q E S p q e k p q ˆ p q s The comple amplitudes of the plane waves ae stoed in the elements of the eigen vectos. Lectue Slide 7 Reaange Tems The field epansion is,,, p ˆ q E S p q e k p q ˆ jk p q p q j We substitute in the epession fo k p, q and facto out. p q j ˆ ˆ E e S p q e j, p q This is just an invese D FFT. We can now wite the electic field epession as j E e FFT Sp, q e Lectue Slide 8 4
5 Visualiing the Data S p, q S p, q on full gid q q p p j e Sp q FFT, E Lectue Slide 9 ocedue ( of ) Step Etact eigen vecto and eshape Step Inset Fouie coefficients into lage gid. s = V(:,m); s = eshape(s,,q); nc = ceil(n/); n = nc - floo(/); n = nc + floo(/); nc = ceil(n/); n = nc - floo(q/); n = nc + floo(q/); sf = eos(n,n); sf(n:n,n:n) = s; Lectue Slide 5
6 ocedue ( of ) Step 3 Calculate invese FFT of S p, q. S p, q on full gid A, on full gid Step 4 Calculate phase tem. q p a = ifft(ifftshift(sf)); a = a/ma(abs(a(:))); phase = ep(-i*(beta()*x + beta()*y)); Step 5 Calculate oveall field..* = E = phase.*a; Lectue Slide Fomulation of Efficient 3D lane Wave Epansion Method Lectue Slide 6
7 Notation Unbolded lettes ae scala quantities. Bold capital lettes ae matices, usuall squae matices. Bold lowe case lettes ae column vectos. a.678 a a a a a a A a a a M N N M M MN Bold lowe case lettes with a vecto sign ae column vectos of column vectos. e e e e a a a Lectue Slide 3 Vecto Mati Opeations Definitions A diagonal mati containing components of A A diagonal mati containing components of A A diagonal mati containing components of A A A i i A i AiN A A A A Vecto Mati Opeations A A A A A A A A A A A A A A A Unless A A B B A, A, A, B, B, and B ae all diagonal. Unless A A B B A, A, A, B, B, and B ae all diagonal. In ou fomulation, onl [ ] and [ ] ae not diagonal matices. Lectue Slide 4 7
8 Epand the Field Into Two Othogonal olaiations We wite the field as the sum of two othogonal polaiations. s s ˆ s ˆ Two polaiation vectos ae chosen fo each spatial hamonic in the plane wave epansion. These must be othogonal to the diection of the hamonic. If these ae chosen coectl, the following equations will be satisfied. ˆ ˆ k ˆ ˆ k K K k ˆ ˆ ˆ ˆ k K K k ˆ ˆ ˆ ˆ k K K Lectue Slide 5 Recipe fo Calculating Othogonal olaiation Vectos Given the wave vecto of the i th spatial hamonic k k aˆ k aˆ k aˆ i, i, i, i We constuct anothe vecto that is in a diffeent diection than k. i v 4k aˆ k aˆ 3k aˆ i, i, i, i We can now calculate two vectos that ae pependicula to the wave vecto. kivi pˆ, i kivi An choice of pˆ, i and pˆ, i can be used as long as pˆ ˆ, i p, i k i. k ˆ i p ehaps thee eists a bette choice than descibed hee., i pˆ, i k pˆ i, i Lectue Slide 6 8
9 Wave Equation with olaiation Epansion We substitute the polaiation epansion into ou mati wave equation. K K s k s K ˆ ˆ ˆ ˆ K s s k s s Lectue Slide 7 ojecting a Vecto Mati Equation Onto Two Othogonal olaiations We can poject one vecto mati onto anothe using ab oj a onto b b b ab oj a onto b a b b b b a Using this equation, we can poject both sides of a vectomati equation onto ou two othogonal polaiations to get two sepaate equations. A B ˆ ˆ A ˆ B ˆ A B ˆ ˆ A ˆ B ˆ A B One 33 mati equation Two mati equations Note: ˆ I ˆ Lectue Slide 8 9
10 oject Ou Vecto Mati Equation Onto Two Othogonal olaiations We appl the esults fom the last slide to ou vecto mati equation. K K s s s s ˆ ˆ ˆ ˆ k oject onto ˆ ˆ K K s s s s ˆ ˆ ˆ ˆ ˆ k oject onto ˆ ˆ K K s s s s ˆ ˆ ˆ ˆ ˆ k We now have two mati equations. Lectue Slide 9 Block Mati Fom We can put ou two pevious equations into block mati fom as follows. ˆ ˆ ˆ ˆ ˆ ˆ K K s s s s ˆ ˆ ˆ ˆ ˆ ˆ K K s s s s k k ˆ ˆ ˆ ˆ K K K K s ˆ k ˆ ˆ ˆ ˆ K K K K s ˆ ˆ ˆ s ˆ ˆ ˆ ˆ s This is a genealied eigen value poblem. It is smalle than the 33 genealied eigen value poblem that we deived peviousl. Lectue Slide
11 Final Fom of Genealied Eigen Value oblems Reduced Eigen Value oblem fo Electic Fields ˆ ˆ ˆ ˆ s s k K K K K A B A s s ˆ ˆ ˆ ˆ K K K K ˆ ˆ ˆ ˆ B ˆ ˆ ˆ ˆ Reduced Eigen Value oblem fo Magnetic Fields ˆ ˆ ˆ ˆ u u k K K K K A B A u u ˆ ˆ ˆ ˆ K K K K ˆ ˆ ˆ ˆ B ˆ ˆ ˆ ˆ Note that both of these eigen value poblems ae valid fo geneal anisotopic media. Lectue Slide Genealied Eigen Value oblem with NO MAGNETIC RESONSE Genealied Eigen Value oblem fo Electic Fields ˆ ˆ ˆ ˆ s s K K A k B A ˆ ˆ ˆ ˆ s s K K ˆ ˆ ˆ ˆ B ˆ ˆ ˆ ˆ Odina Eigen Value oblem fo Magnetic Fields ˆ ˆ ˆ ˆ u u k K K K K A A u u ˆ ˆ ˆ ˆ K K K K This is the most common fomulation because it is an odina eigen value poblem. Note that both of these eigen value poblems ae valid fo geneal anisotopic media. Lectue Slide
12 Recoveing s fom s and s The genealied eigen value poblem fo the electic fields is s s A kb s s Let the eigen vecto and eigen value matices be s s W eigen-vecto mati A kb s s λ eigen-value mati If we etact a single eigen mode, w i and i, we have si wi si Recall that s s ˆ ˆ s. Theefoe, s ˆ ˆ ˆ ˆ si s ˆ ˆ i s is i wi s si s These tems ae now independent of the chosen polaiation vectos. Lectue Slide 3 Fomulation of Efficient D lane Wave Epansion Method Lectue Slide 4
13 Reduction to One Dimension Fo the D devices whee popagation is onl in the diection, the wave has no wave vecto components in the and diections. K and K Mawell s equations educe to Ku jk Ku jk s s jk s Ks jk Ks jk u u jk u Immediatel, we see that s u The field will have no components in the longitudinal diection. Lectue Slide 5 Two Modes Mawell s equations have again decoupled into two distinct sets of equations. Ku K u jk jk jk s s s Ks K s jk jk jk u u u Ku E Mode jk Ks jk E Mode s Ku jk u Ks jk s u Lectue Slide 6 3
14 Final Eigen Value oblems We can substitute one equation into the othe to deive two eigenvalue poblems. E Mode Kμ Ks k u Kε Ku jk o k u s jk s K s K u E Mode Kμ Ks k s u jk Kε Ku o k s u K s K u jk Lectue Slide 7 WEM with Full Anisotopic Mateials Lectue Slide 8 4
15 Anisotopic 3D WEM The mati equations deived fo 3D WEM ae valid fo full anisotopic mateials. In the anisotopic case, howeve, the mateial tensos become block matices composed of convolution matices fo each tenso element individuall. Lectue Slide 9 Anisotopic D WEM In geneal, all field components emain coupled in anisotopic media. Since Mawell s equations will not simplif in this case, thee is almost no numeical advantage to developing a D WEM fo anisotopic media. Instead, use ou anisotopic 3D WEM fo D poblems b using onl on spatial hamonic in the unifom diection. infinite and unifom ais # hamonics along Q # hamonics along R # hamonics along In this case, ou would set R =. Lectue Slide 3 5
16 Incopoating Dispesion Lectue Slide 3 Techniques fo Incopoating Dispesion. Iteative WEM Can onl handle one fequenc at a time.. Do not use k as the eigen value Lectue Slide 3 6
17 Altenate D WEM Fomulation ( of 5) Recall ou stating point fo E and H modes E Mode Ku Ku jk s Ks Ks jk jk u u H Mode Ks Ks jk u Ku Ku We will choose as ou eigen value instead of k. jk jk s s jk jk I G u I G u s E Mode I G s u IG s jk u jk jk I G s I G s u H Mode I G u s IG u jk s G G G G components of wave vecto epansion components of wave vecto epansion Lectue Slide 33 Altenate D WEM Fomulation ( of 5) Net we nomalie ou wave vecto tems b dividing b k. K k K K k K G k G G k G k k Ou govening equations can now be witten as E Mode IG u IG u j s j j I G s u I G s u H Mode j j j I G s I G s u I G u s I G u s Lectue Slide 34 7
18 Altenate D WEM Fomulation (3 of 5) Solving fo s and u, we get u I G s We eliminate these tems b substituting these epessions into ou govening equations. E Mode j IG IG s IG u js IG s j u H Mode j E Mode j H Mode s I G u j j j I G I G u I G s u I G u s Lectue Slide 35 Altenate D WEM Fomulation (4 of 5) We eaange these equations to bing to the ight hand side. E Mode H Mode Gu j IG IG s u j u G s s Gs j IG IG u s j s G u u Lectue Slide 36 8
19 Altenate D WEM Fomulation (5 of 5) Finall, we wite ou equations in block mati fom. E Mode G j IG IG u u j s s G H Mode G j IG IG s s j u u G Lectue Slide 37 Altenate 3D WEM Fomulation ( of ) We choose to be the eigen value. j j Ku IG u s IG u K u s Ku Ku j s Solve fo longitudinal components. Eliminate longitudinal components. j j Ks IG s u Lectue Slide 38 IG s K s u Ks Ks j u jk K s K K s G u u j s j K u K u u j K s K s j K K s jk K s G u u j j j j j Gs K Ku K K u s Gs K K u K Ku s 9
20 Altenate 3D WEM Fomulation ( of ) In block mati fom, we have j j j j j G K K K K s s G j K K K K s s j K K K K G u u u u j K K K K G Lectue Slide 39 Altenate D WEM Fomulation Fo D, ou eigen value poblem is E Mode I G u j s j Ku jk Gs j u s j s G u u G j s s j u u G Ks jk u I G s u s E Mode j j Ks Ku jk s jk Gs j u s j s G u u G j s s j u u G I G u s IG s u u Lectue Slide 4
21 Reduced Bloch Mode Epansion (RBME) Technique Mahmoud I. Hussein, Reduced Bloch mode epansion fo peiodic media band stuctue calculations, oc. R. Soc. A 465, pp , 9. Lectue Slide 4 The oblem Suppose we wish to calculate a photonic band diagam o calculate the full bands thoughout the entie Billouin one. Fo D simulations, tpical mati sie is 4 4. Fo 3D simulations, tpical mati sie is,,. These lage matices must be solved fo each Bloch wave vecto of inteest. This is tpicall s of vectos. a c Run time is long and computation is bogged down. a a a Lectue 4 a
22 Step : Calculate the Full Solution at Onl the Ke oints of Smmet Γ M Notes Ensue that no edundant points ae used o RBME will fail. An asmmetic distibution of ke points will lead to asmmet in the data calculated fom them. You can impove accuac b adding moe ke points, but this will be slowe and less efficient. X 3 N V Γ v Γ v Γ v Γ v Γ 3 N V X v X v X v X v X V M 3 N v M v M v M v M Done? Constuct full WEM Eigen Value oblem Solve full WEM Eigen Value oblem Loop though ke points of smmet Append Tansfomation Mati U Lectue 43 Step : Combine Eigen Vecto Matices Using Lowest Ode Modes 3 N V Γ vγ vγ v Γ vγ 3 N VX vx vx v X vx VM 3 N v M vm v M vm Take the M lowest ode modes fom each eigen vecto mati and constuct a new eigen vecto mati U. M M M U v Γ vγ v X vx v M vm U does not have to be a squae mati. Lectue 44
23 Step 3: efom Gam Schmidt Othonomaliation on the Mati U M M M U v Γ vγ v X vx v M vm M M M U v Γ v Γ v X v X v M v M Lectue 45 Gam Schmidt Othonomaliation U v v v 3 v M vm vm v Step : v v v vvv Step : v v vvv v v v v v v v Step 3: v 3 v v v v v v v This algoithm essentiall just emoves the components of the pevious vectos and then nomalies the amplitude. Lectue 46 3
24 Step 4: Constuct Eigen Value oblem fo the Net Constuct the standad WEM eigen value poblem fo the net, but do not solve it et. A λb E Mode H Mode K K K K A K K K K B E Mode H Mode Lectue 47 Step 5: Reduce Mati Sie b Epanding into Bloch Modes as the Basis efom a linea tansfomation that epesses the field in tems of the new basis fomed fom the educed set of Bloch modes. This damaticall educes the sie of the matices while etaining ecellent accuac because those Bloch modes moe efficientl descibe the fields in the stuctue. A B H U AU H U BU H A U Lectue 48 A U 4
25 Step 6: Solve the Reduced Eigen Value oblem The educed eigen value poblem is solved accoding to A λb V, λ The eigen values ae intepeted the same as in the standad WEM. If the ae needed, the eigen vectos must be tansfomed back into a plane wave basis to be compatible with the standad WEM. V UVU H Lectue 49 Step 7: Repeat ocedue fo Each Bloch Wave Vecto in IBZ Done? Constuct full WEM Eigen Value oblem Tansfom to U Solved Reduced Eigen Value oblem Recod Loop though list of s Lectue 5 5
26 Block Diagam of RBME Dashboad Build Unit Cell Build Convolution Matices Constuct IBZ and Ke oints Geneate List of s Done? Constuct WEM Eigen Value oblem Solve WEM Eigen Value oblem Loop though ke points of smmet Append Tansfomation Mati U Gam Schmidt Done? Constuct WEM Eigen Value oblem Tansfom to U Solved Reduced Eigen Value oblem Recod Loop though list of s Daw Band Diagam Lectue 5 Genealied Smmet Lectue Slide 5 6
27 Mste Bands When Modeling Heagonal Aas with Rectangula Smmet Unit Cell 3 Recipocal Lattice Vectos a c Γ M K Γ Unit Cell 3 Recipocal Lattice Vectos a c Mste nonphsical bands Γ M K Γ Lectue Slide 53 Coect phsical bands Cause of Mste Bands Lattice smmet imposes phase conditions between adjacent unit cells epessed b the Bloch theoem. The highe the smmet, the tighte the phase conditions that ae imposed. Thee ae fewe phase conditions with a ectangula unit cell so thee ae additional allowed states. Lectue Slide 54 7
28 The Fi The fi is to constuct the convolution matices fo the geneal smmet case. The standad FFT can no longe be used and this pocess is moe numeicall intensive. Fo D peiodic functions, the epansion is jptqt jptqt f, apq, e apq, fe, da A p q A Fo 3D peiodic functions, this is jptqtt3 jptqtt3 f apq,, e apq,, f e dv V p q V These must be numeicall evaluated. Lectue Slide 55 8
Transfer Matrix Method
9/6/17 Instucto D. Ramond Rumpf (915) 747 6958 cumpf@utep.edu 5337 Computational lectomagnetics Lectue #4 Tansfe Mati Method Lectue 4These notes ma contain copighted mateial obtained unde fai use ules.
More informationEE 5337 Computational Electromagnetics. Maxwell s Equations
9/15/217 Instucto D. Ramond Rumpf (915) 747 6958 cumpf@utep.edu 5337 Computational lectomagnetics Lectue #2 Mawell s quations Lectue 2These notes ma contain copighted mateial obtained unde fai use ules.
More informationThe Perfectly Matched Layer (PML)
9/13/216 EE 533 Electomagnetic Analsis Using Finite Diffeence Time Domain Lectue #13 The Pefectl Matched Lae (PML) Lectue 13 These notes ma contain copighted mateial obtained unde fai use ules. Distibution
More informationMaxwell s Equations in Fourier Space
Instucto D. Ramond Rumpf (915) 747 6958 cumpf@utep.edu EE 5337 Computational Electomagnetics Lectue #18 Mawell s Equations in Fouie Space Lectue 18 These notes ma contain copighted mateial obtained unde
More informationImplementation of RCWA
Instucto D. Ramond Rumpf (915) 747 6958 cumpf@utep.edu EE 5337 Computational Electomagnetics Lectue # Implementation of RCWA Lectue These notes ma contain copighted mateial obtained unde fai use ules.
More informationFormulation of Rigorous Coupled Wave Analysis
1/8/18 Instucto D. Ramon Rumpf (915) 747 6958 cumpf@utep.eu EE 5337 Computational Electomagnetics Lectue #1 Fomulation of Rigoous Couple Wave Analsis Lectue 1 These notes ma contain copighte mateial obtaine
More informationReading Assignment. Problem Description for Homework #9. Read Chapters 29 and 30.
Reading Assignment Read Chaptes 29 and 30. Poblem Desciption fo Homewok #9 In this homewok, you will solve the inhomogeneous Laplace s equation to calculate the electic scala potential that exists between
More informationSEE LAST PAGE FOR SOME POTENTIALLY USEFUL FORMULAE AND CONSTANTS
Cicle instucto: Moow o Yethiaj Name: MEMORIL UNIVERSITY OF NEWFOUNDLND DEPRTMENT OF PHYSICS ND PHYSICL OCENOGRPHY Final Eam Phsics 5 Winte 3:-5: pil, INSTRUCTIONS:. Do all SIX (6) questions in section
More informationMAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS
The 8 th Intenational Confeence of the Slovenian Society fo Non-Destuctive Testing»pplication of Contempoay Non-Destuctive Testing in Engineeing«Septembe 1-3, 5, Potoož, Slovenia, pp. 17-1 MGNETIC FIELD
More informationPHYS 705: Classical Mechanics. Small Oscillations
PHYS 705: Classical Mechanics Small Oscillations Fomulation of the Poblem Assumptions: V q - A consevative system with depending on position only - The tansfomation equation defining does not dep on time
More informationSupporting Information Wedge Dyakonov Waves and Dyakonov Plasmons in Topological Insulator Bi 2 Se 3 Probed by Electron Beams
Suppoting Infomation Wedge Dyakonov Waves and Dyakonov Plasmons in Topological Insulato Bi 2 Se 3 Pobed by Electon Beams Nahid Talebi, Cigdem Osoy-Keskinboa, Hadj M. Benia, Klaus Ken, Chistoph T. Koch,
More informationChapter Eight Notes N P U1C8S4-6
Chapte Eight Notes N P UC8S-6 Name Peiod Section 8.: Tigonometic Identities An identit is, b definition, an equation that is alwas tue thoughout its domain. B tue thoughout its domain, that is to sa that
More information1D2G - Numerical solution of the neutron diffusion equation
DG - Numeical solution of the neuton diffusion equation Y. Danon Daft: /6/09 Oveview A simple numeical solution of the neuton diffusion equation in one dimension and two enegy goups was implemented. Both
More informationAs is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.
Appendix A Vecto Algeba As is natual, ou Aeospace Stuctues will be descibed in a Euclidean thee-dimensional space R 3. A.1 Vectos A vecto is used to epesent quantities that have both magnitude and diection.
More informationRigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018
Rigid Body Dynamics 2 CSE169: Compute Animation nstucto: Steve Rotenbeg UCSD, Winte 2018 Coss Poduct & Hat Opeato Deivative of a Rotating Vecto Let s say that vecto is otating aound the oigin, maintaining
More informationSimulation of a 2-link Brachiating Robot with Open-Loop Controllers
Simulation of a -link Bachiating Robot with Open-oop Contolles David Uffod Nothwesten Univesit June 009 . Poject Oveview The goal of this poject was to wite a complete simulation of a -link swinging obot
More informationPerfectly Matched Layer
7/2/217 Intucto D. Ramond Rumpf (915) 747 6958 cumpf@utep.edu 5337 Computational lectomagnetic Lectue #9 Pefectl Matched Lae Lectue 9Thee note ma contain copighted mateial obtained unde fai ue ule. Ditibution
More informationworking pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50
woking pages fo Paul Richads class notes; do not copy o ciculate without pemission fom PGR 2004/11/3 10:50 CHAPTER7 Solid angle, 3D integals, Gauss s Theoem, and a Delta Function We define the solid angle,
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationTHE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN
THE INFLUENCE OF THE MAGNETIC NON-LINEARITY ON THE MAGNETOSTATIC SHIELDS DESIGN LIVIU NEAMŢ 1, ALINA NEAMŢ, MIRCEA HORGOŞ 1 Key wods: Magnetostatic shields, Magnetic non-lineaity, Finite element method.
More informationElectromagnetic scattering. Graduate Course Electrical Engineering (Communications) 1 st Semester, Sharif University of Technology
Electomagnetic scatteing Gaduate Couse Electical Engineeing (Communications) 1 st Semeste, 1390-1391 Shaif Univesity of Technology Geneal infomation Infomation about the instucto: Instucto: Behzad Rejaei
More informationPart V: Closed-form solutions to Loop Closure Equations
Pat V: Closed-fom solutions to Loop Closue Equations This section will eview the closed-fom solutions techniques fo loop closue equations. The following thee cases will be consideed. ) Two unknown angles
More informationMath 124B February 02, 2012
Math 24B Febuay 02, 202 Vikto Gigoyan 8 Laplace s equation: popeties We have aleady encounteed Laplace s equation in the context of stationay heat conduction and wave phenomena. Recall that in two spatial
More informationSchool of Electrical and Computer Engineering, Cornell University. ECE 303: Electromagnetic Fields and Waves. Fall 2007
School of Electical and Compute Engineeing, Conell Univesity ECE 303: Electomagnetic Fields and Waves Fall 007 Homewok 8 Due on Oct. 19, 007 by 5:00 PM Reading Assignments: i) Review the lectue notes.
More informationConjugate Gradient Methods. Michael Bader. Summer term 2012
Gadient Methods Outlines Pat I: Quadatic Foms and Steepest Descent Pat II: Gadients Pat III: Summe tem 2012 Pat I: Quadatic Foms and Steepest Descent Outlines Pat I: Quadatic Foms and Steepest Descent
More informationTopic 4a Introduction to Root Finding & Bracketing Methods
/8/18 Couse Instucto D. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: cumpf@utep.edu Topic 4a Intoduction to Root Finding & Backeting Methods EE 4386/531 Computational Methods in EE Outline
More informationLecture 8 - Gauss s Law
Lectue 8 - Gauss s Law A Puzzle... Example Calculate the potential enegy, pe ion, fo an infinite 1D ionic cystal with sepaation a; that is, a ow of equally spaced chages of magnitude e and altenating sign.
More informationModule 9: Electromagnetic Waves-I Lecture 9: Electromagnetic Waves-I
Module 9: Electomagnetic Waves-I Lectue 9: Electomagnetic Waves-I What is light, paticle o wave? Much of ou daily expeience with light, paticulaly the fact that light ays move in staight lines tells us
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationEM Boundary Value Problems
EM Bounday Value Poblems 10/ 9 11/ By Ilekta chistidi & Lee, Seung-Hyun A. Geneal Desciption : Maxwell Equations & Loentz Foce We want to find the equations of motion of chaged paticles. The way to do
More information3. Electromagnetic Waves II
Lectue 3 - Electomagnetic Waves II 9 3. Electomagnetic Waves II Last time, we discussed the following. 1. The popagation of an EM wave though a macoscopic media: We discussed how the wave inteacts with
More informationSection 8.2 Polar Coordinates
Section 8. Pola Coodinates 467 Section 8. Pola Coodinates The coodinate system we ae most familia with is called the Catesian coodinate system, a ectangula plane divided into fou quadants by the hoizontal
More informationrect_patch_cavity.doc Page 1 of 12 Microstrip Antennas- Rectangular Patch Chapter 14 in Antenna Theory, Analysis and Design (4th Edition) by Balanis
ect_patch_cavit.doc Page 1 of 1 Micostip Antennas- Rectangula Patch Chapte 14 in Antenna Theo, Analsis and Design (4th dition) b Balanis Cavit model Micostip antennas esemble dielectic-loaded cavities
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More information2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum
2. Electostatics D. Rakhesh Singh Kshetimayum 1 2.1 Intoduction In this chapte, we will study how to find the electostatic fields fo vaious cases? fo symmetic known chage distibution fo un-symmetic known
More information( ) ( )( ) ˆ. Homework #8. Chapter 27 Magnetic Fields II.
Homewok #8. hapte 7 Magnetic ields. 6 Eplain how ou would modif Gauss s law if scientists discoveed that single, isolated magnetic poles actuall eisted. Detemine the oncept Gauss law fo magnetism now eads
More information18.06 Problem Set 4 Solution
8.6 Poblem Set 4 Solution Total: points Section 3.5. Poblem 2: (Recommended) Find the lagest possible numbe of independent vectos among ) ) ) v = v 4 = v 5 = v 6 = v 2 = v 3 =. Solution (4 points): Since
More informationPhys101 Lectures 30, 31. Wave Motion
Phys0 Lectues 30, 3 Wave Motion Key points: Types of Waves: Tansvese and Longitudinal Mathematical Repesentation of a Taveling Wave The Pinciple of Supeposition Standing Waves; Resonance Ref: -7,8,9,0,,6,,3,6.
More informationdq 1 (5) q 1 where the previously mentioned limit has been taken.
1 Vecto Calculus And Continuum Consevation Equations In Cuvilinea Othogonal Coodinates Robet Maska: Novembe 25, 2008 In ode to ewite the consevation equations(continuit, momentum, eneg) to some cuvilinea
More informationAppendix B The Relativistic Transformation of Forces
Appendix B The Relativistic Tansfomation of oces B. The ou-foce We intoduced the idea of foces in Chapte 3 whee we saw that the change in the fou-momentum pe unit time is given by the expession d d w x
More informationMULTILAYER PERCEPTRONS
Last updated: Nov 26, 2012 MULTILAYER PERCEPTRONS Outline 2 Combining Linea Classifies Leaning Paametes Outline 3 Combining Linea Classifies Leaning Paametes Implementing Logical Relations 4 AND and OR
More informationQuantum theory of angular momentum
Quantum theoy of angula momentum Igo Mazets igo.mazets+e141@tuwien.ac.at (Atominstitut TU Wien, Stadionallee 2, 1020 Wien Time: Fiday, 13:00 14:30 Place: Feihaus, Sem.R. DA gün 06B (exception date 18 Nov.:
More information2. Plane Elasticity Problems
S0 Solid Mechanics Fall 009. Plane lasticity Poblems Main Refeence: Theoy of lasticity by S.P. Timoshenko and J.N. Goodie McGaw-Hill New Yok. Chaptes 3..1 The plane-stess poblem A thin sheet of an isotopic
More informationFlux. Area Vector. Flux of Electric Field. Gauss s Law
Gauss s Law Flux Flux in Physics is used to two distinct ways. The fist meaning is the ate of flow, such as the amount of wate flowing in a ive, i.e. volume pe unit aea pe unit time. O, fo light, it is
More informationPhysics 2A Chapter 10 - Moment of Inertia Fall 2018
Physics Chapte 0 - oment of netia Fall 08 The moment of inetia of a otating object is a measue of its otational inetia in the same way that the mass of an object is a measue of its inetia fo linea motion.
More informationMATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE. We consider second order constant coefficient scalar linear PDEs on R n. These have the form
MATH 220: SECOND ORDER CONSTANT COEFFICIENT PDE ANDRAS VASY We conside second ode constant coefficient scala linea PDEs on R n. These have the fom Lu = f L = a ij xi xj + b i xi + c i whee a ij b i and
More informationB da = 0. Q E da = ε. E da = E dv
lectomagnetic Theo Pof Ruiz, UNC Asheville, doctophs on YouTube Chapte Notes The Maxwell quations in Diffeential Fom 1 The Maxwell quations in Diffeential Fom We will now tansfom the integal fom of the
More information3.6 Applied Optimization
.6 Applied Optimization Section.6 Notes Page In this section we will be looking at wod poblems whee it asks us to maimize o minimize something. Fo all the poblems in this section you will be taking the
More informationSolution to Problem First, the firm minimizes the cost of the inputs: min wl + rk + sf
Econ 0A Poblem Set 4 Solutions ue in class on Tu 4 Novembe. No late Poblem Sets accepted, so! This Poblem set tests the knoledge that ou accumulated mainl in lectues 5 to 9. Some of the mateial ill onl
More informationTheWaveandHelmholtzEquations
TheWaveandHelmholtzEquations Ramani Duaiswami The Univesity of Mayland, College Pak Febuay 3, 2006 Abstact CMSC828D notes (adapted fom mateial witten with Nail Gumeov). Wok in pogess 1 Acoustic Waves 1.1
More informationPHYS 1444 Section 501 Lecture #7
PHYS 1444 Section 51 Lectue #7 Wednesday, Feb. 8, 26 Equi-potential Lines and Sufaces Electic Potential Due to Electic Dipole E detemined fom V Electostatic Potential Enegy of a System of Chages Capacitos
More informationME 210 Applied Mathematics for Mechanical Engineers
Tangent and Ac Length of a Cuve The tangent to a cuve C at a point A on it is defined as the limiting position of the staight line L though A and B, as B appoaches A along the cuve as illustated in the
More informationENGI 4430 Non-Cartesian Coordinates Page xi Fy j Fzk from Cartesian coordinates z to another orthonormal coordinate system u, v, ˆ i ˆ ˆi
ENGI 44 Non-Catesian Coodinates Page 7-7. Conesions between Coodinate Systems In geneal, the conesion of a ecto F F xi Fy j Fzk fom Catesian coodinates x, y, z to anothe othonomal coodinate system u,,
More informationMAC Module 12 Eigenvalues and Eigenvectors
MAC 23 Module 2 Eigenvalues and Eigenvectos Leaning Objectives Upon completing this module, you should be able to:. Solve the eigenvalue poblem by finding the eigenvalues and the coesponding eigenvectos
More informationAppendix A. Appendices. A.1 ɛ ijk and cross products. Vector Operations: δ ij and ɛ ijk
Appendix A Appendices A1 ɛ and coss poducts A11 Vecto Opeations: δ ij and ɛ These ae some notes on the use of the antisymmetic symbol ɛ fo expessing coss poducts This is an extemely poweful tool fo manipulating
More informationScientific Computing II
Scientific Computing II Conjugate Gadient Methods Michael Bade Summe 2014 Conjugate Gadient Methods, Summe 2014 1 Families of Iteative Solves elaxation methods: Jacobi-, Gauss-Seidel-Relaxation,... Ove-Relaxation-Methods
More informationWhen a mass moves because of a force, we can define several types of problem.
Mechanics Lectue 4 3D Foces, gadient opeato, momentum 3D Foces When a mass moves because of a foce, we can define seveal types of poblem. ) When we know the foce F as a function of time t, F=F(t). ) When
More informationPhotonic Crystals: Periodic Surprises in Electromagnetism. Steven G. Johnson MIT
Photonic Cystals: Peiodic Supises in Electomagnetism Steven G. Johnson MIT To Begin: A Catoon in 2d k scatteing planewave E, H ~ e i( k x -wt ) k = w / c = 2p l To Begin: A Catoon in 2d k a planewave E,
More informationMODULE 5a and 5b (Stewart, Sections 12.2, 12.3) INTRO: In MATH 1114 vectors were written either as rows (a1, a2,..., an) or as columns a 1 a. ...
MODULE 5a and 5b (Stewat, Sections 2.2, 2.3) INTRO: In MATH 4 vectos wee witten eithe as ows (a, a2,..., an) o as columns a a 2... a n and the set of all such vectos of fixed length n was called the vecto
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More informationReview: Electrostatics and Magnetostatics
Review: Electostatics and Magnetostatics In the static egime, electomagnetic quantities do not vay as a function of time. We have two main cases: ELECTROSTATICS The electic chages do not change postion
More informationQ. Obtain the Hamiltonian for a one electron atom in the presence of an external magnetic field.
Syed Ashad Hussain Lectue Deatment of Physics Tiua Univesity www.sahussaintu.wodess.com Q. Obtain the Hamiltonian fo a one electon atom in the esence of an extenal magnetic field. To have an idea about
More information7.2.1 Basic relations for Torsion of Circular Members
Section 7. 7. osion In this section, the geomety to be consideed is that of a long slende cicula ba and the load is one which twists the ba. Such poblems ae impotant in the analysis of twisting components,
More informationMCV4U Final Exam Review. 1. Consider the function f (x) Find: f) lim. a) lim. c) lim. d) lim. 3. Consider the function: 4. Evaluate. lim. 5. Evaluate.
MCVU Final Eam Review Answe (o Solution) Pactice Questions Conside the function f () defined b the following gaph Find a) f ( ) c) f ( ) f ( ) d) f ( ) Evaluate the following its a) ( ) c) sin d) π / π
More information( ) ( ) Review of Force. Review of Force. r = =... Example 1. What is the dot product for F r. Solution: Example 2 ( )
: PHYS 55 (Pat, Topic ) Eample Solutions p. Review of Foce Eample ( ) ( ) What is the dot poduct fo F =,,3 and G = 4,5,6? F G = F G + F G + F G = 4 +... = 3 z z Phs55 -: Foce Fields Review of Foce Eample
More informationElectromagnetic Waves
Chapte 32 Electomagnetic Waves PowePoint Lectues fo Univesity Physics, Twelfth Edition Hugh D. Young and Roge A. Feedman Lectues by James Pazun Modified P. Lam 8_11_2008 Topics fo Chapte 32 Maxwell s equations
More informationB. Spherical Wave Propagation
11/8/007 Spheical Wave Popagation notes 1/1 B. Spheical Wave Popagation Evey antenna launches a spheical wave, thus its powe density educes as a function of 1, whee is the distance fom the antenna. We
More informationObjects usually are charged up through the transfer of electrons from one object to the other.
1 Pat 1: Electic Foce 1.1: Review of Vectos Review you vectos! You should know how to convet fom pola fom to component fom and vice vesa add and subtact vectos multiply vectos by scalas Find the esultant
More informationCOORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT
COORDINATE TRANSFORMATIONS - THE JACOBIAN DETERMINANT Link to: phsicspages home page. To leave a comment o epot an eo, please use the auilia blog. Refeence: d Inveno, Ra, Intoducing Einstein s Relativit
More informationMath 2263 Solutions for Spring 2003 Final Exam
Math 6 Solutions fo Sping Final Exam ) A staightfowad appoach to finding the tangent plane to a suface at a point ( x, y, z ) would be to expess the cuve as an explicit function z = f ( x, y ), calculate
More information. Using our polar coordinate conversions, we could write a
504 Chapte 8 Section 8.4.5 Dot Poduct Now that we can add, sutact, and scale vectos, you might e wondeing whethe we can multiply vectos. It tuns out thee ae two diffeent ways to multiply vectos, one which
More informationCHAPTER 25 ELECTRIC POTENTIAL
CHPTE 5 ELECTIC POTENTIL Potential Diffeence and Electic Potential Conside a chaged paticle of chage in a egion of an electic field E. This filed exets an electic foce on the paticle given by F=E. When
More informationGraphs of Sine and Cosine Functions
Gaphs of Sine and Cosine Functions In pevious sections, we defined the tigonometic o cicula functions in tems of the movement of a point aound the cicumfeence of a unit cicle, o the angle fomed by the
More informationFields and Waves I Spring 2005 Homework 4. Due 8 March 2005
Homewok 4 Due 8 Mach 005. Inceasing the Beakdown Voltage: This fist question is a mini design poject. You fist step is to find a commecial cable (coaxial o two wie line) fo which you have the following
More informationSpring 2001 Physics 2048 Test 3 solutions
Sping 001 Physics 048 Test 3 solutions Poblem 1. (Shot Answe: 15 points) a. 1 b. 3 c. 4* d. 9 e. 8 f. 9 *emembe that since KE = ½ mv, KE must be positive Poblem (Estimation Poblem: 15 points) Use momentum-impulse
More informationA Crash Course in (2 2) Matrices
A Cash Couse in ( ) Matices Seveal weeks woth of matix algeba in an hou (Relax, we will only stuy the simplest case, that of matices) Review topics: What is a matix (pl matices)? A matix is a ectangula
More informationChapter 1. Introduction
Chapte 1 Intoduction 1.1 The Natue of Phsics Phsics has developed out of the effots of men and women to eplain ou phsical envionment. Phsics encompasses a emakable vaiet of phenomena: planeta obits adio
More informationRight-handed screw dislocation in an isotropic solid
Dislocation Mechanics Elastic Popeties of Isolated Dislocations Ou study of dislocations to this point has focused on thei geomety and thei ole in accommodating plastic defomation though thei motion. We
More informationPlane Wave Expansion Method (PWEM)
/15/18 Instucto D. Rymond Rumpf (915) 747 6958 cumpf@utep.edu EE 5337 Computtionl Electomgnetics Lectue #19 Plne Wve Expnsion Method (PWEM) Lectue 19 These notes my contin copyighted mteil obtined unde
More informationDOING PHYSICS WITH MATLAB COMPUTATIONAL OPTICS
DOING PHYIC WITH MTLB COMPUTTIONL OPTIC FOUNDTION OF CLR DIFFRCTION THEORY Ian Coope chool of Physics, Univesity of ydney ian.coope@sydney.edu.au DOWNLOD DIRECTORY FOR MTLB CRIPT View document: Numeical
More informationUnit 7: Sources of magnetic field
Unit 7: Souces of magnetic field Oested s expeiment. iot and Savat s law. Magnetic field ceated by a cicula loop Ampèe s law (A.L.). Applications of A.L. Magnetic field ceated by a: Staight cuent-caying
More informationPhysics 2212 GH Quiz #2 Solutions Spring 2016
Physics 2212 GH Quiz #2 Solutions Sping 216 I. 17 points) Thee point chages, each caying a chage Q = +6. nc, ae placed on an equilateal tiangle of side length = 3. mm. An additional point chage, caying
More informationPhysics for Scientists and Engineers
Phsics 111 Sections 003 and 005 Instucto: Pof. Haimin Wang E-mail: haimin@flae.njit.edu Phone: 973-596-5781 Office: 460 Tienan Hall Homepage: http://sola.njit.edu/~haimin Office Hou: 2:30 to 3:50 Monda
More informationGeneral Solution of EM Wave Propagation in Anisotropic Media
Jounal of the Koean Physical Society, Vol. 57, No. 1, July 2010, pp. 55 60 Geneal Solution of EM Wave Popagation in Anisotopic Media Jinyoung Lee Electical and Electonic Engineeing Depatment, Koea Advanced
More information2 Lecture 2: The Bohr atom (1913) and the Schrödinger equation (1925)
1 Lectue 1: The beginnings of quantum physics 1. The Sten-Gelach expeiment. Atomic clocks 3. Planck 1900, blackbody adiation, and E ω 4. Photoelectic effect 5. Electon diffaction though cystals, de Boglie
More information6.4 Period and Frequency for Uniform Circular Motion
6.4 Peiod and Fequency fo Unifom Cicula Motion If the object is constained to move in a cicle and the total tangential foce acting on the total object is zeo, F θ = 0, then (Newton s Second Law), the tangential
More informationCS-184: Computer Graphics. Today
CS-184: Compute Gaphics Lectue #6: 3D Tansfomations and Rotations Pof. James O Bien Univesity of Califonia, Bekeley V2006-F-06-1.0 Today Tansfomations in 3D Rotations Matices Eule angles Eponential maps
More informationMultifrontal sparse QR factorization on the GPU
Multifontal spase QR factoization on the GPU Tim Davis, Sanjay Ranka, Shaanyan Chetlu, Nui Yealan Univesity of Floida Feb 2012 GPU-based Multifontal QR factoization why spase QR? multifontal spase QR in
More informationQualifying Examination Electricity and Magnetism Solutions January 12, 2006
1 Qualifying Examination Electicity and Magnetism Solutions Januay 12, 2006 PROBLEM EA. a. Fist, we conside a unit length of cylinde to find the elationship between the total chage pe unit length λ and
More informationPhysics 161 Fall 2011 Extra Credit 2 Investigating Black Holes - Solutions The Following is Worth 50 Points!!!
Physics 161 Fall 011 Exta Cedit Investigating Black Holes - olutions The Following is Woth 50 Points!!! This exta cedit assignment will investigate vaious popeties of black holes that we didn t have time
More informationQuantum Fourier Transform
Chapte 5 Quantum Fouie Tansfom Many poblems in physics and mathematics ae solved by tansfoming a poblem into some othe poblem with a known solution. Some notable examples ae Laplace tansfom, Legende tansfom,
More informationFresnel Diffraction. monchromatic light source
Fesnel Diffaction Equipment Helium-Neon lase (632.8 nm) on 2 axis tanslation stage, Concave lens (focal length 3.80 cm) mounted on slide holde, iis mounted on slide holde, m optical bench, micoscope slide
More informationPhysics Fall Mechanics, Thermodynamics, Waves, Fluids. Lecture 18: System of Particles II. Slide 18-1
Physics 1501 Fall 2008 Mechanics, Themodynamics, Waves, Fluids Lectue 18: System of Paticles II Slide 18-1 Recap: cente of mass The cente of mass of a composite object o system of paticles is the point
More informationCS 475 / CS 675 Computer Graphics
Pojectos CS 475 / CS 675 Compute Gaphics Lectue 6 : Viewing Cente of Pojection Object Image Plane o Pojection Plane Othogaphic Pojection Pojectos Cente of Pojection? Object Tue sie o shape fo lines Fo
More informationOn the Quasi-inverse of a Non-square Matrix: An Infinite Solution
Applied Mathematical Sciences, Vol 11, 2017, no 27, 1337-1351 HIKARI Ltd, wwwm-hikaicom https://doiog/1012988/ams20177273 On the Quasi-invese of a Non-squae Matix: An Infinite Solution Ruben D Codeo J
More informationHopefully Helpful Hints for Gauss s Law
Hopefully Helpful Hints fo Gauss s Law As befoe, thee ae things you need to know about Gauss s Law. In no paticula ode, they ae: a.) In the context of Gauss s Law, at a diffeential level, the electic flux
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More informationTo Feel a Force Chapter 7 Static equilibrium - torque and friction
To eel a oce Chapte 7 Chapte 7: Static fiction, toque and static equilibium A. Review of foce vectos Between the eath and a small mass, gavitational foces of equal magnitude and opposite diection act on
More informationAn Exact Solution of Navier Stokes Equation
An Exact Solution of Navie Stokes Equation A. Salih Depatment of Aeospace Engineeing Indian Institute of Space Science and Technology, Thiuvananthapuam, Keala, India. July 20 The pincipal difficulty in
More informationNon-Linear Dynamics Homework Solutions Week 2
Non-Linea Dynamics Homewok Solutions Week Chis Small Mach, 7 Please email me at smach9@evegeen.edu with any questions o concens eguading these solutions. Fo the ececises fom section., we sketch all qualitatively
More information