ODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits
|
|
- Kristin Small
- 5 years ago
- Views:
Transcription
1 Solving ODEs
2 General Stuff ODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits need workhorse solvers to deal with them PDEs are more common
3 General Stuff these methods apply to all ODEs higher order ODEs are transformed to systems of first order ODEs, e.g., transforms to dy 1 dx dy 2 dx d 2 y dx 2 + q(x)dy dx = r(x) with y(x) = y 1 (x) and dy/dx = y 2 (x) = y 2 (x) (1) = r(x) q(x)y 2 (x) (2)
4 General Stuff need methods to solve a set of ODEs dy i dx = f i(x, y (1...i) ) equations of more than one independent variable PDEs require far more complex methods next chapter
5 Boundary Conditions need BCs to solve a ODE problem numerically analytical solution integration constants type of BCs determine numerical approach, e.g.: 1. initial value problem: all y i are given at one starting point x s 2. two-point boundary value problem: some y i are given at x s, the others at some point x f first: initial value problems
6 Euler method consider only ODE with i = 1 for simplicity! index y n value of y(x) at x = x n, i.e., the nth point going from x s to x f simplest approach: Euler chose stepsize h for x then write y n+1 = y n + hf (x n, y n ) uses information only at the beginning of an interval expansion in power series first order accuracy: error is O(h 2 )
7 Euler method
8 Euler method not recommended for practical use: Euler Predictor Corrector is better other methods have better accuracy for same step size instability problems (see below)
9 Euler Predictor-Corrector then slope of chord is y 1 = y 0 + d f d x x 0 h (3) y 0 = y 1 d f d x x 1 h (4) y 1 y 0 h = 1 2 {f (x 0, y 0 ) + f (x 1, y 1 )} but don t know y 1 use Euler step to estimate (predictor step) Y 1 = y 0 + f (x 0, y 0 ) h then use cord to advance the equation y 1 = y h{f (x 0, y 0 ) + f (x 1, Y 1 )}
10 Euler PC Method
11 Runge-Kutta method use Euler to take trial step to x n + h/2 use Euler estimate to approximate x and y at midpoint use these to compute full step midpoint method second order Runge-Kutta method
12 Midpoint Method
13 Runge-Kutta method formally: k 1 = hf (x n, y n ) (5) k 2 = hf (x n + h/2, y n + k 1 /2) (6) y n+1 = y n + k 2 + O(h 3 ) (7) second order accurate! needs 2 evals of f () for stepsize h can be extended by using several estimates at start/mid/endpoints accuracy can be increased
14 4th order RK method
15 Runge-Kutta method formally: k 1 = hf (x n, y n ) (8) k 2 = hf (x n + h/2, y n + k 1 /2) (9) k 3 = hf (x n + h/2, y n + k 2 /2) (10) k 4 = hf (x n + h, y n + k 3 ) (11) y n+1 = y n + k 1 /6 + k 3 /3 + k4/6 + O(h 5 ) (12) fourth order accurate better than 2nd order if h is at least 2 times larger for same quality
16 step size control very often f () changes significantly over integration interval to ensure homogeneous quality we need to adapt h needs to be really small in domains where f () has nasty features h can be larger for smooth, slowly varying functions
17 step size control general ODE solver needs to track accuracy to obtain best performance for given quality step-doubling: take each step twice as h and two h/2 steps total of 3 RK steps with 4 function evals each starting point shared 11 evals overhead is 11/ (we reach accuracy for step size h/2!)
18 step size control indicate truncation error difference between the two estimates use this information to adapt h automatically method is 4th order scales as h 5 ( ) h 0 = h 1 1 given a target truncation error 0, use this to estimate h 0 from current 1 and h 1 if 1 > 0 decrease h otherwise increase it
19 Fehlberg method RK of M > 4th order up to M + 2 evals 4th order most efficient? Fehlberg 5th order method with 6 evals these 6 points can also be used to build a 4th order formula combine them to compute truncation error with only 6 evals! embedded RK-Fehlberg method
20 Embedded RKF method
21 Embedded RKF method
22 Modified Midpoint Method go through a (large) step H with h/n sub-steps use an algorithm of the form z 0 = y(x) (13) z 1 = z 0 + hf (x, z 0 ) (14) z m+1 = z m 1 + 2hf (x + mh, z m ) (15) y(z + H) = 1 2 [z n + z n 1 + hf (x + H, z n )] (16) same as midpoint method but first and last points modified midpoint method
23 Modified Midpoint Method truncation error y n y(x + H) = i α i h 2i only even powers in h combining steps will deliver two orders better accuracy per doubling usually adaptive RKF s are better but modified midpoint is very useful for...
24 Richardson extrapolation
25 Bulirsch-Stoer method original Richardson power series limited convergence radius! Bulirsch-Stoer: use modified midpoint method use rational functions in h 2 can take rather large steps H
26 Bulirsch-Stoer method how to subdivide steps best? H/n j Bulirsch-Stoer: n j = 2n j 2 more efficient: n j = 2j (Deuflhard)
27 Bulirsch-Stoer method j max is not determined a priori extrapolation delivers error limits use them to limit steps variation: polynomial extrapolation
28 Bulirsch-Stoer method method works extremely well for well behaved ODEs smooth functions it does not work well for function evals through table lookup or interpolation discontinuous functions internal singular points steep changes of f through an interval
29 Predictor-Corrector Methods write the solution of the ODE as y n+1 = y n + xn+1 x n f (t, y) dt approximate f (t, y) by a polynomial through a number of previous points x n (?), x n 1, x n 2 multistep method result of has the form y n+1 = y n +h(β 0 f (x n+1, y n+1 )+β 1 f (x n, y n )+β 2 f (x n 1, y n 1 )+ )
30 Predictor-Corrector Methods β 0 = 0 explicit, otherwise implicit how to solve? Newton s method multistep solver for stiff problems
31 Predictor-Corrector Methods functional iteration predictor-corrector methods idea: use explicit method to get estimate y n+1 predictor step next: use estimate to correct y n+1 corrector step can be iterated if number of iterations is fixed at the beginning explicit method (why?)
32 Adams-Bashforth-Moulton method overall good stability (explicit!) 3rd order version: predictor step: Adams-Bashforth y n+1 = y n + h 12 corrector step: Adams-Moulton ( 23y n 16y n 1 + 5y n 2 ) + O(h 4 ) y n+1 = y n + h 12 ( 5y n+1 + 8y n y n 1 ) + O(h 4 )
33 predictor-corrector methods have starting/stopping problems (use RK) stepsize control hard use canned routines whenever possible historically very frequently used very good for very smooth functions that are expensive to evaluate
34 ODEINT Useful Workhorse from Num Rec Can be called with either RK4 (Runge-Kutta) OR with BSSTEP (Burlisch-Stoer)
35 My favorite method Hammings Modified Predictor Corrector Method Adaptive Stepsize control Written at IBM in the 1960s 4th order method Uses RK4 to start Will go exactly to b if initial stepsize is (b a)/2 n
Math 225 Scientific Computing II Outline of Lectures
Math 225 Scientific Computing II Outline of Lectures Spring Semester 2003 I. Interpolating polynomials Lagrange formulation of interpolating polynomial Uniqueness of interpolating polynomial of degree
More informationOrdinary Differential Equations
Next: Partial Differential Equations Up: Numerical Analysis for Chemical Previous: Numerical Differentiation and Integration Subsections Runge-Kutta Methods Euler's Method Improvement of Euler's Method
More informationover The idea is to construct an algorithm to solve the IVP ODE (9.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (9.1) over To obtain the solution point we can use the fundamental theorem
More informationx n x n stepnumber k order r error constant C r+1 1/2 5/12 3/8 251/720 abs. stab. interval (α,0) /11-3/10
MATH 573 LECTURE NOTES 77 13.8. Predictor-corrector methods. We consider the Adams methods, obtained from the formula xn+1 xn+1 y(x n+1 y(x n ) = y (x)dx = f(x,y(x))dx x n x n by replacing f by an interpolating
More informationSimulation in Computer Graphics. Particles. Matthias Teschner. Computer Science Department University of Freiburg
Simulation in Computer Graphics Particles Matthias Teschner Computer Science Department University of Freiburg Outline introduction particle motion finite differences system of first order ODEs second
More informationEuler s Methods (a family of Runge- Ku9a methods)
Euler s Methods (a family of Runge- Ku9a methods) ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationover The idea is to construct an algorithm to solve the IVP ODE (8.1)
Runge- Ku(a Methods Review of Heun s Method (Deriva:on from Integra:on) The idea is to construct an algorithm to solve the IVP ODE (8.1) over To obtain the solution point we can use the fundamental theorem
More information16.1 Runge-Kutta Method
70 Chapter 6. Integration of Ordinary Differential Equations CITED REFERENCES AND FURTHER READING: Gear, C.W. 97, Numerical Initial Value Problems in Ordinary Differential Equations (Englewood Cliffs,
More informationMATH2071: LAB 2: Explicit ODE methods
MATH2071: LAB 2: Explicit ODE methods 1 Introduction Introduction Exercise 1 Euler s method review Exercise 2 The Euler Halfstep (RK2) Method Exercise 3 Runge-Kutta Methods Exercise 4 The Midpoint Method
More informationAn Introduction to Flow Visualization (1) Christoph Garth
An Introduction to Flow Visualization (1) Christoph Garth cgarth@ucdavis.edu Motivation What will I be talking about? Classical: Physical experiments to understand flow. 2 Motivation What will I be talking
More informationODE IVP. An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable:
Euler s Methods ODE IVP An Ordinary Differential Equation (ODE) is an equation that contains a function having one independent variable: The equation is coupled with an initial value/condition (i.e., value
More informationMathematics for chemical engineers
Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester 2013-2014 Numerical solution of ordinary differential equations Initial value problem Outline 1 Introduction
More informationHomework Set #2-3, Math 475B
Homework Set #2-3, Math 475B Part I: Matlab In the last semester you learned a number of essential features of MATLAB. 1. In this instance, you will learn to make 3D plots and contour plots of z = f(x,
More informationThere are also formulas of higher order for this situation, but we will refrain from giving them.
136 Chapter 4. Integration of Functions N = 1 2 3 4 (total after N = 4) Figure 4.2.1. Sequential calls to the routine trapzd incorporate the information from previous calls and evaluate the integrand only
More informationChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab
ChE 400: Applied Chemical Engineering Calculations Tutorial 6: Numerical Solution of ODE Using Excel and Matlab Tutorial 6: Numerical Solution of ODE Gerardine G. Botte This handout contains information
More informationOrdinary differential equations solving methods
Radim Hošek, A07237 radhost@students.zcu.cz Ordinary differential equations solving methods Problem: y = y2 (1) y = x y (2) y = sin ( + y 2 ) (3) Where it is possible we try to solve the equations analytically,
More informationNumerical Integration
1 Numerical Integration Jim Van Verth Insomniac Games jim@essentialmath.com Essential mathematics for games and interactive applications Talk Summary Going to talk about: Euler s method subject to errors
More informationModule1: Numerical Solution of Ordinary Differential Equations. Lecture 6. Higher order Runge Kutta Methods
Module1: Numerical Solution of Ordinary Differential Equations Lecture 6 Higher order Runge Kutta Methods Keywords: higher order methods, functional evaluations, accuracy Higher order Runge Kutta Methods
More informationFinite difference methods
Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 / Introduction Finite difference methods
More informationNumerical Resolution of optimal control problems by a Piecewise Linear continuation method
Numerical Resolution of optimal control problems by a Piecewise Linear continuation method Pierre Martinon november 25 Shooting method and homotopy Introduction We would like to solve numerically the
More informationPart I: Theoretical Background and Integration-Based Methods
Large Vector Field Visualization: Theory and Practice Part I: Theoretical Background and Integration-Based Methods Christoph Garth Overview Foundations Time-Varying Vector Fields Numerical Integration
More informationEuler s Method for Approximating Solution Curves
Euler s Method for Approximating Solution Curves As you may have begun to suspect at this point, time constraints will allow us to learn only a few of the many known methods for solving differential equations.
More informationAccuracy Analysis of Charged Particle Trajectory CAE Software
www.integratedsoft.com Accuracy Analysis of Charged Particle Trajectory CAE Software Content Executive Summary... 3 Overview of Charged Particle Beam Analysis... 3 Types of Field Distribution... 4 Simulating
More information(Part - 1) P. Sam Johnson. April 14, Numerical Solution of. Ordinary Differential Equations. (Part - 1) Sam Johnson. NIT Karnataka.
P. April 14, 2015 1/51 Overview We discuss the following important methods of solving ordinary differential equations of first / second order. Picard s method of successive approximations (Method of successive
More informationReview Initial Value Problems Euler s Method Summary
THE EULER METHOD P.V. Johnson School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 INITIAL VALUE PROBLEMS The Problem Posing a Problem 3 EULER S METHOD Method Errors 4 SUMMARY OUTLINE 1 REVIEW 2 INITIAL
More informationLecture 22: Rules for Solving IVP
cs412: introduction to numerical analysis 12/7/10 Lecture 22: Rules for Solving IVP Instructor: Professor Amos Ron Scribes: Yunpeng Li, Pratap Ramamurthy, Mark Cowlishaw, Nathanael Fillmore 1 Review Recall
More information4 Visualization and. Approximation
4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential
More informationComputer Simulations
Computer Simulations A practical approach to simulation Semra Gündüç gunduc@ankara.edu.tr Ankara University Faculty of Engineering, Department of Computer Engineering 2014-2015 Spring Term Ankara University
More informationPrecision Spin Tracking for Electric Dipole Moment Searches
Precision Spin Tracking for Electric Dipole Moment Searches Outline: Spin Evolution in Storage Rings Technique for EDM Measurement and Systematic Errors Integration Methods, Current Status of Our Program
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 16 Lecture 16 May 3, 2018 Numerical solution of systems
More informationDiscontinuous Galerkin Sparse Grid method for Maxwell s equations
Discontinuous Galerkin Sparse Grid method for Maxwell s equations Student: Tianyang Wang Mentor: Dr. Lin Mu, Dr. David L.Green, Dr. Ed D Azevedo, Dr. Kwai Wong Motivation u We consider the Maxwell s equations,
More informationMass-Spring Systems. Last Time?
Mass-Spring Systems Last Time? Implicit Surfaces & Marching Cubes/Tetras Collision Detection & Conservative Bounding Regions Spatial Acceleration Data Structures Octree, k-d tree, BSF tree 1 Today Particle
More informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
More informationThe listing says y(1) = How did the computer know this?
18.03 Class 2, Feb 3, 2010 Numerical Methods [1] How do you know the value of e? [2] Euler's method [3] Sources of error [4] Higher order methods [1] The study of differential equations has three parts:.
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationAsymptotic Error Analysis
Asymptotic Error Analysis Brian Wetton Mathematics Department, UBC www.math.ubc.ca/ wetton PIMS YRC, June 3, 2014 Outline Overview Some History Romberg Integration Cubic Splines - Periodic Case More History:
More informationGPU Implementation of Implicit Runge-Kutta Methods
GPU Implementation of Implicit Runge-Kutta Methods Navchetan Awasthi, Abhijith J Supercomputer Education and Research Centre Indian Institute of Science, Bangalore, India navchetanawasthi@gmail.com, abhijith31792@gmail.com
More informationNumerical Simulation of Dynamic Systems XXIV
Numerical Simulation of Dynamic Systems XXIV Prof. Dr. François E. Cellier Department of Computer Science ETH Zurich May 14, 2013 Introduction Introduction A number of important simulation applications
More informationApplication 2.4 Implementing Euler's Method
Application 2.4 Implementing Euler's Method One's understanding of a numerical algorithm is sharpened by considering its implementation in the form of a calculator or computer program. Figure 2.4.13 in
More informationIntroduction to Programming with Matlab. Jim Dove
Introduction to Programming with Matlab Jim Dove Spring Semester, 2005 Chapter 1 INTRODUCTION 1.1 The Need for Computers Statistical physics or N-body simulations not possible Monte-Carlo Simulations Many
More informationC++ Playground for Numerical Integration Method Developers
C++ Playground for Numerical Integration Method Developers Stepan Orlov Computer Technologies in Engineering dept. Peter the Great St. Petersburg Polytechnic University St. Petersburg, Russian Federation
More information16.2 Adaptive Stepsize Control for Runge-Kutta
708 Chapter 6. Integration of Ordinary Differential Equations call rk4(v,dv,nvar,x,h,v,derivs) if(x+h.eq.x)pause stepsize not significant in rkdumb x=x+h xx(k+)=x Store intermediate steps. do 2 i=,nvar
More informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationData-Driven Modeling. Scientific Computation J. NATHAN KUTZ OXPORD. Methods for Complex Systems & Big Data
Data-Driven Modeling & Scientific Computation Methods for Complex Systems & Big Data J. NATHAN KUTZ Department ofapplied Mathematics University of Washington OXPORD UNIVERSITY PRESS Contents Prolegomenon
More informationStepsize this trip. General step.
724 Chapter 16. Integration of Ordinary Differential Equations ym=vector(1,nvar); yn=vector(1,nvar); h=htot/nstep; for (i=1;i
More informationA Study on Numerical Exact Solution of Euler, Improved Euler and Runge - Kutta Method
A Study on Numerical Exact Solution of Euler, Improved Euler and Runge - Kutta Method Mrs. P. Sampoornam P.K.R Arts College for Women, Gobichettipalayam Abstract: In this paper, I will discuss the Runge
More informationLecture 8: Euler s Methods
Lecture 8: Euler s Methods Forward Euler Method The formula for the forward Euler method is given by equation (8.2) in the lecture note for week 8, as y i+1 = y i + f(x i, y i )h. (1) where f(x i, y i
More informationCommunications in Applied Mathematics and Computational Science
Communications in Applied Mathematics and Computational Science A HYBRID PARAREAL SPECTRAL DEFERRED CORRECTIONS METHOD MICHAEL L. MINION vol. 5 no. 2 2010 mathematical sciences publishers COMM. APP. MATH.
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationNumerical Methods for Hyperbolic and Kinetic Equations
Numerical Methods for Hyperbolic and Kinetic Equations Organizer: G. Puppo Phenomena characterized by conservation (or balance laws) of physical quantities are modelled by hyperbolic and kinetic equations.
More informationMulti-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech
Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated
More informationSuppose we want to find approximate values for the solution of the differential equation
A course in differential equations Chapter 2. Numerical methods of solution We shall discuss the simple numerical method suggested in the last chapter, and explain several more sophisticated ones. We shall
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationThe Jello Cube Assignment 1, CSCI 520. Jernej Barbic, USC
The Jello Cube Assignment 1, CSCI 520 Jernej Barbic, USC 1 The jello cube Undeformed cube Deformed cube The jello cube is elastic, Can be bent, stretched, squeezed,, Without external forces, it eventually
More informationFirst of all, we need to know what it means for a parameterize curve to be differentiable. FACT:
CALCULUS WITH PARAMETERIZED CURVES In calculus I we learned how to differentiate and integrate functions. In the chapter covering the applications of the integral, we learned how to find the length of
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
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 informationSolving systems of nonlinear equations with arbitrary precision
Solving systems of nonlinear equations with arbitrary precision CES Seminar January 2017 Presented by Supervisor Gregor Heiming Dr. Ralf Hannemann-Tamás Systemverfahrenstechnik AVT, RWTH Aachen University
More informationExercises for a Numerical Methods Course
Exercises for a Numerical Methods Course Brian Heinold Department of Mathematics and Computer Science Mount St. Mary s University November 18, 2017 1 / 73 About the class Mix of Math and CS students, mostly
More informationParallel implicit ordinary differential equation solver for cuda. Tomasz M. Kardaś
Parallel implicit ordinary differential equation solver for cuda Tomasz M. Kardaś August 11, 2014 Chapter 1 Parallel Implicit Ordinary Differential Equations Solver A simplest definition of stiffness,
More informationProf. B.S. Thandaveswara. The computation of a flood wave resulting from a dam break basically involves two
41.4 Routing The computation of a flood wave resulting from a dam break basically involves two problems, which may be considered jointly or seperately: 1. Determination of the outflow hydrograph from the
More informationSPH: Towards the simulation of wave-body interactions in extreme seas
SPH: Towards the simulation of wave-body interactions in extreme seas Guillaume Oger, Mathieu Doring, Bertrand Alessandrini, and Pierre Ferrant Fluid Mechanics Laboratory (CNRS UMR6598) Ecole Centrale
More informationModule 2: Single Step Methods Lecture 4: The Euler Method. The Lecture Contains: The Euler Method. Euler's Method (Analytical Interpretations)
The Lecture Contains: The Euler Method Euler's Method (Analytical Interpretations) An Analytical Example file:///g /Numerical_solutions/lecture4/4_1.htm[8/26/2011 11:14:40 AM] We shall now describe methods
More informationAMSC/CMSC 460 Final Exam, Fall 2007
AMSC/CMSC 460 Final Exam, Fall 2007 Show all work. You may leave arithmetic expressions in any form that a calculator could evaluate. By putting your name on this paper, you agree to abide by the university
More informationAPPM/MATH Problem Set 4 Solutions
APPM/MATH 465 Problem Set 4 Solutions This assignment is due by 4pm on Wednesday, October 16th. You may either turn it in to me in class on Monday or in the box outside my office door (ECOT 35). Minimal
More informationMath 133A, September 10: Numerical Methods (Euler & Runge-Kutta) Section 1: Numerical Methods
Math 133A, September 10: Numerical Methods (Euler & Runge-Kutta) Section 1: Numerical Methods We have seen examples of first-order differential equations which can, by one method of another, be solved
More informationDirection Fields; Euler s Method
Direction Fields; Euler s Method It frequently happens that we cannot solve first order systems dy (, ) dx = f xy or corresponding initial value problems in terms of formulas. Remarkably, however, this
More informationBuilding More Efficient Time Parallel Solvers with Applications to Multiscale Modeling and Gyrokinetics
Building More Efficient Time Parallel Solvers with Applications to Multiscale Modeling and Gyrokinetics January 2015 Carl D Lederman AFTC/PA clearance No. 15008, 16 January 2014 1 The Differential Equation
More informationf( x ), or a solution to the equation f( x) 0. You are already familiar with ways of solving
The Bisection Method and Newton s Method. If f( x ) a function, then a number r for which f( r) 0 is called a zero or a root of the function f( x ), or a solution to the equation f( x) 0. You are already
More informationAnnouncements. Midterms graded back at the end of class Help session on Assignment 3 for last ~20 minutes of class. Computer Graphics
Announcements Midterms graded back at the end of class Help session on Assignment 3 for last ~20 minutes of class 1 Scan Conversion Overview of Rendering Scan Conversion Drawing Lines Drawing Polygons
More informationDifferential equations for chemical kinetics
Differential equations for chemical kinetics Atmospheric Modelling course Sampo Smolander sampo.smolander@helsinki.fi 15.11.2010 Sampo Smolander Chemical kinetics 15.11.2010 1 / 24 Chemical equation reactants
More informationData Visualization. Fall 2017
Data Visualization Fall 2017 Vector Fields Vector field v: D R n D is typically 2D planar surface or 2D surface embedded in 3D n = 2 fields tangent to 2D surface n = 3 volumetric fields When visualizing
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationThree-Dimensional Oceanic Flows from Eulerian Velocity Data
Second-year Ph.D. student, Applied Math and Scientific Computing Project Advisor: Kayo Ide Department of Atmospheric and Oceanic Science Center for Scientific Computation and Mathematical Modeling Earth
More informationUSING SPREADSHEETS AND DERIVE TO TEACH DIFFERENTIAL EQUATIONS
USING SPREADSHEETS AND DERIVE TO TEACH DIFFERENTIAL EQUATIONS Kathleen Shannon, Ph.D. Salisbury State University Department of Mathematics and Computer Science Salisbury, MD 21801 KMSHANNON@SAE.SSU.UMD.EDU
More informationdy = dx y dx Project 1 Consider the following differential equations: I. xy II. III.
Project 1 Consider the following differential equations: dy = dy = y dy x + sin( x) y = x I. xy II. III. 1. Graph the direction field in the neighborhood of the origin. Graph approximate integral curves
More informationCONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS
CONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS CHANG YANG AND FRANCIS FILBET Abstract. We develop weighted essentially
More informationNew Very High-Order Upwind Multilayer Compact Schemes with Spectral-Like Resolution for Flow Simulations
New Very High-Order Upwind Multilayer Compact Schemes with Spectral-Lie Resolution for Flow Simulations Zeyu Bai and Xiaolin Zhong University of California, Los Angeles, CA, 995, USA Hypersonic boundary
More information7.1 First-order initial-value problems
7.1 First-order initial-value problems A first-order initial-value problem (IVP) is a first-order ordinary differential equation with a specified value: d dt y t y a y f t, y t That is, this says that
More informationRunge Kutta Methods Optimized For Advection Problems
Runge Kutta Methods Optimized For Advection Problems Brown University December 9, 2009 Introduction Classical Runge Kutta Methods Overview Standard Runge Kutta Schemes for Advection Equation Low Dissipation
More informationSemi-Conservative Schemes for Conservation Laws
Semi-Conservative Schemes for Conservation Laws Rosa Maria Pidatella 1 Gabriella Puppo, 2 Giovanni Russo, 1 Pietro Santagati 1 1 Università degli Studi di Catania, Catania, Italy 2 Università dell Insubria,
More informationComputer Experiments. Designs
Computer Experiments Designs Differences between physical and computer Recall experiments 1. The code is deterministic. There is no random error (measurement error). As a result, no replication is needed.
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationSolving initial value problem by different numerical methods: Practical investigation
Annales Mathematicae et Informaticae 32 (2005) pp. 203 20. Solving initial value problem by different numerical methods: Practical investigation Gábor Geda Department of Computer Science, Eszterházy Károly
More informationNumerical Methods in Physics Lecture 2 Interpolation
Numerical Methods in Physics Pat Scott Department of Physics, Imperial College November 8, 2016 Slides available from http://astro.ic.ac.uk/pscott/ course-webpage-numerical-methods-201617 Outline The problem
More informationComputational Astrophysics 5 Higher-order and AMR schemes
Computational Astrophysics 5 Higher-order and AMR schemes Oscar Agertz Outline - The Godunov Method - Second-order scheme with MUSCL - Slope limiters and TVD schemes - Characteristics tracing and 2D slopes.
More informationAppendix A Numerical Methods
Appendix A Numerical Methods A.1 The Bisection Method In Sect. 5.3 we examined the Newton Raphson method for numerically finding a root of a function. The Newton Raphson method is very effective and fast,
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationFrom the output we see that the largest eigenvalue of B is about (ii) Use the inverse power method to find the smallest eigenvalue.
. Let B be the matrix given by.8635.735.9593.5685.65.735.8984.7439.6447.553.9593.7439.699.635.643.5685.6447.635.056.430.65.553.643.430.5505 (i) Use the power method to find the largest eigenvalue. The
More information9.1 Bracketing and Bisection
350 Chapter 9. Root Finding and Nonlinear Sets of Equations for (i=1;i
More informationModelling and Simulation for Engineers
Unit T7: Modelling and Simulation for Engineers Unit code: F/503/7343 QCF level: 6 Credit value: 15 Aim This unit gives learners the opportunity to develop their understanding of Ordinary Differential
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More information37 Self-Assessment. 38 Two-stage Runge-Kutta Methods. Chapter 10 Runge Kutta Methods
Chapter 10 Runge Kutta Methods In the previous lectures, we have concentrated on multi-step methods. However, another powerful set of methods are known as multi-stage methods. Perhaps the best known of
More informationThe jello cube. Undeformed cube. Deformed cube
The Jello Cube Assignment 1, CSCI 520 Jernej Barbic, USC Undeformed cube The jello cube Deformed cube The jello cube is elastic, Can be bent, stretched, squeezed,, Without external forces, it eventually
More informationInterpolation by Spline Functions
Interpolation by Spline Functions Com S 477/577 Sep 0 007 High-degree polynomials tend to have large oscillations which are not the characteristics of the original data. To yield smooth interpolating curves
More informationPhysics Tutorial 2: Numerical Integration Methods
Physics Tutorial 2: Numerical Integration Methods Summary The concept of numerically solving differential equations is explained, and applied to real time gaming simulation. Some objects are moved in a
More informationFrom Vertices to Fragments: Rasterization. Reading Assignment: Chapter 7. Special memory where pixel colors are stored.
From Vertices to Fragments: Rasterization Reading Assignment: Chapter 7 Frame Buffer Special memory where pixel colors are stored. System Bus CPU Main Memory Graphics Card -- Graphics Processing Unit (GPU)
More informationRecent developments in simulation, optimization and control of flexible multibody systems
Recent developments in simulation, optimization and control of flexible multibody systems Olivier Brüls Department of Aerospace and Mechanical Engineering University of Liège o.bruls@ulg.ac.be Katholieke
More information9.1 Bracketing and Bisection
9.1 Bracketing and Bisection 343 CITED REFERENCES AND FURTHER READING: Stoer, J., and Bulirsch, R. 1980, Introduction to Numerical Analysis (New York: Springer-Verlag), Chapter 5. Acton, F.S. 1970, Numerical
More information