An iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems
|
|
- Phebe McDonald
- 5 years ago
- Views:
Transcription
1 An iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems K. Hertel 1,2, S. Yan 1,2, C. Pflaum 1,2, R. Mittra 3 kai.hertel@fau.de 1 Lehrstuhl für Systemsimulation Department Informatik Friedrich-Alexander-Universität Erlangen-Nürnberg 2 School for Advanced Optical Technologies Friedrich-Alexander-Universität Erlangen-Nürnberg 3 Electromagnetic Communication Laboratory Electrical Engineering Department Pennsylvania State University Copper Mountain March 22 nd / 24
2 Basic Outline 1 Model and Discretization 2 Iterative Solver 3 Extrapolation 4 Numerical Results 2 / 24
3 Overview of the model problem High-frequency Maxwell problem (optical regime, Hz < f < Hz, 0.3µm < λ < 1.1µm resp.) Simple box-shaped domain geometries Both media of positive and negative permittivities Rough interfaces between media Periodic and absorbing boundaries (PML) Figure : Schematic of a multi-junction thin-film silicon solar cell 3 / 24
4 Overview of the numerical scheme Finite Difference Frequency Domain discretization (FDFD) Spatial discretization based on Yee cells Finite Integration Technique (FIT) Explicit leap-frog time stepping scheme Algebraically speaking a Richardson method Prony s extrapolation method Figure : EM field in a solar cell structure 4 / 24
5 Under model assumptions valid for optical applications (linearity, isotropy and non-dispersivity), Maxwell s equations can be stated as follows: Faraday s Law: Ampère s Law: µ t H = E σ H (1) ɛ t E = H σe (2) Gauß s Law for Electric Fields: (ɛe) = 0 (3) Gauß s Law for Magnetic Fields: H = 0 (4) 5 / 24
6 In a time-harmonic setting of (1),(2), we seek a solution: H(ξ, t) = Ĥ(ξ) exp(iωt) E(ξ, t) = Ê(ξ) exp(iωt) Therefore, the problem with source terms SÊ, SĤ of frequency ω can be stated as: Find (Ĥ, Ê) that satisfies: iωµ(ξ)ĥ(ξ) = Ê(ξ) σ (ξ)ĥ(ξ) + µ(ξ)s Ĥ (ξ) iωε(ξ)ê(ξ) = Ĥ(ξ) σ(ξ)ê(ξ) + ε(ξ)s Ê (ξ) with ξ Ω R 3 a box-shaped domain. 6 / 24
7 After discretization in time and space (FDFD, Yee cells, explicit leap-frog), we get the following iterative scheme: µ (α 2 Ĥ k+1 τ Ĥk) = α h Êk σ (Ĥk + α 2 Ĥ k+1) + µsĥ 2 ε ( α 2 Ê k+1 τ Êk) = α h Ĥk+1 σ (Êk + α 2 Ê k+1) + εsê 2 for time steps τ > 0, where α := α(τ) = exp (iω 12 ) τ and ( H k (ξ, kτ) = Ĥk (ξ)α 2k, E ξ, 2k + 1 ) τ = 2 Êk (ξ)α 2k+1 ( H ξ, 2k + 1 ) τ 1 ( ) 2 2 α2k Ĥ k (ξ) + α 2 Ĥ k+1 (ξ) 7 / 24
8 Introduce a local modification to the iteration scheme (THIIM 1 ) to deal with instabilities in media of negative permittivity ε < 0: ε ( α 2 Ê k τ Êk+1) = α h Ĥk+1 σ ( α 2 Ê k + 2 Êk+1) + εsê This yields a sparse system matrix consisting of blocks of the following kind: M := 2µ τσ α 2 (2µ+τσ ) 2τ α(2ε+τσ) h 2µ τσ α 2 (2µ+τσ ) 2τα 2 ε +τσ h 2µ τσ α 2 (2µ+τσ ) 2τ α(2µ+τσ ) h 2ε τσ α 2 (2ε+τσ) 2τ α(2ε+τσ) h 2τ α(2µ+τσ ) h α 2 (2 ε τσ) 2 ε +τσ 2τα 2 ε +τσ h 2τ α(2µ+τσ ) h with eigenpairs (λ m, e m ), e m = 1, λ m < 1 m for 0 < τ sufficiently small and M(τ) = m e m λ m e T m τ 0 Id. 1 Pflaum, Rahimi: An iterative solver for the FDFD method for the simulation of materials with negative permittivity, Num.Lin.Alg.Appl., / 24
9 Define algebraic state variables and load vectors: x k (Ĥk :=, Êk) T (Ĥk ) T =, ε>0 Êk, ( ) T ε<0 Êk, b := τ SĤ, SÊ and restate the discretized time-harmonic Maxwell problem: Find a solution x = (Ĥ, Ê)T R d such that: or equivalently: Solve the system x = Mx + b (Id M)x = b 9 / 24
10 Explicit FDFD iteration scheme: x 0 = 0 x k+1 = Mx k + b = k 1 M κ b κ=0 converges slowly 2. Approximation error: x x k = M(x x k 1 ) = M k (x x 0 ) = M k x = M κ b κ=k 2 Recall the corresponding Neumann series: k=0 M k = (Id M) 1 10 / 24
11 High iteration numbers (typically ) to convergence up to a suitable prescribed tolerance 0.2 Real Imag Number of Iterations Wave Length [nm] Figure : Transient signal and iterations to convergence up to a tolerance of / 24
12 How to speed up convergence? x k+1 x k = M k b = m e m λ k me T m Idea: Use low order model to capture the (few) dominant eigenvalue contributions of the system matrix to the iterates. M e m λ m e T m m {n N: λ n >θ} Problem: How to get an estimate of the dominant eigenpairs (or their effect on given iterates)? Approach: Use an extrapolation method based on component-wise samples of successive iterates Pros: Allows for observation of oscillatory behavior caused by resonant modes. Cons: Decouples the degrees of freedom. 12 / 24
13 Step 1: Establish dependence of iterates on previous iterates x k l s x k l 1 s... x k 1 s u 1 x k l+1 x k s l+1 s... x P s := k 2 s....., u := u 2., r x k l+2 s s :=. x k n 1 s x k n 2 s... x k n l u s l x kn s We solve for u: Equivalent to: P s u = r s l m=0 u m x k t m s = 0 t = l n, where u 0 = 1 (5) 13 / 24
14 Step 2: Exponential basis Represent x k s as linear combination of exponentials z k = exp(ifk) Complex-valued f relates to frequency (Re f) and attenuation ( Im f) Determine roots {z m : p(z m ) = 0} of: p(z) = 1 + l u m z l+1 m (6) m=1 Roots z m = exp(if m ) represent distinct frequencies of the system (as detected in the respective degree of freedom s corresponding to some spatial position ξ Ω). 14 / 24
15 Step 3: Determine amplitudes Q s := z k 1 1 z k z k 1 l z k z k 2 l....., a := z1 kn z2 kn... z kn l a 1 a 2. a l, w s := x k 1 s x k 2 s. x kn s Denote iterates x kt s as a superposition of frequencies: x kt s = l a m exp(if m t) m=1 Solve system to get corresponding amplitudes: Q s a = w s Filter amplitude-frequency pairs to get behavior for t 15 / 24
16 Remarks System generally overdetermined (n 2l) Solve in the least squares sense Leads to singular value decomposition (SVD) Minimizes: Numerically robust P s u r s 2 = inf v R l P s v r s 2 (7) Q s a w s 2 = inf v R l Q s v w s 2 (8) 16 / 24
17 Typical distribution of amplitude and frequency pairs as determined by Prony s method 0.1 Distribution of Frequencies and Amplitudes Amplitude Frequency Figure : The signal s frequency-amplitude pairs (abs) Stabilization: Attenutation a = Im f. Discard frequencies with a > ε, assume that frequencies a ε are / 24
18 Convergence of a single degree of freedom: Order 6 extrapolation error vs iterative error, filter by frequency, 0.97 frequencies used on average 10-1 extrapolation iterated signal Figure : Error: Iterated vs extrapolated signal 18 / 24
19 Convergence of the whole simulation domain: Figure : 2 Error: Iterated vs extrapolated signal 19 / 24
20 Numerical Results Figure : Experimental and simulation results 20 / 24
21 Figure : Absorption for smooth and rough interfaces Figure : Quantum efficiency under different angles of incidence 21 / 24
22 Figure : Surface plasmons Figure : Scattering behavior of rough interfaces 22 / 24
23 Summary Explicit iteration scheme for high-frequency Maxwell problems Stable for negative permittivity Extrapolation reduces error Future work How to improve the extrapolation accuracy? (e.g. by choice of better samples) Can a hierarchical basis approach reduce local inaccuracies? What about the divergence of the extrapolated signal? (a constraint in the original system of PDEs) How to make use of multigrid methods? 23 / 24
24 Thank you for your attention. References: A. Taflove, S. Hagness: Computational Electrodynamics, Artech House, 2005 C. Pflaum, Z. Rahimi: An Iterative Solver for the Finite-Difference Frequency Domain (FDFD) Method for the Simulation of Materials with Negative Permittivity, Numerical Linear Algebra with Applications, 2010, DOI: /nla.746 C. Jandl, K. Hertel, C. Pflaum, H. Stiebig: Simulation of Thin-film Silicon Solar Cells with Integrated AFM Scans for Oblique Incident Waves, Proceedings of PVSEC, 2011 R. Prony: Essai expérimental et analytique sur les lois de la dilatibilite des fluides elastiques et sur celles de la force expansive de la vapeur de l eau et de la vapeur de l alkool a differentes temperatures, J. l Ecole Polytechnique, Paris, 1795 M. van Blaricum, R. Mittra: A Technique for Extracting the Poles and Residues of a System Directly from Its Transient Response, IEEE Transactions on Antennas and Propagation, 1975 R. Kumaresan, D. Tufts: Estimating the Parameters of Exponentially Damped Sinusoids and Pole-Zero Modeling in Noise, IEEE Transactions on Acoustics, Speech, and Signal Processing, / 24
Theoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9
WiSe 202 20.2.202 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg
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 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 informationFINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC
FINITE-DIFFERENCE FREQUENCY-DOMAIN ANALYSIS OF NOVEL PHOTONIC WAVEGUIDES Chin-ping Yu (1) and Hung-chun Chang (2) (1) Graduate Institute of Electro-Optical Engineering, National Taiwan University, Taipei,
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 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 informationFinite Difference Solution of Maxwell s Equations
Chapter 1 Finite Difference Solution of Maxwell s Equations 1.1 Maxwell s Equations The principles of electromagnetism have been deduced from experimental observations. These principles are Faraday s law,
More informationA multivariate generalization of Prony s method
A multivariate generalization of Prony s method Ulrich von der Ohe joint with Stefan Kunis Thomas Peter Tim Römer Institute of Mathematics Osnabrück University SPARS 2015 Signal Processing with Adaptive
More informationElectromagnetically Induced Flows in Water
Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationComparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses
Journal of the Korean Physical Society, Vol. 58, No. 4, April 2011, pp. 1014 1020 Comparison Study of the Band-gap Structure of a 1D-Photonic Crystal by Using TMM and FDTD Analyses Jian-Bo Chen, Yan Shen,
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 informationPart VIII. Interaction with Solids
I with Part VIII I with Solids 214 / 273 vs. long pulse is I with Traditional i physics (ICF ns lasers): heating and creation of long scale-length plasmas Laser reflected at critical density surface Fast
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 informationFinite-Difference Time-Domain and Beam Propagation Methods for Maxwell s Equations
Finite-Difference Time-Domain and Beam Propagation Methods for Maxwell s Equations Wolfgang Freude and Jan Brosi Institute of High-Frequency and Quantum Electronics (IHQ), University of Karlsruhe, Germany
More informationPropagation of Surface Plasmon Polariton in the Single Interface of Gallium Lanthanum Sulfide and Silver
PHOTONIC SENSORS / Vol., No., : 58 6 Propagation of Surface Plasmon Polariton in the Single Interface of Gallium Lanthanum Sulfide and Silver Rakibul Hasan SAGOR, Md. Ghulam SABER *, and Md. Ruhul AMIN
More informationEffective Conductivity Concept for Modeling Conductor Surface Roughness G. Gold, K. Helmreich Institute of Microwaves and Photonics
Effective Conductivity Concept for Modeling Conductor Surface Roughness Institute of Microwaves and Photonics Friedrich-Alexander-University Erlangen-Nuremberg gerald.gold@fau.de klaus.helmreich@fau.de
More informationOverview. 1. What range of ε eff, µ eff parameter space is accessible to simple metamaterial geometries? ``
MURI-Transformational Electromagnetics Innovative use of Metamaterials in Confining, Controlling, and Radiating Intense Microwave Pulses University of New Mexico August 21, 2012 Engineering Dispersive
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 informationNew Scaling Factors of 2-D Isotropic-Dispersion Finite Difference Time Domain (ID-FDTD) Algorithm for Lossy Media
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 56, NO., FEBRUARY 8 63 is shown by squares. Each point is computed by averaging reflection coefficient values calculated for each component of the regular
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Tutorial 14
WiSe 0 9.0.03 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg
More informationII Theory Of Surface Plasmon Resonance (SPR)
II Theory Of Surface Plasmon Resonance (SPR) II.1 Maxwell equations and dielectric constant of metals Surface Plasmons Polaritons (SPP) exist at the interface of a dielectric and a metal whose electrons
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 informationFDTD Simulations of Surface Plasmons Using the Effective Permittivity Applied to the Dispersive Media
American Journal of Electromagnetics and Applications 2017; 5(2): 14-19 http://www.sciencepublishinggroup.com/j/ajea doi: 10.11648/j.ajea.20170502.11 ISSN: 2376-5968 (Print); ISSN: 2376-5984 (Online) Review
More informationA RIGOROUS TWO-DIMENSIONAL FIELD ANALYSIS OF DFB STRUCTURES
Progress In Electromagnetics Research, PIER 22, 197 212, 1999 A RIGOROUS TWO-DIMENSIONAL FIELD ANALYSIS OF DFB STRUCTURES M. Akbari, M. Shahabadi, and K. Schünemann Arbeitsbereich Hochfrequenztechnik Technische
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 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 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 informationarxiv:physics/ v1 [physics.optics] 29 Aug 2005
INTERNATIONAL JOURNAL OF NUMERICAL MODELLING: ELECTRONIC NETWORKS, DEVICES AND FIELDS Int. J. Numer. Model. ; : 6 [Version: /9/8 v.] arxiv:physics/589v [physics.optics] 9 Aug 5 Wave scattering by metamaterial
More informationTHE ADI-FDTD METHOD INCLUDING LUMPED NET- WORKS USING PIECEWISE LINEAR RECURSIVE CON- VOLUTION TECHNIQUE
Progress In Electromagnetics Research M, Vol. 30, 67 77, 203 THE ADI-FDTD METHOD INCLUDING LUMPED NET- WORKS USING PIECEWISE LINEAR RECURSIVE CON- VOLUTION TECHNIQUE Fen Xia, Qing-Xin Chu *, Yong-Dan Kong,
More informationELECTROMAGNETIC FIELDS AND WAVES
ELECTROMAGNETIC FIELDS AND WAVES MAGDY F. ISKANDER Professor of Electrical Engineering University of Utah Englewood Cliffs, New Jersey 07632 CONTENTS PREFACE VECTOR ANALYSIS AND MAXWELL'S EQUATIONS IN
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 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 information12A Reflection & Transmission (Normal Incidence)
12A Reflection & Transmission (Normal Incidence) Topics: Reflection and transmission, boundary conditions, complex exponentials. Summary: Students begin by expressing in exponential notation the boundary
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 informationEXTENSIONS OF THE COMPLEX JACOBI ITERATION TO SIMULATE PHOTONIC WAVELENGTH SCALE COMPONENTS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux Eds c TU Delft, The Netherlands, 2006 EXTENSIONS OF THE COMPLEX JACOBI ITERATION TO SIMULATE PHOTONIC
More informationGreen s functions for planarly layered media
Green s functions for planarly layered media Massachusetts Institute of Technology 6.635 lecture notes Introduction: Green s functions The Green s functions is the solution of the wave equation for a point
More informationHybrid Simulation Method ISSS-10 Banff 2011
Hybrid Simulation Method ISSS-10 Banff 2011 David Burgess Astronomy Unit Queen Mary University of London With thanks to Dietmar Krauss-Varban Space Plasmas: From Sun to Earth Space Plasma Plasma is (mostly)
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationqq k d Chapter 16 Electric and Magnetic Forces Electric charge Electric charges Negative (electron) Positive (proton)
Chapter 16 Electric and Magnetic Forces Electric charge Electric charges Negative (electron) Positive (proton) Electrons and protons in atoms/molecules Ions: atoms/molecules with excess of charge Ions
More informationElectromagnetic Relaxation Time Distribution Inverse Problems in the Time-domain
Electromagnetic Relaxation Time Distribution Inverse Problems in the Time-domain Prof Nathan L Gibson Department of Mathematics Joint Math Meeting Jan 9, 2011 Prof Gibson (OSU) Inverse Problems for Distributions
More informationDispersion of Homogeneous and Inhomogeneous Waves in the Yee Finite-Difference Time-Domain Grid
280 IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 49, NO. 2, FEBRUARY 2001 Dispersion of Homogeneous and Inhomogeneous Waves in the Yee Finite-Difference Time-Domain Grid John B. Schneider,
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 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 information22 Phasor form of Maxwell s equations and damped waves in conducting media
22 Phasor form of Maxwell s equations and damped waves in conducting media When the fields and the sources in Maxwell s equations are all monochromatic functions of time expressed in terms of their phasors,
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 informationACM 106a: Lecture 1 Agenda
1 ACM 106a: Lecture 1 Agenda Introduction to numerical linear algebra Common problems First examples Inexact computation What is this course about? 2 Typical numerical linear algebra problems Systems of
More informationUsing ParExPDE for the numerical solution of bioelectric field problems
Using ParExPDE for the numerical solution of bioelectric field problems Christoph Freundl Friedrich-Alexander-Universität Erlangen-Nürnberg Christoph.Freundl@informatik.uni-erlangen.de Harald Köstler Friedrich-Alexander-Universität
More informationMagnetotelluric (MT) Method
Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to
More information4. The interaction of light with matter
4. The interaction of light with matter The propagation of light through chemical materials is described by a wave equation similar to the one that describes light travel in a vacuum (free space). Again,
More informationInside-out electromagnetic cloaking
Inside-out electromagnetic cloaking Nina A. Zharova 1,2, Ilya V. Shadrivov 1, and Yuri S. Kivshar 1 1 Nonlinear Physics Center, Research School of Physical Sciences and Engineering, Australian National
More informationANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA
Progress In Electromagnetics Research Letters, Vol. 22, 83 93, 2011 ANALYSIS OF DISPERSION RELATION OF PIECEWISE LINEAR RECURSIVE CONVOLUTION FDTD METHOD FOR SPACE-VARYING PLASMA X. Ai Science and Technology
More informationA Novel Design of Photonic Crystal Lens Based on Negative Refractive Index
PIERS ONLINE, VOL. 4, NO. 2, 2008 296 A Novel Design of Photonic Crystal Lens Based on Negative Refractive Index S. Haxha 1 and F. AbdelMalek 2 1 Photonics Group, Department of Electronics, University
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 informationLight-trapping by diffraction gratings in silicon solar cells
Light-trapping by diffraction gratings in silicon solar cells Silicon: a gift of nature Rudolf Morf, Condensed Matter Theory, Paul Scherrer Institute Benefits from improved absorption Diffraction gratings
More informationSUPPLEMENTARY INFORMATION
An effective magnetic field from optically driven phonons T. F. Nova 1 *, A. Cartella 1, A. Cantaluppi 1, M. Först 1, D. Bossini 2 #, R. V. Mikhaylovskiy 2, A.V. Kimel 2, R. Merlin 3 and A. Cavalleri 1,
More informationScalar electromagnetic integral equations
Scalar electromagnetic integral equations Uday K Khankhoje Abstract This brief note derives the two dimensional scalar electromagnetic integral equation starting from Maxwell s equations, and shows how
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 informationA MATLAB GUI FOR SIMULATING THE PROPAGATION OF THE ELECTROMAGNETIC FIELD IN A 2-D INFINITE SPACE
A MATLAB GUI FOR SIMULATING THE PROPAGATION OF THE ELECTROMAGNETIC FIELD IN A 2-D INFINITE SPACE Ioana SĂRĂCUŢ Victor POPESCU Marina Dana ŢOPA Technical University of Cluj-Napoca, G. Bariţiu Street 26-28,
More informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationPublication I Institute of Physics Publishing (IOPP) Reprinted by permission of Institute of Physics Publishing.
Publication I Ilkka Laakso, Sami Ilvonen, and Tero Uusitupa. 7. Performance of convolutional PML absorbing boundary conditions in finite-difference time-domain SAR calculations. Physics in Medicine and
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 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 informationECE 107: Electromagnetism
ECE 107: Electromagnetism Notes Set 1 Instructor: Prof. Vitaliy Lomakin Department of Electrical and Computer Engineering University of California, San Diego, CA 92093 1 Introduction (1) atom Electromagnetism
More informationTransmission Line Model for Rectangular Waveguides accurately incorporating Loss Effects
Transmission Line Model for Rectangular Waveguides accurately incorporating Loss Effects Konstantin Lomakin konstantin.lomakin@fau.de Institute of Microwaves and Photonics Friedrich-Alexander-Universität
More informationIntroduction to Nonlinear Optics
Introduction to Nonlinear Optics Prof. Cleber R. Mendonca http://www.fotonica.ifsc.usp.br Outline Linear optics Introduction to nonlinear optics Second order nonlinearities Third order nonlinearities Two-photon
More informationElectromagnetic Waves
Physics 8 Electromagnetic Waves Overview. The most remarkable conclusion of Maxwell s work on electromagnetism in the 860 s was that waves could exist in the fields themselves, traveling with the speed
More informationMixed Finite Elements Method
Mixed Finite Elements Method A. Ratnani 34, E. Sonnendrücker 34 3 Max-Planck Institut für Plasmaphysik, Garching, Germany 4 Technische Universität München, Garching, Germany Contents Introduction 2. Notations.....................................
More informationNumerical Analysis of Electromagnetic Fields in Multiscale Model
Commun. Theor. Phys. 63 (205) 505 509 Vol. 63, No. 4, April, 205 Numerical Analysis of Electromagnetic Fields in Multiscale Model MA Ji ( ), FANG Guang-You (ྠ), and JI Yi-Cai (Π) Key Laboratory of Electromagnetic
More informationA Spectral Time-Domain Method for Computational Electrodynamics
A Spectral Time-Domain Method for Computational Electrodynamics James V. Lambers Abstract Block Krylov subspace spectral (KSS) methods have previously been applied to the variable-coefficient heat equation
More informationA Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere
A Three-Fluid Approach to Model Coupling of Solar Wind-Magnetosphere-Ionosphere- Thermosphere P. Song Center for Atmospheric Research University of Massachusetts Lowell V. M. Vasyliūnas Max-Planck-Institut
More informationThe Interaction of Light and Matter: α and n
The Interaction of Light and Matter: α and n The interaction of light and matter is what makes life interesting. Everything we see is the result of this interaction. Why is light absorbed or transmitted
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationContour Integral Method for the Simulation of Accelerator Cavities
Contour Integral Method for the Simulation of Accelerator Cavities V. Pham-Xuan, W. Ackermann and H. De Gersem Institut für Theorie Elektromagnetischer Felder DESY meeting (14.11.2017) November 14, 2017
More informationA new method for sensitivity analysis of photonic crystal devices
A new method for sensitivity analysis of photonic crystal devices Georgios Veronis, Robert W. Dutton, and Shanhui Fan Department of Electrical Engineering, Stanford University, Stanford, California 9435
More informationFast and Efficient Modeling of Photonic Crystal Devices
Fast and Efficient Modeling of Photonic Crystal Devices D. Pissoort, E. Michielssen, F. Olyslager, and D.De Zutter Outline 1. Introduction 2. Part I: 2D Photonic Crystal Devices 3. Part II: Photonic Crystal
More informationDynamic Response. Assoc. Prof. Enver Tatlicioglu. Department of Electrical & Electronics Engineering Izmir Institute of Technology.
Dynamic Response Assoc. Prof. Enver Tatlicioglu Department of Electrical & Electronics Engineering Izmir Institute of Technology Chapter 3 Assoc. Prof. Enver Tatlicioglu (EEE@IYTE) EE362 Feedback Control
More informationEFFICIENT SIMULATIONS OF PERIODIC STRUC- TURES WITH OBLIQUE INCIDENCE USING DIRECT SPECTRAL FDTD METHOD
Progress In Electromagnetics Research M, Vol. 17, 101 111, 2011 EFFICIENT SIMULATIONS OF PERIODIC STRUC- TURES WITH OBLIQUE INCIDENCE USING DIRECT SPECTRAL FDTD METHOD Y. J. Zhou, X. Y. Zhou, and T. J.
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 informationOverview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001).
Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) T.Thio et al., Optics Letters 6, 197-1974
More informationS-matrix approach for calculations of the optical properties of metallic-dielectric photonic crystal slabs
S-matrix approach for calculations of the optical properties of metallic-dielectric photonic crystal slabs N. I. Komarevskiy1,2, T. Weiss3, and S. G. Tikhodeev2 1 Faculty of Physics, Lomonosov Moscow State
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 informationLight and Matter. Thursday, 8/31/2006 Physics 158 Peter Beyersdorf. Document info
Light and Matter Thursday, 8/31/2006 Physics 158 Peter Beyersdorf Document info 3. 1 1 Class Outline Common materials used in optics Index of refraction absorption Classical model of light absorption Light
More informationA new Technique for Simulating Semiconductor Laser Resonators
A new Technique for Simulating Semiconductor Laser Resonators Britta Heubeck Christoph Pflaum Department of System Simulation (LSS) University of Erlangen-Nuremberg, Germany NUSOD conference, September
More informationIN TIME-DOMAIN methods for electromagnetic analysis,
1202 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 65, NO. 3, MARCH 2017 Time-Domain Method Having a Naturally Diagonal Mass Matrix Independent of Element Shape for General Electromagnetic Analysis
More informationSatellite Remote Sensing SIO 135/SIO 236. Electromagnetic Radiation and Polarization
Satellite Remote Sensing SIO 135/SIO 236 Electromagnetic Radiation and Polarization 1 Electromagnetic Radiation The first requirement for remote sensing is to have an energy source to illuminate the target.
More informationLecture 10 Light-Matter Interaction Part 4 Surface Polaritons 2. EECS Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.
Lecture 10 Light-Matter Interaction Part 4 Surface Polaritons 2 EECS 598-002 Winter 2006 Nanophotonics and Nano-scale Fabrication P.C.Ku Schedule for the rest of the semester Introduction to light-matter
More informationOn Approximate Cloaking by Nonsingular Transformation Media
On Approximate Cloaking by Nonsingular Transformation Media TING ZHOU MIT Geometric Analysis on Euclidean and Homogeneous Spaces January 9, 2012 Invisibility and Cloaking From Pendry et al s paper J. B.
More informationPrinciples of Mobile Communications
Communication Networks 1 Principles of Mobile Communications University Duisburg-Essen WS 2003/2004 Page 1 N e v e r s t o p t h i n k i n g. Wave Propagation Single- and Multipath Propagation Overview:
More informationChapter Three: Propagation of light waves
Chapter Three Propagation of Light Waves CHAPTER OUTLINE 3.1 Maxwell s Equations 3.2 Physical Significance of Maxwell s Equations 3.3 Properties of Electromagnetic Waves 3.4 Constitutive Relations 3.5
More informationOscillations and Electromagnetic Waves. March 30, 2014 Chapter 31 1
Oscillations and Electromagnetic Waves March 30, 2014 Chapter 31 1 Three Polarizers! Consider the case of unpolarized light with intensity I 0 incident on three polarizers! The first polarizer has a polarizing
More informationExponential Decomposition and Hankel Matrix
Exponential Decomposition and Hankel Matrix Franklin T Luk Department of Computer Science and Engineering, Chinese University of Hong Kong, Shatin, NT, Hong Kong luk@csecuhkeduhk Sanzheng Qiao Department
More informationGeneralized Pencil Of Function Method
Generalized Pencil Of Function Method In recent times, the methodology of approximating a function by a sum of complex exponentials has found applications in many areas of electromagnetics For example
More informationAccelerating Nesterov s Method for Strongly Convex Functions
Accelerating Nesterov s Method for Strongly Convex Functions Hao Chen Xiangrui Meng MATH301, 2011 Outline The Gap 1 The Gap 2 3 Outline The Gap 1 The Gap 2 3 Our talk begins with a tiny gap For any x 0
More informationAn inverse problem for eddy current equations
An inverse problem for eddy current equations Alberto Valli Department of Mathematics, University of Trento, Italy In collaboration with: Ana Alonso Rodríguez, University of Trento, Italy Jessika Camaño,
More informationIMPLEMENTING THE PERFECTLY MATCHED LAYER ABSORBING BOUNDARY CONDITION WITH MIMETIC DIFFERENCING SCHEMES
Progress In Electromagnetics Research, PIER 32, 383 411, 21 IMPLEMENTING THE PERFECTLY MATCHED LAYER ABSORBING BOUNDARY CONDITION WITH MIMETIC DIFFERENCING SCHEMES M. W. Buksas Los Alamos National Laboratory,
More informationFinding damage in Space Shuttle foam
Finding damage in Space Shuttle foam Nathan L. Gibson gibsonn@math.oregonstate.edu In Collaboration with: Prof. H. T. Banks, CRSC Dr. W. P. Winfree, NASA OSU REU p. 1 Motivating Application The particular
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 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 informationLecture 2 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell
Lecture Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell 1. Dispersion Introduction - An electromagnetic wave with an arbitrary wave-shape
More information