Introduction to PDEs and Numerical Methods Tutorial 11. 2D elliptic equations
|
|
- Vincent Simmons
- 6 years ago
- Views:
Transcription
1 Platzalter für Bild, Bild auf Titelfolie inter das Logo einsetzen Introduction to PDEs and Numerical Metods Tutorial 11. 2D elliptic equations Dr. Noemi Friedman,
2 Overview Introduction (classification of PDEs, examples of important PDEs) 1D Poisson equation u x = f(x) Analytical solution (existence, uniqueness of solution, stability) Finite difference approximation of te Poisson eqaution Eigenvalue problem, Fourier transform Elliptic equations ( aa x + b x u x + c x u x = f x ) Analytical solution: existence, and uniqueness of solution, and its bounds FD approximation Consistency, stability, convergence FD approximation of PDEs on non-equidistant grids Higer order FD scemes Wit Taylor expansion Wit Lagrange interpolation Compact scemes Dr. Noemi Friedman PDE tutorial Seite 2
3 Overview Scemes for te numerical approximation of initial value problems Explicit Euler Implicit (backward) Euler Teta-metods Runge-Kutta metods Heat eqaution wit 1d spatial dimension analytical solution (seperation of ) Spatial and time discretisation, consistency and stability analysis Transport and wave equations wit 1d spatial dimension analytical solution numerical soluion Multidimensional problems 2d elliptic equation analytical analysis numerical analysis Dr. Noemi Friedman PDE tutorial Seite 3
4 Seperation of Seperation of - basic steps: 1. Suppose: u x, y = f x g y (ansatz function) 2. Plug-in ansatz function to te original PDE and rearrange in a way tat one side is a function of only x and te oter side is a function of only y: φ x = ψ y Te equality only stands if bot sides is only a scalar (separation constant - σ) φ x = ψ y = σ 3. Solve te two separated ODEs: φ x = σ ψ y = σ Dr. Noemi Friedman PDE tutorial Seite 4
5 Analytical solution of Laplace equation wit separation of Let s solve analyticaly te following problem Laplace equation u(x, y) y u xx + u yy = (LE) Boundary conditions (BC4) u x, = (BC1) u x, = (BC2) u, y = y(y ) (BC3) u x, y = lim x x u(x, y):? Dr. Noemi Friedman PDE tutorial Seite 5
6 Analytical solution of Laplace equation wit separation of 1. Ansatz function u x, y = f x g y 2. Plug-in ansatz function to te original PDE (LE) u xx + u yy = (LE) f x g y + f x g y = (PDE1) Rearrange: f x f x = g y g y (PDE2) = σ f x σσ x = g y + σg y = (ODE1) (ODE2) φ x ψ y seeeeeeeee cccccccc Dr. Noemi Friedman PDE tutorial Seite 6
7 Analytical solution of Laplace equation wit separation of 3. SOLVE (ODE1) and (ODE2) 3a) Let s assume g y + σg y = g y = g y = C 1 y + C 2 (1) (ODE2) were C 1, C 2 : constants Define C 1, C 2 from boundary conditions: u x, = f x g = (BC1) (BC1) is satisfied if f x or if g = (BC 1) would give trivial solution (rearrangement of (PDE1) to (PDE2) would not be possible if f x (BC 1) (1): g = C 1 + C 2 = C 2 = Dr. Noemi Friedman PDE tutorial Seite 7
8 Analytical solution of Laplace equation wit separation of Define C 1 from boundary conditions u x, = f x g = (BC2) (BC2) is satisfied if f x would give trivial solution or if g = (BC 2) (BC 2) (1): g = C 1 = C 1 = doesn t lead to any nontrivial solution Dr. Noemi Friedman PDE tutorial Seite 8
9 Analytical solution of Laplace equation wit separation of 3b) Let s assume g y + σg y = (ODE2) Assuming g y = e λy (2) plug-in (2) to (ODE2) λ 2 e λλ + σe λλ = λ 2 = σ λ = ± σ Te solution is te linear combination of te two possible solutions g y = C 3 e σ y + C 4 e σ y (3) were C 3, C 4 are arbitrary constants Dr. Noemi Friedman PDE tutorial Seite 9
10 Analytical solution of Laplace equation wit separation of Let s assume Define C 3, C 4 from boundary conditions: g = (BC 1) g = (BC 2) see equations (BC 1) and (BC 2) in 3a from (3) using (BC 1) g = C 3 e σ + C 4 e σ = C 3 + C 4 = C 4 = C 3 (4) from (3) using (BC 2) g = C 3 e σ + C 4 e σ = using (4) g = C 3 e σ 1 e σ = C 3 = C 4 = eiter leads to any nontrivial solution Dr. Noemi Friedman PDE tutorial Seite 1
11 Analytical solution of Laplace equation wit separation of Let s assume λ = ±i σ complex roots Te solution: g y = C 5 e i σ y + C 6 e i σ y (5) Dr. Noemi Friedman PDE tutorial Seite 11
12 Analytical solution of Laplace equation wit separation of g y = C 5 e i σ y + C 6 e i σ y = C 5 cos ( σ y) + i sin ( σ y) + C 6 cos σ y i sin ( σ y) = C 5 + C 6 cos σ y + C 5 C 6 i sin( σ y) = A = B were A, B are new constants wit C 5 + C 6 = A C 5 C 6 i = B g y = Acos σ y + B sin( σ y) (6) Let s define A and B from boundary conditions: g = (BC 1) g = (BC 2) see equations (BC 1) and (BC 2) in 3a from (6) using (BC 1) g = Acos σ + B sin( σ ) = Acos σ = A = (7) Dr. Noemi Friedman PDE tutorial Seite 12
13 Analytical solution of Laplace equation wit separation of from equations (6) (BC 2) and (7) g = B sin( σ ) = (8) equations (8) can be satisfied if B = would leave to trivial solution or if sin( σ ) = σ k = kπ were k = 1,2,3 (9) A sequance of solution: g k y = B k sin kπ y Dr. Noemi Friedman PDE tutorial Seite 13
14 Analytical solution of Laplace equation wit separation of SOLVE (ODE1) wit f x σf x = (ODE1) Assuming f x = e κx κ 2 e κκ σe κκ = Caracteristic equation: κ 2 σ = κ = ± σ real roots Te solution: f x = C 7 e σ x + C 8 e σ x (1) were C 7, C 8 are constants Dr. Noemi Friedman PDE tutorial Seite 14
15 Analytical solution of Laplace equation wit separation of Using equation (9) and (1) f k x = D k e kπ x + E k e kπ x (11) And te ansatz: u k x, y = f k x g k y = D k e kπ x + E k e kπ x B k sin kπ y (12) Let s define B k, D k and E k from boundary conditions (BC3) and (BC4). First from (BC4): lim u x, y = lim x x D ke kπ x + E k e kπ x B k sin kπ y = lim e kπ x = x lim x D k e kπ x B k sin kπ y = only if D k = Dr. Noemi Friedman PDE tutorial Seite 15
16 Analytical solution of Laplace equation wit separation of u k x, y = B k E k e kπ x sin kπ y = A ke kπ x sin kπ y A k Te general solution will be te sum of te sequances: u x, y = k=1 u k x, y = A k e kπ x sin kπ k=1 y Let s define A k from boundary conditions (BC3): u, y = u (y) = (y 2 2 /4) (BC3) u, y = A k e kπ sin k=1 kπ y = A k sin kπ y = u (y) k=1 Accordingly te A k coefficients are te coefficients of te Fourier serious of te function u (y) Dr. Noemi Friedman PDE tutorial Seite 16
17 Analytical solution of Laplace equation wit separation of Let s define A k from boundary conditions (BC3): u, y = u y = y(y ) (BC3) u, y = A k e kπ sin k=1 kπ y = A k sin kπ y = u (y) k=1 Accordingly te A k coefficients are te coefficients of te Fourier serious of te function u (y) A k = 2 u (y) sin kπ y dy = 2 ( y(y )) sin kπ y dy u x, y = 2 u (y) sin k=1 kπ y dy e kπ x sin kπ y Excercise: Calculate analyticaly te A k coefficients Dr. Noemi Friedman PDE tutorial Seite 17
18 Analytical solution of Laplace equation wit separation of A k = 2 (y(y )) sin kπ y dy = 2 (y2 y) sin kπ y dy = 2 2 kπ (y2 y) cos kπ y + (2y ) cos kπ y dy = kπ (y2 y) cos kπ y + kπ (2y ) sin kπ y kπ 2 sin kπ y dy 2 kπ (y2 y) cos kπ y + kπ (2y ) sin kπ y + kπ 2 2 cos kπ y = = = 2 1 k 2 A k = 8 2 kπ 3 ii k ii ooo ii k ii eeee Dr. Noemi Friedman PDE tutorial Seite 18
19 Analytical solution of Laplace equation wit separation of u x, y = 8 2 2n 1 3 π 3 e n=1 2n 1 π x sin 2n 1 π y Dr. Noemi Friedman PDE tutorial Seite 19
20 2D Laplace equation - finite difference metod Let s consider te Laplace equation in 2 dimensions: If dx=dy: Dr. Noemi Friedman PDE tutorial Seite 2
21 n = 2D Laplace equation - finite difference metod u n 1,1 u n 2,1 u n 3,1 u n M 1,1 u n 1,2 u n 2,2 u n M 1,2 u n 1,N 1 u n M 1,N 1 u n = u 1,1 u 2,1 u 3,1 u 4,1 u 1,2 u 2,2 u 3,2 u 4,2 u 1,3 u 2,3 u 3,3 u 4,3 l = 1.. N y = lδy (, N) (,1) (1,1) (2,1) (3,1) (M, 1) (,) (1,) (2,) (3,) x = jδx (2,3) (1,2) (2,2) (3,2) j = 1.. M (M, N) (M, 1) wit omogenous Diriclet BC. r 4u j,l + u j 1,l + u j+1,l + u j,l 1 + u j,l+1 = f j,l Dr. Noemi Friedman PDE tutorial Seite 21
22 2D Laplace equation - finite difference metod r 4u j,l + u j 1,l + u j+1,l + u j,l 1 + u j,l+1 = f j,l wit omogenous Diriclet BC r u 1,1 u 2,1 u 3,1 u 4,1 u 1,2 u 2,2 u 3,2 u 4,2 u 1,3 u 2,3 u 3,3 u 4,3 Sparse matrix wit bandwidt: 2M-1 (ere 9) BUT te band itself is sparse, only five diagonals are nonzero Dr. Noemi Friedman PDE tutorial Seite 22
23 2D Laplace equation - finite difference metod r 4u j,l + u j 1,l + u j+1,l + u j,l 1 + u j,l+1 = f j,l B 1 u = f were B 1 = B C C B C C B B = 1 + 4r r r 1 + 4r r r 1 + 4r C = r r r Dr. Noemi Friedman PDE tutorial Seite 23
Introduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 5. Finite difference methods equilibrium equation and iterative solvers Dr. Noemi
More informationIntroduction to PDEs and Numerical Methods 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 informationIntroduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Tutorial 4. Finite difference methods stability, concsistency, convergence Dr. Noemi Friedman,
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
More informationIntroduction to PDEs and Numerical Methods Lecture 7. Solving linear systems
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 7. Solving linear systems Dr. Noemi Friedman, 09.2.205. Reminder: Instationary heat
More informationIntroduction to PDEs and Numerical Methods Lecture 1: Introduction
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 1: Introduction Dr. Noemi Friedman, 28.10.2015. Basic information on the course Course
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen PARAMetric UNCertainties, Budapest STOCHASTIC PROCESSES AND FIELDS Noémi Friedman Institut für Wissenschaftliches Rechnen, wire@tu-bs.de
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More informationFall 2014 MAT 375 Numerical Methods. Numerical Differentiation (Chapter 9)
Fall 2014 MAT 375 Numerical Metods (Capter 9) Idea: Definition of te derivative at x Obviuos approximation: f (x) = lim 0 f (x + ) f (x) f (x) f (x + ) f (x) forward-difference formula? ow good is tis
More informationChapter 8. Numerical Solution of Ordinary Differential Equations. Module No. 2. Predictor-Corrector Methods
Numerical Analysis by Dr. Anita Pal Assistant Professor Department of Matematics National Institute of Tecnology Durgapur Durgapur-7109 email: anita.buie@gmail.com 1 . Capter 8 Numerical Solution of Ordinary
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationExamination paper for TMA4130 Matematikk 4N: SOLUTION
Department of Mathematical Sciences Examination paper for TMA4 Matematikk 4N: SOLUTION Academic contact during examination: Morten Nome Phone: 9 84 97 8 Examination date: December 7 Examination time (from
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
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 informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationMethod of Separation of Variables
MODUE 5: HEAT EQUATION 11 ecture 3 Method of Separation of Variables Separation of variables is one of the oldest technique for solving initial-boundary value problems (IBVP) and applies to problems, where
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 information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationPhysically Based Modeling: Principles and Practice Implicit Methods for Differential Equations
Pysically Based Modeling: Principles and Practice Implicit Metods for Differential Equations David Baraff Robotics Institute Carnegie Mellon University Please note: Tis document is 997 by David Baraff
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationSOLUTION OF EQUATION AND EIGENVALUE PROBLEMS PART A( 2 MARKS)
CHENDU COLLEGE OF ENGINEERING AND TECHNOLOGY (Approved by AICTE New Delhi, Affiliated to Anna University Chennai. Zamin Endathur Village, Madurntakam Taluk, Kancheepuram Dist.-603311.) MA6459 - NUMERICAL
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More information2.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 informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationComputation 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 informationRunge-Kutta methods. With orders of Taylor methods yet without derivatives of f (t, y(t))
Runge-Kutta metods Wit orders of Taylor metods yet witout derivatives of f (t, y(t)) First order Taylor expansion in two variables Teorem: Suppose tat f (t, y) and all its partial derivatives are continuous
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationLecture 21. Numerical differentiation. f ( x+h) f ( x) h h
Lecture Numerical differentiation Introduction We can analytically calculate te derivative of any elementary function, so tere migt seem to be no motivation for calculating derivatives numerically. However
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationNumerical Solution to Parabolic PDE Using Implicit Finite Difference Approach
Numerical Solution to arabolic DE Using Implicit Finite Difference Approac Jon Amoa-Mensa, Francis Oene Boateng, Kwame Bonsu Department of Matematics and Statistics, Sunyani Tecnical University, Sunyani,
More informationModule 4 Linear Boundary value problems
Module 4 Linear Boundary value problems Lecture Content Hours 1 Finite Difference Metods: Diriclet type boundary 1 condition Finite Difference Metods: Mixed boundary condition 1 3 Sooting Metod 1 4 Sooting
More informationNumerical 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 informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
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 informationThird part: finite difference schemes and numerical dispersion
Tird part: finite difference scemes and numerical dispersion BCAM - Basque Center for Applied Matematics Bilbao, Basque Country, Spain BCAM and UPV/EHU courses 2011-2012: Advanced aspects in applied matematics
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationENGI Partial Differentiation Page y f x
ENGI 344 4 Partial Differentiation Page 4-0 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
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 informationIntroduction to PDEs and Numerical Methods Tutorial 10. Finite Element Analysis
Patzhater für Bid, Bid auf Titefoie hinter das Logo einsetzen Introduction to PDEs and Numerica Methods Tutoria. Finite Eement Anaysis Dr. Noemi Friedman, 3..2 FROM STRONG FORM TO WEAK FORM inhomogeneous
More informationDELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Mathematics and Computer Science
DELFT UNIVERSITY OF TECHNOLOGY Faculty of Electrical Engineering, Matematics and Computer Science. ANSWERS OF THE TEST NUMERICAL METHODS FOR DIFFERENTIAL EQUATIONS (WI3097 TU) Tuesday January 9 008, 9:00-:00
More information2.4 Eigenvalue problems
2.4 Eigenvalue problems Associated with the boundary problem (2.1) (Poisson eq.), we call λ an eigenvalue if Lu = λu (2.36) for a nonzero function u C 2 0 ((0, 1)). Recall Lu = u. Then u is called an eigenfunction.
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationNUMERICAL 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 information1. Differential Equations (ODE and PDE)
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
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
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 information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationTwo Step Hybrid Block Method with Two Generalized Off-step Points for Solving Second Ordinary Order Differential Equations Directly
Global Journal of Pure and Applied Matematics. ISSN 0973-768 Volume 2, Number 2 (206), pp. 59-535 Researc India Publications ttp://www.ripublication.com/gjpam.tm Two Step Hybrid Block Metod wit Two Generalized
More informationIn this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation.
Eigen Function Expansion and Applications. In this worksheet we will use the eigenfunction expansion to solve nonhomogeneous equation. a/ The theory. b/ Example: Solving the Euler equation in two ways.
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationDarcy s law in 3-D. K * xx K * yy K * zz
PART 7 Equations of flow Darcy s law in 3-D Specific discarge (vector) is calculated by multiplying te ydraulic conductivity (second-order tensor) by te ydraulic gradient (vector). We obtain a general
More informationIntroduction to Multigrid Method
Introduction to Multigrid Metod Presented by: Bogojeska Jasmina /08/005 JASS, 005, St. Petersburg 1 Te ultimate upsot of MLAT Te amount of computational work sould be proportional to te amount of real
More informationProblem 1. Possible Solution. Let f be the 2π-periodic functions defined by f(x) = cos ( )
Problem Let f be the π-periodic functions defined by f() = cos ( ) hen [ π, π]. Make a draing of the function f for the interval [ 3π, 3π], and compute the Fourier series of f. Use the result to compute
More informationPartial Differential Equations
Next: Using Matlab Up: Numerical Analysis for Chemical Previous: Ordinary Differential Equations Subsections Finite Difference: Elliptic Equations The Laplace Equations Solution Techniques Boundary Conditions
More informationLecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems
Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now
More informationBasic 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 informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
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 informationExercise 19 - OLD EXAM, FDTD
Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms
More informationResearch Article Cubic Spline Iterative Method for Poisson s Equation in Cylindrical Polar Coordinates
International Scolarly Researc Network ISRN Matematical Pysics Volume 202, Article ID 2456, pages doi:0.5402/202/2456 Researc Article Cubic Spline Iterative Metod for Poisson s Equation in Cylindrical
More informationLecture 19: Heat conduction with distributed sources/sinks
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 ecture 19: Heat conduction
More informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationAMS 147 Computational Methods and Applications Lecture 09 Copyright by Hongyun Wang, UCSC. Exact value. Effect of round-off error.
Lecture 09 Copyrigt by Hongyun Wang, UCSC Recap: Te total error in numerical differentiation fl( f ( x + fl( f ( x E T ( = f ( x Numerical result from a computer Exact value = e + f x+ Discretization error
More informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationQuaternion Dynamics, Part 1 Functions, Derivatives, and Integrals. Gary D. Simpson. rev 01 Aug 08, 2016.
Quaternion Dynamics, Part 1 Functions, Derivatives, and Integrals Gary D. Simpson gsim1887@aol.com rev 1 Aug 8, 216 Summary Definitions are presented for "quaternion functions" of a quaternion. Polynomial
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More information1 Separation of Variables
Jim ambers ENERGY 281 Spring Quarter 27-8 ecture 2 Notes 1 Separation of Variables In the previous lecture, we learned how to derive a PDE that describes fluid flow. Now, we will learn a number of analytical
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 informationOrdinary Differential Equations
Ordinary Differential Equations (MA102 Mathematics II) Shyamashree Upadhyay IIT Guwahati Shyamashree Upadhyay ( IIT Guwahati ) Ordinary Differential Equations 1 / 1 Books Shyamashree Upadhyay ( IIT Guwahati
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
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 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 informationM. Croci (Oxford), M. B. Giles (Oxford), P. E. Farrell (Oxford), M. E. Rognes (Simula) MCQMC July 4, 2018
Efficient wite noise sampling and coupling for multilevel Monte Carlo M. Croci (Oxford), M. B. Giles (Oxford), P. E. Farrell (Oxford), M. E. Rognes (Simula) MCQMC2018 - July 4, 2018 Overview Introduction
More information