Mathematics for chemical engineers
|
|
- Meryl Webb
- 5 years ago
- Views:
Transcription
1 Mathematics for chemical engineers Drahoslava Janovská Department of mathematics Winter semester
2 Numerical solution of ordinary differential equations Initial value problem
3 Outline 1 Introduction 2 Euler s method Runge - Kutta methods
4 Introduction Numerical solution of differential equations necessity for engineering applications. We will consider a solution of one differential equation as well as a solution of systems of differential equations. Equations of the higher orders can be transformed into a system of differential equations of the first order. In particular, we will study differential equations of the first order, both initial and boundary value problems.
5 Solution of the initial value problem by a one step method We solve the initial value (Cauchy) problem y = f (x, y), y(x 0 ) = y 0. Uniqueness and existence of a solution: If f (x, y) is continuous in Ω = {(x, y), x x 0 a, a > 0, y y 0 b, b > 0} and if we denote M = max f (x, y), α = min(a, b (x,y) Ω M ), then there exist a solution of the equation y = f (x, y) defined in the interval (x 0 α, x 0 + α). Moreover, if f is a Lipschitz function, i.e. f (x, y 1 ) f (x, y 2 ) L y 1 y 2 L > 0, (x, y 1 ), (x, y 2 ) Ω, the solution is unique. Remark: L... Lipschitz constant
6 Taylor expansion of the solution We are looking for a function y = y(t) such that dy dt = y = f (y, t), y(0) = y 0, 0 < t < T. Let us suppose that we have already computed the solution y n := y(t n) in time t n and we want to find a solution in time t n+1. Let h := t n+1 t n be the corresponding time step. Taylor expansion of y n+1 in the point y n: And y n?... etc. y n+1 = y n + h y n }{{} + h2 2 y n + h3 6 y n +... y = = f (y n, t n) from the differential equation ( ) dy = f dt t 1 + f y y = f t + f y f y = (ft + f fy) + t y (ft + f fy) f = ftt + 2f fyt + ft fy + f f y 2 + f 2 f yy,
7 Euler s method Euler s method Euler s method is based on the first two terms of the Taylor expansion: We choose y(0) := y 0 and construct a sequence y n+1 = y n + hf (y n, t n), i.e.,we replace the function y in the point (y n, t n) by a tangent to the integral curve at this point. Accuracy of calculation: Euler s method is the first order method. Remark The step of the method h can be changed with a particular iteration (adaptive choice of step): y n+1 = y n + h nf (y n, t n).
8 Euler s method Example Solve the following initial value problem by Euler s method. Find the solution in time t = 3, y = 0, 3 y sin(t), y(1) = 2. For simplicity, let us consider n = 4, i.e., h = 0.5. y j+1 = y j + h f (y j, t j ), where f (y j, t j ) = 0, 3y j sin(t j ). t j step 0.5 step
9 Euler s method Stability of Euler s method Let us solve the model problem y = λy, λ a constant (1) Euler s method = y n+1 = y n + λh y n, t.j. y n+1 = y n(1 + λ h), i.e., we obtain y n = y n 1 (1 + λ h) = y n 2 (1 + λ h) 2 =... = y 0 (1 + λ h) n. For λ = λ 1 + iλ 2... imaginary, we have y n = y 0 (1 + λ 1 h + iλ 2 h }{{} )n = y 0 σ n. σ... so called amplification factor The numerical solution is stable (i.e. it will remain limited for the growing (large) n, too), if σ 1.
10 Euler s method Let in (1), λ = λ 1 + iλ 2, λ 1 0, σ = 1 + λ 1 h + iλ 2 h. Then the region of stability for Euler s method is part of the left half of the complex plane, in particular inside of the circle σ 2 = (1 + λ 1 h) 2 + λ 2 2h 2 = 1. For any value of λh in the left half of the complex plane outside of this circle the numerical solution is blowing up, while the exact solution decreases. If we want to have a stable solution, we must reduce the h so, that λh would be inside the circle.
11 Euler s method Implicit (backward) Euler s method The y n+1 occurs in the equation implicitly: y n+1 = y n + h f (y n+1, t n+1 ) Disadvantage: the method is computationally more demanding than explicit Euler s method. Advantage: it is more stable, sometimes linearization of f can be exploited with advantage. Example Let us apply the implicit Euler s method to the model problem (1): We obtain: 1 y n+1 = y n + λhy n+1 = y n+1 = yn, i.e., 1 λh ( ) 2 ( ) n y n = 1 λh y n 1 = y n 2 =... = y 0. 1 λh 1 λh y n = σ n y 0, σ = 1 1 λh.
12 Euler s method θ methods We can define the following one-parameter class of one-step methods, called θ methods: For a given initial approximation y 0, we define y n+1 as a convex combination of f (y n, t n) and f (y n+1, t n+1 ) (θ... parameter): y n+1 = y n + h[(1 θ)f (y n, t n) + θf (y n+1, t n+1 )], n = 0, 1,..., N 1, θ 0, 1. θ = 0 = y n+1 = y n + hf (y n, t n)... (explicit) Euler s method θ = 1 = y n+1 = y n + hf (y n+1, t n+1 )... implicit Euler s method θ = 1 2 = y n+1 = y n h[f (yn, tn) + f (y n+1, t n+1 )]... trapezoidal rule. It can be shown that a θ method is explicit for θ = 0 and it is implicit for 0 < θ 1.
13 Runge - Kutta methods Runge - Kutta methods Runge - Kutta methods (RK) more precise than Euler s methods: explicit: The solution in time t n+1 is computed from the values y n, f (y n, t n) and from f (y, t) enumerated at a point between points t n and t n+1 = better accuracy because we use more information about the function f. implicit: They usually lead to the solution of nonlinear algebraic equations, but the amount of work involved is balanced by better numerical stability.
14 Runge - Kutta methods RK method of the second order Let us have the equation y = f (y, t). In the time step t n+1, we obtain the solution from the equation where y n+1 = y n + γ 1 k 1 + γ 2 k 2, (2) k 1 = hf (y n, t n) k 2 = hf (y n + βk 1, t n + αh), α, β, γ 1, γ 2 R. The constants α, β, γ 1, γ 2 have to be determined in such a way that the method should have the highest order of accuracy possible. To determine the order of accuracy, we exploit the Taylor expansion of y(t n+1 ) y n+1 = y n + h y n }{{} + h2 2 y n }{{} + = f (y n, t n) f t + f f y y n+1 = y n + hf (y n, t n) + h2 2 (ftn + f nf yn )+ (3) and compare the coefficients in (2) and (3).
15 Runge - Kutta methods Taylor series for the function of two variables k 2 = hf (y n + βk 1, t n + αh) = ( ) k 2 = h f (y n, t n) + βk 1 f yn + αhf tn + O(h 2 ). Remark Symbol O (capital O) g(h) = O(h p ) g(h) C h p, C is a constant independent of h. Then y n+1 = y n + (γ 1 + γ 2 )hf n + γ 2 βh 2 f nf yn + γ 2 αh 2 f tn + We compare the result with (3) and obtain three nonlinear equations for 4 unknowns: γ 1 + γ 2 = 1, γ 2 α = 1 2, γ 2β = 1 2. Let α R be a parameter. Then γ 2 = 1 2α, β = α, γ 1 = (1 1 2α ).
16 Runge - Kutta methods We obtain Runge - Kutta methods of the order 2: k 1 = hf (y n, t n) k 2 = hf (y n + αk 1, t n + αh) y n+1 = y n + (1 1 2α )k α k 2. We choose α and get the method. For example α = 1 2 = γ 2 = 1, β = 1 2, γ 1 = 0 = y n+1 = y n + k 2 = y n + hf (y n k 1, t n h). Remark RK method of the 2nd order requires at each step twice quantification of function values.
17 Runge - Kutta methods RK methods of 4th order For solving initial value problem, the RK methods of 4th order are mostly used y n+1 = y n k (k 2 + k 3 ) k 4, where k 1 = hf (y n, t n) k 2 = hf (y n k 1, t n + h 2 ) k 3 = hf (y n k 2, t n + h 2 ) k 4 = hf (y n + k 3, t n + h). At each step, we need to compute function values 4 times. Although laborious, the RK method of 4th order is stable and very accurate. It is easily programmable, because it requires no differentiation, only computation of function values.
18 Runge - Kutta methods Computational error Computational error a) y(t n)... exact solution in time t n, y n... approximate solution in time t n e n = y n y(t n)... global error of the approximation b) The computer works in finite arithmetic: r n = ỹ n y n... rounding error We want to compute f but in fact, we compute the numerical approximation ỹ n. For Euler s method: (1)... the global error of approximation e n, E n is directly proportional to the first power of h (2)... rounding error r n is inversely proportional to the first power of h (3)... total error for Euler s method About rounding errors, we can convince only if we repeat the calculation with the different pracission (double precission,... )
Finite 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 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 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 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 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 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 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 informationODEs occur quite often in physics and astrophysics: Wave Equation in 1-D stellar structure equations hydrostatic equation in atmospheres orbits
Solving ODEs 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
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 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 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 information1. Suppose that the equation F (x, y, z) = 0 implicitly defines each of the three variables x, y, and z as functions of the other two:
Final Solutions. Suppose that the equation F (x, y, z) implicitly defines each of the three variables x, y, and z as functions of the other two: z f(x, y), y g(x, z), x h(y, z). If F is differentiable
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 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 information. Tutorial Class V 3-10/10/2012 First Order Partial Derivatives;...
Tutorial Class V 3-10/10/2012 1 First Order Partial Derivatives; Tutorial Class V 3-10/10/2012 1 First Order Partial Derivatives; 2 Application of Gradient; Tutorial Class V 3-10/10/2012 1 First Order
More informationKevin James. MTHSC 206 Section 14.5 The Chain Rule
MTHSC 206 Section 14.5 The Chain Rule Theorem (Chain Rule - Case 1) Suppose that z = f (x, y) is a differentiable function and that x(t) and y(t) are both differentiable functions as well. Then, dz dt
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 informationIntroduction to Optimization Problems and Methods
Introduction to Optimization Problems and Methods wjch@umich.edu December 10, 2009 Outline 1 Linear Optimization Problem Simplex Method 2 3 Cutting Plane Method 4 Discrete Dynamic Programming Problem Simplex
More informationREVIEW I MATH 254 Calculus IV. Exam I (Friday, April 29) will cover sections
REVIEW I MATH 254 Calculus IV Exam I (Friday, April 29 will cover sections 14.1-8. 1. Functions of multivariables The definition of multivariable functions is similar to that of functions of one variable.
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 6:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_1.htm 1 of 1 6/20/2012 12:24 PM The Lecture deals with: ADI Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_2.htm 1 of 2 6/20/2012
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 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 informationDD2429 Computational Photography :00-19:00
. Examination: DD2429 Computational Photography 202-0-8 4:00-9:00 Each problem gives max 5 points. In order to pass you need about 0-5 points. You are allowed to use the lecture notes and standard list
More informationMath 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 informationMAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations)
MATLAB sessions: Laboratory 3 1 MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations) In this session we look at basic numerical methods to help us understand
More informationInverse and Implicit functions
CHAPTER 3 Inverse and Implicit functions. Inverse Functions and Coordinate Changes Let U R d be a domain. Theorem. (Inverse function theorem). If ϕ : U R d is differentiable at a and Dϕ a is invertible,
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 informationFunctions of Two variables.
Functions of Two variables. Ferdinánd Filip filip.ferdinand@bgk.uni-obuda.hu siva.banki.hu/jegyzetek 27 February 217 Ferdinánd Filip 27 February 217 Functions of Two variables. 1 / 36 Table of contents
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 informationMAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations)
MAT 275 Laboratory 3 Numerical Solutions by Euler and Improved Euler Methods (scalar equations) In this session we look at basic numerical methods to help us understand the fundamentals of numerical approximations.
More informationMath 126 Winter CHECK that your exam contains 8 problems.
Math 126 Winter 2016 Your Name Your Signature Student ID # Quiz Section Professor s Name TA s Name CHECK that your exam contains 8 problems. This exam is closed book. You may use one 8 1 11 sheet of hand-written
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 informationLecture 15. Lecturer: Prof. Sergei Fedotov Calculus and Vectors. Length of a Curve and Parametric Equations
Lecture 15 Lecturer: Prof. Sergei Fedotov 10131 - Calculus and Vectors Length of a Curve and Parametric Equations Sergei Fedotov (University of Manchester) MATH10131 2011 1 / 5 Lecture 15 1 Length of a
More informationMathematical Methods and Modeling Laboratory class. Numerical Integration of Ordinary Differential Equations
Mathematical Methods and Modeling Laboratory class Numerical Integration of Ordinary Differential Equations Exact Solutions of ODEs Cauchy s Initial Value Problem in normal form: Recall: if f is locally
More informationContinuity and Tangent Lines for functions of two variables
Continuity and Tangent Lines for functions of two variables James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University April 4, 2014 Outline 1 Continuity
More informationBackground for Surface Integration
Background for urface Integration 1 urface Integrals We have seen in previous work how to define and compute line integrals in R 2. You should remember the basic surface integrals that we will need to
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 informationSurfaces and Integral Curves
MODULE 1: MATHEMATICAL PRELIMINARIES 16 Lecture 3 Surfaces and Integral Curves In Lecture 3, we recall some geometrical concepts that are essential for understanding the nature of solutions of partial
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 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 informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationParametric and Polar Curves
Chapter 2 Parametric and Polar Curves 2.1 Parametric Equations; Tangent Lines and Arc Length for Parametric Curves Parametric Equations So far we ve described a curve by giving an equation that the coordinates
More informationParametric and Polar Curves
Chapter 2 Parametric and Polar Curves 2.1 Parametric Equations; Tangent Lines and Arc Length for Parametric Curves Parametric Equations So far we ve described a curve by giving an equation that the coordinates
More informationParametric and Polar Curves
Chapter 2 Parametric and Polar Curves 2.1 Parametric Equations; Tangent Lines and Arc Length for Parametric Curves Parametric Equations So far we ve described a curve by giving an equation that the coordinates
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationNumerical Solutions of Differential Equations (1)
Numerical Solutions of Differential Equations (1) Euler method Sources of error Truncation Round-off error Error balance in numerical integration Euler Method Many differential equations don't have an
More informationWinfried Auzinger, Dirk Praetorius (SS2012)
Computer Algebra using Maple Part V: Extensions, some special topics. Didactical aspects Winfried Auzinger, Dirk Praetorius (SS2012) restart: 1 Polynomial equations and systems For polynomial equations
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 information= w. w u. u ; u + w. x x. z z. y y. v + w. . Remark. The formula stated above is very important in the theory of. surface integral.
1 Chain rules 2 Directional derivative 3 Gradient Vector Field 4 Most Rapid Increase 5 Implicit Function Theorem, Implicit Differentiation 6 Lagrange Multiplier 7 Second Derivative Test Theorem Suppose
More information. 1. Chain rules. Directional derivative. Gradient Vector Field. Most Rapid Increase. Implicit Function Theorem, Implicit Differentiation
1 Chain rules 2 Directional derivative 3 Gradient Vector Field 4 Most Rapid Increase 5 Implicit Function Theorem, Implicit Differentiation 6 Lagrange Multiplier 7 Second Derivative Test Theorem Suppose
More informationConvex Optimization Lecture 2
Convex Optimization Lecture 2 Today: Convex Analysis Center-of-mass Algorithm 1 Convex Analysis Convex Sets Definition: A set C R n is convex if for all x, y C and all 0 λ 1, λx + (1 λ)y C Operations that
More informationChapter 6. Curves and Surfaces. 6.1 Graphs as Surfaces
Chapter 6 Curves and Surfaces In Chapter 2 a plane is defined as the zero set of a linear function in R 3. It is expected a surface is the zero set of a differentiable function in R n. To motivate, graphs
More informationJim Lambers MAT 169 Fall Semester Lecture 33 Notes
Jim Lambers MAT 169 Fall Semester 2009-10 Lecture 33 Notes These notes correspond to Section 9.3 in the text. Polar Coordinates Throughout this course, we have denoted a point in the plane by an ordered
More informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationHomework Set 5 (Sections )
Homework Set 5 (Sections.4-.6) 1. Consider the initial value problem (IVP): = 8xx yy, yy(1) = 3 a. Do 10 steps of Euler s Method, using a step size of 0.1. Write the details in the table below. Work to
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 information21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane
12th International Conference on Domain Decomposition Methods Editors: Tony Chan, Takashi Kako, Hideo Kawarada, Olivier Pironneau, c 2001 DDM.org 21. Efficient and fast numerical methods to compute fluid
More informationApplication 7.6A The Runge-Kutta Method for 2-Dimensional Systems
Application 7.6A The Runge-Kutta Method for -Dimensional Systems Figure 7.6. in the text lists TI-85 and BASIC versions of the program RKDIM that implements the Runge-Kutta iteration k (,, ) = f tn xn
More informationEuler s Method and Logistic Growth (BC Only)
Euler s Method Students should be able to: Approximate numerical solutions of differential equations using Euler s method without a calculator. Recognize the method as an recursion formula extension of
More informationCurve and Surface Basics
Curve and Surface Basics Implicit and parametric forms Power basis form Bezier curves Rational Bezier Curves Tensor Product Surfaces ME525x NURBS Curve and Surface Modeling Page 1 Implicit and Parametric
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 informationNumerical Methods on the Image Processing Problems
Numerical Methods on the Image Processing Problems Department of Mathematics and Statistics Mississippi State University December 13, 2006 Objective Develop efficient PDE (partial differential equations)
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 informationA Grid Based Particle Method for Evolution of Open Curves and Surfaces
A Grid Based Particle Method for Evolution of Open Curves and Surfaces Shingyu Leung a,, Hongkai Zhao b a Department of Mathematics, Hong Kong University of Science and Technology, Clear Water Bay, Hong
More informationTruncation Errors. Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4.
Chapter 4: Roundoff and Truncation Errors Applied Numerical Methods with MATLAB for Engineers and Scientists, 2nd ed., Steven C. Chapra, McGraw Hill, 2008, Ch. 4. 1 Outline Errors Accuracy and Precision
More informationDaily WeBWorK, #1. This means the two planes normal vectors must be multiples of each other.
Daily WeBWorK, #1 Consider the ellipsoid x 2 + 3y 2 + z 2 = 11. Find all the points where the tangent plane to this ellipsoid is parallel to the plane 2x + 3y + 2z = 0. In order for the plane tangent to
More informationSparse Optimization Lecture: Proximal Operator/Algorithm and Lagrange Dual
Sparse Optimization Lecture: Proximal Operator/Algorithm and Lagrange Dual Instructor: Wotao Yin July 2013 online discussions on piazza.com Those who complete this lecture will know learn the proximal
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 informationHartley - Zisserman reading club. Part I: Hartley and Zisserman Appendix 6: Part II: Zhengyou Zhang: Presented by Daniel Fontijne
Hartley - Zisserman reading club Part I: Hartley and Zisserman Appendix 6: Iterative estimation methods Part II: Zhengyou Zhang: A Flexible New Technique for Camera Calibration Presented by Daniel Fontijne
More informationMath (Spring 2009): Lecture 5 Planes. Parametric equations of curves and lines
Math 18.02 (Spring 2009): Lecture 5 Planes. Parametric equations of curves and lines February 12 Reading Material: From Simmons: 17.1 and 17.2. Last time: Square Systems. Word problem. How many solutions?
More informationMath 213 Exam 2. Each question is followed by a space to write your answer. Please write your answer neatly in the space provided.
Math 213 Exam 2 Name: Section: Do not remove this answer page you will return the whole exam. You will be allowed two hours to complete this test. No books or notes may be used other than a onepage cheat
More informationEquation of tangent plane: for implicitly defined surfaces section 12.9
Equation of tangent plane: for implicitly defined surfaces section 12.9 Some surfaces are defined implicitly, such as the sphere x 2 + y 2 + z 2 = 1. In general an implicitly defined surface has the equation
More informationPARALLEL METHODS FOR SOLVING PARTIAL DIFFERENTIAL EQUATIONS. Ioana Chiorean
5 Kragujevac J. Math. 25 (2003) 5 18. PARALLEL METHODS FOR SOLVING PARTIAL DIFFERENTIAL EQUATIONS Ioana Chiorean Babeş-Bolyai University, Department of Mathematics, Cluj-Napoca, Romania (Received May 28,
More informationMath 407 Solutions to Exam 3 April 24, Let f = u + iv, and let Γ be parmetrized by Γ = (x (t),y(t)) = x (t) +iy (t) for a t b.
Math 7 Solutions to Exam 3 April, 3 1. (1 Define the following: (a f (z dz, where Γ is a smooth curve with finite length. Γ Let f u + iv, and let Γ be parmetrized by Γ (x (t,y(t x (t +iy (t for a t b.
More informationChapter 6 Some Applications of the Integral
Chapter 6 Some Applications of the Integral More on Area More on Area Integrating the vertical separation gives Riemann Sums of the form More on Area Example Find the area A of the set shaded in Figure
More informationLecture 8: Jointly distributed random variables
Lecture : Jointly distributed random variables Random Vectors and Joint Probability Distributions Definition: Random Vector. An n-dimensional random vector, denoted as Z = (Z, Z,, Z n ), is a function
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 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 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 informationMAT203 OVERVIEW OF CONTENTS AND SAMPLE PROBLEMS
MAT203 OVERVIEW OF CONTENTS AND SAMPLE PROBLEMS MAT203 covers essentially the same material as MAT201, but is more in depth and theoretical. Exam problems are often more sophisticated in scope and difficulty
More informationUpdated: August 24, 2016 Calculus III Section Math 232. Calculus III. Brian Veitch Fall 2015 Northern Illinois University
Updated: August 24, 216 Calculus III Section 1.2 Math 232 Calculus III Brian Veitch Fall 215 Northern Illinois University 1.2 Calculus with Parametric Curves Definition 1: First Derivative of a Parametric
More informationMath 209, Fall 2009 Homework 3
Math 209, Fall 2009 Homework 3 () Find equations of the tangent plane and the normal line to the given surface at the specified point: x 2 + 2y 2 3z 2 = 3, P (2,, ). Solution Using implicit differentiation
More informationMA 114 Worksheet #17: Average value of a function
Spring 2019 MA 114 Worksheet 17 Thursday, 7 March 2019 MA 114 Worksheet #17: Average value of a function 1. Write down the equation for the average value of an integrable function f(x) on [a, b]. 2. Find
More informationTangent & normal vectors The arc-length parameterization of a 2D curve The curvature of a 2D curve
Local Analysis of 2D Curve Patches Topic 4.2: Local analysis of 2D curve patches Representing 2D image curves Estimating differential properties of 2D curves Tangent & normal vectors The arc-length parameterization
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 informationMath 348 Differential Geometry of Curves and Surfaces
Math 348 Differential Geometry of Curves and Surfaces Lecture 3 Curves in Calculus Xinwei Yu Sept. 12, 2017 CAB 527, xinwei2@ualberta.ca Department of Mathematical & Statistical Sciences University of
More informationOpenGL Graphics System. 2D Graphics Primitives. Drawing 2D Graphics Primitives. 2D Graphics Primitives. Mathematical 2D Primitives.
D Graphics Primitives Eye sees Displays - CRT/LCD Frame buffer - Addressable pixel array (D) Graphics processor s main function is to map application model (D) by projection on to D primitives: points,
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 informationSpline Curves. Spline Curves. Prof. Dr. Hans Hagen Algorithmic Geometry WS 2013/2014 1
Spline Curves Prof. Dr. Hans Hagen Algorithmic Geometry WS 2013/2014 1 Problem: In the previous chapter, we have seen that interpolating polynomials, especially those of high degree, tend to produce strong
More informationChapter One. Numerical Algorithms
Chapter One Numerical Algorithms The word algorithm derives from the name of the Persian mathematician (Abu Ja far Muhammad ibn Musa) Al- Khwarizmi who lived from about 790 CE to about 840 CE. He wrote
More informationLAB #8 Numerical Methods
LAB #8 Numerical Methods Goal: The purpose of this lab is to explain how computers numerically approximate solutions to differential equations. Required tools: Matlab routine dfield ; numerical routines
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 information7/28/2011 SECOND HOURLY PRACTICE V Maths 21a, O.Knill, Summer 2011
7/28/2011 SECOND HOURLY PRACTICE V Maths 21a, O.Knill, Summer 2011 Name: Start by printing your name in the above box. Try to answer each question on the same page as the question is asked. If needed,
More informationHyperbola for Curvilinear Interpolation
Applied Mathematical Sciences, Vol. 7, 2013, no. 30, 1477-1481 HIKARI Ltd, www.m-hikari.com Hyperbola for Curvilinear Interpolation G. L. Silver 868 Kristi Lane Los Alamos, NM 87544, USA gsilver@aol.com
More informationKey points. Assume (except for point 4) f : R n R is twice continuously differentiable. strictly above the graph of f except at x
Key points Assume (except for point 4) f : R n R is twice continuously differentiable 1 If Hf is neg def at x, then f attains a strict local max at x iff f(x) = 0 In (1), replace Hf(x) negative definite
More informationCS321 Introduction To Numerical Methods
CS3 Introduction To Numerical Methods Fuhua (Frank) Cheng Department of Computer Science University of Kentucky Lexington KY 456-46 - - Table of Contents Errors and Number Representations 3 Error Types
More informationNumerical Analysis Fall. Numerical Differentiation
Numerical Analysis 5 Fall Numerical Differentiation Differentiation The mathematical definition of a derivative begins with a difference approimation: and as is allowed to approach zero, the difference
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 information(1) Given the following system of linear equations, which depends on a parameter a R, 3x y + 5z = 2 4x + y + (a 2 14)z = a + 2
(1 Given the following system of linear equations, which depends on a parameter a R, x + 2y 3z = 4 3x y + 5z = 2 4x + y + (a 2 14z = a + 2 (a Classify the system of equations depending on the values of
More information