CS205b/CME306. Lecture 4. x v. + t
|
|
- Suzanna Marshall
- 6 years ago
- Views:
Transcription
1 CS05b/CME306 Lecture 4 Time Integration We now consider seeral popular approaches for integrating an ODE Forward Euler Forward Euler eolution takes on the form n+ = n + Because forward Euler is unstable when the system contains pure imaginary eigenalues, it is unstable without damping In particular, a < 0 Forward Euler is stable when λ + < A time step restriction of λ is required As problems get stiffer, ma λ becomes larger, and consequently the time step that can be taken becomes smaller The stability region for forward Euler is shown in Figure a If the system has eigenalues near the imaginary ais, we would be much better off using using a different method such as one of the higher order Runge-Kutta methods Runge-Kutta Methods There are a number of Runge-Kutta methods in common use The Runge-Kutta methods were discussed in more detail in CS05A, so they are only highlighted here The second order Runge- Kutta method is technically linearly unstable without damping since its stability region does not include the imaginary ais Roundoff with a < 0 may be sufficient, and nonlinear effects can sometimes rectify the problem The third order Runge-Kutta method is stable without damping, since its stability region includes part of the imaginary ais It is suitable for system with a = 0 Stability regions for arious Runge-Kutta methods are shown in Figure 3 Backward Euler The need for a time step restriction is true of all eplicit methods By contract, implicit methods often hae a much better time step restriction that makes them more suitable een though they require a linear system sole and are thus much more epensie per step Backward Euler is the simplest of the implicit methods Its stability region Figure a looks ery different from those of the eplicit integrators It takes the form n+ = n + a a n n+
2 a Forward Euler First Order Runge-Kutta b Second Order Runge-Kutta c Third Order Runge-Kutta d Fourth Order Runge-Kutta Figure : Stability regions for Runge-Kutta The scheme is stable with time step size if λ, a comple number in general, is located within the shaded region of the comple plane Note that the third and fourth order schemes hae stability regions that encompass the imaginary ais, making them stable on undamped problems The first and second order schemes do not hae this property
3 a Backward Euler b Trapezoid Rule Figure : Stability regions for two implicit time integration schemes The scheme is stable with time step size if λ, a comple number in general, is located within the shaded region of the comple plane Note that the new alues being obtained occur on both sides of the equation, so a sole is typically necessary at each time step Backward Euler has no time step restriction for stability, but large time steps do result in poor accuracy The sole required for backward Euler can be relatiely easy in the case of linear ODEs For eample, if the only force is drag, F = k and a = k/m Then, n+ = n k n+ /m so that n+ = + k/m n and n+ = n + n+ = n + + k/m n are relatiely easily obtained Note that this damps elocity each time step, and it damps more with larger time steps This is in contrast to forward Euler, where n+ = k/m n Damping here occurs if < m/k = ma λ but amplification occurs for larger time steps More generally, backward Euler requires soling a linear system If the ODE is nonlinear, it will typically be more difficult and epensie to sole and requires that the nonlinear system first be linearized 4 Trapezoid Rule Forward Euler eolution takes on the form n+ = n + a n + a n+ Because forward Euler is unstable when the system contains pure imaginary eigenalues, it is unstable without damping In particular, a < 0 Forward Euler is stable when λ + < A time step restriction of λ is required As problems get stiffer, ma λ becomes larger, and consequently the time step that can be taken becomes smaller The stability region for forward 3
4 = e t = 4 = = = Figure 3: Effects of large time step sizes on numerical simulation of y = y with backward Euler Euler is shown in Figure b If the system has eigenalues near the imaginary ais, we would be much better off using using a different method such as one of the higher order Runge-Kutta methods 5 Newmark Method Newmark methods are the most famous of the multialue methods in computation solids Multialue methods are those that take adantage of the fact that higher order deriaties can be computed at lower accuracy to be more efficient The Newmark method is a two-parameter family of methods: n+ = n + n + βan + βa n+ n+ = n + γa n + γa n+ The specific method depends on the choice of β and γ chosen Constant acceleration equations, β = γ = 0 Implicit acceleration, β = / and γ = n+ = n + n + an n+ = n + a n n+ = n + n + an+ n+ = n + γa n+ The elocity update is the first order backward Euler update The position update can be written as n+ = n + n + n+, which is just the second order accurate trapezoid rule This choice is still only first order accurate oerall There is a theorem that states that second order accuracy is obtained if and only if γ = / The following two methods are second order Newmark methods 4
5 Trapezoid rule, β = /4 and γ = / n+ = n + n + a n + a n+ n+ = n + an + a n+ This method gets its name since the position equation can also be rewritten as a trapezoid rule update as n+ = n + n + n+ Central differencing, β = 0 and γ = / n+ = n + n + an n+ = n + an + a n+ The name of this method comes from the ability of its updates to be epressed as the central differences n+ n = n n+ n + n = a n These equations can be obtained by n+ + n = n+ n n+ n n+ n n+ n = n+ n + n n+ n + n n+ n + n n+ n + n n + n + an = n + an + n + an + = n + n + an + a n = n n + n + an = n + an n an = n n + an a n = n + an + a n = a n n + n + an n + n + an n + an a n 6 Staggering the Position and Velocity Define elocity halfway between time steps so that the update n+/ = n+ n 5
6 is second order accurate Note that n is still defined at the normal times We can then define n by aeraging nearby elocities as n = n/ + n+/ = n+ n, which is also central differencing for elocity The update formula for positions is then just Acceleration is also at grid points so the update a n = a n, n = a n+ = n + n+/ n, n/ + n+/ n+/ = n/ + a n is also second order accurate We also hae the second order accurate central difference formula n+ n + n = a n We can combine these formulas to obtain the eplicit formula for a half time step n+/ n = n/ + n+/ n = n+/ a n + n+/ = n + an which is then used to update positions Note that these two steps are equialent to n+ = n + n + an Finally, elocities are updated another half step using the implicit update The oerall elocity update is n+ = n+/ + n+3/ n+ = n+/ + n+/ + a n+ n+ = n+/ + an+ = n+/ + an+, n+ n+ = n+/ + an+, n+ = n + an, n + an+, n+ which is just the trapezoid rule The dependence of acceleration on elocity is often symmetric eg, damping forces, so that a fast soler like preconditioned conjugate gradient can be employed 6
Lecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationDynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2)
Dynamics ( 동역학 ) Ch. Motion of Translating Bodies (. &.) Motion of Translating Bodies This chapter is usually referred to as Kinematics of Particles. Particles: In dynamics, a particle is a body without
More informationDrive train. Steering System. Figure 1 Vehicle modeled by subsystems
Proceedings of the XII International Symposium on Dynamic Problems of Mechanics (DINAME 7), P.S.Varoto and M.A.Trindade (editors), ABCM, Ilhabela, SP, Brazil, February 6 - March, 7 Wheel Dynamics Georg
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationСollisionless damping of electron waves in non-maxwellian plasma 1
http:/arxi.org/physics/78.748 Сollisionless damping of electron waes in non-mawellian plasma V. N. Soshnio Plasma Physics Dept., All-Russian Institute of Scientific and Technical Information of the Russian
More informationChapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs
Chapter 9 Implicit Methods for Linear and Nonlinear Systems of ODEs In the previous chapter, we investigated stiffness in ODEs. Recall that an ODE is stiff if it exhibits behavior on widelyvarying timescales.
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationA. Idesman. Keywords: time integration, spurious oscillations, numerical dispersion
COMPDYN 0 rd ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Pleris (eds.) Corfu, Greece, -8 May 0 ACCURATE NUMERICAL
More informationStability of Mass-Point Systems
Simulation in Computer Graphics Stability of Mass-Point Systems Matthias Teschner Computer Science Department University of Freiburg Demos surface tension vs. volume preservation distance preservation
More informationQuadratic Equation. ax bx c =. = + + =. Example 2. = + = + = 3 or. The solutions are -7/3 and 1.
Quadratic Equation A quadratic equation is any equation that is equialent to the equation in fmat a c + + = 0 (1.1) where a,, and c are coefficients and a 0. The ariale name is ut the same fmat applies
More informationSemi-implicit Treatment of the Hall Effect in NIMROD Simulations
Semi-implicit Treatment of the Hall Effect in NIMROD Simulations H. Tian and C. R. Soinec Department of Engineering Physics, Uniersity of Wisconsin-Madison Madison, WI 5376 Presented at the 45th Annual
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationDepartment of Physics PHY 111 GENERAL PHYSICS I
EDO UNIVERSITY IYAMHO Department o Physics PHY 111 GENERAL PHYSICS I Instructors: 1. Olayinka, S. A. (Ph.D.) Email: akinola.olayinka@edouniersity.edu.ng Phone: (+234) 8062447411 2. Adekoya, M. A Email:
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationLecture 21: Physical Brownian Motion II
Lecture 21: Physical Brownian Motion II Scribe: Ken Kamrin Department of Mathematics, MIT May 3, 25 Resources An instructie applet illustrating physical Brownian motion can be found at: http://www.phy.ntnu.edu.tw/jaa/gas2d/gas2d.html
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationMath 660 Lecture 4: FDM for evolutionary equations: ODE solvers
Math 660 Lecture 4: FDM for evolutionary equations: ODE solvers Consider the ODE u (t) = f(t, u(t)), u(0) = u 0, where u could be a vector valued function. Any ODE can be reduced to a first order system,
More informationVelocity, Acceleration and Equations of Motion in the Elliptical Coordinate System
Aailable online at www.scholarsresearchlibrary.com Archies of Physics Research, 018, 9 (): 10-16 (http://scholarsresearchlibrary.com/archie.html) ISSN 0976-0970 CODEN (USA): APRRC7 Velocity, Acceleration
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationDYNAMICS. Kinematics of Particles Engineering Dynamics Lecture Note VECTOR MECHANICS FOR ENGINEERS: Eighth Edition CHAPTER
27 The McGraw-Hill Companies, Inc. All rights resered. Eighth E CHAPTER 11 VECTOR MECHANICS FOR ENGINEERS: DYNAMICS Ferdinand P. Beer E. Russell Johnston, Jr. Kinematics of Particles Lecture Notes: J.
More informationChapter 2: 1D Kinematics Tuesday January 13th
Chapter : D Kinematics Tuesday January 3th Motion in a straight line (D Kinematics) Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Short summary Constant acceleration a special
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationThe Full-rank Linear Least Squares Problem
Jim Lambers COS 7 Spring Semeseter 1-11 Lecture 3 Notes The Full-rank Linear Least Squares Problem Gien an m n matrix A, with m n, and an m-ector b, we consider the oerdetermined system of equations Ax
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationIntroduction to Differential Equations
Math0 Lecture # Introduction to Differential Equations Basic definitions Definition : (What is a DE?) A differential equation (DE) is an equation that involves some of the derivatives (or differentials)
More informationTime scales. Notes. Mixed Implicit/Explicit. Newmark Methods [ ] [ ] "!t = O 1 %
Notes Time scales! For using Pixie (the renderer) make sure you type use pixie first! Assignment 1 questions?! [work out]! For position dependence, characteristic time interval is "!t = O 1 % $ ' # K &!
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationSplit explicit methods
Split explicit methods Almut Gassmann Meteorological Institute of the University of Bonn Germany St.Petersburg Summer School 2006 on nonhydrostatic dynamics and fine scale data assimilation Two common
More informationFeb 6, 2013 PHYSICS I Lecture 5
95.141 Feb 6, 213 PHYSICS I Lecture 5 Course website: faculty.uml.edu/pchowdhury/95.141/ www.masteringphysics.com Course: UML95141SPRING213 Lecture Capture h"p://echo36.uml.edu/chowdhury213/physics1spring.html
More informationComputational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations
Computational Methods CMSC/AMSC/MAPL 460 Ordinar differential equations Ramani Duraiswami, Dept. of Computer Science Several slides adapted from Prof. ERIC SANDT, TAMU ODE: Previous class Standard form
More informationNote: the net distance along the path is a scalar quantity its direction is not important so the average speed is also a scalar.
PHY 309 K. Solutions for the first mid-term test /13/014). Problem #1: By definition, aerage speed net distance along the path of motion time. 1) ote: the net distance along the path is a scalar quantity
More informationDO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.
DO PHYSICS ONLINE DISPLACEMENT VELOCITY ACCELERATION The objects that make up space are in motion, we moe, soccer balls moe, the Earth moes, electrons moe, - - -. Motion implies change. The study of the
More informationSimulations of bulk phases. Periodic boundaries. Cubic boxes
Simulations of bulk phases ChE210D Today's lecture: considerations for setting up and running simulations of bulk, isotropic phases (e.g., liquids and gases) Periodic boundaries Cubic boxes In simulations
More informationStability regions of Runge-Kutta methods. Stephan Houben Eindhoven University of Technology
Stability regions of Runge-Kutta methods Stephan Houben Eindhoven University of Technology February 19, 2002 1 Overview of the talk 1. Quick review of some concepts 2. Stability regions 3. L-stability
More informationLecture 17: Ordinary Differential Equation II. First Order (continued)
Lecture 17: Ordinary Differential Equation II. First Order (continued) 1. Key points Maple commands dsolve dsolve[interactive] dsolve(numeric) 2. Linear first order ODE: y' = q(x) - p(x) y In general,
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationDynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems
. Ordinary Differential Equations. An ordinary differential equation (ODE, for short) is an equation involves a dependent variable and its derivatives with respect to an independent variable. When we speak
More informationLecture 4.5 Schemes for Parabolic Type Equations
Lecture 4.5 Schemes for Parabolic Type Equations 1 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationSection 4.3: Quadratic Formula
Objective: Solve quadratic equations using the quadratic formula. In this section we will develop a formula to solve any quadratic equation ab c 0 where a b and c are real numbers and a 0. Solve for this
More informationEuler s Method, cont d
Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example
More informationComputational Methods CMSC/AMSC/MAPL 460. Ordinary differential equations
Computational Methods CMSC/AMSC/MAPL 460 Ordinar differential equations Ramani Duraiswami, Dept. of Computer Science Several slides adapted from Prof. ERIC SANDT, TAMU ODE: Previous class Applications
More informationTHE FIFTH DIMENSION EQUATIONS
JP Journal of Mathematical Sciences Volume 7 Issues 1 & 013 Pages 41-46 013 Ishaan Publishing House This paper is aailable online at http://www.iphsci.com THE FIFTH DIMENSION EQUATIONS Niittytie 1B16 03100
More informationEuler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability.
Euler s Method, Taylor Series Method, Runge Kutta Methods, Multi-Step Methods and Stability. REVIEW: We start with the differential equation dy(t) dt = f (t, y(t)) (1.1) y(0) = y 0 This equation can be
More informationLecture 5: Single Step Methods
Lecture 5: Single Step Methods J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 1 October 2012 () Single Step Methods 1 October 2012 1 / 44 Outline 1
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationSection 6: PRISMATIC BEAMS. Beam Theory
Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam
More informationCS 205b / CME 306. Application Track. Homework 2
CS 205b / CME 306 Application Track Homework 2 1. ALE An Eulerian formulation of conservation of mass uses control volumes that are fixed in space as material flows freely through the control volumes.
More informationIII. Relative Velocity
Adanced Kinematics I. Vector addition/subtraction II. Components III. Relatie Velocity IV. Projectile Motion V. Use of Calculus (nonuniform acceleration) VI. Parametric Equations The student will be able
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 27 Theme of Last Three
More informationGeneral Lorentz Boost Transformations, Acting on Some Important Physical Quantities
General Lorentz Boost Transformations, Acting on Some Important Physical Quantities We are interested in transforming measurements made in a reference frame O into measurements of the same quantities as
More informationBoundary Layer (Reorganization of the Lecture Notes from Professor Anthony Jacobi and Professor Nenad Miljkoic) Consider a steady flow of a Newtonian, Fourier-Biot fluid oer a flat surface with constant
More informationMATHEMATICAL MODELLING AND IDENTIFICATION OF THE FLOW DYNAMICS IN
MATHEMATICAL MOELLING AN IENTIFICATION OF THE FLOW YNAMICS IN MOLTEN GLASS FURNACES Jan Studzinski Systems Research Institute of Polish Academy of Sciences Newelska 6-447 Warsaw, Poland E-mail: studzins@ibspan.waw.pl
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationLesson 2: Kinematics (Sections ) Chapter 2 Motion Along a Line
Lesson : Kinematics (Sections.-.5) Chapter Motion Along a Line In order to specify a position, it is necessary to choose an origin. We talk about the football field is 00 yards from goal line to goal line,
More informationWhat we ll do: Lecture 21. Ordinary Differential Equations (ODEs) Differential Equations. Ordinary Differential Equations
What we ll do: Lecture Ordinary Differential Equations J. Chaudhry Department of Mathematics and Statistics University of New Mexico Review ODEs Single Step Methods Euler s method (st order accurate) Runge-Kutta
More informationLecture #9: Harmonic Oscillator: Creation and Annihilation Operators
5.61 Fall, 013 Lecture #9 Page 1 Last time Lecture #9: Harmonic Oscillator: Creation and Annihilation Operators 1/ Simplified Schrödinger equation: ξ=α x, α =(kμ) 1/ n E +ξ ψ=0 (dimensionless) ξ nω reduced
More informationGeostrophy & Thermal wind
Lecture 10 Geostrophy & Thermal wind 10.1 f and β planes These are planes that are tangent to the earth (taken to be spherical) at a point of interest. The z ais is perpendicular to the plane (anti-parallel
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationA possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k
A possible mechanism to explain wae-particle duality L D HOWE No current affiliation PACS Numbers: 0.50.-r, 03.65.-w, 05.60.-k Abstract The relationship between light speed energy and the kinetic energy
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
More informationMotion in Two and Three Dimensions
PH 1-1D Spring 013 Motion in Two and Three Dimensions Lectures 5,6,7 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter
More informationReading. Particle Systems. What are particle systems? Overview
Readig Particle Systems Required: Witki, Particle System Dyamics, SIGGRAPH 97 course otes o Physically Based Modelig. Optioal Witki ad Baraff, Differetial Equatio Basics, SIGGRAPH 97 course otes o Physically
More informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 18. HW 7 is posted, and will be due in class on 4/25.
Logistics Week 12: Monday, Apr 18 HW 6 is due at 11:59 tonight. HW 7 is posted, and will be due in class on 4/25. The prelim is graded. An analysis and rubric are on CMS. Problem du jour For implicit methods
More informationInvestigation on Ring Valve Motion and Impact Stress in Reciprocating Compressors
Purdue Uniersity Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2010 Inestigation on Ring Vale Motion and Impact Stress in Reciprocating Compressors Yu Wang
More informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
More informationSecond Order Transfer Function Discrete Equations
Second Order Transfer Function Discrete Equations J. Riggs 23 Aug 2017 Transfer Function Equations pg 1 1 Introduction The objective of this paper is to develop the equations for a discrete implementation
More information4-1 Graphing Quadratic Functions
4-1 Graphing Quadratic Functions Quadratic Function in standard form: f() a b c The graph of a quadratic function is a. y intercept Ais of symmetry -coordinate of verte coordinate of verte 1) f ( ) 4 a=
More informationAn example of Lagrangian for a non-holonomic system
Uniersit of North Georgia Nighthaks Open Institutional Repositor Facult Publications Department of Mathematics 9-9-05 An eample of Lagrangian for a non-holonomic sstem Piotr W. Hebda Uniersit of North
More informationMA 266 Review Topics - Exam # 2 (updated)
MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,
More informationHigh Order Well-Balanced Schemes and Applications to Non-Equilibrium Flow with Stiff Source Terms
High Order Well-Balanced Schemes and Applications to Non-Equilibrium Flow with Stiff Source Terms Wei Wang, Chi-Wang Shu, H. C. Yee and Björn Sjögreen February 12, 29 Abstract The stiffness of the source
More informationSupplementary Information Microfluidic quadrupole and floating concentration gradient Mohammad A. Qasaimeh, Thomas Gervais, and David Juncker
Mohammad A. Qasaimeh, Thomas Gerais, and Daid Juncker Supplementary Figure S1 The microfluidic quadrupole (MQ is modeled as two source (Q inj and two drain (Q asp points arranged in the classical quardupolar
More informationMecE 390 Final examination, Winter 2014
MecE 390 Final examination, Winter 2014 Directions: (i) a double-sided 8.5 11 formula sheet is permitted, (ii) no calculators are permitted, (iii) the exam is 80 minutes in duration; please turn your paper
More informationSUPPLEMENTARY MATERIAL. Authors: Alan A. Stocker (1) and Eero P. Simoncelli (2)
SUPPLEMENTARY MATERIAL Authors: Alan A. Stocker () and Eero P. Simoncelli () Affiliations: () Dept. of Psychology, Uniersity of Pennsylania 34 Walnut Street 33C Philadelphia, PA 94-68 U.S.A. () Howard
More information[ ] is a vector of size p.
Lecture 11 Copyright by Hongyun Wang, UCSC Recap: General form of explicit Runger-Kutta methods for solving y = F( y, t) k i = hfy n + i1 j =1 c ij k j, t n + d i h, i = 1,, p + p j =1 b j k j A Runge-Kutta
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) NRiC Chapter 16. ODEs involve derivatives wrt one independent variable, e.g. time t. ODEs can always be reduced to a set of firstorder equations (involving only first
More informationMain points of today s lecture: Example: addition of velocities Trajectories of objects in 2 dimensions: Physic 231 Lecture 5 ( )
Main points of toda s lecture: Eample: addition of elocities Trajectories of objects in dimensions: Phsic 31 Lecture 5 ( ) t g gt t t gt o 1 1 downwards 9.8 m/s g Δ Δ Δ + Δ Motion under Earth s graitational
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 5 Chapter 21 Numerical Differentiation PowerPoints organized by Dr. Michael R. Gustafson II, Duke University 1 All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
More informationChapter 1: Kinematics of Particles
Chapter 1: Kinematics of Particles 1.1 INTRODUCTION Mechanics the state of rest of motion of bodies subjected to the action of forces Static equilibrium of a body that is either at rest or moes with constant
More informationMultistep Methods for IVPs. t 0 < t < T
Multistep Methods for IVPs We are still considering the IVP dy dt = f(t,y) t 0 < t < T y(t 0 ) = y 0 So far we have looked at Euler s method, which was a first order method and Runge Kutta (RK) methods
More informationSome notes about PDEs. -Bill Green Nov. 2015
Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can
More informationMath Numerical Analysis Homework #4 Due End of term. y = 2y t 3y2 t 3, 1 t 2, y(1) = 1. n=(b-a)/h+1; % the number of steps we need to take
Math 32 - Numerical Analysis Homework #4 Due End of term Note: In the following y i is approximation of y(t i ) and f i is f(t i,y i ).. Consider the initial value problem, y = 2y t 3y2 t 3, t 2, y() =.
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 32 Theme of Last Few
More informationSynchronization procedures and light velocity
Synchronization procedures and light elocity Joseph Léy Mail address : 4, square Anatole France 95 St-Germain-lès-orbeil France E-mail : josephley@compusere.com Abstract Although different arguments speak
More information