Interconnect Modeling

Size: px
Start display at page:

Download "Interconnect Modeling"

Transcription

1 Interconnect Modelng Modelng of Interconnects Interconnect R, C and computaton Interconnect models umped RC model Dstrbuted crcut models Hgher-order waveform n dstrbuted RC trees Accuracy and fdelty Prepared by CK 1

2 Interconnect Resstance Computaton Resstance R R Mawell s Equaton of Dfferental Form Mawell s equaton for statc electrc feld: ( ε E) = ρ In Cartesan coordnate system, gradent operator s = v v n y ny z v n ε: permttvty of the feld regon E: strength of electrc feld ρ: charge densty of the feld regon z Prepared by CK

3 Posson s Equaton ( ε E) = ρ Base on E =, the Possson s equaton s ( ε ) = ρ ε: permttvty of the feld regon E: strength of electrc feld ρ: charge densty of the feld regon aplace s Equaton ( ε ) = ρ If assume ρ = 0 and homogenous delectrc, potental satsfes aplace s equaton ( ε ) = ε = 0 = s aplacan operator In Cartesan coordnate system, aplace s equaton s y = 0 z = 0 Prepared by CK 3

4 Eact Etracton of Resstance Solvng aplace s equaton, Ψ = 0 Relaaton method: Potental at every pont s the average potental of ts neghbors Consder square grds [ ( 1, y) ( 1, y) (, y 1) (, 1) ] / 4 (, y) = y Pont on an edge, mrror the potental about the edge, e.g., rght vertcal edge [ ( 1, y) (, y 1) (, 1) ] / 4 (, y) = y Repeatedly replace the potental of all ponts wth the average of ther neghbors Equpotentals and Current Flow Current flow potental Prepared by CK 4

5 Break nes for Varous Shapes Potental Dstrbuton for Parallel and Dagonal Contacts Prepared by CK 5

6 Equvalent Perpendcular Contacts l p = mn(, l ) = l p / r c l p = mn(, l ) = l p / 4 r c Resstance Computaton for Non-rectangular Shapes [Horowtz-Dutton, T-CAD July 83] R = R = R = 4 1 R = 1 All lnes are rectlnear or dagonal (45º) Prepared by CK 6

7 Calculatng the Resstance of a Conve Regon Resstance Computaton for Non-rectangular Shapes (cont d) Rato = Rato = 1 A 1 C 1 1 Rato = E Rato = B D 1 Rato = 1 Shape Rato Resstance A 1 1 A 5 5 B 1.5 B B.55 B 3.83 C C.3 C 3.66 C 4.97 D 1.5 D D.43 D 3.74 E E 1.8 E 3.33 E 4.71 Prepared by CK 7

8 Accuracy of Resstance Etracton Shape Rato Etracted Eact %Error Resstance Resstance A A B B B.55.6 B C C.3.5 C C D 1.5. D D D E E E E Eact resstances are obtaned by solvng aplace s equaton Current Spreadng from a Small Contact Prepared by CK 8

9 Typcal sheet resstance for conductors Sheet Resstance Ω/SQ. Materal Mn. Typcal Ma. Metal[Al] Slcdes Dffuson (n and p) Polyslcon hat s Capactance? Q -Q Smplest model: parallel-plate capactor It has two parallel plates and homogeneous delectrc between them The capactance s C = ε ε: permttvty of delectrc A: area of plate d: dstance between plates The capactance s the capacty to store charge charge at each plate s Q = CV one s postve, the other s negatve A d Prepared by CK 9

10 General Pcture m1 c1 c13 m m3 c3 For multple conductors of any shapes and materals, and n any delectrc, there s a capactance between any two conductors Mutual capactance between m1 and m s C 1 = q 1 /v q 1 s the charge of m1 v 1 = 0 and v 3 = 0 Capactance Matr Capactance s often wrtten as a symmetrc matr m1 c1 c13 m m3 c3 c11 -c1 -c13 C = -c1 c -c3 -c31 -c3 c33 c c ( j ) s the self-capactance for a conductor, m = j= 1 j e.g., c 11 = c 1 c 13 The charge s gven by e.g., q = c v c v c v 1 = c ( v1 v) c13( v1 v3) m mm m T Q = C ( V ) 3 Prepared by CK 10

11 Physcal Desgn Doman Conductors: metal wre, va, polyslcon, substrate Delectrcs: SO,... vctm vctm cross-secton Total cap for a wre delay, power Mutual cap between wres sgnal ntegrty top-vew 3D Model 3D model res wth fnte wdth, thckness and length Compute self-cap and mutual cap 3D model D model D model res wth fnte wdth and thckness, but nfnte length along the 3rd dmenson Compute the unt-length cap (also called cap coeffcent) Prepared by CK 11

12 Appromaton of 3D Structure by a D Model Pck a 3rd dmenson Chop 3D structure nto sectons wth dstngush profles 1 3 3rd dmenson 3rd dmenson Appromaton of 3D structure va D Model Solve each dstngush profle va D model two types of profles wth unt-length cap c 1 and c Type I: 1 Type II: 3 3rd dmenson Total cap s c 1 ( 1 3 )c Prepared by CK 1

13 Appromaton va Quas-3D Model D model effects of sdewalls along the 3rd dmenson a correcton capactance for each sdewall Classfcaton (orthogonal to 3D/D) Numercal method accurate any geometrc structures etremely epensve Formula-based method effcent lmted accuracy and geometrc structures nsght nto dependency on desgn parameters Table lookup accuracy -> numercal method effcency -> analytcal method More fleble handlng of geometrc structures than analytcal method Prepared by CK 13

14 Framework for Numercal Method ❶ Assume voltage [0, 0,, 1,, 0],.e., only conductor has unt voltage ➋ Compute charge q j for every conductor j ➌ Obtan mutual cap c j =q j, and self-cap (sum of mutual cap) ➍ Iterate through steps 1-3 usng dfferent voltage assgnments V m = vctm m1 c1 m c4 c3 m4 m3 C mm = 0 -c c1 c -c3 -c4 0 -c c4 0 0 Solutons to Conductor Charge Dfferental-based use Mawell s equaton n dfferental form Integral-equaton based use Mawell s equaton n ntegral form Prepared by CK 14

15 Capactance Etracton - Electrostatc Analyss For m conductors solve m potental problems for the conductor surface charges Each problem has one conductor at 1V, the rest at zero potental Two Vews of Capactance Can derve C 1 from C 1, C 1, and C Set up a system of charge neutralty equatons C 1 C C 1 1 V V Solve for V and q, then C 1 = q 1 C C 1 C 1 q = 0 q Prepared by CK 15

16 Prepared by CK 16 Volume Methods - Fnte-Dfferences/Fnte-Elements Solve aplace s equaton n the conductor eteror Appromate dervatves by fnte-dfferences Conductors Provde potental boundary condtons Voltage one on conductor 1, zero on conductor -D eample Conductor 1 Conductor m k l j 0 )) ( ) 0.5(( )) ( ) 0.5(( = j k j j k k l m l l m m Volume Methods Generate Sparse Matrces One equaton for each grd node In 3-D, each equaton nvolves at least 7 varables Solve by matr soluton methods Sparse Gaussan Elmnaton Incomplete Cholesky Conjugate-Gradent Method (ICCG) Multgrd methods

17 Two Dfferent Boundary Condtons Closed bo: Overestmates self C, underestmates couplng C Open Bo: Underestmates self C, overestmates couplng C Components of Interconnect Capactance Classfcaton based on profles of nteractng nterconnects Area component Due to overlappng area of two nterconnects on dfferent layers Frngng component Due to sde-walls of an nterconnect and the surface of an nterconnect on a dfferent layer ateral component Due to sde-walls of two nterconnects on the same layer Prepared by CK 17

18 Capactance: The Parallel Plate Model H Electrcal-feld lnes SO t o Substrate C nt = ε o t o S S S S C, wre = = S Parallel-Plate Capactor d (0) = V ( d) = 0 aplace s equaton s = y z = 0 ( ) = k k 1 (0) = k 1 0 k = V d) = k d k 0 ( 1 = ( ) = V d V Prepared by CK 18

19 Parallel-Plate Capactor ( ) = V d V unt-area charge = unt-area cap = d ε parallel-plate cap = ε ( ) = ε V ε A d d Typcal Capactance/Unt Area for 1µm CMOS Interconnect ayer Polyslcon to Substrate Metal1 to Substrate Metal to Substrate Metal3 to Substrate N Dffuson to Substrate (@ 0 Volt) P Dffuson to Substrate (@ 0 Volt) ff/µm ± ± ± ± ± ± 0.06 Prepared by CK 19

20 Frngng Capactance H H - H/ Frngng Component Area Component Frngng and Area Capactances H π C = εo t o to t t o o ln 1 H H H C = ε o to to 0.5 H 1.06 to 0.5 Prepared by CK 0

21 Typcal Frngng Capactance/Unt ength for 1µm CMOS Interconnect ayer Polyslcon to Substrate Metal1 to Substrate Metal to Substrate Metal3 to Substrate ff/µm ± ± ± ± Frngng-Feld Effect Prepared by CK 1

22 Interwre Capactance evel evel 1 SO Substrate Creates crosstalk Couplng Capactances Area component Frngng component C X frngng = c f 1 / 0 / 0 ( e e ) 1 0 measures how frnge capactance vares for ncremental length of the frngng surface ateral component C d 1 X lateral Prepared by CK

23 Interwre Capactance for 1µm CMOS Metal1 to Polyslcon Metal to Polyslcon Metal to Metal1 Area ff/µm Frngng ff/µm Impact of Interwre Capactance Prepared by CK 3

24 Sgnfcance of Couplng Capactance C/Ctotal Mn. Spacng Mn. Spacng Technology (u m) Capactance Crosstalk Assume C X = 5 ff For 55mm overlap of X and Y C C C XY, area XY, frnge XY V X = ff = 1.9 ff = ff C X = C C X = 0.35V XY 5V Prepared by CK 4

25 ma (mm) Couplng Nose 10% Vdd mn. Spacng % Vdd mn. Spacng Technology (u m) Prepared by CK 5

PART 8. Partial Differential Equations PDEs

PART 8. Partial Differential Equations PDEs he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef

More information

V. Electrostatics. Lecture 25: Diffuse double layer structure

V. Electrostatics. Lecture 25: Diffuse double layer structure V. Electrostatcs Lecture 5: Dffuse double layer structure MIT Student Last tme we showed that whenever λ D L the electrolyte has a quas-neutral bulk (or outer ) regon at the geometrcal scale L, where there

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming OPTIMIATION Introducton ngle Varable Unconstraned Optmsaton Multvarable Unconstraned Optmsaton Lnear Programmng Chapter Optmsaton /. Introducton In an engneerng analss, sometmes etremtes, ether mnmum or

More information

Electrical Circuits 2.1 INTRODUCTION CHAPTER

Electrical Circuits 2.1 INTRODUCTION CHAPTER CHAPTE Electrcal Crcuts. INTODUCTION In ths chapter, we brefly revew the three types of basc passve electrcal elements: resstor, nductor and capactor. esstance Elements: Ohm s Law: The voltage drop across

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Lecture #4 Capacitors and Inductors Energy Stored in C and L Equivalent Circuits Thevenin Norton

Lecture #4 Capacitors and Inductors Energy Stored in C and L Equivalent Circuits Thevenin Norton EES ntro. electroncs for S Sprng 003 Lecture : 0/03/03 A.R. Neureuther Verson Date 0/0/03 EES ntroducton to Electroncs for omputer Scence Andrew R. Neureuther Lecture # apactors and nductors Energy Stored

More information

Pattern Classification

Pattern Classification Pattern Classfcaton All materals n these sldes ere taken from Pattern Classfcaton (nd ed) by R. O. Duda, P. E. Hart and D. G. Stork, John Wley & Sons, 000 th the permsson of the authors and the publsher

More information

2.29 Numerical Fluid Mechanics

2.29 Numerical Fluid Mechanics REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit.

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit. Physcs 4B Solutons to Chapter 7 HW Chapter 7: Questons:, 8, 0 Problems:,,, 45, 48,,, 7, 9 Queston 7- (a) no (b) yes (c) all te Queston 7-8 0 μc Queston 7-0, c;, a;, d; 4, b Problem 7- (a) Let be the current

More information

Chapter 2: Electric Energy and Capacitance

Chapter 2: Electric Energy and Capacitance Chapter : Electrc Energy and Capactance Potental One goal of physcs s to dentfy basc forces n our world, such as the electrc force as studed n the prevous lectures. Expermentally, we dscovered that the

More information

Phys102 General Physics II

Phys102 General Physics II Electrc Potental/Energy Phys0 General Physcs II Electrc Potental Topcs Electrc potental energy and electrc potental Equpotental Surace Calculaton o potental rom eld Potental rom a pont charge Potental

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Field and Wave Electromagnetic. Chapter.4

Field and Wave Electromagnetic. Chapter.4 Fel an Wave Electromagnetc Chapter.4 Soluton of electrostatc Problems Posson s s an Laplace s Equatons D = ρ E = E = V D = ε E : Two funamental equatons for electrostatc problem Where, V s scalar electrc

More information

Numerical Solution to Laplace s Equation

Numerical Solution to Laplace s Equation Numercal Soluton to Laplace s Equaton Carleton Unversty Department of Electroncs ELEC 3105 Laboratory Exercse 1 July 8 015 PRE-LABORATORY EXERCISE You need to complete the pre-lab and have the TA sgn off

More information

8.022 (E&M) Lecture 4

8.022 (E&M) Lecture 4 Topcs: 8.0 (E&M) Lecture 4 More applcatons of vector calculus to electrostatcs: Laplacan: Posson and Laplace equaton url: concept and applcatons to electrostatcs Introducton to conductors Last tme Electrc

More information

Computational Electromagnetics in Antenna Analysis and Design

Computational Electromagnetics in Antenna Analysis and Design Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for real-lfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures.

More information

Airflow and Contaminant Simulation with CONTAM

Airflow and Contaminant Simulation with CONTAM Arflow and Contamnant Smulaton wth CONTAM George Walton, NIST CHAMPS Developers Workshop Syracuse Unversty June 19, 2006 Network Analogy Electrc Ppe, Duct & Ar Wre Ppe, Duct, or Openng Juncton Juncton

More information

Finite Wings Steady, incompressible flow

Finite Wings Steady, incompressible flow Steady, ncompressble flow Geometrc propertes of a wng - Fnte thckness much smaller than the span and the chord - Defnton of wng geometry: a) Planform (varaton of chord and sweep angle) b) Secton/Arfol

More information

Supplementary information. Coulomb oscillations in a gate-controlled few-layer graphene quantum dot

Supplementary information. Coulomb oscillations in a gate-controlled few-layer graphene quantum dot Supplementary nformaton Coulomb oscllatons n a gate-controlled few-layer graphene quantum dot Ypu Song * Haonan Xong Wentao Jang Hongy Zhang Xao Xue Cheng Ma Yuln Ma Luyan Sun Hayan Wang and Lumng Duan

More information

Physics 114 Exam 2 Fall 2014 Solutions. Name:

Physics 114 Exam 2 Fall 2014 Solutions. Name: Physcs 114 Exam Fall 014 Name: For gradng purposes (do not wrte here): Queston 1. 1... 3. 3. Problem Answer each of the followng questons. Ponts for each queston are ndcated n red. Unless otherwse ndcated,

More information

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13 CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,

More information

Introduction to circuit analysis. Classification of Materials

Introduction to circuit analysis. Classification of Materials Introducton to crcut analyss OUTLINE Electrcal quanttes Charge Current Voltage Power The deal basc crcut element Sgn conventons Current versus voltage (I-V) graph Readng: 1.2, 1.3,1.6 Lecture 2, Slde 1

More information

Physics 2102 Spring 2007 Lecture 10 Current and Resistance

Physics 2102 Spring 2007 Lecture 10 Current and Resistance esstance Is Futle! Physcs 0 Sprng 007 Jonathan Dowlng Physcs 0 Sprng 007 Lecture 0 Current and esstance Georg Smon Ohm (789-854) What are we gong to learn? A road map lectrc charge lectrc force on other

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Rotor Noise Modeling Kenneth S. Brentner Penn State University

Rotor Noise Modeling Kenneth S. Brentner Penn State University Rotor Nose Modelng Kenneth S. Brentner Penn State Unversty Joby Avaton S4 www.jobyavaton.com 2018 Kenneth S. Brentner. All rghts reserved. 5 th Transformatve Vertcal Flght Workshop, January 18-19, 2018

More information

Electrostatic Potential from Transmembrane Currents

Electrostatic Potential from Transmembrane Currents Electrostatc Potental from Transmembrane Currents Let s assume that the current densty j(r, t) s ohmc;.e., lnearly proportonal to the electrc feld E(r, t): j = σ c (r)e (1) wth conductvty σ c = σ c (r).

More information

Implementation of the Matrix Method

Implementation of the Matrix Method Computatonal Photoncs, Prof. Thomas Pertsch, Abbe School of Photoncs, FSU Jena Computatonal Photoncs Semnar 0 Implementaton of the Matr Method calculaton of the transfer matr calculaton of reflecton and

More information

1. Mean-Field Theory. 2. Bjerrum length

1. Mean-Field Theory. 2. Bjerrum length 1. Mean-Feld Theory Contnuum models lke the Posson-Nernst-Planck equatons are mean-feld approxmatons whch descrbe how dscrete ons are affected by the mean concentratons c and potental φ. Each on mgrates

More information

Grid Generation around a Cylinder by Complex Potential Functions

Grid Generation around a Cylinder by Complex Potential Functions Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around

More information

MAGNETISM MAGNETIC DIPOLES

MAGNETISM MAGNETIC DIPOLES MAGNETISM We now turn to magnetsm. Ths has actually been used for longer than electrcty. People were usng compasses to sal around the Medterranean Sea several hundred years BC. However t was not understood

More information

Numerical Solution to Laplace s Equation

Numerical Solution to Laplace s Equation Numercal Soluton to Laplace s Equaton Carleton Unversty Department of Electroncs ELEC 3105 Laboratory Exercse 1 January 019 PRE-LABORATORY EXERCISE You need to complete the pre-lab and have the TA sgn

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Basic Statistical Analysis and Yield Calculations

Basic Statistical Analysis and Yield Calculations October 17, 007 Basc Statstcal Analyss and Yeld Calculatons Dr. José Ernesto Rayas Sánchez 1 Outlne Sources of desgn-performance uncertanty Desgn and development processes Desgn for manufacturablty A general

More information

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16 STAT 39: MATHEMATICAL COMPUTATIONS I FALL 218 LECTURE 16 1 why teratve methods f we have a lnear system Ax = b where A s very, very large but s ether sparse or structured (eg, banded, Toepltz, banded plus

More information

Solution for singularly perturbed problems via cubic spline in tension

Solution for singularly perturbed problems via cubic spline in tension ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool

More information

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017)

Advanced Circuits Topics - Part 1 by Dr. Colton (Fall 2017) Advanced rcuts Topcs - Part by Dr. olton (Fall 07) Part : Some thngs you should already know from Physcs 0 and 45 These are all thngs that you should have learned n Physcs 0 and/or 45. Ths secton s organzed

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

Dr. Fritz Wilhelm, Physics 230 E:\Excel files\230 lecture\ch26 capacitance.docx 1 of 13 Last saved: 12/27/2008; 8:40 PM. Homework: See website.

Dr. Fritz Wilhelm, Physics 230 E:\Excel files\230 lecture\ch26 capacitance.docx 1 of 13 Last saved: 12/27/2008; 8:40 PM. Homework: See website. Dr. Frtz Wlhelm, Physcs 3 E:\Excel fles\3 lecture\ch6 capactance.docx of 3 Last saved: /7/8; 8:4 PM Homework: See webste. Table of ontents: h.6. Defnton of apactance, 6. alculatng apactance, 6.a Parallel

More information

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018 MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.

( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation. 3 Interpolaton {( y } Gven:,,,,,, [ ] Fnd: y for some Mn, Ma Polynomal Appromaton Theorem (Weerstrass Appromaton Theorem --- estence ε [ ab] f( P( , then there ests a polynomal

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Gouy-Chapman model (1910) The double layer is not as compact as in Helmholtz rigid layer.

Gouy-Chapman model (1910) The double layer is not as compact as in Helmholtz rigid layer. CHE465/865, 6-3, Lecture 1, 7 nd Sep., 6 Gouy-Chapman model (191) The double layer s not as compact as n Helmholtz rgd layer. Consder thermal motons of ons: Tendency to ncrease the entropy and make the

More information

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12 REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal

More information

Sections begin this week. Cancelled Sections: Th Labs begin this week. Attend your only second lab slot this week.

Sections begin this week. Cancelled Sections: Th Labs begin this week. Attend your only second lab slot this week. Announcements Sectons begn ths week Cancelled Sectons: Th 122. Labs begn ths week. Attend your only second lab slot ths week. Cancelled labs: ThF 25. Please check your Lab secton. Homework #1 onlne Due

More information

Aerodynamics. Finite Wings Lifting line theory Glauert s method

Aerodynamics. Finite Wings Lifting line theory Glauert s method α ( y) l Γ( y) r ( y) V c( y) β b 4 V Glauert s method b ( y) + r dy dγ y y dy Soluton procedure that transforms the lftng lne ntegro-dfferental equaton nto a system of algebrac equatons - Restrcted to

More information

Announcements. Lecture #2

Announcements. Lecture #2 Announcements Lectures wll be n 4 LeConte begnnng Frday 8/29 Addtonal dscusson TA Denns Chang (Sectons 101, 105) Offce hours: Mo 2-3 PM; Th 5-6 PM Lab sectons begn Tuesday 9/2 Read Experment #1 onlne Download

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Numerical Investigation of Electroosmotic Flow. in Convergent/Divergent Micronozzle

Numerical Investigation of Electroosmotic Flow. in Convergent/Divergent Micronozzle Appled Mathematcal Scences, Vol. 5, 2011, no. 27, 1317-1323 Numercal Investgaton of Electroosmotc Flow n Convergent/Dvergent Mcronozzle V. Gnanaraj, V. Mohan, B. Vellakannan Thagarajar College of Engneerng

More information

Interconnect Optimization for Deep-Submicron and Giga-Hertz ICs

Interconnect Optimization for Deep-Submicron and Giga-Hertz ICs Interconnect Optmzaton for Deep-Submcron and Gga-Hertz ICs Le He http://cadlab.cs.ucla.edu/~hele UCLA Computer Scence Department Los Angeles, CA 90095 Outlne Background and overvew LR-based STIS optmzaton

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

New Method for Solving Poisson Equation. on Irregular Domains

New Method for Solving Poisson Equation. on Irregular Domains Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad

More information

PHY2049 Exam 2 solutions Fall 2016 Solution:

PHY2049 Exam 2 solutions Fall 2016 Solution: PHY2049 Exam 2 solutons Fall 2016 General strategy: Fnd two resstors, one par at a tme, that are connected ether n SERIES or n PARALLEL; replace these two resstors wth one of an equvalent resstance. Now

More information

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m)

ˆ (0.10 m) E ( N m /C ) 36 ˆj ( j C m) 7.. = = 3 = 4 = 5. The electrc feld s constant everywhere between the plates. Ths s ndcated by the electrc feld vectors, whch are all the same length and n the same drecton. 7.5. Model: The dstances to

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245 Insttuto Tecnológco de Aeronáutca FIITE ELEMETS I Class notes AE-5 Insttuto Tecnológco de Aeronáutca 5. Isoparametrc Elements AE-5 Insttuto Tecnológco de Aeronáutca ISOPARAMETRIC ELEMETS Introducton What

More information

ME 501A Seminar in Engineering Analysis Page 1

ME 501A Seminar in Engineering Analysis Page 1 umercal Solutons of oundary-value Problems n Os ovember 7, 7 umercal Solutons of oundary- Value Problems n Os Larry aretto Mechancal ngneerng 5 Semnar n ngneerng nalyss ovember 7, 7 Outlne Revew stff equaton

More information

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures

Analytical Gradient Evaluation of Cost Functions in. General Field Solvers: A Novel Approach for. Optimization of Microwave Structures IMS 2 Workshop Analytcal Gradent Evaluaton of Cost Functons n General Feld Solvers: A Novel Approach for Optmzaton of Mcrowave Structures P. Harscher, S. Amar* and R. Vahldeck and J. Bornemann* Swss Federal

More information

Numerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I

Numerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I 5 9 Mechancal Engneerng -.30 ME Lab I ME.30 Mechancal Lab I Numercal Methods Volt Sne Seres.5 0.5 SIN(X) 0 3 7 5 9 33 37 4 45 49 53 57 6 65 69 73 77 8 85 89 93 97 0-0.5 Normalzed Squared Functon - 0.07

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Line Drawing and Clipping Week 1, Lecture 2

Line Drawing and Clipping Week 1, Lecture 2 CS 43 Computer Graphcs I Lne Drawng and Clppng Week, Lecture 2 Davd Breen, Wllam Regl and Maxm Peysakhov Geometrc and Intellgent Computng Laboratory Department of Computer Scence Drexel Unversty http://gcl.mcs.drexel.edu

More information

Parallel VLSI CAD Algorithms. Lecture 7 Iterative methods for IC thermal analysis Zhuo Feng

Parallel VLSI CAD Algorithms. Lecture 7 Iterative methods for IC thermal analysis Zhuo Feng Parallel VLSI CAD Algortms Lecture 7 Iteratve metods for IC termal analyss Zuo Feng 7. We ve taled about usng LU to solve a sparse lnear matr t problem Anoter way of dealng wt sparse matr problems s to

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Three-dimensional eddy current analysis by the boundary element method using vector potential

Three-dimensional eddy current analysis by the boundary element method using vector potential Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

The Discretization Process

The Discretization Process FMIA F Moukalled L Mangan M Darwsh An Advanced Introducton wth OpenFOAM and Matlab Ths textbook explores both the theoretcal foundaton of the Fnte Volume Method (FVM) and ts applcatons n Computatonal Flud

More information

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00 ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental

More information

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is. Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

A constant recursive convolution technique for frequency dependent scalar wave equation based FDTD algorithm

A constant recursive convolution technique for frequency dependent scalar wave equation based FDTD algorithm J Comput Electron (213) 12:752 756 DOI 1.17/s1825-13-479-2 A constant recursve convoluton technque for frequency dependent scalar wave equaton bed FDTD algorthm M. Burak Özakın Serkan Aksoy Publshed onlne:

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

CALCU THIRD EDITION. Monfy J. Strauss. Texas Tech University. Gerald L Bradley. Claremont McKenna College. Karl J. Smith. Santa Rosa Junior College

CALCU THIRD EDITION. Monfy J. Strauss. Texas Tech University. Gerald L Bradley. Claremont McKenna College. Karl J. Smith. Santa Rosa Junior College CALCU THIRD EDITION Monfy J. Strauss Texas Tech Unversty Gerald L Bradley Claremont McKenna College Karl J. Smth Santa Rosa Junor College Prentce Hall, Upper Saddle Rver, New Jersey 07458 Preface x 1 Functons

More information

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression 11 MACHINE APPLIED MACHINE LEARNING LEARNING MACHINE LEARNING Gaussan Mture Regresson 22 MACHINE APPLIED MACHINE LEARNING LEARNING Bref summary of last week s lecture 33 MACHINE APPLIED MACHINE LEARNING

More information

Plate Theories for Classical and Laminated plates Weak Formulation and Element Calculations

Plate Theories for Classical and Laminated plates Weak Formulation and Element Calculations Plate heores for Classcal and Lamnated plates Weak Formulaton and Element Calculatons PM Mohte Department of Aerospace Engneerng Indan Insttute of echnolog Kanpur EQIP School on Computatonal Methods n

More information

EXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017

EXAMPLES of THEORETICAL PROBLEMS in the COURSE MMV031 HEAT TRANSFER, version 2017 EXAMPLES of THEORETICAL PROBLEMS n the COURSE MMV03 HEAT TRANSFER, verson 207 a) What s eant by sotropc ateral? b) What s eant by hoogeneous ateral? 2 Defne the theral dffusvty and gve the unts for the

More information

Electrochemistry Thermodynamics

Electrochemistry Thermodynamics CHEM 51 Analytcal Electrochemstry Chapter Oct 5, 016 Electrochemstry Thermodynamcs Bo Zhang Department of Chemstry Unversty of Washngton Seattle, WA 98195 Former SEAC presdent Andy Ewng sellng T-shrts

More information

Mean Field / Variational Approximations

Mean Field / Variational Approximations Mean Feld / Varatonal Appromatons resented by Jose Nuñez 0/24/05 Outlne Introducton Mean Feld Appromaton Structured Mean Feld Weghted Mean Feld Varatonal Methods Introducton roblem: We have dstrbuton but

More information

THEOREMS OF QUANTUM MECHANICS

THEOREMS OF QUANTUM MECHANICS THEOREMS OF QUANTUM MECHANICS In order to develop methods to treat many-electron systems (atoms & molecules), many of the theorems of quantum mechancs are useful. Useful Notaton The matrx element A mn

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE.

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE. !! www.clutchprep.com CONCEPT: ELECTROMAGNETIC INDUCTION A col of wre wth a VOLTAGE across each end wll have a current n t - Wre doesn t HAVE to have voltage source, voltage can be INDUCED V Common ways

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF

More information

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD he Gaussan classfer Nuno Vasconcelos ECE Department, UCSD Bayesan decson theory recall that we have state of the world X observatons g decson functon L[g,y] loss of predctng y wth g Bayes decson rule s

More information

First Year Examination Department of Statistics, University of Florida

First Year Examination Department of Statistics, University of Florida Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve

More information

Finite Difference Method

Finite Difference Method 7/0/07 Instructor r. Ramond Rump (9) 747 698 rcrump@utep.edu EE 337 Computatonal Electromagnetcs (CEM) Lecture #0 Fnte erence Method Lecture 0 These notes ma contan coprghted materal obtaned under ar use

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors. SCALARS AND ECTORS All phscal uanttes n engneerng mechancs are measured usng ether scalars or vectors. Scalar. A scalar s an postve or negatve phscal uantt that can be completel specfed b ts magntude.

More information