ENGI 9420 Lecture Notes 8 - PDEs Page 8.01

Similar documents
ENGI 9420 Lecture Notes 8 - PDEs Page 8.01

Partial Differential Equations

ENGI 4430 PDEs - d Alembert Solutions Page 11.01

ENGI 9420 Engineering Analysis Solutions to Additional Exercises

MA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE

1 Separation of Variables

Lecture6. Partial Differential Equations

ENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations

Partial Differential Equations (PDEs)

Partial Differential Equations Summary

1. Partial differential equations. Chapter 12: Partial Differential Equations. Examples. 2. The one-dimensional wave equation

Chapter 3 Second Order Linear Equations

MATH 251 Final Examination August 14, 2015 FORM A. Name: Student Number: Section:

CHAPTER 10 NOTES DAVID SEAL

Chapter 12 Partial Differential Equations

MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation

Separation of Variables

Math 201 Assignment #11

Method of Separation of Variables

Math Assignment 14

A Guided Tour of the Wave Equation

Differential Equations

Separation of variables in two dimensions. Overview of method: Consider linear, homogeneous equation for u(v 1, v 2 )

Partial Differential Equations

Introduction to the Wave Equation

Wave Equation With Homogeneous Boundary Conditions

Boundary-value Problems in Rectangular Coordinates

Name: Math Homework Set # 5. March 12, 2010

5. Suggestions for the Formula Sheets

Midterm 2: Sample solutions Math 118A, Fall 2013

MA 201: Partial Differential Equations Lecture - 10

Separation of Variables. A. Three Famous PDE s

Math 2930 Worksheet Final Exam Review

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

Homework 7 Solutions

BOUNDARY-VALUE PROBLEMS IN RECTANGULAR COORDINATES

THE METHOD OF SEPARATION OF VARIABLES

Analysis III (BAUG) Assignment 3 Prof. Dr. Alessandro Sisto Due 13th October 2017

Heat Equation, Wave Equation, Properties, External Forcing

Partial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators

Practice Problems For Test 3

ENGI 9420 Lecture Notes 1 - ODEs Page 1.01

Final: Solutions Math 118A, Fall 2013

Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.

MATH 251 Final Examination August 10, 2011 FORM A. Name: Student Number: Section:

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))

Math 251 December 14, 2005 Answer Key to Final Exam. 1 18pt 2 16pt 3 12pt 4 14pt 5 12pt 6 14pt 7 14pt 8 16pt 9 20pt 10 14pt Total 150pt

Branch: Name of the Student: Unit I (Fourier Series) Fourier Series in the interval (0,2 l) Engineering Mathematics Material SUBJECT NAME

Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series

CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE

MATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:

Salmon: Lectures on partial differential equations

Boundary value problems for partial differential equations

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:

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

MATH 251 Final Examination May 3, 2017 FORM A. Name: Student Number: Section:

Introduction of Partial Differential Equations and Boundary Value Problems

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

u tt = a 2 u xx u tt = a 2 (u xx + u yy )

Solving the Heat Equation (Sect. 10.5).

Lecture Notes for Ch 10 Fourier Series and Partial Differential Equations

MATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:

Differential Equations

Partial Differential Equations for Engineering Math 312, Fall 2012

An Introduction to Numerical Methods for Differential Equations. Janet Peterson

CHAPTER 4. Introduction to the. Heat Conduction Model

d Wave Equation. Rectangular membrane.

MA Chapter 10 practice

Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche

Lecture notes: Introduction to Partial Differential Equations

Practice Problems For Test 3

Connection to Laplacian in spherical coordinates (Chapter 13)

LECTURE 19: SEPARATION OF VARIABLES, HEAT CONDUCTION IN A ROD

Advanced Eng. Mathematics

Solutions to Exercises 8.1

Wave Equation Modelling Solutions

MATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012

3150 Review Problems for Final Exam. (1) Find the Fourier series of the 2π-periodic function whose values are given on [0, 2π) by cos(x) 0 x π f(x) =

INTRODUCTION TO PDEs

Math 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 Methods for Polymer Physics Lecture 1: Series Representations of Functions

Chapter 10: Partial Differential Equations

MATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:

Higher-order ordinary differential equations

MATH-UA 263 Partial Differential Equations Recitation Summary

The One Dimensional Heat Equation

c2 2 x2. (1) t = c2 2 u, (2) 2 = 2 x x 2, (3)

Introduction to Partial Differential Equations

Last Update: April 7, 201 0

Plot of temperature u versus x and t for the heat conduction problem of. ln(80/π) = 820 sec. τ = 2500 π 2. + xu t. = 0 3. u xx. + u xt 4.

1. Solve the boundary-value problems or else show that no solutions exist. y (x) = c 1 e 2x + c 2 e 3x. (3)

STABILITY. Phase portraits and local stability

First order Partial Differential equations

Vibrating-string problem

Lecture Notes on PDEs

Instructor s Solutions Manual PARTIAL DIFFERENTIAL EQUATIONS. with FOURIER SERIES and BOUNDARY VALUE PROBLEMS. NAKHLÉ H. ASMAR University of Missouri

LEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.

FINAL EXAM, MATH 353 SUMMER I 2015

14 Separation of Variables Method

Differential equations, comprehensive exam topics and sample questions

Transcription:

ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives with respect to those variables. Most (but not all) physical models in engineering that result in partial differential equations are of at most second order and are often linear. (Some problems such as elastic stresses and bending moments of a beam can be of fourth order). In this course we shall have time to look at only a small subset of second order linear partial differential equations. Sections in this Chapter: 8.0 Major Classifications of Common PDEs 8.0 The Wave Equation Vibrating Finite String 8.03 The Wave Equation Vibrating Infinite String 8.04 d Alembert Solution 8.05 The Maximum-Minimum Principle 8.06 The Heat Equation

ENGI 940 8.0 - Classification of PDEs Page 8.0 8.0 Major Classifications of Common PDEs A general second order linear partial differential equation in two Cartesian variables can be written as u u u u u Axy (, ) + Bxy (, ) + C( xy, ) = f xyu,,,, y y y Three main types arise, based on the value of D = B 4AC (a discriminant): Hyperbolic, wherever (x, y) is such that D > 0; Parabolic, wherever (x, y) is such that D = 0; Elliptic, wherever (x, y) is such that D < 0. Among the most important partial differential equations in engineering are: u The wave equation: = c u u u or its one-dimensional special case = c [which is hyperbolic everywhere] (where u is the displacement and c is the speed of the wave); u The heat (or diffusion) equation: µρ = K u + K u a one-dimensional special case of which is u K u = [which is parabolic everywhere] µρ (where u is the temperature, μ is the specific heat of the medium, ρ is the density and K is the thermal conductivity); The potential (or aplace s) equation: a special case of which is u u y + = 0 u = 0 [which is elliptic everywhere] The complete solution of a PDE requires additional information, in the form of initial conditions (values of the dependent variable and its first partial derivatives at t = 0), boundary conditions (values of the dependent variable on the boundary of the domain) or some combination of these conditions.

ENGI 940 8.0 - Finite String Page 8.03 8.0 The Wave Equation Vibrating Finite String The wave equation is u = c u If u(x, t) is the vertical displacement of a point at location x on a vibrating string at time t, then the governing PDE is u u = c If u(x, y, t) is the vertical displacement of a point at location (x, y) on a vibrating membrane at time t, then the governing PDE is u u u = c u = c + y or, in plane polar coordinates (r, θ), (appropriate for a circular drum), u u u u = c + + r r r r θ 3 For spherically symmetric waves in, the wave equation is u c u u u = r c = + r r r r r r The latter two cases lead to solutions involving Bessel functions. Example 8.0. An elastic string of length is fixed at both ends (x = 0 and x = ). The string is displaced into the form y = f (x) and is released from rest. Find the displacement y(x, t) at all locations on the string (0 < x < ) and at all subsequent times (t > 0). The boundary value problem for the displacement function y(x, t) is: y y = c for 0 < x < and t > 0 Both ends fixed for all time: y(0, t) = y(, t) = 0 for t > 0 Initial configuration of string: y(x, 0) = f (x) for 0 < x < String released from rest: y ( x,0) = 0 for 0 x

ENGI 940 8.0 - Finite String Page 8.04 Example 8.0. (continued) Separation of Variables (or Fourier Method) Attempt a solution of the form y(x, t) = X(x) T(t) Substitute y(x, t) = X(x) T(t) into the PDE: dt d X ( X( xt ) ( t )) = c ( X( xt ) ( t) ) X = c T dt dx dt d X = c T dt X dx The left hand side of this equation is a function of t only. At any instant t it must have the same value at all values of x. Therefore the right hand side, which is a function of x only, must at any one instant have that same value at all values of x. By a similar argument, the right hand side of this equation is a function of x only. At any location x it must have the same value at all times t. Therefore the left hand side, which is a function of t only, must at any one location have that same value at all times t. Thus both sides of this differential equation must be the same absolute constant, which we shall represent for now by k. d X d X = k + kx = 0 X dx dx The general solution of this simple second order ODE is a combination of two exponential functions of x if k < 0, is a linear function of x if k = 0 and is a combination of sine and cosine functions of x if k > 0. It is not possible for both ends of the string to be fixed for all time in the first two cases (unless we admit the trivial solution y(x, t) 0, a string that never moves from its equilibrium position). Therefore k > 0. Replace k by λ (guaranteed positive for all real λ except λ = 0). We now have the pair of ODEs d X dt + λ X = 0 and + λ ct = 0 dx dt The general solutions are X x = A cos λx + B sin λx and T t = C cos λct + Dsin λct respectively, where A, B, C and D are arbitrary constants. ( ) ( ) ( ) ( ) ( ) ( )

ENGI 940 8.0 - Finite String Page 8.05 Example 8.0. (continued) Consider the boundary conditions: y 0, t = X 0 T t = 0 t 0 ( ) ( ) ( ) For a non-trivial solution, this requires X ( ) ( ) ( ) ( ) ( ) 0 = 0 A= 0. y, t = X T t = 0 t 0 X = 0 nπ Bsin = 0 n =, n ( λ ) λ ( ) We now have a solution only for a discrete set of eigenvalues λ n, with corresponding eigenfunctions Xn ( x) nπ x = sin, ( n =,, 3, ) Consider the initial condition: y = X( xt ) ( 0) = 0 x T ( 0) = 0 ( x,0) ( ) ( ) ( ) ( ) T t = Cλc sin λct + Dλc cos λct T 0 = Dλc = 0 D = 0 Therefore our complete solution for y(x, t) is now some linear combination of nπ x nπct yn xt, = Xn xtn t = Cn sin cos, n=,, 3, ( ) ( ) ( ) ( ) There is one condition remaining to be satisfied. The initial configuration of the string is: y(x, 0) = f (x) for 0 < x <. n = nπ x y x = Cn = f x (,0) sin ( ) This is precisely the Fourier sine series expansion of f (x) on [0, ]! From Fourier series theory (Chapter 7), the coefficients C n are nπ u Cn = f ( u) sin du 0 Therefore our complete solution is nπu nπ x nπct y ( x, t) = f ( u) sin du sin cos 0 n =

ENGI 940 8.0 - Finite String Page 8.06 Example 8.0. (continued) This solution is valid for any initial displacement function f (x) that is continuous with a piece-wise continuous derivative on [0, ] with f (0) = f () = 0. nπ x If the initial displacement is itself sinusoidal f ( x) = asin for some n, then the complete solution is a single term from the infinite series, nπ x nπct y( xt, ) = asin cos. Example 8.0. Suppose that the initial configuration is triangular: ( 0 ) ( ) x x y( x,0) = f ( x) = x < x Then the Fourier sine coefficients are nπ u Cn = f ( u) sin du 0 / π sin 0 / n u nπu = u du + ( u) sin du / nπu nπu nπu = cos + sin nπ nπ ( u) nπu nπu + cos sin / 0

ENGI 940 8.0 - Finite String Page 8.07 Example 8.0. (continued) nπ nπ nπ = cos sin 0 0 + ( nπ ) nπ nπ nπ + ( 0 + 0) cos sin ( nπ ) ( k ) ( ) ( n ) ( n ) 4 nπ nπ 0 even = sin But sin = ± odd π k+ sin = ( ), ( k ) Therefore sum over the odd integer values of n only (n = k ). 4 k+ Ck = ( ) k π and (( ) ) k+ ( ) ( k ) k = ( ) π ( ) 4 k x k πct y( xt, ) = sin cos π See the web page "www.engr.mun.ca/~ggeorge/940/demos/ex80.html" for an animation of this solution.

ENGI 940 8.0 - Finite String Page 8.08 Example 8.0. (continued) Some snapshots of the solution are shown here: These graphs were generated from the Fourier series, truncated after the fifth non-zero term.

ENGI 940 8.0 - Finite String Page 8.09 Example 8.0.3 An elastic string of length is fixed at both ends (x = 0 and x = ). The string is initially in its equilibrium state [y(x, 0) = 0 for all x] and is released with the initial velocity y = g( x). Find the displacement y(x, t) at all locations on the string (0 < x < ) ( x,0) and at all subsequent times (t > 0). The boundary value problem for the displacement function y(x, t) is: y y = c for 0 < x < and t > 0 Both ends fixed for all time: y(0, t) = y(, t) = 0 for t > 0 Initial configuration of string: y(x, 0) = 0 for 0 < x < y String released with initial velocity: ( x,0) ( ) = g x for 0 x As before, attempt a solution by the method of the separation of variables. Substitute y(x, t) = X(x) T(t) into the PDE: dt d X ( X( xt ) ( t )) = c ( X( xt ) ( t) ) X = c T dt dx Again, each side must be a negative constant. dt d X = = λ c T dt X dx We now have the pair of ODEs d X dt + λ X = 0 and + λ ct = 0 dx dt The general solutions are X x = A cos λx + B sin λx and T t = C cos λct + Dsin λct ( ) ( ) ( ) ( ) ( ) ( ) respectively, where A, B, C and D are arbitrary constants.

ENGI 940 8.0 - Finite String Page 8.0 Example 8.0.3 (continued) Consider the boundary conditions: y 0, t = X 0 T t = 0 t 0 ( ) ( ) ( ) For a non-trivial solution, this requires X ( ) ( ) ( ) ( ) ( ) 0 = 0 A= 0. y, t = X T t = 0 t 0 X = 0 nπ Bsin = 0 n =, n ( λ ) λ ( ) We now have a solution only for a discrete set of eigenvalues λ n, with corresponding eigenfunctions nπ x Xn ( x) = sin, ( n =,, 3, ) and nπ x yn( x, t) = Xn( x) Tn( t) = sin Tn( t), ( n =,, 3, ) So far, the solution has been identical to Example 8.0.. Consider the initial condition y(x, 0) = 0 : y x,0 = 0 X xt 0 = 0 x T 0 = 0 ( ) ( ) ( ) ( ) The initial value problem for T(t) is now nπ T + λ ct = 0, T( 0) = 0, where λ = the solution to which is nπ ct Tn( t) = Cnsin, ( n ) Our eigenfunctions for y are now nπ x nπct yn xt, = Xn xtn t = Cnsin sin, n ( ) ( ) ( ) ( )

ENGI 940 8.0 - Finite String Page 8. Example 8.0.3 (continued) Differentiate term by term and impose the initial velocity condition: y n π c n sin n π x = C = g( x) ( x,0) n = which is just the Fourier sine series expansion for the function g(x). The coefficients of the expansion are nπc nπu Cn = g ( u) sin du 0 which leads to the complete solution n = nπu nπ x nπct y ( x, t) = g ( u) sin du sin sin πc n 0 This solution is valid for any initial velocity function g(x) that is continuous with a piece-wise continuous derivative on [0, ] with g(0) = g() = 0. The solutions for Examples 8.0. and 8.0.3 may be superposed. et y (, ) et y (, ) Then y( xt, ) y( xt, ) y ( xt, ) xt be the solution for initial displacement f (x) and zero initial velocity. xt be the solution for zero initial displacement and initial velocity g(x). = + satisfies the wave equation (the sum of any two solutions of a linear homogeneous PDE is also a solution), and satisfies the boundary conditions y(0, t) = y(, t) = 0 : ( ) ( ) ( ) ( ) y x,0 = y x,0 + y x,0 = f x + 0, which satisfies the condition for initial displacement f (x). yt ( x,0 ) = yt( x,0 ) + yt( x,0) = 0+ g( x), which satisfies the condition for initial velocity g(x). Therefore the sum of the two solutions is the complete solution for initial displacement f (x) and initial velocity g(x): n = nπu nπ x nπct y ( x, t) = f ( u) sin du sin cos 0 nπu nπ x nπct + g ( u) sin du sin sin πc n 0 n =

ENGI 940 8.0 - Finite String Page 8. Example 8.0.4 An elastic string of length m is fixed at both ends (x = 0 and x = ). The string is initially in the shape of an arc of a parabola [y(x, 0) = x x for 0 < x < ] and is released with the initial velocity y ( x,0) ( 0 ) x x x =. It is known that the wave speed is c = 5 m s. Find the displacement y(x, t) at all locations on the string (0 < x < ) and at all subsequent times (t > 0). In the formula for the complete solution of the wave equation, nπu nπ x nπct y ( x, t) = f ( u) sin du sin cos 0 n = nπu nπ x nπct + g ( u) sin du sin sin πc n 0 n = we know that =, c = 5 and f (x) = g(x) = x x for 0 < x <. Both integrals inside the summations are the same: nπ u ( u u ) sin du = 0 ( ) u u u 3 cos nπu + sin nπu nπ ( nπ) ( nπ) n 0 3 ( ) 0 0 = + 3 + 0 ( nπ) ( nπ) ( nπ ) 3 n ( ( ) ) ( nπ ) 3 ( n ) 0 even = = 4 et (odd n) = k k = ( n odd) The complete solution is 8 sin 5( k ) π t y( xt, ) = 3 3 sin( k ) πx cos5( k ) πt + π ( k ) 5 ( k ) π A Maple file for this solution is available at "www.engr.mun.ca/~ggeorge/940/demos/ex804.mws". 0

ENGI 940 8.03 - Infinite String Page 8.3 8.03 The Wave Equation Vibrating Infinite String Example 8.03. An elastic string of infinite length is displaced into the form y = f (x) and is released from rest. Find the displacement y(x, t) at all locations on the string x and at all subsequent times (t > 0). The boundary value problem for the displacement function y(x, t) is: y y = c for < x < and t > 0 Initial configuration of string: y(x, 0) = f (x) for x String released from rest: y ( x,0) = 0 for x We no longer have the additional boundary conditions of fixed endpoints. However, it is reasonable to insist upon a bounded solution. Separation of Variables (or Fourier Method) Attempt a solution of the form y(x, t) = X(x) T(t) Again we find the linked pair of ordinary differential equations X + ω X = 0 and T + ω ct = 0 If ω = 0 then X(x) = ax + b. However, for a bounded solution, we require a = 0. For other values of ω, X(x) = a cos ωx + b sin ωx, which is bounded for all x and all ω. The ω = 0 case is a special case of this solution. We have a continuum of eigenvalues ω with corresponding eigenfunctions Xω ( x) = aω cos( ωx) + bωsin ( ωx) It then follows that T ( t) = c cosωct + d sinωct ω ω ω Imposing the initial condition of zero velocity, y = X ( x) T ( 0) = X ( x) dωωc = 0 x dω = 0 ( x,0)

ENGI 940 8.03 - Infinite String Page 8.4 Example 8.03. (continued) Therefore we have, for any real ω, a solution of the wave equation and the initial velocity condition, yω ( x, t) = Xω ( x) Tω ( t) = ( aω cos( ωx) + bωsin ( ωx) ) cos( ωct ) [where c ω has been absorbed into the other arbitrary constants a ω and b ω.] The superposition of solutions now leads to an integral, not a discrete sum. 0 ω 0 ω ω (, ) = (, ) = ( cos( ) + sin ( )) cos( ) y xt y xt dω a ωx b ωx ωct dω Imposing the remaining condition, (,0) = ( cos( ω ) + sin ( ω )) ω = ( ) y x aω x bω x d f x 0 But this is just the Fourier integral representation of f (x) on (, ). Therefore a ω and b ω are just the Fourier integral coefficients aω = f ( u) cos( ωu) du and bω f ( u) sin ( ωu) du π = π The complete solution is y ( x, t) = f ( u) cos( ωu) du cos( ωx) 0 π + f ( u) sin ( ωu) du sin ( ωx) cos( ωct ) dω π which, after re-iteration (interchanging the order of integration) is y ( x, t) = ( cos( ωu) cos( ωx) sin ( ωu) sin ( ωx) ) f ( u) cos ( ωct ) dωdu π + 0 y ( x, t) = f ( u) cos( ω( u x) ) cos( ωct ) dωdu π 0 = f ( u) ( cos( ω( u x ct )) cos( ω( u x ct ))) dω du π + + 0 ω = sin ( ω( u x + ct )) sin ( ω( u x ct )) = f ( u) + du π u x + ct u x ct ω = 0

ENGI 940 8.04 - d Alembert Solution Page 8.5 8.04 d Alembert Solution One form of the solution to Example 8.03., ω = sin ( ω( u x + ct )) sin ( ω( u x ct )) y ( x, t) = f ( u) + du π u x + ct u x ct ω = 0 suggests that, in general, one might seek solutions to the wave equation of the form f ( x + ct ) + f ( x ct ) y( xt, ) = f ( r) + f ( s) et r = x + ct and s = x ct, then yrs (, ) = and y y r y s = + = (( f ( r) + 0) + ( 0+ f ( s) ) ), r s = = + = + y y r y s = + = (( f ( r) + 0) c + ( 0+ f ( s) ) ( c) ), r s y y y r y s x x x r x x s x x ( f ( r) f ( s) ) y y r y s = + = ( cf ( r) c cf ( s) ( c) ), r s y y = ( f ( r) + f ( s) ) ( c f ( r) + c f ( s) ) = 0, c c ( + ) + ( ) f x ct f x ct Therefore y( xt, ) = is a solution to the wave equation for all twice differentiable functions f (u). This is part of the d Alembert solution. This d Alembert solution satisfies the initial displacement condition: f ( x+ 0) + f ( x 0) y( x,0) = = f ( x) cf ( x+ ct) cf ( x ct) cf ( x) cf ( x) Also y( xt, ) = = = 0 t t = 0 t = 0 The d Alembert solution therefore satisfies both initial conditions.,

ENGI 940 8.04 - d Alembert Solution Page 8.6 A more general d Alembert solution to the wave equation for an infinitely long string is (, ) y x t ( + ) + ( ) f x ct f x ct = + c x+ ct x ct ( ) g u du This satisfies the wave equation y y = c for < x < and t > 0 and Initial configuration of string: y(x, 0) = f (x) for x and Initial speed of string: y ( x,0) for any twice differentiable functions f (x) and g(x). ( ) = g x for x Physically, this represents two identical waves, moving with speed c in opposite directions along the string. x+ ct = g u du satisfies both initial conditions: c x ct x+ ct x y( x, t) = g( u) du y( x,0) = g( u) du = 0 c x ct c x Proof that y( x, t) ( ) Using a eibnitz differentiation of the integral: x+ ct y = g ( x + ct ) ( x + ct ) g ( x ct ) ( x ct ) + g ( u) du c x ct g ( x + ct ) + g ( x ct ) = ( c g ( x + ct ) + c g ( x ct ) + 0) = c y g( x+ 0) + g( x 0) = = g( x) t t = 0

ENGI 940 8.04 - d Alembert Solution Page 8.7 Example 8.04. An elastic string of infinite length is displaced into the form y = cos π x/ on [, ] only (and y = 0 elsewhere) and is released from rest. Find the displacement y(x, t) at all locations on the string x and at all subsequent times (t > 0). For this solution to the wave equation we have initial conditions π x cos ( x ) y( x,0) = f ( x) = 0 ( otherwise) and y ( x,0) = g ( x ) = 0 The d Alembert solution is f ( x + ct ) + f ( x ct ) x+ ct f ( x + ct ) + f ( x ct ) y( x, t) = + g( u) du = + 0 c x ct π ( x + ct ) cos ( ct x ct ) where f ( x + ct ) = 0 ( otherwise) ( x ct ) π cos ( + ct x + ct ) and f ( x ct ) = 0 ( otherwise) We therefore obtain two waves, each of the form of a single half-period of a cosine function, moving apart from a superposed state at x = 0 at speed c in opposite directions. See the web page "www.engr.mun.ca/~ggeorge/940/demos/ex804.html" for an animation of this solution.

ENGI 940 8.04 - d Alembert Solution Page 8.8 Example 8.04. (continued) Some snapshots of the solution are shown here:

ENGI 940 8.04 - d Alembert Solution Page 8.9 A more general case of a d Alembert solution arises for the homogeneous PDE with constant coefficients u u u A + B + C = 0 y y The characteristic (or auxiliary) equation for this PDE is Aλ + Bλ + C = 0 This leads to the complementary function (which is also the general solution for this homogeneous PDE) u x, y = f y+ λ x + f y+ λ x, ( ) ( ) ( ) where B D and B + λ = λ D = A A and D = B 4AC and f, f are arbitrary twice-differentiable functions of their arguments. λ and λ are the roots (or eigenvalues) of the characteristic equation. In the event of equal roots, the solution changes to uxy (, ) = f( y+ λx) + hxy (, ) f( y+ λx) where h(x, y) is any non-trivial linear function of x and/or y (except y + λx). The wave equation is a special case with y = t, A =, B = 0, C = /c and λ = ± /c.

ENGI 940 8.04 - d Alembert Solution Page 8.0 Example 8.04. u u u 3 + y y = 0 u(x, 0) = x u y (x, 0) = 0 (a) Classify the partial differential equation. (b) Find the value of u at (x, y) = (0, ). (a) Compare this PDE to the standard form u u u A B C y y + + = 0 A =, B = 3, C = D = 9 4 = > 0 Therefore the PDE is hyperbolic everywhere. (b) + 3 ± λ = = or u(x, y) = f (y + x) + g(y + x) u y (x, y) = f '(y + x) + g'(y + x) Boundary conditions: u(x, 0) = f (x) + g(x) = x () and u y (x, 0) = f '(x) + g'(x) = 0 () ( ) ( ) ( ) = f x + g x = x (3) (3) () g'(x) = x g'(x) = x ( ) ( ) ( ) g x = x + k g y+ x = y+ x + k

ENGI 940 8.04 - d Alembert Solution Page 8. Example 8.04. (continued) Also () f (x) = x g(x) = x + ½(x) k = x k f (y + x) = (y + x) k Therefore u(x, y) = f (y + x) + g(y + x) = (y + x) k (y + x) / + k = + + = ( y 4xy x y 4xy 4x ) ( y x ) The complete solution is therefore uxy (, ) = ( y x ) u ( 0,) = ( 0 ) = [It is easy (though tedious) to confirm that uxy (, ) = ( y x ) satisfies the partial u u u differential equation 3 + = 0 together with both initial conditions y y u(x, 0) = x and u y (x, 0) = 0.] [Also note that the arbitrary constants of integration for f and g cancelled each other out. This cancellation happens generally for this d Alembert method of solution.]

ENGI 940 8.04 - d Alembert Solution Page 8. Example 8.04.3 Find the complete solution to u u u 6 5 + = 4, y y u(x, 0) = x +, u y (x, 0) = 4 6x. This PDE is non-homogeneous. For the particular solution, we require a function such that the combination of second partial derivatives resolves to the constant 4. It is reasonable to try a quadratic function of x and y as our particular solution. Try up = ax + bxy + cy up up = ax + by and = bx + cy y up up up = a, = b and = c y y up up up 6 5 + = a 5b+ c = 4 y y We have one condition on three constants, two of which are therefore a free choice. Choose b = 0 and c = a, then 4a = 4 c = a = Therefore a particular solution is up = x + y [But we could have chosen, for example, a = b = 0 and c = 7 instead u = y ] Complementary function: A = 6, B = 5, C = D = 5 4 6 = > 0 Therefore the PDE is hyperbolic everywhere. + 5 ± λ = = or 3 The complementary function is u xy, = f y+ x + g y+ x C ( ) ( ) ( ) 3 and the general solution is u x, y = f y+ x + g y+ x + x + y ( ) ( ) ( ) 3 P 7

ENGI 940 8.04 - d Alembert Solution Page 8.3 Example 8.04.3 (continued) (, ) = ( + ) + ( + ) + + u x y f y x g y x x y 3 u = f ( y+ 3 x) + g ( y+ x) + y y Imposing the two boundary conditions: ( ) ( ) ( ) u x,0 = f x + g x + x = x+ () 3 and u x,0 = f x + g x + 0 = 4 6x () y ( ) ( ) ( ) 3 = f 3 x + g x + x = (3) 3 ( ) ( ) ( ) () (3) ( 3 ) 3 f x x = x ( ) ( ) ( ) 3 3 4 4 6 4 f x = 6x = 8 x f x = 8x ( ) 9 f x = x + k () x g( x) x x f ( 3 x) x x 9 = + = + + k 9 g ( x) = 4( x) + k g( x) = 4x + k But u x, y = f y+ x + g y+ x + x + y ( ) ( ) ( ) 3 3 ( ) ( ) ( ) uxy, = 9 y+ x + k+ 4 y+ x + k+ x + y [again the arbitrary constants cancel - they can be omitted safely.] = 9y 6xy x + 4y + x + + x + y Therefore the complete solution is (, ) = + + 4 6 8 u x y x y xy y

ENGI 940 8.04 - d Alembert Solution Page 8.4 Example 8.04.3 - Alternative Treatment of the Particular Solution Find the complete solution to u u u 6 5 + = 4, y y u(x, 0) = x +, u y (x, 0) = 4 6x. This PDE is non-homogeneous. For the particular solution, we require a function such that the combination of second partial derivatives resolves to the constant 4. It is reasonable to try a quadratic function of x and y as our particular solution. Try u = ax + bxy+ cy P u u ax by and bx cy y P P = + = + up up up = a, = b and = c y y up up up 6 5 + = a 5b+ c = 4 y y We have one condition on three constants, two of which are therefore a free choice. et us leave the free choice unresolved for now. Complementary function: + 5 ± λ = = or 3 The complementary function is C (, ) = ( + ) + ( + ) u xy f y x g y x and the general solution is 3 (, ) = ( + ) + ( + ) + + + u x y f y x g y x ax bxy cy 3

ENGI 940 8.04 - d Alembert Solution Page 8.5 Example 8.04.3 (continued) (, ) = ( + ) + ( + ) + + + u x y f y x g y x ax bxy cy 3 u = f ( y+ 3 x) + g ( y+ x) + bx + cy y Imposing the two boundary conditions: ( ) ( ) ( ) u x,0 = f x + g x + ax + 0+ 0 = x+ (A) and 3 ( ) ( ) ( ) u x,0 = f x + g x + bx + 0 = 4 6x (B) y d dx 3 A = f 3 x + g x + ax = (C) 3 ( ) ( ) ( ) (B) (C) ( 3 ) ( 4 ) 4 6 4 3 f x + b a x = x f x = 3 4a b 6 x = 9 4a b 6 x f x = 9 4a b 6 x ( ) ( ) ( )( ) ( ) ( ) 3 3 ( a b ) 9 4 6 f ( x) = x + k ( 4a b 6) x (A) 9 g( x) = x+ ax f ( 3 x) = x+ ax k 9 g( x) = a+ 4a b 6 ( x) + ( x) k = ( b+ 6 6a) x + 4x+ k But (, ) = ( + 3 ) + ( + ) + + + 9( 4a b 6) ( ) u x y f y x g y x ax bxy cy (, ) uxy = y+ x + k ( )( ) ( ) 3 + b+ 6 6a y+ x + 4 y+ x + k + ax + bxy+ cy [again the arbitrary constants cancel - they can be omitted safely.] ( a b ) ( ) ( )( y xy x b 6 6a y xy x ) 9 4 6 = + + + + + + + 4y+ x+ + ax + bxy+ cy 3 9 4

ENGI 940 8.04 - d Alembert Solution Page 8.6 Example 8.04.3 (continued) ( ) ( ) 9 4 6 4 6 6 = + + + + ( 4a b 6) + ( b+ 6 6a) + a x + + 4 y + x + a 5b+ c 30 = y 6xy + 0x + 4y + x + a 5b+ c 30 4 30 But a 5b+ c = 4 = = 8 a b + b + a + c y a b b a b xy Therefore the complete solution is (, ) = + + 4 6 8 u x y x y xy y ( 34 ( 6) ( 6 6) ) and note how the values of the two free parameters have no effect whatsoever on the solution. In practice, assign values to the free parameters in the particular solution only after one of the two arbitrary functions from the complementary function has been determined - if possible, make that function zero. In the example above, when we found f ( x) = 9( 4a b 6) x, choose a = 0 and b= 6 f ( x) = 0 and, from a 5b+ c = 4, we also have c = 8. (A) ( ) ( ) ( ) g x = x+ ax f x = x+ 0 0 g x = 4x+ 3 The general solution then becomes u x, y = f y+ x + g y+ x + ax + bxy+ cy ( ) ( 3 ) ( ) (, ) = 0 + 4( + ) + + 0+ 6 8 u x y y x xy y as before, but much faster! (, ) = + + 4 6 8 u x y x y xy y

ENGI 940 8.04 - d Alembert Solution Page 8.7 Example 8.04.4 Find the complete solution to u u u + + y y = 0, u = 0 on x = 0, u = x on y =. A =, B =, C = D = 4 4 = 0 Therefore the PDE is parabolic everywhere. ± 0 λ = = or The complementary function (and general solution) is uxy (, ) = f( y x) + hxyg (, ) ( y x) where h(x, y) is any convenient non-trivial linear function of (x, y) except a multiple of (y x). Choosing, arbitrarily, h(x, y) = x, u x, y = f y x + x g y x ( ) ( ) ( ) Imposing the boundary conditions: u(0, y) = 0 f (y) + 0 = 0 Therefore the function f is identically zero, for any argument including (y x). We now have u(x, y) = x g(y x). u(x, ) = x x g( x) = x g( x) = x g(x) = x Therefore u(x, y) = x g(y x) = x ( (y x)) The complete solution is (, ) = ( + ) uxy xx y

ENGI 940 8.04 - d Alembert Solution Page 8.8 Example 8.04.4 - Alternative Treatment of the Complementary Function Find the complete solution to u u u + + y y = 0, u = 0 on x = 0, u = x on y =. A =, B =, C = D = 4 4 = 0 ± 0 λ = = or The complementary function (and general solution) is uxy (, ) = f( y x) + hxyg (, ) ( y x) where h(x, y) is any convenient non-trivial linear function of (x, y) except a multiple of (y x). The most general choice possible is h( x, y) = ax+ by, with the restriction a λb. u ( x, y) = f ( y x) + ( ax + by) g ( y x) u ( 0, y) f ( y) by g ( y) 0 f ( y) by g ( y) [note: choosing b 0 f ( y) 0 u ( x,) = f ( x) + ( ax + b) g ( x) = x = + = = (A) = =, as happened on the previous page] and Using (A), u x, = b x g x + ax+ b g x = x ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) b + bx + ax + b g x = x a+ b x g x = x x g( x) = a+ b (note that the restriction on h(x, y) ensures that a+ b 0 ). y g( y) = a+ b (B) y (A) f ( y) = by a + b The complete solution becomes u ( x, y) ( y x) ( y x) = b( y x) + ( ax + by) a+ b a+ b ( y x) x y+ = ( by + bx + ax + by ) = ( a+ b) x a+ b a + b (, ) = ( + ) uxy xx y

ENGI 940 8.04 - d Alembert Solution Page 8.9 Example 8.04.4 Extension (continued) It is easy to confirm that this solution is correct: (, ) ( ) ( ) ( ) uxy = xx y+ u 0, y = 0 0 y+ = 0 y and u x, = x x + = x x ( ) ( ) and u u u ( x, y) = x xy + x = x y and = x y u u u =, = and = 0 y y u u u + + = + 0 = 0, y y ( xy)

ENGI 940 8.04 - d Alembert Solution Page 8.30 Two-dimensional aplace Equation u u y + = 0 A = C =, B = 0 D = 0 4 < 0 This PDE is elliptic everywhere. 0 ± 4 λ = = ± j The general solution is (, ) = ( ) + ( + ) u x y f y jx g y jx where f and g are any twice-differentiable functions. A function f (x, y) is harmonic if and only if f = 0 everywhere inside a domain Ω. Example 8.04.5 Is u = e x sin y harmonic on? u x u x = e sin y and = e cos y y u x u x = e sin y and = e sin y y u u x x u = + = e sin y e sin y y = 0 x, y Therefore yes, u = e x sin y is harmonic on. ( )

ENGI 940 8.04 - d Alembert Solution Page 8.3 Example 8.04.6 Find the complete solution uxy (, ) to the partial differential equation additional information 3 ( 0, ) u y = y and u x x = 0 = 0 u = 0, given the The PDE is u u y u = + = 0 (which means that the solution uxy (, ) is an harmonic function). A = C = B = D = B AC = < The PDE is elliptic everywhere., 0 4 4 0 A.E.: λ + = 0 λ =± j uc x, y = f y jx + g y+ jx The PDE is homogeneous u xy = C.F.: ( ) ( ) ( ) P.S.: P (, ) 0 G.S.: u( x, y) = f ( y jx) + g( y+ jx) (, ) ( ) ( ) ux x y = j f y jx + jg y+ jx Using the additional information, 3 3 u 0, y = f y + g y = y g y = y f y g y = 3y f y and u 0, y = 0 = j f y + jg y = j f y + 3y f y ( ) ( ) ( ) ( ) ( ) ( ) ( ) x ( ) ( ) ( ) ( ) ( ) ( ) 3 3 f ( y ) 3y f ( y) y f ( y) y 3 3 3 g ( y) y y y = = = = = Therefore the complete solution is u x, y = y jx + y+ jx 3 3 ( ) ( ) ( ) 3 = y 3( jx) y 3 ( jx 3 ) y ( jx 3 + ) + y + 3( jx) y + 3 ( jx ) y + ( jx ) ( ) 3 3 = y + 3 j x y 3 (, ) = 3 uxy y xy Note that the solution is completely real, even though the eigenvalues are not real.

ENGI 940 8.05 - Max-Min Principle Page 8.3 8.05 The Maximum-Minimum Principle et Ω be some finite domain on which a function u(x, y) and its second derivatives are defined. et Ω be the union of the domain with its boundary. et m and M be the minimum and maximum values respectively of u on the boundary of the domain. If u 0 in Ω, then u is subharmonic and u r < M u r M r in Ω ( ) or ( ) If u 0 in Ω, then u is superharmonic and u r > m u r m r in Ω ( ) or ( ) If u = 0 in Ω, then u is harmonic (both subharmonic and superharmonic) and u is either constant on Ω or m < u < M everywhere on Ω. Example 8.05. u = 0 in Ω : x + y < and u(x, y) = on C : x + y =. Find u(x, y) on Ω. u is harmonic on Ω min ( u) C (, ) uxy on Ω max C ( u) But min C (u) = max C (u) = Therefore u(x, y) = everywhere in Ω.

ENGI 940 8.05 - Max-Min Principle Page 8.33 Example 8.05. u = 0 in the square domain Ω : < x < +, < y < +. On the boundary C, on the left and right edges (x = ±), u(x, y) = 4 y, while on the top and bottom edges (y = ±), u(x, y) = x 4. Find bounds on the value of u(x, y) inside the domain Ω. For < y < +, 0 < 4 y < 4. For < x < +, 4 < x 4 < 0. Therefore, on the boundary C of the domain Ω, 4 < u(x, y) < +4 so that m = 4 and M = +4. u(x, y) is harmonic (because u = 0). Therefore, everywhere in Ω, 4 < u(x, y) < +4 Note: u(x, y) = x y is consistent with the boundary condition and u = x 0, u = 0 y u =, u = = u y y u u u = + = 0 y Contours of constant values of u are hyperbolas. A contour map illustrates that 4 < u(x, y) < +4 within the domain is indeed true.

ENGI 940 8.06 - Heat Equation Page 8.34 8.06 The Heat Equation For a material of constant density ρ, constant specific heat μ and constant thermal conductivity K, the partial differential equation governing the temperature u at any location (x, y, z) and any time t is u K = k u, where k = µρ Example 8.06. Heat is conducted along a thin homogeneous bar extending from x = 0 to x =. There is no heat loss from the sides of the bar. The two ends of the bar are maintained at temperatures T (at x = 0) and T (at x = ). The initial temperature throughout the bar at the cross-section x is f (x). Find the temperature at any point in the bar at any subsequent time. The partial differential equation governing the temperature u(x, t) in the bar is u = k u together with the boundary conditions u( 0, t) T u t, = T and the initial condition u(x, 0) = f (x) = and ( ) [Note that if an end of the bar is insulated, instead of being maintained at a constant u u temperature, then the boundary condition changes to ( 0, t) = 0 or ( t, ) = 0.] Attempt a solution by the method of separation of variables. u(x, t) = X(x) T(t) T X XT = k X T = k = c T X Again, when a function of t only equals a function of x only, both functions must equal the same absolute constant. Unfortunately, the two boundary conditions cannot both be satisfied unless T = T = 0. Therefore we need to treat this more general case as a perturbation of the simpler ( T = T = 0) case.

ENGI 940 8.06 - Heat Equation Page 8.35 Example 8.06. (continued) et u(x, t) = v(x, t) + g(x) Substitute this into the PDE: v v ( v( xt, ) + g( x) ) = k v ( xt, ) + g( x) = k + g x This is the standard heat PDE for v if we choose g such that g"(x) = 0. g(x) must therefore be a linear function of x. ( ) ( ) We want the perturbation function g(x) to be such that u( 0, t) = T and u( t, ) = T and v(0, t) = v(, t) = 0 g 0 It follows that T T g( x) = x + T and we now have the simpler problem Therefore g(x) must be the linear function for which ( ) = T and g( ) T =. together with the boundary conditions v(0, t) = v(, t) = 0 and the initial condition v(x, 0) = f (x) g(x) v = k v Now try separation of variables on v(x, t) : v(x, t) = X(x) T(t) T X XT = k X T = = c kt X But v(0, t) = v(, t) = 0 X(0) = X() = 0 This requires c to be a negative constant, say λ. The solution is very similar to that for the wave equation on a finite string with fixed ends nπ (section 8.0). The eigenvalues are λ = and the corresponding eigenfunctions are any non-zero constant multiples of nπ x Xn ( x) = sin

ENGI 940 8.06 - Heat Equation Page 8.36 Example 8.06. (continued) The ODE for T(t) becomes nπ T + kt = 0 whose general solution is nπ kt/ Tn( t) = cne Therefore nπx n π kt vn( xt, ) = Xn( xt ) n( t) = cnsin exp If the initial temperature distribution f (x) g(x) is a simple multiple of sin n π x for nπx n π kt some integer n, then the solution for v is just v( xt, ) = cn sin exp. Otherwise, we must attempt a superposition of solutions. nπx n π kt v( xt, ) = cn sin exp n = nπ x such that v( x,0) = cn sin = f ( x) g( x). n = nπ z The Fourier sine series coefficients are cn = ( f ( z) g ( z) ) sin dz 0 so that the complete solution for v(x, t) is T T nπz nπx n π kt v( x, t) = f ( z) z T sin dz sin exp n = 0 and the complete solution for u(x, t) is T T u( xt, ) = v( xt, ) + x + T Note how this solution can be partitioned into a transient part v(x, t) (which decays to zero as t increases) and a steady-state part g(x) which is the limiting value that the temperature distribution approaches.

ENGI 940 8.06 - Heat Equation Page 8.37 Example 8.06. (continued) As a specific example, let k = 9, T = 00, T = 00, = and f (x) = 45x 40x + 00, (for which f (0) = 00, f () = 00 and f (x) > 0 x). 00 00 Then g( x) = x + 00 = 50 x + 00 The Fourier sine series coefficients are (( nπ z cn = 45z 40z + 00) ( 50z + 00) ) sin dz 0 nπ z cn = 45 ( z z) sin dz 0 After an integration by parts (details omitted here), ( z ) ( ) 6 nπz 8 nπz cn = 45 ( z + z) + 3 cos + sin n π ( nπ) nπ 30 c 3 (( ) n n = ) nπ ( ) The complete solution is z= z= 0 ( ) π π 30 n x 9n t u( xt, ) = 50x + 00 sin exp 3 3 π n = n 4 Some snapshots of the temperature distribution (from the tenth partial sum) from the Maple file at "www.engr.mun.ca/~ggeorge/940/demos/ex806.mws" are shown on the next page. n

ENGI 940 8.06 - Heat Equation Page 8.38 Example 8.06. (continued) The steady state distribution is nearly attained in much less than a second! END OF CHAPTER 8 END OF ENGI. 940!