Part E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
|
|
- Polly Whitehead
- 6 years ago
- Views:
Transcription
1 Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann
2 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 ordinary differential equation (ODEs) Numeric integration using standard methods of numerical analysis utilize the system s special structure, using a central difference quotient f ( t t) f ( t t) f t O t t ( ) ( ) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik
3 Discretization of time: FARADAY s law equidistant time steps ( temporal grid ) h ( m) Transient Fields d h dt h( tm ), tm t0 m t with time step width t 1 () t MC e( t) central DQ for the time derivative h h ( m1) ( m) t d h dt ( m ) 1 M Ce 1 (m 1 ) half time step ( m) h e d h dt ( m ) 1 ( m ) 1 ( m1) h ( m) h t m t m+1/ t m+1 t t m+ t staggered time ais (allocation of quantities) (analogous to staggered spatial allocation in the dual grid system) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 3
4 Transient Fields Discretization of time: AMPERE s law equidistant time steps e ( m ) e( t t m t ), t ( ) 1 m m 1 d e dt (w/o currents!) 1 () t M Ch() t central DQ for the time derivative e ( m ) ( m ) 3 1 ( m) h e t e d h dt d e dt ( m ) 1 ( m ) 1 ( m1) d e dt M Ch ( m1) ( m1) h 1 ( m1) e ( m ) 3 again at full time steps ( m) h t m t m+1/ t m+1 t m+3/ t m+ t Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 4
5 Construction of the time integration scheme d h dt d e dt 1 ( m ) ( m1) ( m) h t ( m1) ( m ) ( m ) h e 3 1 e t M M 1 Ce Ch 1 ( m ) 1 ( m1) Transient Fields h h t M Ce ( m1) ( m) 1 ( m 1 ) m 3 1 m e e t M C ( ) ( ) 1 ( m1) update equations h ( m) h e ( m ) 1 ( m1) h e ( m ) 3 ( m) h t m t m+1/ t m+1 t m+3/ t m+ t leapfrog scheme Etensions Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 5
6 Leapfrog Scheme Special time integration scheme for MAXWELL Grid Equations idea: finite differences on an equidistant time ais motivated by the structure of the ODE system (of MAXWELL s equations) (so far) only the simplest form of the curl equations properties staggered allocation of the quantities on time ais eplicit update equations with high efficiency (no matri inversion, no system of equations must be solved) local accuracy: nd order in the temporal FD approimation global accuracy? stability? proposed for MAXWELL s equations by K.S. YEE 1966 (later referred to as: FDTD = finite difference time domain) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 6
7 Small eample Leapfrog Scheme: Eample 33 grid, PEC boundaries, unit lengths = y = z = 1m only 5 non-zero components: 1, 4 free indeing epected solution: one cavity mode (type E z,110 with f res 106 MHz) e h time domain simulation no current (0) start with e 1 vary time step width t (start with ca. 10 steps per time period) transient results of (normalized) fields h 1 h h 4 e 1 h 3 Matlab... Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 7
8 Leapfrog Scheme: Eample (Scaled) transient fields, varying time step width t t 0 fitted cos function Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 8
9 Leapfrog Scheme: Eample (Scaled) transient fields, varying time step width t t 0 fitted cos function Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 9
10 Leapfrog Scheme: Eample (Scaled) transient fields, varying time step width t 3t 0 fitted cos function Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 10
11 Leapfrog Scheme: Eample (Scaled) transient fields, varying time step width t 4 t 0 fitted cos function Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 11
12 Leapfrog Scheme: Eample (Scaled) transient fields, varying time step width t 4 t 0 semi-logarithmic Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 1
13 Leapfrog Scheme: Eample Impact of the time step width in the leapfrog scheme 1. accuracy very accurate result for (at least) 10 time steps per time period larger time steps: phase error (effective frequency ) insufficient sampling rate: inaccurate results. stability time integration is instable for a too large t! eponential increase (comple frequency) the leapfrog scheme is only conditionally stable (strict) stability limit: t t ma s (eperimental value) calculate a priori? instability: meaningless results Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 13
14 Properties of the discrete solutions Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 14
15 Properties of the discrete solutions Properties of the Solution eample D-grid with absorbing boundaries radiation source (antenna with reference impedance) in the middle time domain: ecitation signal = modulated pulse fields are monitored at several frequencies (FOURIER transform) online demo... Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 15
16 Properties of the Solution Cylindrical wave in a Cartesian FIT grid f 50 MHz D grid - 11 points - symmetry boundaries - = y = 0.5m 1 f f f f f 50 MHz: 6m=1 100 MHz: 3m=6 10 MHz:.5m=5 150 MHz: m=4 00 MHz: 1.5m=3 Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 16
17 Properties of the Solution Cylindrical wave in a Cartesian FIT grid f 150 MHz D grid - 11 points - symmetry boundaries - = y = 0.5m f 50 MHz: 6m=1 4 f 100 MHz: 3m=6 f 10 MHz:.5m=5 f 150 MHz: m=4 f 00 MHz: 1.5m=3 Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 17
18 Grid Solution for Plane Waves Cylindrical wave in a Cartesian FIT grid for fied grid: frequency spatial sampling rate / spatial sampling dominates accuracy (ok: / > 10) discrete wave length depends on direction of propagation f=50 MHz, / 1 f=150 MHz, / 4 Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 18
19 Grid Solution for Plane Waves Wave propagation in grid, simplest case plane waves j( tk r ) j t j( k k y k z) E, H ~ e e e y z (continuous) dispersion relation: k ky kz c (relation of: wave number / frequency = space / time dependency) infinite grid (no specific boundary conditions), equidistant (constant grid steps, y, z), vacuum (constant: 0, 0, s = 0 ) (numerical) dispersion relation in the grid? Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 19
20 Grid Solution for Plane Waves 1) Properties of the spatial discretization impose plane wave for grid quantities jk EH, e r e T y eˆy bˆz eˆ e y e jk r e y T ey T reference point all components can be represented by amplitudes (quantities eˆ, eˆ, eˆ, bˆ, at reference point) phase factors y z jk jk y jkzz y z y T e, T e, T e Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 0
21 Grid Solution for Plane Waves Local derivation of the discrete wave equation eˆ T y eˆ ˆ ˆ ˆ ˆ y ezt eyt z e j b y z eˆy bˆz eˆ T y eˆ ˆ ˆ ˆ ˆ z et ezt e j b z y eˆ eˆ ˆ ˆ ˆ ˆ eyt et y e j b y z 33 version of FARADAY s Law for wave amplitudes ˆ 0 ( T 1) ( 1) ˆ z Ty e b ( T 1) 0 ( 1) ˆ ˆ z T e y j b y ( T 1) ( 1) 0 ˆ ˆ y T e z b z Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 1
22 Grid Solution for Plane Waves Local derivation of the discrete wave equation eˆy eˆ bˆz T y eˆ T y ˆ 0 ( T 1) ( T 1) ˆ z y e b ( T 1) 0 ( T 1) eˆ ˆ j b z y y ( T 1) ( T 1) 0 eˆ ˆ y z b z eˆ Cˆ eˆ j bˆ analogous: insert into each other: ˆ Chˆ j dˆ with ˆ C C dˆ Mˆ eˆ, bˆ Mˆ hˆ ˆ H ˆ M CM Cˆ e ˆ 1 ˆ 1 ˆ eˆ (a local 33- wave equation) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik
23 Solutions of Grid Solution for Plane Waves ˆ M CM Cˆ e ˆ 1 ˆ 1 ˆ eˆ left hand side (spatial discretization): compute eigenvalues ˆ of Aˆ ˆ Mˆ CMˆ Cˆ 1 1 ( Aˆ ˆ ) (...) ˆ 0 1 ˆ static solution sink sinkyy sinkzz y z,3 c plane waves ( polarizations!) MAXWELL s grid equations require: ˆ Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 3
24 Numerical Dispersion Relation Numerical dispersion relation (NDR) of the FIT Gitterdispersionsgleichung k sin k y y k z sin sin z y z c condition for the propagation of plane waves in the grid only valid eactly in the given simplest case dispersion in the grid depends on direction of wave propagation! converges towards continuous dispersion relation (...):, y, z 0 k k k y z c allows a quantitative approimation of the dispersion error (...) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 4
25 Numerical Dispersion Relation Directional wave propagation in the grid k 1 k0 0 sin k c k 1/ k1/ 0 sin k ( ) ( ) c (with y) (both converging) different wave lengths (for identical frequency) k this eplains the results for the cylindrical wave eperiment! Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 5
26 Numerical Dispersion Relation quantitatively: cos k k diskret sin 0 y k Polar plot representation (n=4): diskret kontin. sin cos / / n sin sin / / n n n kontin. n diskret diskret 1 3 cos4 6 n 4 kontin. Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 6
27 Numerical Dispersion: Time Time harmonic wave temporal phase factor Tt j t e M CM Ce 1 1 d dt e discrete 1 st derivative Tt 1 eˆ ê T eˆ t t ref t t ref t ref t t t/ -1 1 discrete nd derivative d 1 eˆ ref e eˆ ref ( Tt 1) (cost j sint 1) dt t t t ref d ˆ ˆ ref t t ref e e ( T T ) e (cost j sint cost j sin t) dt t t 1 1 (sin t ) eˆ (1 cos ) ˆ ref t e ref ( t ) ( t ) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 7
28 Numerical Dispersion Relation Numerical dispersion relation (time domain) k sin k y y k z 1 t sin sin sin z y z c t now a complete description of the result above (in the simplest case ) identical epressions for space- and time-dependence converges for step widths 0 (later: temporal error < spatial error) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 8
29 NDR and Stability Stability eperimental result: stable for t < t ma can the stability limit be computed a priori? stable: f(t) = ep(j t) instable: f(t) = ep(j t) ep(at) = ep(j t) only real eigenfrequencies comple eigenfrequency Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 9
30 Stability NDR and Stability the scheme is stable, if all (eigen)frequencies are real in the simplest case: modelled by the numerical grid dispersion relation: eigenfrequency, wave number, step widths k sin k y y k z 1 t sin sin sin z y z c t condition for real frequencies real sin( t ) 1 k, k, k y z worst-case approimation of the left-hand side (wave numbers) sin( k ) 1 l.h.s y z Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 30
31 NDR and Stability Stability ( ) t sin y z c t sin t t c 1 y z! 1 1 t c y z COURANT - FRIEDRICHS - LEVY (CFL) criterion = stability limit (maimum stable time step width) for FIT + leapfrog scheme (or FDTD, resp.) conditionally stable time integration scheme Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 31
32 CFL criterion Leapfrog: Stable Time Step maimum stable time step for the time integration t ma depends on spatial grid step width (and material) for y z : tma c 3 3 only valid in the simples case; otherwise (heuristically): 1 1 t c y z t min i i y z i i i i (worst-case for all grid cells) can easily be evaluated a-priori! 1 9 yields here: tma s (eperimentally: s) c 3m Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 3
33 Leapfrog: Stable Time Step CFL criterion worst-case approimation for all grid cells in the infinite, homogeneous grid t min i i y z i i i i all possible (real) wave numbers k, k y, k z real situation with losses, boundaries, etc. (strictly speaking:) the derivation is no longer valid inhomogeneous, finite grids: worst-case situation often only locally deviation in the demo eample: small grid, limited range of wave numbers k, k y, k z (typically the estimation is much closer) a superior (eact!) formula is t ma = / ma 1 1 ma CC (eigenvalue with largest magnitude of A M CM C ) Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 33
34 Interpretation of TD Eigenvalues 1. Time continuous / frequency domain / t 0 1 rot E 1 rot H t H E t 1 rot H H 1 rot E t E H / t Ht () Im j E/ t Et () FIT Re Sign: operations still missing ( ) already shown: real frequencies for lossless FIT systems the time-continuous system itself is stable Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 34
35 Interpretation of TD Eigenvalues. Time discrete: t > 0 (e.g., leapfrog) each time step advances the phase of each eigensolution by we already know: frequencies only approimate t TD TD eigenvalue of the time stepping scheme j j TD t e e TD the scheme is stable, if j TD t e 1 TD TD ( m1) h the stability of the time stepping scheme depends on how it approimates the eigenfrequencies ( m) h Im ( m1) ( m1/ ) e ( m) Re leapfrog ( ): j ( t) ( t) e 1 t 1 4 t t ma ma Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 35
36 TD Stability fully continuous system (MAXWELL) is stable with real eigenfrequencies phys FIT with C C etc., Others? (FD, FE, ) T FIT semi-discrete system (timecontinuous MAXWELL Grid Eqs.) is stable with real eigenfrequencies Leapfrog with limited time step Others? j TD t e 1, TD fully discrete time-domain system is stable with eigenvalues within unit circle Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 36
37 Accuracy and Stability Accuracy and stability stability criterion determines ma. time step the ratio time step / time period also determines accuracy space: sampling of the wave length by grid spatial error time: sampling of the time period by time ais temporal error coupled by t t ma c quantitative relation between spatial and temporal approimation? is the choice of the maimum stable reasonable w.r.t. accuracy? Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 39
38 Accuracy and stability Accuracy and Stability Eample (1D): plane wave propagating in -direction sin k 1 sin t numerical grid dispersion: c t analyze accuracy: compare FIT results with continuous values k fied cont. = k c vs. FIT (eample: cavity), or fied k,cont. = / c vs. k,fit (eample: cyl. wave) vary the grid step -- sampling rates -- (n steps per wave length) k n n k vary the time step T t m m (m steps per period) t n m Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 40
39 insert: k Accuracy and Stability n n m FIT n n FIT k c n c m m sin m sin si sin with n m T t k c,, cont. relative frequency error: FIT cont. cont. cont. sin m sin m n n 1 for further analysis: series epansion: consider separately: 3 5 sin a a a O( a ) m (spatial error: only due to > 0 ) n (temporal error: only due to t > 0 ) O 4 space 6n n 1 O 4 time 6m m Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 41
40 Accuracy and Stability 1 O 4 space 6n n 1 O 4 time 6m m quantitatively: error < % for (6 n ) 0.0 n lines per lambda rule spatial and temporal errors compensate each other (partially) is it possible to take the best choice m=n? apply CFL criterion c 3 3 t m 3 n c 3 t the stability limit yields errors in the same magnitude, m=n only possible in the (academic) 1D case ( magic time step ) no (considerable) additional effort by stability requirement! Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 4
41 The leapfrog scheme is NOT efficient, if Limitations of Leapfrog the CFL criterion requires a very small time step very small grid steps (misalignments, small geometric details related to ) very small frequencies (e.g. 50 Hz, related to =6.000 km; object size 1m1m1m, = 1cm, t s, T 10 8 t ) good conductors with small skin depth ( use special models) the computational domain is electrically very large (related to ) very high frequencies and large objects: > 10 grid steps per required due to dispersion errors e.g. plane wave (1 GHz) on airplane (30m=100 ): 1000 lines per spatial dimension = 10 9 grid points in 3D! Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 44
42 The leapfrog scheme is VERY efficient Limitations of Leapfrog if the size of the structure is in the range of a few wavelength many technical applications work with or below the fundamental resonance microwave filter, antennas, resonators, connectors, waveguiding structures & transitions,... if one is interested in broadband results ecitation in time domain by pulses with broadband spectral content broadband results (such as S-parameter) within one simulation run details... Rolf Schuhmann TU Berlin, FG Theoretische Elektrotechnik 45
Limitations of the Leap-Frog Scheme
Limitations of the Leap-Frog Scheme 1. Structures with high quality factor (resonators, filter,...): ω Ws frequency stored energy The quality factor Q is defined as Q= = Pv losses long decay time of resonant
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 informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
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
More informationApproximate Analytical Calculation of the Skin Effect in Rectangular Wires
Approimate Analytical Calculation of the Skin Effect in Rectangular Wires Dieter Gerling University of Federal Defense Munich, 85579 Neubiberg, Germany Dieter.Gerling@unibw.de Abstract- It is well-known
More informationThe equations so far... Gauss Law for E Fields. Gauss Law for B Fields. B da. inside. d dt. n C 3/28/2018
The equations so far... Gauss Law for E Fields E da S n 1 Q inside Gauss Law for B Fields B da S n C Faraday s Law d E dl dt S B da n Ampere s Law B dl I C 3/8/18 1 Ampere s Law B dl I inside _ path No
More informationPeriodic Structures in FDTD
EE 5303 Electromagnetic Analsis Using Finite Difference Time Domain Lecture #19 Periodic Structures in FDTD Lecture 19 These notes ma contain coprighted material obtained under fair use rules. Distribution
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 informationElectromagnetic Waves & Polarization
Course Instructor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@utep.edu EE 4347 Applied Electromagnetics Topic 3a Electromagnetic Waves & Polarization Electromagnetic These
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 informationFinite difference modelling, Fourier analysis, and stability
Finite difference modelling, Fourier analysis, and stability Peter M. Manning and Gary F. Margrave ABSTRACT This paper uses Fourier analysis to present conclusions about stability and dispersion in finite
More informationSecondary School Mathematics & Science Competition. Mathematics. Date: 1 st May, 2013
Secondary School Mathematics & Science Competition Mathematics Date: 1 st May, 01 Time allowed: 1 hour 15 minutes 1. Write your Name (both in English and Chinese), Name of School, Form, Date, Se, Language,
More informationPolynomial Chaos Approach for Maxwell s Equations in Dispersive Media
Polynomial Chaos Approach for Maxwell s Equations in Dispersive Media Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 15, 2013 Prof. Gibson (OSU) PC-FDTD
More information1-D Implementation of Maxwell s Equations in MATLAB to Study the Effect of Absorption Using PML
1-D Implementation of Mawell s Equations in MATLAB to Study the Effect of Absorption Using PML 1 Vikas Rathi, Pranav 2 K. Shrivastava, 3 Hemant S Pokhariya ECE Department Graphic Era University Dehradun,
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 informationCoupling Impedance of Ferrite Devices Description of Simulation Approach
Coupling Impedance of Ferrite Devices Description of Simulation Approach 09 May 2012 TU Darmstadt Fachbereich 18 Institut Theorie Elektromagnetischer Felder Uwe Niedermayer 1 Content Coupling Impedance
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 informationClass 30: Outline. Hour 1: Traveling & Standing Waves. Hour 2: Electromagnetic (EM) Waves P30-
Class 30: Outline Hour 1: Traveling & Standing Waves Hour : Electromagnetic (EM) Waves P30-1 Last Time: Traveling Waves P30- Amplitude (y 0 ) Traveling Sine Wave Now consider f(x) = y = y 0 sin(kx): π
More informationPseudospectral 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 informationTime-harmonic form Phasor form. =0 (1.11d)
Chapter 2 Wave in an Unbounded Medium Maxwell s Equations Time-harmonic form Phasor form (Real quantity) (complex quantity) B E = Eˆ = jωbˆ (1.11 a) t D= ρ Dˆ = ρ (1.11 b) D H = J + Hˆ = Jˆ+ jωdˆ ( 1.11
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More information( z) ( ) ( )( ) ω ω. Wave equation. Transmission line formulas. = v. Helmholtz equation. Exponential Equation. Trig Formulas = Γ. cos sin 1 1+Γ = VSWR
Wave equation 1 u tu v u(, t f ( vt + g( + vt Helmholt equation U + ku jk U Ae + Be + jk Eponential Equation γ e + e + γ + γ Trig Formulas sin( + y sin cos y+ sin y cos cos( + y cos cos y sin sin y + cos
More informationEnergy conserving coupling through small apertures in an infinite perfect conducting screen
doi:1.519/ars-13-7-15 Author(s) 15. CC Attribution 3. License. Energy conserving coupling through small apertures in an infinite perfect conducting screen J. Petzold,. Tkachenko, and R. Vick Chair of Electromagnetic
More informationElectromagnetic Waves
Electromagnetic Waves Our discussion on dynamic electromagnetic field is incomplete. I H E An AC current induces a magnetic field, which is also AC and thus induces an AC electric field. H dl Edl J ds
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 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 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 informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Pacet The Group Velocity Phys 375 Overview and Motivation: Last time we looed at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
More informationUniversity of Alberta ENGM 541: Modeling and Simulation of Engineering Systems Laboratory #5
University of Alberta ENGM 54: Modeling and Simulation of Engineering Systems Laboratory #5 M.G. Lipsett, Updated 00 Integration Methods with Higher-Order Truncation Errors with MATLAB MATLAB is capable
More informationElectromagnetic waves in free space
Waveguide notes 018 Electromagnetic waves in free space We start with Maxwell s equations for an LIH medum in the case that the source terms are both zero. = =0 =0 = = Take the curl of Faraday s law, then
More informationCERN Accelerator School. RF Cavities. Erk Jensen CERN BE-RF
CERN Accelerator School RF Cavities Erk Jensen CERN BE-RF CERN Accelerator School, Varna 010 - "Introduction to Accelerator Physics" What is a cavity? 3-Sept-010 CAS Varna/Bulgaria 010- RF Cavities Lorentz
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 informationMath 1160 Final Review (Sponsored by The Learning Center) cos xcsc tan. 2 x. . Make the trigonometric substitution into
Math 60 Final Review (Sponsored by The Learning Center). Simplify cot csc csc. Prove the following identities: cos csc csc sin. Let 7sin simplify.. Prove: tan y csc y cos y sec y cos y cos sin y cos csc
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Packet The Group Velocity Phys 375 Overview and Motivation: Last time we looked at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
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 informationCourse Updates. 2) This week: Electromagnetic Waves +
Course Updates http://www.phys.hawaii.edu/~varner/phys272-spr1/physics272.html Reminders: 1) Assignment #11 due Wednesday 2) This week: Electromagnetic Waves + 3) In the home stretch [review schedule]
More informationAn Explicit and Unconditionally Stable FDTD Method for Electromagnetic Analysis
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 1 An Explicit and Unconditionally Stable FDTD Method for Electromagnetic Analysis Md. Gaffar and Dan Jiao, Senior Member, IEEE Abstract In this paper,
More informationAPPM 1360 Final Exam Spring 2016
APPM 36 Final Eam Spring 6. 8 points) State whether each of the following quantities converge or diverge. Eplain your reasoning. a) The sequence a, a, a 3,... where a n ln8n) lnn + ) n!) b) ln d c) arctan
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationECE 6340 Intermediate EM Waves. Fall Prof. David R. Jackson Dept. of ECE. Notes 17
ECE 634 Intermediate EM Waves Fall 16 Prof. David R. Jacson Dept. of ECE Notes 17 1 General Plane Waves General form of plane wave: E( xz,, ) = Eψ ( xz,, ) where ψ ( xz,, ) = e j( xx+ + zz) The wavenumber
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 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 informationECE 6341 Spring 2016 HW 2
ECE 6341 Spring 216 HW 2 Assigned problems: 1-6 9-11 13-15 1) Assume that a TEN models a layered structure where the direction (the direction perpendicular to the layers) is the direction that the transmission
More informationThe Quantum Theory of Atoms and Molecules
The Quantum Theory of Atoms and Molecules The postulates of quantum mechanics Dr Grant Ritchie The postulates.. 1. Associated with any particle moving in a conservative field of force is a wave function,
More informationDOING PHYSICS WITH MATLAB
DOING PHYSICS WITH MATLAB ELECTROMAGNETISM USING THE FDTD METHOD [1D] Propagation of Electromagnetic Waves Matlab Download Director ft_3.m ft_sources.m Download and run the script ft_3.m. Carefull inspect
More informationApplications of Time Domain Vector Potential Formulation to 3-D Electromagnetic Problems
Applications of Time Domain Vector Potential Formulation to 3-D Electromagnetic Problems F. De Flaviis, M. G. Noro, R. E. Diaz, G. Franceschetti and N. G. Alexopoulos Department of Electrical Engineering
More informationClass 11: Thermal radiation
Class : Thermal radiation By analyzing the results from a number of eperiments, Planck found the energy density of the radiation emitted by a black body in wavelength interval (, d + was well described
More informationTHERE are two main aspects to the interaction of electromagnetic
334 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 46, NO. 3, MARCH 1998 FDTD Analysis of Wave Propagation in Nonlinear Absorbing and Gain Media Amit S. Nagra, Student Member, IEEE, and Robert A.
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 7
ECE 6341 Spring 216 Prof. David R. Jackson ECE Dept. Notes 7 1 Two-ayer Stripline Structure h 2 h 1 ε, µ r2 r2 ε, µ r1 r1 Goal: Derive a transcendental equation for the wavenumber k of the TM modes of
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 information1 Lectures 10 and 11: resonance cavities
1 1 Lectures 10 and 11: resonance cavities We now analyze cavities that consist of a waveguide of length h, terminated by perfectly conducting plates at both ends. The coordinate system is oriented such
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 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 informationComparison of Unit Cell Geometry for Bloch Wave Analysis in Two Dimensional Periodic Beam Structures
Clemson University TigerPrints All Theses Theses 8-018 Comparison of Unit Cell Geometry for Bloch Wave Analysis in Two Dimensional Periodic Beam Structures Likitha Marneni Clemson University, lmarnen@g.clemson.edu
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 informationCartesian Coordinates
Cartesian Coordinates Daniel S. Weile Department of Electrical and Computer Engineering University of Delaware ELEG 648 Cartesian Coordinates Outline Outline Separation of Variables Away from sources,
More informationTASK A. TRANSMISSION LINE AND DISCONTINUITIES
TASK A. TRANSMISSION LINE AND DISCONTINUITIES Task A. Transmission Line and Discontinuities... 1 A.I. TEM Transmission Line... A.I.1. Circuit Representation of a Uniform Transmission Line... A.I.. Time
More informationBragg waveguides. Leaky modes 1D Bragg reflector waveguides
Outline Bragg waveguides These are covered etensively in Nanophotonics. Let s look closely at the nature of the modes that are common to all such structures. Leaky modes 1D Bragg reflector waveguides 142
More informationM445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
More informationSystem Identification & Parameter Estimation
System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
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 informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
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 information3.5 Continuity of a Function One Sided Continuity Intermediate Value Theorem... 23
Chapter 3 Limit and Continuity Contents 3. Definition of Limit 3 3.2 Basic Limit Theorems 8 3.3 One sided Limit 4 3.4 Infinite Limit, Limit at infinity and Asymptotes 5 3.4. Infinite Limit and Vertical
More informationVibrating Strings and Heat Flow
Vibrating Strings and Heat Flow Consider an infinite vibrating string Assume that the -ais is the equilibrium position of the string and that the tension in the string at rest in equilibrium is τ Let u(,
More informationJapanese Journal of Applied Physics
Japanese Journal of Applied Physics Efficient finite-difference time-domain computation of resonant frequencies of rectangular Lamé mode resonators via a combination of the symmetry boundary conditions
More informationMultidimensional partitions of unity and Gaussian terrains
and Gaussian terrains Richard A. Bale, Jeff P. Grossman, Gary F. Margrave, and Michael P. Lamoureu ABSTRACT Partitions of unity play an important rôle as amplitude-preserving windows for nonstationary
More informationThe Schrödinger Equation in One Dimension
The Schrödinger Equation in One Dimension Introduction We have defined a comple wave function Ψ(, t) for a particle and interpreted it such that Ψ ( r, t d gives the probability that the particle is at
More informationComputational Modeling of Oscillatory Motion with Finite Difference Time Domain Method. Yue Wu
Computational Modeling of Oscillatory Motion with Finite Difference Time Domain Method by Yue Wu A Report Submitted in Partial Fulfillment of the Requirements for the Degree of MASTER OF ENGINEERING in
More informationTHE total-field/scattered-field (TFSF) boundary, first proposed
454 IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS, VOL. 5, 2006 Analytic Field Propagation TFSF Boundary for FDTD Problems Involving Planar Interfaces: Lossy Material and Evanescent Fields Kakhkhor Abdijalilov
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 informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More informationLecture 4.5 Schemes for Parabolic Type Equations
Lecture 4.5 Schemes for Parabolic Type Equations 1 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed
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 informationLight-induced spiral mass transport in azo-polymer films under vortex-beam illumination Supplementary Information
Light-induced spiral mass transport in azo-polymer films under vorte-beam illumination Supplementary Information Antonio Ambrosio a),1, Lorenzo Marrucci 1, Fabio Borbone, Antonio Roviello and Pasqualino
More informationFourier spectral methods for solving some nonlinear partial differential equations
Int. J. Open Problems Compt. Math., Vol., No., June ISSN 99-; Copyright ICSRS Publication, www.i-csrs.org Fourier spectral methods for solving some nonlinear partial differential equations Hany N. HASSAN
More informationLecture 16 February 25, 2016
MTH 262/CME 372: pplied Fourier nalysis and Winter 2016 Elements of Modern Signal Processing Lecture 16 February 25, 2016 Prof. Emmanuel Candes Scribe: Carlos. Sing-Long, Edited by E. Bates 1 Outline genda:
More informationIntroduction to the FDTD method
Introduction to the FDTD method Ilkka Laakso Department of Electrical Engineering and Automa8on Tfy- 99.3227 4.11.2015 Contents Principle of FDTD Deriva8on Basic proper8es Stability Dispersion Boundary
More information1. Waves and Particles 2. Interference of Waves 3. Wave Nature of Light
1. Waves and Particles 2. Interference of Waves 3. Wave Nature of Light 1. Double-Slit Eperiment reading: Chapter 22 2. Single-Slit Diffraction reading: Chapter 22 3. Diffraction Grating reading: Chapter
More informationCourse Secretary: Christine Berber O3.095, phone x-6351,
IMPRS: Ultrafast Source Technologies Franz X. Kärtner (Umit Demirbas) & Thorsten Uphues, Bldg. 99, O3.097 & Room 6/3 Email & phone: franz.kaertner@cfel.de, 040 8998 6350 thorsten.uphues@cfel.de, 040 8998
More informationEngineering Electromagnetics
Nathan Ida Engineering Electromagnetics With 821 Illustrations Springer Contents Preface vu Vector Algebra 1 1.1 Introduction 1 1.2 Scalars and Vectors 2 1.3 Products of Vectors 13 1.4 Definition of Fields
More informationComposition of and the Transformation of Functions
1 3 Specific Outcome Demonstrate an understanding of operations on, and compositions of, functions. Demonstrate an understanding of the effects of horizontal and vertical translations on the graphs of
More informationPre-Calculus and Trigonometry Capacity Matrix
Pre-Calculus and Capacity Matri Review Polynomials A1.1.4 A1.2.5 Add, subtract, multiply and simplify polynomials and rational epressions Solve polynomial equations and equations involving rational epressions
More informationLecture 2 Review of Maxwell s Equations and the EM Constitutive Parameters
Lecture 2 Review of Maxwell s Equations and the EM Constitutive Parameters Optional Reading: Steer Appendix D, or Pozar Section 1.2,1.6, or any text on Engineering Electromagnetics (e.g., Hayt/Buck) Time-domain
More information( ) ( ) ( ), ( 0 ), ( 0)
. (a Find the eigenvalues and eigenfunctions of problem: (b The differential equation ( ( ( =, - =, =. (8% - - = has one basis solution =. Show that the other basis solution on the interval - < < is =
More informationPHYS Concept Tests Fall 2009
PHYS 30 Concept Tests Fall 009 In classical mechanics, given the state (i.e. position and velocity) of a particle at a certain time instant, the state of the particle at a later time A) cannot be determined
More informationLeast Squares with Examples in Signal Processing 1. 2 Overdetermined equations. 1 Notation. The sum of squares of x is denoted by x 2 2, i.e.
Least Squares with Eamples in Signal Processing Ivan Selesnick March 7, 3 NYU-Poly These notes address (approimate) solutions to linear equations by least squares We deal with the easy case wherein the
More informationElectrical and optical properties of materials
Electrical and optical properties of materials John JL Morton Part 4: Mawell s Equations We have already used Mawell s equations for electromagnetism, and in many ways they are simply a reformulation (or
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationA Computer Controlled Microwave Spectrometer for Dielectric Relaxation Studies of Polar Molecules
384 A Computer Controlled Microave Spectrometer for Dielectric Relaation Studies of Polar Molecules Jai N. Dahiya*, Santaneel Ghosh and J. A. Roberts+ Department of Physics and Engineering Physics, Southeast
More information2. Higher-order Linear ODE s
2. Higher-order Linear ODE s 2A. Second-order Linear ODE s: General Properties 2A-1. On the right below is an abbreviated form of the ODE on the left: (*) y + p()y + q()y = r() Ly = r() ; where L is the
More informationPhysics 3323, Fall 2014 Problem Set 13 due Friday, Dec 5, 2014
Physics 333, Fall 014 Problem Set 13 due Friday, Dec 5, 014 Reading: Finish Griffiths Ch. 9, and 10..1, 10.3, and 11.1.1-1. Reflecting on polarizations Griffiths 9.15 (3rd ed.: 9.14). In writing (9.76)
More informationY m = y n e 2πi(m 1)(n 1)/N (1) Y m e 2πi(m 1)(n 1)/N (2) m=1
The Discrete Fourier Transform (Bretherton notes): 1 Definition Let y n, n = 1,..., N be a sequence of N possibly comple values. The discrete Fourier transform () of this sequence is the sequence Y m,
More information2014 Mathematics. Advanced Higher. Finalised Marking Instructions
0 Mathematics Advanced Higher Finalised ing Instructions Scottish Qualifications Authority 0 The information in this publication may be reproduced to support SQA qualifications only on a noncommercial
More informationUniversity of Ljubljana Faculty of Mathematics and Physics. Photons in a double quantum resonator
University of Ljubljana Faculty of Mathematics and Physics Photons in a double quantum resonator Author : Jernej Györköš 28010210 Mentor: Norma Mankoč Borštnik Date : 22.3.2006 Abstract The interaction
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 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 informationResearch on Ground Penetrating Radar Migration Imaging Technology
Sensors & Transducers, Vol. 80, Issue 0, October 04, pp. 5-55 Sensors & Transducers 04 by IFSA Publishing, S. L. http://www.sensorsportal.com Research on Ground Penetrating Radar Migration Imaging Technology
More informationThe Simple Harmonic Oscillator
The Simple Harmonic Oscillator Michael Fowler, University of Virginia Einstein s Solution of the Specific Heat Puzzle The simple harmonic oscillator, a nonrelativistic particle in a potential ½C, is a
More informationUNIT I ELECTROSTATIC FIELDS
UNIT I ELECTROSTATIC FIELDS 1) Define electric potential and potential difference. 2) Name few applications of gauss law in electrostatics. 3) State point form of Ohm s Law. 4) State Divergence Theorem.
More information1. The first observations of the small-scale ETG-mode drift wave turbulence performed in FT-2 research tokamak.
Correlative Enhanced Scattering in the upper hybrid resonance for study of micro and meso-scale wave phenomena in low-temperature and tokamak plasmas Evgeniy Gusakov Ioffe Institute 6 Politekhnicheskaya,
More information