Partial Differential Equations II

Size: px
Start display at page:

Download "Partial Differential Equations II"

Transcription

1 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

2 Almost Done! Homework 7: 12/2 (two days late!) Homework 8: 12/9 (optional) Section: 12/6 (final review) Final exam: 12/12, 12:15pm (Gates B03) Go to office hours! CS 205A: Mathematical Methods Partial Differential Equations II 2 / 28

3 Course Reviews On Axess! Additional comments: CS 205A: Mathematical Methods Partial Differential Equations II 3 / 28

4 Request for Help CS 205A notes your help! Textbook Review text Write reference implementations Solidify your CS205A knowledge CS 205A: Mathematical Methods Partial Differential Equations II 4 / 28

5 Final Exam Cumulative Similar format to midterm Two sheets of notes CS 205A: Mathematical Methods Partial Differential Equations II 5 / 28

6 This Week Couple relationships between derivatives. Pressure gradient determining fluid flow Image operators using x and y derivatives Partial Differential Equations (PDE) CS 205A: Mathematical Methods Partial Differential Equations II 6 / 28

7 Boundary Value Problems Dirichlet conditions: Value of f( x) on Ω Neumann conditions: Derivatives of f( x) on Ω Mixed or Robin conditions: Combination CS 205A: Mathematical Methods Partial Differential Equations II 7 / 28

8 Second-Order Model Equation ij a ij f x i x j + i b i f x i + cf = 0 ( A + b + c)f = 0 CS 205A: Mathematical Methods Partial Differential Equations II 8 / 28

9 Classification of Second-Order PDE ( A + b + c)f = 0 If A is positive or negative definite, system is elliptic. If A is positive or negative semidefinite, the system is parabolic. If A has only one eigenvalue of different sign from the rest, the system is hyperbolic. If A satisfies none of the criteria, the system is ultrahyperbolic. CS 205A: Mathematical Methods Partial Differential Equations II 9 / 28

10 Derivative Operator Matrix h 2 w = L 1 y Dirichlet CS 205A: Mathematical Methods Partial Differential Equations II 10 / 28

11 What About First Derivative? Potential for asymmetry at boundary Centered differences: Fencepost problem Possible resolution: Imitate leapfrog CS 205A: Mathematical Methods Partial Differential Equations II 11 / 28

12 Fencepost Problem CS 205A: Mathematical Methods Partial Differential Equations II 12 / 28

13 Big Idea Derivatives : Functions :: Matrices : Vectors CS 205A: Mathematical Methods Partial Differential Equations II 13 / 28

14 Elliptic PDE Lf = g L y = b CS 205A: Mathematical Methods Partial Differential Equations II 14 / 28

15 Elliptic PDE Lf = g L y = b Example: Laplace s equation on a line CS 205A: Mathematical Methods Partial Differential Equations II 14 / 28

16 Common Theme Elliptic PDE Positive definite matrix CS 205A: Mathematical Methods Partial Differential Equations II 15 / 28

17 Common Theme Elliptic PDE Positive definite matrix L = D D, D = CS 205A: Mathematical Methods Partial Differential Equations II 15 / 28

18 Common Theme Elliptic PDE Positive definite matrix L = D D, D = Review: Name two ways to solve. CS 205A: Mathematical Methods Partial Differential Equations II 15 / 28

19 Time Dependence Choice: 1. Treat t separate from x ( semidiscrete ) 2. Treat all variables democratically ( fully discrete ) CS 205A: Mathematical Methods Partial Differential Equations II 16 / 28

20 Semidiscrete Heat Equation f t = f xx CS 205A: Mathematical Methods Partial Differential Equations II 17 / 28

21 Semidiscrete Heat Equation f t = f xx f t = Lf CS 205A: Mathematical Methods Partial Differential Equations II 17 / 28

22 Semidiscrete Heat Equation f t = f xx f t = Lf Stability for elliptic spatial operator (parabolic PDE) CS 205A: Mathematical Methods Partial Differential Equations II 17 / 28

23 Semidiscrete Time Stepping Left with a multivariable ODE problem! Forward/backward Euler, RK, and friends CS 205A: Mathematical Methods Partial Differential Equations II 18 / 28

24 Semidiscrete Time Stepping Left with a multivariable ODE problem! Forward/backward Euler, RK, and friends Implicit vs. explicit (vs. symplectic) CS 205A: Mathematical Methods Partial Differential Equations II 18 / 28

25 Semidiscrete Time Stepping Left with a multivariable ODE problem! Forward/backward Euler, RK, and friends Implicit vs. explicit (vs. symplectic) Alternative: Eigenvector methods (low-frequency approximation) CS 205A: Mathematical Methods Partial Differential Equations II 18 / 28

26 Fully Discrete PDE Discretize x and t simultaneously Can create larger linear algebra problems Philosophical point: What is fully discrete? CS 205A: Mathematical Methods Partial Differential Equations II 19 / 28

27 Reminders Review Numerical PDEs Gradient Domain Inpainting Gradient Domain Inpainting CS 205A: Mathematical Methods Partial Differential Equations II 20 / 28 Fluids

28 Gradient Domain Pipeline for image I(x, y): 1. Compute gradient: v(x, y) = I(x, y) 2. Edit: v v 3. Reconstruct: g =? v CS 205A: Mathematical Methods Partial Differential Equations II 21 / 28

29 Gradient Domain Reconstruction min g Ω g v 2 2 da CS 205A: Mathematical Methods Partial Differential Equations II 22 / 28

30 Gradient Domain Reconstruction min g Ω g v 2 2 da 2 g = v Elliptic! CS 205A: Mathematical Methods Partial Differential Equations II 22 / 28

31 Incompressible Navier-Stokes ( v ρ t ) + v v = p + µ 2 v + f t [0, ): Time v(t) : Ω R 3 : Velocity ρ(t) : Ω R: Density p(t) : Ω R: Pressure f(t) : Ω R 3 : External forces (e.g. gravity) CS 205A: Mathematical Methods Partial Differential Equations II 23 / 28

32 Lagrangian vs. Eulerian Lagrangian: Track parcels of fluid Eulerian: Fluid flows past a point in space CS 205A: Mathematical Methods Partial Differential Equations II 24 / 28

33 Marker-and-Cell (MAC) Grid CS 205A: Mathematical Methods Partial Differential Equations II 25 / 28

34 Splitting for Incompressible Flow u = 0 (divergence-free) ρ t + u ρ = 0 (density advection) u t + u u + p ρ = g (velocity advection) CS 205A: Mathematical Methods Partial Differential Equations II 26 / 28

35 Steps for Flow (on board) 1. Adjust t 2. Advect velocity 3. Apply forces 4. Solve for pressure: p ρ divergence-free projection 5. Advect density = u; CS 205A: Mathematical Methods Partial Differential Equations II 27 / 28

36 Semilagrangian Advection ecmwf.int/newsevents/training/rcourse_notes/numerical_methods/numerical_methods/numerical_methods6.html Next CS 205A: Mathematical Methods Partial Differential Equations II 28 / 28

Partial Differential Equations

Partial Differential Equations Chapter 14 Partial Differential Equations Our intuition for ordinary differential equations generally stems from the time evolution of physical systems. Equations like Newton s second law determining the

More information

Ordinary Differential Equations II

Ordinary Differential Equations II Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No

More information

Ordinary Differential Equations II

Ordinary 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 information

Math background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids

Math background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure

More information

PDE Solvers for Fluid Flow

PDE Solvers for Fluid Flow PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic

More information

Course Information Course Overview Study Skills Background Material. Introduction. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics

Course Information Course Overview Study Skills Background Material. Introduction. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Introduction CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James CS 205A: Mathematical Methods Introduction 1 / 16 Instructor Prof. Doug James Office: Gates 363 Telephone: (650)

More information

Ordinary Differential Equations I

Ordinary Differential Equations I Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 35

More information

Partial Differential Equations

Partial 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 information

Module 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 1: Finite Difference Method

Module 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 1: Finite Difference Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture1/1_1.htm 1 of 1 6/19/2012 4:29 PM The Lecture deals with: Classification of Partial Differential Equations Boundary and Initial Conditions Finite Differences

More information

7 Hyperbolic Differential Equations

7 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 information

Finite 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 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO

More information

Basic Aspects of Discretization

Basic Aspects of Discretization Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)

More information

Ordinary Differential Equations I

Ordinary Differential Equations I Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 27 Theme of Last Three

More information

Scientific Computing: An Introductory Survey

Scientific Computing: An Introductory Survey Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.

More information

Ordinary Differential Equations I

Ordinary 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 information

Partitioned Methods for Multifield Problems

Partitioned Methods for Multifield Problems C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible

More information

INTRODUCTION TO PDEs

INTRODUCTION TO PDEs INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial

More information

Final Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University

Final Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University Final Examination CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University The exam runs for 3 hours. The exam contains eight problems. You must complete the first

More information

An Overview of Fluid Animation. Christopher Batty March 11, 2014

An Overview of Fluid Animation. Christopher Batty March 11, 2014 An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.

More information

Linear Algebra Review

Linear Algebra Review Linear Algebra Review CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Linear Algebra Review 1 / 16 Midterm Exam Tuesday Feb

More information

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation

A Study on Numerical Solution to the Incompressible Navier-Stokes Equation A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field

More information

Fluid Animation. Christopher Batty November 17, 2011

Fluid Animation. Christopher Batty November 17, 2011 Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces

More information

Introduction to Partial Differential Equations

Introduction to Partial Differential Equations Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define

More information

PARTIAL DIFFERENTIAL EQUATIONS. MTH 5230, Fall 2007, MW 6:30 pm - 7:45 pm. George M. Skurla Hall 116

PARTIAL DIFFERENTIAL EQUATIONS. MTH 5230, Fall 2007, MW 6:30 pm - 7:45 pm. George M. Skurla Hall 116 PARTIAL DIFFERENTIAL EQUATIONS MTH 5230, Fall 2007, MW 6:30 pm - 7:45 pm George M. Skurla Hall 116 Ugur G. Abdulla Office Hours: S311, TR 2-3 pm COURSE DESCRIPTION The course presents partial diffrential

More information

Block-Structured Adaptive Mesh Refinement

Block-Structured Adaptive Mesh Refinement Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell

More information

Numerical Methods for Engineers and Scientists

Numerical Methods for Engineers and Scientists Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL

More information

d v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =

d v 2 v = d v d t i n where in and rot denote the inertial (absolute) and rotating frames. Equation of motion F = Governing equations of fluid dynamics under the influence of Earth rotation (Navier-Stokes Equations in rotating frame) Recap: From kinematic consideration, d v i n d t i n = d v rot d t r o t 2 v rot

More information

Differential equations, comprehensive exam topics and sample questions

Differential equations, comprehensive exam topics and sample questions Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions

More information

Webster s horn model on Bernoulli flow

Webster s horn model on Bernoulli flow Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular

More information

Index. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems

Index. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,

More information

Announce Statistical Motivation ODE Theory Spectral Embedding Properties Spectral Theorem Other. Eigenproblems I

Announce Statistical Motivation ODE Theory Spectral Embedding Properties Spectral Theorem Other. Eigenproblems I Eigenproblems I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Eigenproblems I 1 / 33 Announcements Homework 1: Due tonight

More information

Course Syllabus: Continuum Mechanics - ME 212A

Course Syllabus: Continuum Mechanics - ME 212A Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester

More information

Lecture Introduction

Lecture Introduction Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary

More information

Elliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II

Elliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods

More information

Numerical Solutions to PDE s

Numerical Solutions to PDE s Introduction Numerical Solutions to PDE s Mathematical Modelling Week 5 Kurt Bryan Let s start by recalling a simple numerical scheme for solving ODE s. Suppose we have an ODE u (t) = f(t, u(t)) for some

More information

Physics-Based Animation

Physics-Based Animation CSCI 5980/8980: Special Topics in Computer Science Physics-Based Animation 13 Fluid simulation with grids October 20, 2015 Today Presentation schedule Fluid simulation with grids Course feedback survey

More information

Numerical Solutions of Partial Differential Equations

Numerical Solutions of Partial Differential Equations Numerical Solutions of Partial Differential Equations Dr. Xiaozhou Li xiaozhouli@uestc.edu.cn School of Mathematical Sciences University of Electronic Science and Technology of China Introduction Overview

More information

Structure of the Comprehensive Examination in the ME Department. For circulation to students

Structure of the Comprehensive Examination in the ME Department. For circulation to students Structure of the Comprehensive Examination in the ME Department For circulation to students i. The qualifying exams will be held up to 3 times every year. ii. Generally, the qualifying examination will

More information

ME Computational Fluid Mechanics Lecture 5

ME 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 information

Before we look at numerical methods, it is important to understand the types of equations we will be dealing with.

Before we look at numerical methods, it is important to understand the types of equations we will be dealing with. Chapter 1. Partial Differential Equations (PDEs) Required Readings: Chapter of Tannehill et al (text book) Chapter 1 of Lapidus and Pinder (Numerical Solution of Partial Differential Equations in Science

More information

Math 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 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 information

Computational Fluid Dynamics-1(CFDI)

Computational Fluid Dynamics-1(CFDI) بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to

More information

Finite Difference Methods for

Finite 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 information

Math 330 (Section 7699 ): Fall 2015 Syllabus

Math 330 (Section 7699 ): Fall 2015 Syllabus College of Staten Island, City University of New York (CUNY) Math 330 (Section 7699 ): Fall 2015 Syllabus Instructor: Joseph Maher Applied Mathematical Analysis I Office: 1S-222 Phone: (718) 982-3623 Email:

More information

Classification of partial differential equations and their solution characteristics

Classification of partial differential equations and their solution characteristics 9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.

More information

Linear Systems and LU

Linear Systems and LU Linear Systems and LU CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Linear Systems and LU 1 / 48 Homework 1. Homework 0:

More information

Additive Manufacturing Module 8

Additive Manufacturing Module 8 Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p

More information

Announce Statistical Motivation Properties Spectral Theorem Other ODE Theory Spectral Embedding. Eigenproblems I

Announce Statistical Motivation Properties Spectral Theorem Other ODE Theory Spectral Embedding. Eigenproblems I Eigenproblems I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Eigenproblems I 1 / 33 Announcements Homework 1: Due tonight

More information

Numerical Methods for Partial Differential Equations: an Overview.

Numerical Methods for Partial Differential Equations: an Overview. Numerical Methods for Partial Differential Equations: an Overview math652_spring2009@colorstate PDEs are mathematical models of physical phenomena Heat conduction Wave motion PDEs are mathematical models

More information

Introduction to Computational Fluid Dynamics

Introduction to Computational Fluid Dynamics AML2506 Biomechanics and Flow Simulation Day Introduction to Computational Fluid Dynamics Session Speaker Dr. M. D. Deshpande M.S. Ramaiah School of Advanced Studies - Bangalore 1 Session Objectives At

More information

Chapter 5 Types of Governing Equations. Chapter 5: Governing Equations

Chapter 5 Types of Governing Equations. Chapter 5: Governing Equations Chapter 5 Types of Governing Equations Types of Governing Equations (1) Physical Classification-1 Equilibrium problems: (1) They are problems in which a solution of a given PDE is desired in a closed domain

More information

First, Second, and Third Order Finite-Volume Schemes for Diffusion

First, Second, and Third Order Finite-Volume Schemes for Diffusion First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.

More information

Time-adaptive methods for the incompressible Navier-Stokes equations

Time-adaptive methods for the incompressible Navier-Stokes equations Time-adaptive methods for the incompressible Navier-Stokes equations Joachim Rang, Thorsten Grahs, Justin Wiegmann, 29.09.2016 Contents Introduction Diagonally implicit Runge Kutta methods Results with

More information

Diffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics

Diffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods

More information

Introduction to PDEs and Numerical Methods Lecture 1: Introduction

Introduction 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 information

FBE / Advanced Topics in Engineering Mathematics. Date Fall Credits 3 credits Course Title

FBE / Advanced Topics in Engineering Mathematics. Date Fall Credits 3 credits Course Title Date Fall 2015-2016 Credits 3 credits Course Title Course Number Math 501 Advanced Topics in Engineering Mathematics Pre-requisite (s) None Co-requisite (s) None Hours 36 Out of Class 90 Work Hours Place

More information

Length Learning Objectives Learning Objectives Assessment

Length Learning Objectives Learning Objectives Assessment Universidade Federal Fluminense PGMEC Course: Advanced Computational Fluid Dynamics Coordinator: Vassilis Theofilis Academic Year: 2018, 2 nd Semester Length: 60hrs (48hrs classroom and 12hrs tutorials)

More information

Computation Fluid Dynamics

Computation Fluid Dynamics Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand

More information

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference

More information

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9

2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9 Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers

More information

Numerical Methods for Partial Differential Equations CAAM 452. Spring 2005

Numerical Methods for Partial Differential Equations CAAM 452. Spring 2005 Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks

More information

Table of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard

Table of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm

More information

SOLVING ELLIPTIC PDES

SOLVING ELLIPTIC PDES university-logo SOLVING ELLIPTIC PDES School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 POISSON S EQUATION Equation and Boundary Conditions Solving the Model Problem 3 THE LINEAR ALGEBRA PROBLEM

More information

M.Sc. in Meteorology. Numerical Weather Prediction

M.Sc. in Meteorology. Numerical Weather Prediction M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section

More information

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned

More information

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,

More information

Conservation of Mass (Eulerian Framework)

Conservation of Mass (Eulerian Framework) CME306 / CS205B Homework 1 (Theory) Conservation of Mass (Eulerian Framework) 1. In an Eulerian framework, the strong form of Conservation of Mass takes the form below. Please briefly explain the three

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Final Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Spring 2016), Stanford University

Final Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Spring 2016), Stanford University Final Examination CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Spring 2016), Stanford University The exam runs for 3 hours. The exam contains seven problems. You must complete the

More information

Index. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2

Index. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2 Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

Problem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015

Problem 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 information

Streamline calculations. Lecture note 2

Streamline calculations. Lecture note 2 Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

Numerical methods for the Navier- Stokes equations

Numerical methods for the Navier- Stokes equations Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:

More information

Introduction to PDEs and Numerical Methods Tutorial 1: Overview of essential linear algebra and analysis

Introduction 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 information

AMSC/MATH 673, CLASSICAL METHODS IN PDE, FALL Required text: Evans, Partial Differential Equations second edition

AMSC/MATH 673, CLASSICAL METHODS IN PDE, FALL Required text: Evans, Partial Differential Equations second edition AMSC/MATH 673, CLASSICAL METHODS IN PDE, FALL 2018. MWF 2:00pm - 2:50pm MTH 0407 Instructor: M. Machedon Office: MTH 3311 e-mail: mxm@math.umd.edu Required text: Evans, Partial Differential Equations second

More information

x j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2

x j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2 Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological

More information

Finite Difference Methods (FDMs) 2

Finite 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 information

arxiv: v3 [math.na] 18 Feb 2019

arxiv: v3 [math.na] 18 Feb 2019 Constraint-consistent Runge-Kutta methods for one-dimensional incompressible multiphase flow B. Sanderse a, A.E.P. Veldman b a Centrum Wiskunde & Informatica (CWI), Amsterdam, The Netherlands b Bernoulli

More information

Intro to Research Computing with Python: Partial Differential Equations

Intro to Research Computing with Python: Partial Differential Equations Intro to Research Computing with Python: Partial Differential Equations Erik Spence SciNet HPC Consortium 28 November 2013 Erik Spence (SciNet HPC Consortium) PDEs 28 November 2013 1 / 23 Today s class

More information

Engineering. Mathematics. GATE 2019 and ESE 2019 Prelims. For. Comprehensive Theory with Solved Examples

Engineering. Mathematics. GATE 2019 and ESE 2019 Prelims. For. Comprehensive Theory with Solved Examples Thoroughly Revised and Updated Engineering Mathematics For GATE 2019 and ESE 2019 Prelims Comprehensive Theory with Solved Examples Including Previous Solved Questions of GATE (2003-2018) and ESE-Prelims

More information

Engineering Mathematics

Engineering Mathematics Thoroughly Revised and Updated Engineering Mathematics For GATE 2017 and ESE 2017 Prelims Note: ESE Mains Electrical Engineering also covered Publications Publications MADE EASY Publications Corporate

More information

Announcements Monday, November 13

Announcements Monday, November 13 Announcements Monday, November 13 The third midterm is on this Friday, November 17. The exam covers 3.1, 3.2, 5.1, 5.2, 5.3, and 5.5. About half the problems will be conceptual, and the other half computational.

More information

Numerical Analysis of Differential Equations Numerical Solution of Parabolic Equations

Numerical 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 information

Preliminary Examination, Numerical Analysis, August 2016

Preliminary Examination, Numerical Analysis, August 2016 Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any

More information

Reading: P1-P20 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering)

Reading: P1-P20 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering) Chapter 1. Partial Differential Equations Reading: P1-P0 of Durran, Chapter 1 of Lapidus and Pinder (Numerical solution of Partial Differential Equations in Science and Engineering) Before even looking

More information

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

PH.D. PRELIMINARY EXAMINATION MATHEMATICS UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

More information

Physically Based Simulations (on the GPU)

Physically Based Simulations (on the GPU) Physically Based Simulations (on the GPU) (some material from slides of Mark Harris) CS535 Fall 2014 Daniel G. Aliaga Department of Computer Science Purdue University Simulating the world Floating point

More information

Modeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan

Modeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods

More information

Module 7: The Laplace Equation

Module 7: The Laplace Equation Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

Some notes about PDEs. -Bill Green Nov. 2015

Some notes about PDEs. -Bill Green Nov. 2015 Some notes about PDEs -Bill Green Nov. 2015 Partial differential equations (PDEs) are all BVPs, with the same issues about specifying boundary conditions etc. Because they are multi-dimensional, they can

More information

Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice

Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice 1 Lecture Notes, HCI, 4.1.211 Chapter 2 Gradient Descent and Implementation Solving the Euler-Lagrange Equations in Practice Bastian Goldlücke Computer Vision Group Technical University of Munich 2 Bastian

More information

McGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA

McGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper

More information

Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction

Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de

More information

NUMERICAL METHODS FOR ENGINEERING APPLICATION

NUMERICAL METHODS FOR ENGINEERING APPLICATION NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto

More information

AS MATHEMATICS HOMEWORK C1

AS MATHEMATICS HOMEWORK C1 Student Teacher AS MATHEMATICS HOMEWORK C September 05 City and Islington Sixth Form College Mathematics Department www.candimaths.uk HOMEWORK INTRODUCTION You should attempt all the questions. If you

More information

Introduction to PDEs and Numerical Methods: Exam 1

Introduction 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 information