Roberto B. Armenta and Costas D. Sarris
|
|
- Alexandra Bates
- 6 years ago
- Views:
Transcription
1 A Method for Introducing Nonuniform Grids into the FDTD Solution of the Transmission-Line Equations Based on the Renormalization of the Per-Unit-Length Parameters Roberto B. Armenta and Costas D. Sarris Edward S. Rogers Sr. Department of Computer and Electrical Engineering, University of Toronto, ON M5S 3G4, Canada. Abstract: one of the challenging aspects of using the finite-difference time-domain (FDTD) method to solve the transmission-line (TL) equations is to choose a discretization grid with an adequate spatial resolution. In the standard formulation of the algorithm, a uniformly spaced discretization grid is used; nonetheless, different regions of the TL may require different levels of spatial resolution. To address this issue, this paper presents a nonuniform gridding scheme with an attractive feature: it preserves the simplicity of the standard FDTD update equations. It is shown that all the details of the nonuniform grid can be absorbed into effective per-unit-length parameters through a coordinate system transformation. By using the proposed step summing coordinate system transformation, a nonuniformly spaced grid can be mapped onto a uniformly spaced grid where a standard set of FDTD update equations can be applied. Keywords: finite-difference time-domain method and transmission lines.. Introduction When the finite-difference time-domain (FDTD) method is used to sove the transmission-line (TL) equations, a discretization grid with an adequate spatial resolution level must be carefully chosen. In the standard (2,p)-version of the algorithm, a uniform discretization grid is used. However, when dealing with position dependent per-unit-length parameters, different regions of the TL may require very different resolution levels. In such case, using a uniformly spaced discretization grid produces an unnecessary computational overhead. In the past, two different alternatives have been proposed to address this issue: the first one is to use a wavelet based multiresolution timedomain method (e.g., []-[2]); and the second one is to incorporate a nonuniform grid (e.g., [3]-[4]) into the original (2,2)-algorithm [5]. While many useful methods for introducing a nonuniform grid have been proposed before (see [6] for a summary), all of them introduce cumbersome modifications into the standard FDTD update equations. This paper proposes a method for introducing a nonuniform grid that preserves the simplicity of the standard (2,p)-algorithm on a uniform grid. It is shown that all the details of the nonuniform grid can be absorbed into effective or renormalized per-unit-length parameters through a one-dimensional coordinate system transformation. By using The notation (2,p) means second order accurate in time and p-th order accurate in space.
2 a coordinate system transformation, a nonuniformly spaced grid can be mapped onto a uniformly spaced grid where a (2,p)-set of standard FDTD update equations can be applied. The rest of this paper is divided into three sections: Section 2 briefly reviews the (2,p)-FDTD procedure for solving the TL equations; Section 3 introduces a one-dimensional coordinate system transformation that can be used to map a nonuniformly spaced grid onto a uniformly spaced one; and, in Section 4, these two are put together with a simple example. 2. The Standard (2,p)-FDTD Algorithm on a Uniform Grid For a general TL, the position and time dependent voltage and current functions V (z, t) and I(z, t) are related by V (z, t) I(z, t) = G(z) V (z, t) C(z) C(z) = R(z) I(z, t) L(z) L(z) I(z, t), () z V (z, t), (2) z where G(z) is the per-unit-length conductance, R(z) is the per-unit-length resistance, C(z) is the per-unit-length capacitance and L(z) is the per-unit-length inductance of the TL. In the standard (2,p)-FDTD algorithm, () and (2) are discretized on a uniformly spaced grid using second order accurate centered difference approximations in time and p-th order accurate centered difference approximations in space at points ((n+/2) t, i z) and (n t, (i+/2) z) for all n = 0,, 2,..., N T and i = 0,, 2,..., N Z. When a uniform discretization grid is used, the cell size z, which determines the spatial resolution of the grid, is usually selected so that λ min 20 z λ min 0 where λ min = z 0+L Z min {λ(z)} ; λ(z) = z=z 0 f max L(z)C(z). (3) In the above, f max is the maximum frequency excited by the source, and L Z is the z-axis length of the TL. Also, in the case of the (2,4)-algorithm, the time step t is picked according to λ min f max t/ z 6/7 (4) in order to satisfy stability requirements. Similar formulas apply to other high-order methods. Whenever z and t are picked according to (3) and (4), the uniform grid will be unnecessarily fine in regions where λ(z) is much larger than λ min. This well-known problem can be addressed by introducing a nonuniform discretization grid. Nonuniform grids are usually introduced in three steps. In step, a piecewise constant approximation of λ(z) is introduced. For this purpose, the TL s z-axis domain is divided into M nonoverlapping sections where the m-th section is contained within the interval [z m, z m ]. Then, λ(z) is approximated by a constant value λ m min = inside the m-th section. In step 2, all the sections where λ m min λ max = zm min z=z m {λ(z)} (5) M max m= {λm min} (6)
3 are put together into a single entity called the main grid region, which receives a coarse uniform discretization grid. Once the main grid region has been defined, finer grids are assigned to the remaining TL sections, which are called the refined grid regions. Finally, in step 3, a special set of FDTD update equations is generated to deal with the refined grid regions. The alternative nonuniform gridding approach proposed here consists of keeping the first two steps and replacing the last one with a coordinate transformation that maps the nonuniform grid generated in step 2 from (t, z)-space to (t, z)-space. If chosen as described in the next section, the coordinate transformation can be used to map a nonuniform grid in (t, z)-space to a uniform grid in (t, z)-space, where the problem can be solved using a standard (2,p)-set of update equations. 3. The Step Summing Coordinate System Transformation To define a coordinate transformation z z, we need to specify a function z = f( z). Substituting z = f( z) into () and (2) yields V ( z, t) I( z, t) = Ḡ( z) V ( z, t) C( z) C( z) = R( z) I( z, t) L( z) L( z) I( z, t), (7) z V ( z, t), (8) z Ḡ( z) = G( z)f ( z), C( z) = C( z)f ( z), R( z) = R( z)f ( z) and L( z) = L( z)f ( z). Observe that (7) and (8) have the same form as () and (2). This property which is known as metric invariance implies that (7) and (8) can be solved numerically in the (t, z)-domain by using a standard (2,p)-set of equations provided that the renormalized per-unit-length parameters are used in the computations. So, now the crucial question is: what choice of f( z) will achieve the desired mapping? First of all, f( z) must satisfy two key properties:. it must offer a one-to-one mapping z z, and 2. it must be continuous with continuous first derivatives. The first property is required so that the voltages and currents are uniquely defined in the (t, z)- domain. The second property is required so that no artificial discontinuities are introduced into the renormalized per-unit-length parameters (which could cause spurious reflections). The set of functions that satisfy the above two properties is very large; nonetheless, a simple coordinate transformation that has been found to be very versatile is presented next. For convenience, let us work with the inverse relation z = f (z) instead of z = f( z). Since a one-to-one mapping is required, this does not cause problems. Consider the coordinate transformation z = f (z) = z + ϕ(z), ϕ(z) = a m = λ max λ m min M a m d m ϕ m (z) where (9a) m= 0, d m = z m z m + 20 z; (9b)
4 ϕ m (z) = 2 + ( ) π 2 cos (z (z m + 0 z)), z m 0 z z z m + 0 z d m, z > z m + 0 z. (9c) Parameters λ max and λ m min were previously defined in (5) and (6). Also, z represents the cell size of the uniform grid in (t, z)-space. This parameter is usually chosen so that λ max 20 z λ max 0. (0) Once z has been chosen, the total number of cells N Z which must be the same in both spatial domains can be determined from { f (z M ) f } (z 0 ) N Z = round. () z To understand how the proposed coordinate system transformation achieves the desired mapping, observe that (9) yields z = z in the main grid region (i.e. where a m = 0). This implies that the main grid will have identical cell sizes in both spatial domains. On the other hand, in the refined grid regions (i.e. where a m > 0), the coordinate transformation adds a step function ϕ m (z) that effectively contracts distances (by a factor /(a m +)) as we go from the z-domain to the z-domain. This spatial contraction is what allocates a higher density of grid points in the refined grid regions even though a uniform grid is used in the (t, z)-space. When solving a TL problem using (9), three important practical issues must be kept in mind. First, in order to apply the transformation, we must know z = f( z) at a uniformly spaced set of points z i = i z for all i =, 2,..., N Z. While it is possible to analytically invert (9) to obtain z i = f( z i ), the easiest way to obtain the z i s is through an interpolation routine that computes them from the data of a plot of f (z) vs. z. Second, in order to renormalize the per-unit-length parameters, it is necessary to obtain the renormalization factor f ( z i ). Once the z i s have been obtained by interpolation, it is possible to obtain f ( z i ) from f ( z i ) = d dz f (z). (2) z=zi A closed form expression for the derivative of f (z) can be easily obtained from (9). The third and most important issue is the stability of the resulting FDTD scheme in the (t, z)-domain. As the per-unit-length parameters are renormalized, the stability criterion will be different from that of the conventional (2,p)-FDTD scheme. It can be shown that, for the (2,4)-set of update equations, the proposed time stepping procedure is stable as long as 7λ max f max t 6 z S, S = z 0+L Z min z=z 0 { } d dz f. (3) (z) Similar formulas will apply for higher order cases. Moreover, since renormalizing the per-unitlength parameters does not affect the dual structure of the TL equations, the proposed approach does not suffer from late time instability problems [7].
5 4. Transmssion-Line Example: a Bragg Reflector For demonstration purposes, consider the TL structure shown in Fig.. The TL is lossless with piecewise constant per-unit-length parameters given in Table. This structure is the TL equivalent of a Bragg reflector with two dielectric slabs. The goal in this example is to compute the reflection coefficient at z = cm assuming that absorbing boundaries are present at z = z 0 and z = z 5. This problem was solved in the (t, z)-domain by using a (2,4)-set of update equations and a uniform discretization grid with the following parameters: z =.268 mm = λ max /.82, λ max = c 0 /f max, f max = 20 GHz, t = ps, N Z = 424, N T = Both z and N Z were picked according to (0) and (). Observe that if we were to solve the problem using a uniform discretization grid in the (t, z)-domain with z = λ max /.82, the condition in (3) would be violated (since λ max = 5λ min ). To excite the TL, a Gaussian pulse of the form ( ( ) ) t 2 t0 f S (t) = exp, t 0 = 3T S, T S = 2f max T S was added at every time step to the voltage node located at z = cm. The voltage waveform was sampled at the node located at z = cm in order to compute the reflection coefficient from a discrete Fourier transform of the incident and reflected pulses. The amplitude of the computed reflection coefficient is shown in Fig. 2 together with the corresponding analytic solution cm.268 cm cm.268 cm cm TL Conductors z 0 = 0 z z 2 z 3 z 4 z 5 = L Z Figure : TL example under consideration. The TL has been divided into five sections (labeled with underlined numbers) with per-unit-length parameters given in Table. Table : Per-unit-length parameters of the TL example. z m L m C m µ 0 ε 0 2 µ 0 25ε 0 3 µ 0 ε 0 4 µ 0 25ε 0 5 µ 0 ε 0 m = Section Number L m = Per-Unit-Length Inductance C m = Per-Unit-Length Capacitance µ 0 = 4π0 7 H/m ε 0 = F/m c 0 = / µ 0 ε 0
6 Reflection Coefficient Amplitude (db) 5 0!0!20!30!40!50 FDTD Solution Analytic Solution! Frequency (GHz) (a) Error (db)!0!20!30!40!50!60!70!80!90!00!0! Frequency (GHz) (b) Figure 2: Plots of (a) the computed reflection coefficient amplitude and (b) the error in the estimate of the reflection coefficient amplitude. 5. Conclusions In summary, a method for introducing a nonuniform grid into the FDTD solution of the TL equations has been presented. A unique feature of this method is that it uses a coordinate transformation to introduce localized grid refinements. This allows the method to preserve the simplicity of the standard (2,p)-FDTD update equations. The presented coordinate transformation automatically selects an appropriate grid density for each region of the TL based on the TL s material properties. Finally, while the presented computed results were obtained using a (2,4)-set of update equations, the proposed approach easily incorporates other high-order FDTD schemes. References [] S. Grivet-Talocia and F. Canavero Wavelet-Based High-Order Adaptive Modeling of Lossy Interconnects, IEEE Trans. Electromagn. Compat., vol. 43, pp , 200. [2] S. Barmada and M. Raugi, Transient Numerical Solutions of Nonuniform MTL Equations with Nonlinear Loads by Wavelet Expansion in Time or Space Domain, IEEE Trans. Circuit Sys. I, vol. 47, pp.78-90, [3] P. Monk, Subgridding FDTD Schemes, ACES Jrnl., vol., pp , 996. [4] K. Xiao, D. J. Pommerenke, and J. L. Drewniak, A Three-Dimensional FDTD Subgridding Method with Separate Spatial and Temporal Subgridding Interfaces, in Proc. Intl. Symp. Electromagn. Compat., EMC 05, pp , [5] J. A. Roden, C. R. Paul, W. T. Smith and S. D. Gedney, Finite-Difference, Time-Domain Analysis of Lossy Transmission Lines, IEEE Trans. Electromagn. Compat., vol. 38, pp. 5-24, 996. [6] A. Taflove and S. Hagness (ed.), Computational Electrodynamics: the Finite-Difference Time- Domain Method. Boston: Artech House, 2005, ch.. [7] R. Schuhmann and T. Weiland, A Stable Interpolation Technique for FDTD on Non- Orhtogonal Grids, Intl. Jrnl. Numer. Model., vol., pp , 998.
A 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 informationPeriodic FDTD Characterization of Guiding and Radiation Properties of Negative Refractive Index Transmission Line Metamaterials
Periodic FDTD Characterization of Guiding and Radiation Properties of Negative Refractive Index Transmission Line Metamaterials Costas D. Sarris The Edward S. Rogers Sr. Department of Electrical and Computer
More informationlimit of the time-step decreases as more resolutions are added requires the use of an eective multitime-stepping algorithm, that will maintain the req
Invited to the Session: "Wavelet and TLM Modeling Techniques" organized by Dr. W. J. R. Hoefer in ACES 2000 COMPUTATIONAL OPTIMIZATION OF MRTD HAAR-BASED ADAPTIVE SCHEMES USED FOR THE DESIGN OF RF PACKAGING
More informationCombined FDTD/Macromodel Simulation of Interconnected Digital Devices
Combined FDTD/Macromodel Simulation of Interconnected Digital Devices S. Grivet-Talocia, I. S. Stievano, I. A. Maio, F. G. Canavero Dip. Elettronica, Politecnico di Torino, Torino, Italy (E-mail: grivet@polito.it)
More informationMulti-transmission Lines Loaded by Linear and Nonlinear Lumped Elements: FDTD Approach
Journal of Electrical Engineering 5 (2017) 67-73 doi: 10.17265/2328-2223/2017.02.002 D DAVID PUBLISHING Multi-transmission Lines Loaded by Linear and Nonlinear Lumped Elements: FDTD Approach Ismail ALAOUI
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 informationStatistical approach in complex circuits by a wavelet based Thevenin s theorem
Proceedings of the 11th WSEAS International Conference on CIRCUITS, Agios Nikolaos, Crete Island, Greece, July 23-25, 2007 185 Statistical approach in complex circuits by a wavelet based Thevenin s theorem
More informationNONUNIFORM multiconductor transmission
1 Nonuniform Multiconductor Transmission Line Analysis by a Two-Step Perturbation Technique Mykola Chernobryvko, Daniël De Zutter, Fellow, IEEE and Dries Vande Ginste, Senior Member, IEEE. Abstract A two-step
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 informationExtension of the FDTD Huygens Subgridding to Frequency Dependent Media
Noname manuscript No. (will be inserted by the editor) Extension of the FDTD Huygens Subgridding to Frequency Dependent Media Fumie Costen Jean-Pierre Bérenger eceived: date / Accepted: date Abstract A
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 informationSTUDY OF LOSS EFFECT OF TRANSMISSION LINES AND VALIDITY OF A SPICE MODEL IN ELECTROMAG- NETIC TOPOLOGY
Progress In Electromagnetics Research, PIER 90, 89 103, 2009 STUDY OF LOSS EFFECT OF TRANSMISSION LINES AND VALIDITY OF A SPICE MODEL IN ELECTROMAG- NETIC TOPOLOGY H. Xie, J. Wang, R. Fan, andy. Liu Department
More information/$ IEEE
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 55, NO. 9, SEPTEMBER 2008 937 Analytical Stability Condition of the Latency Insertion Method for Nonuniform GLC Circuits Subramanian N.
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 informationImplementation of a Transmission Line Model with the PEEC Method for Lightning Surge Analysis
Implementation of a Transmission Line Model with the PEEC Method for Lightning Surge Analysis PEERAWUT YUTTHAGOWITH Department of Electrical Engineering, Faculty of Engineering King Mongkut s Institute
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 informationB. H. Jung Department of Information and Communication Engineering Hoseo University Asan, Chungnam , Korea
Progress In Electromagnetics Research, PIER 77, 111 120, 2007 ANALYSIS OF TRANSIENT ELECTROMAGNETIC SCATTERING WITH PLANE WAVE INCIDENCE USING MOD-FDM B. H. Jung Department of Information and Communication
More informationEECS 117 Lecture 3: Transmission Line Junctions / Time Harmonic Excitation
EECS 117 Lecture 3: Transmission Line Junctions / Time Harmonic Excitation Prof. Niknejad University of California, Berkeley University of California, Berkeley EECS 117 Lecture 3 p. 1/23 Transmission Line
More informationSPICE MODELS FOR RADIATED AND CONDUCTED SUSCEPTIBILITY ANALYSES OF MULTICONDUCTOR SHIELDED CABLES
Progress In Electromagnetics Research, PIER 103, 241 257, 2010 SPICE MODELS FOR RADIATED AND CONDUCTED SUSCEPTIBILITY ANALYSES OF MULTICONDUCTOR SHIELDED CABLES H. Xie, J. Wang, R. Fan, and Y. Liu Department
More informationANALYSIS OF NON-UNIFORM TRANSMISSION LINE USING THE DISCRETIZATION THEORY
JAE, VOL. 5, NO., 3 ANALYSIS OF NON-UNIFORM TRANSMISSION LINE USING THE DISCRETIZATION THEORY L.Kumar*, H. Parathasarathy**, D.Upadhyay** Netaji Subhas Institute of Technology, Dwarka, University Of Delhi,
More informationECE357H1F ELECTROMAGNETIC FIELDS FINAL EXAM. 28 April Examiner: Prof. Sean V. Hum. Duration: hours
UNIVERSITY OF TORONTO FACULTY OF APPLIED SCIENCE AND ENGINEERING The Edward S. Rogers Sr. Department of Electrical and Computer Engineering ECE357H1F ELECTROMAGNETIC FIELDS FINAL EXAM 28 April 15 Examiner:
More informationCLOSED FORM SOLUTIONS FOR NONUNIFORM TRANSMISSION LINES
Progress In Electromagnetics Research B, Vol. 2, 243 258, 2008 CLOSED FORM SOLUTIONS FOR NONUNIFORM TRANSMISSION LINES M. Khalaj-Amirhosseini College of Electrical Engineering Iran University of Science
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 informationTransient Analysis of Interconnects by Means of Time-Domain Scattering Parameters
Transient Analysis of Interconnects by Means of Time-Domain Scattering Parameters Wojciech Bandurski, Poznań Univ. of Technology 60-965 Poznań, Piotrowo 3a, Poland, bandursk@zpe.iee.put.poznan.pl INTRODUCTION
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 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 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 informationTransmission Lines. Plane wave propagating in air Y unguided wave propagation. Transmission lines / waveguides Y. guided wave propagation
Transmission Lines Transmission lines and waveguides may be defined as devices used to guide energy from one point to another (from a source to a load). Transmission lines can consist of a set of conductors,
More informationEQUIVALENT CIRCUIT MODEL FOR ANALYSIS OF INHOMOGENEOUS GRATINGS
Progress In Electromagnetics Research, PIER 69, 21 34, 2007 EQUIVALENT CIRCUIT MODEL FOR ANALYSIS OF INHOMOGENEOUS GRATINGS M. Khalaj-Amirhosseini College of Electrical Engineering Iran University of Science
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 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 informationNon-Sinusoidal Waves on (Mostly Lossless)Transmission Lines
Non-Sinusoidal Waves on (Mostly Lossless)Transmission Lines Don Estreich Salazar 21C Adjunct Professor Engineering Science October 212 https://www.iol.unh.edu/services/testing/sas/tools.php 1 Outline of
More informationCircuit Representation of TL s A uniform TL may be modeled by the following circuit representation:
TRANSMISSION LINE THEORY (TEM Line) A uniform transmission line is defined as the one whose dimensions and electrical properties are identical at all planes transverse to the direction of propaation. Circuit
More informationOne-Dimensional Numerical Solution of the Maxwell-Minkowski Equations
Tamkang Journal of Science and Engineering, Vol. 12, No. 2, pp. 161168 (2009) 161 One-Dimensional Numerical Solution of the Maxwell-Minkowski Equations Mingtsu Ho 1 and Yao-Han Chen 2 1 Department of Electronic
More informationTransducer design simulation using finite element method
Transducer design simulation using finite element method Wenwu Cao Whitaker Center for Medical Ultrasonic Transducer Engineering Department of Mathematics and Materials Research Laboratory The Pennsylvania
More informationECE357H1S ELECTROMAGNETIC FIELDS TERM TEST 1. 8 February 2016, 19:00 20:00. Examiner: Prof. Sean V. Hum
UNIVERSITY OF TORONTO FACULTY OF APPLIED SCIENCE AND ENGINEERING The Edward S. Rogers Sr. Department of Electrical and Computer Engineering ECE57HS ELECTROMAGNETIC FIELDS TERM TEST 8 February 6, 9:00 :00
More informationI. INTRODUCTION II. NONLINEAR AUXILIARY DIFFERENTIAL EQUATION ADI-FDTD: FORMULATION
JOURNAL OF LIGHTWAVE TECHNOLOGY, VOL. 29, NO. 7, APRIL 1, 2011 1003 Time-Domain Modeling of Nonlinear Optical Structures With Extended Stability FDTD Schemes Dongying Li, Student Member, IEEE, and Costas
More informationModel-Order Reduction of High-Speed Interconnects: Challenges and Opportunities
Model-Order Reduction of High-Speed Interconnects: Challenges and Opportunities Michel Nakhla Carleton University Canada Model Reduction for Complex Dynamical Systems Berlin 2010 EMI Delay Crosstalk Reflection
More informationTwo-Layer Network Equivalent for Electromagnetic Transients
1328 IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 18, NO. 4, OCTOBER 2003 Two-Layer Network Equivalent for Electromagnetic Transients Mohamed Abdel-Rahman, Member, IEEE, Adam Semlyen, Life Fellow, IEEE, and
More informationTHE FINITE-DIFFERENCE TIME-DOMAIN (FDTD) METHOD
THE FINITE-DIFFERENCE TIME-DOMAIN (FDTD) METHOD Applied to the analysis of crosstalk between parallel striplines Edward Chan June 10, 1996 0. Introduction 1. The FDTD method 1.1 Maxwell s equations 1.2
More informationNumerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material
Chapter 30 Numerical Technique for Electromagnetic Field Computation Including High Contrast Composite Material Hiroshi Maeda Additional information is available at the end of the chapter http://dx.doi.org/10.5772/50555
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 informationSpectral Domain Analysis of Open Planar Transmission Lines
Mikrotalasna revija Novembar 4. Spectral Domain Analysis of Open Planar Transmission Lines Ján Zehentner, Jan Mrkvica, Jan Macháč Abstract The paper presents a new code calculating the basic characteristics
More informationModeling of Kerr non-linear photonic components with mode expansion
Modeling of Kerr non-linear photonic components with mode expansion Björn Maes (bjorn.maes@intec.ugent.be), Peter Bienstman and Roel Baets Department of Information Technology, Ghent University IMEC, St.-Pietersnieuwstraat
More informationINTEGRAL PML ABSORBING BOUNDARY CONDITIONS FOR THE HIGH-ORDER M24 FDTD ALGORITHM
Progress In Electromagnetics Research, PIER 76, 141 15, 007 INTEGRAL PML ABSORBING BOUNDARY CONDITIONS FOR THE HIGH-ORDER M4 FDTD ALGORITHM A. M. Shreim andm. F. Hadi Electrical Engineering Kuwait University
More informationThin-Slot/Thin-Layer Subcell FDTD Algorithms for EM Penetration through Apertures
EMG 23(2) #6315 Electromagnetics, 23:119 133, 2003 Copyright 2003 Taylor & Francis 0272-6343/03 $12.00 +.00 DOI: 10.1080/02726340390159469 Thin-Slot/Thin-Layer Subcell FDTD Algorithms for EM Penetration
More informationTransient analysis of spectrally asymmetric magnetic photonic crystals with ferromagnetic losses
PHYSICAL REVIEW B 7, 16507 006 Transient analysis of spectrally asymmetric magnetic photonic crystals with ferromagnetic losses K.-Y. Jung, B. Donderici, and F. L. Teixeira ElectroScience Laboratory and
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 informationEELE 3332 Electromagnetic II Chapter 11. Transmission Lines. Islamic University of Gaza Electrical Engineering Department Dr.
EEE 333 Electromagnetic II Chapter 11 Transmission ines Islamic University of Gaza Electrical Engineering Department Dr. Talal Skaik 1 1 11.1 Introduction Wave propagation in unbounded media is used in
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 informationProgress In Electromagnetics Research B, Vol. 56, , 2013
Progress In Electromagnetics Research B, Vol. 56, 409 427, 2013 A COMBINED FDTD/TLM TIME DOMAIN METHOD TO SOLVE EFFICIENTLY ELECTROMAGNETIC PROB- LEMS Nathanael Muot 1, Christophe Girard 1, Xavier Ferrieres
More informationEffects from the Thin Metallic Substrate Sandwiched in Planar Multilayer Microstrip Lines
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 115 Effects from the Thin Metallic Substrate Sandwiched in Planar Multilayer Microstrip Lines L. Zhang and J. M. Song Iowa
More informationPRACTICE NO. PD-AP-1309 PREFERRED PAGE 1 OF 5 RELIABILITY PRACTICES ANALYSIS OF RADIATED EMI FROM ESD EVENTS CAUSED BY SPACE CHARGING
PREFERRED PAGE 1 OF 5 RELIABILITY PRACTICES ANALYSIS OF RADIATED EMI FROM ESD EVENTS Practice: Modeling is utilized for the analysis of conducted and radiated electromagnetic interference (EMI) caused
More informationAnalyses of characteristics of spatial lightning electromagnetic field at different horizontal distances and observation point heights
Journal of Physics: Conference Series Analyses of characteristics of spatial lightning electromagnetic field at different horizontal distances and observation point heights To cite this article: Xiaojia
More informationand Ee = E ; 0 they are separated by a dielectric material having u = io-s S/m, µ, = µ, 0
602 CHAPTER 11 TRANSMISSION LINES 11.10 Two identical pulses each of magnitude 12 V and width 2 µs are incident at t = 0 on a lossless transmission line of length 400 m terminated with a load. If the two
More informationMICROSTRIP NONUNIFORM IMPEDANCE RESONATORS
Progress In Electromagnetics Research, PIER 67, 329 339, 2007 MICROSTRIP NONUNIFORM IMPEDANCE RESONATORS M. Khalaj-Amirhosseini College of Electrical Engineering Iran University of Science and Technology
More informationPulses in transmission lines
Pulses in transmission lines Physics 401, Fall 013 Eugene V. Colla Definition Distributed parameters networ Pulses in transmission line Wave equation and wave propagation eflections. esistive load Thévenin's
More informationKimmo Silvonen, Transmission lines, ver
Kimmo Silvonen, Transmission lines, ver. 13.10.2008 1 1 Basic Theory The increasing operating and clock frequencies require transmission line theory to be considered more and more often! 1.1 Some practical
More informationContents. Transmission Lines The Smith Chart Vector Network Analyser (VNA) ü structure ü calibration ü operation. Measurements
Contents Transmission Lines The Smith Chart Vector Network Analyser (VNA) ü structure ü calibration ü operation Measurements Göran Jönsson, EIT 2017-05-12 Vector Network Analysis 2 Waves on Lines If the
More informationStability and Passivity of the Super Node Algorithm for EM Modeling of IC s
Stability and Passivity of the Super Node Algorithm for EM Modeling of IC s M.V. Ugryumova and W.H.A. Schilders Abstract The super node algorithm performs model order reduction based on physical principles.
More informationLecture Outline 9/27/2017. EE 4347 Applied Electromagnetics. Topic 4a
9/7/17 Course Instructor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@utep.edu EE 4347 Applied Electromagnetics Topic 4a Transmission Lines Transmission These Lines notes may
More informationModeling and optimization of noise coupling in TSV-based 3D ICs
LETTER IEICE Electronics Express, Vol.11, No.20, 1 7 Modeling and optimization of noise coupling in TSV-based 3D ICs Yingbo Zhao, Yintang Yang, and Gang Dong a) School of Microelectronics, Xidian University,
More informationECE 546 Lecture 04 Resistance, Capacitance, Inductance
ECE 546 Lecture 04 Resistance, Capacitance, Inductance Spring 2018 Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jschutt@emlab.uiuc.edu ECE 546 Jose Schutt Aine 1 What is
More informationSørensen, Morten; Bonev, Ivan Bonev; Franek, Ondrej; Pedersen, Gert F.; Ebert, Hans
Downloaded from vbn.aau.dk on: januar 31, 219 Aalborg Universitet How to handle a Huygens' box inside an enclosure Sørensen, Morten; Bonev, Ivan Bonev; Franek, Ondrej; Pedersen, Gert F.; Ebert, Hans Published
More informationPower Absorption of Near Field of Elementary Radiators in Proximity of a Composite Layer
Power Absorption of Near Field of Elementary Radiators in Proximity of a Composite Layer M. Y. Koledintseva, P. C. Ravva, J. Y. Huang, and J. L. Drewniak University of Missouri-Rolla, USA M. Sabirov, V.
More informationSome observation on developments in lossy TLM
Some observation on developments in lossy TLM Donard de Cogan School of Information Systems, UEA, Norwich N4 7TJ (UK) Xiang Gui Department of Electrical and Computer Engineering, The University of Calgary,
More informationThe Finite-Difference Time-Domain Method for Electromagnetics with MATLAB Simulations
The Finite-Difference Time-Domain Method for Electromagnetics with MATLAB Simulations Atef Z. Elsherbeni and Veysel Demir SciTech Publishing, Inc Raleigh, NC scitechpublishing.com Contents Preface Author
More informationANALYSIS OF NONUNIFORM TRANSMISSION LINES USING THE EQUIVALENT SOURCES
Progress In Electromagnetics Research, PIER 7, 95 7, 27 ANALYSIS O NONUNIORM TRANSMISSION LINES USING THE EQUIVALENT SOURCES M. Khalaj-Amirhosseini College of Electrical Engineering Iran University of
More informationIntroduction to RF Design. RF Electronics Spring, 2016 Robert R. Krchnavek Rowan University
Introduction to RF Design RF Electronics Spring, 2016 Robert R. Krchnavek Rowan University Objectives Understand why RF design is different from lowfrequency design. Develop RF models of passive components.
More informationECE 546 Lecture 13 Scattering Parameters
ECE 546 Lecture 3 Scattering Parameters Spring 08 Jose E. Schutt-Aine Electrical & Computer Engineering University of Illinois jesa@illinois.edu ECE 546 Jose Schutt Aine Transfer Function Representation
More informationJoint time-frequency and finite-difference time-domain analysis of precursor fields in dispersive media
Joint time-frequency and finite-difference time-domain analysis of precursor fields in dispersive media Reza Safian, Costas D. Sarris, and Mohammad Mojahedi The Edward S. Rogers, Sr. Department of Electrical
More informationMicrowave Phase Shift Using Ferrite Filled Waveguide Below Cutoff
Microwave Phase Shift Using Ferrite Filled Waveguide Below Cutoff CHARLES R. BOYD, JR. Microwave Applications Group, Santa Maria, California, U. S. A. ABSTRACT Unlike conventional waveguides, lossless
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 informationECE 391 supplemental notes - #11. Adding a Lumped Series Element
ECE 391 supplemental notes - #11 Adding a umped Series Element Consider the following T-line circuit: Z R,1! Z,2! Z z in,1 = r in,1 + jx in,1 Z in,1 = z in,1 Z,1 z = Z Z,2 zin,2 = r in,2 + jx in,2 z,1
More informationTransmission Line Basics II - Class 6
Transmission Line Basics II - Class 6 Prerequisite Reading assignment: CH2 Acknowledgements: Intel Bus Boot Camp: Michael Leddige Agenda 2 The Transmission Line Concept Transmission line equivalent circuits
More informationTC 412 Microwave Communications. Lecture 6 Transmission lines problems and microstrip lines
TC 412 Microwave Communications Lecture 6 Transmission lines problems and microstrip lines RS 1 Review Input impedance for finite length line Quarter wavelength line Half wavelength line Smith chart A
More informationName. Section. Short Answer Questions. 1. (20 Pts) 2. (10 Pts) 3. (5 Pts) 4. (10 Pts) 5. (10 Pts) Regular Questions. 6. (25 Pts) 7.
Name Section Short Answer Questions 1. (20 Pts) 2. (10 Pts) 3. (5 Pts). (10 Pts) 5. (10 Pts) Regular Questions 6. (25 Pts) 7. (20 Pts) Notes: 1. Please read over all questions before you begin your work.
More informationScattering Parameters
Berkeley Scattering Parameters Prof. Ali M. Niknejad U.C. Berkeley Copyright c 2016 by Ali M. Niknejad September 7, 2017 1 / 57 Scattering Parameters 2 / 57 Scattering Matrix Voltages and currents are
More informationWIRE ANTENNA MODEL FOR TRANSIENT ANALYSIS OF SIMPLE GROUNDINGSYSTEMS, PART I: THE VERTICAL GROUNDING ELECTRODE
Progress In Electromagnetics Research, PIER 64, 149 166, 2006 WIRE ANTENNA MODEL FOR TRANSIENT ANALYSIS OF SIMPLE GROUNDINGSYSTEMS, PART I: THE VERTICAL GROUNDING ELECTRODE D. Poljak and V. Doric Department
More informationElectromagnetic Modeling and Signal Integrity Simulation of Power/Ground Networks in High Speed Digital Packages and Printed Circuit Boards
Electromagnetic Modeling and Signal Integrity Simulation of Power/Ground Networks in High Speed Digital Packages and Printed Circuit Boards Frank Y. Yuan Viewlogic Systems Group, Inc. 385 Del Norte Road
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 informationSTUDY ON THE STABILITY AND NUMERICAL ERROR OF THE FOUR-STAGES SPLIT-STEP FDTD METHOD INCLUDING LUMPED INDUCTORS
Progress In Electromagnetics Research B, Vol. 44, 117 135, 2012 STUDY ON THE STABILITY AND NUMERICAL ERROR OF THE FOUR-STAGES SPLIT-STEP FDTD METHOD INCLUDING LUMPED INDUCTORS Y.-D. Kong 1, Q.-X. Chu 1,
More informationContents. Transmission Lines The Smith Chart Vector Network Analyser (VNA) ü structure ü calibration ü operation. Measurements
Contents Transmission Lines The Smith Chart Vector Network Analyser (VNA) ü structure ü calibration ü operation Measurements Göran Jönsson, EIT 2015-04-27 Vector Network Analysis 2 Waves on Lines If the
More informationLarge omnidirectional band gaps in metallodielectric photonic crystals
PHYSICAL REVIEW B VOLUME, NUMBER 16 15 OCTOBER 1996-II Large omnidirectional band gaps in metallodielectric photonic crystals Shanhui Fan, Pierre R. Villeneuve, and J. D. Joannopoulos Department of Physics,
More informationProgress In Electromagnetics Research, Vol. 106, 1 16, 2010
Progress In Electromagnetics Research, Vol. 106, 1 16, 2010 TRANSIENT RESPONSES OF COAXIAL CABLES IN AN ELECTRICALLY LARGE CABIN WITH SLOTS AND WINDOWNS ILLUMINATED BY AN ELECTROMAG- NETIC PULSE J. Wang
More informationCapacitance of Dielectric Coated Cylinder of Finite Axial Length and Truncated Cone Isolated in Free Space
394 IEEE TRANSACTIONS ON ELECTROMAGNETIC COMPATIBILITY, VOL. 44, NO., MAY Capacitance of Dielectric Coated Cylinder of Finite Axial Length and Truncated Cone Isolated in Free Space S. B. Chakrabarty, Sukhendu
More informationElectrical Package Design TKK 2009 Lecture 2
Electrical Package Design TKK 2009 Lecture 2 James E. Morris Dept of Electrical & Computer Engineering Portland State University i Electrical Package Design Lecture topics A: Introduction CMOS; R, L, &
More informationALGORITHM FOR ACCURATE CAPACITANCE MATRIX MEASUREMENTS OF THE MULTICONDUCTOR STRUCTURE FOR VLSI INTERCONNECTIONS
ALGORITHM FOR ACCURATE CAPACITANCE MATRIX MEASUREMENTS OF THE MULTICONDUCTOR STRUCTURE FOR VLSI INTERCONNECTIONS Lech ZNAMIROWSKI Olgierd A. PALUSINSKI DEPARTMENT OF AUTOMATIC CONTROL, DEPARTMENT OF ELECTRICAL
More informationCURRENT- and voltage-driven common-mode coupling
18 IEEE TRANSACTIONS ON ELECTROMAGNETIC COMPATIBILITY, VOL 51, NO 1, FEBRUARY 2009 Cavity Model for the Slot Radiation of an Enclosure Excited by Printed Circuit Board Traces With Different Loads Christian
More informationEfficient Optimization of In-Package Decoupling Capacitors for I/O Power Integrity
1 Efficient Optimization of In-Package Decoupling Capacitors for I/O Power Integrity Jun Chen and Lei He Electrical Engineering Department, University of California, Los Angeles Abstract With high integration
More informationElectromagnetic Field Interaction with
Electromagnetic Field Interaction with Transmission Lines From classical theory to HF radiation effects Edited by F Rachidi sc S Tkachenko WITPRESS Southampton, Boston Contents Preface xv PART I: CLASSICAL
More informationTransmission Lines. Author: Michael Leddige
Transmission Lines Author: Michael Leddige 1 Contents PCB Transmission line structures Equivalent Circuits and Key Parameters Lossless Transmission Line Analysis Driving Reflections Systems Reactive Elements
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 informationIdentification and Classification of High Impedance Faults using Wavelet Multiresolution Analysis
92 NATIONAL POWER SYSTEMS CONFERENCE, NPSC 2002 Identification Classification of High Impedance Faults using Wavelet Multiresolution Analysis D. Cha N. K. Kishore A. K. Sinha Abstract: This paper presents
More informationTransmission Lines in the Frequency Domain
Berkeley Transmission Lines in the Frequency Domain Prof. Ali M. Niknejad U.C. Berkeley Copyright c 2016 by Ali M. Niknejad August 30, 2017 1 / 38 Why Sinusoidal Steady-State? 2 / 38 Time Harmonic Steady-State
More informationElimination of Conductor Foil Roughness Effects in Characterization of Dielectric Properties of Printed Circuit Boards 14 TH1
Elimination of Conductor Foil Roughness Effects in Characterization of Dielectric Properties of Printed Circuit Boards 14 TH1 Marina Koledintseva, Aleksei Rakov, Alexei Koledintsev, James Drewniak (Missouri
More informationCONTRIBUTION TO CALCULATING THE IMPEDANCE OF GROUNDING ELECTRODES USING CIRCUIT EQUIVALENTS. Andrijana Kuhar, Leonid Grcev
FACTA UNIVERSITATIS Series: Electronics and Energetics Vol. 29, N o 4, December 2016, pp. 721-732 DOI: 10.2298/FUEE1604721K CONTRIBUTION TO CALCULATING THE IMPEDANCE OF GROUNDING ELECTRODES USING CIRCUIT
More informationTHE PROPOSAL to use multiple-input multiple-output
IEEE TRANSACTIONS ON ELECTROMAGNETIC COMPATIBILITY, VOL 49, NO 4, NOVEMBER 2007 901 Modeling the Terminal Response of a Bundle of Twisted-Wire Pairs Excited by a Plane Wave Roberto B Armenta, Member, IEEE,
More informationPhysics-Based CAD Models for the Analysis of Vias in Parallel-Plate Environments
IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES, VOL. 49, NO. 10, OCTOBER 2001 1697 Physics-Based CAD Models for the Analysis of Vias in Parallel-Plate Environments Ramesh Abhari, Student Member,
More informationANTENNAS and MICROWAVES ENGINEERING (650427)
Philadelphia University Faculty of Engineering Communication and Electronics Engineering ANTENNAS and MICROWAVES ENGINEERING (65427) Part 2 Dr. Omar R Daoud 1 General Considerations It is a two-port network
More information