FDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
|
|
- Barrie Dorsey
- 5 years ago
- Views:
Transcription
1 FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c 2 2 ( t) ( x) t n +1 n +1 n 1 i + 1 i i 1 x 2 n+ 1 c t n n 2 n n 1 i i i i i 2 ( + 1 1) 21 ( γ ) H = H + H + H H x γ 26
2 FDTD for 1D wave equation (cont d) Properties o Fully explicit second-order accurate scheme o Two previous time steps are required to be stored o The analysis is recursive in time o Each field is found at every time step at every node t n +1 n +1 n 1 i + 1 i i 1 x 27
3 FDTD for 1D Maxwell s equations Equation: Notations: o Redefine: 2ΔΔ ΔΔ o Central difference: H 1 E E 1 H = ; = t µ x t ε x H Ht ( + t2) Ht ( t2) t t o Involved time steps: t = nδδ, t = nδδ + ΔΔ/2 o Involved spatial steps: x = iδx, x = iδx + Δx/2 o Involved fields: Hn ( ti, x+ x2) = H En ( t+ t 2, i x) = E n i+ 12 n+ 12 i t n +1 n +1 2 n i 12 i i + 12 x H E H 28
4 FDTD for 1D Maxwell s equations (cont d) Update equations: H = H t 1 + E x µ E E = E t 1 + H x ε H ( ) n+ 1 n n+ 12 n+ 12 i+ 12 i+ 12 i+ 1 i ( ) n+ 12 n 12 n n i i i i n +1 n +1 2 n Yee algorithm: o The fields are given in different space and time locations o Interleaved E & H grids o Leapfrog in time o Explicit stepping no matrix inversions, no iterations o The solution is more robust and straightforward o E & H boundary conditions are enforced automatically t i 12 i i + 12 x H E H 29
5 Numerical dispersion of FDTD Consider a time harmonic wave H = Continuous dispersion relation is linear: n j( n t ki x) Numerical wave: Hi = e ω Use the FDTD update equation γ ( ) 21 ( γ ) e = γ e + e + e 2 ω t 2 2 k x sin = γ sin 2 2 jω t 2 jk x jk x 2 jω t j( t kx) e ω ω = ck Properties: o In general the phase velocity depends on the frequency numerical dispersion! o ΔΔ, Δx 0 ω 2 = cc 2 - weak dispersion o Δx = cδt, γ = 1 magic time step with no dispersion (for the 1D case) ( + 1 1) 21 ( γ ) H = H + H + H n 1 2 n n 2 n 1 i i i i 30
6 Stability of FDTD Consider H = e j(ωω kk) and H i n = e j(ωωωω kkkk) Allow ω to be complex, i.e. ω = ω + jω The propagation is stable if ω 0 (wave~e ω t ) For ω < 0 (wave~e ω t ) the propagation is unstable Define jω Ω= ω t 2, K = k x 2, g = e and use sin Ω γ sin K = g + 2 jγg sin K 1 = 0 g = jk sin K ± 1 γ sin K Stable solution: 2 2 g 1 1 γ sin K 0 γ 1 Courant stability criterion for 1D case: t x c 31
7 Stability of FDTD (cont d) Physical meaning: the wave should not pass through more than one cell in a single time step Under the magic time step the scheme is exact How to choose ΔΔ and Δt If no fine geometrical features are present: ΔΔ < λ min 10 to satisfy Nyquist sampling criterion If fine geometrical features are present: ΔΔ is determined by the features Choose Δt < ΔΔ/c (note that for fine geometrical features, ΔΔ is small and the time step is small as well!) 32
8 Source modeling Initial conditions: Define E i 0, H i 1 for all i. Hard source o Specify fields at certain locations independent of surrounding cells, i.e. replace the Yee equations, e.g. n = F(t n E z,is o Easy to implement but lead to reflections Soft source o Sources that allow field updates, e.g. current sources o Equations are updated and no reflections but harder to implement Plane wave injection: Scattered field / total field approach 33
9 Time dependence of sources Time harmonic sources o Simple sin functions leads to discontinuities and a relatively broadband spectrum o Use sin with a smooth transition function (g(t)) Ft () = Fsin( ω n tgt ) () 0 0 Pulsed sources o Gaussian: Contains DC a problem o Time derivative of a Gaussian no DC o Modulated Gaussian 34
10 Domain truncation: Perfectly matched layer (PML) Add an additional domain where a truncation is needed Fill the domain with a matching lossy material with σ & σ e satisfying σ m σ µ 1 µ 2 = & = µ 1 ε1 ε1 ε2 No reflections are obtained because computational domain PML domain µ µ 1 + σ ( jωµ ) η = = η = 1 ( ) 1 2 m ε1 ε 2 + σ jωε 2 This idea can be extended to 2D and 3D but the lossy materials should be anisotropic to match all directions Other absorbing boundary conditions can be used as well 35
11 36
12 37
13 38
14 39
15 40
16 41
17 Outline Introduction Maxwell s equations Finite Difference Time Domain method Finite element method Integral equation methods 42
18 Integral equation procedure Formulate an integral equation Expand the unknowns in terms of basis functions Define a discrete matrix representation of the continuous integral equation Solve the matrix equation Post-process to derive the required parameters 43
19 Advantages Only the structure of interest is discretized Outgoing conditions are satisfied automatically No numerical dispersion Can be formulated in time and frequency Can be defined for surfaces and volumes Can be formulated for complex backgrounds Can be higher order accurate 44
20 Challenges Computational cost of simple methods is ~O(N 2 ) Need for fast methods (~O(N) or O(NlogN)) for N degrees of freedom Different formulations are required for different problem types, e.g. volumes vs. surfaces 45
21 Equation of electrostatics Electrostatics t = ω = 0 two uncoupled sets of (static) equations are obtained for E& D and H& B Equation of electrostatics E = 0 D = ρ Boundary conditions nˆ ( E E ) = or nˆ 12 ( D2 D1) = ρs Constitutive relation D S C = εe E dl = D ds= 0 Q 46
22 Electric scalar potential Potential due to various charge distributions Point charge Many discrete charges Various continuous charge distributions in free space 47
23 Capacitance Consider a capacitor consisting of two conductors with a voltage Opposite charges appear on the conductors Capacitance Alternative expression Relation to resistance 48
24 Numerical capacitance extraction (1) Integral equation method Given V = V 0 2 Two conductors Voltage on the conductors is Find the capacitance Formulation C V 0 2and V 0 2 Potential on the surfaces is given by V 0 2 and V Potential everywhere is calculated as V () r = ρ s ( r ) ds 4 πε 0 r r S surface charge where ρ s is an unknown surface charge distribution Green's function Equate the potentials on the surfaces S (1) (2) + S 1 V0 2; r S1 ρ ( ) 4 s r ds = πε 0 r r V 2; S surface charge r Green's function distribution (unknown) 0 2 V = V 0 2 (1) surface S (2) surface S integral equation 49
25 Numerical capacitance extraction (2) Matrix equation V = V 0 2 Discretize the problem Divide the surfaces into N = N + N S small patches, N (1), N (2) are numbers S S (1) (2) of patches on surfaces S, S Replace the integration by summation Enforce the IE at the location of charges N ρ ( ) 0 2; s r ds V rm S1 = 4 πε r r V 2; r S (1) (2) n= 1 S 0 m 0 m 2 n S V = V 0 2 patch S n (1) surface S of area s (2) surface S n Obtain a matrix equation N 0 2; N ρs ( r ) ds V r 0 m S1 = ZmnQn = Vm ZQ = V 4 πε r r V 2; r S n= 1 S 0 m 0 m 2 n= 1 n matrix equation 50
26 Numerical capacitance extraction (3) Matrix equation and solution Matrix equation Z Z : N N matrix; Z = Q: N vector; Q ( ) s V : N vector; V Solution ZQ mn m n 1 ds 4 πε 0 rm rn mn s 4 1 n πε 0 r Sn m r Zmn ; rm = rn 2ε0 π s n = ρs rn n n m V0 2; rm S1 = V 2; r S Q 1 Z V Extracted capacitance = = V 0 m 2 total charge on surface 1 C = Q (1) S V C = 1 N S (1) n= 1 V 0 0 Q n ; r r 51
27 52
28 53
29 Outline Introduction Maxwell s equations Finite Difference Time Domain method Finite element method Integral equation methods 54
30 Basic FEM Steps 1. Discretization/subdivision of solution domain 1-D: 2-D: 3-D: 2. Selection of interpolation schemes Linear or higher-order polynomials 3. Formulation of the system of equations Using either the Ritz or Galerkin method: Formulate elemental equations and assemble 4. Solution of the system of equations Using either a direct or iterative method 55
31 Advantages Works for general structures Can be formulated in time and frequency Can be higher order accurate Low cost (~O(N)) 56
32 Boundary-Value Problem Ω ˆn Γ 57
33 Equivalent Variational Problem 58
34 Equivalent Variational Problem 59
35 Equivalent Variational Problem 60
36 Equivalent Variational Problem 61
37 Basic FEM Steps 1. Discretization/subdivision of solution domain 1-D: 2-D: 3-D: 2. Selection of interpolation schemes Linear or higher-order polynomials 3. Formulation of the system of equations Using either the Ritz or Galerkin method: Formulate elemental equations and assemble 4. Solution of the system of equations Using either a direct or iterative method 62
38 FEM Analysis Domain subdivision Step 1: Domain Discretization 63
39 FEM Analysis Domain subdivision Step 1: Domain Discretization 64
40 FEM Analysis Domain subdivision Step 1: Domain Discretization 65
41 FEM Analysis Domain subdivision Step 1: Domain Discretization 66
42 FEM Analysis Element interpolation Step 2: Element Interpolation 3 1 e 2 67
43 FEM Analysis Element interpolation 68
44 FEM Analysis Element formulation Step 3: Formulation of the System of Equations A. Elemental equations 69
45 FEM Analysis Element formulation Elemental functional: 70
46 FEM Analysis Element formulation Integration formula: 71
47 FEM Analysis Element formulation Use matrix notation: 72
48 FEM Analysis Assembly B. Assembly 73
49 FEM Analysis Assembly Carry out the summation: Apply the stationarity condition: 74
50 FEM Analysis Assembly How to carry out the summation? 75
51 FEM Analysis Assembly Example:
52 FEM Analysis 1. Start from a null matrix and add in the first element: 77
53 FEM Analysis 2. Add in the second element: 78
54 FEM Analysis 3. Add in the third element: 79
55 FEM Analysis 4. Add in the fourth element: 80
56 FEM Analysis 5. Follow a similar procedure: 81
57 FEM Analysis Apply BC C. Impose the Dirichlet Boundary Condition: Approach #1: To impose, simply set: To maintain symmetry, set: 82
58 FEM Analysis Apply BC After imposing,, : Remains Symmetric! 83
59 FEM Analysis Apply BC Can be made smaller: Worthwhile when there are many prescribed nodes. 84
60 FEM Analysis Apply BC Approach #2 (Simple one): To impose, simply set: After imposing,, : 85
Pseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationECE 107: Electromagnetism
ECE 107: Electromagnetism Set 7: Dynamic fields Instructor: Prof. Vitaliy Lomakin Department of Electrical and Computer Engineering University of California, San Diego, CA 92093 1 Maxwell s equations Maxwell
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More informationA Time Domain Approach to Power Integrity for Printed Circuit Boards
A Time Domain Approach to Power Integrity for Printed Circuit Boards N. L. Mattey 1*, G. Edwards 2 and R. J. Hood 2 1 Electrical & Optical Systems Research Division, Faculty of Engineering, University
More informationChapter 4. Two-Dimensional Finite Element Analysis
Chapter 4. Two-Dimensional Finite Element Analysis general two-dimensional boundary-value problems 4.1 The Boundary-Value Problem 2nd-order differential equation to consider α α β φ Ω (4.1) Laplace, Poisson
More informationFDFD. The Finite-Difference Frequency-Domain Method. Hans-Dieter Lang
FDFD The Finite-Difference Frequency-Domain Method Hans-Dieter Lang Friday, December 14, 212 ECE 1252 Computational Electrodynamics Course Project Presentation University of Toronto H.-D. Lang FDFD 1/18
More informationAntennas and Propagation. Chapter 2: Basic Electromagnetic Analysis
Antennas and Propagation : Basic Electromagnetic Analysis Outline Vector Potentials, Wave Equation Far-field Radiation Duality/Reciprocity Transmission Lines Antennas and Propagation Slide 2 Antenna Theory
More informationDivergent Fields, Charge, and Capacitance in FDTD Simulations
Divergent Fields, Charge, and Capacitance in FDTD Simulations Christopher L. Wagner and John B. Schneider August 2, 1998 Abstract Finite-difference time-domain (FDTD) grids are often described as being
More informationField and Wave Electromagnetic
Field and Wave Electromagnetic Chapter7 The time varying fields and Maxwell s equation Introduction () Time static fields ) Electrostatic E =, id= ρ, D= εe ) Magnetostatic ib=, H = J, H = B μ note) E and
More informationA Nested Dissection Parallel Direct Solver. for Simulations of 3D DC/AC Resistivity. Measurements. Maciej Paszyński (1,2)
A Nested Dissection Parallel Direct Solver for Simulations of 3D DC/AC Resistivity Measurements Maciej Paszyński (1,2) David Pardo (2), Carlos Torres-Verdín (2) (1) Department of Computer Science, AGH
More informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationSimulation and Numerical Modeling of a Rectangular Patch Antenna Using Finite Difference Time Domain (FDTD) Method
Journal of Computer Science and Information Technology June 2014, Vol. 2, No. 2, pp. 01-08 ISSN: 2334-2366 (Print), 2334-2374 (Online) Copyright The Author(s). 2014. All Rights Reserved. Published by American
More informationTHE FINITE-DIFFERENCE TIME-DOMAIN (FDTD) METHOD PART IV
Numerical Techniques in Electromagnetics ECE 757 THE FINITE-DIFFERENCE TIME-DOMAIN (FDTD) METHOD PART IV The Perfectly Matched Layer (PML) Absorbing Boundary Condition Nikolova 2009 1 1. The need for good
More informationECE 497 JS Lecture - 13 Projects
ECE 497 JS Lecture - 13 Projects Spring 2004 Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jose@emlab.uiuc.edu 1 ECE 497 JS - Projects All projects should be accompanied
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationHigh-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials
Progress In Electromagnetics Research Letters, Vol. 77, 103 107, 01 High-Order FDTD with Exponential Time Differencing Algorithm for Modeling Wave Propagation in Debye Dispersive Materials Wei-Jun Chen
More informationChapter 4 Reflection and Transmission of Waves
4-1 Chapter 4 Reflection and Transmission of Waves ECE 3317 Dr. Stuart Long www.bridgat.com www.ranamok.com Boundary Conditions 4- -The convention is that is the outward pointing normal at the boundary
More informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More information6.013 Recitation 11. Quasistatic Electric and Magnetic Fields in Devices and Circuit Elements
6.013 Recitation 11 Quasistatic Electric and Magnetic Fields in Devices and Circuit Elements A. Introduction The behavior of most electric devices depends on static or slowly varying (quasistatic 1 ) electric
More informationModel Order Reduction and Stability Enforcement of Finite-Difference Time-Domain Equations Beyond the CFL Limit. Xihao Li
Model Order Reduction and Stability Enforcement of Finite-Difference Time-Domain Equations Beyond the CFL Limit by Xihao Li A thesis submitted in conformity with the requirements for the degree of Master
More informationCOLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE
Progress In Electromagnetics Research M, Vol. 3, 239 252, 213 COLLOCATED SIBC-FDTD METHOD FOR COATED CONDUCTORS AT OBLIQUE INCIDENCE Lijuan Shi 1, 3, Lixia Yang 2, *, Hui Ma 2, and Jianning Ding 3 1 School
More informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
More informationINSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad Electronics and Communicaton Engineering
INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 00 04 Electronics and Communicaton Engineering Question Bank Course Name : Electromagnetic Theory and Transmission Lines (EMTL) Course Code :
More informationELE 3310 Tutorial 10. Maxwell s Equations & Plane Waves
ELE 3310 Tutorial 10 Mawell s Equations & Plane Waves Mawell s Equations Differential Form Integral Form Faraday s law Ampere s law Gauss s law No isolated magnetic charge E H D B B D J + ρ 0 C C E r dl
More informationComparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer
Comparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer V. A. Bokil a and M. W. Buksas b Center for Research in Scientific Computation a North Carolina
More informationChapter 1. Introduction to Electrostatics
Chapter. Introduction to Electrostatics. Electric charge, Coulomb s Law, and Electric field Electric charge Fundamental and characteristic property of the elementary particles There are two and only two
More informationUniform Plane Waves. Ranga Rodrigo. University of Moratuwa. November 7, 2008
Uniform Plane Waves Ranga Rodrigo University of Moratuwa November 7, 2008 Ranga Rodrigo (University of Moratuwa) Uniform Plane Waves November 7, 2008 1 / 51 Summary of Last Week s Lecture Basic Relations
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 informationLecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations
Poisson s and Laplace s Equations Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay We will spend some time in looking at the mathematical foundations of electrostatics.
More informationFOR most applications, the finite-difference time-domain
ECE 1252 INTRODUCTION TO COMPUTATIONAL ELECTRODYNAMICS 1 The Finite-Difference Frequency-Domain Method Hans-Dieter Lang, Student-Member, IEEE Abstract The finite-difference frequency-domain (FDFD) method
More informationGeneral review: - a) Dot Product
General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationECE 6340 Fall Homework 2. Please do the following problems (you may do the others for practice if you wish): Probs. 1, 2, 3, 4, 5, 6, 7, 10, 12
ECE 634 Fall 16 Homework Please do the following problems (you may do the others for practice if you wish: Probs. 1,, 3, 4, 5, 6, 7, 1, 1 1 Consider two parallel infinite wires in free space each carrying
More informationENGIN 211, Engineering Math. Fourier Series and Transform
ENGIN 11, Engineering Math Fourier Series and ransform 1 Periodic Functions and Harmonics f(t) Period: a a+ t Frequency: f = 1 Angular velocity (or angular frequency): ω = ππ = π Such a periodic function
More informationTECHNO INDIA BATANAGAR
TECHNO INDIA BATANAGAR ( DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING) QUESTION BANK- 2018 1.Vector Calculus Assistant Professor 9432183958.mukherjee@tib.edu.in 1. When the operator operates on
More informationFinite Element Method (FEM)
Finite Element Method (FEM) The finite element method (FEM) is the oldest numerical technique applied to engineering problems. FEM itself is not rigorous, but when combined with integral equation techniques
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 1
ECE 6341 Spring 16 Prof. David R. Jackson ECE Dept. Notes 1 1 Fields in a Source-Free Region Sources Source-free homogeneous region ( ε, µ ) ( EH, ) Note: For a lossy region, we replace ε ε c ( / ) εc
More informationModul 3. Finite-difference time-domain (FDTD)
Modul 3 Finite-difference time-domain (FDTD) based on Dennis Sullivan, A Brief Introduction to The Finite-Difference Time-Domain (FDTD) Method http://www.mrc.uidaho.edu/~dennis/ece538-files/intro(fdtd).doc
More informationEE243 Advanced Electromagnetic Theory Lec #3: Electrostatics (Apps., Form),
EE4 Advanced Electromagnetic Theory Lec #: Electrostatics Apps., Form, Electrostatic Boundary Conditions Energy, Force and Capacitance Electrostatic Boundary Conditions on Φ Image Solutions Eample Green
More informationPerfectly Matched Layer (PML) for Computational Electromagnetics
Perfectly Matched Layer (PML) for Computational Electromagnetics Copyright 2007 by Morgan & Claypool All rights reserved. No part of this publication may be reproduced, stored in a retrieval system, or
More informationProgress In Electromagnetics Research M, Vol. 13, 29 40, 2010
Progress In Electromagnetics Research M, Vol. 13, 9 40, 010 SHIFT-OPERATOR FINITE DIFFERENCE TIME DO- MAIN ANALYSIS OF CHIRAL MEDIUM A. M. Attiya Electrical Engineering Department King Saud University,
More informationBasics of Wave Propagation
Basics of Wave Propagation S. R. Zinka zinka@hyderabad.bits-pilani.ac.in Department of Electrical & Electronics Engineering BITS Pilani, Hyderbad Campus May 7, 2015 Outline 1 Time Harmonic Fields 2 Helmholtz
More informationDevelopment of a stable coupling of the Yee scheme with linear current
Development of a stable coupling of the Yee scheme with linear current Martin Campos Pinto (LJLL), Bruno Després (LJLL) Stéphane Heuraux (IJL), Filipe Da Silva (IPFN+IST) 15 octobre 2013 Munich: Pereverzev
More informationElectromagnetic Modeling and Simulation
Electromagnetic Modeling and Simulation Erin Bela and Erik Hortsch Department of Mathematics Research Experiences for Undergraduates April 7, 2011 Bela and Hortsch (OSU) EM REU 2010 1 / 45 Maxwell s Equations
More informationCHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L
CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit
More informationElectromagnetism. 1 ENGN6521 / ENGN4521: Embedded Wireless
Electromagnetism 1 ENGN6521 / ENGN4521: Embedded Wireless Radio Spectrum use for Communications 2 ENGN6521 / ENGN4521: Embedded Wireless 3 ENGN6521 / ENGN4521: Embedded Wireless Electromagnetism I Gauss
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More informationChapter 7. Time-Varying Fields and Maxwell s Equation
Chapter 7. Time-Varying Fields and Maxwell s Equation Electrostatic & Time Varying Fields Electrostatic fields E, D B, H =J D H 1 E B In the electrostatic model, electric field and magnetic fields are
More informationEnergy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media
Energy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media Yingda Cheng Michigan State University Computational Aspects of Time Dependent Electromagnetic Wave Problems
More informationA Single-Field Finite-Difference Time-Domain Formulations for Electromagnetic Simulations
Syracuse University SURFACE Electrical Engineering and Computer Science - Dissertations College of Engineering and Computer Science 2011 A Single-Field Finite-Difference Time-Domain Formulations for Electromagnetic
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationPlane Waves and Planar Boundaries in FDTD Simulations
Plane Waves and Planar Boundaries in FDTD Simulations School of Electrical Engineering and Computer Science John B. Schneider Elec. Eng. & Comp. Sci. Washington State Univ. Pullman, WA Robert J. Kruhlak
More informationWorked Examples Set 2
Worked Examples Set 2 Q.1. Application of Maxwell s eqns. [Griffiths Problem 7.42] In a perfect conductor the conductivity σ is infinite, so from Ohm s law J = σe, E = 0. Any net charge must be on the
More informationECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 7
ECE 634 Intermediate EM Waves Fall 16 Prof. David R. Jackson Dept. of ECE Notes 7 1 TEM Transmission Line conductors 4 parameters C capacitance/length [F/m] L inductance/length [H/m] R resistance/length
More informationEECS 117 Lecture 19: Faraday s Law and Maxwell s Eq.
University of California, Berkeley EECS 117 Lecture 19 p. 1/2 EECS 117 Lecture 19: Faraday s Law and Maxwell s Eq. Prof. Niknejad University of California, Berkeley University of California, Berkeley EECS
More informationAdvanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III
Advanced Engineering Electromagnetics, ECE750 LECTURE 11 THE FDTD METHOD PART III 1 11. Yee s discrete algorithm Maxwell s equations are discretized using central FDs. We set the magnetic loss equal to
More informationEfficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques
1 Efficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques From Resonators to Open Domains Rob Remis Delft University of Technology November 3, 2017 ICERM Brown University 2 Thanks
More informationUSAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION
October 4-6, 2007 - Chiinu, Rep.Moldova USAGE OF NUMERICAL METHODS FOR ELECTROMAGNETIC SHIELDS OPTIMIZATION Ionu- P. NICA, Valeriu Gh. DAVID, /tefan URSACHE Gh. Asachi Technical University Iai, Faculty
More informationSimulation of Electromagnetic Fields: The Finite-Difference Time-Domain (FDTD) Method and Its Applications
Simulation of Electromagnetic Fields: The Finite-Difference Time-Domain (FDTD) Method and Its Applications Veysel Demir, Ph.D. demir@ceet.niu.edu Department of Electrical Engineering, Northern Illinois
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
More informationChapter 7. Time-Varying Fields and Maxwell s Equations
Chapter 7. Time-arying Fields and Maxwell s Equations Electrostatic & Time arying Fields Electrostatic fields E, D B, H =J D H 1 E B In the electrostatic model, electric field and magnetic fields are not
More information3.3 Capacitance, relative permittivity & dielectrics 4
3.3 Capacitance, relative permittivity & dielectrics 4 +Q d E Gaussian surface Voltage, V Q Fig. 3.2. Parallel plate capacitor with the plates separated by a distance d which have been charged by a power
More informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
More informationThe Cooper Union Department of Electrical Engineering ECE135 Engineering Electromagnetics Exam II April 12, Z T E = η/ cos θ, Z T M = η cos θ
The Cooper Union Department of Electrical Engineering ECE135 Engineering Electromagnetics Exam II April 12, 2012 Time: 2 hours. Closed book, closed notes. Calculator provided. For oblique incidence of
More informationELE3310: Basic ElectroMagnetic Theory
A summary for the final examination EE Department The Chinese University of Hong Kong November 2008 Outline Mathematics 1 Mathematics Vectors and products Differential operators Integrals 2 Integral expressions
More informationPlane Waves GATE Problems (Part I)
Plane Waves GATE Problems (Part I). A plane electromagnetic wave traveling along the + z direction, has its electric field given by E x = cos(ωt) and E y = cos(ω + 90 0 ) the wave is (a) linearly polarized
More informationIEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 3, JUNE
IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 28, NO. 3, JUNE 2000 841 CTLSS An Advanced Electromagnetic Simulation Tool for Designing High-Power Microwave Sources Simon J. Cooke, Member, IEEE, Alfred A. Mondelli,
More informationPROCEEDINGS OF SPIE. FDTD method and models in optical education. Xiaogang Lin, Nan Wan, Lingdong Weng, Hao Zhu, Jihe Du
PROCEEDINGS OF SPIE SPIEDigitalLibrary.org/conference-proceedings-of-spie FDTD method and models in optical education Xiaogang Lin, Nan Wan, Lingdong Weng, Hao Zhu, Jihe Du Xiaogang Lin, Nan Wan, Lingdong
More informationElectromagnetic Field Theory (EMT)
Electromagnetic Field Theory (EMT) Lecture # 9 1) Coulomb s Law and Field Intensity 2) Electric Fields Due to Continuous Charge Distributions Line Charge Surface Charge Volume Charge Coulomb's Law Coulomb's
More information1 Electromagnetic concepts useful for radar applications
Electromagnetic concepts useful for radar applications The scattering of electromagnetic waves by precipitation particles and their propagation through precipitation media are of fundamental importance
More informationInterpolation and Deformations A short cookbook
Interpolation and Deformations A short cookbook 600.445 Fall 2000; Updated: 0 October 2002 Linear Interpolation p ρ 2 2 = [ 40 30 20] = 20 T p ρ = [ 0 5 20] = 5 p 3 = ρ3 =? 0?? T [ 20 20 20] T 2 600.445
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationSOLUTIONS Aug 2016 exam TFY4102
SOLUTIONS Aug 2016 exam TFY4102 1) In a perfectly ELASTIC collision between two perfectly rigid objects A) the momentum of each object is conserved. B) the kinetic energy of each object is conserved. C)
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationGeneralized Analysis of Stability and Numerical Dispersion in the Discrete-Convolution FDTD Method
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 48, NO. 6, JUNE 2000 887 Generalized Analysis of Stability and Numerical Dispersion in the Discrete-Convolution FDTD Method William A. Beck, Member,
More informationStability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media
Contemporary Mathematics Volume 586 013 http://dx.doi.org/.90/conm/586/11666 Stability and dispersion analysis of high order FDTD methods for Maxwell s equations in dispersive media V. A. Bokil and N.
More informationSpectral Elements. Introduction recalling the elastic wave equation
Spectral Elements Introduction recalling the elastic wave equation The spectral-element method: General concept domain mapping from space-continuous to space-discrete time extrapolation Gauss-Lobatto-Legendre
More informationAcoustic Wave Equation
Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral
More informationA thesis submitted to the University of Manchester for the degree of Doctor of Philosophy in the Faculty of Science and Engineering
INTRODUCTION OF THE DEBYE MEDIA TO THE FILTERED FINITE-DIFFERENCE TIME-DOMAIN METHOD WITH COMPLEX-FREQUENCY-SHIFTED PERFECTLY MATCHED LAYER ABSORBING BOUNDARY CONDITIONS A thesis submitted to the University
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationCHAPTER 4 ANALYSIS AND DESIGN OF THE DUAL INVERTED-F ANTENNA
CHAPTER 4 ANALYSIS AND DESIGN OF THE DUAL INVERTED-F ANTENNA 4.1. Introduction The previous chapter presented the Inverted-F Antenna (IFA) and its variations as antenna designs suitable for use in hand-held
More informationPhysics GRE: Electromagnetism. G. J. Loges 1. University of Rochester Dept. of Physics & Astronomy. xkcd.com/567/
Physics GRE: Electromagnetism G. J. Loges University of Rochester Dept. of Physics & stronomy xkcd.com/567/ c Gregory Loges, 206 Contents Electrostatics 2 Magnetostatics 2 3 Method of Images 3 4 Lorentz
More informationLecture 3. Electric Field Flux, Gauss Law
Lecture 3. Electric Field Flux, Gauss Law Attention: the list of unregistered iclickers will be posted on our Web page after this lecture. From the concept of electric field flux to the calculation of
More informationTransduction Based on Changes in the Energy Stored in an Electrical Field
Lecture 6-1 Transduction Based on Changes in the Energy Stored in an Electrical Field Electric Field and Forces Suppose a charged fixed q 1 in a space, an exploring charge q is moving toward the fixed
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationPublication II Wiley Periodicals. Reprinted by permission of John Wiley & Sons.
Publication II Ilkka Laakso and Tero Uusitupa. 2008. Alternative approach for modeling material interfaces in FDTD. Microwave and Optical Technology Letters, volume 50, number 5, pages 1211-1214. 2008
More informationCONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING FDTD METHODS FOR THE 3D MAXWELL S EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 5 Number 6 Pages 747 784 c 08 Institute for Scientific Computing and Information CONSTRUCTION AND ANALYSIS OF WEIGHTED SEQUENTIAL SPLITTING
More informationModeling of Overhead Power Lines for Broadband PLC Applications.
Modeling of Overhead Power Lines for Broadband PLC Applications. T. A. Papadopoulos, G. K. Papagiannis, D. P. Labridis Power Systems Laboratory Dept. of Electrical & Computer Engineering Aristotle University
More informationANALOGY BETWEEN ELECTROSTATIC FIELD AND HEAT TRANSFER SIMPLE MODELS
ANALOGY BETWEEN ELECTROSTATIC FIELD AND HEAT TRANSFER SIMPLE MODELS PROF. RNDR ING. MILOSLAV KOŠEK, CSC. ING. JAN VODOLAN. Abstract: The perfect analogy allows solving of heat transfer problems by a lot
More informationarxiv:physics/ v1 19 Feb 2004
Finite Difference Time Domain (FDTD) Simulations of Electromagnetic Wave Propagation Using a Spreadsheet David W. Ward and Keith A. Nelson Department of Chemistry Massachusetts Institute of Technology,
More informationAntenna Theory (Engineering 9816) Course Notes. Winter 2016
Antenna Theory (Engineering 9816) Course Notes Winter 2016 by E.W. Gill, Ph.D., P.Eng. Unit 1 Electromagnetics Review (Mostly) 1.1 Introduction Antennas act as transducers associated with the region of
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationCHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD
CHAPTER II MATHEMATICAL BACKGROUND OF THE BOUNDARY ELEMENT METHOD For the second chapter in the thesis, we start with surveying the mathematical background which is used directly in the Boundary Element
More informationECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 1
EE 6340 Intermediate EM Waves Fall 2016 Prof. David R. Jackson Dept. of EE Notes 1 1 Maxwell s Equations E D rt 2, V/m, rt, Wb/m T ( ) [ ] ( ) ( ) 2 rt, /m, H ( rt, ) [ A/m] B E = t (Faraday's Law) D H
More informationLecture 7. Capacitors and Electric Field Energy. Last lecture review: Electrostatic potential
Lecture 7. Capacitors and Electric Field Energy Last lecture review: Electrostatic potential V r = U r q Q Iclicker question The figure shows cross sections through two equipotential surfaces. In both
More informationECE 6340 Intermediate EM Waves. Fall 2016 Prof. David R. Jackson Dept. of ECE. Notes 15
ECE 634 Intermediate EM Waves Fall 6 Prof. David R. Jackson Dept. of ECE Notes 5 Attenuation Formula Waveguiding system (WG or TL): S z Waveguiding system Exyz (,, ) = E( xye, ) = E( xye, ) e γz jβz αz
More informationSUB-NANOSECOND ELECTROMAGNETIC-MICROMA- GNETIC DYNAMIC SIMULATIONS USING THE FINITE- DIFFERENCE TIME-DOMAIN METHOD
Progress In Electromagnetics Research B, Vol. 15, 1 29, 2009 SUB-NANOSECOND ELECTROMAGNETIC-MICROMA- GNETIC DYNAMIC SIMULATIONS USING THE FINITE- DIFFERENCE TIME-DOMAIN METHOD M. M. Aziz School of Engineering,
More informationSolving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation
Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University
More information