Finite difference methods
|
|
- Marcus Eugene Simon
- 5 years ago
- Views:
Transcription
1 Finite difference methods Siltanen/Railo/Kaarnioja Spring 8 Applications of matrix computations Applications of matrix computations Finite difference methods Spring 8 /
2 Introduction Finite difference methods (FDM) are numerical methods for solving (partial) differential equations, where (partial) derivatives are approximated by finite differences. Example For < h, we have the approximations (i) f f (x + h) f (x) (x), h (ii) f f (x) f (x h) (x), h (iii) f f (x + h) f (x h) (x), and h (iv) f f (x + h) f (x) + f (x h) (x) h. Applications of matrix computations Finite difference methods Spring 8 /
3 Let a = x < x < < x n = b be a uniform partition of the interval [a, b], where x k = a + hk, k {,,..., n}, and h = (b a)/n. The formulae on the previous page can be used to derive discrete approximations for (i) [f (x ), f (x ),..., f (x n )] T = Af, (ii) [f (x ), f (x ),..., f (x n )] T = Bf, (iii) [f (x ), f (x ),..., f (x n )] T = Cf, and (iv) [f (x ), f (x ),..., f (x n )] T = Df using the point values f = [f (x ), f (x ),..., f (x n )] T. Proof. Homework! :) Note: By using a suitable quadrature rule (e.g., the left-, mid-, or right-point rule or the trapezoidal rule), one can obtain a discretized approximation for the antiderivative as well! Applications of matrix computations Finite difference methods Spring 8 3 /
4 Finite difference methods for the equation y (t) = f (t, y(t)) Applications of matrix computations Finite difference methods Spring 8 /
5 Let s consider the IVP { y (t) = f (t, y(t)), t, y() = y. To ensure (local) unique solvability, it is assumed that f is continuous and that f (t, ) is Lipschitz-continuous [Picard Lindelöf theorem]. How do we get rid of the derivative? By integration: t y(t) = y + f (τ, y(τ)) dτ. Bad news: We have just turned our differential equation into an integral equation, which we still need to solve for y(t). Good news: We know how to do numerical integration! Applications of matrix computations Finite difference methods Spring 8 /
6 We are interested in solving the IVP for t [, T ], T >. Let = t < t < < t n = T be an equispaced partition of the interval [, T ], t k = hk, k {,,..., n}, and h = T /n. Let us denote y(t k ) = y k. We can recast the integral equation on the previous page as y k+ = y k + t k+ t k f (t, y(t)) dt for k =,,..., n. Let s try different quadratures for approximating the integral above! Remark. We omit the detailed error analysis of the resulting formulae and focus on some heuristics of the solution methods. Applications of matrix computations Finite difference methods Spring 8 6 /
7 One-point Riemann sum approximations over [t k, t k+ ] Left-point rule: y k+ y k + hf (t k, y k ). (Euler s method) Right-point rule: y k+ y k + hf (t k+, y k+ ). (Implicit Euler s method) Midpoint rule: y k+ y k + hf (t k + h, y k + h f (t k, y k )). (Runge Kutta method) Applications of matrix computations Finite difference methods Spring 8 7 /
8 Two-point trapezoidal rule over [t k, t k+ ] y k+ y k + h (f (t k, y k ) + f (t k+, y k+ )) y k + h f (t k, y k ) + h f (t k + h, y k + hf (t k, y k )). Applications of matrix computations Finite difference methods Spring 8 8 /
9 Two-point trapezoidal rule over [t k, t k+ ] y k+ y k + h (f (t k, y k ) + f (t k+, y k+ )) y k + h f (t k, y k ) + h f (t k + h, y k + hf (t k, y k )). Applications of matrix computations Finite difference methods Spring 8 8 /
10 Two-point trapezoidal rule over [t k, t k+ ] k = hf (t k, y k ), y k+ y k + k + h f (t k + h, y k + k ). Applications of matrix computations Finite difference methods Spring 8 8 /
11 Two-point trapezoidal rule over [t k, t k+ ] k = hf (t k, y k ), y k+ y k + k + h f (t k + h, y k + k ). Applications of matrix computations Finite difference methods Spring 8 8 /
12 Two-point trapezoidal rule over [t k, t k+ ] k = hf (t k, y k ), k = hf (t k + h, y k + k ), y k+ y k + k + k. This is the second order Runge Kutta method (RK). Applications of matrix computations Finite difference methods Spring 8 8 /
13 There is a similar relationship between the famous fourth order Runge Kutta method (RK) and Simpson s rule: k = hf (t k, y k ), k = hf (t k + h, y k + k ), k 3 = hf (t k + h, y k + k ), k = hf (t k+, y k + k 3 ), y k+ y k + 6 (k + k + k 3 + k ). The RK rule can be derived by first applying Simpson s rule and then balancing the midstep (here y k+/ = y(t k+/ ) and t k+/ = t k + h/) f (t k+/, y k+/ ) = f (t k+/, y k+/ ) + f (t k+/, y k+/ ) by using two different methods: f (t k+/, y k+/ ) f (t k + h, y k + k ), f (t k+/, y k+/ ) f (t k + h, y k + k ). (left-point rule) (right-point rule) Applications of matrix computations Finite difference methods Spring 8 9 /
14 In summary There exists a wide range of literature on FD methods for first order IVPs and the theory extends nicely to first order systems of ODEs as well. The examples on pg. 7 highlight a couple of canonical cases: (A) Explicit methods, (B) Implicit methods. While simple and intuitive, explicit methods are never unconditionally stable. For example, when working with [stiff equations], certain implicit methods are able to capture the asymptotics of the solution regardless of the step size. Another class of FD methods is represented by Runge Kutta methods, where the FD solution is adjusted based on point evaluations at additional intermediate points between the discretization mesh. Applications of matrix computations Finite difference methods Spring 8 /
15 For an introduction into the approximation and stability theory of first order FD methods, I highly recommend the monograph [Griffiths and Higham]! Applications of matrix computations Finite difference methods Spring 8 /
16 Finite difference methods for Poisson s equation u = f Applications of matrix computations Finite difference methods Spring 8 /
17 Let denote the Laplacian. u(x, y) = u(x, y) + u(x, y) x y In the following, we discuss the solution of Poisson s equation u = f in a bounded domain Ω R with the Dirichlet boundary condition u = g on Ω. where f and g are given functions. Applications of matrix computations Finite difference methods Spring 8 3 /
18 We need to discretize Poisson s equation in order to find a suitable matrix approximation Au = m, which we can then solve numerically. We can obtain an approximate formula for u near (x, y) by developing Taylor expansions of u(x, y) with respect to x and y variables, respectively. Assuming sufficient smoothness, we obtain u(x + h, y) = u(x, y) + h x u(x, y) + h x u(x, y) + 6 h3 3 x u(x, y) + O(h ), u(x h, y) = u(x, y) h x u(x, y) + h x u(x, y) 6 h3 3 x u(x, y) + O(h ). The sum of these is u(x + h, y) + u(x h, y) = u(x, y) + h x u(x, y) + O(h ) Applications of matrix computations Finite difference methods Spring 8 /
19 Developing the Taylor expansion w.r.t. the x variable yielded u(x + h, y) + u(x h, y) = u(x, y) + h x u(x, y) + O(h ). () Proceeding analogously with the other component: u(x, y + h) = u(x, y) + h y u(x, y) + h y u(x, y) + 6 h3 3 y u(x, y) + O(h ), u(x, y h) = u(x, y) h y u(x, y) + h y u(x, y) 6 h3 3 y u(x, y) + O(h ) u(x, y + h) + u(x, y h) = u(x, y) + h y u(x, y) + O(h ). () By summing () and () together, we obtain u(x + h, y) + u(x h, y) + u(x, y + h) + u(x, y h) = u(x, y) + h u(x, y) + O(h ) u(x, y) = u(x + h, y) + u(x h, y) + u(x, y + h) + u(x, y h) u(x, y) h + O(h ). Applications of matrix computations Finite difference methods Spring 8 /
20 Five-point stencil (quincunx): u (x,y+h) u (x-h,y) u (x,y) u (x+h,y) u (x,y-h) u(x, y) u(x + h, y) + u(x h, y) + u(x, y + h) + u(x, y h) u(x, y) h Applications of matrix computations Finite difference methods Spring 8 6 /
21 Recovering missing pixels During the lecture, we discussed the example presented in Poisson FD v.pdf (see course page). Figure: Image by Samuli Siltanen. Applications of matrix computations Finite difference methods Spring 8 7 /
22 For a given (polygonal) computational domain, we could attempt to create a quincunx mesh where, for each individual mesh element (x, y), the neighbouring coordinates (x ± h, y) and (x, y ± h) are either Other mesh elements, Boundary points Applications of matrix computations Finite difference methods Spring 8 8 /
23 unknowns u(x i) boundary values g i Example Let s solve the Dirichlet problem { u = in Ω, u Ω = g on Ω, where the computational domain Ω is the L-shaped domain on the left, the discretization mesh is denoted by black dots, and the boundary values of g are given in red at the discretization nodes lying on the boundary. Applications of matrix computations Finite difference methods Spring 8 9 /
24 Let s give a numbering to the unknowns! unknowns u(xi) boundary values gi 6 Applications of matrix computations Finite difference methods Spring 8 /
25 u(x ) u(x ) + u(x 3 ) u(x ) h [ ] u(x ) h u(x ) = h u(x 3 ) 6 Applications of matrix computations Finite difference methods Spring 8 /
26 u(x ) u(x ) + u(x ) u(x ) h [ ] u(x ) u(x ) h u(x 3 ) = [ ] h 6 u(x ) 6 Applications of matrix computations Finite difference methods Spring 8 /
27 u(x 3 ) u(x )+u(x )+u(x ) u(x 3 ) h u(x ) u(x ) h u(x 3 ) u(x ) = u(x ) h 6 Applications of matrix computations Finite difference methods Spring 8 /
28 m := [, 6,,,,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
29 m := [, 6,,, 3,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
30 m := [, 6,,, 3,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
31 m := [, 6,,, 3,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
32 m := [, 6,,, 3,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
33 m := [, 6,,, 3,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
34 =: A m := [, 6,,, 3,,,,,,,,,,, ] T u(x ) h A. = h m u(x 6 ) Applications of matrix computations Finite difference methods Spring 8 /
35 =: A m := [, 6,,, 3,,,,,, 6,,,,, ] T u(x ) h A. = h m u(x 6 ) Applications of matrix computations Finite difference methods Spring 8 /
36 =: A m := [, 6,,, 3,,,,,, 6,,,,, ] T u(x ) h A. = h m u(x 6 ) Applications of matrix computations Finite difference methods Spring 8 /
37 m := [, 6,,, 3,,,,,, 6,,,,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
38 m := [, 6,,, 3,,,,,, 6,,, 3,, ] T u(x ) h A. = h m u(x 6 ) =: A Applications of matrix computations Finite difference methods Spring 8 /
39 m := =: A [, 6,,, 3,,,,,, 6,,, 3,, ] T u(x ) h A. = h m u(x 6 ) Applications of matrix computations Finite difference methods Spring 8 /
40 m := =: A [, 6,,, 3,,,,,, 6,,, 3,, ] T u(x ) h A. = h m u(x 6 ) Applications of matrix computations Finite difference methods Spring 8 /
41 m = = A [, 6,,, 3,,,,,, 6,,, 3,, ] T u(x ) h A. = h m u(x 6 ) Applications of matrix computations Finite difference methods Spring 8 /
42 Numerical solution: >> A\m should give [u(x ),..., u(x 6 )] T = [, 3, 3,,,,,,,, 3,,,, 3, ] T. Compare this with the analytical solution u(x, y) = x y. (They are exactly the same!) Applications of matrix computations Finite difference methods Spring 8 /
43 Numerical experiment Let s return to the isospectral drums discussed earlier during the course Figure: Isospectral drum shapes in D [Gordon, Webb, and Wolpert]. Two drums with clamped boundaries give the same sound if they have the same set of (Dirichlet) eigenvalues λ satisfying u = λu in D, u D =, D R bounded domain. Applications of matrix computations Finite difference methods Spring 8 3 /
44 The quincunx pattern has extremely important applications in statistics as well! Applications of matrix computations Finite difference methods Spring 8 /
45 Bibliography C. Gordon, D. Webb, and S. Wolpert. Isospectral plane domains and surfaces via Riemannian orbifolds. Inventiones Mathematicae ():, 99. DOI:.7/BF33 D. F. Griffiths and D. J. Higham. Numerical Methods for Ordinary Differential Equations. Springer,. Applications of matrix computations Finite difference methods Spring 8 /
Mathematics 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 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 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 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 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 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 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 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 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 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 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 informationLet be a function. We say, is a plane curve given by the. Let a curve be given by function where is differentiable with continuous.
Module 8 : Applications of Integration - II Lecture 22 : Arc Length of a Plane Curve [Section 221] Objectives In this section you will learn the following : How to find the length of a plane curve 221
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 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 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 informationQuasilinear First-Order PDEs
MODULE 2: FIRST-ORDER PARTIAL DIFFERENTIAL EQUATIONS 16 Lecture 3 Quasilinear First-Order PDEs A first order quasilinear PDE is of the form a(x, y, z) + b(x, y, z) x y = c(x, y, z). (1) Such equations
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 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 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 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 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 informationChapter 13. Boundary Value Problems for Partial Differential Equations* Linz 2002/ page
Chapter 13 Boundary Value Problems for Partial Differential Equations* E lliptic equations constitute the third category of partial differential equations. As a prototype, we take the Poisson equation
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 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 informationIndex. C m (Ω), 141 L 2 (Ω) space, 143 p-th order, 17
Bibliography [1] J. Adams, P. Swarztrauber, and R. Sweet. Fishpack: Efficient Fortran subprograms for the solution of separable elliptic partial differential equations. http://www.netlib.org/fishpack/.
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 informationCOMPUTER VISION > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE
COMPUTER VISION 2017-2018 > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE OUTLINE Optical flow Lucas-Kanade Horn-Schunck Applications of optical flow Optical flow tracking Histograms of oriented flow Assignment
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 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 informationPROPERTIES OF NATURAL ELEMENT COORDINATES ON ANY POLYHEDRON
PROPRTIS OF NATURAL LMNT COORDINATS ON ANY POLYHDRON P. Milbradt and T. Fröbel Institute of Computer Science in Civil ngineering, Univercity of Hanover, 3067, Hanover, Germany; PH (+49) 5-76-5757; FAX
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 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 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 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 informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
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 informationSurface Parameterization
Surface Parameterization A Tutorial and Survey Michael Floater and Kai Hormann Presented by Afra Zomorodian CS 468 10/19/5 1 Problem 1-1 mapping from domain to surface Original application: Texture mapping
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 informationGeometric Modeling Assignment 3: Discrete Differential Quantities
Geometric Modeling Assignment : Discrete Differential Quantities Acknowledgements: Julian Panetta, Olga Diamanti Assignment (Optional) Topic: Discrete Differential Quantities with libigl Vertex Normals,
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 informationShape Modeling and Geometry Processing
252-0538-00L, Spring 2018 Shape Modeling and Geometry Processing Discrete Differential Geometry Differential Geometry Motivation Formalize geometric properties of shapes Roi Poranne # 2 Differential Geometry
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
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 informationMETU Mechanical Engineering Department ME 582 Finite Element Analysis in Thermofluids Spring 2018 (Dr. C. Sert) Handout 12 COMSOL 1 Tutorial 3
METU Mechanical Engineering Department ME 582 Finite Element Analysis in Thermofluids Spring 2018 (Dr. C. Sert) Handout 12 COMSOL 1 Tutorial 3 In this third COMSOL tutorial we ll solve Example 6 of Handout
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationAnnouncements. Edges. Last Lecture. Gradients: Numerical Derivatives f(x) Edge Detection, Lines. Intro Computer Vision. CSE 152 Lecture 10
Announcements Assignment 2 due Tuesday, May 4. Edge Detection, Lines Midterm: Thursday, May 6. Introduction to Computer Vision CSE 152 Lecture 10 Edges Last Lecture 1. Object boundaries 2. Surface normal
More informationUsing Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs).
Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs). D.I. Lanlege, Ph.D. * ; U.M. Garba, B.Sc.; and A. Aluebho, B.Sc. Department
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 informationDigital Geometry Processing Parameterization I
Problem Definition Given a surface (mesh) S in R 3 and a domain find a bective F: S Typical Domains Cutting to a Disk disk = genus zero + boundary sphere = closed genus zero Creates artificial boundary
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 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 informationCOMPUTING AND DATA ANALYSIS WITH EXCEL. Numerical integration techniques
COMPUTING AND DATA ANALYSIS WITH EXCEL Numerical integration techniques Outline 1 Quadrature in one dimension Mid-point method Trapezium method Simpson s methods Uniform random number generation in Excel,
More informationIntensity Transformations and Spatial Filtering
77 Chapter 3 Intensity Transformations and Spatial Filtering Spatial domain refers to the image plane itself, and image processing methods in this category are based on direct manipulation of pixels in
More informationChapter 6 Visualization Techniques for Vector Fields
Chapter 6 Visualization Techniques for Vector Fields 6.1 Introduction 6.2 Vector Glyphs 6.3 Particle Advection 6.4 Streamlines 6.5 Line Integral Convolution 6.6 Vector Topology 6.7 References 2006 Burkhard
More informationAn ε-uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems
An -Uniform Initial Value Technique For Convection-Diffusion Singularly Perturbed Problems RKBawa and Vinod Kumar Abstract In this paper, we have proposed an -uniform initial value technique for singularly
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 informationOn a nested refinement of anisotropic tetrahedral grids under Hessian metrics
On a nested refinement of anisotropic tetrahedral grids under Hessian metrics Shangyou Zhang Abstract Anisotropic grids, having drastically different grid sizes in different directions, are efficient and
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 informationLecture 6: Chain rule, Mean Value Theorem, Tangent Plane
Lecture 6: Chain rule, Mean Value Theorem, Tangent Plane Rafikul Alam Department of Mathematics IIT Guwahati Chain rule Theorem-A: Let x : R R n be differentiable at t 0 and f : R n R be differentiable
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationAP Calculus BC Course Description
AP Calculus BC Course Description COURSE OUTLINE: The following topics define the AP Calculus BC course as it is taught over three trimesters, each consisting of twelve week grading periods. Limits and
More information本讲内容 5.1 概述 5.2 FEM 5.3 FDM 5.4 DEM. 5.5 Case study
5 地下工程数值计算方法 本讲内容 5.1 概述 5.2 FEM 5.3 FDM 5.4 DEM 5.5 Case study Agenda 4.1 Introduction 4.2 Finite Element Method (FEM) 4.3 Finite Difference Method (FDM) 4.4 Boundary Element Method (BEM) 4.5 Discrete
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 informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationNumerical Methods for Differential Equations Contents Review of numerical integration methods Rectangular Rule Trapezoidal Rule Simpson s Rule How to
Numerical Methods for Differential Equations Contents Review of numerical integration methods Rectangular Rule Trapezoidal Rule Simpson s Rule How to make a connect-the-dots graphic Numerical Methods for
More informationThursday 1/8 1 * syllabus, Introduction * preview reading assgt., chapter 1 (modeling) * HW review chapters 1, 2, & 3
Topics and Syllabus Class # Text Reading I. NUMERICAL ANALYSIS CHAPRA AND CANALE A. INTRODUCTION AND MATLAB REVIEW :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::Week
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationProblem #3 Daily Lessons and Assessments for AP* Calculus AB, A Complete Course Page Mark Sparks 2012
Problem # Daily Lessons and Assessments for AP* Calculus AB, A Complete Course Page 490 Mark Sparks 01 Finding Anti-derivatives of Polynomial-Type Functions If you had to explain to someone how to find
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 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 information4.2 and 4.6 filled in notes.notebook. December 08, Integration. Copyright Cengage Learning. All rights reserved.
4 Integration Copyright Cengage Learning. All rights reserved. 1 4.2 Area Copyright Cengage Learning. All rights reserved. 2 Objectives Use sigma notation to write and evaluate a sum. Understand the concept
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 informationTreewidth and graph minors
Treewidth and graph minors Lectures 9 and 10, December 29, 2011, January 5, 2012 We shall touch upon the theory of Graph Minors by Robertson and Seymour. This theory gives a very general condition under
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 informationCloth Simulation. Tanja Munz. Master of Science Computer Animation and Visual Effects. CGI Techniques Report
Cloth Simulation CGI Techniques Report Tanja Munz Master of Science Computer Animation and Visual Effects 21st November, 2014 Abstract Cloth simulation is a wide and popular area of research. First papers
More informationMath 485, Graph Theory: Homework #3
Math 485, Graph Theory: Homework #3 Stephen G Simpson Due Monday, October 26, 2009 The assignment consists of Exercises 2129, 2135, 2137, 2218, 238, 2310, 2313, 2314, 2315 in the West textbook, plus the
More informationTrees. 3. (Minimally Connected) G is connected and deleting any of its edges gives rise to a disconnected graph.
Trees 1 Introduction Trees are very special kind of (undirected) graphs. Formally speaking, a tree is a connected graph that is acyclic. 1 This definition has some drawbacks: given a graph it is not trivial
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationComputer Vision I. Announcement. Corners. Edges. Numerical Derivatives f(x) Edge and Corner Detection. CSE252A Lecture 11
Announcement Edge and Corner Detection Slides are posted HW due Friday CSE5A Lecture 11 Edges Corners Edge is Where Change Occurs: 1-D Change is measured by derivative in 1D Numerical Derivatives f(x)
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
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 informationMaximizing an interpolating quadratic
Week 11: Monday, Apr 9 Maximizing an interpolating quadratic Suppose that a function f is evaluated on a reasonably fine, uniform mesh {x i } n i=0 with spacing h = x i+1 x i. How can we find any local
More informationVon Neumann Analysis for Higher Order Methods
1. Introduction Von Neumann Analysis for Higher Order Methods Von Neumann analysis is a widely used method to study how an initial wave is propagated with certain numerical schemes for a linear wave equation
More informationPTE 519 Lecture Note Finite Difference Approximation (Model)
PTE 519 Lecture Note 3 3.0 Finite Difference Approximation (Model) In this section of the lecture material, the focus is to define the terminology and to summarize the basic facts. The basic idea of any
More informationNumerical Methods for Solving Optimal Control Problems
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Masters Theses Graduate School 5-2015 Numerical Methods for Solving Optimal Control Problems Garrett Robert Rose University
More information= f (a, b) + (hf x + kf y ) (a,b) +
Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals
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 information7.8. Approximate Integration
7.8. Approimate Integration Recall that when we originally defined the definite integral of a function interval, that we did so with the limit of a Riemann sum: on an where is the number of intervals that
More informationFall CSCI 420: Computer Graphics. 4.2 Splines. Hao Li.
Fall 2014 CSCI 420: Computer Graphics 4.2 Splines Hao Li http://cs420.hao-li.com 1 Roller coaster Next programming assignment involves creating a 3D roller coaster animation We must model the 3D curve
More informationNodal Basis Functions for Serendipity Finite Elements
Nodal Basis Functions for Serendipity Finite Elements Andrew Gillette Department of Mathematics University of Arizona joint work with Michael Floater (University of Oslo) Andrew Gillette - U. Arizona Nodal
More informationIn what follows, we will focus on Voronoi diagrams in Euclidean space. Later, we will generalize to other distance spaces.
Voronoi Diagrams 4 A city builds a set of post offices, and now needs to determine which houses will be served by which office. It would be wasteful for a postman to go out of their way to make a delivery
More informationValidated Solution of Initial Value Problem for Ordinary Differential Equations based on Explicit and Implicit Runge-Kutta Schemes
Validated Solution of Initial Value Problem for Ordinary Differential Equations based on Explicit and Implicit Runge-Kutta Schemes Julien Alexandre Dit Sandretto, Alexandre Chapoutot To cite this version:
More informationSemester Final Report
CSUMS SemesterFinalReport InLaTex AnnKimball 5/20/2009 ThisreportisageneralsummaryoftheaccumulationofknowledgethatIhavegatheredthroughoutthis semester. I was able to get a birds eye view of many different
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 information1. Two double lectures about deformable contours. 4. The transparencies define the exam requirements. 1. Matlab demonstration
Practical information INF 5300 Deformable contours, I An introduction 1. Two double lectures about deformable contours. 2. The lectures are based on articles, references will be given during the course.
More information