Methods for Ordinary Differential Equations. Jacob White
|
|
- Silvia Houston
- 6 years ago
- Views:
Transcription
1 Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy
2 Outine Initia Vaue probem exampes Signa propagation (circuits with capacitors). Space frame dynamics (struts and masses). Chemica reaction dynamics. Investigate the simpe finite-difference methods Forward-Euer, Backward-Euer, Trap Rue. Look at the approximations and agorithms Examine properties experimentay. Anayze Convergence for Forward-Euer
3 Appication Probems Signa Transmission in an Integrated Circuit Logic Gate Signa Wire Wire has resistance Wire and ground pane form a capacitor Ground Pane Logic Gate Meta Wires carry signas from gate to gate. How ong is the signa deayed?
4 Appication Probems Signa Transmission in an Integrated Circuit Circuit Mode resistor capacitor Cut the wire into sections. Mode wire resistance with resistors. Mode wire-pane capacitance with capacitors. Constructing the Mode
5 Appication Probems Osciations in a Space Frame What is the osciation ampitude?
6 Appication Probems Osciations in a Space Frame Simpified Structure Struts Bots Ground Load Exampe Simpified for Iustration
7 Appication Probems Osciations in a Space Frame Modeing with Struts, Joints and Point Masses Point Mass Strut Constructing the Mode Repace Meta Beams with Struts. Repace cargo with point mass.
8 Appication Probems Chemica Reaction Dynamics Crucibe Reagent Strange green stuff How fast is product produced? Does it expode?
9 i C1 Appication Probems v1 i v R2 2 i R1 i R3 C 1 R 2 R 1 R 3 C 2 Signa Transmission in an Integrated Circuit Constitutive Equations dvc ic = C dt 1 ir = vr R A 2x2 Exampe Noda Equations Yieds 2x2 System dv i C Conservation Laws i i i + + = C R R i i i 1 + C1 0 dt R1 R2 R 2 v1 0 C = 2 dv v 2 + dt R2 R3 R 2 + = C R R
10 Appication Probems Signa Transmission in an Integrated Circuit A 2x2 Exampe Let C = C = 1, R = R = 10, R = 1 dv C1 0 dt R1 R2 R 2 v1 0 C = dv v dt R2 R3 R Eigenvaues and Eigenvectors eigenvectors dx = dt A A = Eigenvaues x
11 Eigendecomposition: An Aside on Eigenanaysis () Consider an ODE: dx t = Axt (), x(0) = x dt λ 0 0 A E E E E E 1 = 1 2 n E n 0 0 λ n E 1 Change of variab e s : Ey( t) = xt ( ) yt ( ) = E xt ( ) dey() t Substituting: = AEy(), t Ey(0) = x0 dt λ dy() t 1 1 Mutipy by E : = E AEy() t = 0 0 dt 0 0 λ n 0 yt () 1
12 An Aside on Eigenanaysis Continued From ast side: dy() t dt λ1 0 0 = λ n yt () Decouped Equations! dyi () t λit Decouping: = λiyi() t yi() t = e y(0) dt () Steps for soving dx t = Axt (), x(0) = x0 dt 1) Determine E, λ 1 λ1t 2) Compute y(0) = E x e ) Compute yt ( ) = 0 0 y(0) λnt 0 0 e 4) xt () = Eyt ()
13 Appication Probems Signa Transmission in an Integrated Circuit A 2x2 Exampe v 1(0) = 1 v 2(0) = 0 Notice two time scae behavior v 1 and v 2 come together quicky (fast eigenmode). v 1 and v 2 decay to zero sowy (sow eigenmode).
14 Appication Probems Struts, Joints and point mass exampe A 2x2 Exampe y = y0 + u Define v as veocity (du/dt) to yied a 2x2 System f s f m Constitutive Equations 2 du fm = M dt 2 y y EA 0 c fs = EAc = u y0 y0 dv EAc M 0 dt 0 v y = du u 1 0 dt Conservation Law f s + f = m 0
15 Appication Probems Struts, Joints and point mass exampe A 2x2 Exampe Let M = 1, EA c = 1 y dv EAc M 0 dt 0 v y = du u 1 0 dt 0 Eigenvaues and Eigenvectors i A = i i 0 i i i eigenvectors Eigenvaues dx = x dt A
16 0.5 1 v (0) = 1 Appication Probems Struts, Joints and point mass exampe A 2x2 Exampe u (0) = Note the system has imaginary eigenvaues Persistent Osciation Veocity, v, peaks when dispacement, u, is zero.
17 Appication Probems Amount of reactant Chemica Reaction Exampe A 2x2 Exampe = R, the temperature = T dt T R dt = + More reactant causes temperature to rise, higher temperatures increases heat dissipation causing temperature to fa dr R 4T dt = + Higher temperatures raises reaction rates, increased reactant interferes with reaction and sows rate.
18 Appication Probems Chemica Reaction Exampe A 2x2 Exampe dt dt 1 1 T dx 1 1 = dr 4 1 R = dt 4 1 dt A Eigenvaues and Eigenvectors A = eigenvectors Eigenvaues x
19 Appication Probems R (0) = 0 Chemica Reaction Exampe A 2x2 Exampe T (0) = Note the system has a positive eigenvaue Soutions grow exponentiay with time.
20 First - Discretize Time Basic Concepts Second - Represent x(t) using vaues at t i ˆx 2 3 ˆx xˆ 4 ˆx x( t ) ˆx 1 0 t1 t2 t3 d xt dt Third - Approximate ( using ) the discrete tl t Approx. so n ˆ ˆ ˆ ˆ d x x x x Exampe: xt ( ) or dt t t 0 t t1 t2 tl t tl Exact so n xˆ 's = T
21 Basic Concepts Forward Euer Approximation x d sope = xt ( ) dt xt+ 1 xt sope = t t + 1 ( ) ( ) t t d xt ( 1) xt ( ) xt ( ) Axt + = ( ) dt t or xt ( ) xt ( ) + taxt ( ) + 1 = xt ( ) xt ( ) + taxt ( ) ( ) + 1
22 Basic Concepts Forward Euer Agorithm xt ( ) x = x(0) + tax 0 1 ˆ 1 ( ) x xt ( ˆ ˆ ˆ 2) x = x + tax xt ( ) xˆ = xˆ + taxˆ L L L 1 L 1 1 taxˆ tax(0) t1 t2 t 3 t
23 x Finite Difference d sope = xt ( + 1) dt xt+ 1 xt sope = t t + 1 ( ) ( ) t t Basic Concepts Backward Euer Approximation d xt ( 1) xt ( ) xt ( 1) Axt + + = ( + 1) dt t or xt ( ) xt ( ) + taxt ( ) = xt ( ) xt ( ) + taxt ( ) ( )
24 Basic Concepts Backward Euer Agorithm Sove with Gaussian Eimination xt ( ˆ ˆ 1) x = x(0) + tax 1 1 I ta xˆ = x 1 [ ] (0) xt ( ˆ ˆ 2) x = [ I ta] x L xt ( ) xˆ = [ I ta] xˆ L L 1 2 taxˆ 1 taxˆ x t1 t2 t
25 1 d d ( xt ( + 1 ) + xt ( )) 2 dt dt 1 = ( Ax ( t+ 1 ) + Ax ( t )) 2 xt ( + 1) xt ( ) t 1 xt ( + 1) xt ( ) + taxt ( ( + 1) + xt ( )) 2 Basic Concepts x Trapezoida Rue d sope = xt ( ) dt d sope = xt ( + 1) dt xt+ 1 xt sope = ( ) ( ) t t 1 1 = (( xt+ 1) taxt ()) (() xt + taxt ( + 1)) 2 2
26 Basic Concepts Trapezoida Rue Agorithm Sove with Gaussian Eimination t x( t1) xˆ = x(0) + Ax(0) + Axˆ 2 t t I A x = I + A x 2 2 ( ) ˆ (0) x 1 2 t t 1 xt ( 2) xˆ = I 2 A I+ 2 A xˆ 1 L ( ) ˆ t t xtl x = I A I+ A xˆ 2 2 L 1 t1 t2 t
27 Basic Concepts Numerica Integration View d xt xt xt t A d dt + 1 () = A xt ( ) ( + 1) = ( ) + x ( τ ) t t t + 1 Ax( τ ) dτ t 2 ( ) tax( t + 1) tax( t ) τ Ax( t) + Axt ( ) Trap BE FE t t + 1
28 Basic Concepts Summary Trap Rue, Forward-Euer, Backward-Euer Are a one-step methods xˆ is computed using ony xˆ, not xˆ, xˆ, etc. Forward-Euer is simpest No equation soution expicit method. Boxcar approximation to integra Backward-Euer is more expensive Equation soution each step impicit method Trapezoida Rue might be more accurate Equation soution each step impicit method Trapezoida approximation to integra
29 4 Numerica Experiments Unstabe Reaction t = 0.1 Backward-Euer Exact Soution Trap rue Forward-Euer FE and BE resuts have arger errors than Trap Rue, and the errors grow with time.
30 Numerica Experiments Unstabe Reaction-Error Pots x Forward Euer Backward Euer Trap Rue A methods have errors which grow exponentiay
31 M ax Finite Difference Backward-Euer Numerica Experiments Unstabe Reaction-Convergence 10-2 E rr 10-4 Trap rue Forward-Euer o r Timestep For FE and BE, Error t For Trap, Error ( t) 2
32 4 2 t = 0.1 Numerica Experiments Osciating Strut and Mass Forward-Euer Why does FE resut grow, BE resut decay and the Trap rue preserve osciations Trap rue Backward-Euer
33 Numerica Experiments Two timescae RC Circuit sma t Backward-Euer Computed Soution arge t With Backward-Euer it is easy to use sma timesteps for the fast dynamics and then switch to arge timesteps for the sow decay
34 Numerica Experiments Two timescae RC Circuit Forward-Euer Computed Soution The Forward-Euer is accurate for sma timesteps, but goes unstabe when the timestep is enarged
35 Convergence Numerica Experiments Summary Did the computed soution approach the exact soution? Why did the trap rue approach faster than BE or FE? Energy Preservation Why did BE produce a decaying osciation? Why did FE produce a growing osciation? Why did trap rue maintain osciation ampitude? Two timeconstant (stiff) probems Why did FE go unstabe when the timestep increased? We wi focus on convergence today
36 Convergence Anaysis Convergence Definition Definition: A finite-difference method for soving initia vaue probems on [0,T] is said to be convergent if given any A and any initia condition T 0, t ( ) max xˆ x t 0 as t 0 xˆ computed with t t xˆ computed with 2 x exact
37 Convergence Anaysis Order-p convergence Definition: A finite-difference method for soving initia vaue probems on [0,T] is said to be order p convergent if given any A and any initia condition max T 0, t for a t ess than a given ( ) ( ) xˆ x t C t t Forward- and Backward-Euer are order 1 convergent Trapezoida Rue is order 2 convergent 0 p
38 Convergence Anaysis Two Conditions for Convergence 1) Loca Condition: One step errors are sma (consistency) Typicay verified using Tayor Series 2) Goba Condition: The singe step errors do not grow too quicky (stabiity) A one-step methods are stabe in this sense.
39 Convergence Anaysis Consistency Definition Definition: A one-step method for soving initia vaue probems on an interva [0,T] is said to be consistent if for any A and any initia condition 1 ˆ ( ) x x t t 0 as t 0
40 Convergence Anaysis Consistency for Forward Euer Forward-Euer definition xˆ 1 = x 0 + tax 0 τ [ 0, t ] ( ) ( ) Expanding in t about zero yieds 2 2 dx 0 t d x x( t) = x( 0) + t + dt 2 dt ( ) ( ) ( ) d Noting that x (0) = Ax (0) and subtracting dt 2 2 τ ˆ 1 x x t ( ) t ( ) ( ) 2 d x dt 2 2 τ Proves the theorem if derivatives of x are bounded
41 Forward-Euer definition xˆ + = xˆ + taxˆ 1 Convergence Anaysis Convergence Anaysis for Forward Euer Expanding in t about t yieds (( 1) ) ( ) ( ) x + t = x t + tax t + e where e is the "one-step" error bounded by ( ) 2 e C t, where C = 0.5 max 2 d x τ [0, T ] 2 dt τ ( )
42 Convergence Anaysis Convergence Anaysis for Forward Euer Continued Subtracting the previous side equations ( ) (( 1) ) ( ) ˆ ( ) xˆ + x + t = I + ta x x t + e 1 Define the "Goba" error E + = I + ta E + e ( ) 1 ˆ ( ) E x x t Taking norms and using the bound on + 1 E I + ta E + C t ( ) ( ) ( ) 1 t A E C( t) e
43 Convergence Anaysis A hepfu bound on difference equations A emma bounding difference equation soutions ( ) If u 1 + ε u + b, u = 0, ε > 0 Then u e ε ε b To prove, first write u as a power series and sum 1 1 ( 1+ ε ) ( 1 ) j u + ε b = b 1 ( 1 ε ) j= 0 + j
44 Convergence Anaysis A hepfu bound on difference equations cont. ε To finish, note (1 + ε) e (1 + ε) Mapping the goba error equation to the emma + 1 E 1+ t A E + C ( t )2 ε b ( ε) ( ) j ( ε) e u b = b b 1 1+ j e ε ε ε ε ε
45 max C [ 0, L] A Convergence Anaysis Back to Forward Euer Convergence anaysis. Appying the emma and canceing terms 1 E 1+ t A E + C ( t )2 ε b Finay noting that t T, C( t) 2 AT E e t t A e t A
46 Convergence Anaysis Observations about the forward-euer anaysis. max C AT [ 0, L] E e t A forward-euer is order 1 convergent The bound grows exponentiay with time interva C is reated to the soution second derivative The bound grows exponentiay fast with norm(a).
47 12 Convergence Anaysis Exact and forward-euer(fe) Pots for Unstabe Reaction Texact Rexact RFE TFE Forward-Euer Errors appear to grow with time
48 E r r o r Rexact-RFE Time Convergence Anaysis forward-euer errors for soving reaction equation. Texact - TFE Note error grows exponentiay with time, as bound predicts
49 1 Convergence Anaysis Exact and forward-euer(fe) Pots for Circuit. 0.8 v1fe 0.6 v1exact 0.4 v2fe Forward-Euer Errors don t aways grow with time v2exact
50 E r r o r Finite Difference v1exact - v1fe Convergence Anaysis forward-euer errors for soving circuit equation v2exact-v2fe Time Error does not aways grow exponentiay with time! Bound is conservative
51 Summary Initia Vaue probem exampes Signa propagation (two time scaes). Space frame dynamics (osciator). Chemica reaction dynamics (unstabe system). Looked at the simpe finite-difference methods Forward-Euer, Backward-Euer, Trap Rue. Look at the approximations and agorithms Experiments generated many questions Anayzed Convergence for Forward-Euer Many more questions to answer, some next time
Introduction to Simulation - Lecture 13. Convergence of Multistep Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 13 Convergence of Mutistep Methods Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Loca truncation
More informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationMethods for Computing Periodic Steady-State Jacob White
Introduction to Simulation - Lecture 15 Methods for Computing Periodic Steady-State Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Periodic Steady-state problems Application
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More information1. Measurements and error calculus
EV 1 Measurements and error cacuus 11 Introduction The goa of this aboratory course is to introduce the notions of carrying out an experiment, acquiring and writing up the data, and finay anayzing the
More informationModule 22: Simple Harmonic Oscillation and Torque
Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationLECTURE 10. The world of pendula
LECTURE 0 The word of pendua For the next few ectures we are going to ook at the word of the pane penduum (Figure 0.). In a previous probem set we showed that we coud use the Euer- Lagrange method to derive
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 7
Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationChapter 5. Wave equation. 5.1 Physical derivation
Chapter 5 Wave equation In this chapter, we discuss the wave equation u tt a 2 u = f, (5.1) where a > is a constant. We wi discover that soutions of the wave equation behave in a different way comparing
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationDemonstration of Ohm s Law Electromotive force (EMF), internal resistance and potential difference Power and Energy Applications of Ohm s Law
Lesson 4 Demonstration of Ohm s Law Eectromotive force (EMF), interna resistance and potentia difference Power and Energy Appications of Ohm s Law esistors in Series and Parae Ces in series and Parae Kirchhoff
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationWave Equation Dirichlet Boundary Conditions
Wave Equation Dirichet Boundary Conditions u tt x, t = c u xx x, t, < x 1 u, t =, u, t = ux, = fx u t x, = gx Look for simpe soutions in the form ux, t = XxT t Substituting into 13 and dividing
More informationExpectation-Maximization for Estimating Parameters for a Mixture of Poissons
Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating
More information1 Heat Equation Dirichlet Boundary Conditions
Chapter 3 Heat Exampes in Rectanges Heat Equation Dirichet Boundary Conditions u t (x, t) = ku xx (x, t), < x (.) u(, t) =, u(, t) = u(x, ) = f(x). Separate Variabes Look for simpe soutions in the
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationInput-to-state stability for a class of Lurie systems
Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More information1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2
MATEC Web of Conferences 159, 01005 (018) IJCAET & ISAMPE 017 https://doi.org/10.1051/matecconf/01815901005 Dynamic Response of RC Cantiever Beam by Equivaent Singe Degree of Freedom Method on Eastic Anaysis
More informationA Fictitious Time Integration Method for a One-Dimensional Hyperbolic Boundary Value Problem
Journa o mathematics and computer science 14 (15) 87-96 A Fictitious ime Integration Method or a One-Dimensiona Hyperboic Boundary Vaue Probem Mir Saad Hashemi 1,*, Maryam Sariri 1 1 Department o Mathematics,
More informationNOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS
NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs
More informationFourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form
Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationChemical Kinetics Part 2
Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationIntroduction to Simulation - Lecture 16. Methods for Computing Periodic Steady-State - Part II Jacob White
Inroducion o Simuaion - ecure 6 Mehods for Compuing Periodic Seady-Sae - Par II Jacob Whie hanks o Deepak Ramaswamy, Micha Rewienski, and Karen Veroy Ouine hree Mehods so far ime inegraion uni seady-sae
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationChemical Kinetics Part 2. Chapter 16
Chemica Kinetics Part 2 Chapter 16 Integrated Rate Laws The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationInduction and Inductance
Induction and Inductance How we generate E by B, and the passive component inductor in a circuit. 1. A review of emf and the magnetic fux. 2. Faraday s Law of Induction 3. Lentz Law 4. Inductance and inductor
More informationLecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling
Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationTwo Kinds of Parabolic Equation algorithms in the Computational Electromagnetics
Avaiabe onine at www.sciencedirect.com Procedia Engineering 9 (0) 45 49 0 Internationa Workshop on Information and Eectronics Engineering (IWIEE) Two Kinds of Paraboic Equation agorithms in the Computationa
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationTHE ROLE OF ENERGY IMBALANCE MANAGEMENT ON POWER MARKET STABILITY
Proceedings of HICSS-31, Big Isand of Hawaii, January 6-9, 1998, Voume III, pp. 4-8. THE ROLE OF ENERGY IMBALANCE MANAGEMENT ON POWER MARKET STABILITY Fernando L. Avarado Department of Eectrica and C.
More informationNonlinear oscillation of a dielectric elastomer balloon
Research Artice Received: 28 May 2009 Revised: 31 Juy 2009 Accepted: 15 October 2009 Pubished onine in Wiey Interscience: 1 February 2010 (www.interscience.wiey.com) DOI 10.1002/pi.2767 Noninear osciation
More informationCSE464 Coupling Calculations via Odd/Even Modes Spring, David M. Zar (Based on notes by Fred U. Rosenberger)
SE464 ouping acuations ia Odd/Een Modes Spring, 009 aid M. Zar (Based on notes by Fred U. Rosenberger) dzar@wust.edu - aid M. Zar - 3/4/009 ouped Lines Here are two couped ines I I ( )/ I ( I )/ I ( Odd,
More informationWavelet Galerkin Solution for Boundary Value Problems
Internationa Journa of Engineering Research and Deveopment e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Voume, Issue 5 (May 24), PP.2-3 Waveet Gaerkin Soution for Boundary Vaue Probems D. Pate, M.K.
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationFourier series. Part - A
Fourier series Part - A 1.Define Dirichet s conditions Ans: A function defined in c x c + can be expanded as an infinite trigonometric series of the form a + a n cos nx n 1 + b n sin nx, provided i) f
More informationFinite element method for structural dynamic and stability analyses
Finite eement method for structura dynamic and stabiity anayses Modue-9 Structura stabiity anaysis Lecture-33 Dynamic anaysis of stabiity and anaysis of time varying systems Prof C S Manohar Department
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationMA 201: Partial Differential Equations Lecture - 11
MA 201: Partia Differentia Equations Lecture - 11 Heat Equation Heat conduction in a thin rod The IBVP under consideration consists of: The governing equation: u t = αu xx, (1) where α is the therma diffusivity.
More informationMultigrid Method for Elliptic Control Problems
J OHANNES KEPLER UNIVERSITÄT LINZ Netzwerk f ür Forschung, L ehre und Praxis Mutigrid Method for Eiptic Contro Probems MASTERARBEIT zur Erangung des akademischen Grades MASTER OF SCIENCE in der Studienrichtung
More informationVTU-NPTEL-NMEICT Project
MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid
More informationForces of Friction. through a viscous medium, there will be a resistance to the motion. and its environment
Forces of Friction When an object is in motion on a surface or through a viscous medium, there wi be a resistance to the motion This is due to the interactions between the object and its environment This
More informationLecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series
Lecture Notes for Math 251: ODE and PDE. Lecture 32: 1.2 Fourier Series Shawn D. Ryan Spring 212 Last Time: We studied the heat equation and the method of Separation of Variabes. We then used Separation
More information3.10 Implications of Redundancy
118 IB Structures 2008-9 3.10 Impications of Redundancy An important aspect of redundant structures is that it is possibe to have interna forces within the structure, with no externa oading being appied.
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More informationInternational Journal of Advance Engineering and Research Development
Scientific Journa of Impact Factor (SJIF): 4.4 Internationa Journa of Advance Engineering and Research Deveopment Voume 3, Issue 3, March -206 e-issn (O): 2348-4470 p-issn (P): 2348-6406 Study and comparison
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationSolution Set Seven. 1 Goldstein Components of Torque Along Principal Axes Components of Torque Along Cartesian Axes...
: Soution Set Seven Northwestern University, Cassica Mechanics Cassica Mechanics, Third Ed.- Godstein November 8, 25 Contents Godstein 5.8. 2. Components of Torque Aong Principa Axes.......................
More informationAST 418/518 Instrumentation and Statistics
AST 418/518 Instrumentation and Statistics Cass Website: http://ircamera.as.arizona.edu/astr_518 Cass Texts: Practica Statistics for Astronomers, J.V. Wa, and C.R. Jenkins, Second Edition. Measuring the
More informationIntroduction to Simulation - Lecture 9. Multidimensional Newton Methods. Jacob White
Introduction to Simulation - Lecture 9 Multidimensional Newton Methods Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Michal Rewienski, and Karen Veroy Outline Quick Review of -D Newton Convergence
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations
.615, MHD Theory of Fusion ystems Prof. Freidberg Lecture : The Moment Equations Botzmann-Maxwe Equations 1. Reca that the genera couped Botzmann-Maxwe equations can be written as f q + v + E + v B f =
More informationMode in Output Participation Factors for Linear Systems
2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong
More informationResearch Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method
Science and Technoogy of Nucear Instaations Voume 213, Artice ID 261327, 6 pages http://dx.doi.org/1.1155/213/261327 Research Artice Soution of Point Reactor Neutron Kinetics Equations with Temperature
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationMultilayer Kerceptron
Mutiayer Kerceptron Zotán Szabó, András Lőrincz Department of Information Systems, Facuty of Informatics Eötvös Loránd University Pázmány Péter sétány 1/C H-1117, Budapest, Hungary e-mai: szzoi@csetehu,
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationMAS 315 Waves 1 of 8 Answers to Examples Sheet 1. To solve the three problems, we use the methods of 1.3 (with necessary changes in notation).
MAS 35 Waves of 8 Answers to Exampes Sheet. From.) and.5), the genera soution of φ xx = c φ yy is φ = fx cy) + gx + cy). Put c = : the genera soution of φ xx = φ yy is therefore φ = fx y) + gx + y) ) To
More informationLaplace s Equation FEM Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy, Jaime Peraire and Tony Patera
Introduction to Simulation - Lecture 19 Laplace s Equation FEM Methods Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy, Jaime Peraire and Tony Patera Outline for Poisson Equation
More informationCopyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationMath 1600 Lecture 5, Section 2, 15 Sep 2014
1 of 6 Math 1600 Lecture 5, Section 2, 15 Sep 2014 Announcements: Continue reading Section 1.3 and aso the Exporation on cross products for next cass. Work through recommended homework questions. Quiz
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More information8 Digifl'.11 Cth:uits and devices
8 Digif'. Cth:uits and devices 8. Introduction In anaog eectronics, votage is a continuous variabe. This is usefu because most physica quantities we encounter are continuous: sound eves, ight intensity,
More informationIntroduction to Riemann Solvers
CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid
More informationT.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA
ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationEquilibrium of Heterogeneous Congestion Control Protocols
Equiibrium of Heterogeneous Congestion Contro Protocos Ao Tang Jiantao Wang Steven H. Low EAS Division, Caifornia Institute of Technoogy Mung Chiang EE Department, Princeton University Abstract When heterogeneous
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),
More informationWalking C H A P T E R 5
C H A P T E R 5 Waking q Practica egged ocomotion is one of the fundamenta unsoved probems in robotics. Many chaenges are in mechanica design - a waking robot must carry a of it s actuators and power,
More informationReliability: Theory & Applications No.3, September 2006
REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate
More informationInternational Journal of Mass Spectrometry
Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg
More informationTraffic data collection
Chapter 32 Traffic data coection 32.1 Overview Unike many other discipines of the engineering, the situations that are interesting to a traffic engineer cannot be reproduced in a aboratory. Even if road
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationComputational studies of discrete breathers. Sergej Flach MPIPKS Dresden January 2003
Computationa studies of discrete breathers Sergej Fach MPIPKS Dresden January 2003 CONTENT: 0. A bit on numerics of soving ODEs 1. How to observe breathers in simpe numerica runs 2. Obtaining breathers
More informationIndirect Optimal Control of Dynamical Systems
Computationa Mathematics and Mathematica Physics, Vo. 44, No. 3, 24, pp. 48 439. Transated from Zhurna Vychisite noi Matematiki i Matematicheskoi Fiziki, Vo. 44, No. 3, 24, pp. 444 466. Origina Russian
More informationHydrodynamic Instability of Liquid Films on Moving Fibers
Journa of Cooid and Interface Science 215, 381 396 (1999) Artice ID jcis.1999.6258, avaiabe onine at http://www.ideaibrary.com on Hydrodynamic Instabiity of Liquid Fims on Moving Fibers Konstantin G. Kornev*
More information6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7
6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the
More informationSteepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1
Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rues 1 R.J. Marks II, S. Oh, P. Arabshahi Λ, T.P. Caude, J.J. Choi, B.G. Song Λ Λ Dept. of Eectrica Engineering Boeing Computer Services University
More information