Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
|
|
- Baldwin Arnold
- 6 years ago
- Views:
Transcription
1 Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler forward, Euler backward, the Theta method, and their stability Dr. Noemi Friedman,
2 Overview of the course Introduction (definition of PDEs, classification, basic math, introductory examples of PDEs) Analytical solution of elementary PDEs (Fourier series/transform, separation of variables, Green s function) Numerical solutions of PDEs: Finite difference method Finite element method PDE lecture Seite 2
3 Overview of this lecture Finite difference operators The heat equation Analytical solution Semidiscretization- Method of Lines (spatial discretization) Time discretization Euler foward Euler backward The theta-method Consisteny, stability, convergence PDE lecture Seite 3
4 We only analyse now the hom. equation - why can we do that? d u = Au(t) + f dt Instead of analyzing stability of the inhomogenous case, we discretize the homogenous one. We can do that, because of the following. Let s take a stationary function u 0 for which the equation: d dt u 0 = Au 0 + f = 0 holds. And let s suppose, that the solution can be written in the form u t = u 0 + v(t) r.h.s: Au t + f = Au 0 + f + Av t = Av t l.h.s: = 0 d dt u t = d dt u 0 + v t = d v t dt d v t = Av t dt Dr. Noemi Friedman PDE lecture Seite 4
5 Heat equation Comparism of the solutions: analytical and method of lines Homogeneous equation (no internal source term): d u = Au(t) dt Analytical solution with hom. BC: u x, t = j d j e ω j 2 βt sin ω j x ω j = jπ l j = 0,,2.. as t u 0 With method of lines: u t = i c i e λ i(t t 0 ) v i λ i = 2β2 h 2 cos iπ N > 0 < as t u Dr. Noemi Friedman PDE lecture Seite 5
6 2.) Euler forward - explicit Euler forward: approximate time derivate with the forward difference u j (t) t = β2 h 2 u j 2u j + u j+ + O(h 2 ) u j,n = u j,n+ u j,n + O() u n = u n+ u n = Au n u j,n+ u j,n = β2 h 2 u j,n 2u j,n + u j+,n + O h 2 + O() Stable? (does it give decaying solution?) B = β2 h Dr. Noemi Friedman PDE lecture Seite 6
7 2.) Euler forward d-tour on eigenvalues and eigenvectors Stable? (does it give decaying solution?) u n+ = Bu n u n+2 = B 2 u n u n+3 = B 3 u n u n+k = B k u n Some discussion about egienvalues and eigenvectors: Av i = λ i v i AV = DV V AV = D V AV is a diagonal matrix, with the eigenvalues in the diagonal. D = V = λ λ 2 λ n λ n v v 2 v n v n Dr. Noemi Friedman PDE lecture Seite 7
8 Eigenvalues, eigenvectors change of basis (coordinate system) How to write an arbitrary vector, u, in a new basis v (), v (2), v (3) What is V AV? u = u e + u 2 e 2 + u 3 e 3 u = u v () + u 2 v (2) + u 3 v (3) v () () () () = v e + v 2 e2 + v 3 e3 v (2) (2) (2) (2) = v e + v 2 e2 + v 3 e3 v (3) (3) (3) (3) = v e + v 2 e2 + v 3 e3 u = u () () () v e + v 2 e2 + v 3 e3 + + u 2 (2) (2) (2) v e + v 2 e2 + v 3 e3 + + u 3 (3) (3) (3) v e + v 2 e2 + v 3 e3 + v (3) v (2) e 2 u e e 3 u v () u = e u v () + u2 v (2) + u3 v (3) u 2 e 2 () (2) (3) u v 2 + u2 v 2 + u3 v 2 e 3 () (2) (3) u v 3 + u2 v 3 + u3 v Dr. Noemi Friedman PDE lecture Seite 8 u 3
9 Eigenvalues, eigenvectors change of basis (coordinate system) Some discussion about egienvalues and eigenvectors: v () v (2) v (3) u u u 2 u 3 = v () v 2 () v 3 () v (2) v 2 (2) v 3 (2) v (3) v 2 (3) v 3 (2) u u 2 u 3 u = e u v () + u2 v (2) + u3 v (3) u 2 e 2 () (2) (3) u v 2 + u2 v 2 + u3 v 2 e 3 () (2) (3) u v 3 + u2 v 3 + u3 v u = V u u = V u u 3 Let s consider now the respresentation of an operator in the new basis: Au = b AV u = V b V AV is the representationa of the A operator in the new basis defined by v (), v (2), v (3) V AV u = V V b V AV u = b Dr. Noemi Friedman PDE lecture Seite 9
10 2.) Euler forward d-tour on eigenvalues and eigenvectors V AV is the representationa of the A operator in the new basis defined by v (), v (2), v (3) If v (), v (2), v (3) are the eigenvectors of A, this representation is a diagonal matrix. Let s go back to the problem of the Euler forward method used for solving the heat equation Stable? (does it give decaying solution?) u n+ = Bu n u n+2 = B 2 u n u n+k = B k u n Let s check instead in the basis defined by the eigenvectors: u i = V u i u n+ = V BV u n u n+2 = V B 2 V u n u n+k = V B k V u n V B k V = V BB BBV = V BVV BVV VV BVV BV = D k = λ k λ 2 k D D λ n k Dr. Noemi Friedman PDE lecture Seite 0
11 2.) Euler forward stability analysis Stable? (does it give decaying solution?) B = β2 h When for all the eigenvalues λ i < then response is decaying with time. β2 a = 2 h 2 λ i = a + 2b cos iπ N b = β2 h 2 = 2 β2 h 2 cos iπ N 2 β2 h 2 cos iπ N < i =.. N Dr. Noemi Friedman PDE lecture Seite
12 2.) Euler forward stability analysis Stable? (does it give decaying solution?) When for all the eigenvalues λ i 2 β2 h 2 cos iπ N < i =.. N < then response is decaying with time. 2 β2 h 2 cos iπ N < 2 β2 h 2 cos iπ N < allways satisfied < 2 β2 h 2 cos iπ N 0< 2 Max value: 2 2 β2 h 2 cos iπ N > 2 β2 h 2 cos iπ N < 2 cos iπ N < h2 β 2 2 < h2 β 2 > Dr. Noemi Friedman PDE lecture Seite 2 < h2 2β 2 Method is not unconditionally stable, only stable if criteria is satisfied!
13 2.) Euler forward summary Euler forward: approximate time derivate with the forward difference u j (t) t = β2 h 2 u j 2u j + u j+ + O(h 2 ) u j,n = u j,n+ u j,n + O() u n = u n+ u n = Au n u j,n+ u j,n = β2 h 2 u j,n 2u j,n + u j+,n + O h 2 + O() Stable? (does it give decaying solution?) the absolut values of the eigenvalues of matrix B can not be greater than one. λ j = 2 β2 h 2 cos N π N λ j < h2 2β Dr. Noemi Friedman PDE lecture Seite 3
14 3.) Euler backward (implicit) derivation Euler backward: approximate time derivate with the backward difference u j (t) t = β2 h 2 u j 2u j + u j+ + O(h 2 ) u j,n+ = u j,n+ u j,n + O() u n+ = u n+ u n = Au n+ u j,n+ u j,n = β2 h 2 u j,n+ 2u j,n+ + u j+,n+ + O h 2 + O() u n = u n+ Au n+ = (I A)u n+ (Stability criteria: unconditionaly stable, to be shown later) Dr. Noemi Friedman PDE lecture Seite 4
15 4.) Theta method (implicit) derivation Theta method (Crank-Nicolson) u j (t) t = β2 h 2 u j 2u j + u j+ + O(h 2 ) u j,n+θ = u j,n+ u j,n + O( p ) u j,n = u j,n+ u j,n u n = u n+ u n + O() = Au n Euler f. u j,n+ u j,n = θ β2 h 2 u j,n+ 2u j,n+ + u j+,n+ + + θ β2 h 2 u j,n 2u j,n + u j+,n + O h 2 + O( p ) u j,n+ = u j,n+ u j,n u n+ = u n+ u n + O() = Au n+ Euler b. u n+θ = u n+ u n = θau n+ + θ Au n (I θa)u n+ =(I + θ A)u n B B 2 θ = 0 Euler forward θ = Euler backward Dr. Noemi Friedman PDE lecture Seite 5
16 4.) Theta method (implicit) stability analysis Stability criteria: B = (I θa) B u n+ = B 2 u n B 2 = (I + θ A) ) B and B 2 are both tridiagonal sym. matrices B and B 2 have the same eigenvectors v i = v i v i N v i k = sin ikπ N 2) The representation in the basis defined by the eigenvectors V B V u n+ = V B 2 V u n D λ (B ) λ 2 (B ) D 2 λ (B 2 ) λ 2 (B 2 ) u n+ = u n λ n (B ) λ n (B 2 ) Dr. Noemi Friedman PDE lecture Seite 6
17 4.) Theta method (implicit) stability analysis Stability criteria: B u n+ = B 2 u n B = (I θa) a b u n+ = λ (B 2 ) λ (B ) λ 2 (B 2 ) λ 2 (B ) u n B = b a b b a λ n (B 2 ) λ n (B ) B 2 = (I + θ A) a 2 b 2 λ i B = a + 2b cos iπ N = + 2β2 h 2 θ cos iπ N λ i B 2 = a 2 + 2b 2 cos iπ N = 2β2 h 2 θ cos iπ N B 2 = b 2 a 2 b 2 b 2 a Dr. Noemi Friedman PDE lecture Seite 7
18 u n+ = Numerical solution of the heat equation 4.) Theta method (implicit) stability analysis λ (B 2 ) λ (B ) λ 2 (B 2 ) λ 2 (B ) λ n (B 2 ) λ n (B ) u n 2r 0< C i 2 2β λ i (B 2 ) 2 λ i (B ) = h 2 θ cos iπ N + 2β2 h 2 θ cos iπ N r β2 h 2 2r 0< C i 2 < < 2r θ C i + 2rθC i < C i cos iπ N Dr. Noemi Friedman PDE lecture Seite 8
19 4.) Theta method (implicit) stability analysis ) Right hand side of the equation: < 2r θ C i + 2rθC i > 0 2rθC i < 2r θ C i / 2 2rθC i < 2rC i + 2rθC i /+2rθC i 2 < 2rC i + 4rθC i > r 2θ C i /: ( 2) r β2 h 2 2) Left hand side of the equation: 2r θ C i + 2rθC i < > 0 2r θ C i < + 2rθC i 2rC i > 0 unconditionally satisfied 2θ 0 (θ 0.5) 2θ > 0 (θ < 0.5) > r 2θ C i < 0 unconditionally satisfied C i 2θ > r h 2 β 2 C i 2θ > max 2 max Dr. Noemi Friedman PDE lecture Seite 9 h 2 2β 2 2θ >
20 Summary: Instationary heat equation what to solve? Stability checking from eigenvalue analysis: Method of lines Euler forward method find the eigenvalues (λ j ) and eigenvectors (v j ) of matrix A u n+ = I + A B u n u n+ = Bu n Euler backward method u n = I A B u n+ Theta method B u n+ = u n Solve system of equations for u n+ I θa u n+ = I + θ A u n B θ u n+ = B 2θ u n B θ B 2θ Dr. Noemi Friedman PDE lecture Seite 20
Introduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers Dr. Noemi
More informationIntroduction to PDEs and Numerical Methods Lecture 7. Solving linear systems
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 7. Solving linear systems Dr. Noemi Friedman, 09.2.205. Reminder: Instationary heat
More informationIntroduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence Dr. Noemi Friedman,
More informationIntroduction to PDEs and Numerical Methods Tutorial 1: Overview of essential linear algebra and analysis
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 1: Overview of essential linear algebra and analysis Dr. Noemi Friedman, 25.10.201.
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen PARAMetric UNCertainties, Budapest STOCHASTIC PROCESSES AND FIELDS Noémi Friedman Institut für Wissenschaftliches Rechnen, wire@tu-bs.de
More informationIntroduction to PDEs and Numerical Methods Tutorial 11. 2D elliptic equations
Platzalter für Bild, Bild auf Titelfolie inter das Logo einsetzen Introduction to PDEs and Numerical Metods Tutorial 11. 2D elliptic equations Dr. Noemi Friedman, 3. 1. 215. Overview Introduction (classification
More informationIntroduction to PDEs and Numerical Methods Lecture 1: Introduction
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 1: Introduction Dr. Noemi Friedman, 28.10.2015. Basic information on the course Course
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationIntroduction to PDEs and Numerical Methods Tutorial 10. Finite Element Analysis
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Introduction to PDEs and Numerica Methods Tutoria. Finite Eement Anaysis Dr. Noemi Friedman, 3..2 FROM STRONG FORM TO WEAK FORM inhomogeneous
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationChapter 18. Remarks on partial differential equations
Chapter 8. Remarks on partial differential equations If we try to analyze heat flow or vibration in a continuous system such as a building or an airplane, we arrive at a kind of infinite system of ordinary
More informationSTABILITY FOR PARABOLIC SOLVERS
Review STABILITY FOR PARABOLIC SOLVERS School of Mathematics Semester 1 2008 OUTLINE Review 1 REVIEW 2 STABILITY: EXPLICIT METHOD Explicit Method as a Matrix Equation Growing Errors Stability Constraint
More informationMath 3191 Applied Linear Algebra
Math 9 Applied Linear Algebra Lecture 9: Diagonalization Stephen Billups University of Colorado at Denver Math 9Applied Linear Algebra p./9 Section. Diagonalization The goal here is to develop a useful
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 informationNumerical methods for PDEs FEM - abstract formulation, the Galerkin method
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Numerica methods for PDEs FEM - abstract formuation, the Gaerkin method Dr. Noemi Friedman Contents of the course Fundamentas of functiona
More informationNear-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Near-wall Reynolds stress modelling for RANS and hybrid RANS/LES methods Axel Probst (now at: C 2 A 2 S 2 E, DLR Göttingen) René Cécora,
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationUNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH 2 IMPLICIT-EXPLICIT NUMERICAL METHOD
UNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH IMPLICIT-EXPLICIT NUMERICAL METHOD ANDREW JORGENSON Abstract. Systems of nonlinear partial differential equations modeling turbulent fluid flow
More informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationDepartment of Mathematics California State University, Los Angeles Master s Degree Comprehensive Examination in. NUMERICAL ANALYSIS Spring 2015
Department of Mathematics California State University, Los Angeles Master s Degree Comprehensive Examination in NUMERICAL ANALYSIS Spring 2015 Instructions: Do exactly two problems from Part A AND two
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationHomework 3 Solutions Math 309, Fall 2015
Homework 3 Solutions Math 39, Fall 25 782 One easily checks that the only eigenvalue of the coefficient matrix is λ To find the associated eigenvector, we have 4 2 v v 8 4 (up to scalar multiplication)
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationOn Schemes for Exponential Decay
On Schemes for Exponential Decay Hans Petter Langtangen 1,2 Center for Biomedical Computing, Simula Research Laboratory 1 Department of Informatics, University of Oslo 2 Sep 24, 2015 1.0 0.8 numerical
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationPartial Differential Equations II
Partial Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Partial Differential Equations II 1 / 28 Almost Done! Homework
More information1.6: 16, 20, 24, 27, 28
.6: 6, 2, 24, 27, 28 6) If A is positive definite, then A is positive definite. The proof of the above statement can easily be shown for the following 2 2 matrix, a b A = b c If that matrix is positive
More informationSpace-time kinetics. Alain Hébert. Institut de génie nucléaire École Polytechnique de Montréal.
Space-time kinetics Alain Hébert alain.hebert@polymtl.ca Institut de génie nucléaire École Polytechnique de Montréal ENE613: Week 1 Space-time kinetics 1/24 Content (week 1) 1 Space-time kinetics equations
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationFinite Difference Methods (FDMs) 2
Finite Difference Methods (FDMs) 2 Time- dependent PDEs A partial differential equation of the form (15.1) where A, B, and C are constants, is called quasilinear. There are three types of quasilinear equations:
More informationUNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH 2 NUMERICAL METHOD
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 5, Number 3, Pages 7 87 c 04 Institute for Scientific Computing and Information UNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More information1 Last time: least-squares problems
MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that
More informationJordan Canonical Form
Jordan Canonical Form Massoud Malek Jordan normal form or Jordan canonical form (named in honor of Camille Jordan) shows that by changing the basis, a given square matrix M can be transformed into a certain
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationProperties of Linear Transformations from R n to R m
Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
More informationNumerical Solution of Initial Boundary Value Problems. Jan Nordström Division of Computational Mathematics Department of Mathematics
Numerical Solution of Initial Boundary Value Problems Jan Nordström Division of Computational Mathematics Department of Mathematics Overview Material: Notes + GUS + GKO + HANDOUTS Schedule: 5 lectures
More informationMatrices related to linear transformations
Math 4326 Fall 207 Matrices related to linear transformations We have encountered several ways in which matrices relate to linear transformations. In this note, I summarize the important facts and formulas
More informationLecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems
Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now
More informationHTS Cable Integration into Rural Networks with Renewable Energy Resources
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen HTS Cable Integration into Rural Networks with Renewable Energy Resources Dr.-Ing. Nasser Hemdan Outline Introduction Objectives Network
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationProof by induction ME 8
Proof by induction ME 8 n Let f ( n) 9, where n. f () 9 8, which is divisible by 8. f ( n) is divisible by 8 when n =. Assume that for n =, f ( ) 9 is divisible by 8 for. f ( ) 9 9.9 9(9 ) f ( ) f ( )
More informationBeam Propagation Method Solution to the Seminar Tasks
Beam Propagation Method Solution to the Seminar Tasks Matthias Zilk The task was to implement a 1D beam propagation method (BPM) that solves the equation z v(xz) = i 2 [ 2k x 2 + (x) k 2 ik2 v(x, z) =
More informationStudy guide: Generalizations of exponential decay models. Extension to a variable coecient; Crank-Nicolson
Extension to a variable coecient; Forward and Backward Euler Study guide: Generalizations of exponential decay models Hans Petter Langtangen 1,2 Center for Biomedical Computing, Simula Research Laboratory
More informationMath 211. Lecture #6. Linear Equations. September 9, 2002
1 Math 211 Lecture #6 Linear Equations September 9, 2002 2 Air Resistance 2 Air Resistance Acts in the direction opposite to the velocity. 2 Air Resistance Acts in the direction opposite to the velocity.
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationStudy guide: Generalizations of exponential decay models
Study guide: Generalizations of exponential decay models Hans Petter Langtangen 1,2 Center for Biomedical Computing, Simula Research Laboratory 1 Department of Informatics, University of Oslo 2 Sep 13,
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More information15 n=0. zz = re jθ re jθ = r 2. (b) For division and multiplication, it is handy to use the polar representation: z = rejθ. = z 1z 2.
Professor Fearing EECS0/Problem Set v.0 Fall 06 Due at 4 pm, Fri. Sep. in HW box under stairs (st floor Cory) Reading: EE6AB notes. This problem set should be review of material from EE6AB. (Please note,
More informationUNCONDITIONAL STABILITY OF A CRANK-NICOLSON/ADAMS-BASHFORTH 2 IMPLICIT/EXPLICIT METHOD FOR ORDINARY DIFFERENTIAL EQUATIONS
UNCONDITIONAL STABILITY OF A CRANK-NICOLSON/ADAMS-BASHFORTH IMPLICIT/EXPLICIT METHOD FOR ORDINARY DIFFERENTIAL EQUATIONS by Andrew D. Jorgenson B.S., Gonzaga University, 009 M.A., University of Pittsburgh,
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationProperties of Transformations
6. - 6.4 Properties of Transformations P. Danziger Transformations from R n R m. General Transformations A general transformation maps vectors in R n to vectors in R m. We write T : R n R m to indicate
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationGalerkin Finite Element Model for Heat Transfer
Galerkin Finite Element Model for Heat Transfer Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents 1 Notation remarks 1 2 Local differential
More information3D Coordinate Transformations. Tuesday September 8 th 2015
3D Coordinate Transformations Tuesday September 8 th 25 CS 4 Ross Beveridge & Bruce Draper Questions / Practice (from last week I messed up!) Write a matrix to rotate a set of 2D points about the origin
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationInstitute of Dynamics and Vibrations Field of Expertise, Tribology
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Institute of Dynamics and Vibrations Field of Expertise, Tribology Prof. Dr.-Ing. habil. G.-P. Ostermeyer Our Approach to Solving Your
More informationRemark By definition, an eigenvector must be a nonzero vector, but eigenvalue could be zero.
Sec 6 Eigenvalues and Eigenvectors Definition An eigenvector of an n n matrix A is a nonzero vector x such that A x λ x for some scalar λ A scalar λ is called an eigenvalue of A if there is a nontrivial
More informationMath113: Linear Algebra. Beifang Chen
Math3: Linear Algebra Beifang Chen Spring 26 Contents Systems of Linear Equations 3 Systems of Linear Equations 3 Linear Systems 3 2 Geometric Interpretation 3 3 Matrices of Linear Systems 4 4 Elementary
More informationSystem of Linear Equations
Math 20F Linear Algebra Lecture 2 1 System of Linear Equations Slide 1 Definition 1 Fix a set of numbers a ij, b i, where i = 1,, m and j = 1,, n A system of m linear equations in n variables x j, is given
More informationNational Taiwan University
National Taiwan University Meshless Methods for Scientific Computing (Advisor: C.S. Chen, D.L. oung) Final Project Department: Mechanical Engineering Student: Kai-Nung Cheng SID: D9956 Date: Jan. 8 The
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 informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationSolutions to practice questions for the final
Math A UC Davis, Winter Prof. Dan Romik Solutions to practice questions for the final. You are given the linear system of equations x + 4x + x 3 + x 4 = 8 x + x + x 3 = 5 x x + x 3 x 4 = x + x + x 4 =
More informationExamples of solutions to the examination in Computational Physics (FYTN03), October 24, 2016
Examples of solutions to the examination in Computational Physics (FYTN3), October 4, 6 Allowed calculation aids: TEFYMA or Physics Handbook. The examination consists of eight problems each worth five
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationTime-dependent variational forms
Time-dependent variational forms Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Oct 30, 2015 PRELIMINARY VERSION
More informationLecture 10: Fourier Sine Series
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without eplicit written permission from the copyright owner. ecture : Fourier Sine Series
More informationCheat Sheet for MATH461
Cheat Sheet for MATH46 Here is the stuff you really need to remember for the exams Linear systems Ax = b Problem: We consider a linear system of m equations for n unknowns x,,x n : For a given matrix A
More informationA Fast Fourier transform based direct solver for the Helmholtz problem. AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017
A Fast Fourier transform based direct solver for the Helmholtz problem Jari Toivanen Monika Wolfmayr AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017 Content 1 Model problem 2 Discretization 3 Fast solver
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 informationImplicit Scheme for the Heat Equation
Implicit Scheme for the Heat Equation Implicit scheme for the one-dimensional heat equation Once again we consider the one-dimensional heat equation where we seek a u(x, t) satisfying u t = νu xx + f(x,
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationUsing pseudo-parabolic and fractional equations for option pricing in jump diffusion models
Using pseudo-parabolic and fractional equations for option pricing in jump diffusion models Andrey Itkin 1 1 HAP Capital LLC and NYU Poly joined work with Peter Carr) 6-th World Bachelier Congress, Toronto,
More informationMATH SOLUTIONS TO PRACTICE MIDTERM LECTURE 1, SUMMER Given vector spaces V and W, V W is the vector space given by
MATH 110 - SOLUTIONS TO PRACTICE MIDTERM LECTURE 1, SUMMER 2009 GSI: SANTIAGO CAÑEZ 1. Given vector spaces V and W, V W is the vector space given by V W = {(v, w) v V and w W }, with addition and scalar
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
More informationWith this expanded version of what we mean by a solution to an equation we can solve equations that previously had no solution.
M 74 An introduction to Complex Numbers. 1 Solving equations Throughout the calculus sequence we have limited our discussion to real valued solutions to equations. We know the equation x 1 = 0 has distinct
More informationPartial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction
Page of 8 Partial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction FILE:Chap 3 Partial Differential Equations-V6. Original: May 7, 05 Revised: Dec 9, 06, Feb 0, 07,
More informationLecture 19: Heat conduction with distributed sources/sinks
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 ecture 19: Heat conduction
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationAssignment on iterative solution methods and preconditioning
Division of Scientific Computing, Department of Information Technology, Uppsala University Numerical Linear Algebra October-November, 2018 Assignment on iterative solution methods and preconditioning 1.
More information