1. Differential Equations (ODE and PDE)
|
|
- Bridget Maxwell
- 6 years ago
- Views:
Transcription
1 1. Differential Equations (ODE and PDE) 1.1. Ordinary Differential Equations (ODE) So far we have dealt with Ordinary Differential Equations (ODE): involve derivatives with respect to only one variable if the sought after function depends on more than one variable, the extra variables were treated as constants! examples u (t) = λu(t) (radioactive decay) ( ) p (t) = r 1 p(t) p(t) (population model, logistic growth) k (ṗ(t) ) ( ) f(p(t), q(t)) p(t) = (predator-prey) q(t) g(p(t), q(t)) q(t) usually, initial conditions were required Page 1 of 19
2 1.2. (PDE) involve (partial) derivatives with respect to more than one variable, and the sought after function depends on more than one variable boundary conditions and/or initial conditions required (depending on type of PDE) Examples: 1D heat equation: 2 u(x, t) = u(x, t) t x2 also written as u t (x, t) = u xx (x, t) 2D Laplace equation: 2 2 u(x, y) + u(x, y) = x2 y2 also written as u xx (x, y) + u yy (x, y) = Page 2 of 19
3 2. choosing PDE instead of ODE in modelling is often a result of adding spatial dimensions: in population modelling: ODE-modelling assumes homogeneous distribution of the population in space switch to PDE if assumption is no longer applicable heat conduction: localising hot spots requires spatial resolution (i.e. PDE) time is often treated as a special dimension, we distinguish stationary problems: no time-dependence unsteady/instationary problems: time-dependence (perhaps, but not necessarily, with a stationary limit) Page 3 of 19
4 Heat Conduction models propagation of heat within a given object examples: a heated wire (1D) a metal plate, heated or cooled at its boundaries (2D) water cooling in a reactor (3D problem) air conditioning: where to place ventilations for heating or cooling boiling water in a pot function of interest: temperature T T (x, t), T (x, y), T (x, y, t) or T (x, y, z, t) or... heat propagation depends on the material s heat conductivity. Page 4 of 19
5 3.2. Derivation of the Heat Equation basic idea for modelling: analyse conservation of energy given an arbitrary part D of our computational domain the change of heat energy (per time) is a result of transfer of heat energy across D s surface, heat sources and drains in D (external influences) ρct dv = q dv + k T n ds t D D density ρ, specific heat c, and heat conductivity k are material parameters heat sources and drains are modelled in term q n is outer normal vector on surface element dv/ds according to Gauß theorem: k T n ds = k T dv D D D Page 5 of 19
6 leads to integral equation for arbitrary domain D: ρct t q k T dv = D hence, the integrand has to be identically : T t = κ T + q ρc, κ := k ρc κ > is called the thermal diffusion coefficient (since the Laplace operator models a (heat) diffusion process) Heat Equation: For vanishing external influence q =, we finally get the heat equation: T t = κ T alternate notation T t = κ ( ) 2 T x + 2 T 2 y + 2 T 2 z 2 Page 6 of 19
7 3.3. Boundary and Initial Conditions the heat equation needs boundary or initial-boundary conditions to provide a unique solution Dirichlet boundary conditions: fix T on (part of) the boundary Neumann boundary conditions: T (x, y, z) = ϕ(x, y, z) fix T s normal derivative on (part of) the boundary: T (x, y, z) = ϕ(x, y, z) n pure Dirichlet conditions and mixtures of Dirichlet and Neumann conditions lead to a unique solution. pure Neumann b.c. determines solution up to a constant (if T is a solution, T + c solves the PDE, as well) Page 7 of 19
8 Physical interpretation of boundary conditions: Dirichlet conditions: the temperature T is prescribed itself along (part of) the boundary (some defined heating or cooling) Neumann conditions: the temperature flux across (part of) the boundary is prescribed special case: homogeneous Neumann conditions: complete isolation, no orthogonal transport of heat into or out of the domain For time-dependent heat equation(s): initial conditions required, e.g. for t = Analytical solution: can be given for simple configurations (separation of variables/fourier s method, see exercises) heat equation is one of few examples of PDE, where general statements concerning existence and uniqueness of solutions are possible. Page 8 of 19
9 4. Poisson s Equation in One Dimension We remember the general heat conduction equation, T t = κ T + q ρc. Assuming a steady state problem (T t = ) on the unit interval (i.e. in the 1D case), we get Poisson s equation: or, using an alternate notation: T xx (x) = f(x), x (, 1) u xx (x) = f(x) or even u (x) = f(x) f(x) combines all material parameters and the source term q. Further, we assume homogeneous Dirichlet boundary conditions: u() = u(1) = Page 9 of 19
10 5. Recalling a fundamental theorem of calculus, we may state that x for a suitable constant c 1, and similar u(x) = c 1 + u (ξ) dξ (1) x u (x) = c 2 + u (ζ) dζ (2) Now, suppose that u(x) satisfies the PDE ( u = f(x)), then: u (x) = c 2 f(ζ) dζ ξ Inserting this result into equation (1), we get u(x) = c 1 + c 2 x x ξ f(ζ) dζ dξ Page 1 of 19
11 To simplify the equation, we define ξ F (ξ) := f(ζ) dζ and integrate by parts to get x ξ x x f(ζ) dζ dξ = F (ξ)dξ = [ξf (ξ)] x ξf (ξ)dξ x x x = xf (x) ξf(ξ)dξ = xf(ξ)dξ ξf(ξ)dξ = x (x ξ)f(ξ)dξ Thus, we get as a preliminary result that x u(x) = c 1 + c 2 x (x ξ)f(ξ)dξ (3) Page 11 of 19
12 5.1. Considering Boundary Conditions So far c 1 and c 2 are arbitrary constants. equation (3) solves the PDE in (, 1) the boundary conditions u() = u(1) = have not been taken into account yet. We therefore compute c 1 and c 2 to satisfy u() = u(1) =. u() = leads to u() = c 1 = u(1) = leads to 1 u(x) = c 2 (1 ξ)f(ξ)dξ = or 1 c 2 = (1 ξ)f(ξ)dξ Page 12 of 19
13 5.2. Green s Function As a final result, we get 1 x u(x) = x (1 ξ)f(ξ)dξ (x ξ)f(ξ)dξ being the solution of our boundary value problem. Examples: f(x) = 1, then u(x) = x 1 (1 ξ)dξ x (x ξ)dξ = 1 x(1 x) 2 f(x) = x (see homework) f(x) = e x, then (see Maple worksheet) u(x) = 1 e x + (e 1)x Page 13 of 19
14 Green s Function: Introduce the function G(x, ξ) := { ξ(1 x) if ξ x x(1 ξ) if x ξ 1 then our solution u(x) can be simplified to 1 u(x) = G(x, ξ)f(ξ) dξ (4) Properties: G is continuous G is symmetric: G(x, y) = G(y, x) G(, ξ) = G(1, ξ) = G(x, ) = G(x, 1) = G is a piecewise linear function of x for fixed ξ (and vice versa!) G is positive: G(x, ξ) for all x, ξ [, 1] Page 14 of 19
15 5.3. A Maximum Principle For any continuous function f C([, 1]), let f = sup f(x) x [,1] (maximum norm) Maximum principle (proposition): If u is our solution given by equation (4), then u 1 8 f if f is continuous on [, 1], i.e. f C([, 1]). Proof: 1 1 u(x) G(x, ξ) f(ξ) dξ f G(x, ξ)dξ A short computation leads to Page 15 of 19 u 1 8 f
16 5.4. Uniqueness of the Solution Consider two boundary value problems u xx (x) = f 1 (x) and u xx (x) = f 2 (x), u() = u(1) =, with their respective solutions u 1 (x) and u 2 (x). Due to the linearity of the problem, u 1 (x) u 2 (x) is then a solution of the boundary value problem u xx (x) = f 1 (x) f 2 (x). If f 1 = f 2 this means that, if the solutions u 1 and u 2 were actually different, the maximum principle would require that u 1 u f 1 f 2 = Hence, u 1 and u 2 must be identical. That means, the solution to our 1D boundary value problem is unique! Page 16 of 19
17 6. A Finite Difference Approximation Similar to numerical Methods for ODE, we will approximate the solution on a grid of n equidistant grid points (meshsize h := 1 n ): x i := ih, i =,..., n; compute approximations u i u(x i ) for i =,..., n; replace the second derivative by an appropriate difference term u (x i ) = u i+1 2u i + u i 1 h 2 where the error term E h (x) is O(h 2 ); + E h (x) Thus, we get the following system equations: u i+1 2u i + u i 1 h 2 = f(x i ) for i = 1,..., n 1; (5) with boundary conditions u = u() = and u n = u(1) =. Page 17 of 19
18 6.1. Setting up a Linear System of Equations The linear system of equations (5) can also be written in matrixvector-notation: Au = b where u = (u 1,..., u n 1 ) T is the vector of unknowns; b = (b 1,..., b n 1 ) T, with component b i given by b i := h 2 f(x i ) for i = 1,..., n 1; the system matrix A is given by A = Page 18 of 19
19 6.2. Properties of the Linear System of Equations the matrix A is diagonal dominant, which means that for matrix elements a ij the equation a ii i j a ij holds for every i = 1,..., n 1 as a consequence, all eigenvalues of A are positive, and A is non-singular. therefore, the linear system of equations has a unique solution Page 19 of 19
Part I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationScientific Computing I
Scientific Computing I Session 11: Basics of Partial Differential Equations Tobias Weinzierl Winter 2010/2011 Session 11: Basics of Partial Differential Equations, Winter 2010/2011 1 Motivation Instead
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 6. Solve the completely inhomogeneous diffusion problem on the half-line
Strauss PDEs 2e: Section 3.3 - Exercise 2 Page of 6 Exercise 2 Solve the completely inhomogeneous diffusion problem on the half-line v t kv xx = f(x, t) for < x
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationMathematical Methods - Lecture 9
Mathematical Methods - Lecture 9 Yuliya Tarabalka Inria Sophia-Antipolis Méditerranée, Titane team, http://www-sop.inria.fr/members/yuliya.tarabalka/ Tel.: +33 (0)4 92 38 77 09 email: yuliya.tarabalka@inria.fr
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationApplications of the Maximum Principle
Jim Lambers MAT 606 Spring Semester 2015-16 Lecture 26 Notes These notes correspond to Sections 7.4-7.6 in the text. Applications of the Maximum Principle The maximum principle for Laplace s equation is
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationMath 5587 Lecture 2. Jeff Calder. August 31, Initial/boundary conditions and well-posedness
Math 5587 Lecture 2 Jeff Calder August 31, 2016 1 Initial/boundary conditions and well-posedness 1.1 ODE vs PDE Recall that the general solutions of ODEs involve a number of arbitrary constants. Example
More informationIntroduction 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 informationClass Meeting # 2: The Diffusion (aka Heat) Equation
MATH 8.52 COURSE NOTES - CLASS MEETING # 2 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 2: The Diffusion (aka Heat) Equation The heat equation for a function u(, x (.0.). Introduction
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationCHAPTER 1 Introduction to Differential Equations 1 CHAPTER 2 First-Order Equations 29
Contents PREFACE xiii CHAPTER 1 Introduction to Differential Equations 1 1.1 Introduction to Differential Equations: Vocabulary... 2 Exercises 1.1 10 1.2 A Graphical Approach to Solutions: Slope Fields
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationMATH-UA 263 Partial Differential Equations Recitation Summary
MATH-UA 263 Partial Differential Equations Recitation Summary Yuanxun (Bill) Bao Office Hour: Wednesday 2-4pm, WWH 1003 Email: yxb201@nyu.edu 1 February 2, 2018 Topics: verifying solution to a PDE, dispersion
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationLECTURE # 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 informationINTRODUCTION 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 informationIntroduction to Partial Differential Equation - I. Quick overview
Introduction to Partial Differential Equation - I. Quick overview To help explain the correspondence between a PDE and a real world phenomenon, we will use t to denote time and (x, y, z) to denote the
More informationReading: 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 informationChapter 12 Partial Differential Equations
Chapter 12 Partial Differential Equations Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 12.1 Basic Concepts of PDEs Partial Differential Equation A
More informationMATH3203 Lecture 1 Mathematical Modelling and ODEs
MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1
More informationClassification 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 informationNumerical 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 informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationNumerical Analysis and Simulation of Partial Differential Equations
Numerical Analysis and Simulation of Partial Differential Equations Roland Pulch Institut für Mathematik und Informatik Ernst-Moritz-Arndt-Universität Greifswald Contents: 1. Examples and Classification
More informationPreliminary 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 informationAn Introduction to Numerical Methods for Differential Equations. Janet Peterson
An Introduction to Numerical Methods for Differential Equations Janet Peterson Fall 2015 2 Chapter 1 Introduction Differential equations arise in many disciplines such as engineering, mathematics, sciences
More informationCS 542G: The Poisson Problem, Finite Differences
CS 542G: The Poisson Problem, Finite Differences Robert Bridson November 10, 2008 1 The Poisson Problem At the end last time, we noticed that the gravitational potential has a zero Laplacian except at
More informationElectrodynamics PHY712. Lecture 3 Electrostatic potentials and fields. Reference: Chap. 1 in J. D. Jackson s textbook.
Electrodynamics PHY712 Lecture 3 Electrostatic potentials and fields Reference: Chap. 1 in J. D. Jackson s textbook. 1. Poisson and Laplace Equations 2. Green s Theorem 3. One-dimensional examples 1 Poisson
More information17 Source Problems for Heat and Wave IB- VPs
17 Source Problems for Heat and Wave IB- VPs We have mostly dealt with homogeneous equations, homogeneous b.c.s in this course so far. Recall that if we have non-homogeneous b.c.s, then we want to first
More informationA Computational Approach to Study a Logistic Equation
Communications in MathematicalAnalysis Volume 1, Number 2, pp. 75 84, 2006 ISSN 0973-3841 2006 Research India Publications A Computational Approach to Study a Logistic Equation G. A. Afrouzi and S. Khademloo
More informationStaple or bind all pages together. DO NOT dog ear pages as a method to bind.
Math 3337 Homework Instructions: Staple or bind all pages together. DO NOT dog ear pages as a method to bind. Hand-drawn sketches should be neat, clear, of reasonable size, with axis and tick marks appropriately
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationQuestion 9: PDEs Given the function f(x, y), consider the problem: = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1. x 2 u. u(x, 0) = u(x, 1) = 0 for 0 x 1
Question 9: PDEs Given the function f(x, y), consider the problem: 2 u x 2 u = f(x, y) 2 y2 for 0 < x < 1 and 0 < x < 1 u(x, 0) = u(x, 1) = 0 for 0 x 1 u(0, y) = u(1, y) = 0 for 0 y 1. a. Discuss how you
More informationModule 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 informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationMATH 220 solution to homework 5
MATH 220 solution to homework 5 Problem. (i Define E(t = k(t + p(t = then E (t = 2 = 2 = 2 u t u tt + u x u xt dx u 2 t + u 2 x dx, u t u xx + u x u xt dx x [u tu x ] dx. Because f and g are compactly
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationMath Partial Differential Equations 1
Math 9 - Partial Differential Equations Homework 5 and Answers. The one-dimensional shallow water equations are h t + (hv) x, v t + ( v + h) x, or equivalently for classical solutions, h t + (hv) x, (hv)
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
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 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 informationMath 251 December 14, 2005 Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt
Math 251 December 14, 2005 Final Exam Name Section There are 10 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning of each question
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationPartial Differential Equations Summary
Partial Differential Equations Summary 1. The heat equation Many physical processes are governed by partial differential equations. temperature of a rod. In this chapter, we will examine exactly that.
More informationQ ( q(m, t 0 ) n) S t.
THE HEAT EQUATION The main equations that we will be dealing with are the heat equation, the wave equation, and the potential equation. We use simple physical principles to show how these equations are
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 15 Heat with a source So far we considered homogeneous wave and heat equations and the associated initial value problems on the whole line, as
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationBOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES
1 BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES 1.1 Separable Partial Differential Equations 1. Classical PDEs and Boundary-Value Problems 1.3 Heat Equation 1.4 Wave Equation 1.5 Laplace s Equation
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationSome Aspects of Solutions of Partial Differential Equations
Some Aspects of Solutions of Partial Differential Equations K. Sakthivel Department of Mathematics Indian Institute of Space Science & Technology(IIST) Trivandrum - 695 547, Kerala Sakthivel@iist.ac.in
More informationSeparation of variables
Separation of variables Idea: Transform a PDE of 2 variables into a pair of ODEs Example : Find the general solution of u x u y = 0 Step. Assume that u(x,y) = G(x)H(y), i.e., u can be written as the product
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationDiffusion - The Heat Equation
Chapter 6 Diffusion - The Heat Equation 6.1 Goal Understand how to model a simple diffusion process and apply it to derive the heat equation in one dimension. We begin with the fundamental conservation
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations MATH 365 Ordinary Differential Equations J Robert Buchanan Department of Mathematics Fall 2018 Objectives Many physical problems involve a number of separate
More informationMath 124A October 11, 2011
Math 14A October 11, 11 Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This corresponds to a string of infinite length. Although
More informationMATH 173: PRACTICE MIDTERM SOLUTIONS
MATH 73: PACTICE MIDTEM SOLUTIONS This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve all of them. Write your solutions to problems and in blue book #, and your
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationFOURIER TRANSFORMS. 1. Fourier series 1.1. The trigonometric system. The sequence of functions
FOURIER TRANSFORMS. Fourier series.. The trigonometric system. The sequence of functions, cos x, sin x,..., cos nx, sin nx,... is called the trigonometric system. These functions have period π. The trigonometric
More informationCalculus C (ordinary differential equations)
Calculus C (ordinary differential equations) Lesson 9: Matrix exponential of a symmetric matrix Coefficient matrices with a full set of eigenvectors Solving linear ODE s by power series Solutions to linear
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationModeling using conservation laws. Let u(x, t) = density (heat, momentum, probability,...) so that. u dx = amount in region R Ω. R
Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...) so that u dx = amount in region R Ω. R Modeling using conservation laws Let u(x, t) = density (heat, momentum, probability,...)
More information18 Green s function for the Poisson equation
8 Green s function for the Poisson equation Now we have some experience working with Green s functions in dimension, therefore, we are ready to see how Green s functions can be obtained in dimensions 2
More informationM.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 informationLecture 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 informationBoundary value problems
1 Introduction Boundary value problems Lecture 5 We have found that the electric potential is a solution of the partial differential equation; 2 V = ρ/ǫ 0 The above is Poisson s equation where ρ is the
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
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 informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationGetting 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 informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationDiffusion of a density in a static fluid
Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More information6 Non-homogeneous Heat Problems
6 Non-homogeneous Heat Problems Up to this point all the problems we have considered for the heat or wave equation we what we call homogeneous problems. This means that for an interval < x < l the problems
More information