ECE 497 JS Lecture - 14 Projects: FDTD & LVDS

Similar documents
ECE 598 JS Lecture 08 Lossy Transmission Lines

Transmission Characteristics of Overhead MV/PLC Channels

Module 4: General Formulation of Electric Circuit Theory

ECE 546 Lecture 02 Review of Electromagnetics

Design and Simulation of Dc-Dc Voltage Converters Using Matlab/Simulink

Calculation of Ramp Response of Lossy Transmission Lines Using Two-port Network Functions

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

ECE 497 JS Lecture - 13 Projects

Computational modeling techniques

ELECTROMAGNETIC CHARACTERISTICS OF CON- FORMAL DIPOLE ANTENNAS OVER A PEC SPHERE

Department of Electrical Engineering, University of Waterloo. Introduction

MODAL ANALYSIS OF A TWO-DIMENSIONAL DIELECTRIC GRATING SLAB EXCITED BY AN OBLIQUELY INCIDENT PLANE WAVE

Bicycle Generator Dump Load Control Circuit: An Op Amp Comparator with Hysteresis

Evaluation of Eigenmodes of Dielectric Waveguides by a Numerical Method Based on the BPM

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

the results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

An RWG Basis Function Based Near- to Far-Field Transformation for 2D Triangular-Grid Finite Difference Time Domain Method

TOPPER SAMPLE PAPER 2 Class XII- Physics

Laboratory #2: Introduction to Microstripline Transmission Lines, Reflection and Transmission Coefficients, and S-Parameters

RADIATED EM FIELDS FROM A ROTATING CURRENT- CARRYING CIRCULAR CYLINDER: 2-DIMENSIONAL NUMERICAL SIMULATION

POLARISATION VISUAL PHYSICS ONLINE. View video on polarisation of light

ZVS Boost Converter. (a) (b) Fig 6.29 (a) Quasi-resonant boost converter with M-type switch. (b) Equivalent circuit.

Transduction Based on Changes in the Energy Stored in an Electrical Field

Direct Monte Carlo Simulation of Time- Dependent Problems

EDA Engineering Design & Analysis Ltd

Power Flow in Electromagnetic Waves. The time-dependent power flow density of an electromagnetic wave is given by the instantaneous Poynting vector

Current/voltage-mode third order quadrature oscillator employing two multiple outputs CCIIs and grounded capacitors

ENSC Discrete Time Systems. Project Outline. Semester

Building to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.

Synchronous Motor V-Curves

A Novel Isolated Buck-Boost Converter

PHYS College Physics II Final Examination Review

ECE 2100 Circuit Analysis

Modelling of NOLM Demultiplexers Employing Optical Soliton Control Pulse

ERRATA for RF and Microwave Design: A Systems Approach, Second Edition, First Printing by Michael Steer, 2013.

Least Squares Optimal Filtering with Multirate Observations

King Fahd University of Petroleum and Minerals. Electrical Engineering Department EE 420. Fiber Optics Communication.

CBSE Board Class XII Physics Set 1 Board Paper 2008 (Solution)

Schedule. Time Varying electromagnetic fields (1 Week) 6.1 Overview 6.2 Faraday s law (6.2.1 only) 6.3 Maxwell s equations

ELECTROSTATIC FIELDS IN MATERIAL MEDIA

Determining the Accuracy of Modal Parameter Estimation Methods

( ) ( ) ( ) ( ) ( z) ( )

NGSS High School Physics Domain Model

Guaranteeing Reliability with Vibration Simulation and Testing. Dr. Nathan Blattau

Kinetics of Particles. Chapter 3

THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES

4) What is the magnitude of the net electric field at the center of the square?

Verification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance

(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f

ECEN620: Network Theory Broadband Circuit Design Fall 2012

Chapter 32. Maxwell s Equations and Electromagnetic Waves

Applying Kirchoff s law on the primary circuit. V = - e1 V+ e1 = 0 V.D. e.m.f. From the secondary circuit e2 = v2. K e. Equivalent circuit :

1.1 The main transmission network of Eskom The classical two generator model 11

Optical property of Few-Mode Fiber with Non-uniform Refractive Index for Cylindrical Vector Beam Generation

Micro and Smart Systems

Progress In Electromagnetics Research M, Vol. 9, 9 20, 2009

Propagation of High Speed Digital Signals in Printed Circuit Board Systems - Phase I. Michael B. Steer Jeff Kastan Mark S. Basel Sasan H.

ECEN620: Network Theory Broadband Circuit Design Fall 2014

MICROWAVE COMMUNICATIONS AND RADAR

Lecture 02 CSE 40547/60547 Computing at the Nanoscale

Lab 11 LRC Circuits, Damped Forced Harmonic Motion

8 th Grade Math: Pre-Algebra

Collective effects in Beam Dynamics

Formal Uncertainty Assessment in Aquarius Salinity Retrieval Algorithm

Three charges, all with a charge of 10 C are situated as shown (each grid line is separated by 1 meter).

A Novel Electro-thermal Simulation Approach to Power IGBT Modules for Automotive Traction Applications

Technical Bulletin. Generation Interconnection Procedures. Revisions to Cluster 4, Phase 1 Study Methodology

Enhancing Performance of MLP/RBF Neural Classifiers via an Multivariate Data Distribution Scheme

Convergence Aid Techniques for Compact Modelling with Verilog-A

ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES

ME 458: Electronic Packaging and Materials Electromagnetic part Lecture note

Math 302 Learning Objectives

David HORN and Irit OPHER. School of Physics and Astronomy. Raymond and Beverly Sackler Faculty of Exact Sciences

Springs in parallel. Two springs in parallel. Springs in series. Springs in series. Resonance. Forced vibrations and resonance. 2 C. 2 1/2 m.

Linearization of the Output of a Wheatstone Bridge for Single Active Sensor. Madhu Mohan N., Geetha T., Sankaran P. and Jagadeesh Kumar V.

ADP EARTH ELECTRODES FOR GROUNDING NEMP TYPE LIGHTNING. David P. Mi Hard Georgia Institute of Technology, Atlanta, Georgia

Simulation of Line Outage Distribution Factors (L.O.D.F) Calculation for N-Buses System

Design and Analysis of Gas Turbine Blade by Potential Flow Approach

COUPLED THERMO-MECHANICAL ANALYSES OF DYNAMICALLY LOADED RUBBER CYLINDERS 1

ENG2410 Digital Design Sequential Circuits: Part A

Investigation of a Single-Point Nonlinearity Indicator in One-Dimensional Propagation. 2 Theory

Series and Parallel Resonances

Ray tracing equations in transversely isotropic media Cosmin Macesanu and Faruq Akbar, Seimax Technologies, Inc.

I. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is

ecommons University of Dayton Monish Ranjan Chatterjee University of Dayton, Fathi H.A. Mohamed University of Dayton

1. Transformer A transformer is used to obtain the approximate output voltage of the power supply. The output of the transformer is still AC.

ScholarlyCommons. University of Pennsylvania. Andrea Alù University of Pennsylvania. Nader Engheta University of Pennsylvania,

ECE 546 Lecture 13 Scattering Parameters

PYROMETRIC METHOD OF TEMPERATURE MEASUREMENT WITH COMPENSATION FOR SOLAR RADIATION

Unit/Comments. %; Equal Output Loads %; Equal Output Loads Full Load

Supplementary Course Notes Adding and Subtracting AC Voltages and Currents

High penetration of renewable resources and the impact on power system stability. Dharshana Muthumuni

Chapter 23 Electromagnetic Waves Lecture 14

Increasing Heat Transfer in Microchannels with Surface Acoustic Waves*

GENESIS Structural Optimization for ANSYS Mechanical

Supplementary Course Notes Adding and Subtracting AC Voltages and Currents

EXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE

ECE 545 Project Deliverables

Two-Dimensional COMSOL Simulation of Heavy-Oil Recovery by Electromagnetic Heating

Transcription:

ECE 497 JS Lecture - 14 Prjects: FDTD & LVDS Spring 2004 Jse E. Schutt-Aine Electrical & Cmputer Engineering University f Illinis jse@emlab.uiuc.edu 1

ECE 497 JS - Prjects All prjects shuld be accmpanied with a shrt paper (3-5 pages) Paper surveys shuld be abut 10-15 pages. P1. Write a prgram that simulates transients n a unifrm lssless line. P2. Write a mment methd cde t calculate the capacitance per unit length f a single micrstrip line. P3. Write a prgram that predicts the TDR respnse f a device frm the measured s parameters. P4. Write an FDTD prgram t calculate the frequency dependence f a micrstrip line P5. Develp an IBIS mdel fr a CMOS differential amplifier P6. Write a single TL prgram that will accept IBIS mdels at its terminatins P7. TBD n pwer distributin P8. LVDS versus single-ended design P9. Paper survey n related subjects 2

MOTIVATION Mdel signal cupling and distrtin in high speed circuit Predict prpagatin characteristics beynd TEM Extract frequency dependence f circuit parameters Full wave simulatin is necessary 3

Full-Wave Methds E r = r B t H r = r J + D r t B r = 0 D r = ρ v FDTD: Discretize equatins and slve with apprpriate bundary cnditins 4

Finite Difference Time Dmain Yee Algrithm E n n 1 c x (i, j,k)= E t x + ε y H z n 1/ 2 n (i, j, k) H 1/2 z (i, j 1, k) c t ε z H y n 1/ 2 n 1/ (i, j, k) H 2 y (i, j,k 1) ( ) ( ) Ex Ey x z y Hz Ez Hx Ey Ex Ez Hy Hy n +1/ H 2 n 1/2 c x (i, j,k)= H t x µ y E z n (i, j + 1, k) E n z (i, j,k) ( ) + c t µ z E y n (i, j,k + 1) E n y (i, j,k) ( ) Hx Ex Hz Ez Ey Ex Ey 5

2D-FDTD E y E x H z y x E n x i + 1 2,j = E x n 1 i + 1 2,j + t ε y H z n 1/2 i + 1 2,j+ 1 2 H z n 1/2 i + 1 2,j 1 2 E n y i, j + 1 2 = E y n 1 i, j + 1 2 t ε x H z n 1/2 i + 1 2,j+ 1 2 H z n 1/2 i 1 2,j+ 1 2 n+1/2 H 1 z i + 2,j+ 1 2 = H z n 1/2 i + 1 2,j+ 1 2 + t µ y E x n i + 1 2,j+1 E x n i + 1 2,j t - µ x E n 1 y i +1, j + 2 E n 1 x i, j + 2 6

2D-PML Frmulatin Simulatin Medium ε E x t ε E y t µ H z t = H z y = H z x = E x y E y x y ε E x t ε E y t µ H z t PML Medium +σe x = H z y +σe y = H z x +σ*h z = E x y E y x N reflectin frm PML interface x σ ε = σ * µ 7

PML-FDTD Frmulatin Perfectly Matched Layer (PML) artificial lssy medium reflectinless absrptin f EM waves independent f frequency r angle f incidence characterized by electrical cnductivity and magnetic cnductivity σ ε = σ µ The abve relatin ensures that the wave impedance f the PML medium is matched t that f the adjacent physical medium 8

PML-FDTD Frmulatin Mdified Maxwell s equatins: and where µ H sxz t +σ x H sxz = E y x ε E si +σ t i E si = iî H E = E sx +E sy +E sz and H = H sx +H sy +H sz i = x,y,z The abve equatins cntain twelve scalar equatins with twelve split field unknwns. Fr example, µ H sxz t +σ x H sxz = E y x 9

PML-FDTD Frmulatin The FDTD implementatin f these scalar equatins n a Yee grid is straightfrward. where n H sxz i, j,k = α n 1 mh sxz i, j,k β m n 1 E y i +1,j,k E y n 1 i, j,k x α m = µ t σ x 2 β m = µ t +σ x 2 / 1 µ t +σ x 2 10

CIRCUIT MODEL EM fields are simulated by the FDTD, where the PML is used as ABC The fields at sample psitins are recrded The current and vltage n the line are caculated by the fllwing integratins i(z,t) = H dl c v(z,t) = E dl [ ] After FFT, get the vltage and current matrices V(z,ω), [ I(z,ω) ] Calculate the circuit parameters L(ω), C(ω) 11

CIRCUIT MODEL Quasi-TEM mde f prpagatin in MTL s cnsisting f M lssless lines L(ω), C(ω) are the M M per unit length inductance and capacitance parameter matrices, respectively. V(z,ω), I(z,ω) are the line vltage and current vectrs. The circuit parameter matrices can be caculated by L(ω) = 1 jω ( d dz [ V(z,ω) ]I(z,ω) [ ] 1 ) C(ω)= 1 jω ( d dz [ I(z,ω) ]V(z,ω) [ ] 1 ) [ ] [ I(z,ω) ] dz d V(z,ω)= jωl(ω)i(z,ω) dz d I(z,ω)= jωc(ω)v(z,ω) V(z,ω), are the line vltage and current vectr matrices (M excitatins) 12

Numerical Results Verificatin f the FDTD-PML cde The cnductivities in PML are chsen with parablic prfiles The fields are excited with a single frequency surce (18 GHz) Nrmalized currents 1 0.5 0-0.5-1 210 3 410 3 610 3 Time steps 13

Numerical Results Cupled lssless transmissin lines w s t z h y Parameters: w=s=0.3 mm, t=0.05 mm and h=0.25 mm 14

. Cupled Lines s/h=0.8 0.107 0.106 s/h=1.0 s/h=1.2-0.01-0.011-0.012 C 11 (nf/m) 0.105 0.104 0.103 C 12 (nf/m) -0.013-0.014-0.015-0.016 s/h=0.8 s/h=1. s/h=1.2 0.102 0 5 10 15 20 25 Frequency (GHz) w s w -0.017 0 5 10 15 20 25 Frequency (GHz) h Parameters: w=0.3mm, h=0.25mm. Dielectric cnstant is 4.5. s/h=0.8 L 11 (nh/m) 360 355 350 345 340 s/h=0.8 s/h=1.0 s/h=1.2 L 12 (nh/m) 105 100 95 90 85 80 75 s/h=1.0 s/h=1.2 335 0 5 10 15 20 25 Frequency (GHz) 70 65 0 5 10 15 20 25 Fr equency( GHz) 15

. EMI Simulatins 1 2 W 3 4 h w=1.2mm, h=1mm, ε r = 4.5 Prpagating signal (db) -60-70 -80-90 -100 1 2 3 4 Prpagating signal (db) -40-50 -60-70 -80 1 2 3 4-110 0 2 4 6 8 10 Frequency (GHz) -90 0 10 20 30 40 50 Frequency (GHz) 16

2D-FDTD fr Intercnnects MOTIVATION * Take advantage f single mde prpagatin * Reduce cmputatinal dmain * Frmulate FDTD prblem in transverse directin STRATEGY E(x,y,z) E(x, y) H(x,y,z) H(x, y) z jβ z e jβ zz Use βz as input and btain transverse field prfiles fr different values f βz 17

. 2D-FDTD W h Parameters: w=h=1mm. The effective dielectric cnstant is 9.8. 0.19 C 600 L Capacitance (nf/m) 0.185 0.18 0.175 0.17 0.165 2D-FDTD 3D-FDTD Inductance (nh/m) 550 500 450 400 2D-FDTD 3D-FDTD 0.16 0 5 10 15 20 25 Frequency (GHz) 350 0 5 10 15 20 25 Frequency (GHz) 18

FDTD/PML Summary FDTD-PML prvides an accurate methd fr the extractin f intercnnect characteristics. Radiatin effects can be simulated accurately using FDTD-PML. 2D-FDTD is an efficient methd fr extracting the frequency dependence f intercnnects. Implementatin f absrbing bundary cnditin is critical t the accuracy f the methd. 19

References [1] K. S. Yee, Numerical slutin f initial bundary value prblems invlving Maxwell s equatins in istrpic media, IEEE trans. Antennas Prpagat., vl 14, pp. 302-307, May 1966. [2] J. P. Berenger, " A perfectly matched layer fr the absrptin f electrmagnetic waves", J. Cmputatinal Physics, vl 144, pp. 185-200, Oct. 1994 [3] R. Mittra, W. D. Becker, and P. H. Harms, " A general purpse Maxwell slver fr the extractin f equivalent circuits f electric package cmpnent fr circuit simulatin." IEEE Trans. Circuits Syst. I, vl.39, pp964-973, Nv. 1992. [4] T. Dhaene, S. Criel, and D. D. Zutter, " Analysis and mdeling f cupled dispersive intercnnectin lines," IEEE Trans. MTT., vl40, pp2103-2105, Nv.1992. [5] J. Zha and Z. F. Li, " A time-dmain full-wave extractin methd f frequency-dependent equivalent circuit parameters f multicnductr intercnnectin lines", IEEE Trans. MTT, vl 45, pp23-31, Jan. 1997 [6] T. Daehne and D. De Zutter, " CAD-riented general circuit descriptin f unifrm cupled lssy dispersive waveguide structures, " IEEE Trans. MTT., vl. 40, pp1545-1559, July 1992. 20

Lw-Vltage Differential Signaling (LVDS) Definitin: Methd t cmmunicate data using a very lw vltage swing (abut 350mV) differentially ver tw PCB traces r a balanced cable Criteria fr high-perfrmance cmmunicatin - Bandwidth - Lw Pwer - Lw Nise Slutin exists fr very shrt and very lng distances; hwever fr bard-t-bard r bx-t-bx, this is a challenge 21

Why LVDS? 1. Differential transmissin is less susceptible t cmmn mde nise 2. Cnsequently they can use lwer vltage swings 3. In PC bard (micrstrip) dd-mde prpagatin is faster 22

LVDS Attributes fr EMI 1. Lw utput vltage swing 2. Slw edge rates 3. Odd-mde peratin (magnetic fields cancel) 4. Sft utput crner transitins 23

LVDS Driver and Receiver - Majrity f current flws acrss 100-hm resistr - Switching changes the directin f current - Lgic state determined by current directin 24

Differential Signaling Technlgies RS-422 PECL LVDS Differential Driver Output Vltage ±2 t ±5V ±600-1000 mv ±250-450 mv Receiver Input Threshld ±200 mv ±200-300mV ±100 mv Data Rate <30Mbps >400Mbps >400Mbps Supply Current Quad Driver (n lad, static) 60 ma (max) 32-65mA (max) 8.0mA Supply Current Quad Receiver (n lad, static) 23mA (max) 40mA (max) 15mA (max) Prpagatin Delay f Driver 11ns (max) 4.5ns (max) 1.7ns (max) Prpagatin Delay f Receiver 30ns (max) 7.0ns (max) 2.7ns (max) Pulse Skew (Driver r Receiver) N/A 500ps (max) 400ps (max) 25

LVDS Standard Maximum Switching Speed Depends n line driver Depends n selected media (type and length) LVDS Saves Pwer Pwer dissipated in lad is small LVDS devices are in CMOS=>lw static pwer Lwers system pwer thrugh current-mde Design Practices Matching is critical Preserve balance 26

Prject P8 - LVDS Design Perfrm simulatin experiments t cmpare LVDS and single-ended signaling with CMOS driver and receiver. Emphasize speed, pwer and nise issues t validate the use f LVDS. http://jsa6.ece.uiuc.edu/prjects/p6 http://www.natinal.cm/appinf/lvds 27