LINEAR SECOND-ORDER EQUATIONS

Similar documents
Chapter 3 Second Order Linear Equations

Strauss PDEs 2e: Section Exercise 1 Page 1 of 6

Classification of Second order PDEs

LECTURE NOTES IN PARTIAL DIFFERENTIAL EQUATIONS. Fourth Edition, February by Tadeusz STYŠ. University of Botswana

Department of mathematics MA201 Mathematics III

Extra Problems and Examples

Math 342 Partial Differential Equations «Viktor Grigoryan

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

First order Partial Differential equations

Partial Differential Equations

Math 124A October 11, 2011

INTRODUCTION TO PDEs

Numerical Methods for PDEs

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u

Scientific Computing: An Introductory Survey

2. Method of Characteristics

MATH443 PARTIAL DIFFERENTIAL EQUATIONS Second Midterm Exam-Solutions. December 6, 2017, Wednesday 10:40-12:30, SA-Z02

(1) u (t) = f(t, u(t)), 0 t a.

Introduction of Partial Differential Equations and Boundary Value Problems

PARTIAL DIFFERENTIAL EQUATIONS MA 3132 LECTURE NOTES

MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20

Homework 3 solutions Math 136 Gyu Eun Lee 2016 April 15. R = b a

First order wave equations. Transport equation is conservation law with J = cu, u t + cu x = 0, < x <.

Math 220A - Fall 2002 Homework 5 Solutions

Hyperbolic PDEs. Chapter 6

Waves and characteristics: Overview 5-1

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES

M.Sc. in Meteorology. Numerical Weather Prediction

An Introduction to Numerical Methods for Differential Equations. Janet Peterson

Math 342 Partial Differential Equations «Viktor Grigoryan

The double layer potential

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

New Symmetries and Particular Solutions for 2D Black-Scholes Model

First Order Partial Differential Equations: a simple approach for beginners

Introduction to Partial Differential Equations

Numerical Integration of Linear and Nonlinear Wave Equations

FOURIER METHODS AND DISTRIBUTIONS: SOLUTIONS

Lecture Notes on Partial Dierential Equations (PDE)/ MaSc 221+MaSc 225

Tutorial 2. Introduction to numerical schemes

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

Lecture No 1 Introduction to Diffusion equations The heat equat

ANALYTICAL METHOD FOR SOLVING THE ONE-DIMENSIONAL WAVE EQUATION WITH MOVING BOUNDARY

Summer 2017 MATH Solution to Exercise 5

u xx + u yy = 0. (5.1)

Chapter 3. Second Order Linear PDEs

Numerical Solutions to Partial Differential Equations

Grid Generation and Applications

Partial Differential Equations

Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps

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

MATH 220: Problem Set 3 Solutions

Week 2: First Order Semi-Linear PDEs

UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH

Lesson 41. where the linear operator x xy y. to enforce the one to one correspondence of this transformation, we assume

MATH 220: MIDTERM OCTOBER 29, 2015

Propagation of discontinuities in solutions of First Order Partial Differential Equations

Lie Theory of Differential Equations and Computer Algebra

1 Quasi-Linear Equations and the Cauchy Kowalewski Theorem

Math512 PDE Homework 2

Partial Differential Equations, Winter 2015

Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche

Sobolev Spaces. Chapter 10

MAT 372 K.T.D.D. Final Sınavın Çözümleri N. Course. Question 1 (Canonical Forms). Consider the second order partial differential equation

Generalized evolutionary equations and their invariant solutions

Chapter 6. Partial Differential Equations. 6.1 Introduction

Grid Generation and Applications

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

Optimal Transportation. Nonlinear Partial Differential Equations

On a class of second-order PDEs admitting partner symmetries

A non-local problem with integral conditions for hyperbolic equations

Symmetry Reductions of (2+1) dimensional Equal Width. Wave Equation

Mathematical Methods - Lecture 9

5CCM211A & 6CCM211B Partial Differential Equations and Complex Variables

A review of stability and dynamical behaviors of differential equations:

HUYGENS PRINCIPLE FOR THE DIRICHLET BOUNDARY VALUE PROBLEM FOR THE WAVE EQUATION

Solving linear and nonlinear partial differential equations by the method of characteristics

PDE (Math 4163) Spring 2016

ECONOMETRIC AND GEOMETRIC ANALYSIS OF COBB-DOUGLAS AND CES PRODUCTION FUNCTIONS

Final: Solutions Math 118A, Fall 2013

ISSN: [Engida* et al., 6(4): April, 2017] Impact Factor: 4.116

Green s Functions and Distributions

Applied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.

MATH 173: PRACTICE MIDTERM SOLUTIONS

Module 2: First-Order Partial Differential Equations

Salmon: Lectures on partial differential equations

Numerical Methods in Fluid Dynamics MPO 662

Classification of cubics up to affine transformations

0.3.4 Burgers Equation and Nonlinear Wave

MATH 205C: STATIONARY PHASE LEMMA

The second-order 1D wave equation

The linear equations can be classificed into the following cases, from easier to more difficult: 1. Linear: u y. u x

SELF-SIMILAR SOLUTIONS FOR THE 2-D BURGERS SYSTEM IN INFINITE SUBSONIC CHANNELS

Oblique derivative problems for elliptic and parabolic equations, Lecture II

Math 4263 Homework Set 1

i=1 α i. Given an m-times continuously

LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Symmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota

SINC PACK, and Separation of Variables

Transcription:

LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y) in Ω. We will use the 1-1 transformation ξ = ξ(x, y), η = η(x, y) to transform the equation to one of three canonical forms determined by the sign of the discriminant B 2 AC. With w(ξ, η) = u(x(ξ, η), y(ξ, η)), the transformed equation takes the form aw ξξ + 2bw ξη + cw ηη + dw ξ + ew η + fw + g = 0, in which a = Aξ 2 x + 2Bξ x ξ y + Cξ 2 y, c = Aη 2 x + 2Bη x η y + Cη 2 y, and b = Aξ x η x + B(ξ x η y + ξ y η x ) + Cξ y η y. An important connection between the coefficients is b 2 ac = (B 2 AC)J 2 where J = ξ x η y ξ y η x ( = 0 in Ω) is the Jacobian of the 1-1 transformation. The sign of the discriminant is thus an invariant of the transformation. 1

The pde is classified according to the sign of the discriminant as following: (1) hyperbolic if B 2 AC > 0, canonical form u ξη +... = 0, (2) parabolic if B 2 AC = 0, canonical form u ξξ +... = 0, (3) elliptic if B 2 AC < 0, canonical form u ξξ + u ηη +... = 0, where... stands for terms containing at most first-order partials. The hyperbolic type To turn the equation into the canonical hyperbolic form, it is necessary to choose ξ and η that make the coefficients a and c both zero. A and C are not simultaneously zero or the equation is already in canonical form. Suppose A is not zero, then a is zero if one of the following equations is satisfied: ξ x + m ± ξ y = 0 roots are -m ± where m ± = (B ± B 2 AC)/A are two different real-value functions. The same equations, applied to η, can also make c = 0. One can pick either one of these equations for ξ and the another one for η. Therefore, there are two characteristic equations dy dx = m ± 2

that describe two families of characteristics. In terms of ξ and η, the canonical hyperbolic equation u ξη +... = 0 can be obtained. An alternate canonical form With ξ = ξ + η, η = ξ η, the canonical form of the hyperbolic type can be transformed to an alternate form: u ξ ξ u η η +... = 0. The simplest hyperbolic equation The equation u ξη = 0 has the general solution u = f(ξ) + g(η) where f, g are arbitrary differentiable functions. Solve u xx + 2u xy 3u yy = 0. u = f(y 3x) + g(y + x) Equation containing an inhomogeneous term Consider u ξη = h(ξ, η). The general solution is u = ( hdη ) dξ + f(ξ) + g(η) where f, g are arbitrary differentiable functions. 3

The parabolic type Suppose the discriminant B 2 AC = 0 in Ω, there is only one characteristic equation η x + B A η y = 0 (or dy/dx = B/A) which makes c = 0. As B 2 AC = 0 B/A = C/B Bη x + Cη y Aη x + Bη y = 0, b = (Aη x + Bη y )ξ x + (Bη x + Cη y )ξ y = 0. Thus, the equation can be written as u ξξ +... = 0. Note that as long as the Jacobian of the transformation is non-zero, ξ can be quite arbitrarily chosen. We normally adopt the simple choice ξ = x. Solve u xx + 4u xy + 4u yy = 0 η = y 2x, ξ = x, u = xf(y 2x) + g(y 2x) The elliptic type As B 2 AC < 0, there are no characteristics. While a general approach to reduce the equation is in principle possible, it can be cumbersome. We ll proceed by considering only cases with constant coefficients. 4

Equation with constant coefficients In such a case, conversion to canonical form can follow a simple procedure. Suppose that A 0 so that the equation can be divided throughout by A. (1) hyperbolic case Completing square gives ( x + B A y) 2 u ( A y) 2 u +... = 0 where = B 2 AC is the discriminant. By setting ξ = x + B A y and η = A y, the pde can be converted to the canonical form: Since u ξξ u ηη +... = 0. ξ = x ξ x + y ξ y and η = x η x + y η y, comparison of coefficients yields x y = 1, ξ ξ = B A, x y = 0, and η η = A. The transformation equations can be chosen as x = ξ and y = B A ξ + A η. 5

(2) parabolic case Rewrite the equation in the form: u ξξ ( A )2 u ηη +... = 0. As = 0, there is no constraint on x η and y η. Setting them to 0 and 1 respectively produce the transformation x = ξ, y = B ξ + η. The resultant equation is A u ξξ +... = 0. (3) elliptic case Use the transformation: x = ξ, y = B A ξ + η to get A u ξξ + u ηη +... = 0. Find a transformation to transform the equation to its canonical form. u xx + 2u xy + 3u yy + 4u = 0 ξ = x, η = (y x)/ 2 The second-order Cauchy problem Let Γ be a curve on which the values of u and its derivative u n in the direction normal to Γ are given. It is possible to find a unique solution if Γ is non-characteristic. The 6

Cauchy-Kowaleski theorem states that a unique solution exists in a neighborhood of a point of Γ if the coefficients of the pde, the functions specifying the solution and its normal derivative along Γ are analytic (possess power series expansions). Consider the following pde with the Cauchy data defined on the x-axis. u xx + u yy = 0 u(x, 0) = φ(x), u y (x, 0) = ψ(x) Use the pde and the Cauchy data to calculate the first few terms of the Taylor expansion of u about (x, 0): u(x, y) = u(x, 0) + u y (x, 0)y + 1 2! u yy(x, 0)y 2 + 1 3! u yyy(x, 0)y 3 + 1 4! u yyyy(x, 0)y 4 +... u(x,y) = φ + ψy 1 2! φ y 2 1 3! ψ y 3 + 1 4! φ(4) y 4 +... The Cauchy-Kowalewski theorem is a local existence theorem. It only asserts that a solution exists in a neighborhood of the point, not in all space. A peculiar feature of the theorem is that the type of the differential equation (wether it is hyperbolic, parabolic, or elliptic) is irrelevant. Problem Set 2 7